Changeset 1837 for trunk/MagicSoft/Mars

Timestamp:

03/19/03 12:41:28 (22 years ago)

Author:

moralejo

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (2 diffs)
• macros/CT1collarea.C (modified) (3 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -1 +1 @@
                                                  -*-*- END OF LINE -*-*-
+
+ 2003/03/19: Abelardo Moralejo
+
+    * macros/CT1collarea.C:
+      - Added filter to cut in hadronness (now available in new version of 
+        root CT1 Monte Carlo file from Wolfgang).
 
  2003/03/19: Robert Wagner
@@ -6 +12 @@
       - changed MH3::New() such that constructor according to Thomas'
         previous changes is supported
-
 
 

trunk/MagicSoft/Mars/macros/CT1collarea.C

@@ -24 +24 @@
 
 
-void CT1collarea(TString filename="MC_ON1.root", TString outname="")
+void CT1collarea(TString filename="MC_ON2.root", TString outname="")
 { 
     //
@@ -40 +40 @@
     //
     MReadMarsFile reader("Events", filename);
+    reader.DisableAutoScheme();
+    //    reader.EnableBranch("fHadronness");
+
     MMcCT1CollectionAreaCalc effi;
+    MHMcCT1CollectionArea* collarea = new MHMcCT1CollectionArea("MHMcCT1CollectionArea","",30,2.,5.);
+    parlist.AddToList(collarea);
 
     tasklist.AddToList(&reader);   
+
+    MF filterhadrons("MHadronness.fHadronness<0.25");
+    tasklist.AddToList(&filterhadrons);
+    effi.SetFilter(&filterhadrons);
 
     tasklist.AddToList(&effi);
@@ -75 +84 @@
         return;
 
-    TFile f(outname, "recreate");
-    if (!f)
-        return;
-    parlist.FindObject("MHMcCT1CollectionArea")->Write();
+    TFile f("area.root","recreate");
+    collarea->GetHist()->Write();
+    collarea->GetHAll()->Write();
+    collarea->GetHSel()->Write();
+    f.Close();
 }