Changeset 2801

Timestamp:

01/14/04 14:03:26 (21 years ago)

Author:

moralejo

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• manalysis/MMcCalibrationCalc.cc (modified) (1 diff)

trunk/MagicSoft/Mars/Changelog

@@ -4 +4 @@
 
                                                  -*-*- END OF LINE -*-*-
+
+ 2004/01/14: Abelardo Moralejo
+
+   * manalysis/MMcCalibrationCalc.cc:
+     - set pedestal rms per FADC slice to 0.01 counts in case it is 0.
+       This can happen (depending on camera simulation parameters), and
+       rms = 0 later resulted in the pixel being eliminated in the 
+       image cleaning.
 
  2004/01/14: Thomas Bretz

trunk/MagicSoft/Mars/manalysis/MMcCalibrationCalc.cc

@@ -216 +216 @@
           fSignalCam->GetNumUsedHiGainFADCSlices()*fHeaderFadc->GetPedestal(i);
 
+        //
+        // In some cases, depending on the camera simulation parameters, one can have
+        // very little or no noise in the FADC. In the case the rms of pedestal is zero, 
+        // the pixel will be cleaned out later in the image cleaning. To avoid this problem,
+        // we set a default value of 0.01 ADC counts for the RMS per slice:
+        //
+
         const Float_t pedestrms  = sigpix.IsLoGainUsed()? 
-          sqrt((Double_t)(fSignalCam->GetNumUsedLoGainFADCSlices())) * fHeaderFadc->GetPedestalRmsLow(i) : 
-          sqrt((Double_t)(fSignalCam->GetNumUsedHiGainFADCSlices())) * fHeaderFadc->GetPedestalRmsHigh(i);
+          sqrt((Double_t)(fSignalCam->GetNumUsedLoGainFADCSlices())) * TMath::Max(0.01, fHeaderFadc->GetPedestalRmsLow(i)) : 
+          sqrt((Double_t)(fSignalCam->GetNumUsedHiGainFADCSlices())) * TMath::Max(0.01, fHeaderFadc->GetPedestalRmsHigh(i));
 
         calpix.SetPedestal(pedestmean, pedestrms);