Changeset 4368 for trunk/MagicSoft/Mars/msignal

Timestamp:

06/30/04 17:55:26 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/msignal

Files:

• MExtractor.cc (modified) (10 diffs)
• MExtractor.h (modified) (1 diff)

trunk/MagicSoft/Mars/msignal/MExtractor.cc

@@ -25 +25 @@
 //////////////////////////////////////////////////////////////////////////////
 //
-//   MExtractor
+// MExtractor
+// ==========
 //
 //   Base class for the signal extractors, used the functions 
@@ -38 +39 @@
 //   - fSqrtLoGainSamples
 //
+//   The signal extractor classes can be setup from an environmental
+//   setup file. For more information see ReadEnv and MEvtLoop::ReadEnv
+//
 // Input Containers:
 //   MRawEvtData
@@ -71 +75 @@
 
 const Byte_t MExtractor::fgSaturationLimit = 254;
+
 // --------------------------------------------------------------------------
 //
@@ -88 +93 @@
       fSaturationLimit(fgSaturationLimit)
 {
-
     fName  = name  ? name  : "MExtractor";
     fTitle = title ? title : "Base class for signal extractors";
@@ -99 +103 @@
 void MExtractor::SetRange(Byte_t hifirst, Byte_t hilast, Byte_t lofirst, Byte_t lolast)
 {
-  
-  fHiGainFirst = hifirst;
-  fHiGainLast  = hilast;
-
-  fLoGainFirst = lofirst;
-  fLoGainLast  = lolast;
-
-}
-
+    fHiGainFirst = hifirst;
+    fHiGainLast  = hilast;
+
+    fLoGainFirst = lofirst;
+    fLoGainLast  = lolast;
+}
 
 // --------------------------------------------------------------------------
