Changeset 1462


Ignore:
Timestamp:
Aug 15, 2014 12:45:04 PM (8 years ago)
Author:
vondreele
Message:

add sig-2 to TOF profile coeff.
new d-spacing plot option

Location:
trunk
Files:
7 edited

Legend:

Unmodified
Added
Removed

  • trunk/GSASII.py

    r1461 r1462  
    980980            elif 'T' in DataType:
    981981                names = ['Type','fltPath','2-theta','difC','difA', 'difB','Zero','alpha','beta-0','beta-1',
    982                     'beta-q','sig-0','sig-1','sig-q','X', 'Y','Azimuth',]
    983                 codes = [0,0,0,0,0, 0,0,0,0,0, 0,0,0,0,0, 0,0,]
     982                    'beta-q','sig-0','sig-1','sig-2','sig-q', 'X','Y','Azimuth',]
     983                codes = [0,0,0,0,0, 0,0,0,0,0, 0,0,0,0,0, 0,0,0,]
    984984                azm = 0.
    985985                if 'INS  1DETAZM' in Iparm:
    … …  
    998998                    s = Iparm['INS  1PRCF 1'].split()
    999999                    if abs(pfType) == 1:
    1000                         data.extend([G2IO.sfloat(s[1]),G2IO.sfloat(s[2]),G2IO.sfloat(s[3])])
     1000                        data.extend([G2IO.sfloat(s[1]),G2IO.sfloat(s[2]),G2IO.sfloat(s[3])]) #alpha, beta-0, beta-1
    10011001                        s = Iparm['INS  1PRCF 2'].split()
    1002                         data.extend([0.0,0.0,G2IO.sfloat(s[1]),0.0,0.0,0.0,azm])
     1002                        data.extend([0.0,0.0,G2IO.sfloat(s[1]),G2IO.sfloat(s[2]),0.0,0.0,0.0,azm])    #beta-q, sig-0, sig-1, sig-2, sig-q, X, Y
    10031003                    elif abs(pfType) in [3,4,5]:
    1004                         data.extend([G2IO.sfloat(s[0]),G2IO.sfloat(s[1]),G2IO.sfloat(s[2])])
     1004                        data.extend([G2IO.sfloat(s[0]),G2IO.sfloat(s[1]),G2IO.sfloat(s[2])]) #alpha, beta-0, beta-1
    10051005                        if abs(pfType) == 4:
    1006                             data.extend([0.0,0.0,G2IO.sfloat(s[3]),0.0,0.0,0.0,azm])
     1006                            data.extend([0.0,0.0,G2IO.sfloat(s[3]),0.0,0.0,0.0,0.0,azm])    #beta-q, sig-0, sig-1, sig-2, sig-q, X, Y
    10071007                        else:
    10081008                            s = Iparm['INS  1PRCF 2'].split()
    1009                             data.extend([0.0,0.0,G2IO.sfloat(s[0]),0.0,0.0,0.0,azm])                       
     1009                            data.extend([0.0,0.0,G2IO.sfloat(s[0]),G2IO.sfloat(s[1]),0.0,0.0,0.0,azm])    #beta-q, sig-0, sig-1, sig-2, sig-q, X, Y                       
    10101010                else:
    10111011                    s = Iparm['INS  1PRCF1 '].split()
    … …  
    10131013                    s = Iparm['INS  1PRCF11'].split()
    10141014                    if abs(pfType) == 1:
    1015                         data.extend([G2IO.sfloat(s[1]),G2IO.sfloat(s[2]),G2IO.sfloat(s[3])])
     1015                        data.extend([G2IO.sfloat(s[1]),G2IO.sfloat(s[2]),G2IO.sfloat(s[3])]) #alpha, beta-0, beta-1
    10161016                        s = Iparm['INS  1PRCF12'].split()
    1017                         data.extend([0.0,0.0,G2IO.sfloat(s[1]),0.0,0.0,0.0,azm])
     1017                        data.extend([0.0,0.0,G2IO.sfloat(s[1]),G2IO.sfloat(s[2]),0.0,0.0,0.0,azm])    #beta-q, sig-0, sig-1, sig-2, sig-q, X, Y
    10181018                    elif abs(pfType) in [3,4,5]:
    1019                         data.extend([G2IO.sfloat(s[0]),G2IO.sfloat(s[1]),G2IO.sfloat(s[2])])
     1019                        data.extend([G2IO.sfloat(s[0]),G2IO.sfloat(s[1]),G2IO.sfloat(s[2])]) #alpha, beta-0, beta-1
    10201020                        if abs(pfType) == 4:
    1021                             data.extend([0.0,0.0,G2IO.sfloat(s[3]),0.0,0.0,0.0,azm])
     1021                            data.extend([0.0,0.0,G2IO.sfloat(s[3]),0.0,0.0,0.0,0.0,azm])    #beta-q, sig-0, sig-1, sig-2, sig-q, X, Y
    10221022                        else:
    10231023                            s = Iparm['INS  1PRCF12'].split()
    1024                             data.extend([0.0,0.0,G2IO.sfloat(s[0]),0.0,0.0,0.0,azm])                       
     1024                            data.extend([0.0,0.0,G2IO.sfloat(s[0]),G2IO.sfloat(s[1]),0.0,0.0,0.0,azm])    #beta-q, sig-0, sig-1, sig-2, sig-q, X, Y                       
    10251025                Inst1 = G2IO.makeInstDict(names,data,codes)
    10261026                Inst2 = {}
    … …  
    19171917        self.logPlot = False
    19181918        self.qPlot = False
     1919        self.dPlot = False
    19191920        self.sqPlot = False
    19201921        self.SqrtPlot = False

  • trunk/GSASIImath.py

    r1453 r1462  
    20492049   
    20502050    '''
    2051     return ins['sig-0']+ins['sig-1']*dsp**2+ins['sig-q']*dsp
     2051    return ins['sig-0']+ins['sig-1']*dsp**2+ins['sig-2']*dsp**4+ins['sig-q']*dsp
    20522052   
    20532053def getTOFsigDeriv(dsp):
    … …  
    20592059   
    20602060    '''
    2061     return 1.0,dsp**2,dsp
     2061    return 1.0,dsp**2,dsp**4,dsp
    20622062   
    20632063def getTOFgamma(ins,dsp):
    … …  
    21622162            dsp = pos/Parms['difC'][1]
    21632163        if 'Pdabc' in Parms2:
    2164             for x in ['sig-0','sig-1','sig-q','X','Y']:
     2164            for x in ['sig-0','sig-1','sig-2','sig-q','X','Y']:
    21652165                ins[x] = Parms[x][ind]
    21662166            Pdabc = Parms2['Pdabc'].T
    … …  
    21682168            bet = np.interp(dsp,Pdabc[0],Pdabc[2])
    21692169        else:
    2170             for x in ['alpha','beta-0','beta-1','beta-q','sig-0','sig-1','sig-q','X','Y']:
     2170            for x in ['alpha','beta-0','beta-1','beta-q','sig-0','sig-1','sig-2','sig-q','X','Y']:
    21712171                ins[x] = Parms[x][ind]
    21722172            alp = getTOFalpha(ins,dsp)

  • trunk/GSASIIplot.py

    r1460 r1462  
    910910                newPlot = True
    911911                G2frame.qPlot = not G2frame.qPlot
     912                G2frame.dPlot = False
    912913            elif 'SASD' in plottype:
    913914                newPlot = True
    914                 G2frame.sqPlot = not G2frame.sqPlot       
     915                G2frame.sqPlot = not G2frame.sqPlot
     916        elif event.key == 't' and 'PWDR' in plottype:
     917            G2frame.dPlot = not G2frame.dPlot
     918            G2frame.qPlot = False
     919            newPlot = True     
    915920        elif event.key == 'm':
    916921            G2frame.SqrtPlot = False
    … …  
    943948                if 'C' in Parms['Type'][0]:
    944949                    wave = G2mth.getWave(Parms)
    945                     if G2frame.qPlot and 'PWDR' in plottype:
     950                    if (G2frame.qPlot or G2frame.dPlot) and 'PWDR' in plottype:
    946951                        try:
    947952                            xpos = 2.0*asind(xpos*wave/(4*math.pi))
    … …  
    10241029            if ind.all() != [0] and ObsLine[0].get_label() in str(pick):                                    #picked a data point
    10251030                data = G2frame.PatternTree.GetItemPyData(G2frame.PickId)
    1026                 if 'C' in Parms['Type'][0]:                            #CW data - TOF later in an elif
     1031                if 'C' in Parms['Type'][0]:
    10271032                    if G2frame.qPlot:                              #qplot - convert back to 2-theta
    10281033                        xy[0] = 2.0*asind(xy[0]*wave/(4*math.pi))
     1034                    elif G2frame.dPlot:                            #dplot - convert back to 2-theta
     1035                        xy[0] = 2.0*asind(wave/(2.*xy[0]))
     1036                elif 'T' in Parms['Type'][0]:
     1037                    if G2frame.qPlot:                              #qplot - convert back to TOF
     1038                        xy[0] = 2.0*math.pi*difC/xy[0]
     1039                    elif G2frame.dPlot:                            #dplot - convert back to TOF
     1040                        xy[0] *= difC                   
    10291041                XY = G2mth.setPeakparms(Parms,Parms2,xy[0],xy[1])
    10301042                data['peaks'].append(XY)
    … …  
    10411053                    if G2frame.qPlot and 'PWDR' in plottype:                              #qplot - convert back to 2-theta
    10421054                        xy[0] = 2.0*asind(xy[0]*wave/(4*math.pi))
     1055                    elif G2frame.dPlot:                            #dplot - convert back to 2-theta
     1056                        xy[0] = 2.0*asind(wave/(2.*xy[0]))
     1057                elif 'T' in Parms['Type'][0]:
     1058                    if G2frame.qPlot:                              #qplot - convert back to TOF
     1059                        xy[0] = 2.0*math.pi*difC/xy[0]
     1060                    elif G2frame.dPlot:                            #dplot - convert back to TOF
     1061                        xy[0] *= difC                   
    10431062                if G2frame.ifGetExclude:
    10441063                    excl = [0,0]
    … …  
    11071126                    else:
    11081127                        data[id][id2] = 2*math.pi*Parms['difC'][1]/xpos
     1128                elif G2frame.dPlot and 'PWDR' in plottype:
     1129                    if 'C' in Parms['Type'][0]:
     1130                        data[id][id2] = 2.0*asind(wave/(2*xpos))
     1131                    else:
     1132                        data[id][id2] *= Parms['difC'][1]
    11091133                else:
    11101134                    data[id][id2] = xpos
    … …  
    11231147                else:
    11241148                    if G2frame.qPlot:
    1125                         data['peaks'][lineNo-2][0] = 2.0*asind(wave*xpos/(4*math.pi))
     1149                        if 'C' in Parms['Type'][0]:
     1150                            data['peaks'][lineNo-2][0] = 2.0*asind(wave*xpos/(4*math.pi))
     1151                        else:
     1152                            data['peaks'][lineNo-2][0] = 2*math.pi*Parms['difC'][1]/xpos
     1153                    elif G2frame.dPlot:
     1154                        if 'C' in Parms['Type'][0]:
     1155                            data['peaks'][lineNo-2][0] = 2.0*asind(wave/(2*xpos))
     1156                        else:
     1157                            data['peaks'][lineNo-2][0] *= Parms['difC'][1]
    11261158                    else:
    11271159                        data['peaks'][lineNo-2][0] = xpos
    … …  
    11901222                if G2frame.SinglePlot:
    11911223                    Page.Choice = (' key press','n: log(I) off',
    1192                         'c: contour on','q: toggle q plot','m: toggle multidata plot','w: toggle divide by sig','+: no selection')
     1224                        'c: contour on','q: toggle q plot','t: toggle d-spacing plot',
     1225                            'm: toggle multidata plot','w: toggle divide by sig','+: no selection')
    11931226                else:
    11941227                    Page.Choice = (' key press','n: log(I) off',
    11951228                        'd: offset down','l: offset left','r: offset right','u: offset up','o: reset offset',
    1196                         'c: contour on','q: toggle q plot','m: toggle multidata plot','w: toggle divide by sig','+: no selection')
     1229                        'c: contour on','q: toggle q plot','t: toggle d-spacing plot',
     1230                        'm: toggle multidata plot','w: toggle divide by sig','+: no selection')
    11971231            elif 'SASD' in plottype:
    11981232                if G2frame.SinglePlot:
    … …  
    12081242                    Page.Choice = (' key press',
    12091243                        'b: toggle subtract background','n: log(I) on','s: toggle sqrt plot','c: contour on',
    1210                         'q: toggle q plot','m: toggle multidata plot','w: toggle divide by sig','+: no selection')
     1244                        'q: toggle q plot','t: toggle d-spacing plot','m: toggle multidata plot',
     1245                        'w: toggle divide by sig','+: no selection')
    12111246                else:
    12121247                    Page.Choice = (' key press','l: offset left','r: offset right','d: offset down',
    12131248                        'u: offset up','o: reset offset','b: toggle subtract background','n: log(I) on','c: contour on',
    1214                         'q: toggle q plot','m: toggle multidata plot','w: toggle divide by sig','+: no selection')
     1249                        'q: toggle q plot','t: toggle d-spacing plot','m: toggle multidata plot',
     1250                        'w: toggle divide by sig','+: no selection')
    12151251            elif 'SASD' in plottype:
    12161252                if G2frame.SinglePlot:
    … …  
    12831319    if G2frame.qPlot or 'SASD' in plottype:
    12841320        Plot.set_xlabel(r'$Q, \AA^{-1}$',fontsize=16)
     1321    elif G2frame.dPlot and 'PWDR' in plottype:
     1322        Plot.set_xlabel(r'$d, \AA$',fontsize=16)
    12851323    else:
    12861324        if 'C' in ParmList[0]['Type'][0]:       
    … …  
    13391377            else:
    13401378                X = 2*np.pi*Parms['difC'][1]/(xye[0]-Zero)
     1379        elif G2frame.dPlot and 'PWDR' in plottype:
     1380            Id = G2gd.GetPatternTreeItemId(G2frame,G2frame.root, Pattern[2])
     1381            if 'C' in Parms['Type'][0]:
     1382                X = wave/(2.*npsind(xye[0]-Zero)/2.0)
     1383            else:
     1384                X = (xye[0]-Zero)/Parms['difC'][1]
    13411385        else:
    13421386            X = xye[0]-Zero
    … …  
    13631407                else:
    13641408                    limits = 2*np.pi*difC/limits
     1409            elif G2frame.dPlot and 'PWDR' in plottype:
     1410                if 'C' in Parms['Type'][0]:
     1411                    limits = wave/(2*npsind(limits/2.0))
     1412                else:
     1413                    limits /= difC               
    13651414            Lines.append(Plot.axvline(limits[1][0],color='g',dashes=(5,5),picker=3.))   
    13661415            Lines.append(Plot.axvline(limits[1][1],color='r',dashes=(5,5),picker=3.))
    … …  
    14821531                            else:
    14831532                                Lines.append(Plot.axvline(2*math.pi*difC/item[0],color=colors[N%6],picker=2.))                               
     1533                        if G2frame.dPlot:
     1534                            if 'C' in Parms['Type'][0]:
     1535                                Lines.append(Plot.axvline(wave/(2*sind(item[0]/2.)),color=colors[N%6],picker=2.))
     1536                            else:
     1537                                Lines.append(Plot.axvline(item[0]/difC,color=colors[N%6],picker=2.))                               
    14841538                        else:
    14851539                            Lines.append(Plot.axvline(item[0],color=colors[N%6],picker=2.))
    … …  
    15141568            for peak in peaks:
    15151569                if G2frame.qPlot:
    1516                     Plot.axvline(4*np.pi*sind(peak[0]/2.0)/wave,color='b')
     1570                    if 'C' in Parms['Type'][0]:
     1571                        Plot.axvline(4*np.pi*sind(peak[0]/2.0)/wave,color='b')
     1572                    else:
     1573                        Plot.axvline(2*np.pi*Parms['difC'][1]/peak[0],color='b')
     1574                if G2frame.dPlot:
     1575                    if 'C' in Parms['Type'][0]:
     1576                        Plot.axvline(wave/(2*sind(peak[0]/2.0)),color='b')
     1577                    else:
     1578                        Plot.avline(peak[0]/Parms['difC'][1],color='b')
    15171579                else:
    15181580                    Plot.axvline(peak[0],color='b')
    15191581            for hkl in G2frame.HKL:
    15201582                if G2frame.qPlot:
    1521                     Plot.axvline(4*np.pi*sind(hkl[5]/2.0)/wave,color='r',dashes=(5,5))
     1583                    if 'C' in Parms['Type'][0]:
     1584                        Plot.axvline(4*np.pi*sind(hkl[5]/2.0)/wave,color='b')
     1585                    else:
     1586                        Plot.axvline(2*np.pi*Parms['difC'][1]/hkl[5],color='b')
     1587                if G2frame.dPlot:
     1588                    if 'C' in Parms['Type'][0]:
     1589                        Plot.axvline(wave/(2*sind(hkl[5]/2.0)),color='b')
     1590                    else:
     1591                        Plot.avline(hkl[5]/Parms['difC'][1],color='b')
    15221592                else:
    15231593                    Plot.axvline(hkl[5],color='r',dashes=(5,5))
    … …  
    15371607                if G2frame.qPlot:
    15381608                    Plot.plot(2*np.pi/peak.T[0],pos,refColors[pId%6]+'|',mew=1,ms=8,picker=3.,label=phase)
     1609                elif G2frame.dPlot:
     1610                    Plot.plot(peak.T[0],pos,refColors[pId%6]+'|',mew=1,ms=8,picker=3.,label=phase)
    15391611                else:
    15401612                    Plot.plot(peak.T[1],pos,refColors[pId%6]+'|',mew=1,ms=8,picker=3.,label=phase)
    … …  
    22452317                '%.3f'%(2*theta)
    22462318            G2frame.G2plotNB.Delete('Peak Widths')
    2247     else:
     2319    else:   #'T'OF
    22482320        Plot.set_title('Instrument and sample peak coefficients')
    22492321        Plot.set_xlabel(r'$Q, \AA^{-1}$',fontsize=14)

  • trunk/GSASIIpwd.py

    r1459 r1462  
    480480    'needs a doc string'
    481481    sig = lambda Th,U,V,W: 1.17741*math.sqrt(max(0.001,U*tand(Th)**2+V*tand(Th)+W))*math.pi/180.
    482     sigTOF = lambda dsp,S0,S1,Sq:  S0+S1*dsp**2+Sq*dsp
     482    sigTOF = lambda dsp,S0,S1,S2,Sq:  S0+S1*dsp**2+S2*dsp**4+Sq*dsp
    483483    gam = lambda Th,X,Y: (X/cosd(Th)+Y*tand(Th))*math.pi/180.
    484484    gamTOF = lambda dsp,X,Y: X*dsp+Y*dsp**2
    … …  
    488488    else:
    489489        dsp = pos/Inst['difC'][0]
    490         s = sigTOF(dsp,Inst['sig-0'][1],Inst['sig-1'][1],Inst['sig-q'][1])
     490        s = sigTOF(dsp,Inst['sig-0'][1],Inst['sig-1'][1],Inst['sig-2'][1],Inst['sig-q'][1])
    491491        g = gamTOF(dsp,Inst['X'][1],Inst['Y'][1])
    492492    return getgamFW(g,s)
    … …  
    10131013                if sigName in varyList:
    10141014                    sig = parmDict[sigName]
    1015                     dsds0 = dsds1 = dsds2 = 0
     1015                    dsds0 = dsds1 = dsds2 = dsds3 = 0
    10161016                else:
    10171017                    sig = G2mth.getTOFsig(parmDict,dsp)
    1018                     dsds0,dsds1,dsds2 = G2mth.getTOFsigDeriv(dsp)
     1018                    dsds0,dsds1,dsds2,dsds3 = G2mth.getTOFsigDeriv(dsp)
    10191019                gamName = 'gam'+str(iPeak)
    10201020                if gamName in varyList:
    … …  
    10631063                if 'sig-1' in varyList:
    10641064                    dMdv[varyList.index('sig-1')] += dsds1*dervDict['sig']
     1065                if 'sig-2' in varyList:
     1066                    dMdv[varyList.index('sig-2')] += dsds2*dervDict['sig']
    10651067                if 'sig-q' in varyList:
    1066                     dMdv[varyList.index('sig-q')] += dsds2*dervDict['sig']
     1068                    dMdv[varyList.index('sig-q')] += dsds3*dervDict['sig']
    10671069                if 'X' in varyList:
    10681070                    dMdv[varyList.index('X')] += dsdX*dervDict['gam']
    … …  
    13021304            insVals.append(Inst[parm][1])
    13031305            if parm in ['U','V','W','X','Y','SH/L','I(L2)/I(L1)','alpha',
    1304                 'beta-0','beta-1','beta-q','sig-0','sig-1','sig-q',] and Inst[parm][2]:
     1306                'beta-0','beta-1','beta-q','sig-0','sig-1','sig-2','sig-q',] and Inst[parm][2]:
    13051307                    insVary.append(parm)
    13061308        instDict = dict(zip(insNames,insVals))
    … …  
    13441346        for parm in Inst:
    13451347            if parm in  ['U','V','W','X','Y','SH/L','I(L2)/I(L1)','alpha',
    1346                 'beta-0','beta-1','beta-q','sig-0','sig-1','sig-q',]:
     1348                'beta-0','beta-1','beta-q','sig-0','sig-1','sig-2','sig-q',]:
    13471349                ptlbls += "%s" % (parm.center(12))
    13481350                ptstr += ptfmt % (Inst[parm][1])

  • trunk/GSASIIpwdGUI.py

    r1461 r1462  
    929929        for key in keys:
    930930            if key in ['Type','U','V','W','X','Y','SH/L','I(L2)/I(L1)','alpha',
    931                 'beta-0','beta-1','beta-q','sig-0','sig-1','sig-q','Polariz.',
     931                'beta-0','beta-1','beta-q','sig-0','sig-1','sig-2','sig-q','Polariz.',
    932932                'Lam','Azimuth','2-theta','fltPath','difC','difA','difB','Zero','Lam1','Lam2']:
    933933                good.append(key)
    … …  
    10071007                File.close()
    10081008                Inst,Inst2 = G2frame.PatternTree.GetItemPyData(G2gd.GetPatternTreeItemId(G2frame,G2frame.PatternId,'Instrument Parameters'))
    1009                 inst = G2IO.makeInstDict(newItems,newVals,len(newVals)*[False,])
    1010                 G2frame.PatternTree.SetItemPyData(G2gd.GetPatternTreeItemId(G2frame,G2frame.PatternId,'Instrument Parameters'),[inst,Inst2])
     1009                data = G2IO.makeInstDict(newItems,newVals,len(newVals)*[False,])
     1010                G2frame.PatternTree.SetItemPyData(G2gd.GetPatternTreeItemId(G2frame,G2frame.PatternId,'Instrument Parameters'),[data,Inst2])
    10111011                RefreshInstrumentGrid(event,doAnyway=True)          #to get peaks updated
    10121012                UpdateInstrumentGrid(G2frame,data)
    … …  
    12731273                refFlgElem.append(None)                   
    12741274                if 'Pdabc' in Inst2:
    1275                     Items = ['sig-0','sig-1','X','Y']
     1275                    Items = ['sig-0','sig-1','sig-2','sig-q','X','Y']
    12761276                    subSizer.Add(wx.StaticText(G2frame.dataDisplay,-1,'  difC: '),0,WACV)
    12771277                    txt = '%8.2f'%(insVal['difC'])
    … …  
    12831283                    subSizer.Add(wx.StaticText(G2frame.dataDisplay,-1,'  alpha, beta: fixed by table'),0,WACV)
    12841284                else:
    1285                     Items = ['difC','difA','difB','Zero','alpha','beta-0','beta-1','beta-q','sig-0','sig-1','sig-q','','X','Y']
     1285                    Items = ['difC','difA','difB','Zero','alpha','beta-0','beta-1','beta-q','sig-0','sig-1','sig-2','sig-q','X','Y']
    12861286                mainSizer.Add((5,5),0)
    12871287                mainSizer.Add(subSizer)
    … …  
    13931393                insDef['Azimuth'] = 0.0
    13941394                insRef['Azimuth'] = False
    1395         if 'T' in insVal['Type']:
    1396             if 'difB' not in insVal:
    1397                 insVal['difB'] = 0.0
    1398                 insDef['difB'] = 0.0
    1399                 insRef['difB'] = False
     1395#        if 'T' in insVal['Type']:
     1396#            if 'difB' not in insVal:
     1397#                insVal['difB'] = 0.0
     1398#                insDef['difB'] = 0.0
     1399#                insRef['difB'] = False
    14001400    #end of patch
    14011401    if 'P' in insVal['Type']:                   #powder data menu commands

  • trunk/GSASIIstrMain.py

    r1461 r1462  
    4242   
    4343ateln2 = 8.0*math.log(2.0)
    44 DEBUG = False
     44DEBUG = True
    4545
    4646def RefineCore(Controls,Histograms,Phases,restraintDict,rigidbodyDict,parmDict,varyList,

  • trunk/GSASIIstrMath.py

    r1461 r1462  
    11851185        for iPO in dIdPO:
    11861186            dIdPO[iPO] /= POcorr
    1187     AbsCorr = GetAbsorb(refl,hfx,calcControls,parmDict)
    1188     dFdAb = GetAbsorbDerv(refl,hfx,calcControls,parmDict)/AbsCorr
     1187    dFdAb = GetAbsorbDerv(refl,hfx,calcControls,parmDict)*refl[14]/refl[16] #wave/abs corr
    11891188    dFdEx = GetPwdrExtDerv(refl,pfx,phfx,hfx,calcControls,parmDict)
    11901189    return dIdsh,dIdsp,dIdPola,dIdPO,dFdODF,dFdSA,dFdAb,dFdEx
    … …  
    16121611               
    16131612    def GetReflSigGamTOF(refl,G,GB,phfx,calcControls,parmDict):
    1614         sig = parmDict[hfx+'sig-0']+parmDict[hfx+'sig-1']*refl[4]**2+parmDict[hfx+'sig-q']*refl[4]
     1613        sig = parmDict[hfx+'sig-0']+parmDict[hfx+'sig-1']*refl[4]**2+   \
     1614            parmDict[hfx+'sig-2']*refl[4]**4+parmDict[hfx+'sig-q']*refl[4]
    16151615        gam = parmDict[hfx+'X']*refl[4]+parmDict[hfx+'Y']*refl[4]**2
    16161616        Ssig,Sgam = GetSampleSigGam(refl,0.0,G,GB,hfx,phfx,calcControls,parmDict)
    … …  
    18841884                    hfx+'alpha':[1./refl[4],'alp'],hfx+'beta-0':[1.0,'bet'],hfx+'beta-1':[1./refl[4]**4,'bet'],
    18851885                    hfx+'beta-q':[1./refl[4],'bet'],hfx+'sig-0':[1.0,'sig'],hfx+'sig-1':[refl[4]**2,'sig'],
    1886                     hfx+'sig-q':[refl[4],'sig'],hfx+'Absorption':[dFdAb,'int'],hfx+'Extinction':[dFdEx,'int'],}
     1886                    hfx+'sig-2':[refl[4]**4,'sig'],hfx+'sig-q':[refl[4],'sig'],
     1887                    hfx+'Absorption':[dFdAb,'int'],hfx+'Extinction':[dFdEx,'int'],}
    18871888            for name in names:
    18881889                item = names[name]
Note: See TracChangeset for help on using the changeset viewer.