Changeset 2543 for trunk/MagicSoft/Mars/macros/CT1Analysis.C

Timestamp:

11/20/03 09:13:31 (21 years ago)

Author:

wittek

Message:

*** empty log message ***

File:

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

trunk/MagicSoft/Mars/macros/CT1Analysis.C

@@ -286 +286 @@
 
 
+
     // Job E_XX : extended version of E_XX (including flux plots)  
     //  - select g/h separation method XX
@@ -297 +298 @@
 
 
-    Bool_t JobE_XX  = kFALSE;  
+    Bool_t JobE_XX  = kTRUE;  
+    Bool_t CCollArea= kFALSE;  // calculate eff. collection areas
     Bool_t OEEst    = kFALSE;  // optimize energy estimator
-    Bool_t WEX      = kFALSE;  // update root file  ?
-    Bool_t WRobert  = kFALSE;  // write out Robert's file  ?
+    Bool_t WEX      = kTRUE;  // update root file  ?
+    Bool_t WRobert  = kTRUE;  // write out Robert's file  ?
+
 
 
@@ -2527 +2530 @@
 
 
-
   //---------------------------------------------------------------------
   // Job E_XX
@@ -2547 +2549 @@
 
     gLog << "" << endl;
-    gLog << "Macro CT1Analysis : JobE_XX, OEEst, WEX = " 
+    gLog << "Macro CT1Analysis : JobE_XX, CCollArea, OEEst, WEX = " 
          << (JobE_XX ? "kTRUE" : "kFALSE")  << ",  " 
+         << (CCollArea?"kTRUE" : "kFALSE")  << ",  " 
          << (OEEst ?   "kTRUE" : "kFALSE")  << ",  " 
          << (WEX     ? "kTRUE" : "kFALSE")  << endl;
@@ -2554 +2557 @@
 
     // type of data to be analysed
-    TString typeData = "ON";
+    //TString typeData = "ON";
     //TString typeData = "OFF";
-    //TString typeData = "MC";
+    TString typeData = "MC";
     gLog << "typeData = " << typeData << endl;
 
@@ -2567 +2570 @@
     // and definition of final selections
 
-    TString XX("SC");
-    //TString XX("RF");
+    //TString XX("SC");
+    TString XX("RF");
     TString fhadronnessName("Had");
     fhadronnessName += XX;
@@ -2574 +2577 @@
 
     // maximum values of the hadronness, |ALPHA| and DIST
-    Float_t maxhadronness   = 0.4;
+    Float_t maxhadronness   = 0.23;
     Float_t maxalpha        = 20.0;
     Float_t maxdist         = 10.0;
@@ -2638 +2641 @@
 
     //====================================================================
+
+    MHMcCT1CollectionArea collarea;
+    collarea.SetEaxis(MHMcCT1CollectionArea::kLinear);
+
+    MParList  parlist;
+    InitBinnings(&parlist);
+
+  if (CCollArea)
+  {
     gLog << "-----------------------------------------------" << endl;
     gLog << "Start calculation of effective collection areas" << endl;
-    MParList  parlist;
+
+
     MTaskList tasklist;
 
@@ -2655 +2668 @@
     MContinue conthadrons(&cuthadrons);
 
-    MHMcCT1CollectionArea collarea;
-    collarea.SetEaxis(MHMcCT1CollectionArea::kLinear);
 
     MFillH filler("MHMcCT1CollectionArea", "MMcEvt");
@@ -2665 +2676 @@
 
     parlist.AddToList(&tasklist);
-    InitBinnings(&parlist);
+
     parlist.AddToList(&collarea);
 
@@ -2697 +2708 @@
     collarea.CalcEfficiency();
     collarea.DrawClone();
+
+
+
+    //---------------------------------------------
+    // Write histograms to a file 
+    //
+
+    TFile f(collareaName, "RECREATE");
+    //collarea.GetHist()->Write();
+    //collarea.GetHAll()->Write();
+    //collarea.GetHSel()->Write();
+    collarea.Write();
+
+    f.Close();
+
+    gLog << "Collection area plots written onto file " << collareaName << endl;
+
+    gLog << "Calculation of effective collection areas done" << endl;
+    gLog << "-----------------------------------------------" << endl;    
+    //------------------------------------------------------------------
+  }
+
+  if (!CCollArea)
+  {
+    gLog << "-----------------------------------------------" << endl;
+    gLog << "Read in effective collection areas from file " 
+         << collareaName << endl;
+
+    TFile collfile(collareaName);
+    collfile.ls();
+    collarea.Read("MHMcCT1CollectionArea");
+    collarea.DrawClone();
+
+    gLog << "Effective collection areas were read in from file " 
+         << collareaName << endl;
+    gLog << "-----------------------------------------------" << endl;    
+  }
+
 
     // save binnings for call to CT1EEst
@@ -2712 +2761 @@
         }
 
-
-    //---------------------------------------------
-    // Write histograms to a file 
-    //
-
-    TFile f(collareaName, "RECREATE");
-    collarea.GetHist()->Write();
-    collarea.GetHAll()->Write();
-    collarea.GetHSel()->Write();
-    f.Close();
-
-    gLog << "Collection area plots written onto file " << collareaName << endl;
-
-    gLog << "Calculation of effective collection areas done" << endl;
-    gLog << "-----------------------------------------------" << endl;    
-    //------------------------------------------------------------------
-
-
+    //-------------------------------------
     TString fHilName    = "MHillas"; 
     TString fHilNameExt = "MHillasExt";