Changeset 8151 for trunk/MagicSoft/Mars/mpedestal/MPedCalcFromLoGain.cc

Timestamp:

10/24/06 08:58:13 (18 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mpedestal/MPedCalcFromLoGain.cc (modified) (16 diffs)

trunk/MagicSoft/Mars/mpedestal/MPedCalcFromLoGain.cc

@@ -1 +1 @@
 /* ======================================================================== *\
+! $Name: not supported by cvs2svn $:$Id: MPedCalcFromLoGain.cc,v 1.34 2006-10-24 07:58:13 tbretz Exp $
+! --------------------------------------------------------------------------
 !
 ! *
@@ -24 +26 @@
 !   Author(s): Nepomuk Otte 10/2004 <mailto:otte@mppmu.mpg.de>
 !
-!   Copyright: MAGIC Software Development, 2000-2004
+!   Copyright: MAGIC Software Development, 2000-2006
 !
 !
@@ -142 +144 @@
 #include "MRawRunHeader.h"  
 #include "MRawEvtPixelIter.h"
-#include "MRawEvtData.h"
 
 #include "MPedestalPix.h"
@@ -151 +152 @@
 
 #include "MExtractPedestal.h"
-#include "MExtractTimeAndCharge.h"
+#include "MPedestalSubtractedEvt.h"
 
 ClassImp(MPedCalcFromLoGain);
@@ -157 +158 @@
 using namespace std;
 
-const UShort_t MPedCalcFromLoGain::fgCheckWinFirst   =  0;
-const UShort_t MPedCalcFromLoGain::fgCheckWinLast    = 29;
 const UShort_t MPedCalcFromLoGain::fgExtractWinFirst = 17;
 const UShort_t MPedCalcFromLoGain::fgExtractWinSize  =  6;
-const UShort_t MPedCalcFromLoGain::fgMaxSignalVar    = 40;
+const UInt_t   MPedCalcFromLoGain::fgNumDump         = 500;
 
 // --------------------------------------------------------------------------
@@ -167 +166 @@
 // Default constructor: 
 //
-// Sets:
-// - all pointers to NULL
-//
 // Calls: 
-// - AddToBranchList("fHiGainPixId");
-// - AddToBranchList("fHiGainFadcSamples");
-// - SetCheckRange(fgCheckWinFirst, fgCheckWinLast, fgExtractWinFirst, fgExtractWinSize)
+// - SetExtractWindow(fgExtractWinFirst, fgExtractWinSize)
 //
 MPedCalcFromLoGain::MPedCalcFromLoGain(const char *name, const char *title)
@@ -180 +174 @@
     fTitle = title ? title : "Task to calculate pedestals from lo-gains";
 
-    SetCheckRange(fgCheckWinFirst, fgCheckWinLast);
     SetExtractWindow(fgExtractWinFirst, fgExtractWinSize);
-
-    SetMaxSignalVar(fgMaxSignalVar);
-}
-
+    SetPedestalUpdate(kTRUE);
+    SetNumDump();
+}
+
+// --------------------------------------------------------------------------
+//
+// Call MExtractPedestl::ResetArrays aand reset fNumAventsUsed and
+// fTotalCounter
+//
 void MPedCalcFromLoGain::ResetArrays()
 {
@@ -192 +190 @@
     fNumEventsUsed.Reset();
     fTotalCounter.Reset();
-}
-
-// --------------------------------------------------------------------------
-//
-// SetCheckRange: 
-// 
-// Exits, if the first argument is smaller than 0
-// Exits, if the the last argument is smaller than the first
-//
-Bool_t MPedCalcFromLoGain::SetCheckRange(UShort_t chfirst, UShort_t chlast)
-{
-
-  Bool_t rc = kTRUE;
-  
-  if (chlast<=chfirst)
-    {
-      *fLog << warn << GetDescriptor();
-      *fLog << " - WARNING: Last slice in SetCheckRange smaller than first slice... set to first+2" << endl;
-      chlast = chfirst+1;
-      rc = kFALSE;
-    }
-
-  fCheckWinFirst = chfirst;
-  fCheckWinLast  = chlast;
-
-  return rc;
 }
 
@@ -231 +203 @@
 Bool_t MPedCalcFromLoGain::ReInit(MParList *pList)
 {
-
-  const UShort_t hisamples = fRunHeader->GetNumSamplesHiGain();
-  const UShort_t losamples = fRunHeader->GetNumSamplesLoGain();
-  
-  fSlices.Set(hisamples+losamples);
-  
-  UShort_t lastavailable   = hisamples+losamples-1;
-
-  if (fExtractor)
-    fExtractWinLast = fExtractWinFirst + fExtractor->GetWindowSizeHiGain() - 1;
-  
-  // If the size is not yet set, set the size
-  if (fSumx.GetSize()==0)
-    {
-      const Int_t npixels  = fPedestalsOut->GetSize();
-      fNumEventsUsed.Set(npixels);
-      fTotalCounter.Set(npixels);
-    }
-  
-  if (fCheckWinLast > lastavailable) //changed to override check
-    {
-      *fLog << warn << GetDescriptor();
-      *fLog << " - WARNING: Last Check Window slice out of range...adjusting to last available slice ";
-      *fLog << lastavailable << endl;
-      
-      fCheckWinLast = lastavailable;
-    }
-  
-  if (fExtractWinLast > lastavailable) 
-    {
-      if (fExtractor)
-        {
-          *fLog << err << GetDescriptor();
-          *fLog << " - ERROR: Selected Last Extraction Window: " 
-                << fExtractWinFirst + fExtractor->GetWindowSizeHiGain()-1 
-                << "ranges out of range: " << lastavailable-1 << endl;
-          return kFALSE;
-        }
-      else
-        {
-          const UShort_t diff = fExtractWinLast - lastavailable;
-          *fLog << warn << GetDescriptor();
-          *fLog << " - WARNING: Selected Extract Window ranges out of range...adjusting to last available slice ";
-          *fLog << lastavailable << endl;
-          
-          fExtractWinLast -= diff;
-          fExtractWinSize -= diff;
-        }
-    }
-  
-  return MExtractPedestal::ReInit(pList);
+    // If the size is not yet set, set the size
+    if (fSumx.GetSize()==0)
+    {
+        const Int_t npixels = fPedestalsOut->GetSize();
+        fNumEventsUsed.Set(npixels);
+        fTotalCounter.Set(npixels);
+    }
+
+    if (!MExtractPedestal::ReInit(pList))
+        return kFALSE;
+
+    return kTRUE;
 }
 
@@ -292 +225 @@
 Int_t MPedCalcFromLoGain::Calc()
 {
-    // This is the workaround to put hi- and lo-gains together
-    const Int_t nhigain = fRunHeader->GetNumSamplesHiGain();
-    const Int_t nlogain = fRunHeader->GetNumSamplesLoGain();
-
-    Byte_t *slices = fSlices.GetArray();
-
     // Real Process
     MRawEvtPixelIter pixel(fRawEvt);
@@ -303 +230 @@
     while (pixel.Next())
     {
-        // This is the fast workaround to put hi- and lo-gains together
-        memcpy(slices,         pixel.GetHiGainSamples(), nhigain);
-        memcpy(slices+nhigain, pixel.GetLoGainSamples(), nlogain);
-
-        // Start 'real' work
         const UInt_t idx = pixel.GetPixelId();
+
+        if (!CheckVariation(idx))
+            continue;
+
+        //extract pedestal
+        UInt_t ab[2];
+        const Float_t sum = fExtractor ?
+            CalcExtractor(pixel, pixel.GetNumHiGainSamples()) :
+            CalcSums(pixel, pixel.GetNumHiGainSamples(), ab[0], ab[1]);
 
         const UInt_t aidx   = (*fGeom)[idx].GetAidx();
         const UInt_t sector = (*fGeom)[idx].GetSector();
-
-        UShort_t max = 0;
-        UShort_t min = (UShort_t)-1;
-
-        // Find the maximum and minimum signal per slice in the high gain window
-        for (Byte_t *slice=slices+fCheckWinFirst; slice<=slices+fCheckWinLast; slice++)
-        {
-            if (*slice > max)
-                max = *slice;
-            if (*slice < min)
-                min = *slice;
-        }
-      
-        // If the maximum in the high gain window is smaller than
-        if (max-min>=fMaxSignalVar || max>=250)
-            continue;
-
-        Float_t sum  = 0.;
-
-        //extract pedestal
-        if (fExtractor)
-            CalcExtractor(pixel, sum, (*fPedestalsIn)[idx]);
-        else
-        {
-            UInt_t  sumi = 0;
-            for(Byte_t *slice=slices+fExtractWinFirst; slice<=slices+fExtractWinLast; slice++)
-                sumi += *slice;
-            sum = (Float_t)sumi;
-        }
 
         const Float_t sqrsum = sum*sum;
@@ -352 +254 @@
 
         if (fIntermediateStorage)
-          (*fPedestalsInter)[idx].Set(sum,0.,0.,fNumEventsUsed[idx]);
+            (*fPedestalsInter)[idx].Set(sum, 0, 0, fNumEventsUsed[idx]);
 
         fNumEventsUsed[idx]   ++;
@@ -358 +260 @@
         fSectorFilled [sector]++;
 
-        if (!fExtractor)
+        if (!fExtractor && pixel.IsABFlagValid())
         {
-            //
-            // Calculate the amplitude of the 150MHz "AB" noise
-            //
-            const UShort_t abFlag = (fExtractWinFirst + pixel.HasABFlag()) & 0x1;
-
-            for (Byte_t *slice=slices+fExtractWinFirst; slice<=slices+fExtractWinLast; slice+=2)
-            {
-                const UShort_t ab0 = *(slice + abFlag);
-                const UShort_t ab1 = *(slice - abFlag + 1);
-
-                fSumAB0[idx]        += ab0;
-                fSumAB1[idx]        += ab1;
-                fAreaSumAB0[aidx]   += ab0;
-                fAreaSumAB1[aidx]   += ab1;
-                fSectorSumAB0[aidx] += ab0;
-                fSectorSumAB1[aidx] += ab1;
-            }
+            fSumAB0[idx]        += ab[0];
+            fSumAB1[idx]        += ab[1];
+            fAreaSumAB0[aidx]   += ab[0];
+            fAreaSumAB1[aidx]   += ab[1];
+            fSectorSumAB0[aidx] += ab[0];
+            fSectorSumAB1[aidx] += ab[1];
         }
 
@@ -425 +316 @@
             CalcAreaResults(fAreaFilled[aidx], fAreaValid[aidx], aidx);
 
-}
-
-void MPedCalcFromLoGain::CalcExtractor(const MRawEvtPixelIter &pixel, Float_t &sum, MPedestalPix &ped)
-{
-    Byte_t  sat  = 0;
-    Byte_t *logain = pixel.GetLoGainSamples() + fExtractWinFirst;
-
-    const Bool_t logainabflag = (pixel.HasABFlag() + pixel.GetNumHiGainSamples()) & 0x1;
-
-    Float_t dummy;
-    fExtractor->FindTimeAndChargeHiGain(logain,logain,sum,dummy,dummy,dummy,sat,ped,logainabflag);
 }
 
@@ -469 +349 @@
 }
 
-
 // --------------------------------------------------------------------------
 //
 //  The following resources are available:
-//    FirstCheckWindowSlice:  0
-//    LastCheckWindowSlice:  29
-//    MaxSignalVar:          40
 //
 Int_t MPedCalcFromLoGain::ReadEnv(const TEnv &env, TString prefix, Bool_t print)
@@ -481 +357 @@
     Bool_t rc=kFALSE;
 
-    // Find resources for CheckWindow
-    Int_t fs = fCheckWinFirst;
-    Int_t ls = fCheckWinLast;
-    if (IsEnvDefined(env, prefix, "CheckWinFirst", print))
-    {
-        fs = GetEnvValue(env, prefix, "CheckWinFirst", fs);
+    // find resource for pedestal update
+    if (IsEnvDefined(env, prefix, "PedestalUpdate", print))
+    {
+        SetPedestalUpdate(GetEnvValue(env, prefix, "PedestalUpdate", fPedestalUpdate));
         rc = kTRUE;
     }
-    if (IsEnvDefined(env, prefix, "CheckWinLast", print))
-    {
-        ls = GetEnvValue(env, prefix, "CheckWinLast", ls);
+
+    // find resource for numdump
+    if (IsEnvDefined(env, prefix, "NumDump", print))
+    {
+        const Int_t num = GetEnvValue(env, prefix, "NumDump", -1);
+        if (num<=0)
+        {
+            *fLog << err << GetDescriptor() << ": ERROR - NumDump invalid!" << endl;
+            return kERROR;
+        }
+
+        SetNumDump(num);
         rc = kTRUE;
     }
 
-    SetCheckRange(fs,ls);
-
-    // find resource for maximum signal variation
-    if (IsEnvDefined(env, prefix, "MaxSignalVar", print))
-    {
-        SetMaxSignalVar(GetEnvValue(env, prefix, "MaxSignalVar", fMaxSignalVar));
+    // find resource for numeventsdump
+    if (IsEnvDefined(env, prefix, "NumEventsDump", print))
+    {
+        SetNumEventsDump(GetEnvValue(env, prefix, "NumEventsDump", (Int_t)fNumEventsDump));
         rc = kTRUE;
     }
 
-    return MExtractPedestal::ReadEnv(env,prefix,print) ? kTRUE : rc;
+    // find resource for numeventsdump
+    if (IsEnvDefined(env, prefix, "NumAreasDump", print))
+    {
+        SetNumAreasDump(GetEnvValue(env, prefix, "NumAreasDump", (Int_t)fNumAreasDump));
+        rc = kTRUE;
+    }
+
+    // find resource for numeventsdump
+    if (IsEnvDefined(env, prefix, "NumSectorsDump", print))
+    {
+        SetNumSectorsDump(GetEnvValue(env, prefix, "NumSectorsDump", (Int_t)fNumSectorsDump));
+        rc = kTRUE;
+    }
+
+    return rc;
 }
 
 void MPedCalcFromLoGain::Print(Option_t *o) const
 {
-
     MExtractPedestal::Print(o);
 
-    const Int_t last = fExtractor 
-      ? fExtractWinFirst + fExtractor->GetWindowSizeHiGain() -1
-      : fExtractWinLast;
-
-    *fLog << "ExtractWindow from slice " << fExtractWinFirst << " to " << last << " incl." << endl;
-    *fLog << "Max.allowed signal variation: " << fMaxSignalVar << endl;
-    *fLog << "CheckWindow from slice " << fCheckWinFirst   << " to " << fCheckWinLast << " incl." << endl;
-}
+    *fLog << "Pedestal Update is            " << (fPedestalUpdate?"on":"off") << endl;
+    if (fPedestalUpdate)
+    {
+        *fLog << "Num evts for pedestal   calc: " << fNumEventsDump << endl;
+        *fLog << "Num evts for avg.areas  calc: " << fNumAreasDump << endl;
+        *fLog << "Num evts for avg.sector calc: " << fNumSectorsDump << endl;
+    }
+}