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

Timestamp:

08/02/06 10:33:31 (18 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

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

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

@@ -118 +118 @@
     }
 
+    MHCamera *difflo = (MHCamera*)arr.FindObjectInCanvas("DiffLo;avg", "MHCamera", "DiffLo");
+    if (!difflo)
+    {
+        cout << "WARNING - Reading of DiffLo;avg failed." << endl;
+        return 2;
+    }
+    MHCamera *diffhi = (MHCamera*)arr.FindObjectInCanvas("DiffHi;avg", "MHCamera", "DiffHi");
+    if (!diffhi)
+    {
+        cout << "WARNING - Reading of DiffHi;avg failed." << endl;
+        return 2;
+    }
+
+    MHCamera *pulhi = (MHCamera*)arr.FindObjectInCanvas("PulseHi;avg", "TH1D", "PulseHi");
+    if (!pulhi)
+    {
+        cout << "WARNING - Reading of PulseHi failed." << endl;
+        return 2;
+    }
+    MHCamera *pullo = (MHCamera*)arr.FindObjectInCanvas("PulseLo;avg", "TH1D", "PulseLo");
+    if (!pullo)
+    {
+        cout << "WARNING - Reading of PulseLo failed." << endl;
+        return 2;
+    }
+
     TArrayI inner(1);
     inner[0] = 0;
@@ -173 +199 @@
 
     meanpul = TMath::Nint(meanpul*100)/100.;
-    rmspul  = TMath::Nint(rmspul*100)/100.;
+    rmspul  = TMath::Nint(rmspul *100)/100.;
+
+
+    Double_t meanhi  = TMath::Nint(pulhi->GetMean()*100.)/100.;
+    Double_t rmshi   = TMath::Nint(pulhi->GetRMS() *100.)/100.;
+
+    Double_t meanlo  = TMath::Nint(pullo->GetMean()*100.)/100.;
+    Double_t rmslo   = TMath::Nint(pullo->GetRMS() *100.)/100.;
+
+    pullo->Add(pullo, pulhi, 1, -1);
+    pullo->ResetBit(MHCamera::kProfile);
+
+    Double_t meanoff = TMath::Nint(pullo->GetMean()*100.)/100.;
+    Double_t rmsoff  = TMath::Nint(pullo->GetRMS() *100.)/100.;
 
     //get sequence number from the filename
@@ -194 +233 @@
     TString rmsextpulpos =Form("%6.2f", rmsextpul);
 
+    TString meanpulhi    =Form("%6.2f", meanhi);
+    TString rmspulhi     =Form("%6.2f", rmshi);
+
+    TString meanpullo    =Form("%6.2f", meanlo);
+    TString rmspullo     =Form("%6.2f", rmslo);
+
+    TString meanpuloff   =Form("%6.2f", meanoff);
+    TString rmspuloff    =Form("%6.2f", rmsoff);
+
     if (meanextpul<0 && rmsextpul<0)
     {
@@ -200 +248 @@
     }
 
+    // *****************************************************
+
+    // *****************************************************
+
+
     cout << "Sequence #" << seq << endl;
-    cout << "  Mean Ped RMS inner [phe] " << meanrmsinner << endl;
-    cout << "  Mean Ped RMS outer [phe] " << meanrmsouter << endl;
-    cout << "  Mean Signal  inner [phe] " << meansiginner << endl;
-    cout << "  Mean Signal  outer [phe] " << meansigouter << 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;
+    cout << "  Mean Ped RMS inner [phe] " << meanrmsinner  << endl;
+    cout << "  Mean Ped RMS outer [phe] " << meanrmsouter  << endl;
+    cout << "  Mean Signal  inner [phe] " << meansiginner  << endl;
+    cout << "  Mean Signal  outer [phe] " << meansigouter  << endl;
+    cout << "  Mean calibrated PulsePos " << meanpulpos    << " +- " << rmspulpos    << endl;
+    cout << "  Mean extracted  PulsePos " << meanextpulpos << " +- " << rmsextpulpos << endl;
+    cout << "  Mean ext.HiGain PulsePos " << meanpulhi     << " +- " << rmspulhi     << endl;
+    cout << "  Mean ext.LoGain PulsePos " << meanpullo     << " +- " << rmspullo     << endl;
+    cout << "  Lo-Hi gain offset:       " << meannpuloff   << " +- " << rmspuloff    << endl;
+    cout << endl;
 
     //build query and insert information into the database
@@ -219 +274 @@
                          " fPulsePosMean=%s,      fPulsePosRms=%s,     "
                          " fPulsePosCheckMean=%s, fPulsePosCheckRms=%s "
+                         " fPulsePosHiMean=%s,    fPulsePosHiRms=%s,   "
+                         " fPulsePosLoMean=%s,    fPulsePosLoRms=%s,   "
+                         " fPulsePosOffMean=%s,   fPulsePosOffRms=%s,  "
                          " WHERE fSequenceFirst='%d' ",
                          meanrmsinner.Data(),  meanrmsouter.Data(),
@@ -224 +282 @@
                          meanpulpos.Data(),    rmspulpos.Data(),
                          meanextpulpos.Data(), rmsextpulpos.Data(),
+                         meanpulhi.Data(),     rmspulhi.Data(),
+                         meanpullo.Data(),     rmspullo.Data(),
+                         meanpuloff.Data(),    rmspuloff.Data(),
                          seq);