Changeset 2640 for trunk/MagicSoft/Mars/manalysis/MExtractSignal.cc

Timestamp:

12/11/03 13:58:40 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/manalysis/MExtractSignal.cc (modified) (5 diffs)

trunk/MagicSoft/Mars/manalysis/MExtractSignal.cc

@@ -57 +57 @@
 //
 MExtractSignal::MExtractSignal(const char *name, const char *title)
-  : fSaturationLimit(255),
-    fConversionHiLo(10.)
+  : fBefore(2), fAfter(4), fSaturationLimit(254), fConversionHiLo(10.)
 {
 
@@ -123 +122 @@
     }
 
-    fNumHiGainSamples =  fRunHeader->GetNumSamplesHiGain();
-    fNumLoGainSamples =  fRunHeader->GetNumSamplesLoGain();
-    fSqrtHiGainSamples =  TMath::Sqrt((float)fNumHiGainSamples);
-    fSqrtLoGainSamples =  TMath::Sqrt((float)fNumLoGainSamples);
+    //fNumHiGainSamples =  fRunHeader->GetNumSamplesHiGain();
+    //fNumLoGainSamples =  fRunHeader->GetNumSamplesLoGain();
+    fNumHiGainSamples =  fAfter+fBefore;
+    fNumLoGainSamples =  fAfter+fBefore;
+
+    fSqrtHiGainSamples = TMath::Sqrt((float)fNumHiGainSamples); 
+    fSqrtLoGainSamples = TMath::Sqrt((float)fNumLoGainSamples); 
 
     //
@@ -156 +158 @@
         UShort_t satHi = 0;
         UShort_t satLo = 0;
-        Float_t  sum   = 0.;
-        Float_t sumerr = 0.;
+
+        Byte_t *ptr   = pixel.GetHiGainSamples();
+        Byte_t *max   = ptr + pixel.GetIdxMaxHiGainSample();
+        Byte_t *end   = ptr + fRawEvt->GetNumHiGainSamples();
+        Byte_t *first;
+        Byte_t *last;
+
+        UInt_t  sumHi   = 0;
+        UInt_t  sumLo   = 0;
 
         const Int_t pixid = pixel.GetPixelId();
@@ -165 +174 @@
         const Float_t pedes   = ped.GetPedestal();
         const Float_t pedrms  = ped.GetPedestalRms();
+
+        if (max-3 < ptr)
+          {
+            first = ptr+3;
+            last  = ptr+3+fAfter+fBefore;
+          }
+        else if (max+fAfter-1 > end)
+          {
+            first = end-fAfter-fBefore-1;
+            last  = end-1;
+          }
+        else
+          {
+            first = max-fBefore;
+            last  = max+fAfter;
+          }
+
+        for (ptr=first;ptr<last;ptr++)
+          {
+            sumHi += *ptr;
+            if (*ptr >= fSaturationLimit)
+              satHi++;
+          }
+
+        ptr   = pixel.GetLoGainSamples();
+        max   = ptr+pixel.GetIdxMaxLoGainSample();
+        end   = ptr+fRawEvt->GetNumLoGainSamples();
+
+        if (max-4 < ptr)
+          {
+            first = ptr+4;
+            last  = ptr+4+fAfter+fBefore;
+          }
+        else if (max+fAfter > end)
+          {
+            first = end-fAfter-fBefore;
+            last  = end;
+          }
+        else
+          {
+            first = max-fBefore;
+            last  = max+fAfter;
+          }
+
+        for (ptr=first;ptr<last;ptr++)
+          {
+            sumLo += *ptr;
+            if (*ptr >= fSaturationLimit)
+              {
+                *fLog << err << dbginf 
+                      << "Warning: Saturation of Lo Gain reached in pixel: " 
+                      << pixid << " " << "   sum = " << sumLo << endl;
+                satLo++;
+              }
+          }
 
         //
@@ -171 +235 @@
         //
 
-        UInt_t max  = pixel.GetMaxHiGainSample();
-
-        if (max <= fSaturationLimit)  // take Hi Gain, no saturation
-          {
-            sum = (float)pixel.GetSumHiGainSamples() - pedes*fNumHiGainSamples;  
-            sumerr = pedrms*fSqrtHiGainSamples; 
-          }
-        else                          // Lo Gain
-          {
-
-            satHi++;
-
-            sum = (float)pixel.GetSumLoGainSamples() - pedes*fNumLoGainSamples ;
-            sum *= fConversionHiLo; 
-
-            sumerr = pedrms*fSqrtLoGainSamples; 
-
-            max = pixel.GetMaxLoGainSample();
-
-            if (max > fSaturationLimit)
-              {
-                *fLog << err << dbginf 
-                      << "Warning: Saturation of Lo Gain reached in pixel: " 
-                      << pixid << " " << "   sum = " << sum << endl;
-                satLo++;
-              }
-
-          }
-
-        pix.SetGainSaturation((satHi),(satHi),(satLo));
-        pix.SetExtractedSignal(sum,sumerr);
+        pix.SetExtractedSignal((float)sumHi - pedes*(float)fNumHiGainSamples,
+                               pedrms*fSqrtHiGainSamples,
+                               ((float)sumLo - pedes*(float)fNumLoGainSamples)*fConversionHiLo,
+                               pedrms*fSqrtLoGainSamples
+                               );
+
+        pix.SetGainSaturation((satHi),satHi,satLo);
 
       } /* while (pixel.Next()) */
 
+
+    fSignals->SetNumUsedFADCSlices(fNumHiGainSamples);
     fSignals->SetReadyToSave();