Changeset 5505

Timestamp:

11/29/04 18:46:05 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• mpedestal/MPedCalcFromLoGain.cc (modified) (3 diffs)
• mpedestal/MPedCalcPedRun.cc (modified) (3 diffs)
• msignal/MExtractTimeAndChargeDigitalFilter.cc (modified) (9 diffs)
• msignal/MExtractTimeAndChargeDigitalFilter.h (modified) (2 diffs)
• msignal/MExtractTimeAndChargeSlidingWindow.cc (modified) (1 diff)

trunk/MagicSoft/Mars/Changelog

@@ -34 +34 @@
      - fixed documentation of data members
      - Introduced setter and getter for name of weights file
+     - modified way to use the NoiseCalculation(). Now more compatible
+       with MPedCalcFromLoGain
 
    * msignal/MExtractFixedWindowPeakSearch.cc

trunk/MagicSoft/Mars/mpedestal/MPedCalcFromLoGain.cc

@@ -137 +137 @@
 
 #include "MExtractTimeAndCharge.h"
+#include "MExtractTimeAndChargeDigitalFilter.h"
 
 #include "MParList.h"
@@ -245 +246 @@
       SetExtractWindow(fExtractor->GetHiGainFirst(),(Int_t)fExtractor->GetNumHiGainSamples());
       lastavailable = losamples-1;
+      if (fExtractor->InheritsFrom("MExtractTimeAndChargeDigitalFilter"))
+        {
+          MExtractTimeAndChargeDigitalFilter* df = (MExtractTimeAndChargeDigitalFilter*)fExtractor;
+          if (!df->IsWeightsSet())
+            if (!df->ReadWeightsFile(df->GetNameWeightsFile()))
+              return kFALSE;
+        }
     }
 
@@ -438 +446 @@
   const Bool_t logainabflag = (pixel->HasABFlag() + pixel->GetNumHiGainSamples()) & 0x1;
   Float_t dummy;
-  fExtractor->FindTimeAndChargeLoGain(logain,sum,dummy,dummy,dummy,sat,ped,logainabflag);
+  fExtractor->FindTimeAndChargeHiGain(logain,logain,sum,dummy,dummy,dummy,sat,ped,logainabflag);
 }
 

trunk/MagicSoft/Mars/mpedestal/MPedCalcPedRun.cc

@@ -122 +122 @@
 
 #include "MExtractTimeAndCharge.h"
+#include "MExtractTimeAndChargeDigitalFilter.h"
 
 #include "MParList.h"
@@ -204 +205 @@
 
   if (fExtractor)
-    SetExtractWindow(fExtractor->GetHiGainFirst(),(Int_t)fExtractor->GetNumHiGainSamples());
-
+    {
+      SetExtractWindow(fExtractor->GetHiGainFirst(),(Int_t)fExtractor->GetNumHiGainSamples());
+      if (fExtractor->InheritsFrom("MExtractTimeAndChargeDigitalFilter"))
+        {
+          MExtractTimeAndChargeDigitalFilter* df = (MExtractTimeAndChargeDigitalFilter*)fExtractor;
+          if (!df->IsWeightsSet())
+            if (!df->ReadWeightsFile(df->GetNameWeightsFile()))
+              return kFALSE;
+        }
+    }
+  
   MExtractPedestal::ReInit(pList);
 
@@ -318 +328 @@
 }
 
