Changeset 5999 for trunk/MagicSoft/Mars/mtemp/mifae/macros

Timestamp:

01/25/05 17:00:22 (20 years ago)

Author:

rico

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mtemp/mifae/macros/mmcCleaning.C (modified) (9 diffs)

trunk/MagicSoft/Mars/mtemp/mifae/macros/mmcCleaning.C

@@ -6 +6 @@
 
 
-void mmcCleaning()
+void mmcCleaning(TString NoiseFilename, TString NoiseOutFilename)
 {
   //
@@ -15 +15 @@
 
   // ------------- user change -----------------
-  Char_t*  NonoiseFilename = "/local_disk/jrico/mc/Gammas/Gamma_zbin0_0_7_1000to1009_w0_nonoise.root";  // File to be used for the calibration (must be a camera file without added noise)
-  Char_t*  NoiseFilename   = "/local_disk/jrico/mc/Gammas/Gamma_zbin0_0_7_1000to1009_w0.root"; // File to be analyzed
+  Char_t*  NonoiseFilename = "/local_disk/jrico/mc/Gamma_zbin0_0_7_1000to1009_w0.root";
+  //  Char_t*  NoiseFilename   = "/mnt/wdflix/root_0.73mirror/wuerzburg/gammas/Gamma_zbin3_0_7_1160to1169_w0.root"; // File to be analyzed
 
-  Char_t* NoiseOutFilename   = "/local_disk/jrico/mc/Gammas/prueba.root";   // Output file name
+  //  Char_t* NoiseOutFilename   = "/mnt/users/jrico/prueba.root";   // Output file name
 
-  // Signal extractor
+  const Int_t wsize=2;
+  // Signal sigextract
   // (other extraction methods can be used)
-  MExtractFixedWindowPeakSearch sigextract;
-
+  MExtractor* sigextract = new MExtractTimeAndChargeSpline();
+  ((MExtractTimeAndChargeSpline*)sigextract)->SetTimeType(MExtractTimeAndChargeSpline::kHalfMaximum);
+  ((MExtractTimeAndChargeSpline*)sigextract)->SetChargeType(MExtractTimeAndChargeSpline::kIntegral);      
+  ((MExtractTimeAndChargeSpline*)sigextract)->SetRiseTime((Float_t)wsize*0.25); 
+  ((MExtractTimeAndChargeSpline*)sigextract)->SetFallTime((Float_t)wsize*0.75);
+  
   // Define FADC slices to be integrated in high and low gain:
-  sigextract.SetRange(1, 14, 2, 14);
-  sigextract.SetWindows(6,6,4);
+  sigextract->SetRange(1, 11, 2, 12);
 
   // Defines when to switch to low gain
-  sigextract.SetSaturationLimit(240);
+  sigextract->SetSaturationLimit(240);
 
   // ---------------------------------------------------------------------
@@ -82 +86 @@
   calib.SetCalibrationMode(MCalibrate::kFfactor);
 
-  MImgCleanStd clean(2.5,2.0);
+  MImgCleanStd clean(4.0,3.5);
   clean.SetCleanRings(1);
   //  clean.SetRemoveSingle(kFALSE);
@@ -97 +101 @@
   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(sigextract);
   tlist.AddToList(&mccalibupdate);
   tlist.AddToList(&calib);
@@ -118 +122 @@
 
   MEvtLoop evtloop;
-  evtloop.SetProgressBar(&bar);
+  //  evtloop.SetProgressBar(&bar);
   evtloop.SetParList(&plist);
 
@@ -133 +137 @@
   MIslands      isl;
   MArrivalTimeCam timecam;      
+  MTopology     topology;
   plist.AddToList(&isl);
   plist.AddToList(&timecam);
+  plist.AddToList(&topology);
 
   MArrivalTimeCalc2 timecalc;
@@ -143 +149 @@
   islclean.SetInputName("MIslands");
   islclean.SetMethod(1);
+  MTopologyCalc     topcalc;
 
   MReadMarsFile read2("Events");
@@ -152 +159 @@
   tlist.AddToListBefore(&timecalc,&mccalibupdate,"All");
   tlist.AddToListBefore(&islandcalc,&hcalc,"All");
-  tlist.AddToListBefore(&islclean,&hcalc,"All");
+  tlist.AddToListBefore(&topcalc,&hcalc,"All");
 
   MWriteRootFile write(NoiseOutFilename); // Writes output
+  write.AddContainer("MRawRunHeader", "RunHeaders");
+  write.AddContainer("MMcRunHeader",  "RunHeaders");
+  write.AddContainer("MSrcPosCam",    "RunHeaders");
+  write.AddContainer("MGeomCam",      "RunHeaders");
+  write.AddContainer("MMcConfigRunHeader",  "RunHeaders");
+  write.AddContainer("MMcCorsikaRunHeader", "RunHeaders");
+  write.AddContainer("MMcFadcHeader",  "RunHeaders");
+  write.AddContainer("MMcTrigHeader",  "RunHeaders");
+
   write.AddContainer("MRawRunHeader", "Events");
   write.AddContainer("MMcEvt",        "Events");
@@ -160 +176 @@
   write.AddContainer("MHillas",       "Events");
   write.AddContainer("MHillasSrc",    "Events");
+  write.AddContainer("MHillasExt",    "Events");
+  write.AddContainer("MConcentration","Events");
   write.AddContainer("MImagePar",     "Events");
   write.AddContainer("MNewImagePar",  "Events");
   write.AddContainer("MIslands",      "Events");
+  write.AddContainer("MPointingPos",  "Events");
+  write.AddContainer("MTopology",      "Events");
 
   tlist.RemoveFromList(&mccalibcalc);