Changeset 1791


Ignore:
Timestamp:
Apr 16, 2015 10:17:53 AM (7 years ago)
Author:
vondreele
Message:

use sample angles in sp.harm. PO refinement for non Bragg-Brentano expts.
export powder sig, gam & FWHM in deg - several places

Location:
trunk
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASII.py

    r1785 r1791  
    30673067                                    if len(peak) == 8:  #CW
    30683068                                        FWHM = G2pwd.getgamFW(peak[6],peak[4])
    3069                                         file.write("%10.5f %12.2f %10.3f %10.3f %10.3f \n" % \
    3070                                             (peak[0],peak[2],peak[4],peak[6],FWHM))
     3069                                        file.write("%10.5f %12.2f %10.5f %10.5f %10.5f \n" % \
     3070                                            (peak[0],peak[2],np.sqrt(max(0.0001,peak[4]))/100.,peak[6]/100.,FWHM/100.)) #convert to deg
    30713071                                    else:               #TOF - more cols
    30723072                                        FWHM = G2pwd.getgamFW(peak[10],peak[8])
  • trunk/GSASIIstrMath.py

    r1790 r1791  
    12681268        IFCoup = True
    12691269    else:
    1270         Gangls = [0.,0.,0.,parmDict[hfx+'Azimuth']]
     1270        Gangls = [parmDict[hfx+'Phi'],parmDict[hfx+'Chi'],parmDict[hfx+'Omega'],parmDict[hfx+'Azimuth']]
    12711271        IFCoup = False
    12721272    phi,beta = G2lat.CrsAng(H,cell,SGData)
     
    12941294        IFCoup = True
    12951295    else:
    1296         Gangls = [0.,0.,0.,parmDict[hfx+'Azimuth']]
     1296        Gangls = [parmDict[hfx+'Phi'],parmDict[hfx+'Chi'],parmDict[hfx+'Omega'],parmDict[hfx+'Azimuth']]
    12971297        IFCoup = False
    12981298    phi,beta = G2lat.CrsAng(H,cell,SGData)
  • trunk/exports/G2export_csv.py

    r1773 r1791  
    250250            if phasDict.get('Super',False):
    251251                WriteList(self,("h","k","l","m",tname,"F_obs","F_calc","phase","mult","sig","gam","FWHM","phase #"))
    252                 fmt = "{:.0f},{:.0f},{:.0f},{:.0f},{:.3f},{:.3f},{:.3f},{:.2f},{:.0f},{:.3f},{:.3f},{:.3f},{:d}"
     252                if 'T' in phasDict['Type']:
     253                    fmt = "{:.0f},{:.0f},{:.0f},{:.0f},{:.3f},{:.3f},{:.3f},{:.2f},{:.0f},{:.3f},{:.3f},{:.3f},{:d}"
     254                else:
     255                    fmt = "{:.0f},{:.0f},{:.0f},{:.0f},{:.3f},{:.3f},{:.3f},{:.2f},{:.0f},{:.5f},{:.5f},{:.5f},{:d}"
    253256                refList = phasDict['RefList']
    254257                for refItem in refList:
    255258                    h,k,l,m,mult,dsp,pos,sig,gam,Fobs,Fcalc,phase,Icorr = refItem[:13]
    256259                    FWHM = G2pwd.getgamFW(gam,sig)
    257                     self.Write(fmt.format(h,k,l,m,pos,Fobs,Fcalc,phase,mult,sig,gam,FWHM,i))               
     260                    if 'T' in phasDict['Type']:
     261                        self.Write(fmt.format(h,k,l,m,pos,Fobs,Fcalc,phase,mult,sig,gam,FWHM,i))
     262                    else:        #convert to deg       
     263                        self.Write(fmt.format(h,k,l,m,pos,Fobs,Fcalc,phase,mult,    \
     264                            np.sqrt(max(sig,0.0001))/100.,gam/100.,FWHM/100.,i))
    258265            else:
    259266                WriteList(self,("h","k","l",tname,"F_obs","F_calc","phase","mult","sig","gam","FWHM","phase #"))
    260                 fmt = "{:.0f},{:.0f},{:.0f},{:.3f},{:.3f},{:.3f},{:.2f},{:.0f},{:.3f},{:.3f},{:.3f},{:d}"
     267                if 'T' in phasDict['Type']:
     268                    fmt = "{:.0f},{:.0f},{:.0f},{:.3f},{:.3f},{:.3f},{:.2f},{:.0f},{:.3f},{:.3f},{:.3f},{:d}"
     269                else:
     270                    fmt = "{:.0f},{:.0f},{:.0f},{:.3f},{:.3f},{:.3f},{:.2f},{:.0f},{:.5f},{:.5f},{:.5f},{:d}"
    261271                refList = phasDict['RefList']
    262272                for refItem in refList:
    263273                    h,k,l,mult,dsp,pos,sig,gam,Fobs,Fcalc,phase,Icorr = refItem[:12]
    264274                    FWHM = G2pwd.getgamFW(gam,sig)
    265                     self.Write(fmt.format(h,k,l,pos,Fobs,Fcalc,phase,mult,sig,gam,FWHM,i))
     275                    if 'T' in phasDict['Type']:
     276                        self.Write(fmt.format(h,k,l,pos,Fobs,Fcalc,phase,mult,sig,gam,FWHM,i))
     277                    else:        #convert to deg       
     278                        self.Write(fmt.format(h,k,l,pos,Fobs,Fcalc,phase,mult,  \
     279                            np.sqrt(max(sig,0.0001))/100.,gam/100.,FWHM/100.,i))
    266280        self.CloseFile()
    267281        print(str(hist)+'reflections written to file '+str(self.fullpath))
  • trunk/exports/G2export_examples.py

    r1773 r1791  
    2121'''
    2222import os.path
     23import numpy as np
    2324import GSASIIpath
    2425GSASIIpath.SetVersionNumber("$Revision$")
     
    197198                hklfmt = "{:.0f},{:.0f},{:.0f},{:.0f}"
    198199                hfmt = "{:>10s} {:>8s} {:>12s} {:>12s} {:>7s} {:>6s} {:>8s} {:>8s} {:>8s}"
    199                 fmt = "{:>10s} {:8.3f} {:12.3f} {:12.3f} {:7.2f} {:6.0f} {:8.3f} {:8.3f} {:8.3f}"
     200                if 'T' in phasDict['Type']:
     201                    fmt = "{:>10s} {:8.3f} {:12.3f} {:12.3f} {:7.2f} {:6.0f} {:8.3f} {:8.3f} {:8.3f}"
     202                else:
     203                    fmt = "{:>10s} {:8.3f} {:12.3f} {:12.3f} {:7.2f} {:6.0f} {:8.5f} {:8.5f} {:8.5f}"
    200204                self.Write(hfmt.format("h,k,l,m",tname,"F_obs","F_calc","phase","mult","sig","gam","FWHM"))
    201205                self.Write(87*'=')
     
    204208                    h,k,l,m,mult,dsp,pos,sig,gam,Fobs,Fcalc,phase = refItem[:12]
    205209                    FWHM = G2pwd.getgamFW(gam,sig)
    206                     self.Write(fmt.format(hklfmt.format(h,k,l,m),pos,Fobs,Fcalc,phase,mult,sig,gam,FWHM))
     210                    if 'T' in phasDict['Type']:
     211                        self.Write(fmt.format(hklfmt.format(h,k,l,m),pos,Fobs,Fcalc,phase,mult,sig,gam,FWHM))
     212                    else:
     213                        self.Write(fmt.format(hklfmt.format(h,k,l,m),pos,Fobs,Fcalc,phase,mult, \
     214                            np.sqrt(max(sig,0.0001))/100.,gam/100.,FWHM/100.))
    207215            else:
    208216                self.Write(85*'=')
    209217                hklfmt = "{:.0f},{:.0f},{:.0f}"
    210218                hfmt = "{:>8s} {:>8s} {:>12s} {:>12s} {:>7s} {:>6s} {:>8s} {:>8s} {:>8s}"
    211                 fmt = "{:>8s} {:8.3f} {:12.3f} {:12.3f} {:7.2f} {:6.0f} {:8.3f} {:8.3f} {:8.3f}"
     219                if 'T' in phasDict['Type']:
     220                    fmt = "{:>8s} {:8.3f} {:12.3f} {:12.3f} {:7.2f} {:6.0f} {:8.3f} {:8.3f} {:8.3f}"
     221                else:
     222                    fmt = "{:>8s} {:8.3f} {:12.3f} {:12.3f} {:7.2f} {:6.0f} {:8.5f} {:8.5f} {:8.5f}"
    212223                self.Write(hfmt.format("h,k,l",tname,"F_obs","F_calc","phase","mult","sig","gam","FWHM"))
    213224                self.Write(85*'=')
     
    216227                    h,k,l,mult,dsp,pos,sig,gam,Fobs,Fcalc,phase = refItem[:11]
    217228                    FWHM = G2pwd.getgamFW(gam,sig)
    218                     self.Write(fmt.format(hklfmt.format(h,k,l),pos,Fobs,Fcalc,phase,mult,sig,gam,FWHM))
     229                    if 'T' in phasDict['Type']:
     230                        self.Write(fmt.format(hklfmt.format(h,k,l),pos,Fobs,Fcalc,phase,mult,sig,gam,FWHM))
     231                    else:
     232                        self.Write(fmt.format(hklfmt.format(h,k,l),pos,Fobs,Fcalc,phase,mult,   \
     233                            np.sqrt(max(sig,0.0001))/100.,gam/100.,FWHM/100.))
    219234        self.CloseFile()
    220235        print(str(hist)+'reflections written to file '+str(self.fullpath))                       
Note: See TracChangeset for help on using the changeset viewer.