Changeset 4868 for trunk/MagicSoft/Mars/macros/ONOFFAnalysis.C

Timestamp:

09/06/04 08:31:05 (20 years ago)

Author:

wittek

Message:

*** empty log message ***

File:

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

trunk/MagicSoft/Mars/macros/ONOFFAnalysis.C

@@ -172 +172 @@
       //TString tag = "040126";
       //TString tag = "040127";
-      TString tag = "040215";
+      //TString tag = "040215";
+      TString tag = "040422";
 
       //const char *offfile = "~magican/ct1test/wittek/offdata.preproc"; 
@@ -184 +185 @@
       // 15 Feb 04
       //const char *offfile  = "*.OFF";
-      const char *offfile  = "174*.OFF";
-
+      //const char *offfile  = "174*.OFF";
+      const char *offfile  = "*.OffMrk421*";
 
       //const char *onfile  = "~magican/ct1test/wittek/mkn421_on.preproc"; 
@@ -199 +200 @@
       // 15 Feb 04
       //const char *onfile  = "*.ON";
-      const char *onfile  = "174*.ON";
+      //const char *onfile  = "174*.ON";
+      const char *onfile  = "*.Mrk421*";
 
 
@@ -231 +233 @@
      //if (tag == "040215")
      //  TString inPath = "/mnt/magicserv01/scratch/David/CalibratedData/Crab/2004_02_15/";
-     if (tag == "040215")   
+     if (tag == "040215")  
+     { 
        TString inPathON = "/.magic/magicserv01/scratch/Daniel/CalibratedData/MispointingTest/2004_02_15ForgottenData/finesam20/";
        TString inPathOFF = "/.magic/magicserv01/scratch/David/CalibratedData/Mkn421/2004_02_15/";
        TString inPathMC = "/.magic/magicserv01/scratch/David/MCData/MCApril2004/Data/gammas_highnoise/";
+     }
+
+     if (tag == "040422")  
+     { 
+       TString inPathON  = "/.magic/data03a/mazin/calibrdata/2004_04_22/";
+       TString inPathOFF = "/.magic/data03a/mazin/calibrdata/2004_04_22/";
+       TString inPathMC  = "/.magic/magicserv01/scratch/David/MCData/MCApril2004/Data/gammas_highnoise/";
+     }
 
       // path for output from Mars
       //TString outPath = "~wittek/datacrab_feb04/";
      //TString outPath = "~wittek/datacrab_01march04/";
-      TString outPath = "~wittek/datacrab_27april04/";
+     // TString outPath = "~wittek/datacrab_27april04/";
+     TString outPath = "/mnt/data03a/wittek/";
+     outPath += tag;
+     outPath += "/";
 
       //-----------------------------------------------
@@ -256 +270 @@
 
     Bool_t JobA    = kTRUE;  
-    Bool_t GPadON  = kFALSE;    // \  generate padding histograms 
+    Bool_t GPadON  = kTRUE;    // \  generate padding histograms 
     Bool_t GPadOFF = kFALSE;    //  | and write them onto disk
     Bool_t GPadMC  = kFALSE;    // /
     Bool_t Merge   = kFALSE;   // merge padding histograms
                                // and write them onto disk
-    Bool_t Wout    = kTRUE;   // read in merged padding histograms and
+    Bool_t Wout    = kFALSE;   // read in merged padding histograms and
                               // write out root file of padded data
                               // (ON1.root, OFF1.root or MC1.root) 
@@ -477 +491 @@
     MGeomApply        apply;
 
-    MSourcePosfromStarPos sourcefromstaron;
-    sourcefromstaron.AddFile("~wittek/datacrab_26feb04/positions040215.txt", 0);
-
-      MBlindPixelsCalc2 blindon;
-      //blindon.SetUseBlindPixels();
-      blindon.SetCheckPedestalRms();
+    TString pedphotcamname("MPedPhotCamFromData");
+
+    MBadPixelsCalc badcalcon;
+    badcalcon.SetNamePedPhotContainer(pedphotcamname);
+    badcalcon.SetPedestalLevel(2.0);
+
+    MBadPixelsTreat badtreaton;
+    badtreaton.SetNamePedPhotCam(pedphotcamname);
+    badtreaton.SetUseInterpolation();
+    badtreaton.SetProcessPedestalEvt();
+    badtreaton.SetProcessTimes();
 
       MFSelBasic selbasicon;
       MContinue contbasicon(&selbasicon);
 
-      MHBlindPixels blindON("BlindPixelsON");
-      MFillH fillblindON("BlindPixelsON[MHBlindPixels]", "MBlindPixels");
-      fillblindON.SetName("FillBlind");
+      MHBadPixels badON("BadPixelsON");
+      MFillH fillbadON("BadPixelsON[MHBadPixels]", "MBadPixelsCam");
+      fillbadON.SetName("FillBad");
 
       MSigmabarCalc sigbarcalcon;
 
       MHSigmaTheta sigthON("SigmaThetaON");
-      MFillH fillsigthetaON ("SigmaThetaON[MHSigmaTheta]", "MPointingPos");
+      sigthON.SetNamePedPhotCam(pedphotcamname);
+      MFillH fillsigthetaON ("SigmaThetaON[MHSigmaTheta]", "");
       fillsigthetaON.SetName("FillSigTheta");    
  
@@ -503 +523 @@
       InitBinnings(&pliston);
       pliston.AddToList(&observon);
-      pliston.AddToList(&blindON);
+      pliston.AddToList(&badON);
       pliston.AddToList(&sigthON);
 
@@ -512 +532 @@
       tliston.AddToList(&readON);
       tliston.AddToList(&apply);
-      tliston.AddToList(&sourcefromstaron);
-
-      tliston.AddToList(&blindon);
+
+      tliston.AddToList(&badcalcon);
+      tliston.AddToList(&badtreaton);
 
       tliston.AddToList(&contbasicon);
-      tliston.AddToList(&fillblindON);
+      tliston.AddToList(&fillbadON);
       tliston.AddToList(&sigbarcalcon);
       tliston.AddToList(&fillsigthetaON);
@@ -533 +553 @@
       tliston.PrintStatistics(0, kTRUE);
 
-      blindON.DrawClone();
+      badON.DrawClone();
       sigthON.DrawClone();
 
@@ -541 +561 @@
       TH3D *diffpixthon = sigthON.GetDiffPixTheta();
 
-      TH2D *blindidthon = blindON.GetBlindId();
-      TH2D *blindnthon  = blindON.GetBlindN();
+      TH2D *badidthon = badON.GetBadId();
+      TH2D *badnthon  = baddON.GetBadN();
 
       gLog << "" << endl;
@@ -552 +572 @@
       sigpixthon->Write();
       diffpixthon->Write();
-      blindidthon->Write();
-      blindnthon->Write();
+      badidthon->Write();
+      badnthon->Write();
       filewriteon.Close();