Changeset 3697 for trunk/MagicSoft/Mars/mcalib/MHCalibrationChargeCam.cc

Timestamp:

04/09/04 17:05:05 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MHCalibrationChargeCam.cc (modified) (18 diffs)

trunk/MagicSoft/Mars/mcalib/MHCalibrationChargeCam.cc

@@ -144 +144 @@
 const Float_t MHCalibrationChargeCam::fgNumHiGainSaturationLimit = 0.01;
 const Float_t MHCalibrationChargeCam::fgNumLoGainSaturationLimit = 0.005;
+const Float_t MHCalibrationChargeCam::fgTimeLowerLimit           = 1.;
+const Float_t MHCalibrationChargeCam::fgTimeUpperLimit           = 2.;
 // --------------------------------------------------------------------------
 //
@@ -154 +156 @@
 // - fNumHiGainSaturationLimit to fgNumHiGainSaturationLimit
 // - fNumLoGainSaturationLimit to fgNumLoGainSaturationLimit
+// - fTimeLowerLimit           to fgTimeLowerLimit 
+// - fTimeUpperLimit           to fgTimeUpperLimit 
 //
 MHCalibrationChargeCam::MHCalibrationChargeCam(const char *name, const char *title)
-    : fRawEvt(NULL)
-{
-    fName  = name  ? name  : "MHCalibrationChargeCam";
-    fTitle = title ? title : "Class to fill the calibration histograms ";
-
-    SetNumHiGainSaturationLimit();
-    SetNumLoGainSaturationLimit();
+    : fRawEvt(NULL), fSig(NULL)
+{
+  fName  = name  ? name  : "MHCalibrationChargeCam";
+  fTitle = title ? title : "Class to fill the calibration histograms ";
+  
+  SetNumHiGainSaturationLimit();
+  SetNumLoGainSaturationLimit();
+  SetTimeLowerLimit();
+  SetTimeUpperLimit();
 }
 
@@ -169 +175 @@
 // Gets the pointers to:
 // - MRawEvtData
+// - MExtractedSignalCam
 //
 Bool_t MHCalibrationChargeCam::SetupHists(const MParList *pList)
@@ -180 +187 @@
   }
 
+  MExtractedSignalCam *fSig = (MExtractedSignalCam*)pList->FindObject("MExtractedSignalCam");
+  if (!fSig)
+  {
+      *fLog << err << "MExtractedSignalCam not found... abort." << endl;
+      return kFALSE;
+  }
+
   return kTRUE;
 }
@@ -187 +201 @@
 // Gets or creates the pointers to:
 // - MCalibrationChargeCam
-//
-// Searches pointer to:
-// - MExtractedSignalCam
+// - MBadPixelsCam
+//
+// Initializes the number of used FADC slices from MExtractedSignalCam 
+// into MCalibrationChargeCam and test for changes in that variable
 //
 // Initializes, if empty to MGeomCam::GetNumPixels():
@@ -229 +244 @@
     }
 
-  MExtractedSignalCam *signal = (MExtractedSignalCam*)pList->FindObject("MExtractedSignalCam");
-  if (!signal)
-  {
-      *fLog << err << "MExtractedSignalCam not found... abort." << endl;
+  const Float_t numhigain = fSig->GetNumUsedHiGainFADCSlices();
+  const Float_t numlogain = fSig->GetNumUsedLoGainFADCSlices();  
+
+  if (fCam->GetNumHiGainFADCSlices() == 0.)
+    fCam->SetNumHiGainFADCSlices ( numhigain );
+  else if (fCam->GetNumHiGainFADCSlices() != numhigain)
+    {
+      *fLog << err << GetDescriptor() 
+            << ": Number of High Gain FADC extraction slices has changes, abort..." << endl;
       return kFALSE;
-  }
+    }
+  
+  if (fCam->GetNumLoGainFADCSlices() == 0.)
+    fCam->SetNumLoGainFADCSlices ( numlogain );
+  else if (fCam->GetNumLoGainFADCSlices() != numlogain)
+    {
+      *fLog << err << GetDescriptor() 
+            << ": Number of Low Gain FADC extraction slices has changes, abort..." << endl;
+      return kFALSE;
+    }
 
   const Int_t npixels  = fGeom->GetNumPixels();
@@ -532 +561 @@
 {
 
+  const Byte_t firsthi = fSig->GetFirstUsedSliceHiGain();
+  const Byte_t lasthi  = fSig->GetLastUsedSliceHiGain();
+  const Byte_t firstlo = fSig->GetFirstUsedSliceLoGain();
+  const Byte_t lastlo  = fSig->GetLastUsedSliceLoGain();
+
+
   for (Int_t i=0; i<fHiGainArray->GetSize(); i++)
     {
@@ -537 +572 @@
       MHCalibrationChargePix &histhi = (MHCalibrationChargePix&)(*this)[i];
       MCalibrationChargePix  &pix    = (MCalibrationChargePix&)(*fCam)[i];
+      MBadPixelsPix          &bad    = (*fBadPixels)[i];
       
       if (histhi.IsExcluded())
@@ -548 +584 @@
         }
 
-      pix.SetAbsTimeMean ( histhi.GetAbsTimeMean());
-      pix.SetAbsTimeRms  ( histhi.GetAbsTimeRms() );
+      FinalizeAbsTimes(histhi, pix, bad, firsthi, lasthi);
     }
 
