Changeset 7895 for trunk/MagicSoft/Mars/msignal

Timestamp:

08/21/06 14:38:52 (18 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/msignal

Files:

• MExtractTimeAndChargeDigitalFilter.cc (modified) (7 diffs)
• MExtractTimeAndChargeSpline.cc (modified) (2 diffs)

trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeDigitalFilter.cc

@@ -239 +239 @@
     // here for speed reason. The low-gain has a changing extraction offset,
     // so it will be added at every event (in FindTimeAndChargeLoGain)
-    fTimeShiftHiGain = 0.5 + 1./fBinningResolutionHiGain + fHiGainFirst;
+    fTimeShiftHiGain = 0.5 + 1./fBinningResolutionHiGain;
     fTimeShiftLoGain = 0.5 + 1./fBinningResolutionLoGain;
 
@@ -404 +404 @@
   
   Float_t time_sum  = 0.;
-  Float_t fmax      = 0.;
+  Float_t fmax      = -FLT_MAX;
   Float_t ftime_max = 0.;
-  Int_t   max_p     = 0;
+  Int_t   max_p     = -1;
   
   //
@@ -436 +436 @@
     } /*   for (Int_t i=0;i<range-fWindowSizeHiGain+1;i++) */
 
+  if (max_p<0)
+  {
+      sum = 0;
+      time = -1;
+      dsum = -1;
+      dtime = -1;
+      return;
+  }
+
   if (fmax==0)
+  {
+      sum = 0;
+      time = -1;
+      dsum = 0;
+      dtime = 0;
       return;
+  }
 
   ftime_max        /= fmax;
@@ -479 +494 @@
 
   if (sum == 0)
+  {
+      time = -1;
+      dsum = 0;
+      dtime = 0;
       return;
+  }
 
   // here, the first high-gain slice is already included in the fTimeShiftHiGain
@@ -544 +564 @@
 
   Float_t time_sum  = 0.;
-  Float_t fmax      = 0.;
+  Float_t fmax      = -FLT_MAX;
   Float_t ftime_max = 0.;
-  Int_t   max_p     = 0;
+  Int_t   max_p     = -1;
   
   //
@@ -576 +596 @@
     } /*   for (Int_t i=0;i<range-fWindowSizeLoGain+1;i++) */
 
-  time = 0;
+  if (max_p<0)
+  {
+      sum = 0;
+      time = -1;
+      dsum = -1;
+      dtime = -1;
+      return;
+  }
+
   if (fmax==0)
+  {
+      sum = 0;
+      time = -1;
+      dsum = 0;
+      dtime = 0;
       return;
+  }
 
   ftime_max        /= fmax;
@@ -617 +651 @@
 
   if (sum == 0)
+  {
+      time = -1;
+      dsum = 0;
+      dtime = 0;
       return;
+  }
 
   time = max_p + fTimeShiftLoGain + (Float_t)fLoGainFirst /* this shifts the time to the start of the rising edge */

trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeSpline.cc

@@ -374 +374 @@
 
   sat = 0;
+  dsum = 0; // In all cases the extracted signal is valid
 
   const Float_t pedes  = ped.GetPedestal();
@@ -741 +742 @@
   Int_t  max    = -9999;
 
+  dsum = 0; // In all cases the extracted signal is valid
+
   //
   // Check for saturation in all other slices