Changeset 1203


Ignore:
Timestamp:
Jan 20, 2014 1:07:50 PM (9 years ago)
Author:
vondreele
Message:

show q in powder/sasd plots
work on absorption corr. for sasd data

Location:
trunk
Files:
7 edited

Legend:

Unmodified
Added
Removed

  • trunk/GSASIIIO.py

    r1199 r1203  
    858858                Id = item
    859859            item, cookie = G2frame.PatternTree.GetNextChild(G2frame.root, cookie)
    860         if 'PWDR' in name:
    861             parms = ['PXC',data['wavelength'],0.0,0.99,1.0,-0.10,0.4,0.30,1.0,0.0001,Azms[i]]    #set polarization for synchrotron radiation!
    862         elif 'SASD' in name:   
    863             parms = ['LXC',data['wavelength'],0.0,Azms[i]]           
    864         Y = G2frame.Integrate[0][i]
    865         W = 1./Y                    #probably not true
    866860        Sample = G2pdG.SetDefaultSample()
    867861        Sample['Gonio. radius'] = data['distance']
    … …  
    869863        Sample['Chi'] = data['GonioAngles'][1]
    870864        Sample['Phi'] = data['GonioAngles'][2]
     865        if 'PWDR' in name:
     866            parms = ['PXC',data['wavelength'],0.0,0.99,1.0,-0.10,0.4,0.30,1.0,0.0001,Azms[i]]    #set polarization for synchrotron radiation!
     867        elif 'SASD' in name:   
     868            parms = ['LXC',data['wavelength'],0.0,Azms[i]]
     869        Y = G2frame.Integrate[0][i]
     870        W = 1./Y                    #probably not true
    871871        if Id:
    872872            G2frame.PatternTree.SetItemPyData(G2gd.GetPatternTreeItemId(G2frame,Id, 'Comments'),Comments)                   

  • trunk/GSASIIgrid.py

    r1202 r1203  
    23032303        self.PrefillDataMenu(self.SampleMenu,helpType='Sample Parameters')
    23042304        self.SampleEdit = wx.Menu(title='')
    2305         self.SampleMenu.Append(menu=self.SampleEdit, title='File')
     2305        self.SampleMenu.Append(menu=self.SampleEdit, title='Command')
    23062306        self.SampleEdit.Append(id=wxID_SAMPLELOAD, kind=wx.ITEM_NORMAL,text='Load',
    23072307            help='Load sample parameters from file')

  • trunk/GSASIIimage.py

    r1199 r1203  
    797797    muT = data['SampleAbs'][0]
    798798    if 'SASD' in data['type']:
    799         muT = np.log(muT)
     799        muT = -np.log(muT)/2.       #Transmission to 1/2 thickness muT
    800800    NST = np.zeros(shape=(numAzms,numChans),order='F',dtype=np.float32)
    801801    H0 = np.zeros(shape=(numAzms,numChans),order='F',dtype=np.float32)

  • trunk/GSASIIimgGUI.py

    r1199 r1203  
    278278            data['type'] = typeSel.GetValue()[:4]
    279279            if 'SASD' in data['type'] and not data['SampleAbs'][0]:
    280                 data['SampleAbs'][0] = 1.0  #switch from muT to trans!
     280                data['SampleAbs'][0] = 1.0  #switch from muT=0 to trans=1!
    281281            wx.CallAfter(UpdateImageControls,G2frame,data,masks)
    282282   

  • trunk/GSASIIphsGUI.py

    r1183 r1203  
    48584858                    ref[6] = FWHM*(xdata[1][indx]-xdata[4][indx])*cosd(pos/2.)**3/dx
    48594859                    Lorenz = 1./(2.*sind(xdata[0][indx]/2.)**2*cosd(xdata[0][indx]/2.))           #Lorentz correction
    4860                     pola,dIdPola = G2pwd.Polarization(Inst['Polariz.'][1],xdata[0][indx],Inst['Azimuth'][1])
     4860                    pola = 1.0
     4861                    if 'X' in Inst['Type']:
     4862                        pola,dIdPola = G2pwd.Polarization(Inst['Polariz.'][1],xdata[0][indx],Inst['Azimuth'][1])
    48614863                    ref[6] /= (Lorenz*pola*ref[3])
    48624864                except IndexError:

  • trunk/GSASIIplot.py

    r1201 r1203  
    508508                    if abs(xpos) > 0.:                  #avoid possible singularity at beam center
    509509                        dsp = wave/(2.*sind(abs(xpos)/2.0))
     510                        q = 2.*np.pi/dsp
    510511                    if G2frame.Contour:
    511                         G2frame.G2plotNB.status.SetStatusText('2-theta =%9.3f d =%9.5f pattern ID =%5d'%(xpos,dsp,int(ypos)),1)
     512                        G2frame.G2plotNB.status.SetStatusText('2-theta =%9.3f d =%9.5f q = %9.5f pattern ID =%5d'%(xpos,dsp,q,int(ypos)),1)
    512513                    else:
    513                         G2frame.G2plotNB.status.SetStatusText('2-theta =%9.3f d =%9.5f Intensity =%9.1f'%(xpos,dsp,ypos),1)
     514                        G2frame.G2plotNB.status.SetStatusText('2-theta =%9.3f d =%9.5f q = %9.5f Intensity =%9.2f'%(xpos,dsp,q,ypos),1)
    514515                else:       #TOF neutrons
    515516                    dsp = 0.0
    516517                    difC = Parms['difC'][1]
    517518                    dsp = xpos/difC             #rough approx.!
     519                    q = 2.*np.pi/dsp
    518520                    if G2frame.Contour:
    519                         G2frame.G2plotNB.status.SetStatusText('TOF =%9.3f d =%9.5f pattern ID =%5d'%(xpos,dsp,int(ypos)),1)
     521                        G2frame.G2plotNB.status.SetStatusText('TOF =%9.3f d =%9.5f q = %9.5f pattern ID =%5d'%(xpos,dsp,q,int(ypos)),1)
    520522                    else:
    521                         G2frame.G2plotNB.status.SetStatusText('TOF =%9.3f d =%9.5f Intensity =%9.1f'%(xpos,dsp,ypos),1)
     523                        G2frame.G2plotNB.status.SetStatusText('TOF =%9.3f d =%9.5f q = %9.5f Intensity =%9.2f'%(xpos,dsp,q,ypos),1)
    522524                if G2frame.itemPicked:
    523525                    Page.canvas.SetToolTipString('%9.3f'%(xpos))
    … …  
    706708        Sample = G2frame.PatternTree.GetItemPyData(G2gd.GetPatternTreeItemId(G2frame,G2frame.PatternId, 'Sample Parameters'))
    707709        ParmList = [Parms,]
     710        SampleList = [Sample,]
    708711        Title = Pattern[-1]
    709712    else:       
    … …  
    751754    Plot.set_title(Title)
    752755    if G2frame.qPlot:
    753         Plot.set_xlabel(r'$Q, \AA^{-1}$',fontsize=14)
     756        Plot.set_xlabel(r'$Q, \AA^{-1}$',fontsize=16)
    754757    else:
    755758        if 'C' in ParmList[0]['Type'][0]:       
    756             Plot.set_xlabel(r'$\mathsf{2\theta}$',fontsize=14)
     759            Plot.set_xlabel(r'$\mathsf{2\theta}$',fontsize=16)
    757760        else:
    758             Plot.set_xlabel(r'TOF, $\mathsf{\mu}$s',fontsize=14)           
     761            Plot.set_xlabel(r'$TOF, \mathsf{\mu}$s',fontsize=16)           
    759762    if G2frame.Weight:
    760         Plot.set_ylabel(r'$\mathsf{I/\sigma(I)}$',fontsize=14)
     763        Plot.set_ylabel(r'$\mathsf{I/\sigma(I)}$',fontsize=16)
    761764    else:
    762765        if 'C' in ParmList[0]['Type'][0]:
    763766            if 'PWDR' in plottype:
    764                 Plot.set_ylabel('Intensity',fontsize=14)
     767                Plot.set_ylabel('$Intensity$',fontsize=16)
    765768            elif 'SASD' in plottype:
    766                 Plot.set_ylabel('Intensity, cm-1',fontsize=14)
     769                Plot.set_ylabel('$Intensity, cm^{-1}$',fontsize=16)
    767770        else:
    768             Plot.set_ylabel('Normalized intensity',fontsize=14)
     771            Plot.set_ylabel('Normalized intensity',fontsize=16)
    769772    if G2frame.Contour:
    770773        ContourZ = []

  • trunk/GSASIIpwdGUI.py

    r1201 r1203  
    6464    'Needs a doc string'
    6565    return {
     66        'InstrName':'',
    6667        'ranId':ran.randint(0,sys.maxint),
    6768        'Scale':[1.0,True],'Type':'Debye-Scherrer','Absorption':[0.0,False],
    … …  
    11421143        finally:
    11431144            dlg.Destroy()
    1144                                                
    1145        
     1145                                                       
    11461146    def OnSampleLoad(event):
    11471147        '''Loads sample parameters from a G2 .samprm file
Note: See TracChangeset for help on using the changeset viewer.