Changeset 9062 for trunk/MagicSoft/Mars/datacenter/macros/fillcalib.C

Timestamp:

07/31/08 17:41:29 (16 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/datacenter/macros/fillcalib.C (modified) (10 diffs)

trunk/MagicSoft/Mars/datacenter/macros/fillcalib.C

@@ -61 +61 @@
 #include <iostream>
 
-#include <TEnv.h>
-#include <TRegexp.h>
-
 #include <TH1.h>
 
@@ -74 +71 @@
 #include "MStatusArray.h"
 #include "MHCamera.h"
+#include "MSequence.h"
 #include "MGeomCamMagic.h"
 #include "MBadPixelsCam.h"
@@ -145 +143 @@
     if (!h)
     {
-        cout << "WARNING - Could not find histogram HRelTimeHiGainArea0." << endl;
+        cout << "ERROR - Could not find histogram HRelTimeHiGainArea0 in Time." << endl;
         return 2;
     }
@@ -156 +154 @@
     if (!h)
     {
-        cout << "WARNING - Could not find histogram HRelTimeHiGainArea1." << endl;
+        cout << "ERROR - Could not find histogram HRelTimeHiGainArea1 in Time." << endl;
         return 2;
     }
@@ -167 +165 @@
     if (!c)
     {
-        cout << "WARNING - Could not find MHCamera TotalConv." << endl;
+        cout << "ERROR - Could not find MHCamera TotalConv in Conversion." << endl;
         return 2;
     }
@@ -182 +180 @@
     TString meanconvouter=Form("%6.3f", meanconvo);
 
-    //Getting sequ# from filename
-    TString sequence = fname(TRegexp("calib[0-9]+[.]root$"));
-    if (sequence.IsNull())
-    {
-        cout << "WARNING - Could get sequence# from filename: " << fname << endl;
-        return 2;
-    }
-
-    Int_t seq = atoi(sequence.Data()+5);
+    //Getting relative charge rms
+    c = (MHCamera*)arr.FindObjectInCanvas("RMSperMean", "MHCamera", "FitCharge");
+    if (!c)
+    {
+        cout << "ERROR - Could not find MHCamera RMSperMean in FitCharge." << endl;
+        return 2;
+    }
+
+    Stat_t relrmsi = c->GetMedianSectors(TArrayI(6, s0), inner);
+    Stat_t relrmso = c->GetMedianSectors(TArrayI(6, s0), outer);
+    TString relrmsinner=Form("%6.3f", relrmsi);
+    TString relrmsouter=Form("%6.3f", relrmso);
+
+    //Getting relative charge rms
+    c = (MHCamera*)arr.FindObjectInCanvas("NumPhes", "MHCamera", "FitCharge");
+    if (!c)
+    {
+        cout << "ERROR - Could not find MHCamera NumPhes in FitCharge." << endl;
+        return 2;
+    }
+
+    Stat_t numphei = c->GetMedianSectors(TArrayI(6, s0), inner);
+    Stat_t numpheo = c->GetMedianSectors(TArrayI(6, s0), outer);
+    TString numpheinner=Form("%5.1f", numphei);
+    TString numpheouter=Form("%5.1f", numpheo);
+
+    MSequence seq;
+    if (seq.Read("sequence[0-9]{8}[.]txt|MSequence")<=0)
+    {
+        cout << "ERROR - Could not find sequence in file: " << fname << endl;
+        return 2;
+    }
+    if (!seq.IsValid())
+    {
+        cout << "ERROR - Sequence read from file inavlid: " << fname << endl;
+        return 2;
+    }
 
     //getting the ratio of calibration events used
@@ -196 +222 @@
     if (!h)
     {
-        cout << "WARNING - Could not find histogram ArrTime;avg." << endl;
-        return 2;
-    }
-
-    UInt_t nevts = h->GetEntries();
+        cout << "ERROR - Could not find histogram ArrTime;avg." << endl;
+        return 2;
+    }
+
+    UInt_t nevts = TMath::Nint(h->GetEntries());
 
     TString query;
     query = Form("SELECT SUM(fNumEvents) FROM RunData "
-                 "LEFT JOIN RunType ON RunType.fRunTypeKEY=RunData.fRunTypeKEY "
-                 "WHERE fSequenceFirst=%d AND RunType.fRunTypeName='Calibration'",
-                 seq);
+                 "LEFT JOIN RunType USING (fRunTypeKEY) "
+                 "WHERE fSequenceFirst=%d AND fTelescopeNumber=%d AND "
+                 "RunType.fRunTypeName='Calibration'",
+                 seq.GetSequence(), seq.GetTelescope());
 
     TSQLResult *res = serv.Query(query);
