Changeset 4903 for trunk/MagicSoft/Mars/mcalib/MHCalibrationRelTimeCam.cc

Timestamp:

09/09/04 17:37:26 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MHCalibrationRelTimeCam.cc (modified) (19 diffs)

trunk/MagicSoft/Mars/mcalib/MHCalibrationRelTimeCam.cc

@@ -27 +27 @@
 //                                                                        
 // Fills the extracted relative arrival times of MArrivalTimeCam into 
-// the MHGausEvents-classes MHCalibrationRelTimePix for every:
+// the MHCalibrationPix-classes MHCalibrationRelTimePix for every:
 //
 // - Pixel, stored in the TObjArray's MHCalibrationCam::fHiGainArray  
@@ -52 +52 @@
 // the fit is declared valid.
 // Otherwise, the fit is repeated within ranges of the previous mean 
-// +- MHGausEvents::fPickupLimit (default: 5) sigma (see MHGausEvents::RepeatFit())
+// +- MHCalibrationPix::fPickupLimit (default: 5) sigma (see MHCalibrationPix::RepeatFit())
 // In case this does not make the fit valid, the histogram means and RMS's are 
-// taken directly (see MHGausEvents::BypassFit()) and the following flags are set:
+// taken directly (see MHCalibrationPix::BypassFit()) and the following flags are set:
 // - MBadPixelsPix::SetUncalibrated( MBadPixelsPix::kRelTimeNotFitted ) and
 // - MBadPixelsPix::SetUnsuitable(   MBadPixelsPix::kUnreliableRun    ) 
 // 
-// Outliers of more than MHGausEvents::fPickupLimit (default: 5) sigmas 
-// from the mean are counted as Pickup events (stored in MHGausEvents::fPickup) 
+// Outliers of more than MHCalibrationPix::fPickupLimit (default: 5) sigmas 
+// from the mean are counted as Pickup events (stored in MHCalibrationPix::fPickup) 
 //
 // The class also fills arrays with the signal vs. event number, creates a fourier 
@@ -88 +88 @@
 #include "MHCalibrationRelTimePix.h"
 
+#include "MHCalibrationPix.h"
+
 #include "MLog.h"
 #include "MLogManip.h"
@@ -93 +95 @@
 #include "MParList.h"
 
+#include "MCalibrationIntensityRelTimeCam.h"
+
 #include "MCalibrationRelTimeCam.h"
+#include "MCalibrationRelTimePix.h"
 #include "MCalibrationPix.h"
 
@@ -104 +109 @@
 #include "MBadPixelsCam.h"
 #include "MBadPixelsPix.h"
-
-#include "MHGausEvents.h"
 
 ClassImp(MHCalibrationRelTimeCam);
@@ -163 +166 @@
 {
 
-  
-  fCam = (MCalibrationCam*)pList->FindObject("MCalibrationRelTimeCam");
-  if (!fCam)
-    {
-      fCam = (MCalibrationCam*)pList->FindCreateObj(AddSerialNumber("MCalibrationRelTimeCam"));
+  fIntensCam = (MCalibrationIntensityCam*)pList->FindObject(AddSerialNumber("MCalibrationIntensityRelTimeCam"));
+  if (fIntensCam)
+    *fLog << inf << "Found MCalibrationIntensityRelTimeCam ... " << endl;
+  else
+    {
+      fCam = (MCalibrationCam*)pList->FindObject(AddSerialNumber("MCalibrationRelTimeCam"));
       if (!fCam)
         {
-          gLog << err << "Cannot find nor create MCalibrationRelTimeCam ... abort." << endl;
-          return kFALSE;
-        }
-      else 
-        fCam->Init(*fGeom);
-    }
-  
+          fCam = (MCalibrationCam*)pList->FindCreateObj(AddSerialNumber("MCalibrationRelTimeCam"));
+          if (!fCam)
+            {
+              *fLog << err << "Cannot find nor create MCalibrationRelTimeCam ... abort." << endl;
+              return kFALSE;
+            }
+          fCam->Init(*fGeom);
+        }
+    }
 
   MArrivalTimeCam *signal = (MArrivalTimeCam*)pList->FindObject("MArrivalTimeCam");