@@ -572 +607 @@
       
       if (pix.IsHiGainSaturation())
-        {
-          pix.SetAbsTimeMean     ( histlo.GetAbsTimeMean());
-          pix.SetAbsTimeRms      ( histlo.GetAbsTimeRms() );
-        }           
+        FinalizeAbsTimes(histlo, pix, bad, firstlo, lastlo);
     }
 
@@ -583 +615 @@
       MHCalibrationChargePix &histhi = (MHCalibrationChargePix&)GetAverageHiGainArea(j);      
       MCalibrationChargePix  &pix    = (MCalibrationChargePix&)fCam->GetAverageArea(j);
+      MBadPixelsPix          &bad    = fCam->GetAverageBadArea(j);
       
       if (histhi.GetSaturated() > fNumHiGainSaturationLimit*histhi.GetHGausHist()->GetEntries())
@@ -591 +624 @@
         }
 
-      pix.SetAbsTimeMean ( histhi.GetAbsTimeMean());
-      pix.SetAbsTimeRms  ( histhi.GetAbsTimeRms() );
-    }
+      FinalizeAbsTimes(histhi, pix, bad, firsthi, lasthi);
+   }
   
   for (Int_t j=0; j<fAverageLoGainAreas->GetSize(); j++)
@@ -600 +632 @@
       MHCalibrationChargePix &histlo = (MHCalibrationChargePix&)GetAverageLoGainArea(j);      
       MCalibrationChargePix  &pix    = (MCalibrationChargePix&)fCam->GetAverageArea(j);
-      
+      MBadPixelsPix          &bad    = fCam->GetAverageBadArea(j);      
+
       if (histlo.GetSaturated() > fNumLoGainSaturationLimit*histlo.GetHGausHist()->GetEntries())
         {
@@ -609 +642 @@
 
       if (pix.IsHiGainSaturation())
-        {
-          pix.SetAbsTimeMean  ( histlo.GetAbsTimeMean());
-          pix.SetAbsTimeRms   ( histlo.GetAbsTimeRms() );
-        }           
+        FinalizeAbsTimes(histlo, pix, bad, firstlo, lastlo);
     }
 
@@ -620 +650 @@
       MHCalibrationChargePix &histhi = (MHCalibrationChargePix&)GetAverageHiGainSector(j);      
       MCalibrationChargePix  &pix    = (MCalibrationChargePix&)fCam->GetAverageSector(j);
-      
+      MBadPixelsPix          &bad    = fCam->GetAverageBadSector(j);      
+
       if (histhi.GetSaturated() > fNumHiGainSaturationLimit*histhi.GetHGausHist()->GetEntries())
         {
@@ -628 +659 @@
         }
 
-      pix.SetAbsTimeMean ( histhi.GetAbsTimeMean());
-      pix.SetAbsTimeRms  ( histhi.GetAbsTimeRms() );
+      FinalizeAbsTimes(histhi, pix, bad, firsthi, lasthi);
     }
   
@@ -648 +678 @@
 
       if (pix.IsHiGainSaturation())
-        {
-          pix.SetAbsTimeMean  ( histlo.GetAbsTimeMean());
-          pix.SetAbsTimeRms   ( histlo.GetAbsTimeRms() );
-        }           
+        FinalizeAbsTimes(histlo, pix, bad, firstlo, lastlo);
     }
   
@@ -670 +697 @@
 }
 
+// --------------------------------------------------------------------------------
+//
+// Fill the absolute time results into MCalibrationChargePix
+//
+// Check absolute time validity:
+// - Mean arrival time is at least fTimeLowerLimit slices from the lower edge 
+// - Mean arrival time is at least fUpperLimit     slices from the upper edge
+//
+void MHCalibrationChargeCam::FinalizeAbsTimes(MHCalibrationChargePix &hist, MCalibrationChargePix &pix, MBadPixelsPix &bad, 
+                                              Byte_t first, Byte_t last)
+{
+  
+  const Float_t mean = hist.GetAbsTimeMean();
+  const Float_t rms  = hist.GetAbsTimeRms();
+
+  pix.SetAbsTimeMean ( mean );
+  pix.SetAbsTimeRms  ( rms  );
+  
+  const Float_t lowerlimit = (Float_t)first + fTimeLowerLimit;
+  const Float_t upperlimit = (Float_t)last  + fTimeUpperLimit;  
+
+  if ( mean < lowerlimit)
+    {
+      *fLog << warn << GetDescriptor() 
+            << Form("%s%3.1f%s%2i%s%4i",": Mean ArrivalTime: ",mean," smaller than ",fTimeLowerLimit,
+                    " FADC slices from lower edge in pixel ",hist.GetPixId()) << endl;
+      bad.SetUncalibrated( MBadPixelsPix::kMeanTimeInFirstBin );
+    }
+  
+  if ( mean  > upperlimit )
+    {
+      *fLog << warn << GetDescriptor() 
+            << Form("%s%3.1f%s%2i%s%4i",": Mean ArrivalTime: ",mean," greater than ",fTimeUpperLimit,
+                    " FADC slices from upper edge in pixel ",hist.GetPixId()) << endl;
+      bad.SetUncalibrated( MBadPixelsPix::kMeanTimeInLast2Bins );
+    }
+}
+
 // --------------------------------------------------------------------------
 //