Changeset 998


Ignore:
Timestamp:
Jul 18, 2013 8:53:20 AM (9 years ago)
Author:
vondreele
Message:

make widths correct in Export HKLs from PWD data

Location:
trunk
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASII.py

    r997 r998  
    21602160                                    file.write('%s \n'%(' h  k  l  m  2-theta wid F**2'))               
    21612161                                    for peak in peaks:
    2162                                         FWHM = G2pwd.getgamFW(peak[7],peak[6])
     2162                                        FWHM = G2pwd.getgamFW(peak[7],peak[6])/50.      #to get delta-2-theta in deg.
    21632163                                        file.write(" %3d %3d %3d %3d %10.5f %10.5f %10.3f \n" % \
    21642164                                            (int(peak[0]),int(peak[1]),int(peak[2]),int(peak[3]),peak[5],FWHM,peak[8]))
  • trunk/GSASIIplot.py

    r993 r998  
    13151315            s = 1.17741*np.sqrt(data[4])*np.pi/18000.
    13161316            g = data[6]*np.pi/18000.
    1317             G = gamFW(g,s)
     1317            G = G2pwd.getgamFW(g,s)
    13181318            Y = s/nptand(X/2.)
    13191319            Z = g/nptand(X/2.)
     
    13261326            sf = 1.17741*np.sqrt(fit[4])*np.pi/18000.
    13271327            gf = fit[6]*np.pi/18000.
    1328             Gf = gamFW(gf,sf)
     1328            Gf = G2pwd.getgamFW(gf,sf)
    13291329            Yf = sf/nptand(X/2.)
    13301330            Zf = gf/nptand(X/2.)
     
    13461346                    s = 0.01
    13471347                g = peak[6]*math.pi/18000.
    1348                 G = gamFW(g,s)
     1348                G = G2pwd.getgamFW(g,s)
    13491349                Y.append(s/tand(peak[0]/2.))
    13501350                Z.append(g/tand(peak[0]/2.))
  • trunk/GSASIIpwd.py

    r997 r998  
    449449def getgamFW(g,s):
    450450    'needs a doc string'
    451     gamFW = lambda s,g: math.exp(math.log(s**5+2.69269*s**4*g+2.42843*s**3*g**2+4.47163*s**2*g**3+0.07842*s*g**4+g**5)/5.)
     451    gamFW = lambda s,g: np.exp(np.log(s**5+2.69269*s**4*g+2.42843*s**3*g**2+4.47163*s**2*g**3+0.07842*s*g**4+g**5)/5.)
    452452    return gamFW(g,s)
    453453               
Note: See TracChangeset for help on using the changeset viewer.