Changeset 1462 for trunk/GSASIIpwd.py


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

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • 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])
Note: See TracChangeset for help on using the changeset viewer.