Changeset 3032 for trunk/MagicSoft/Mars/manalysis/MArrivalTimeCalc.cc

Timestamp:

02/05/04 17:37:02 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/manalysis/MArrivalTimeCalc.cc (modified) (4 diffs)

trunk/MagicSoft/Mars/manalysis/MArrivalTimeCalc.cc

@@ -61 +61 @@
 // Default constructor.
 //
+// Initialize step size by default to 0.03 time slices == 100 ps.
+// 
 MArrivalTimeCalc::MArrivalTimeCalc(const char *name, const char *title)
-    : fStepSize(0.1)
+    : fStepSize(0.03) 
 {
 
@@ -80 +82 @@
 //  - MArrivalTime
 //
+
 Int_t MArrivalTimeCalc::PreProcess(MParList *pList)
 {
+
+
     fRawEvt = (MRawEvtData*)pList->FindObject(AddSerialNumber("MRawEvtData"));
     if (!fRawEvt)
@@ -186 +191 @@
 
   //
-  // Now find the maximum evaluating the spline function at every 1/10 time slice
+  // Now find the half maximum (!) 
+  // evaluating the spline function at every fStepSize time slice
   //
   Double_t abscissa=0;
@@ -201 +207 @@
           maxAb  = abscissa;
         }
+      // make step size a bit bigger first
       abscissa += fStepSize;
-    }
-
-  return (Float_t)maxAb;
-}
-
+      //      abscissa += fStepSize;
+    }
+
+  //
+  // another (much smaller) loop to move back from the maximum 
+  //
+  Double_t halfMaxAb = 0;
+
+  abscissa = maxAb;
+  
+  while (abscissa > 0)
+    {
+      
+      const Double_t swap = spline.Eval(abscissa);
+      
+      if (swap < maxOrd/2.)
+        {
+          halfMaxAb = abscissa;
+          break;
+        }
+      
+      abscissa -= fStepSize;
+    }
+  
+
+  //  return (Float_t)maxAb;
+  return (Float_t)halfMaxAb;
+}
+