Changeset 16671 for trunk/Mars/fact

Timestamp:

06/03/13 22:20:46 (11 years ago)

Author:

Daniela Dorner

Message:

implemented fOnTimeAfterCuts, made macro compile,fixed some small bugs

File:

• trunk/Mars/fact/processing/numevents.C (modified) (9 diffs)

trunk/Mars/fact/processing/numevents.C

@@ +1 @@
+#include <iostream>
+
+#include <TCanvas.h>
+#include <TString.h>
+#include <TTree.h>
+#include <TFile.h>
+#include <TSystem.h>
+#include <TSQLResult.h>
+#include <TSQLRow.h>
+
+#include "MSQLMagic.h"
+#include "MAlphaFitter.h"
+#include "MHThetaSq.h"
+#include "MStatusArray.h"
+
+
 using namespace std;
 
@@ -16 +32 @@
         return 0;
     }
+    serv.SetIsDummy(dummy);
 
     //select runs, where star finished, format the list nicely
@@ -40 +57 @@
 
     //allocate variables to use in the loop
-    TString night;
+    TString night2;
     TString runid;
     TString run;
@@ -59 +76 @@
     Float_t fNumExcEvts         = 0;
     Float_t fNumIslandsMean     = 0;
+    Float_t fOnTimeAfterCuts    = 0;
 
     //loop over the data files
     TSQLRow *row=0;
-    Int_t counter=0;
 
     /*
@@ -77 +94 @@
     {
         //use the results from the query
-        night=(*row)[4];
+        night2=(*row)[4];
         runid=(*row)[5];
         run=(*row)[0];
@@ -88 +105 @@
             ganymed_fname=inpath+"/ganymed_run/"+year+"/"+month+"/"+day+"/"+run+"-analysis.root";
         else
-            ganymed_fname=inpath+"/ganymed_night/"+night+"-analysis.root";
+            ganymed_fname=inpath+"/ganymed_night/"+night2+"-analysis.root";
         //cout << "gf => " << ganymed_fname << endl;
 
@@ -150 +167 @@
             continue;
         }
+
+        TH1D *vstime = (TH1D*)arr.FindObjectInCanvas("OnTime",  "TH1D", "OnTime");
+        if (!vstime)
+        {
+            cout << "WARNING - Reading of Theta failed." << endl;
+            continue;
+        }
+        fOnTimeAfterCuts = vstime->Integral();
 
         //Get number of events after quality cuts
@@ -204 +229 @@
 
         //inserting or updating the information in the database
-        vars = Form("fRunID=%s, fNight=%s,", runid.Data(), night.Data());
+        vars = Form("fRunID=%s, fNight=%s,", runid.Data(), night2.Data());
         vars += Form(" fNumEvtsAfterQualCuts=%d, "
                     " fNumEvtsAfterBgCuts=%d, "
                     " fNumBgEvts=%6.1f, "
                     " fNumSigEvts=%6.1f, "
-                    " fNumExcEvts=%6.1f ",
+                    " fNumExcEvts=%6.1f, "
+                    " fOnTimeAfterCuts=%7.2f ",
                     fNumEvtsAfterQualCuts,
                     fNumEvtsAfterBgCuts,
                     fNumBgEvts,
                     fNumSigEvts,
-                    fNumExcEvts
+                    fNumExcEvts,
+                    fOnTimeAfterCuts
                    );
         if (!pernight)
-            vars += Form(", fNumIslandsMean=%6.2f ",
-                         ", fNumEvtsAfterCleaning=%d ",
+            vars += Form(", fNumIslandsMean=%6.2f, fNumEvtsAfterCleaning=%d ",
                          fNumIslandsMean,
                          fNumEvtsAfterCleaning
@@ -225 +251 @@
 
         if (pernight)
-            where = Form("fNight=%s", night.Data());
+            where = Form("fNight=%s", night2.Data());
         else
-            where = Form("fRunID=%s AND fNight=%s", runid.Data(), night.Data());
+            where = Form("fRunID=%s AND fNight=%s", runid.Data(), night2.Data());
 
         cout << "vars: " << vars << endl;