Changeset 5743 for trunk/MagicSoft

Timestamp:

01/08/05 12:32:08 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeDigitalFilter.cc (modified) (5 diffs)

trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeDigitalFilter.cc

@@ -260 +260 @@
       end = logain + fHiLoLast;
 
-      ids = fHiGainFirst;
+      ids = fHiGainFirst + range;
       sample = fHiGainSignal.GetArray();
       while (logain<end)
@@ -266 +266 @@
 
           *sample++ = (Float_t)*logain - pedmean[(ids++ + abflag) & 0x1];
-          range++;
 
           if (*logain > max)
@@ -276 +275 @@
           if (*logain++ >= fSaturationLimit)
             sat++;
-        }
-    }
-
-  //
-  // allow one saturated slice 
+
+          range++;
+        }
+    }
+
+  //
+  // allow no saturated slice 
   //
   if (sat > 0)
@@ -291 +292 @@
   if (IsNoiseCalculation())
     {
-      if (fRandomIter == 10)
+      if (fRandomIter == fBinningResolution)
         fRandomIter = 0;
       for (Int_t ids=0; ids < fWindowSizeHiGain; ids++)
-      {
+        {
           const Int_t   idx = fBinningResolutionHiGain*ids + fRandomIter;
           sum              += fAmpWeightsHiGain [idx]*fHiGainSignal[ids]; 
@@ -425 +426 @@
 
   //
-  // Slide with a window of size fWindowSizeHiGain over the sample 
+  // Slide with a window of size fWindowSizeLoGain over the sample 
   // and multiply the entries with the corresponding weights
   //
   if (IsNoiseCalculation())
     {
+      if (fRandomIter == fBinningResolutionLoGain)
+        fRandomIter = 0;
       for (Int_t ids=0; ids < fWindowSizeLoGain; ids++)
         {
-          const Int_t   idx = fBinningResolutionHiGain*ids + fRandomIter;
+          const Int_t   idx = fBinningResolutionLoGain*ids + fRandomIter;
           sum              += fAmpWeightsLoGain [idx]*fLoGainSignal[ids]; 
         }