Changeset 13380 for fact/tools/pyscripts/pyfact/pyfact.py

Timestamp:

04/18/12 17:07:40 (13 years ago)

Author:

neise

Message:

new option in RawData.__init__() one can set do_calibration=False. and the amplitude calibration comes without np.roll() ... thx to TPK for the idea

File:

• fact/tools/pyscripts/pyfact/pyfact.py (modified) (4 diffs)

fact/tools/pyscripts/pyfact/pyfact.py

@@ -83 +83 @@
                 user_action_calib=lambda acal_data, data, blm, tom, gm, scells, nroi: None,
                 baseline_file_name='',
-                return_dict = None):
+                return_dict = None,
+                do_calibration = True):
         """ initialize object
         
@@ -107 +108 @@
             print 'and some time later, the option will not be supported anymore'
             return_dict = False
-
         self.return_dict = return_dict
         
+        self.do_calibration = do_calibration
+
         self.data_file_name = data_file_name
         self.calib_file_name = calib_file_name
@@ -198 +200 @@
             raise StopIteration
         else:
-            self.calibrate_drs_amplitude()
+            if self.do_calibration == True:
+                self.calibrate_drs_amplitude()
 
         #print 'nevents = ', self.nevents, 'event_id = ', self.event_id.value
@@ -229 +232 @@
         gm  = np.reshape(self.gm,  (self.npix, self.Ngm) )
         
+        # double the calibration constants
+        blm = np.hstack((blm,blm))
+        gm = np.hstack((gm,gm))
+        
         for pixel in range( self.npix ):
             # rotate the pixel baseline mean to the Data startCell
-            blm_pixel = np.roll( blm[pixel,:], -self.start_cells[pixel] )
+            #blm_pixel = np.roll( blm[pixel,:], -self.start_cells[pixel] )
+            blm_pixel = blm[pixel, self.start_cells[pixel]:self.start_cells[pixel]+self.nroi]
             # rotate the pixel gain mean to the Data startCell
-            gm_pixel = np.roll( gm[pixel,:], -self.start_cells[pixel] )
+            #gm_pixel = np.roll( gm[pixel,:], -self.start_cells[pixel] )
+            gm_pixel = gm[pixel, self.start_cells[pixel]:self.start_cells[pixel]+self.nroi]
             # the 'trigger offset mean' does not need to be rolled
             # on the contrary, it seems there is an offset in the DRS data,
             # which is related to its distance to the startCell, not to its 
             # distance to the beginning of the physical pipeline in the DRS chip
-            tom_pixel = tom[pixel,:]
+            tom_pixel = tom[pixel,0:self.nroi ]
             
-            acal_data[pixel,:] -= blm_pixel[0:self.nroi]
-            acal_data[pixel,:] -= tom_pixel[0:self.nroi]
-            acal_data[pixel,:] /= gm_pixel[0:self.nroi]
+            acal_data[pixel,:] -= blm_pixel
+            acal_data[pixel,:] -= tom_pixel
+            acal_data[pixel,:] /= gm_pixel
             
         self.acal_data = acal_data * 1907.35