Changeset 7844 for trunk/MagicSoft/Mars/msignal

Timestamp:

08/04/06 10:16:57 (18 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/msignal

Files:

• MExtractTimeAndCharge.cc (modified) (4 diffs)
• MExtractTimeAndChargeDigitalFilter.cc (modified) (6 diffs)
• MExtractTimeAndChargeSpline.cc (modified) (3 diffs)

trunk/MagicSoft/Mars/msignal/MExtractTimeAndCharge.cc

@@ -179 +179 @@
 
   while (pixel.Next())
-    {
+  {
+      // COPY HERE PRODUCING ARRAY WITH SAMPLES
+
+      /*
+       const MPedestalPix  &ped = (*fPedestals)[pixidx];
+
+       const Float_t pedes  = ped.GetPedestal();
+       const Float_t ABoffs = ped.GetPedestalABoffset();
+
+       const Float_t pedmean[2] = { pedes + ABoffs, pedes - ABoffs };
+
+       const Int_t num = pixel.GetNumHiGainSamples()+pixel.GetNumLoGainSamples();
+
+       MArrayF sample(num);
+
+       const Int_t abflag = pixel.HasABFlag() ? 1 : 0;
+       const Int_t ids0 = fHiGainFirst + abflag;
+
+       Int_t ids  = ids0;
+
+       Float_t  null   = 0;      // Starting value for maxpos
+       Float_t *maxpos = &null;  // Position of maximum
+       Float_t *sat1   = 0;      // First saturating slice
+       Float_t *sat2   = 0;      // Last saturating slice
+
+       const Float_t *beg = sample.GetArray();
+       const Float_t *end = beg + num;
+
+       Float_t *ptr = beg;
+       while (ptr<end)
+       {
+          *sample++ = (Float_t)*ptr - pedmean[ids++&0x1];
+
+          // This extraction must be done independant for lo- and hi-gain
+          // if (*ptr > *maxpos)
+          //     maxpos = ptr;
+          //
+          // if (*ptr >= fSaturationLimit)
+          // {
+          //    sat2 = ptr;
+          //    if (!sat1)
+          //       sat1 = ptr;
+          // }
+          //
+          // ptr++;
+        }
+       */
+
       //
       // Find signal in hi- and lo-gain
@@ -210 +257 @@
       // Adapt the low-gain extraction range from the obtained high-gain time
       //
-      if (pixel.HasLoGain() && (fMaxBinContent > fLoGainSwitch) )
+
+      // IN THIS CASE THE PIXEL SHOULD BE MARKED BAD!!!!
+      // MEANS: Hi gain has saturated, but it is too early to extract
+      // the lo-gain properly
+      // THIS produces pulse positions ~= -1
+      // The signal might be handled in MCalibrateData, but hwat's about
+      // the arrival times in MCalibrateRelTime
+      if (sathi && fMaxBinContent<=fLoGainSwitch)
+          deltasumlo=deltasumhi=deltatimelo=deltatimehi=-1;
+
+      // FIXME: What to do with the pixel if it saturates too early???
+      if (pixel.HasLoGain() && (fMaxBinContent > fLoGainSwitch /*|| sathi>0*/) )
       {
           deltasumlo  = 0; // make logain of MExtractedSignalPix valid
@@ -217 +275 @@
           fLoGainFirstSave = fLoGainFirst;
 
-          const Float_t pos = sathi==0 ? timehi : sathi;
+          // sathi is the number (not index!) of the first saturating slice
+          // 0 indicates that no saturation was found
+          // FIMXME: Is 3 an accurate value?
+
+          const Float_t pos = sathi==0 ? timehi : (int)(sathi)-3;
 
           if (pos>-fLoGainStartShift)
@@ -225 +287 @@
               fLoGainFirst = fLoGainFirstSave;
 
-          if ( fLoGainFirst <= fLoGainLast-fWindowSizeLoGain)
+          if (fLoGainFirst <= fLoGainLast-fWindowSizeLoGain)
           {
               const Bool_t logainabflag = (higainabflag + pixel.GetNumHiGainSamples()) & 0x1;

trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeDigitalFilter.cc

@@ -24 +24 @@
 !
 \* ======================================================================== */
+
 //////////////////////////////////////////////////////////////////////////////
 //
@@ -88 +89 @@
 const Int_t  MExtractTimeAndChargeDigitalFilter::fgSignalStartBinHiGain    =  4;
 const Int_t  MExtractTimeAndChargeDigitalFilter::fgSignalStartBinLoGain    =  4;
-const Float_t MExtractTimeAndChargeDigitalFilter::fgOffsetLoGain           =  1.05;
+const Float_t MExtractTimeAndChargeDigitalFilter::fgOffsetLoGain           =  0.95;
 const Float_t MExtractTimeAndChargeDigitalFilter::fgLoGainStartShift       = -1.8;
 
@@ -329 +330 @@
   const Byte_t *end = ptr + range;
   Byte_t *p     = ptr;
-  Byte_t maxpos = 0;
 
   //
@@ -352 +352 @@
       if (*p > fMaxBinContent)
         {
-          maxpos = p-ptr;
+          Byte_t maxpos = p-ptr;
+
+          // range-fWindowSizeHiGain+1 == fHiLoLast isn't it?
           if (maxpos > 1 && maxpos < (range - fWindowSizeHiGain + 1))
             fMaxBinContent = *p;
@@ -359 +361 @@
       if (*p++ >= fSaturationLimit)
         if (!sat)
-          sat = ids-3;
+          sat = ids-fHiGainFirst;
     }
 
@@ -391 +393 @@
           if (*logain++ >= fSaturationLimit)
             if (!sat)
-              sat = ids-3;
+              sat = ids-fHiGainFirst;
         }
     }

trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeSpline.cc

@@ -160 +160 @@
 const Float_t MExtractTimeAndChargeSpline::fgLoGainStretch    = 1.5;
 const Float_t MExtractTimeAndChargeSpline::fgOffsetLoGain     = 1.39;  // 5 ns
-const Float_t MExtractTimeAndChargeSpline::fgLoGainStartShift = -1.8;  
+const Float_t MExtractTimeAndChargeSpline::fgLoGainStartShift = -2.4;
 
 // --------------------------------------------------------------------------
@@ -404 +404 @@
       if (*p++ >= fSaturationLimit)
         if (!sat)
-          sat = ids-3;
+          sat = ids-fHiGainFirst;
       
     }
@@ -426 +426 @@
           if (*logain++ >= fSaturationLimit)
             if (!sat)
-              sat = ids-3;
+              sat = ids-fHiGainFirst;
 
           range++;