@@ -228 +255 @@
     delete res;
 
-    cout << "Sequence #" << seq << endl;
-    cout << "  Unsuitable: (i/o) " << Form("%3d %3d", (int)unsin, (int)unsout) << endl; // Unbrauchbar
-    cout << "  Unreliable: (i/o) " << Form("%3d %3d", (int)unrin, (int)unrout) << endl; // Unzuverlaessig
-    cout << "  Isolated:   (i/o) " << Form("%3d %3d", (int)isoin, (int)isoout) << endl; // Isolated (unbrauchbar)
-    cout << "  Max.Cluster:      " << Form("%3d", (int)clumax)                 << endl; // Max Cluster
-    cout << "  Arr Time inner:    " << meaninner << " +- " << rmsinner  << endl;
-    cout << "  Arr Time outer:    " << meanouter << " +- " << rmsouter  << endl;
-    cout << "  Mean Conv inner:   " << meanconvinner << endl;
-    cout << "  Mean Conv outer:   " << meanconvouter << endl;
-    cout << "  Ratio Calib Evts:  " << ratiocal << endl;
+    cout << "Sequence M" << seq.GetTelescope() << ":" << seq.GetSequence() << endl;
+    cout << "  Unsuitable: (i/o)  " << Form("%3d %3d", (int)unsin, (int)unsout) << endl; // Unbrauchbar
+    cout << "  Unreliable: (i/o)  " << Form("%3d %3d", (int)unrin, (int)unrout) << endl; // Unzuverlaessig
+    cout << "  Isolated:   (i/o)  " << Form("%3d %3d", (int)isoin, (int)isoout) << endl; // Isolated (unbrauchbar)
+    cout << "  Max.Cluster:       " << Form("%3d", (int)clumax)                 << endl; // Max Cluster
+    cout << "  Arr Time inner:     " << meaninner << " +- " << rmsinner  << endl;
+    cout << "  Arr Time outer:     " << meanouter << " +- " << rmsouter  << endl;
+    cout << "  Mean Conv inner:     " << meanconvinner << endl;
+    cout << "  Mean Conv outer:     " << meanconvouter << endl;
+    cout << "  Rel.charge rms in:   " << relrmsinner << endl;
+    cout << "  Rel.charge rms out:  " << relrmsouter << endl;
+    cout << "  Med. num phe inner: " << numpheinner << endl;
+    cout << "  Med. num phe outer: " << numpheouter << endl;
+    cout << "  Ratio Calib Evts:   " << ratiocal << endl;
+
+    // FIXME: Fill calibration charge
 
     //inserting or updating the information in the database
     TString vars =
-        Form(" fUnsuitableInner=%d, "
+        Form(" fSequenceFirst=%d, "
+             " fTelescopeNumber=%d, "
+             " fUnsuitableInner=%d, "
              " fUnsuitableOuter=%d, "
              " fUnreliableInner=%d, "
@@ -254 +289 @@
              " fConvFactorInner=%s, "
              " fConvFactorOuter=%s, "
-             " fRatioCalEvents=%s ",
+             " fRatioCalEvents=%s, "
+             " fRelChargeRmsInner=%s, "
+             " fRelChargeRmsOuter=%s, "
+             " fMedNumPheInner=%s, "
+             " fMedNumPheOuter=%s ",
+             seq.GetSequence(), seq.GetTelescope(),
              (int)unsin, (int)unsout, (int)unrin,
              (int)unrout, (int)isoin, (int)isoout, (int)clumax,
@@ -260 +300 @@
              meanouter.Data(), rmsouter.Data(),
              meanconvinner.Data(), meanconvouter.Data(),
-             ratiocal.Data());
-
-    return serv.InsertUpdate("Calibration", "fSequenceFirst",
-                             Form("%d", seq), vars.Data()) ? 1 : 2;
+             ratiocal.Data(),
+             relrmsinner.Data(), relrmsouter.Data(),
+             numpheinner.Data(), numpheouter.Data()
+            );
+
+    TString where = Form("fSequenceFirst=%d AND fTelescopeNumber=%d",
+                         seq.GetSequence(), seq.GetTelescope());
+
+    return serv.InsertUpdate("Calibration", vars, where) ? 1 : 2;
 }
 
 int fillcalib(TString fname, Bool_t dummy=kTRUE)
 {
-    TEnv env("sql.rc");
-
-    MSQLMagic serv(env);
+    MSQLMagic serv("sql.rc");
     if (!serv.IsConnected())
     {