Changeset 3339 for trunk/MagicSoft/Mars/macros

Timestamp:

02/27/04 06:52:39 (21 years ago)

Author:

wittek

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/ONOFFAnalysis.C (modified) (12 diffs)

trunk/MagicSoft/Mars/macros/ONOFFAnalysis.C

@@ -53 +53 @@
         //.....
         MBinning *binsb = new MBinning("BinningSigmabar");
-        binsb->SetEdges( 100,  0.0,  5.0);
+        binsb->SetEdges( 100,  0.0,  50.0);
         plist->AddToList(binsb);
 
@@ -80 +80 @@
 
         MBinning *binsdiff = new MBinning("BinningDiffsigma2");
-        binsdiff->SetEdges(100, -10.0, 15.0);
+        binsdiff->SetEdges(100, -300.0, 500.0);
         plist->AddToList(binsdiff);
 
@@ -181 +181 @@
       //const char *onfile  = "MCerPhot_output";
       //const char *onfile  = "*.ON";
-      //const char *onfile  = "12*.ON";
-      const char *onfile  = "12410.ON";
+      const char *onfile  = "12*.ON";
+      //const char *onfile  = "1238*.ON";
 
      const char *mcfile  = "/data/MAGIC/mc_eth/magLQE_3/gh/0/0/G_M0_00_0_550*.root";
@@ -191 +191 @@
       // path for input for Mars
       //TString inPath = "/.magic/magicserv01/scratch/";
-     TString inPath = "/mnt/data17a/hbartko/";
-     //TString inPath = "~wittek/datacrab_feb04/";
+      //TString inPath = "/mnt/data17a/hbartko/";
+      //TString inPath = "~wittek/datacrab_feb04/";
+      //TString inPath = "~wittek/datacrab_26feb04/";
+     TString inPath = "/.magic/magicserv01/scratch/calibrated/";
+     //TString inPath = "/.magic/magicserv01/scratch/David/CalibratedRuns/";
 
       // path for output from Mars
-      TString outPath = "~wittek/datacrab_feb04/";
+      //TString outPath = "~wittek/datacrab_feb04/";
+      TString outPath = "~wittek/datacrab_26feb04/";
 
       //-----------------------------------------------
@@ -371 +375 @@
     outNameImage += "Hillas";
     outNameImage += typeInput;
-    outNameImage += "1e.root";
+    outNameImage += "1.root";
     gLog << "padded data to be written onto : " << outNameImage << endl;
 
@@ -701 +705 @@
     sourcefromstar.SetSourceAndStarPosition("Crab", 22, 0, 52, 5, 34, 32,
                                         "Zeta-Tau", 21, 8, 33, 5, 37, 38.7);
-    sourcefromstar.AddFile("~wittek/datacrab_feb04/positions.4.txt", 0);
-
-    MBlindPixelCalc blindbeforepad;
+    if (typeInput == "ON")
+      sourcefromstar.AddFile("~wittek/datacrab_26feb04/positionsOn.4.txt", 0);
+    else if (typeInput == "OFF")
+      sourcefromstar.AddFile("~wittek/datacrab_26feb04/positionsOff.txt", 0);
+
+    //MBlindPixelCalc blindbeforepad;
     //blindbeforepad.SetUseBlindPixels();
-    blindbeforepad.SetName("BlindBeforePadding");
-
-    MBlindPixelCalc blind;
+    //blindbeforepad.SetName("BlindBeforePadding");
+
+    //MBlindPixelCalc blind;
     //blind.SetUseBlindPixels();
+    //blind.SetUseInterpolation();
+    //blind.SetName("BlindAfterPadding");
+
+    MSigmabarCalc sigbar;
+
+    MBadPixelCalcRms blind;
+    //blind.SetUseBlindPixels();
+    blind.SetUseInterpolation();
+    //blind.SetCheckPedestalRMS();
     blind.SetName("BlindAfterPadding");
 
@@ -715 +731 @@
     contbasic.SetName("SelBasic");
 
-    //MFillH fillblind("BlindPixels[MHBlindPixels]", "MBlindPixels");
-    //fillblind.SetName("HBlind");
+    MFillH fillblind("BlindPixels[MHBlindPixels]", "MBlindPixels");
+    fillblind.SetName("HBlind");
 
     MSigmabarCalc sigbarcalc;
@@ -759 +775 @@
     selstandard.SetHillasName(fHilName);
     selstandard.SetImgParName(fImgParName);
-    selstandard.SetCuts(200, 6, 600, 0.4, 1.1, 0.0, 0.0);
+    selstandard.SetCuts(200, 4, 600, 0.2, 1.3, 0.0, 0.0);
     MContinue contstandard(&selstandard);
     contstandard.SetName("SelStandard");
@@ -766 +782 @@
       MWriteRootFile write(outNameImage);
 
-      //write.AddContainer("MRawRunHeader", "RunHeaders");
+      write.AddContainer("MRawRunHeader", "RunHeaders");
       //write.AddContainer("MMcRunHeader",  "RunHeaders", kFALSE);
       //write.AddContainer("MTime",         "Events");
@@ -803 +819 @@
     //  tliston.AddToList(&pointcorr);
 
+    tliston.AddToList(&sigbar);
     tliston.AddToList(&blind);
     tliston.AddToList(&contbasic);
 
-    //tliston.AddToList(&fillblind);
-    tliston.AddToList(&sigbarcalc);
+    tliston.AddToList(&fillblind);
+    //tliston.AddToList(&sigbarcalc);
     tliston.AddToList(&fillsigtheta);
     tliston.AddToList(&clean);
@@ -847 +864 @@
 
     pliston.FindObject("SigmaTheta", "MHSigmaTheta")->DrawClone();
-    //pliston.FindObject("BlindPixels", "MHBlindPixels")->DrawClone();
+    pliston.FindObject("BlindPixels", "MHBlindPixels")->DrawClone();
 
     pliston.FindObject("MHHillas")->DrawClone();
@@ -855 +872 @@
     pliston.FindObject("MHStarMap")->DrawClone();
 
-    DeleteBinnings(&pliston);
+    //DeleteBinnings(&pliston);
 
       gLog << "End of padding" << endl;