Changeset 728


Ignore:
Timestamp:
Aug 22, 2012 1:09:12 PM (10 years ago)
Author:
vondreele
Message:

trap sqrt error from bad U,V,W

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIplot.py

    r727 r728  
    12301230    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.)
    12311231#    gamFW2 = lambda s,g: math.sqrt(s**2+(0.4654996*g)**2)+.5345004*g  #Ubaldo Bafile - private communication
    1232     for theta in thetas:
    1233         X.append(4.0*math.pi*sind(theta)/lam)              #q
    1234         s = sig(theta,GU,GV,GW)
    1235         g = gam(theta,LX,LY)
    1236         G = gamFW(g,s)
    1237         Y.append(s/tand(theta))
    1238         Z.append(g/tand(theta))
    1239         W.append(G/tand(theta))
    1240     Plot.set_title('Instrument and sample peak widths')
    1241     Plot.set_ylabel(r'$\Delta q/q, \Delta d/d$',fontsize=14)
    1242     Plot.set_xlabel(r'$q, \AA^{-1}$',fontsize=14)
    1243     Plot.plot(X,Y,color='r',label='Gaussian')
    1244     Plot.plot(X,Z,color='g',label='Lorentzian')
    1245     Plot.plot(X,W,color='b',label='G+L')
    1246     X = []
    1247     Y = []
    1248     Z = []
    1249     W = []
    1250     V = []
    1251     for peak in peaks:
    1252         X.append(4.0*math.pi*sind(peak[0]/2.0)/lam)
    1253         try:
    1254             s = 1.17741*math.sqrt(peak[4])*math.pi/18000.
    1255         except ValueError:
    1256             s = 0.01
    1257         g = peak[6]*math.pi/18000.
    1258         G = gamFW(g,s)
    1259         Y.append(s/tand(peak[0]/2.))
    1260         Z.append(g/tand(peak[0]/2.))
    1261         W.append(G/tand(peak[0]/2.))
    1262     Plot.plot(X,Y,'+',color='r',label='G peak')
    1263     Plot.plot(X,Z,'+',color='g',label='L peak')
    1264     Plot.plot(X,W,'+',color='b',label='G+L peak')
    1265     Plot.legend(loc='best')
    1266     Page.canvas.draw()
     1232    try:
     1233        for theta in thetas:
     1234            X.append(4.0*math.pi*sind(theta)/lam)              #q
     1235            s = sig(theta,GU,GV,GW)
     1236            g = gam(theta,LX,LY)
     1237            G = gamFW(g,s)
     1238            Y.append(s/tand(theta))
     1239            Z.append(g/tand(theta))
     1240            W.append(G/tand(theta))
     1241        Plot.set_title('Instrument and sample peak widths')
     1242        Plot.set_ylabel(r'$\Delta q/q, \Delta d/d$',fontsize=14)
     1243        Plot.set_xlabel(r'$q, \AA^{-1}$',fontsize=14)
     1244        Plot.plot(X,Y,color='r',label='Gaussian')
     1245        Plot.plot(X,Z,color='g',label='Lorentzian')
     1246        Plot.plot(X,W,color='b',label='G+L')
     1247        X = []
     1248        Y = []
     1249        Z = []
     1250        W = []
     1251        V = []
     1252        for peak in peaks:
     1253            X.append(4.0*math.pi*sind(peak[0]/2.0)/lam)
     1254            try:
     1255                s = 1.17741*math.sqrt(peak[4])*math.pi/18000.
     1256            except ValueError:
     1257                s = 0.01
     1258            g = peak[6]*math.pi/18000.
     1259            G = gamFW(g,s)
     1260            Y.append(s/tand(peak[0]/2.))
     1261            Z.append(g/tand(peak[0]/2.))
     1262            W.append(G/tand(peak[0]/2.))
     1263        Plot.plot(X,Y,'+',color='r',label='G peak')
     1264        Plot.plot(X,Z,'+',color='g',label='L peak')
     1265        Plot.plot(X,W,'+',color='b',label='G+L peak')
     1266        Plot.legend(loc='best')
     1267        Page.canvas.draw()
     1268    except ValueError:
     1269        print '**** ERROR - default U,V,W profile coefficients yield sqrt of negative value at 2theta =', \
     1270            '%.3f'%(2*theta)
     1271        G2frame.G2plotNB.Delete('Peak Widths')
     1272
    12671273   
    12681274################################################################################
Note: See TracChangeset for help on using the changeset viewer.