Changeset 797


Ignore:
Timestamp:
Nov 8, 2012 4:05:35 PM (10 years ago)
Author:
vondreele
Message:

some fixes to TOF stuff
some fixes to SXC stuff

Location:
trunk
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASII.py

    r796 r797  
    874874                    Ibeg = np.searchsorted(rd.powderdata[0],Iparm2['Tminmax'][0])
    875875                    Ifin = np.searchsorted(rd.powderdata[0],Iparm2['Tminmax'][1])
    876                     print Ibeg,Ifin,Iparm2['Tminmax']
    877876                    rd.powderdata[0] = rd.powderdata[0][Ibeg:Ifin]
    878877                    YI,WYI = G2pwd.calcIncident(Iparm2,rd.powderdata[0])
    879878                    rd.powderdata[1] = rd.powderdata[1][Ibeg:Ifin]/YI
    880                     rd.powderdata[2] = rd.powderdata[2][Ibeg:Ifin]+(rd.powderdata[1]**2+WYI)
    881                     rd.powderdata[2] /= YI**2
     879                    var = 1./rd.powderdata[2][Ibeg:Ifin]
     880                    var += rd.powderdata[1]**2+WYI
     881                    var /= YI**2
     882                    rd.powderdata[2] = 1./var
    882883                rd.powderdata[3] = np.zeros_like(rd.powderdata[0])                                       
    883884                rd.powderdata[4] = np.zeros_like(rd.powderdata[0])                                       
  • trunk/GSASIIIO.py

    r796 r797  
    12671267        if Wave is not None:
    12681268            HistWave = Wave
    1269         self.Parameters = [(HistType,HistWave),[HistType,HistWave],  # overwrite entire list
    1270             ['Type','Lam']]
     1269        self.Parameters = [{'Type':[HistType,HistType],'Lam':[HistWave,HistWave]},{}]  # overwrite entire list
    12711270           
    12721271    def UpdateControls(self,Type='Fosq',FcalcPresent=False):
  • trunk/GSASIImath.py

    r796 r797  
    10021002            for x in ['alpha','beta-0','beta-0','var-inst','X','Y']:
    10031003                ins[x] = Parms[x][1]
    1004             alp = ins['alpha']*dsp
     1004            alp = ins['alpha']/dsp
    10051005            bet = ins['beta-0']+ins['beta-0']/dsp**4
    10061006        sig = ins['var-inst']*dsp**2
  • trunk/GSASIIphsGUI.py

    r796 r797  
    919919            kFactor.Bind(wx.EVT_KILL_FOCUS,OnkFactor)
    920920            line2Sizer.Add(kFactor,0,wx.ALIGN_CENTER_VERTICAL)
    921             line2Sizer.Add(wx.StaticText(dataDisplay,label=' k-Max (<10.0): '),0,wx.ALIGN_CENTER_VERTICAL)
     921            line2Sizer.Add(wx.StaticText(dataDisplay,label=' k-Max (>=10.0): '),0,wx.ALIGN_CENTER_VERTICAL)
    922922            kMax = wx.TextCtrl(dataDisplay,value='%.0f'%(Flip['k-Max']),style=wx.TE_PROCESS_ENTER)
    923923            kMax.Bind(wx.EVT_TEXT_ENTER,OnkMax)       
  • trunk/GSASIIplot.py

    r796 r797  
    192192            dlg.Destroy()
    193193   
    194 
    195194################################################################################
    196195##### PlotSngl
     
    12241223    Seen when "Instrument Parameters" chosen from powder pattern data tree
    12251224    '''
     1225    sig = lambda Th,U,V,W: 1.17741*math.sqrt(U*tand(Th)**2+V*tand(Th)+W)*math.pi/18000.
     1226    gam = lambda Th,X,Y: (X/cosd(Th)+Y*tand(Th))*math.pi/18000.
     1227    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.)
     1228#    gamFW2 = lambda s,g: math.sqrt(s**2+(0.4654996*g)**2)+.5345004*g  #Ubaldo Bafile - private communication
    12261229    PatternId = G2frame.PatternId
    12271230    limitID = G2gd.GetPatternTreeItemId(G2frame,PatternId, 'Limits')
     
    12441247    else:
    12451248        difC = Parms['difC'][0]
     1249        difA = Parms['difA'][0]
     1250        Zero = Parms['Zero'][0]
    12461251        alp = Parms['alpha'][0]
    12471252        bet0 = Parms['beta-0'][0]
     
    12701275    Page.canvas.SetToolTipString('')
    12711276    colors=['b','g','r','c','m','k']
    1272     Xmin,Xmax = limits[1]
    1273     Xmin = min(0.5,max(Xmin,1))
    1274     Xmin /= 2
    1275     Xmax /= 2
    1276     nPts = 100
    1277     delt = (Xmax-Xmin)/nPts
    1278     thetas = []
    1279     for i in range(nPts):
    1280         thetas.append(Xmin+i*delt)
    12811277    X = []
    12821278    Y = []
    12831279    Z = []
    12841280    W = []
    1285     sig = lambda Th,U,V,W: 1.17741*math.sqrt(U*tand(Th)**2+V*tand(Th)+W)*math.pi/18000.
    1286     gam = lambda Th,X,Y: (X/cosd(Th)+Y*tand(Th))*math.pi/18000.
    1287     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.)
    1288 #    gamFW2 = lambda s,g: math.sqrt(s**2+(0.4654996*g)**2)+.5345004*g  #Ubaldo Bafile - private communication
    1289     Plot.set_title('Instrument and sample peak widths')
    1290     Plot.set_ylabel(r'$\Delta q/q, \Delta d/d$',fontsize=14)
    1291     Plot.set_xlabel(r'$q, \AA^{-1}$',fontsize=14)
    12921281    if 'C' in Parms['Type'][0]:
     1282        Plot.set_title('Instrument and sample peak widths')
     1283        Plot.set_xlabel(r'$q, \AA^{-1}$',fontsize=14)
     1284        Plot.set_ylabel(r'$\Delta q/q, \Delta d/d$',fontsize=14)
    12931285        try:
     1286            Xmin,Xmax = limits[1]
     1287            Xmin = min(0.5,max(Xmin,1))
     1288            Xmin /= 2
     1289            Xmax /= 2
     1290            nPts = 100
     1291            delt = (Xmax-Xmin)/nPts
     1292            thetas = []
     1293            for i in range(nPts):
     1294                thetas.append(Xmin+i*delt)
    12941295            for theta in thetas:
    12951296                X.append(4.0*math.pi*sind(theta)/lam)              #q
     
    13291330            G2frame.G2plotNB.Delete('Peak Widths')
    13301331    else:
    1331         pass    #for TOF peak parms
     1332        Plot.set_title('Instrument and sample peak coefficients')
     1333        Plot.set_xlabel(r'$TOF, \mu s$',fontsize=14)
     1334        Plot.set_ylabel(r'$\alpha, \beta, \Delta T$',fontsize=14)
     1335        Xmin,Xmax = limits[1]
     1336        if 'Pabc' in Parms2:
     1337            Pabc = Parms2['Pabc']
     1338            print Pabc.shape,Pabc[0]
     1339           
     1340        else:           
     1341            T = np.linspace(Xmin,Xmax,num=101,endpoint=True)
     1342            ds = T/difC
     1343            A = alp/ds
     1344            B = bet0+bet1/ds**4
     1345            D = difA*ds**2+Zero
     1346            Plot.plot(T,A,color='r',label='Alpha')
     1347            Plot.plot(T,B,color='g',label='Beta')
     1348            Plot.plot(T,D,color='b',label='Delta-T')
     1349        T = []
     1350        A = []
     1351        B = []
     1352        W = []
     1353        V = []
     1354        for peak in peaks:
     1355            T.append(peak[0])
     1356            A.append(peak[4])
     1357            B.append(peak[6])
     1358            W.append(peak[0]/difC)
     1359           
     1360       
     1361        Plot.plot(T,A,'+',color='r',label='Alpha peak')
     1362        Plot.plot(T,B,'+',color='g',label='Beta peak')
     1363        Plot.plot(T,W,'+',color='k',label='T/difC')
     1364        Plot.legend(loc='best')
     1365        Page.canvas.draw()
    13321366
    13331367   
  • trunk/GSASIIpwd.py

    r796 r797  
    931931                    alp = parmDict[alpName]
    932932                else:
    933                     alp = alp0*dsp
     933                    alp = alp0/dsp
    934934                betName = 'bet'+str(iPeak)
    935935                if betName in varyList:
     
    10981098                    alp = parmDict[alpName]
    10991099                else:
    1100                     alp = alp0*dsp
    1101                     dada0 = dsp
     1100                    alp = alp0/dsp
     1101                    dada0 = 1./dsp
    11021102                betName = 'bet'+str(iPeak)
    11031103                if betName in varyList:
     
    13891389                    peak[2*j] = parmDict[parName]
    13901390                elif 'alpha' in parName:
    1391                     peak[2*j] = parmDict['alpha']*dsp
     1391                    peak[2*j] = parmDict['alpha']/dsp
    13921392                elif 'beta' in parName:
    13931393                    peak[2*j] = parmDict['beta-0']+parmDict['beta-1']/dsp**4
Note: See TracChangeset for help on using the changeset viewer.