Changeset 7d6351e in sasview for sansdataloader/src/sans/dataloader/readers/hfir1d_reader.py
- Timestamp:
- Apr 27, 2012 11:22:31 AM (12 years ago)
- Branches:
- master, ESS_GUI, ESS_GUI_Docs, ESS_GUI_batch_fitting, ESS_GUI_bumps_abstraction, ESS_GUI_iss1116, ESS_GUI_iss879, ESS_GUI_iss959, ESS_GUI_opencl, ESS_GUI_ordering, ESS_GUI_sync_sascalc, costrafo411, magnetic_scatt, release-4.1.1, release-4.1.2, release-4.2.2, release_4.0.1, ticket-1009, ticket-1094-headless, ticket-1242-2d-resolution, ticket-1243, ticket-1249, ticket885, unittest-saveload
- Children:
- 4a96b8b
- Parents:
- f60a8c2
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
sansdataloader/src/sans/dataloader/readers/hfir1d_reader.py
rad8034f r7d6351e 1 1 """ 2 HFIR 1D 4-column data reader 3 """ 2 4 ##################################################################### 3 5 #This software was developed by the University of Tennessee as part of the 4 6 #Distributed Data Analysis of Neutron Scattering Experiments (DANSE) 5 #project funded by the US National Science Foundation. 7 #project funded by the US National Science Foundation. 6 8 #See the license text in license.txt 7 9 #copyright 2008, University of Tennessee 8 10 ###################################################################### 9 10 11 import numpy 11 12 import os … … 24 25 """ 25 26 ## File type 26 type_name = "HFIR 1D" 27 type_name = "HFIR 1D" 27 28 ## Wildcards 28 29 type = ["HFIR 1D files (*.d1d)|*.d1d"] 29 30 ## List of allowed extensions 30 ext = ['.d1d'] 31 ext = ['.d1d'] 31 32 32 33 def read(self, path): … … 42 43 """ 43 44 if os.path.isfile(path): 44 basename 45 basename = os.path.basename(path) 45 46 root, extension = os.path.splitext(basename) 46 47 if extension.lower() in self.ext: 47 48 try: 48 input_f = 49 except 49 input_f = open(path,'r') 50 except: 50 51 raise RuntimeError, "hfir1d_reader: cannot open %s" % path 51 52 buff = input_f.read() 52 53 lines = buff.split('\n') 53 x 54 y 54 x = numpy.zeros(0) 55 y = numpy.zeros(0) 55 56 dx = numpy.zeros(0) 56 57 dy = numpy.zeros(0) … … 84 85 85 86 if data_conv_i is not None: 86 _y = data_conv_i(_y, units=output.y_unit) 87 _dy = data_conv_i(_dy, units=output.y_unit) 87 _y = data_conv_i(_y, units=output.y_unit) 88 _dy = data_conv_i(_dy, units=output.y_unit) 88 89 89 x = numpy.append(x, _x)90 y 91 dx 92 dy = numpy.append(dy, _dy)90 x = numpy.append(x, _x) 91 y = numpy.append(y, _y) 92 dx = numpy.append(dx, _dx) 93 dy = numpy.append(dy, _dy) 93 94 except: 94 95 # Couldn't parse this line, skip it 95 96 pass 96 97 97 98 98 # Sanity check 99 99 if not len(y) == len(dy): … … 120 120 output.yaxis("\\rm{Intensity}", output.y_unit) 121 121 else: 122 output.yaxis("\\rm{Intensity}", "cm^{-1}")122 output.yaxis("\\rm{Intensity}", "cm^{-1}") 123 123 124 124 # Store loading process information 125 output.meta_data['loader'] = self.type_name 125 output.meta_data['loader'] = self.type_name 126 126 return output 127 127 else: 128 128 raise RuntimeError, "%s is not a file" % path 129 129 return None 130 131 if __name__ == "__main__":132 reader = Reader()133 print reader.read("../test/S2-30dq.d1d")134 135 136
Note: See TracChangeset
for help on using the changeset viewer.