Changeset 5293 for trunk/MagicSoft/Mars/msignal

Timestamp:

10/20/04 11:40:32 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeSpline.cc (modified) (4 diffs)

trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeSpline.cc

@@ -341 +341 @@
 // --------------------------------------------------------------------------
 //
-// Calculates the arrival time for each pixel 
+// Calculates the arrival time and charge for each pixel 
 //
 void MExtractTimeAndChargeSpline::FindTimeAndChargeHiGain(Byte_t *first, Byte_t *logain, Float_t &sum, Float_t &dsum, 
@@ -443 +443 @@
 
   for (Int_t k=range-2;k>=0;k--)
-    fHiGainSecondDeriv[k] = (fHiGainSecondDeriv[k]*fHiGainSecondDeriv[k+1] + fHiGainFirstDeriv[k])/6.;
+    fHiGainSecondDeriv[k] = fHiGainSecondDeriv[k]*fHiGainSecondDeriv[k+1] + fHiGainFirstDeriv[k];
+  for (Int_t k=range-2;k>=0;k--)
+    fHiGainSecondDeriv[k] /= 6.;
   
   //
@@ -711 +713 @@
 // --------------------------------------------------------------------------
 //
-// Calculates the arrival time for each pixel 
+// Calculates the arrival time and charge for each pixel 
 //
 void MExtractTimeAndChargeSpline::FindTimeAndChargeLoGain(Byte_t *first, Float_t &sum, Float_t &dsum, 
@@ -791 +793 @@
   fLoGainSecondDeriv[range-1] = 0.;
   for (Int_t k=range-2;k>=0;k--)
-    fLoGainSecondDeriv[k] = (fLoGainSecondDeriv[k]*fLoGainSecondDeriv[k+1] + fLoGainFirstDeriv[k])/6.;
+    fLoGainSecondDeriv[k] = fLoGainSecondDeriv[k]*fLoGainSecondDeriv[k+1] + fLoGainFirstDeriv[k];
+  for (Int_t k=range-2;k>=0;k--)
+    fLoGainSecondDeriv[k] /= 6.;
   
   //