Changeset 3852 for trunk/MagicSoft/Mars/macros

Timestamp:

04/27/04 18:44:34 (21 years ago)

Author:

moralejo

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/macros

Files:

• mccalibrate.C (modified) (7 diffs)
• starmc2.C (modified) (3 diffs)

trunk/MagicSoft/Mars/macros/mccalibrate.C

@@ -28 +28 @@
 //  MMCALIBRATE - Calibration of MC data
 //
-//  This macro is a version of the standard procedure to convert raw MC data 
-//  into calibrated data (photons per pixel) 
+//  This macro converts raw MC data into calibrated data (photons per pixel) 
 //
 /////////////////////////////////////////////////////////////////////////////
@@ -43 +42 @@
   //
 
-
   // ------------- user change -----------------
   TString* CalibrationFilename;
-  CalibrationFilename = new TString("../../gammas_nonoise/Gamma_zbin0_0*.root");
-  // File to be used for the calibration (must be a camera file without added noise)
-
-  Char_t* AnalysisFilename = "Gamma_zbin0*.root";  // File to be analyzed
-
-  Char_t* OutFilename      = "calibrated_data.root";       // Output file name
-
-
-  Int_t BinsHigh[2] = {5, 10}; // First and last FADC bin of the range to be integrated,
-  Int_t BinsLow[2]  = {5, 10}; // for high and low gain respectively.
-
-  // -------------------------------------------
-
+  CalibrationFilename = new TString("../../gammas_nonoise/Gamma_zbin0_0*.root");  // File to be used for the calibration (must be a camera file without added noise)
+
+  Char_t* AnalysisFilename = "Proton_zbin0_0*.root";  // File to be analyzed
+
+  Char_t* OutFilename      = "calibrated_data.root";  // Output file name
+
+
+  MExtractSignal    sigextract;  
+  // (other extraction methods can be used)
+
+  sigextract.SetSaturationLimit(240);
+  // Defines when to switch to low gain
+
+  // Define FADC slices to be integrated in high and low gain:
+  sigextract.SetRange(5, 10, 5, 10);
+
+  // ---------------------------------------------------------------------
   //
   // Create a empty Parameter List and an empty Task List
@@ -74 +76 @@
 
   MBadPixelsCam badpix;
-  plist.AddToList(&badpix);
-  
-
+  plist.AddToList(&badpix);  // Not used for now.
+ 
   //
   // Now setup the tasks and tasklist:
@@ -88 +89 @@
   read.DisableAutoScheme();
 
-  MGeomApply geom; // Reads in geometry from MC file and sets the right sizes for
+  MGeomApply geom; 
+  // Reads in geometry from MC file and sets the right sizes for
   // several parameter containers.
 
   MMcPedestalCopy   pcopy; 
-  // Copies pedestal data from the MC file run fadc header to the MPedestalCam container.
-
-  MExtractSignal    sigextract;
-  sigextract.SetSaturationLimit(240);
-
-  // Define ADC slices to be integrated in high and low gain:
-  sigextract.SetRange(BinsHigh[0], BinsHigh[1], BinsLow[0], BinsLow[1]);
+  // Copies pedestal data from the MC file run fadc header to the 
+  // MPedestalCam container.
+
+  MPointingPosCalc pointcalc;
+  // Creates MPointingPos object and fill it with the telescope orientation
+  // information taken from MMcEvt.
 
   MMcCalibrationUpdate  mccalibupdate;
 
-  MCalibrate calib; // Transforms signals from ADC counts into photons.
+  MCalibrate calib; 
+  // MCalibrate transforms signals from ADC counts into photons. In the first
+  // loop it applies a "dummy" calibration supplied by MMcCalibrationUpdate, just 
+  // to equalize inner and outer pixels. At the end of the first loop, in the
+  // PostProcess of MMcCalibrationCalc (see below) the true calibration constants
+  // are calculated.
+
   calib.SetCalibrationMode(MCalibrate::kFfactor);
 
-  //    MBlindPixelCalc   blind;
-  //    blind.SetUseInterpolation();
-
+  MImgCleanStd clean;
   //
   // Applies tail cuts to image. Since the calibration is performed on 
@@ -114 +119 @@
   // be rejected).
   //
-  MImgCleanStd      clean;
-
-
-  MHillasCalc       hcalc; // Calculates Hillas parameters not dependent on source position.
-
-  MMcCalibrationCalc mccalibcalc;
+
+  MHillasCalc hcalc; // Calculates Hillas parameters not dependent on source position.
+
+  MMcCalibrationCalc mccalibcalc; 
+  // Calculates calibration constants to convert from ADC counts to photons.
+  
 
   tlist.AddToList(&read);
   tlist.AddToList(&geom);
   tlist.AddToList(&pcopy);
-
+  tlist.AddToList(&pointcalc);
   tlist.AddToList(&sigextract);
   tlist.AddToList(&mccalibupdate);
   tlist.AddToList(&calib);
   tlist.AddToList(&clean);
-  //    tlist.AddToList(&blind);
   tlist.AddToList(&hcalc);
 
@@ -150 +154 @@
   write.AddContainer("MMcEvt",        "Events");
   write.AddContainer("MMcTrig",       "Events");
+  write.AddContainer("MPointingPos",  "Events");
   write.AddContainer("MRawEvtHeader", "Events");
   write.AddContainer("MCerPhotEvt",   "Events");
@@ -169 +174 @@
       if (!evtloop.Eventloop())
         return;
-      mccalibcalc->GetHist()->Write();
+      mccalibcalc->GetHistADC2PhotEl()->Write();
+      mccalibcalc->GetHistPhot2PhotEl()->Write();
+      // Writes out the histograms used for calibration.
     }
 
   //
-  // Second loop: analysis loop
+  // Second loop: apply calibration factors to MC events in the 
+  // file to be anlyzed:
   //
 

trunk/MagicSoft/Mars/macros/starmc2.C

@@ -38 +38 @@
 void starmc2()
 {
-  Char_t* AnalysisFilename = "Calibrated_run.root"; // File to be analyzed
+  Char_t* AnalysisFilename = "calibrated_data.root"; // File to be analyzed
   Char_t* OutFilename      = "star.root";        // Output file name
 
@@ -77 +77 @@
   tlist.AddToList(&read);
   tlist.AddToList(&clean);
-  //    tlist.AddToList(&blind);
   tlist.AddToList(&hcalc);
   tlist.AddToList(&scalc); // Calculates Source-dependent Hillas parameters
@@ -89 +88 @@
   write.AddContainer("MSrcPosCam",    "RunHeaders");
   write.AddContainer("MMcEvt",        "Events", kFALSE);
+  write.AddContainer("MPointingPos",  "Events");
   write.AddContainer("MHillas",       "Events");
   write.AddContainer("MHillasExt",    "Events");