Changeset 4752 for trunk/MagicSoft/Mars/msignal

Timestamp:

08/27/04 11:36:24 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/msignal/MExtractTimeFastSpline.cc (modified) (9 diffs)

trunk/MagicSoft/Mars/msignal/MExtractTimeFastSpline.cc

@@ -200 +200 @@
       pp = fHiGainSecondDeriv[i-1] + 4.;
       fHiGainSecondDeriv[i] = -1.0/pp;
-      fHiGainFirstDeriv [i] = *(p+1) - 2.* *(p) + *(p-1);
-      fHiGainFirstDeriv [i] = (6.0*fHiGainFirstDeriv[i]-fHiGainFirstDeriv[i-1])/pp;
+      const Double_t deriv = *(p+1) - 2.* *(p) + *(p-1);
+      fHiGainFirstDeriv [i] = (6.0*deriv-fHiGainFirstDeriv[i-1])/pp;
     }
 
@@ -232 +232 @@
   // interval maxpos+1.
   //
+
+  Float_t higainklo = fHiGainSecondDeriv[klo];
+  Float_t higainkhi = fHiGainSecondDeriv[khi];
   while (x<upper-0.3)
     {
@@ -241 +244 @@
       y = a*klocont
         + b*khicont
-        + (a*a*a-a)*fHiGainSecondDeriv[klo] 
-        + (b*b*b-b)*fHiGainSecondDeriv[khi];
+        + (a*a*a-a)*higainklo
+        + (b*b*b-b)*higainkhi;
       
       if (y > abmax)
@@ -265 +268 @@
       khicont = (Float_t)*(first+khi);
 
+      higainklo = fHiGainSecondDeriv[klo];
+      higainkhi = fHiGainSecondDeriv[khi];
       while (x<upper-0.3)
         {
@@ -274 +279 @@
           y = a* klocont
             + b* khicont
-            + (a*a*a-a)*fHiGainSecondDeriv[klo] 
-            + (b*b*b-b)*fHiGainSecondDeriv[khi];
+            + (a*a*a-a)*higainklo
+            + (b*b*b-b)*higainkhi;
           
           if (y > abmax)
@@ -295 +300 @@
   step  = 0.04; // step size of 83 ps 
 
+  higainklo = fHiGainSecondDeriv[klo];
+  higainkhi = fHiGainSecondDeriv[khi];
   while (x<up)
     {
@@ -304 +311 @@
       y = a* klocont
         + b* khicont
-        + (a*a*a-a)*fHiGainSecondDeriv[klo] 
-        + (b*b*b-b)*fHiGainSecondDeriv[khi];
+        + (a*a*a-a)*higainklo
+        + (b*b*b-b)*higainkhi;
       
       if (y > abmax)
@@ -329 +336 @@
   b     = x - lower;
 
+  higainklo = fHiGainSecondDeriv[klo];
+  higainkhi = fHiGainSecondDeriv[khi];
   while (x>lo)
     {
@@ -338 +347 @@
       y = a* klocont
         + b* khicont
-        + (a*a*a-a)*fHiGainSecondDeriv[klo] 
-        + (b*b*b-b)*fHiGainSecondDeriv[khi];
+        + (a*a*a-a)*higainklo
+        + (b*b*b-b)*higainkhi;
       
       if (y > abmax)