Changeset 4992 for trunk/MagicSoft/Mars/mtemp/mpisa/macros

Timestamp:

09/14/04 12:36:12 (20 years ago)

Author:

stamerra

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mtemp/mpisa/macros/readIPR.C (modified) (4 diffs)

trunk/MagicSoft/Mars/mtemp/mpisa/macros/readIPR.C

@@ -45 +45 @@
     plist.AddToList(&geom);
 
-    // First Task: Read ROOT file with REPORT data
-    MReportFileRead read(fname);
-    read.SetHasNoHeader();
-    // The Trigger container is added to the list
-    read.AddToList("MReportTrigger");
+    // First Task: Read ROOT file with Trigger-REPORT data
+    MReadTree read("Trigger", fname);
+    read.DisableAutoScheme();
+
     tlist.AddToList(&read);
 
+    // Create the container for the IPRs
+    MTriggerIPR ipr;
+    plist.AddToList(&ipr);
 
-    MHCamEvent IPRcam("IPRcam","IPRs");
-    plist.AddToList(&IPRcam);
-    // create a task to fill a histogram
-    //    MFillH fill("MHCamEvent", "MTriggerIPR");
-
-    MFillH fillIPR(&IPRcam, "MTriggerIPR");
+    // Create the histo to display the IPRs
+    MHCamEvent IPRhist("IPRhist","IPRs");
+    plist.AddToList(&IPRhist);
+        
+    // create a task to fill a histogram from the container
+    MFillH fillIPR(&IPRhist, "MTriggerIPR");
     tlist.AddToList(&fillIPR);
 
@@ -68 +70 @@
 
     //
-    // Execute your analysis
+    // Execute the analysis
     //
+    if (!evtloop.PreProcess())
+      return;
+    
     if (!evtloop.Eventloop())
-        return;
+      return;
 
     tlist.PrintStatistics();
 
-    //
+    //-----------------------------------
     // Now display the result of the loop
     //
 
-    //    MHCamEvent &h2 = *(MHCamEvent*)plist->FindObject("MHCamEvent");
-    //    MHCamera &h = *(MHCamera*)h2.GetHistByName("sum");
-    MHCamera &h = *(MHCamera*)IPRcam.GetHistByName("sum");
+    // create a MHCamera histo where the sum off all events is filled
+    MHCamera &h = *(MHCamera*)IPRhist.GetHistByName("sum");
 
     TCanvas *c = MH::MakeDefCanvas();
@@ -92 +96 @@
     disp3->SetCamContent(h, 2);
 
-    disp1->SetYTitle("I [nA]");
-    disp2->SetYTitle("\\sigma_{I} [\\mu A]");
-    disp3->SetYTitle("\\sigma_{I} [%]");
+    disp1->SetYTitle("Rate [Hz]");
+    disp2->SetYTitle("\\sigma_{Rate} [Hz]");
+    disp3->SetYTitle("\\sigma_{Rate} [%]");
     disp1->SetName("IPRs;avg");
     disp2->SetName("IPRs;err");
@@ -139 +143 @@
 
     c->SaveAs(fname(0, fname.Last('.')+1) + "ps");
-    c->SaveAs(fname(0, fname.Last('.')+1) + "root");
+    //c->SaveAs(fname(0, fname.Last('.')+1) + "root");
 }