Changeset 1376 for trunk/MagicSoft/Mars/macros

Timestamp:

07/04/02 18:40:55 (22 years ago)

Author:

bigongia

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/macros/trigrate.C (modified) (2 diffs)

trunk/MagicSoft/Mars/macros/trigrate.C

@@ -20 +20 @@
 !   Copyright: MAGIC Software Development, 2000-2001
 !
+!   Modified 4/7/2002, Abelardo Moralejo:
+!   Added one optional input parameter: a camera .root file containing
+!   pure NSB events. One such file is generated running the camera over an
+!   "empty" reflector file, with the NSB option on, and MINPHOT=0. The nsb
+!   files must contain the same trigger conditions as the proton file.
+!   If no nsb file is supplied, the macro will assume no triggers from 
+!   pure NSB fluctuations.
 !
 \* ======================================================================== */
 
 
-void trigrate(int dim=0, char *filename = "data/camera.root" )
+void trigrate(int dim=0, char *filename = "data/camera.root", 
+              char *nsbfile = NULL)
 {
-    // This macro has two input parameter:
+    // The dim parameter has three possible values:
     // dim : = 0 -> root file with 1 trigger condition.
     //       > 0 -> number of trigger condition to be analised 
@@ -90 +98 @@
     tasklist.AddToList(&reader);
 
-    Float_t BgR[10]={0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
+ 
+    // Now we have to build the BgR array, containing the number 
+    // of triggers (may be more than 1 trigger/event!) from the
+    // total_nsb_events simulated in the nsbfile (3rd input parameter).
+    // If no nsbfile is supplied, we assume no triggers from NSB
+
+    Float_t total_nsb_events = 5.e4; // some default value.
+
+    Float_t* BgR = new float[num];
+    for(Int_t i = 0; i < num; i++)
+      BgR[i] = 0.;
+
+    //    Float_t* BgR = new float[36];
+    //    for(Int_t i = 0; i < 36; i++)
+    //      BgR[i] = 0.;
+
+    if (nsbfile)
+      {
+        nsb_events = new TChain("Events");
+        nsb_events->Add(nsbfile);
+
+        //      nsb_events = new TFile(nsbfile);
+        //      nsb_events->cd();
+
+        h = new TH1F("h","",5,.5,5.5);
+
+        Char_t plot[256];
+        for (Int_t i = from; i <= to; i++)
+          {
+            // For each trigger condition (i)...
+            // histogram Number of L1 triggers:
+            sprintf(plot, "MMcTrig;%d.fNumFirstLevel>>h", i);
+            nsb_events->Draw(plot,"", "goff");
+            // Get total number of L1 triggers from histogram:
+            if (dim > 0)
+              BgR[i-1] = h->Integral()*h->GetMean();
+            else
+              BgR[0] = h->Integral()*h->GetMean();
+
+            total_nsb_events = (Float_t) h->GetEntries();
+          }
+        h->Delete();
+        nsb_events->Delete();
+        //      nsb_headers = new TChain("RunHeaders");
+        //      nsb_headers->Add(nsbfile);
+      }
+
     cout << "Number of Trigger conditions: " << num << endl;
 
-    MMcTriggerRateCalc rate(dim, 14, BgR, 100000);
+    MMcTriggerRateCalc rate(dim, kPROTON, BgR, total_nsb_events);
     tasklist.AddToList(&rate);