Source Code for Package reflectometry.reduction.examples.e3a12

 1  """ 
 2  Sample data from NG-1: 
 3   
 4      e3a12 rem (71g) cofe2o4,cofe(10nm),ta 
 5  """ 
 6   
 7  import numpy, os 
 8  from reflectometry.reduction.data import Data,PolarizedData 
 9   
10  PATH = os.path.dirname(os.path.realpath(__file__)) 
11   
12 -def loadABCD(file): 
13      data = PolarizedData(xlabel="Qz",xunits="invA",vlabel="counts") 
14      for ext,d in [('A',data.pp), ('B',data.pm), ('C',data.mp), ('D',data.mm)]: 
15          A = numpy.loadtxt(os.path.join(PATH,file+ext)) 
16          d.x,d.v,d.dv = A[:,0],A[:,1],A[:,2] 
17      return data 
18   
19   
20 -def slits(): 
21      data = loadABCD('e3a12042.slit') 
22      data.set(xlabel="slit1",xunits="mm",vlabel="counts") 
23      for ext,d in [('A',data.pp), ('B',data.pm), ('C',data.mp), ('D',data.mm)]: 
24          A = numpy.loadtxt(os.path.join(PATH,'e3a12042.slit'+ext)) 
25          d.x,d.v,d.dv = A[:,0],A[:,1],A[:,2] 
26      return data 
27   
28 -def spec(): 
29      data = PolarizedData(xlabel="Qz",xunits="invA",vlabel="counts") 
30      for ext,d in [('A',data.pp), ('B',data.pm), ('C',data.mp), ('D',data.mm)]: 
31          A = numpy.loadtxt(os.path.join(PATH,'e3a12026.spec'+ext)) 
32          d.x,d.v,d.dv = A[:,0],A[:,1],A[:,2] 
33      return data 
34   
35   
36  if __name__ == "__main__": 
37      print slits() 
38