@@ -219 +225 @@
       {
         (*fAverageHiGainAreas)[j] = 
-          new MHCalibrationRelTimePix("MHCalibrationRelTimeAverageAreaHiGain",
+          new MHCalibrationRelTimePix("RelTimeAverageHiGainArea",
                                       "Average Rel. Arr. Times Hi-Gain Area Idx ");
 
-        GetAverageHiGainArea(j).GetHGausHist()->SetTitle("Rel. Arr. Times HiGain Area Idx ");
-        GetAverageHiGainArea(j).SetNbins(fAverageNbins);
-
-        InitHists(GetAverageHiGainArea(j),fCam->GetAverageBadArea(j),j);
+        MHCalibrationRelTimePix &hist = (MHCalibrationRelTimePix&)GetAverageHiGainArea(j);
+
+        hist.SetNbins(fAverageNbins);
+        hist.GetHGausHist()->SetTitle("Rel. Arr. Times average HiGain Area Idx ");
+
+        if (fGeom->InheritsFrom("MGeomCamMagic"))
+          {
+            hist.GetHGausHist()->SetTitle(Form("%s%s%s","Signal averaged on event-by-event basis ",
+                                               j==0 ? "Inner Pixels " : "Outer Pixels ","High Gain Runs: "));
+            hist.InitBins();
+            hist.SetEventFrequency(fPulserFrequency);
+          }
+        else
+          {
+            hist.GetHGausHist()->SetTitle("Signal averaged on event-by-event basis High Gain Area Idx ");
+            InitHists(hist,fIntensCam ? fIntensCam->GetAverageBadArea(j) : fCam->GetAverageBadArea(j),j);
+          }
       }
   }
@@ -236 +255 @@
       {
         (*fAverageLoGainAreas)[j] = 
-          new MHCalibrationRelTimePix("MHCalibrationRelTimeAverageAreaLoGain",
+          new MHCalibrationRelTimePix("RelTimeAverageAreaLoGain",
                                       "Average Rel. Arr. Times Lo-Gain Area Idx ");
 
-        GetAverageLoGainArea(j).GetHGausHist()->SetTitle("Rel. Arr. Times LoGain Area Idx ");
-        GetAverageLoGainArea(j).SetNbins(fAverageNbins);
-
-        InitHists(GetAverageLoGainArea(j),fCam->GetAverageBadArea(j),j);
+        MHCalibrationRelTimePix &hist = (MHCalibrationRelTimePix&)GetAverageLoGainArea(j);
+        hist.SetNbins(fAverageNbins);
+
+        if (fGeom->InheritsFrom("MGeomCamMagic"))
+          {
+            hist.GetHGausHist()->SetTitle(Form("%s%s%s","Rel. Arr. Times averaged on event-by-event basis ",
+                                               j==0 ? "Inner Pixels " : "Outer Pixels ","Low Gain Runs: "));
+            hist.InitBins();
+            hist.SetEventFrequency(fPulserFrequency);
+          }
+        else
+          {
+            hist.GetHGausHist()->SetTitle("Rel. Arr. Times averaged on event-by-event basis Low Gain Area Idx ");
+            InitHists(hist,fIntensCam ? fIntensCam->GetAverageBadArea(j) : fCam->GetAverageBadArea(j),j);
+          }
       }
   }
