Changeset 13400 for fact/tools/pyscripts/sandbox/kraehenb/CalFitsTest.py

Timestamp:

04/20/12 14:41:09 (13 years ago)

Author:

kraehenb

Message:

Completed the fast implementation to read calibrated data: calfits.h, CalFitsTest.py.

File:

• fact/tools/pyscripts/sandbox/kraehenb/CalFitsTest.py (modified) (1 diff)

fact/tools/pyscripts/sandbox/kraehenb/CalFitsTest.py

@@ +1 @@
+#!/usr/bin/python -tt
+# ********************************
+# Test script for the CalFits class
+# 
+# written by Thomas Kraehenbuehl, ETH Zurich
+# tpk@phys.ethz.ch, +41 44 633 3973
+# April 2012
+# ********************************
+
 datafilename = '/fact/raw/2012/04/17/20120417_004.fits.gz'
 calibfilename = '/fact/raw/2012/04/17/20120417_003.drs.fits.gz'
-#data_filename =  '/fact/raw/2012/01/31/20120131_018.fits.gz'
-#calib_filename = '/fact/raw/2012/01/31/20120131_014.drs.fits.gz'
 
-DoCalFits=True
+from ROOT import gSystem
+gSystem.Load("calfits_h.so")
+from ROOT import *
+print "Testing object creation: "
+caltest = CalFits(datafilename,calibfilename)
 
-if DoCalFits:
-        from ROOT import gSystem
-        gSystem.Load("CalFits_h.so")
-        from ROOT import *
-        print "Testing object creation: "
-        caltest = CalFits(datafilename,calibfilename)
-        print "Testing the python access to the function: "
-        print caltest.data_npix
+print "Common variables:"
+print "ROI: ", caltest.nroi
+print "#Pix: ", caltest.npix
+print "Number of events: ", caltest.nevents
+print
 
-        print "Calibfile ROI: ", caltest.calib_nroi
-        print "BaselineMean: ", caltest.calibfile.GetN("BaselineMean")
-        print "GainMean: ", caltest.calibfile.GetN("GainMean")
-        print "TriggerOffsetMean: ", caltest.calibfile.GetN("TriggerOffsetMean")
-        print "Datafile ROI: ", caltest.data_nroi
-        print "Data: ", caltest.datafile.GetN("Data")
-        print "StartCellData: ", caltest.datafile.GetN("StartCellData")
-        print "Direct datafile access: ", caltest.datafile.GetN("StartCellData")
+print "Information per Event:"
+caltest.GetCalEvent()
+print "Event ID: ", caltest.event_id
+print "Trigger type: ", caltest.event_triggertype
+print "Uncalibrated data: ", caltest.event_data
+print "Calibrated data: ", caltest.event_caldata
+print "Board times: ", caltest.event_boardtimes
+print "Trigger offsets: ", caltest.event_offset
+print
 
-#       print caltest.GetCalEvent()
-        for i in range(1000):
-                caltest.GetCalEvent()
-                print caltest.data_eventid
-        print 
+print "Examples of other information"
+print "Calibfile ROI: ", caltest.calib_nroi
+print "Column size BaselineMean: ", caltest.calibfile.GetN("BaselineMean")
+print "Datafile ROI: ", caltest.data_nroi
+print "Data: ", caltest.datafile.GetN("Data")
+print "StartCellData: ", caltest.datafile.GetN("StartCellData")
+print "Direct datafile access: ", caltest.datafile.GetN("StartCellData")
+print
+print "Columns of the datafile: "
+caltest.datafile.PrintColumns()
 
-        del caltest
+while caltest.GetCalEvent():
+    print caltest.event_id, caltest.event_triggertype, caltest.event_caldata[0]
+print
 
-else:
-        from pyfact import RawData
-        pyfacttest = RawData(datafilename, calibfilename, return_dict = True)
-        for i,event in enumerate(pyfacttest):
-                if i>0:
-                        break
-                print event['acal_data'][0][0:5]
+del caltest