Changeset 8585 for trunk/MagicSoft/Mars/msignal

Timestamp:

06/19/07 12:17:36 (17 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/msignal/MExtractor.cc (modified) (4 diffs)

trunk/MagicSoft/Mars/msignal/MExtractor.cc

@@ -136 +136 @@
     : fResolutionPerPheHiGain(0), fResolutionPerPheLoGain(0),
     fSignals(NULL), fRawEvt(NULL), fRunHeader(NULL), fSignal(NULL),
-      /*fHiLoLast(0),*/ fNumHiGainSamples(0), fNumLoGainSamples(0)
+    fNumHiGainSamples(0), fNumLoGainSamples(0)
 {
     fName  = name  ? name  : "MExtractor";
@@ -251 +251 @@
         fLoGainFirst=0;
         fLoGainLast =0;
-        if (fSaturationLimit>fRunHeader->GetMax())
-        {
-            *fLog << " and saturation limit to " << fRunHeader->GetMax();
-            fSaturationLimit=fRunHeader->GetMax();
-        }
         *fLog << "." << endl;
     }
 
-/*
-  const Int_t logainsamples = fRunHeader->GetNumSamplesLoGain();
-
-  Int_t lastdesired;
-  Int_t lastavailable;
-
-  if (logainsamples)
-    {
-
-      lastdesired   = (Int_t)(fLoGainLast);
-      lastavailable = logainsamples-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);
-      }
-    }
-  else
-    SetRange(fHiGainFirst, fHiGainLast, 0,0);      
-  
-  const Int_t higainsamples = fRunHeader->GetNumSamplesHiGain();
-
-  if (higainsamples <= 0)
-    {
-      *fLog << err << GetDescriptor();
-      *fLog << " - ERROR: Number of available High-Gain Slices is smaller than or equal zero!" << endl;
-      return kFALSE;
-    }
-
-  lastdesired   = (Int_t)fHiGainLast;
-  lastavailable = higainsamples-1;
-  
-  if (lastdesired > lastavailable)
-    {
-      const Int_t diff = lastdesired - lastavailable;
-      
-      *fLog << endl;
-      *fLog << inf << GetDescriptor() << ": Selected Hi Gain FADC Window [";
-      *fLog << Form("%2i,%2i", (int)fHiGainFirst,lastdesired);
-      *fLog << "] ranges out of the available limits: [0," << Form("%2i", lastavailable) << "]" << endl;
-      *fLog << inf << 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;
 }
@@ -329 +265 @@
 {
     return kERROR;
-/*
-  MRawEvtPixelIter pixel(fRawEvt);
-
-  while (pixel.Next())
-    {
-      Float_t sumhi = 0.;
-      Byte_t sathi = 0;
-
-      FindSignalHiGain(pixel.GetHiGainSamples()+fHiGainFirst, pixel.GetLoGainSamples(), sumhi, sathi);
-
-      Float_t sumlo  = 0.;
-      Byte_t  satlo  = 0;
-
-      if (pixel.HasLoGain())
-        FindSignalLoGain(pixel.GetLoGainSamples()+fLoGainFirst, sumlo, satlo);
-
-      const Int_t pixid = pixel.GetPixelId();
-      
-      const MPedestalPix  &ped = (*fPedestals)[pixid]; 
-      MExtractedSignalPix &pix = (*fSignals)[pixid];
-      
-      const Float_t pedes  = ped.GetPedestal();
-      const Float_t pedrms = ped.GetPedestalRms();
-
-      pix.SetExtractedSignal(sumhi - pedes*fNumHiGainSamples, pedrms*fSqrtHiGainSamples,
-                             sumlo - pedes*fNumLoGainSamples, pedrms*fSqrtLoGainSamples);
-      
-      pix.SetGainSaturation(sathi, satlo);
-      
-    }
-
-    fSignals->SetReadyToSave();
-
-    return kTRUE;
-    */
 }
 
@@ -461 +362 @@
 
     *fLog << " Hi Gain Range:      " << Form("%2d %2d", fHiGainFirst, fHiGainLast) << endl;
+    *fLog << " Saturation Lim:     " << Form("%3d", fSaturationLimit) << endl;
     if (HasLoGain())
     {
         *fLog << " Lo Gain Range:      " << Form("%2d %2d", fLoGainFirst, fLoGainLast) << endl;
-        *fLog << " Saturation Lim:     " << Form("%3d", fSaturationLimit) << endl;
         *fLog << " Num Samples Hi/Lo:  " << Form("%2.1f %2.1f", fNumHiGainSamples, fNumLoGainSamples) << endl;
     }