@@ -253 +283 @@
         {
           (*fAverageHiGainSectors)[j] = 
-            new MHCalibrationRelTimePix("MHCalibrationRelTimeAverageSectorHiGain",
-                                        "Average Rel. Arr. Times Hi-Gain Sector ");
+            new MHCalibrationRelTimePix("RelTimeAverageSectorHiGain",
+                                        "Average HiGain Rel. Arr. Times Sector ");
           
-          GetAverageHiGainSector(j).GetHGausHist()->SetTitle("Rel. Arr. Times HiGain Sector ");
-          GetAverageHiGainSector(j).SetNbins(fAverageNbins);
+          MHCalibrationRelTimePix &hist = (MHCalibrationRelTimePix&)GetAverageHiGainSector(j);
+
+          hist.GetHGausHist()->SetTitle("Rel. Arr. Times averaged on event-by-event basis HiGain Sector ");
+          hist.SetNbins(fAverageNbins);
           
-          InitHists(GetAverageHiGainSector(j),fCam->GetAverageBadSector(j),j);
+          InitHists(hist,fIntensCam ? fIntensCam->GetAverageBadSector(j) : fCam->GetAverageBadSector(j),j);
       }
   }
@@ -270 +302 @@
         {
           (*fAverageLoGainSectors)[j] = 
-            new MHCalibrationRelTimePix("MHCalibrationRelTimeAverageSectorLoGain",
-                                        "Average Rel. Arr. Times Lo-Gain Sector ");
+            new MHCalibrationRelTimePix("RelTimeAverageSectorLoGain",
+                                        "Average LoGain Rel. Arr. Times Sector ");
           
-          GetAverageLoGainSector(j).GetHGausHist()->SetTitle("Rel. Arr. Times LoGain Sector ");
-          GetAverageLoGainSector(j).SetNbins(fAverageNbins);
+          MHCalibrationRelTimePix &hist = (MHCalibrationRelTimePix&)GetAverageLoGainSector(j);
+
+          hist.GetHGausHist()->SetTitle("Rel. Arr. Times averaged on event-by-event basis LoGain Sector ");
+          hist.SetNbins(fAverageNbins);
           
-          InitHists(GetAverageLoGainSector(j),fCam->GetAverageBadSector(j),j);
-          
+          InitHists(hist,fIntensCam ? fIntensCam->GetAverageBadSector(j) : fCam->GetAverageBadSector(j),j);
         }
     }
