Changeset 12362 for fact/tools/rootmacros

Timestamp:

11/02/11 13:50:54 (13 years ago)

Author:

kraehenb

Message:

New version of the calscope-script which is easier to adapt. The included file TPKplotevent was dropped for the FCalibrateEvent.C, which takes an event and returns an array with the DRS-calibrated pipeline.

Location:

fact/tools/rootmacros

Files:

• FCalibrateEvent.c (added)
• calscope.C (modified) (4 diffs)

fact/tools/rootmacros/calscope.C

@@ -8 +8 @@
 #define HAVE_ZLIB
 #include "fits.h"
-#include "TPKplotevent.c"
 #include "FOpenDataFile.c"
 #include "FOpenCalibFile.c"
+#include "FCalibrateEvent.c"
 
-int calscope(const char *name = "../raw/20110916_025.fits", const char *drsname = "../raw/20110916_024.drs.fits", size_t eventnr = 0, size_t pixelnr = 0)
+int calscope(const char *name, const char *drsname, size_t eventnr, size_t pixelnr)
 {
 //******************************************************************************
@@ -20 +20 @@
 //Example call in ROOT:
 //root [74] .x calscope.C++("20110804_024.fits","20110804_023.drs.fits",10,1348)
-//T. Krähenbühl, August 2011, tpk@phys.ethz.ch
+//T. Kr�henb�hl, August 2011, tpk@phys.ethz.ch
 //******************************************************************************
 
@@ -64 +64 @@
         }
         
-        
 //-------------------------------------------
 //Create the title
@@ -72 +71 @@
         
 //-------------------------------------------
-//Get the event
+//Get the calibrated event
 //-------------------------------------------
+        vector<float> calevent(data_px*data_roi); //Vector for the calibrated event
+
         cout << "--------------------- Data --------------------" << endl;
         datafile.GetRow(eventnr);
         cout << "Event number: " << data_num << endl;
-                
+        
+        FCalibrateEvent(data, data_offset, drs_basemean, drs_gainmean, drs_triggeroffsetmean, calevent, data_px, data_roi);
+        
 //-------------------------------------------
 //Draw the data
 //-------------------------------------------
-        TPKplotevent(title, data, data_offset, drs_basemean, drs_gainmean, drs_triggeroffsetmean, data_roi, pixelnr);
+        TCanvas *canv = new TCanvas( "canv", "Mean values of the first event", 100, 10, 700, 500 );
+        TProfile *pix = new TProfile("pix", title, 1024, -0.5, 1023.5);
+        
+        for (UInt_t k=0; k<data_roi; k++)
+        {
+                pix->Fill(k,calevent[pixelnr*data_roi+k]);
+        }
+        
+        pix->Draw();
+        canv->Modified();
+        canv->Update();
+        
         return 0;
 }