Changeset 2132 for trunk/GSASIIplot.py


Ignore:
Timestamp:
Jan 22, 2016 10:10:43 AM (6 years ago)
Author:
vondreele
Message:

more rationalization of FWHM calculations

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIplot.py

    r2127 r2132  
    26632663            Z = np.ones_like(X)
    26642664            data = G2mth.setPeakparms(Parms,Parms2,X,Z)
    2665             s = 1.17741*np.sqrt(data[4])*np.pi/18000.
    2666             g = data[6]*np.pi/36000.
    2667             G = G2pwd.getgamFW(g,s)
     2665            s = np.sqrt(data[4])*np.pi/18000.   #var -> sig(radians)
     2666            g = data[6]*np.pi/18000.    #centideg -> radians
     2667            G = G2pwd.getgamFW(g,s)/2.  #delt-theta
    26682668            Y = s/nptand(X/2.)
    26692669            Z = g/nptand(X/2.)
     
    26742674           
    26752675            fit = G2mth.setPeakparms(Parms,Parms2,X,Z,useFit=True)
    2676             sf = 1.17741*np.sqrt(fit[4])*np.pi/18000.
    2677             gf = fit[6]*np.pi/36000.
    2678             Gf = G2pwd.getgamFW(gf,sf)
     2676            sf = np.sqrt(fit[4])*np.pi/18000.
     2677            gf = fit[6]*np.pi/18000.
     2678            Gf = G2pwd.getgamFW(gf,sf)/2.
    26792679            Yf = sf/nptand(X/2.)
    26802680            Zf = gf/nptand(X/2.)
     
    26922692                X.append(4.0*math.pi*sind(peak[0]/2.0)/lam)
    26932693                try:
    2694                     s = 1.17741*math.sqrt(peak[4])*math.pi/18000.
     2694                    s = math.sqrt(peak[4])*math.pi/18000.
    26952695                except ValueError:
    26962696                    s = 0.01
    2697                 g = peak[6]*math.pi/36000.
    2698                 G = G2pwd.getgamFW(g,s)
     2697                g = peak[6]*math.pi/18000.
     2698                G = G2pwd.getgamFW(g,s)/2.
    26992699                Y.append(s/tand(peak[0]/2.))
    27002700                Z.append(g/tand(peak[0]/2.))
Note: See TracChangeset for help on using the changeset viewer.