@@ -338 +371 @@
     {
 
-      MHGausEvents &histhi = (*this)[i];
-      MHGausEvents &histlo = (*this)(i);
+      MHCalibrationPix &histhi = (*this)[i];
+      MHCalibrationPix &histlo = (*this)(i);
 
       if (histhi.IsExcluded())
@@ -372 +405 @@
   for (Int_t j=0; j<nareas; j++)
     {
-      MHGausEvents &histhi = GetAverageHiGainArea(j);
+      MHCalibrationPix &histhi = GetAverageHiGainArea(j);
       histhi.FillHistAndArray(fNumareahi[j] == 0 ? 0. : fSumareahi[j]/fNumareahi[j]);
 
       if (fLoGain)
         {
-          MHGausEvents &histlo = GetAverageLoGainArea(j);
+          MHCalibrationPix &histlo = GetAverageLoGainArea(j);
           histlo.FillHistAndArray(fNumarealo[j] == 0 ? 0. : fSumarealo[j]/fNumarealo[j]);
         }
@@ -384 +417 @@
   for (Int_t j=0; j<nsectors; j++)
     {
-      MHGausEvents &histhi = GetAverageHiGainSector(j);
+      MHCalibrationPix &histhi = GetAverageHiGainSector(j);
       histhi.FillHistAndArray(fNumsectorhi[j] == 0 ? 0. : fSumsectorhi[j]/fNumsectorhi[j]);
 
       if (fLoGain)
         {
-          MHGausEvents &histlo = GetAverageLoGainSector(j);
+          MHCalibrationPix &histlo = GetAverageLoGainSector(j);
           histlo.FillHistAndArray(fNumsectorlo[j] == 0 ? 0. : fSumsectorlo[j]/fNumsectorlo[j]);
         }
@@ -415 +448 @@
         continue;
       
+      MCalibrationRelTimePix  &pix    = fIntensCam 
+        ? (MCalibrationRelTimePix&)(*fIntensCam)[i] 
+        : (MCalibrationRelTimePix&)(*fCam)[i];
+
       if (histhi.GetSaturated() > fNumHiGainSaturationLimit*histhi.GetHGausHist()->GetEntries())
         {
-          (*fCam)[i].SetHiGainSaturation();
+          pix.SetHiGainSaturation();
           histhi.SetExcluded();
         }
@@ -423 +460 @@
         if (fLoGain)
           (*this)(i).SetExcluded();
-      
-    }
-
+
+      Stat_t overflow = histhi.GetHGausHist()->GetBinContent(histhi.GetHGausHist()->GetNbinsX()+1);
+      if (overflow > 0.1)
+        {
+          *fLog << warn << GetDescriptor()
+                << ": HiGain Histogram Overflow occurred " << overflow 
+                << " times in pixel: " << i << " (without saturation!) " << endl;
+          //          bad.SetUncalibrated( MBadPixelsPix::kHiGainOverFlow ); 
+        }
+
+      overflow = histhi.GetHGausHist()->GetBinContent(0);
+      if (overflow > 0.1)
+        {
+          *fLog << warn << GetDescriptor()
+                << ": HiGain Histogram Underflow occurred " << overflow
+                << " times in pixel: " << i << " (without saturation!) " << endl;
+          //          bad.SetUncalibrated( MBadPixelsPix::kHiGainOverFlow ); 
+        }
+    }
 
   for (Int_t j=0; j<fAverageHiGainAreas->GetSize(); j++)
@@ -434 +487 @@
       if (histhi.GetSaturated() > fNumHiGainSaturationLimit*histhi.GetHGausHist()->GetEntries())
         {
-          fCam->GetAverageArea(j).SetHiGainSaturation();
+          MCalibrationRelTimePix  &pix    = fIntensCam 
+            ? (MCalibrationRelTimePix&)fIntensCam->GetAverageArea(j)
+            : (MCalibrationRelTimePix&)fCam->GetAverageArea(j);
+          pix.SetHiGainSaturation();
           histhi.SetExcluded();
         }
@@ -447 +503 @@
       
       MHCalibrationRelTimePix &histhi = (MHCalibrationRelTimePix&)GetAverageHiGainSector(j);      
+      MCalibrationRelTimePix  &pix    = fIntensCam
+        ? (MCalibrationRelTimePix&)fIntensCam->GetAverageSector(j)  
+        : (MCalibrationRelTimePix&)fCam->GetAverageSector(j);
 
       if (histhi.GetSaturated() > fNumHiGainSaturationLimit*histhi.GetHGausHist()->GetEntries())
         {
-          fCam->GetAverageSector(j).SetHiGainSaturation();
+          pix.SetHiGainSaturation();
           histhi.SetExcluded();
         }
@@ -458 +517 @@
     }
 
-  FitHiGainArrays((*fCam),*fBadPixels,
+  FitHiGainArrays(fIntensCam ? (MCalibrationCam&)(*fIntensCam->GetCam()) : (MCalibrationCam&)(*fCam),
+                  *fBadPixels,
                   MBadPixelsPix::kRelTimeNotFitted,
                   MBadPixelsPix::kRelTimeOscillating);
   
   if (fLoGain)        
-    FitLoGainArrays((*fCam),*fBadPixels,
+    FitLoGainArrays(fIntensCam ? (MCalibrationCam&)(*fIntensCam->GetCam()) : (MCalibrationCam&)(*fCam),
+                    *fBadPixels,
                     MBadPixelsPix::kRelTimeNotFitted,
                     MBadPixelsPix::kRelTimeOscillating);
@@ -557 +618 @@
 // --------------------------------------------------------------------------
 //
-// Calls MHGausEvents::DrawClone() for pixel idx
+// Calls MHCalibrationPix::DrawClone() for pixel idx
 //
 void MHCalibrationRelTimeCam::DrawPixelContent(Int_t idx) const