Changeset 7188 for trunk/MagicSoft/Mars/macros/starmc.C

Timestamp:

07/13/05 19:06:26 (19 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/starmc.C (modified) (7 diffs)

trunk/MagicSoft/Mars/macros/starmc.C

@@ -52 +52 @@
   // differences in gain of outer pixels)
   //
-  CalibrationFilename = new TString("/data1/magic/mc_data/root/Period025/gammas_nonoise/Gamma_zbin0_*.root");
+  CalibrationFilename = new TString("nonoise/Gamma_zbin0_0_7_1000to1009_w0.root");
   // File to be used in the calibration (must be a camera file without added noise)
 
-
-  Char_t* AnalysisFilename = "Gamma_zbin1_0_*.root";  // File to be analyzed
-
+  Char_t* AnalysisFilename = "Gamma_zbin0_0_*.root";  // File to be analyzed
 
 
@@ -76 +74 @@
   // USER CHANGE: tail cuts for image analysis
 
-  Float_t CleanLev[2] = {7., 5.}; 
-  MImgCleanStd  clean(CleanLev[0], CleanLev[1]); // Applies tail cuts to image.
-  clean.SetMethod(MImgCleanStd::kAbsolute);      // In absolute units (phot or phe- as chosen below)
-
+  Float_t CleanLev[2] = {10., 5.}; // Tail cuts for the analysis loop
+  MImgCleanStd clean2(CleanLev[0], CleanLev[1]); // Applies tail cuts to image.
+  clean2.SetMethod(MImgCleanStd::kAbsolute); 
 
   //  USER CHANGE: signal extraction
-  //
-  //  MExtractFixedWindowPeakSearch sigextract;
-  //  sigextract.SetWindows(6, 6, 4);
 
   //  MExtractTimeAndChargeDigitalFilter sigextract;
@@ -91 +85 @@
 
   MExtractTimeAndChargeSpline sigextract;
+  sigextract.SetChargeType(MExtractTimeAndChargeSpline::kIntegral);
   sigextract.SetRiseTimeHiGain(0.5);
   sigextract.SetFallTimeHiGain(0.5);
@@ -96 +91 @@
   // USER CHANGE: high to low gain ratio. DEPENDS on EXTRACTOR!!
   // One should calculate somewhere else what this factor is for each extractor!
-  Float_t hi2low = 12.; // tentative value for spline with risetime 0.5, fall time 0.5
+
+  //  Float_t hi2low = 10.83; 
+  // value for digital filter, HG window 4, LG window 6
+
+  Float_t hi2low = 11.28; 
+  // value for spline with risetime 0.5, fall time 0.5, low gain stretch 1.5
 
 
@@ -108 +108 @@
 
 
+  MImgCleanStd  clean(0.,0.); 
+  // All pixels containing any photon will be accepted. This is what we want
+  // for the calibration loop (since no noise is added)
+  clean.SetMethod(MImgCleanStd::kAbsolute); 
+  // In absolute units (phot or phe- as chosen below)
+
   MSrcPosCam src;
   //
-  // ONLY FOR WOBBLE MODE!!
-  //   src.SetX(120.);  // units: mm
+  // ONLY FOR WOBBLE MODE!! : set the rigt source position on camera!
+  //   src.SetX(120.);   // units: mm. Value for MC w+ files
+  //   src.SetX(-120.);  // units: mm. Value for MC w- files
 
   src.SetReadyToSave();
@@ -244 +251 @@
   //
 
+
   MProgressBar bar;
   bar.SetWindowName("Calibrating...");
@@ -273 +281 @@
   tlist.RemoveFromList(&read);
 
+  tlist.AddToListBefore(&clean2, &clean);
+  tlist.RemoveFromList(&clean);
+
   //
   // Analyzed only the desired fraction of events, skip the rest: