Changeset 267


Ignore:
Timestamp:
Apr 20, 2011 1:09:53 PM (11 years ago)
Author:
vondreele
Message:

removed azm rotation - didn't work
fix azimuthal polarization

Location:
trunk
Files:
7 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASII.py

    r265 r267  
    474474                            Data['outChannels'] = 2500
    475475                            Data['outAzimuths'] = 1
    476                             Data['azmthRotate'] = 0.0
    477476                            Data['fullIntegrate'] = False
    478477                            Data['setRings'] = False
     
    545544            mainSizer.Add(topLabl,0,wx.ALIGN_CENTER_VERTICAL|wx.LEFT,10)
    546545            mainSizer.Add((10,10),1)
    547             dataGridSizer = wx.FlexGridSizer(rows=len(data),cols=1,hgap=2,vgap=2)
     546            dataGridSizer = wx.FlexGridSizer(rows=len(data),cols=2,hgap=2,vgap=2)
    548547            for id,item in enumerate(self.data):
    549548                ckbox = wx.CheckBox(panel,id,item[1])
     
    12751274                            'Container Bkg.':{'Name':'','Mult':-1.0,'Add':0.0},'ElList':{},
    12761275                            'Geometry':'Cylinder','Diam':1.0,'Pack':0.50,'Form Vol':10.0,
    1277                             'DetType':'Image plate','ObliqCoeff':0.2,'Ruland':0.003,'QScaleLim':[0,100],
     1276                            'DetType':'Image plate','ObliqCoeff':0.2,'Ruland':0.025,'QScaleLim':[0,100],
    12781277                            'Lorch':True,}
    12791278                        self.PatternTree.SetItemPyData(self.PatternTree.AppendItem(Id,text='PDF Controls'),Data)
  • trunk/GSASIIIO.py

    r265 r267  
    804804                Id = item
    805805            item, cookie = self.PatternTree.GetNextChild(self.root, cookie)
    806         parms = ['PXC',data['wavelength'],0.0,0.95,1.0,-1.0,0.3,0.0,1.0,0.0,azm]    #set polarization for synchrotron radiation!
     806        parms = ['PXC',data['wavelength'],0.0,0.99,1.0,-1.0,0.3,0.0,1.0,0.0,azm]    #set polarization for synchrotron radiation!
    807807        Y = self.Integrate[0][i]
    808808        W = 1./Y                    #probably not true
  • trunk/GSASIIimage.py

    r265 r267  
    293293    LRazim = data['LRazimuth']
    294294    azmthoff = data['azmthOff']
    295     azmRot = data['azmthRotate']
    296     Full = data['fullIntegrate']
    297295    dx = np.array(x-cent[0],dtype=np.float32)
    298296    dy = np.array(y-cent[1],dtype=np.float32)
     
    303301    dsp = wave/(2.*npsind(tth/2.))
    304302    azm = (npatan2d(dx,-dy)+azmthoff+720.)%360.
    305     if Full:
    306         azm = (azm+azmRot+720.)%360.
    307     else:
    308         azm = np.where(azm<LRazim[0],azm+360.,azm)
     303    azm = np.where(azm<LRazim[0],azm+360.,azm)
    309304    return tth,azm,dsp
    310305   
     
    565560    numAzms = data['outAzimuths']
    566561    numChans = data['outChannels']
    567     azmRot = data['azmthRotate']
    568     Full = data['fullIntegrate']
    569562    Dtth = (LUtth[1]-LUtth[0])/numChans
    570563    Dazm = (LRazm[1]-LRazm[0])/numAzms
     
    611604        if Dazm:       
    612605            H1 = np.array([azm for azm in np.linspace(LRazm[0],LRazm[1],numAzms+1)])
    613             if Full:
    614                 H1 = H1+azmRot
    615606        else:
    616607            H1 = LRazm
  • trunk/GSASIIimgGUI.py

    r265 r267  
    179179            data['fullIntegrate'] = True
    180180            self.Lazim.SetEditable(False)           
    181             self.Razim.SetEditable(False)           
    182         G2plt.PlotExposedImage(self,event=event)
    183        
    184     def OnAzmRot(event):
    185         try:
    186             azmrot = float(azmRot.GetValue())
    187             data['azmthRotate'] = azmrot
    188         except ValueError:
    189             pass
    190         azmRot.SetValue("%.2f" % (data['azmthRotate']))          #reset in case of error         
    191         G2plt.PlotExposedImage(self,event=event)
    192                
     181            self.Razim.SetEditable(False)
     182        UpdateImageControls(self,data,masks)           
     183        G2plt.PlotExposedImage(self,event=event)
     184       
    193185    def OnSetDefault(event):
    194186        import copy
     
    258250                            id = G2gd.GetPatternTreeItemId(self, self.root, name)
    259251                            Npix,imagefile = self.PatternTree.GetItemPyData(id)
    260                             print imagefile
    261252                            image = G2IO.GetImageData(self,imagefile,True)
    262253                            Id = G2gd.GetPatternTreeItemId(self,id, 'Image Controls')
     
    373364    if 'azmthOff' not in data:
    374365        data['azmthOff'] = 0.0
    375     if 'azmthRotate' not in data:
    376         data['azmthRotate'] = 0.0       
    377366    #end fix
    378367   
     
    548537    rotSel.SetBackgroundColour(VERY_LIGHT_GREY)
    549538    dataSizer.Add(rotSel,0,wx.ALIGN_CENTER_VERTICAL)
    550     dataSizer.Add(wx.StaticText(parent=self.dataDisplay,label=' Bin offset'),0,
    551         wx.ALIGN_CENTER_VERTICAL)
    552     azmRot = wx.TextCtrl(parent=self.dataDisplay,value=("%.2f" % (data['azmthRotate'])),style=wx.TE_PROCESS_ENTER)
    553     azmRot.Bind(wx.EVT_TEXT_ENTER,OnAzmRot)
    554     azmRot.Bind(wx.EVT_KILL_FOCUS,OnAzmRot)
    555     dataSizer.Add(azmRot,0,wx.ALIGN_CENTER_VERTICAL)
    556539    setDefault = wx.CheckBox(parent=self.dataDisplay,label='Use as default for all images?')
    557540    dataSizer.Add(setDefault,0,wx.ALIGN_CENTER_VERTICAL)
  • trunk/GSASIIpeak.py

    r265 r267  
    166166#    return (Pola*npcosd(Azm)**2+(1.-Pola)*npsind(Azm)**2)*npcosd(Tth)**2+ \
    167167#        Pola*npsind(Azm)**2+(1.-Pola)*npcosd(Azm)**2
    168 #    return (Pola*npcosd(Azm)**2+npsind(Azm)**2)*npcosd(Tth)**2+   \
    169 #        Pola*npsind(Azm)**2+npcosd(Azm)**2
    170     return Pola*npcosd(Tth)**2+1.0
     168    return ((1.0-Pola)*npcosd(Azm)**2+Pola*npsind(Azm)**2)*npcosd(Tth)**2+   \
     169        (1.0-Pola)*npsind(Azm)**2+Pola*npcosd(Azm)**2
    171170   
    172171def Oblique(ObCoeff,Tth):
     
    672671    dt = (Tth[1]-Tth[0])
    673672    xydata['IofQ'][1][1] /= Absorb(data['Geometry'],Abs,data['Diam'],Tth)
    674     pola = Polarization(inst['Polariz.'],Tth,Azm=inst['Azimuth'])
    675     auxPlot.append([Tth,pola,'Pola'])     
    676     xydata['IofQ'][1][1] /= pola
     673    xydata['IofQ'][1][1] /= Polarization(inst['Polariz.'],Tth,Azm=inst['Azimuth'])
    677674    xydata['IofQ'][1][1] *= Oblique(data['ObliqCoeff'],Tth)
    678675    XY = xydata['IofQ'][1]   
     
    706703    Q = xydata['SofQ'][1][0]
    707704    ruland = Ruland(data['Ruland'],wave,Q,CF)
    708     auxPlot.append([Q,ruland,'Ruland'])     
     705#    auxPlot.append([Q,ruland,'Ruland'])     
    709706    CF *= ruland
    710     auxPlot.append([Q,CF,'Compton*Ruland'])
    711     scale =       np.sum((FFSq+CF)[minQ:maxQ])/np.sum(xydata['SofQ'][1][1][minQ:maxQ])
     707#    auxPlot.append([Q,CF,'Compton*Ruland'])
     708    scale = np.sum((FFSq+CF)[minQ:maxQ])/np.sum(xydata['SofQ'][1][1][minQ:maxQ])
    712709    xydata['SofQ'][1][1] *= scale
    713710    xydata['SofQ'][1][1] -= CF
    714711    xydata['SofQ'][1][1] = xydata['SofQ'][1][1]/SqFF
     712   
    715713
    716714    xydata['FofQ'] = copy.deepcopy(xydata['SofQ'])
  • trunk/GSASIIplot.py

    r265 r267  
    11631163        Page.canvas.SetToolTipString('')
    11641164        sizexy = Data['size']
    1165         azmRot = 0.
    1166         if Data['fullIntegrate']:
    1167             azmRot = Data['azmthRotate']
    11681165        if event.xdata and event.ydata:                 #avoid out of frame errors
    11691166            Page.canvas.SetCursor(wx.CROSS_CURSOR)
     
    11811178                    tth,azm = G2img.GetTthAzm(event.xdata,event.ydata,Data)
    11821179                    if 'line3' in  str(item) or 'line4' in str(item) and not Data['fullIntegrate']:
    1183                         Page.canvas.SetToolTipString('%6d deg'%(azm-azmRot))
     1180                        Page.canvas.SetToolTipString('%6d deg'%(azm))
    11841181                    elif 'line1' in  str(item) or 'line2' in str(item):
    11851182                        Page.canvas.SetToolTipString('%8.3fdeg'%(tth))                           
     
    11981195                else:
    11991196                    self.G2plotNB.status.SetFields(\
    1200                         ['','Detector 2-th =%9.3fdeg, dsp =%9.3fA, Q = %6.5fA-1, azm = %7.2fdeg, I = %6d'%(tth,dsp,Q,azm-azmRot,Int)])
     1197                        ['','Detector 2-th =%9.3fdeg, dsp =%9.3fA, Q = %6.5fA-1, azm = %7.2fdeg, I = %6d'%(tth,dsp,Q,azm,Int)])
    12011198
    12021199    def OnImPlotKeyPress(event):
     
    14941491            if Data['fullIntegrate']:
    14951492                Azm = np.array(range(0,361))
    1496                 AzmRot = Data['azmthRotate']
    14971493            else:
    14981494                Azm = np.array(range(LRAzim[0],LRAzim[1]+1))-AzmthOff
    1499                 AzmRot = 0.0
    15001495            if ellI:
    15011496                xyI = []
     
    15171512            if Nazm > 1:
    15181513                for i in range(Nazm):
    1519                     cake = (LRAzim[0]+i*delAzm+AzmRot+720)%360
     1514                    if Data['fullIntegrate']:
     1515                        cake = (LRAzim[0]+i*delAzm+720)%360
     1516                    else:
     1517                        cake = LRAzim[0]+i*delAzm
    15201518                    ind = np.searchsorted(Azm,cake)
    15211519                    Plot.plot([arcxI[ind],arcxO[ind]],[arcyI[ind],arcyO[ind]],color='k',dashes=(5,5))
  • trunk/ImageCalibrants.py

    r251 r267  
    1515'NaCl even hkl only':(2,(2.8201,2.8201,2.8201,90,90,90),0,(0.5,20,10)),
    1616'Ag behenate':(6,(1.0,1.0,58.380,90,90,90),0,(7.0,5,1)),
    17 'Spun Si grating':(6,(1.0,1.0,2777.78,90,90,90),2,(200.,5,1))
     17'Spun Si grating':(6,(1.0,1.0,2777.78,90,90,90),2,(200.,5,1)),
     18'Pt   @ 298K':(0,(3.9231,3.9231,3.9231,90,90,90),0,(0.5,5,1)),
    1819    }
Note: See TracChangeset for help on using the changeset viewer.