Changeset 6685 for trunk/MagicSoft/Mars/mhcalib/MHCalibrationPulseTimeCam.cc

Timestamp:

02/28/05 17:54:50 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mhcalib/MHCalibrationPulseTimeCam.cc (modified) (3 diffs)

trunk/MagicSoft/Mars/mhcalib/MHCalibrationPulseTimeCam.cc

@@ -149 +149 @@
 const Axis_t  MHCalibrationPulseTimeCam::fgHiGainLast  =  19.5;
 const Float_t MHCalibrationPulseTimeCam::fgProbLimit   = 0.0001;
-const Float_t MHCalibrationPulseTimeCam::fgLowerSignalLimit = 100.;
+const Int_t   MHCalibrationPulseTimeCam::fgLowerSignalLimit = 100;
 const TString MHCalibrationPulseTimeCam::gsHistName    = "PulseTime";
 const TString MHCalibrationPulseTimeCam::gsHistTitle   = "Extracted Times";
@@ -348 +348 @@
         continue;
 
-      const Float_t sighi = pix.GetExtractedSignalHiGain();
-
-      if (sighi < fLowerSignalLimit)
+      Byte_t *start = pixel.GetHiGainSamples();
+      Byte_t *end   = start + pixel.GetNumHiGainSamples();
+      Byte_t *p     = start;
+      Byte_t max    = 0;
+      Int_t  maxpos = 0;
+
+      while (p < end)
+        {
+          if (*p++ > max)
+            maxpos = p-start-1;
+        }
+
+      start   = pixel.GetLoGainSamples();
+      end     = start + pixel.GetNumLoGainSamples();
+      p       = start;
+
+      while (p < end)
+        {
+          if (*p++ > max)
+            maxpos = p-start+pixel.GetNumHiGainSamples() - 1;
+        }
+
+      if ((Int_t)max < fLowerSignalLimit)
         continue;
 
-      const Byte_t slicehi = pixel.GetMaxHiGainSample();
-      const Byte_t slicelo = pixel.GetMaxLoGainSample();
-
-      Float_t time = 0.;
-
-      if (slicehi > slicelo)
-        time = (Float_t)pixel.GetIdxMaxHiGainSample();
-      else
-        time = (Float_t)pixel.GetIdxMaxLoGainSample();
-
+      const Float_t time = (Float_t)maxpos;
       histhi.FillHist(time);
         
@@ -423 +434 @@
   *fLog << endl;
 
-  MCalibrationCam *pulsecam  = fIntensCam ? fIntensCam->GetCam() : fCam;
-  MBadPixelsCam   *badcam    = fIntensBad ? fIntensBad->GetCam() : fBadPixels;
-
+  MCalibrationCam *calcam  = fIntensCam ? fIntensCam->GetCam() : fCam;
   //
   // Perform the fitting for the High Gain (done in MHCalibrationCam)
   // 
-  FitHiGainArrays(*pulsecam, *badcam, 
-                  MBadPixelsPix::kHiGainNotFitted,
-                  MBadPixelsPix::kHiGainOscillating);
+  for (Int_t i=0; i<fHiGainArray->GetSize(); i++)
+    {
+      
+      MHCalibrationPix &hist = (*this)[i];
+      
+      if (hist.IsExcluded())
+        continue;
+      
+      MCalibrationPix &pix    = (*calcam)[i];
+      CalcHists(hist,pix);
+    }
+
+  if (!IsAverageing())
+    return kTRUE;
+  
+  for (Int_t j=0; j<fAverageHiGainAreas->GetSize(); j++)
+    {
+      
+      MHCalibrationPix &hist = GetAverageHiGainArea(j);      
+      MCalibrationPix  &pix  = calcam->GetAverageArea(j);
+      CalcHists(hist,pix);
+  }
+  
+  for (Int_t j=0; j<fAverageHiGainSectors->GetSize(); j++)
+    {
+      MHCalibrationPix &hist = GetAverageHiGainSector(j);      
+      MCalibrationPix  &pix  = calcam->GetAverageSector(j);
+      CalcHists(hist,pix);
+    }
   
   return kTRUE;
 }
+
+void MHCalibrationPulseTimeCam::CalcHists(MHCalibrationPix &hist, MCalibrationPix &pix)
+{
+
+
+  if (hist.IsEmpty() || hist.IsOnlyOverflow() || hist.IsOnlyUnderflow())
+    {
+      *fLog << warn << GetDescriptor()
+            << ": Only overflow or underflow in hi-gain pixel: " << pix.GetPixId() << endl;
+      return;
+    }  
+
+  hist.BypassFit();
+
+  pix.SetHiGainMean       ( hist.GetMean()      );
+  pix.SetHiGainMeanVar    ( hist.GetMeanErr() * hist.GetMeanErr()   );
+  pix.SetHiGainRms        ( hist.GetHistRms()   );
+  pix.SetHiGainSigma      ( hist.GetSigma()     );
+  pix.SetHiGainSigmaVar   ( hist.GetSigmaErr()* hist.GetSigmaErr()  );
+
+  if (IsDebug())
+    {
+      *fLog << dbginf << GetDescriptor() << ": ID " << GetName() 
+            << " Mean: "         << hist.GetMean    ()
+            << " MeanErr: "      << hist.GetMeanErr ()
+            << " MeanSigma: "    << hist.GetSigma   ()
+            << " MeanSigmaErr: " << hist.GetSigmaErr()
+            << " Prob: "         << hist.GetProb    ()
+            << endl;
+    }
+  return;
+}
+
 
 // --------------------------------------------------------------------------