Changeset 7428 for trunk/MagicSoft/Mars/datacenter/macros/fillsignal.C

Timestamp:

11/25/05 11:29:14 (19 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

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

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

@@ -65 +65 @@
 #include "MHCamera.h"
 
+#include "MCalibrationPulseTimeCam.h"
+#include "MCalibrationPix.h"
+
 using namespace std;
 
@@ -76 +79 @@
     }
 
-    MStatusArray *arr = (MStatusArray*)file.Get("MStatusDisplay");
-    if (!arr)
+    Float_t meanextpul = -1;
+    Float_t rmsextpul  = -1;
+
+    MCalibrationPulseTimeCam *pt;
+    file.GetObject("MCalibrationPulseTimeCam", pt);
+    if (pt)
+    {
+        meanextpul = pt->GetAverageArea(0).GetHiGainMean();
+        rmsextpul  = pt->GetAverageArea(0).GetHiGainRms();
+
+        meanextpul = TMath::Nint(meanextpul*100)/100.;
+        rmsextpul  = TMath::Nint(rmsextpul*100)/100.;
+    }
+
+    MStatusArray arr;
+    if (arr.Read()<=0)
     {
         cout << "ERROR - Reading of MStatusDisplay failed." << endl;
@@ -83 +100 @@
     }
 
-    MHCamera *cam = (MHCamera*)arr->FindObjectInCanvas("PedRMS;avg", "MHCamera", "PedRMS");
+    MHCamera *cam = (MHCamera*)arr.FindObjectInCanvas("PedRMS;avg", "MHCamera", "PedRMS");
     if (!cam)
     {
@@ -90 +107 @@
     }
 
-    MHCamera *pul = (MHCamera*)arr->FindObjectInCanvas("PulsePos;avg", "MHCamera", "PulsePos");
+    MHCamera *pul = (MHCamera*)arr.FindObjectInCanvas("PulsePos;avg", "MHCamera", "PulsePos");
     if (!pul)
     {
@@ -119 +136 @@
     meanrmso = TMath::Nint(meanrmso*100)/100.;
 
-    cam = (MHCamera*)arr->FindObjectInCanvas("Interp'd;avg", "MHCamera", "Interp'd");
+    cam = (MHCamera*)arr.FindObjectInCanvas("Interp'd;avg", "MHCamera", "Interp'd");
     if (!cam)
     {
@@ -163 +180 @@
 
     Int_t seq = atoi(sequence.Data()+6);
-    TString meanrmsinner=Form("%6.2f", meanrmsi);
-    TString meanrmsouter=Form("%6.2f", meanrmso);
-    TString meansiginner=Form("%6.2f", meansigi);
-    TString meansigouter=Form("%6.2f", meansigo);
-    TString meanpulpos  =Form("%6.2f", meanpul);
-    TString rmspulpos   =Form("%6.2f", rmspul);
+    TString meanrmsinner =Form("%6.2f", meanrmsi);
+    TString meanrmsouter =Form("%6.2f", meanrmso);
+    TString meansiginner =Form("%6.2f", meansigi);
+    TString meansigouter =Form("%6.2f", meansigo);
+    TString meanpulpos   =Form("%6.2f", meanpul);
+    TString rmspulpos    =Form("%6.2f", rmspul);
+    TString meanextpulpos=Form("%6.2f", meanextpul);
+    TString rmsextpulpos =Form("%6.2f", rmsextpul);
+
+    if (meanextpul<0 && rmsextpul<0)
+    {
+        meanextpulpos = "NULL";
+        rmsextpulpos  = "NULL";
+    }
 
     cout << "Sequence #" << seq << endl;
@@ -175 +200 @@
     cout << "  Mean Signal  inner [phe] " << meansiginner << endl;
     cout << "  Mean Signal  outer [phe] " << meansigouter << endl;
-    cout << "  Mean PulsePos            " << meanpulpos   << endl;
-    cout << "  Rms  PulsePos            " << rmspulpos    << endl;
+    cout << "  Mean calibrated PulsePos " << meanpulpos   << endl;
+    cout << "  Rms  calibrated PulsePos " << rmspulpos    << endl;
+    cout << "  Mean extracted  PulsePos " << meanextpulpos   << endl;
+    cout << "  Rms  extracted  PulsePos " << rmsextpulpos    << endl;
 
     TString query = Form("UPDATE Calibration SET "
-                         " fMeanPedRmsInner=%s, fMeanPedRmsOuter=%s, "
-                         " fMeanSignalInner=%s, fMeanSignalOuter=%s, "
-                         " fPulsePosMean=%s,    fPulsePosRms=%s      "
+                         " fMeanPedRmsInner=%s,   fMeanPedRmsOuter=%s, "
+                         " fMeanSignalInner=%s,   fMeanSignalOuter=%s, "
+                         " fPulsePosMean=%s,      fPulsePosRms=%s,     "
+                         " fPulsePosCheckMean=%s, fPulsePosCheckRms=%s "
                          " WHERE fSequenceFirst='%d' ",
-                         meanrmsinner.Data(), meanrmsouter.Data(),
-                         meansiginner.Data(), meansigouter.Data(),
-                         meanpulpos.Data(),   rmspulpos.Data(),
+                         meanrmsinner.Data(),  meanrmsouter.Data(),
+                         meansiginner.Data(),  meansigouter.Data(),
+                         meanpulpos.Data(),    rmspulpos.Data(),
+                         meanextpulpos.Data(), rmsextpulpos.Data(),
                          seq);