Changeset 6366 for trunk/MagicSoft/Mars/msignal

Timestamp:

02/10/05 19:28:16 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/msignal

Files:

• MExtractTimeAndChargeSpline.cc (modified) (3 diffs)
• MExtractTimeAndChargeSpline.h (modified) (1 diff)

trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeSpline.cc

@@ -174 +174 @@
 //
 MExtractTimeAndChargeSpline::MExtractTimeAndChargeSpline(const char *name, const char *title) 
-    : fAbMax(0.), fAbMaxPos(0.), fHalfMax(0.), fRandomIter(0)
+    : fAbMax(0.), fAbMaxPos(0.), fHalfMax(0.), 
+      fRiseTimeHiGain(fgRiseTimeHiGain), fFallTimeHiGain(fgFallTimeHiGain),
+      fRandomIter(0)
 {
 
@@ -181 +183 @@
 
   SetResolution();
-
-  SetRiseTimeHiGain();
-  SetFallTimeHiGain();
-  
+  SetLoGainStretch();
+  SetOffsetLoGain(fgOffsetLoGain);
+
   SetChargeType();
-  SetLoGainStretch();
-
-  SetOffsetLoGain(fgOffsetLoGain);
-
   SetRange(fgHiGainFirst, fgHiGainLast, fgLoGainFirst, fgLoGainLast);
 }
@@ -229 +226 @@
   SETBIT(fFlags,typ);
 
+  if (IsExtractionType(kAmplitude))
+    {
+      fNumHiGainSamples = 1.;
+      fNumLoGainSamples = fLoGainLast ? 1. : 0.; 
+      fSqrtHiGainSamples = 1.;
+      fSqrtLoGainSamples = 1.;
+      fWindowSizeHiGain  = 1;
+      fWindowSizeLoGain  = 1;
+      fRiseTimeHiGain    = 0.5;
+      
+      return;
+    }
+
+  if (IsExtractionType(kIntegral))
+    {
+
+      fNumHiGainSamples  = fRiseTimeHiGain + fFallTimeHiGain;
+      fNumLoGainSamples  = fLoGainLast ? fRiseTimeLoGain + fFallTimeLoGain : 0.;
+      //      fNumLoGainSamples  *= 0.75;      
+
+      fSqrtHiGainSamples = TMath::Sqrt(fNumHiGainSamples);
+      fSqrtLoGainSamples = TMath::Sqrt(fNumLoGainSamples);
+      fWindowSizeHiGain  = (Int_t)(fRiseTimeHiGain + fFallTimeHiGain);
+      fWindowSizeLoGain  = (Int_t)(fRiseTimeLoGain + fFallTimeLoGain);
+      //      fNumLoGainSamples  *= 0.75;      
+    }
 }
 

trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeSpline.h

@@ -76 +76 @@
   void SetRange      ( Byte_t hifirst=0, Byte_t hilast=0, Byte_t lofirst=0, Byte_t lolast=0 );  
   void SetResolution ( const Float_t f=fgResolution  )  { fResolution  = f;  }
-  void SetRiseTimeHiGain( const Float_t f=fgRiseTimeHiGain    )  { fRiseTimeHiGain = f;  }
-  void SetFallTimeHiGain( const Float_t f=fgFallTimeHiGain    )  { fFallTimeHiGain = f;  }
+  void SetRiseTimeHiGain( const Float_t f=fgRiseTimeHiGain    )
+    {
+      fRiseTimeHiGain    = f;
+      fRiseTimeLoGain    = f*fLoGainStretch;
+      fSqrtHiGainSamples = TMath::Sqrt(fNumHiGainSamples);
+      fWindowSizeHiGain  = (Int_t)(fRiseTimeHiGain + fFallTimeHiGain);
+    }
+  void SetFallTimeHiGain( const Float_t f=fgFallTimeHiGain    )
+    {
+      fFallTimeHiGain    = f;
+      fFallTimeLoGain    = f*fLoGainStretch;
+      fNumHiGainSamples  = fRiseTimeHiGain + fFallTimeHiGain;
+      fNumLoGainSamples  = fLoGainLast ? fRiseTimeLoGain + fFallTimeLoGain : 0.;
+      fSqrtLoGainSamples = TMath::Sqrt(fNumLoGainSamples);
+      fWindowSizeLoGain  = (Int_t)(fRiseTimeLoGain + fFallTimeLoGain);
+    }
+
   void SetLoGainStretch ( const Float_t f=fgLoGainStretch    )  { fLoGainStretch = f;   }