Changeset 2489 for trunk/MagicSoft/Mars/macros

Timestamp:

11/08/03 09:20:10 (21 years ago)

Author:

galante

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/triglvl2.C (modified) (6 diffs)

trunk/MagicSoft/Mars/macros/triglvl2.C

@@ -38 +38 @@
    
     tasklist.AddToList(&reader);
-
     MGeomCamMagic geocam;
     parlist.AddToList(&geocam);
@@ -55 +54 @@
     // Set the number of Next Neighbourhoods that define a compact pixel
     //
-    cell.SetCompactNN(2);
+    cell.SetCompactNN(3);
     //  USER:       --^
 
@@ -61 +60 @@
     //  A filter to select events using the L2T parameters is created
     //
-
     // MF lvl2filter("MMcTriggerLvl2.fPseudoSize > 25 && MMcTriggerLvl2.fPseudoSize < 31");
-    MF lvl2filter("MMcTriggerLvl2.fPseudoSize > 15");
+    //    MF lvl2filter("MMcTriggerLvl2.fSizeBiggerCell > 34");
+    MF lvl2filter("MMcTriggerLvl2.fSizeBiggerCell > 18");
  
     //
@@ -90 +89 @@
     MFilterList flist;
     //flist.AddToList(&energyfilter);
-    //flist.AddToList(&lvl2filter);
+    flist.AddToList(&lvl2filter);
     //   flist.AddToList(&selector);
     //   flist.AddToList(&eslope);
     //   flist.AddToList(&fTrig);
 
-
-    
     //
     // The task to calculate the L2T parameter is added to the task list
@@ -107 +104 @@
     //
     tasklist.AddToList(&flist);
+
+
+    MContinue L1offline ("MMcTriggerLvl2.fLutPseudoSize < 6");
+    //    tasklist.AddToList(&L1offline);
 
     //
@@ -176 +177 @@
     // filled and can be displayd
     //
-    parlist.FindObject("MHMcCollectionArea")->DrawClone();
-
-    parlist.FindObject("MHMcTriggerLvl2")->DrawClone("sbc");
-    parlist.FindObject("MHMcTriggerLvl2")->DrawClone("lps");
-    parlist.FindObject("MHMcTriggerLvl2")->DrawClone();
+    
+
+    //parlist.FindObject("MHMcTriggerLvl2")->DrawClone("sbc");
+    //parlist.FindObject("MHMcTriggerLvl2")->DrawClone("lps");
+    //parlist.FindObject("MHMcTriggerLvl2")->DrawClone();
     //parlist.FindObject("MHMcTriggerLvl2")->DrawClone("energy");
-    parlist.FindObject("MHMcTriggerLvl2")->DrawClone("size-energy");
-
+    //parlist.FindObject("MHMcTriggerLvl2")->DrawClone("lut-energy");
+    //parlist.FindObject("MHMcCollectionArea")->DrawClone();
     // Returns histogram of the class MHMcTriggerLvl2  
-
+    
     MHMcTriggerLvl2 *htrig = (MHMcTriggerLvl2 *)parlist.FindObject("MHMcTriggerLvl2");
-    TH1F *h1 = (TH1F *)htrig->GetHistByName("fHistPseudoSize");
-    TH2D *h2 = htrig->GetHistSizeBiggerCellEnergy();
-    TH2D *h3 = (TH2D *)htrig->GetHistByName("fHistLutPseudoSizeEnergy");
-    //csbce = new TCanvas();
-    //h2->DrawClone();
-
-    hfile = new TFile("HistFileLUT.root", "RECREATE");
-    h1->Write();
-    h2->Write();
-    h3->Write();
+    MHMcCollectionArea *collarea = (MHMcCollectionArea *)parlist.FindObject("MHMcCollectionArea");
+    TH1F *hps = (TH1F *)htrig->GetHistByName("fHistPseudoSize");
+    TH1F *hlps = (TH1F *)htrig->GetHistByName("fHistLutPseudoSize");
+    TH1F *hsbc = (TH1F *)htrig->GetHistSizeBiggerCell();
+    TH1F *hpsn = (TH1F *)htrig->GetHistByName("fHistPseudoSizeNorm");
+    TH1F *hlpsn = (TH1F *)htrig->GetHistByName("fHistLutPseudoSizeNorm");
+    TH1F *hsbcn = (TH1F *)htrig->GetHistSizeBiggerCellNorm();
+    TH2D *hlpse = (TH2D *)htrig->GetHistByName("fHistLutPseudoSizeEnergy");
+    TH2D *hpse = (TH2D *)htrig->GetHistByName("fHistPseudoSizeEnergy");
+    TH2D *hsbce = (TH2D *)htrig->GetHistSizeBiggerCellEnergy();
+    TH1D *hcollarea = collarea->GetHist();
+
+    // Please set rootfile as the directory where you want to store
+    // your file.root containing histograms
+    rootfile = new TString(filename);
+    if (rootfile->Index("/*.root",1)>=0){
+      rootfile->Resize(rootfile->Index("/*.root",1));
+      rootfile[0]+="_sbcmag18_CNN3.root";
+    }
+      
+    hfile = new TFile((char *)rootfile[0], "RECREATE");
+    hlps->Write();
+    hps->Write();
+    hsbc->Write();
+    hlpsn->Write();
+    hpsn->Write();
+    hsbcn->Write();
+    hlpse->Write();
+    hpse->Write();
+    hsbce->Write();
+    hcollarea->Write();
     hfile->Close();
+    cout << "Histograms stored into " << rootfile[0] << endl;
+    
 
 }