+
 void MPedCalcPedRun::CalcExtractor( MRawEvtPixelIter *pixel, Float_t &sum, MPedestalPix &ped)
 {

trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeDigitalFilter.cc

@@ -91 +91 @@
 //
 MExtractTimeAndChargeDigitalFilter::MExtractTimeAndChargeDigitalFilter(const char *name, const char *title) 
-    : fWeightsSet(kFALSE)
+    : fWeightsSet(kFALSE), fRandomIter(0)
 {
     fName  = name  ? name  : "MExtractTimeAndChargeDigitalFilter";
@@ -188 +188 @@
 Int_t MExtractTimeAndChargeDigitalFilter::PreProcess(MParList *pList)
 {
-
   if (!fWeightsSet)
     if (!ReadWeightsFile(fNameWeightsFile))
@@ -223 +222 @@
 void MExtractTimeAndChargeDigitalFilter::InitArrays()
 {
-  fHiGainSignal.Set(fHiGainLast - fHiGainFirst + 1 + fHiLoLast);
-  fLoGainSignal.Set(fLoGainLast - fLoGainFirst + 1);
+  Int_t range = (Int_t)(fHiGainLast - fHiGainFirst + 1 + fHiLoLast);
+  fHiGainSignal.Set(range);
+  range = (Int_t)(fLoGainLast - fLoGainFirst + 1);
+  fLoGainSignal.Set(range);
 }
 
@@ -237 +238 @@
 
   Int_t range = fHiGainLast - fHiGainFirst + 1;
+  
   const Byte_t *end = ptr + range;
   Byte_t *p     = ptr;
@@ -305 +307 @@
     return;
 
+  //
+  // Slide with a window of size fWindowSizeHiGain over the sample 
+  // and multiply the entries with the corresponding weights
+  //
+  if (IsNoiseCalculation())
+    {
+      if (fRandomIter == 10)
+        fRandomIter = 0;
+      for (Int_t ids=0; ids < fWindowSizeHiGain; ids++)
+        {
+          const Int_t   idx = fBinningResolutionHiGain*ids + fRandomIter;
+          sum              += fAmpWeightsHiGain [idx]*fHiGainSignal[ids]; 
+        }
+      fRandomIter++;
+      return;
+    }
+  
   Float_t time_sum  = 0.;
   Float_t fmax      = 0.;
@@ -330 +349 @@
       }
 
-      //  if (IsNoiseCalculation()) return;
-
       if (sum>fmax)
       {
@@ -373 +390 @@
   {
     const Int_t   idx = fBinningResolutionHiGain*sample + fBinningResolutionHalfHiGain + t_iter;
-    const Int_t   ids = IsNoiseCalculation() ? sample : max_p + sample;
+    const Int_t   ids = max_p + sample;
     
-    const Float_t pex = ids < 0 ? 0. : 
-      ( ids > range ? 0. : fHiGainSignal[ids]);
+    const Float_t pex = ids < 0 ? 0. : ( ids > range ? 0. : fHiGainSignal[ids]);
     sum              += fAmpWeightsHiGain [idx]*pex; 
     time_sum         += fTimeWeightsHiGain[idx]*pex;
@@ -436 +452 @@
     }
   
+  //
+  // Slide with a window of size fWindowSizeHiGain over the sample 
+  // and multiply the entries with the corresponding weights
+  //
+  if (IsNoiseCalculation())
+    {
+      for (Int_t ids=0; ids < fWindowSizeLoGain; ids++)
+        {
+          const Int_t   idx = fBinningResolutionHiGain*ids + fRandomIter;
+          sum              += fAmpWeightsLoGain [idx]*fLoGainSignal[ids]; 
+        }
+      return;
+    }
+
   Float_t time_sum  = 0.;
   Float_t fmax      = 0.;
@@ -503 +533 @@
   {
     const Int_t   idx = fBinningResolutionLoGain*sample + fBinningResolutionHalfLoGain + t_iter;
-    const Int_t   ids = (IsNoiseCalculation()) ? sample : max_p + sample;;
-    const Float_t pex = ids < 0 ? 0. : 
-      ( ids > range ? 0. : fLoGainSignal[ids]);
+    const Int_t   ids = max_p + sample;;
+    const Float_t pex = ids < 0 ? 0. : ( ids > range ? 0. : fLoGainSignal[ids]);
     sum              += fAmpWeightsLoGain [idx]*pex; 
     time_sum         += fTimeWeightsLoGain[idx]*pex;

trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeDigitalFilter.h

@@ -50 +50 @@
   TString fNameWeightsFile;                     // Name of the weights file
   Bool_t  fWeightsSet;                          // Flag if weights have alreayd been set
-  
+  UInt_t  fRandomIter;                          // Counter used to randomize weights for noise calculation
+
   Bool_t ReInit(MParList *pList);
   Int_t PreProcess(MParList *pList);
@@ -87 +88 @@
 
   void Print(Option_t *o="") const;
-  
+
+  Bool_t IsWeightsSet() const  { return fWeightsSet;  }
+
   void FindTimeAndChargeHiGain(Byte_t *firstused, Byte_t *logain, Float_t &sum, Float_t &dsum,
                                Float_t &time, Float_t &dtime,

trunk/MagicSoft/Mars/msignal/MExtractTimeAndChargeSlidingWindow.cc

@@ -187 +187 @@
 void MExtractTimeAndChargeSlidingWindow::InitArrays()
 {
-  fHiGainSignal.Set(fHiGainLast - fHiGainFirst + 1 + fHiLoLast);
-  fLoGainSignal.Set(fLoGainLast - fLoGainFirst + 1);
+  Int_t range = (Int_t)(fHiGainLast - fHiGainFirst + 1 + fHiLoLast);
+  fHiGainSignal.Set(range);
+  range = (Int_t)(fLoGainLast - fLoGainFirst + 1);
+  fLoGainSignal.Set(range);
 }