Changeset 6166 for trunk/MagicSoft/Mars

Timestamp:

02/01/05 14:16:45 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeDigitalFilterPeakSearch.cc (modified) (10 diffs)

trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeDigitalFilterPeakSearch.cc

@@ -81 +81 @@
 const Byte_t MExtractTimeAndChargeDigitalFilterPeakSearch::fgLoGainFirst             =  0;
 const Byte_t MExtractTimeAndChargeDigitalFilterPeakSearch::fgLoGainLast              = 14;
-const Byte_t MExtractTimeAndChargeDigitalFilterPeakSearch::fgOffsetLeftFromPeak      =  1;
+const Byte_t MExtractTimeAndChargeDigitalFilterPeakSearch::fgOffsetLeftFromPeak      =  3;
 const Byte_t MExtractTimeAndChargeDigitalFilterPeakSearch::fgOffsetRightFromPeak     =  3;
 const Byte_t MExtractTimeAndChargeDigitalFilterPeakSearch::fgPeakSearchWindowSize    =  2;
@@ -102 +102 @@
     fTitle = title ? title : "Digital Filter";
 
+    SetRange(fgHiGainFirst, fgHiGainLast, fgLoGainFirst, fgLoGainLast);
+
     SetOffsetLeftFromPeak();
     SetOffsetRightFromPeak();
@@ -137 +139 @@
   fHiGainOutOfRangeLeft  = 0;
   fHiGainOutOfRangeRight = 0;
-  fLoGainOutOfRangeLeft   = 0;
-  fLoGainOutOfRangeRight  = 0;
+  fLoGainOutOfRangeLeft  = 0;
+  fLoGainOutOfRangeRight = 0;
   
   return kTRUE;
@@ -152 +154 @@
 // saturated slices.
 //
-void MExtractTimeAndChargeDigitalFilterPeakSearch::FindPeak(Byte_t *ptr, Byte_t &startslice, Int_t &max, 
+void MExtractTimeAndChargeDigitalFilterPeakSearch::FindPeak(Byte_t *ptr, Byte_t *logain, Byte_t &startslice, Int_t &max, 
                                                             Int_t &sat, Byte_t &satpos) const
 {
 
-  const Byte_t *end = ptr + fHiGainLast - fHiGainFirst + 1;
+  Byte_t *end = ptr + fHiGainLast - fHiGainFirst + 1;
 
   sat = 0;
@@ -208 +210 @@
     }
 
+  if (!fHiLoLast)
+    return;
+
+  Byte_t *l = logain;
+
+  while (++p < end && l < logain+fHiLoLast)
+    {
+      sum += *l - *p;
+      if (sum > max)
+        {
+          max = sum;
+          startslice = p-ptr+1;
+        }
+
+      if (*l++ >= fSaturationLimit)
+        {
+          if (sat == 0)
+            satpos = l-logain + fHiGainLast - fHiGainFirst;
+          sat++;
+        }
+    }
+  
+  if (fHiLoLast <= fPeakSearchWindowSize)
+    return;
+
+  while (l < logain+fHiLoLast)
+    {
+      if (*l++ >= fSaturationLimit)
+        {
+          if (sat == 0)
+            satpos = l-logain+fHiGainLast-fHiGainFirst;
+          sat++;
+        }
+
+      sum += *l - *(l-fPeakSearchWindowSize);
+      if (sum > max)
+        {
+          max = sum;
+          startslice = l-logain + fHiGainLast - fHiGainFirst - fPeakSearchWindowSize + 1;
+        }
+    }
+
   return;
 }
@@ -239 +283 @@
   Byte_t loGainFirstsave = fLoGainFirst;
   Byte_t loGainLastsave  = fLoGainLast; 
+  Byte_t hiLoLastsave    = fHiLoLast;
 
   Byte_t higainfirst     = fHiGainFirst;
@@ -248 +293 @@
       sat = 0;
 
-      FindPeak(pixel.GetHiGainSamples()+fHiGainFirst, startslice, sumhi, sat, satpos);
+      FindPeak(pixel.GetHiGainSamples()+fHiGainFirst, pixel.GetLoGainSamples(), startslice, sumhi, sat, satpos);
 
       if (sumhi > maxsumhi && sat == 0)
@@ -279 +324 @@
   // Shift the last slice to the right:
   //
-  if (higainfirst + fOffsetRightFromPeak + fWindowSizeHiGain <= hiGainLastsave)
-    fHiGainLast  = higainfirst + fOffsetRightFromPeak + fWindowSizeHiGain;
+  const Byte_t rlim = higainfirst + fOffsetRightFromPeak + fWindowSizeHiGain;
+  if (rlim <= hiGainLastsave+fHiLoLast)
+    if (rlim > hiGainLastsave)
+      fHiLoLast   = rlim - fHiGainLast;
+    else
+      {
+        fHiLoLast   = 0;
+        fHiGainLast = rlim;
+      }
   else 
     fHiGainOutOfRangeRight++;
       
 
-  if ( fHiGainFirst+(Int_t)fOffsetLoGain >= fLoGainFirst ) 
-    fLoGainFirst = fHiGainFirst + (Int_t)fOffsetLoGain;
+  const Byte_t llim = fHiGainFirst + (Int_t)fOffsetLoGain;
+  if ( llim >= fLoGainFirst ) 
+    fLoGainFirst = llim;
   else
     fLoGainOutOfRangeLeft++;
@@ -302 +355 @@
   while (pixel.Next())
     {
+
       //
       // Find signal in hi- and lo-gain
@@ -384 +438 @@
   fLoGainFirst  = loGainFirstsave;
   fLoGainLast   = loGainLastsave ; 
+  fHiLoLast     = hiLoLastsave;
 
   return kTRUE;