Changeset 4284 for trunk

Timestamp:

06/08/04 18:09:34 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/msignal

Files:

• MExtractTimeFastSpline.cc (modified) (5 diffs)
• MExtractTimeFastSpline.h (modified) (3 diffs)

trunk/MagicSoft/Mars/msignal/MExtractTimeFastSpline.cc

@@ -53 +53 @@
 const Byte_t  MExtractTimeFastSpline::fgLoGainLast   = 14;
 const Float_t MExtractTimeFastSpline::fgResolution   = 0.003;
+const Float_t MExtractTimeFastSpline::fgRiseTime     = 1.5;
 // --------------------------------------------------------------------------
 //
@@ -72 +73 @@
 
   SetResolution();
+  SetRiseTime  ();
   SetRange(fgHiGainFirst, fgHiGainLast, fgLoGainFirst, fgLoGainLast);
 
@@ -144 +146 @@
 
     }
+
 }
 
@@ -154 +157 @@
                                         Byte_t &sat, const MPedestalPix &ped) const
 {
-  
+
   const Int_t range = fHiGainLast - fHiGainFirst + 1;
   const Byte_t *end = first + range;
@@ -352 +355 @@
   // 
   klo   = maxpos;
-  while (klo > maxpos-4)
+  while (klo > maxpos-fStartBeforeMax)
     {
       if (*(first+klo) < (Byte_t)halfmax)

trunk/MagicSoft/Mars/msignal/MExtractTimeFastSpline.h

@@ -15 +15 @@
   static const Byte_t  fgLoGainLast;     // Default for fLoGainLast   (now set to: 14)
   static const Float_t fgResolution;     // Default for fResolution   (now set to: 0.003)
+  static const Float_t fgRiseTime  ;     // Default for fRiseTime     (now set to: 1.5  C  
 
   Float_t *fHiGainFirstDeriv;
@@ -22 +23 @@
 
   Float_t fResolution;                        // The time resolution in FADC units
+  Float_t fRiseTime  ;                        // The rise time of the pulse 
+  Byte_t  fStartBeforeMax;                    // Slices to start searching for the halfmax before max
   
   void FindTimeHiGain(Byte_t *first, Float_t &time, Float_t &dtime, Byte_t &sat, const MPedestalPix &ped) const;
@@ -33 +36 @@
   void SetRange(Byte_t hifirst=0, Byte_t hilast=0, Byte_t lofirst=0, Byte_t lolast=0);      
   void SetResolution(Float_t f=fgResolution)     { fResolution = f;  }
+  void SetRiseTime  (Float_t f=fgRiseTime  )     { fRiseTime   = f; fStartBeforeMax = (Int_t)(2*fRiseTime); }  
   
   ClassDef(MExtractTimeFastSpline, 0)   // Task to Extract the Arrival Times using a Fast Spline