Changeset 1412


Ignore:
Timestamp:
Jul 6, 2014 9:21:55 AM (7 years ago)
Author:
vondreele
Message:

fix theta vs 2-theta error in peak fitting; gave wrong sig,gam vs U,V,W,X,Y now all match OK

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIpwd.py

    r1411 r1412  
    791791            try:
    792792                pos = parmDict['pos'+str(iPeak)]
    793                 theta = (pos-parmDict['Zero'])/2.0
     793                tth = (pos-parmDict['Zero'])
    794794                intens = parmDict['int'+str(iPeak)]
    795795                sigName = 'sig'+str(iPeak)
     
    797797                    sig = parmDict[sigName]
    798798                else:
    799                     sig = G2mth.getCWsig(parmDict,theta)
     799                    sig = G2mth.getCWsig(parmDict,tth)
    800800                sig = max(sig,0.001)          #avoid neg sigma
    801801                gamName = 'gam'+str(iPeak)
     
    803803                    gam = parmDict[gamName]
    804804                else:
    805                     gam = G2mth.getCWgam(parmDict,theta)
     805                    gam = G2mth.getCWgam(parmDict,tth)
    806806                gam = max(gam,0.001)             #avoid neg gamma
    807807                Wd,fmin,fmax = getWidthsCW(pos,sig,gam,shl)
     
    912912            try:
    913913                pos = parmDict['pos'+str(iPeak)]
    914                 theta = (pos-parmDict['Zero'])/2.0
     914                tth = (pos-parmDict['Zero'])
    915915                intens = parmDict['int'+str(iPeak)]
    916916                sigName = 'sig'+str(iPeak)
    917                 tanth = tand(theta)
    918                 costh = cosd(theta)
    919917                if sigName in varyList:
    920918                    sig = parmDict[sigName]
    921919                    dsdU = dsdV = dsdW = 0
    922920                else:
    923                     sig = G2mth.getCWsig(parmDict,theta)
    924                     dsdU,dsdV,dsdW = G2mth.getCWsigDeriv(theta)
     921                    sig = G2mth.getCWsig(parmDict,tth)
     922                    dsdU,dsdV,dsdW = G2mth.getCWsigDeriv(tth)
    925923                sig = max(sig,0.001)          #avoid neg sigma
    926924                gamName = 'gam'+str(iPeak)
     
    929927                    dgdX = dgdY = 0
    930928                else:
    931                     gam = G2mth.getCWgam(parmDict,theta)
    932                     dgdX,dgdY = G2mth.getCWgamDeriv(theta)
     929                    gam = G2mth.getCWgam(parmDict,tth)
     930                    dgdX,dgdY = G2mth.getCWgamDeriv(tth)
    933931                gam = max(gam,0.001)             #avoid neg gamma
    934932                Wd,fmin,fmax = getWidthsCW(pos,sig,gam,shl)
Note: See TracChangeset for help on using the changeset viewer.