@@ -123 +124 @@
 Int_t MExtractor::PreProcess(MParList *pList)
 {
-
     fRawEvt = (MRawEvtData*)pList->FindObject(AddSerialNumber("MRawEvtData"));
     if (!fRawEvt)
@@ -178 +178 @@
 Bool_t MExtractor::ReInit(MParList *pList)
 {
-  
-  Int_t lastdesired   = (Int_t)(fLoGainLast);
-  Int_t lastavailable = (Int_t)fRunHeader->GetNumSamplesLoGain()-1;
-
-  if (lastavailable < 0)
-    *fLog << warn << GetDescriptor() 
-          << ": WARNING: Number of available Low-Gain Slices is smaller than or equal zero!" << endl;
-  
+    Int_t lastdesired   = (Int_t)(fLoGainLast);
+    Int_t lastavailable = (Int_t)fRunHeader->GetNumSamplesLoGain()-1;
+
+    if (lastavailable < 0)
+        *fLog << warn << GetDescriptor() << " - WARNING: Number of available Low-Gain Slices is smaller than or equal zero!" << endl;
+
+    if (lastdesired > lastavailable)
+    {
+        const Int_t diff = lastdesired - lastavailable;
+
+        *fLog << endl;
+        *fLog << warn << GetDescriptor() << ": Selected Lo Gain FADC Window [";
+        *fLog << Form("%2i,%2i", (int)fLoGainFirst, lastdesired);
+        *fLog << "] ranges out of the available limits: [0," << Form("%2i", lastavailable) << "]" << endl;
+        *fLog << GetDescriptor() << ": Will reduce the upper edge to " << (int)(fLoGainLast - diff) << endl;
+        SetRange(fHiGainFirst, fHiGainLast, fLoGainFirst, fLoGainLast-diff);
+    }
+
+    lastdesired   = (Int_t)fHiGainLast;
+    lastavailable = (Int_t)fRunHeader->GetNumSamplesHiGain()-1;
+
+    if (lastavailable < 0)
+    {
+        *fLog << err << GetDescriptor();
+        *fLog << " - ERROR: Number of available High-Gain Slices is smaller than or equal zero!" << endl;
+        return kFALSE;
+    }
+
   if (lastdesired > lastavailable)
-    {
+  {
       const Int_t diff = lastdesired - lastavailable;
-      *fLog << endl; 
-      *fLog << warn << GetDescriptor()
-            << Form("%s%2i%s%2i%s%2i%s",": Selected Lo Gain FADC Window [",
-                    (int)fLoGainFirst,",",lastdesired,
-                    "] ranges out of the available limits: [0,",lastavailable,"].") << endl;
-      *fLog << GetDescriptor() << ": Will reduce the upper edge to " << (int)(fLoGainLast - diff) << endl;
-      SetRange(fHiGainFirst, fHiGainLast, fLoGainFirst, fLoGainLast-diff);
-    }
-
-  lastdesired   = (Int_t)fHiGainLast;
-  lastavailable = (Int_t)fRunHeader->GetNumSamplesHiGain()-1;
-  
-  if (lastavailable < 0)
-    {
-      *fLog << err << GetDescriptor() 
-            << ": ERROR: Number of available High-Gain Slices is smaller than or equal zero!" << endl;
-      return kFALSE;
-    }
-
-  if (lastdesired > lastavailable)
-    {
-      const Int_t diff = lastdesired - lastavailable;
+
       *fLog << endl;
-      *fLog << warn << GetDescriptor()
-            << Form("%s%2i%s%2i%s%2i%s",": Selected Hi Gain FADC Window [",
-                    (int)fHiGainFirst,",",lastdesired,
-                    "] ranges out of the available limits: [0,",lastavailable,"].") << endl;
-      *fLog << warn << GetDescriptor() 
-            << Form("%s%2i%s",": Will use ",diff," samples from the Low-Gain for the High-Gain extraction")
-            << endl;
+      *fLog << warn << GetDescriptor() << ": Selected Hi Gain FADC Window [";
+      *fLog << Form("%2i,%2i%s%2i%s", (int)fHiGainFirst,",",lastdesired);
+      *fLog << "] ranges out of the available limits: [0," << Form("%2i", lastavailable) << "]" << endl;
+      *fLog << warn << GetDescriptor() << ": Will use ";
+      *fLog << Form("%2i", diff) << " samples from the Low-Gain for the High-Gain extraction";
+      *fLog << endl;
+
       fHiGainLast -= diff;
       fHiLoLast    = diff;
-    }
+  }
 
   return kTRUE;
-}
-
-
-
-void MExtractor::FindSignalHiGain(Byte_t *firstused, Byte_t *logain, Float_t &sum, Byte_t &sat) const
-{
-  return;
-}
-
-void MExtractor::FindSignalLoGain(Byte_t *firstused, Float_t &sum, Byte_t &sat) const
-{
-  return;
 }
 
@@ -245 +232 @@
 Int_t MExtractor::Process()
 {
-
   MRawEvtPixelIter pixel(fRawEvt);
   fSignals->Clear();
@@ -290 +276 @@
 void MExtractor::StreamPrimitive(ofstream &out) const
 {
-
   out << "   " << ClassName() << " " << GetUniqueName() << "(\"";
   out << "\"" << fName << "\", \"" << fTitle << "\");" << endl;
   
   if (fSaturationLimit!=fgSaturationLimit)
-    {
+  {
       out << "   " << GetUniqueName() << ".SetSaturationLimit(";
       out << (int)fSaturationLimit << ");" << endl;
-    }
+  }
   
   out << "   " << GetUniqueName() << ".SetRange(";
@@ -308 +293 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Read the setup from a TEnv, eg:
+//   MJPedestal.MExtractor.HiGainFirst: 5
+//   MJPedestal.MExtractor.LoGainFirst: 5
+//   MJPedestal.MExtractor.HiGainLast:  10
+//   MJPedestal.MExtractor.LoGainLast:  10
+//   MJPedestal.MExtractor.SaturationLimit: 88
+//
+Bool_t MExtractor::ReadEnv(const TEnv &env, TString prefix, Bool_t print)
+{
+    Byte_t hf = fHiGainFirst;
+    Byte_t lf = fLoGainFirst;
+    Byte_t hl = fHiGainLast;
+    Byte_t ll = fLoGainLast;
+
+    if (IsEnvDefined(env, prefix, "HiGainFirst", print))
+        hf = GetEnvValue(env, prefix, "HiGainFirst", hf);
+    if (IsEnvDefined(env, prefix, "LoGainFirst", print))
+        lf = GetEnvValue(env, prefix, "LoGainFirst", lf);
+
+    if (IsEnvDefined(env, prefix, "HiGainLast", print))
+        hl = GetEnvValue(env, prefix, "HiGainLast", hl);
+    if (IsEnvDefined(env, prefix, "LoGainLast", print))
+        ll = GetEnvValue(env, prefix, "LoGainLast", ll);
+
+    SetRange(hf, hl, lf, ll);
+
+    if (IsEnvDefined(env, prefix, "SaturationLimit", print))
+        SetSaturationLimit(GetEnvValue(env, prefix, "SaturationLimit", fSaturationLimit));
+
+    return kTRUE;
+}

trunk/MagicSoft/Mars/msignal/MExtractor.h

@@ -47 +47 @@
     Byte_t  fSaturationLimit;
 
-    virtual void FindSignalHiGain(Byte_t *firstused, Byte_t *lowgain, Float_t &sum, Byte_t &sat) const;
-    virtual void FindSignalLoGain(Byte_t *firstused, Float_t &sum, Byte_t &sat) const;
+    virtual void FindSignalHiGain(Byte_t *firstused, Byte_t *lowgain, Float_t &sum, Byte_t &sat) const { }
+    virtual void FindSignalLoGain(Byte_t *firstused, Float_t &sum, Byte_t &sat) const { }
 
-    virtual Int_t  PreProcess(MParList *pList);
-    virtual Bool_t ReInit(MParList *pList);
-    virtual Int_t  Process();
-    virtual void   StreamPrimitive(ofstream &out) const;
+    Int_t  PreProcess(MParList *pList);
+    Bool_t ReInit(MParList *pList);
+    Int_t  Process();
+    void   StreamPrimitive(ofstream &out) const;
+    Bool_t ReadEnv(const TEnv &env, TString prefix, Bool_t print);
 
 public:
-
     MExtractor(const char *name=NULL, const char *title=NULL);
 
-    Byte_t  GetHiGainFirst()      const      { return fHiGainFirst;      }
-    Byte_t  GetHiGainLast ()      const      { return fHiGainLast ;      }
-    Byte_t  GetLoGainFirst()      const      { return fLoGainFirst;      }
-    Byte_t  GetLoGainLast ()      const      { return fLoGainLast ;      }
-    Float_t GetNumHiGainSamples() const      { return fNumHiGainSamples; }
-    Float_t GetNumLoGainSamples() const      { return fNumLoGainSamples; }
+    Byte_t  GetHiGainFirst()      const { return fHiGainFirst;      }
+    Byte_t  GetHiGainLast ()      const { return fHiGainLast ;      }
+    Byte_t  GetLoGainFirst()      const { return fLoGainFirst;      }
+    Byte_t  GetLoGainLast ()      const { return fLoGainLast ;      }
+    Float_t GetNumHiGainSamples() const { return fNumHiGainSamples; }
+    Float_t GetNumLoGainSamples() const { return fNumLoGainSamples; }
     
     virtual void SetRange  (Byte_t hifirst=0, Byte_t hilast=0, Byte_t lofirst=0, Byte_t lolast=0);