Changeset 2784


Ignore:
Timestamp:
Apr 15, 2017 10:46:00 PM (5 years ago)
Author:
vondreele
Message:

fix REFD plotting
put neutron reflectometry reader back in import routine

Location:
trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIplot.py

    r2777 r2784  
    20132013                Y = xye[1]+offsetY*N*Ymax/100.0
    20142014        elif plottype in ['SASD','REFD']:
    2015             B = xye[5]
     2015            if plottype == 'SASD':
     2016                B = xye[5]
     2017            else:
     2018                B = np.zeros_like(xye[5])
    20162019            if G2frame.plotStyle['sqPlot']:
    20172020                Y = xye[1]*Sample['Scale'][0]*(1.05)**(offsetY*N)*X**4
  • trunk/imports/G2rfd_xye.py

    r2783 r2784  
    122122        return True
    123123
    124 class txt_nmXRayReaderClass(G2IO.ImportReflectometryData):
    125     'Routines to import X-ray q REFD data from a .xrfd or .xdat file, q in nm-1'
     124class txt_NeutronReaderClass(G2IO.ImportReflectometryData):
     125    'Routines to import neutron q REFD data from a .nrfd or .ndat file'
    126126    def __init__(self):
    127127        super(self.__class__,self).__init__( # fancy way to self-reference
    128             extensionlist=('.xrfd','.xdat'),
     128            extensionlist=('.nrfd','.ndat'),
    129129            strictExtension=False,
    130             formatName = 'q (nm-1) step X-ray QRE data',
    131             longFormatName = 'q (nm-1) stepped X-ray text data file in Q,R,E order; E optional'
     130            formatName = 'q (A-1) step neutron QRE data',
     131            longFormatName = 'q (A-1) stepped neutron text data file in Q,R,E order; E optional'
    132132            )
    133133
     
    176176                try:
    177177                    data = [float(val) for val in vals]
    178                     x.append(float(data[0])/10.)        #convert nm-1 to A-1
     178                    x.append(float(data[0]))
    179179                    f = float(data[1])
    180180                    if f <= 0.0:
    181                         x.pop()
     181                        del x[-1]
    182182                        continue
    183183                    elif len(vals) > 2:
     
    203203                    pass
    204204        self.instdict['wave'] = wave
    205         self.instdict['type'] = 'RXC'
     205        self.instdict['type'] = 'RNC'
    206206        x = np.array(x)
    207207        self.reflectometrydata = [
     
    211211            np.zeros(N), # calc. intensities (zero)
    212212            np.zeros(N), # obs-calc profiles
    213             np.array(sq), # fix bkg
     213            np.array(sq), # Q FWHM
    214214            ]
    215215        self.reflectometryentry[0] = filename
     
    220220
    221221        return True
     222
Note: See TracChangeset for help on using the changeset viewer.