Changeset 3124 for trunk/MagicSoft/Mars/macros

Timestamp:

02/12/04 19:54:27 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/calibration.C (modified) (4 diffs)

trunk/MagicSoft/Mars/macros/calibration.C

@@ -26 +26 @@
 //const TString calfile = "/mnt/users/mdoro/Mars/Data/20040201_1441*_C_OffMrk421-1_E.root";
 
-const TString pedfile = "/mnt/Data/rootdata/CrabNebula/2004_01_27/20040126_12386_P_Cab-On_E.root";
-const TString calfile = "/mnt/Data/rootdata/CrabNebula/2004_01_27/20040126_12525_C_Cab-On_E.root";
+const TString pedfile = "/mnt/Data/rootdata/CrabNebula/2004_02_10/20040210_14607_P_CrabNebula_E.root";
+const TString calfile = "/mnt/Data/rootdata/CrabNebula/2004_02_10/20040210_14608_C_CrabNebula_E.root";
+
+//const TString pedfile = "/mnt/Data/rootdata/CrabNebula/2004_01_26/20040125_10412_P_Crab-On_E.root";
+//const TString calfile = "/mnt/Data/rootdata/CrabNebula/2004_01_26/20040125_1041*_C_Crab-On_E.root";
 
 //const TString pedfile = "/mnt/Data/rootdata/Miscellaneous/2003_12_19/20031218_03522_P_Park_E.root";
@@ -220 +223 @@
     MCalibrationCalc     calcalc;
     
+    //
+    // Set the range (other than default) 
+    // of FADC slices for the blind pixel
+    //
+    // calcalc.SetBlindPixelRange(10,25);
+
+    //
+    // Set the cut upon which a superposition of the blind pixel 
+    // FADC slices will be filled into the SinglePHE histogram
+    //
+    calcalc.SetBlindPixelSinglePheCut(500);
+
+    //
+    // Skip the HiGain vs. LoGain calibration
+    // 
+    calcalc.SkipHiLoGainCalibration();
+
+    //
+    // As long, as we don't have digital modules, 
+    // we have to set the color of the pulser LED by hand
+    //
+    calcalc.SetPulserColor(MCalibrationCalc::kECT1);
+
+    //
+    // In case, we want to exclude a pre-defined list of bad pixels:
+    // (This is a preliminary feature)
+    //
+    // calcalc.ExcludePixelsFromAsciiFile("badpixels.dat");
+    
+    //
+    // In case, you want to skip the blind pixel method: 
+    // (NOT RECOMMENDED!!!)
+    //
+    // calcalc.SkipBlindPixelFit();
+
+    //
+    // In case, you want to skip the quality checks
+    // (NOT RECOMMENDED!!!)
+    //
+    // calcalc.SkipQualityChecks();
+
+    //
+    // In case, we want to apply another fit function to the 
+    // blind pixel 
+    //
+    MCalibrationBlindPix *bp = calcam.GetBlindPixel();
+//    bp->ChangeFitFunc(MHCalibrationBlindPixel::kEPolya);
+    bp->ChangeFitFunc(MHCalibrationBlindPixel::kEPoisson5);
+
     // 
     // Apply a filter against cosmics
@@ -226 +278 @@
     MFCosmics            cosmics;
     MContinue            cont(&cosmics);
-
-    //
-    // Skip the HiGain vs. LoGain calibration
-    // 
-    calcalc.SkipHiLoGainCalibration();
-
-    //
-    // Making the step size a bit bigger, gives us 
-    // faster results
-    // 
-    //    timecalc.SetStepSize(0.05);
-
-    //
-    // As long, as we don't have digital modules, 
-    // we have to set the color of the pulser LED by hand
-    //
-    calcalc.SetPulserColor(MCalibrationCalc::kECT1);
-
-    //
-    // In case, we want to exclude a pre-defined list of bad pixels:
-    // (This is a preliminary feature)
-    //
-   calcalc.ExcludePixelsFromAsciiFile("badpixels.dat");
-    
-    //
-    // In case, you want to skip the blind pixel method: 
-    // (NOT RECOMMENDED!!!)
-    //
-    // calcalc.SkipBlindPixelFit();
-
-    //
-    // In case, you want to skip the quality checks
-    // (NOT RECOMMENDED!!!)
-    //
-    // calcalc.SkipQualityChecks();
-
-    //
-    // In case, we want to apply another fit function to the 
-    // blind pixel 
-    //
-    MCalibrationBlindPix *bp = calcam.GetBlindPixel();
-//    bp->ChangeFitFunc(MHCalibrationBlindPixel::kEPolya);
-//    bp->ChangeFitFunc(MHCalibrationBlindPixel::kEPoisson4);
 
     tlist2.AddToList(&read2);
@@ -298 +307 @@
     // print the most important results of all pixels
     //
-    calcam.Print();
+    //    calcam.Print();
 
     //
     // just one example how to get the plots of individual pixels
     //
-    calcam[543].DrawClone();
+    //    calcam[563].DrawClone();
+    //    calcam[564].DrawClone();
 
     // Create histograms to display