Changeset 5116 for trunk/GSASIIphsGUI.py


Ignore:
Timestamp:
Dec 21, 2021 2:46:05 PM (9 months ago)
Author:
vondreele
Message:

fixes to EDX calcs for LeBail?; start on Pawley for EDX.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIphsGUI.py

    r5114 r5116  
    1297212972        '''Set Pawley parameters and optionally recompute
    1297312973        '''
    12974         #GSASIIpath.IPyBreak()
    1297512974       
    1297612975        def DisablePawleyOpts(*args):
     
    1311813117        Inst = G2frame.GPXtree.GetItemPyData(G2gd.GetGPXtreeItemId(G2frame,PatternId,'Instrument Parameters'))[0]
    1311913118        Sample = G2frame.GPXtree.GetItemPyData(G2gd.GetGPXtreeItemId(G2frame,PatternId,'Sample Parameters'))
    13120         wave = G2mth.getWave(Inst)
    1312113119        const = 9.e-2/(np.pi*Sample['Gonio. radius'])                  #shifts in microns
    1312213120        gconst = 2.35482 # sqrt(8 ln 2)
     
    1312713125        try:
    1312813126            for ref in Refs:
    13129                 pos = 2.0*asind(wave/(2.0*ref[4+im]))
     13127                pos = G2lat.Dsp2pos(Inst,ref[4+im])
    1313013128                if 'Bragg' in Sample['Type']:
    1313113129                    pos -= const*(4.*Sample['Shift'][0]*cosd(pos/2.0)+ \
    1313213130                        Sample['Transparency'][0]*sind(pos)*100.0)            #trans(=1/mueff) in cm
    13133                 else:               #Debye-Scherrer - simple but maybe not right - +Layer Disp from DData?
     13131                elif 'E' not in Inst['Type'][0]:               #Debye-Scherrer - simple but maybe not right - +Layer Disp from DData?
    1313413132                    pos -= const*(Sample['DisplaceX'][0]*cosd(pos)+Sample['DisplaceY'][0]*sind(pos))
    1313513133                indx = np.searchsorted(xdata[0],pos)
     
    1314113139                    # then divide by Icorr to get F^2.
    1314213140                    ref[6+im] = (xdata[1][indx]-xdata[4][indx])*gconst/(FWHM*np.sqrt(np.pi))  #Area of Gaussian is height * FWHM * sqrt(pi)
    13143                     Lorenz = 1./(2.*sind(xdata[0][indx]/2.)**2*cosd(xdata[0][indx]/2.))           #Lorentz correction
    13144                     pola = 1.0
    13145                     if 'X' in Inst['Type']:
    13146                         pola,dIdPola = G2pwd.Polarization(Inst['Polariz.'][1],xdata[0][indx],Inst['Azimuth'][1])
     13141                    if 'E' not in Inst['Type'][0]:
     13142                        Lorenz = 1./(2.*sind(xdata[0][indx]/2.)**2*cosd(xdata[0][indx]/2.))           #Lorentz correction
     13143                        pola = 1.0
     13144                        if 'X' in Inst['Type']:
     13145                            pola,dIdPola = G2pwd.Polarization(Inst['Polariz.'][1],xdata[0][indx],Inst['Azimuth'][1])
     13146                        else:
     13147                            pola = 1.0
     13148                    # Include histo scale and volume in calculation
     13149                        ref[6+im] /= (Sample['Scale'][0] * Vst * Lorenz * pola * ref[3+im])
    1314713150                    else:
    13148                         pola = 1.0
    13149                     # Include histo scale and volume in calculation
    13150                     ref[6+im] /= (Sample['Scale'][0] * Vst * Lorenz * pola * ref[3+im])
     13151                        ref[6+im] /= (Sample['Scale'][0] * ref[3+im])
    1315113152                except IndexError:
    1315213153                    pass
Note: See TracChangeset for help on using the changeset viewer.