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

Timestamp:

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

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/msignal/MExtractTimeAndCharge.cc (modified) (4 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;