Changeset 8564 for trunk/MagicSoft/Mars/mpedestal/MExtractPedestal.cc

Timestamp:

06/16/07 22:59:26 (17 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mpedestal/MExtractPedestal.cc (modified) (14 diffs)

trunk/MagicSoft/Mars/mpedestal/MExtractPedestal.cc

@@ -1 +1 @@
 /* ======================================================================== *\
-! $Name: not supported by cvs2svn $:$Id: MExtractPedestal.cc,v 1.31 2007-05-16 13:56:17 tbretz Exp $
+! $Name: not supported by cvs2svn $:$Id: MExtractPedestal.cc,v 1.32 2007-06-16 21:59:19 tbretz Exp $
 ! --------------------------------------------------------------------------
 !
@@ -179 +179 @@
 const TString  MExtractPedestal::fgNameRawEvtData  = "MRawEvtData";
 
-const UShort_t MExtractPedestal::fgCheckWinFirst   =  0;
-const UShort_t MExtractPedestal::fgCheckWinLast    = 29;
-const UShort_t MExtractPedestal::fgMaxSignalVar    = 40;
+const UShort_t MExtractPedestal::fgCheckWinFirst   =   0;
+const UShort_t MExtractPedestal::fgCheckWinLast    =  29;
+const UShort_t MExtractPedestal::fgMaxSignalVar    =  40;
+const UShort_t MExtractPedestal::fgMaxSignalAbs    = 250;
 
 // --------------------------------------------------------------------------
@@ -210 +211 @@
     SetCheckRange(fgCheckWinFirst, fgCheckWinLast);
     SetMaxSignalVar(fgMaxSignalVar);
+    SetMaxSignalAbs(fgMaxSignalAbs);
 
     Clear();
@@ -598 +600 @@
 //
 // Check whether the signal variation between fCheckWinFirst and fCheckWinLast
-// exceeds fMaxSignalVar or the signal is greater than 250
+// exceeds fMaxSignalVar or the signal is greater than fMaxSignalAbs
 //
 Bool_t MExtractPedestal::CheckVariation(UInt_t idx) const
 {
     // This is the fast workaround to put hi- and lo-gains together
-    Byte_t *slices = fSignal->GetSamplesRaw(idx);
+    USample_t *slices = fSignal->GetSamplesRaw(idx);
 
     // Start 'real' work
-    UInt_t max = 0;
-    UInt_t min = (UInt_t)-1;
+    USample_t max = 0;
+    USample_t min = (USample_t)-1;
 
     // Find the maximum and minimum signal per slice in the high gain window
-    for (Byte_t *slice=slices+fCheckWinFirst; slice<=slices+fCheckWinLast; slice++)
+    for (USample_t *slice=slices+fCheckWinFirst; slice<=slices+fCheckWinLast; slice++)
     {
         if (*slice > max)
@@ -620 +622 @@
     // If the maximum in the high gain window is smaller than
     // FIXME: Precompiled value!
-    return max-min<fMaxSignalVar && max<250;
+    return max-min<fMaxSignalVar && max<fMaxSignalAbs;
 }
 
@@ -648 +650 @@
 
     // Do some handling if maxpos is last slice?
-    const Int_t maxposhi = fSignal->GetMaxPos(idx, start, start+range-1);
+    const Int_t maxposhi = fRandomCalculation ? 0 : fSignal->GetMaxPos(idx, start, start+range-1);
 
     const Float_t *sig = fSignal->GetSamples(idx);
@@ -659 +661 @@
                                          dummy[0], dummy[1], dummy[2],
                                          0, maxposhi);
+
     return sum;
 }
@@ -672 +675 @@
     const Int_t first = fExtractWinFirst+offset;
 
-    Byte_t *ptr = fSignal->GetSamplesRaw(pixel.GetPixelId())+first;
-    Byte_t *end = ptr + fExtractWinSize;
+    USample_t *ptr = fSignal->GetSamplesRaw(pixel.GetPixelId())+first;
+    USample_t *end = ptr + fExtractWinSize;
 
     Int_t abflag = pixel.HasABFlag() + first;
@@ -786 +789 @@
     // multiplies by this number
 
+    // scale to 256
+    const UInt_t scale = fExtractor ? 1 : fRunHeader->GetScale();
+
+    ped    /= scale;
+    abOffs /= scale;
+
     // 5. Calculate the RMS from the Variance:
-    const Double_t rms = var<0 ? 0 : TMath::Sqrt(var);
+    const Double_t rms = var<0 ? 0 : TMath::Sqrt(var)/scale;
 
     // abOffs contains only half of the signal as ped.
@@ -834 +843 @@
     // multiplies by this number
 
+    // scale to 256
+    const UInt_t scale = fExtractor ? 1 : fRunHeader->GetScale();
+
     // 5. Scale the mean, variance and AB-noise to the number of pixels:
-    ped    /= napix;
-    var    /= napix;
-    abOffs /= napix;
+    ped    /= napix*scale;
+    abOffs /= napix*scale;
 
     // 6. Calculate the RMS from the Variance:
-    const Double_t rms = var<0 ? 0 : TMath::Sqrt(var);
+    const Double_t rms = var<0 ? 0 : TMath::Sqrt(var)/scale;
 
     // abOffs contains only half of the signal as ped.
@@ -887 +898 @@
     // multiplies by this number
 
+    // scale to 256
+    const UInt_t scale = fExtractor ? 1 : fRunHeader->GetScale();
+
     // 5. Scale the mean, variance and AB-noise to the number of pixels:
-    ped    /= nspix;
-    var    /= nspix;
-    abOffs /= nspix;
+    ped    /= nspix*scale;
+    abOffs /= nspix*scale;
 
     // 6. Calculate the RMS from the Variance:
-    const Double_t rms = var<0 ? 0 : TMath::Sqrt(var);
+    const Double_t rms = var<0 ? 0 : TMath::Sqrt(var)/scale;
 
     // abOffs contains only half of the signal as ped.
@@ -947 +960 @@
     *fLog << "CheckWindow from slice " << fCheckWinFirst   << " to " << fCheckWinLast << " incl." << endl;
     *fLog << "Max.allowed signal variation: " << fMaxSignalVar << endl;
+    *fLog << "Max.allowed signal absolute:  " << fMaxSignalAbs << endl;
 }
 
@@ -1020 +1034 @@
     }
 
+    if (IsEnvDefined(env, prefix, "MaxSignalAbs", print))
+    {
+        SetMaxSignalAbs(GetEnvValue(env, prefix, "MaxSignalAbs", fMaxSignalAbs));
+        rc = kTRUE;
+    }
+
+    // find resource for MPedestalCam
     if (IsEnvDefined(env, prefix, "NamePedestalCamInter", print))
     {
@@ -1034 +1055 @@
     return rc;
 }
-