Changeset 8519 for trunk/MagicSoft/Mars/msignal/MExtractTimeAndCharge.cc

Timestamp:

05/16/07 14:56:18 (18 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/msignal/MExtractTimeAndCharge.cc (modified) (6 diffs)

trunk/MagicSoft/Mars/msignal/MExtractTimeAndCharge.cc

@@ -1 +1 @@
 /* ======================================================================== *\
-! $Name: not supported by cvs2svn $:$Id: MExtractTimeAndCharge.cc,v 1.64 2007-03-04 13:25:07 tbretz Exp $
+! $Name: not supported by cvs2svn $:$Id: MExtractTimeAndCharge.cc,v 1.65 2007-05-16 13:56:17 tbretz Exp $
 ! --------------------------------------------------------------------------
 !
@@ -189 +189 @@
 // order polynomial interpolation.
 //
-Double_t MExtractTimeAndCharge::GetSaturationTime(Int_t sat0, const Float_t *sig, Int_t maxcont) const
+Double_t MExtractTimeAndCharge::GetSaturationTime(Int_t sat0, const Float_t *sig, Int_t maxpos) const
 {
     const Int_t p = sat0>1 ? sat0-2 : sat0-1;
     if (sat0<=0)
         return 0;
+
+    const Float_t &maxcont = sig[maxpos];
     if (sat0==1)
         return sig[0]>maxcont/2 ? 0 : 0.5;
@@ -228 +230 @@
         const Float_t *sig = fSignal->GetSamples(pixidx);
 
-        Int_t sathi0 = fHiGainFirst;  // First slice to extract and first saturating slice
-        Int_t sathi1 = fHiGainLast;   // Last  slice to extract and last saturating slice
-
-        UInt_t maxcont;
-        Int_t maxposhi = fSignal->GetMax(pixidx, sathi0, sathi1, maxcont);
+        const UInt_t maxcont  = fSignal->GetRawMaxVal(pixidx, fHiGainFirst, fHiGainLast);
+        const  Int_t maxposhi = fSignal->GetMaxPos(pixidx, fHiGainFirst, fHiGainLast);
+
         // Would it be better to take lastsat-firstsat?
+        Int_t sathi0   = fHiGainFirst;  // First slice to extract and first saturating slice
+        Int_t sathi1   = fHiGainLast;   // Last  slice to extract and last saturating slice
         Int_t numsathi = fSignal->GetSaturation(pixidx, fSaturationLimit, sathi0, sathi1);
 
@@ -253 +255 @@
         if (numsathi>1)
         {
-            timehi = GetSaturationTime(sathi0, sig, maxcont)-fHiGainFirst;
+            timehi = GetSaturationTime(sathi0, sig, fHiGainFirst+maxposhi)-fHiGainFirst;
             deltatimehi = 0;
         }
@@ -335 +337 @@
                 first = numh+fLoGainFirst;
              */
+            // Would it be better to take lastsat-firstsat?
+            const Int_t maxposlo = fSignal->GetMaxPos(pixidx, first, last);
+
             Int_t satlo0 = first;   // First slice to extract and first saturating slice
             Int_t satlo1 = last;    // Last  slice to extract and last saturating slice
-
-            // Would it be better to take lastsat-firstsat?
-            UInt_t maxlo;
-            Int_t maxposlo = fSignal->GetMax(pixidx, satlo0, satlo1, maxlo);
             numsatlo = fSignal->GetSaturation(pixidx, fSaturationLimit, satlo0, satlo1);
 
@@ -359 +360 @@
             if (numsatlo>1)
             {
-                timelo = GetSaturationTime(satlo0, sig, maxlo)-numh-first;
+                timelo = GetSaturationTime(satlo0, sig, first+maxposlo)-numh-first;
                 deltatimelo = 0;
             }