Changeset 988 for trunk/MagicSoft/Mars/macros/threshold.C

Timestamp:

10/24/01 15:12:21 (23 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/threshold.C (modified) (7 diffs)

trunk/MagicSoft/Mars/macros/threshold.C

@@ -24 +24 @@
 
 
-void threshold(char* filename="data/CrabNebula_dnsb_09_loop.root")
+void threshold(char* filename="data/camera.root")
 {
     //
@@ -30 +30 @@
     // MMcThresholdCalc and shows the results.
     //
-    MParList  parlist;
+    MParList parlist;
 
     MTaskList tasklist;
@@ -43 +43 @@
     //      taking only the trigger information from MMcTrig
     //
-    const UInt_t numtriggerconditions = 4;
+    const UInt_t numtrigcond = 0;
+
+    UInt_t from = numtrigcond>0 ?           1 : -numtrigcond;
+    UInt_t to   = numtrigcond>0 ? numtrigcond : -numtrigcond;
+
+    Int_t dim = to-from+1;
 
     //
@@ -49 +54 @@
     // and store the histograms in an TObjArray
     //
-    TObjArray *hists = new TObjArray(MParList::CreateObjList("MHMcEnergy", numtriggerconditions));
+    TObjArray hists(MParList::CreateObjList("MHMcEnergy", from, to));
+    hists.SetOwner();
 
     //
     // Check if the list really contains the right number of histograms
     //
-    if (hists->GetEntriesFast() != numtriggerconditions)
+    if (hists.GetEntriesFast() != dim)
         return;
 
@@ -60 +66 @@
     // Add the histograms to the paramater list.
     //
-    parlist.AddToList(hists);
+    parlist.AddToList(&hists);
 
     //
@@ -71 +77 @@
     //      like one dimension MMcThresholdCalc
     //
-    MReadTree        read("Events;7", filename);
-    MMcThresholdCalc calc(numtriggerconditions);
+    MReadTree read("Events", filename);
+    read.UseLeaf("fEnergy");
+    read.UseLeaf("fNumFirstLevel");
+
+    MMcThresholdCalc calc(numtrigcond);
 
     tasklist.AddToList(&read);
@@ -90 +99 @@
     // Now you can display the results
     //
-    for (UInt_t i=0; i<numtriggerconditions; i++)
-        ((*hists)[i])->Draw();
+    TIter Next(&hists);
+    TObject *obj;
+    while ((obj=Next()))
+        obj->DrawClone();
 }