Changeset 3112 for trunk/MagicSoft/Mars/mfilter

Timestamp:

02/12/04 11:34:32 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mfilter

Files:

• MFCosmics.cc (modified) (7 diffs)
• MFCosmics.h (modified) (2 diffs)

trunk/MagicSoft/Mars/mfilter/MFCosmics.cc

@@ -48 +48 @@
 //
 //  Output Containers:
+//   -/-
 //
 //////////////////////////////////////////////////////////////////////////////
@@ -77 +78 @@
       fRawEvt(NULL), fMaxEmptyPixels(230)
 {
-
     fName  = name  ? name  : "MFCosmics";
     fTitle = title ? title : "Filter to reject cosmics";
@@ -91 +91 @@
 Int_t MFCosmics::PreProcess(MParList *pList)
 {
-
     fRawEvt = (MRawEvtData*)pList->FindObject("MRawEvtData");
     if (!fRawEvt)
     {
-      *fLog << err << dbginf << "MRawEvtData not found... aborting." << endl;
+      *fLog << err << "MRawEvtData not found... aborting." << endl;
       return kFALSE;
     }
@@ -101 +100 @@
     fPedestals = (MPedestalCam*)pList->FindObject("MPedestalCam");
     if (!fPedestals)
-      {
-        *fLog << err << dbginf << "Cannot find MPedestalCam ... aborting" << endl;
+    {
+        *fLog << err << "MPedestalCam not found... aborting." << endl;
         return kFALSE;
-      }
+    }
 
     fSignals = (MExtractedSignalCam*)pList->FindObject("MExtractedSignalCam");
     if (!fSignals)
-      {
-        *fLog << err << dbginf << "Cannot find MExtractedSignalCam ... aborting" << endl;
+    {
+        *fLog << err << "MExtractedSignalCam not found... aborting." << endl;
         return kFALSE;
-      }
+    }
 
     memset(fCut, 0, sizeof(fCut));
@@ -124 +123 @@
 //  - MExtractedSignalCam
 //
-Bool_t MFCosmics::ReInit(MParList *pList )
-{
- 
+Bool_t MFCosmics::ReInit(MParList *pList)
+{
     fSqrtHiGainSamples = TMath::Sqrt((Float_t) fSignals->GetNumUsedHiGainFADCSlices());
 
@@ -140 +138 @@
 Int_t MFCosmics::Process()
 {
-
-  fResult = CosmicsRejection();
-  
-  fCut[fResult ? 0 : 1]++;
-  return kTRUE;
+    fResult = CosmicsRejection();
+
+    fCut[fResult ? 0 : 1]++;
+    return kTRUE;
 }
 
@@ -158 +155 @@
 // the outer pixels have some defect).
 //
-Bool_t MFCosmics::CosmicsRejection()
-{
-  
-  MRawEvtPixelIter pixel(fRawEvt);
-  
-  Int_t cosmicpix = 0;
-      
-  //
-  // Create a first loop to sort out the cosmics ...
-  // 
-  while (pixel.Next())
-    {
-      
-      const UInt_t idx = pixel.GetPixelId();
-      
-      MExtractedSignalPix &sig =  (*fSignals)[idx];
-      MPedestalPix        &ped =  (*fPedestals)[idx];
-      const Float_t pedrms      = ped.GetPedestalRms()*fSqrtHiGainSamples;
-      const Float_t sumhi       = sig.GetExtractedSignalHiGain();
-
-      //
-      // We consider a pixel as presumably due to cosmics 
-      // if its sum of FADC slices is lower than 3 pedestal RMS
-      //
-      if (sumhi < 3.*pedrms )  
-        cosmicpix++;
-    }
-  
-
-  //
-  // If the camera contains more than fMaxEmptyPixels
-  // presumed pixels due to cosmics, then the event is discarted. 
-  //
-  if (cosmicpix > fMaxEmptyPixels)
-    return kTRUE;
-
-  return kFALSE;
+Bool_t MFCosmics::CosmicsRejection() const
+{
+    MRawEvtPixelIter pixel(fRawEvt);
+
+    Int_t cosmicpix = 0;
+
+    //
+    // Create a first loop to sort out the cosmics ...
+    //
+    while (pixel.Next())
+    {
+        const UInt_t idx = pixel.GetPixelId();
+
+        MExtractedSignalPix &sig =  (*fSignals)[idx];
+        MPedestalPix        &ped =  (*fPedestals)[idx];
+
+        const Float_t pedrms = ped.GetPedestalRms()*fSqrtHiGainSamples;
+        const Float_t sumhi  = sig.GetExtractedSignalHiGain();
+
+        //
+        // We consider a pixel as presumably due to cosmics
+        // if its sum of FADC slices is lower than 3 pedestal RMS
+        //
+        if (sumhi < 3.*pedrms )
+            cosmicpix++;
+    }
+
+    //
+    // If the camera contains more than fMaxEmptyPixels
+    // presumed pixels due to cosmics, then the event is discarted.
+    //
+    return cosmicpix > fMaxEmptyPixels;
 }
 
 Int_t MFCosmics::PostProcess()
 {
-
-  if (GetNumExecutions()==0)
-    return kTRUE;
-  
-  *fLog << inf << endl;
-  *fLog << GetDescriptor() << " execution statistics:" << endl;
-  *fLog << dec << setfill(' ');
-  
-  *fLog << " " << setw(7) << fCut[1] << " (" << setw(3) ;
-  *fLog << (int)(fCut[1]*100/GetNumExecutions()) ;
-  *fLog << "%) Evts skipped due to: Cosmics Rejection applied " ;
-  *fLog << " (with fMaxEmptyPixels = " << fMaxEmptyPixels << ")" << endl;
-
-  *fLog << " " << setw(7) << fCut[0] << " (" << setw(3) ;
-  *fLog << (int)(fCut[0]*100/GetNumExecutions()) ;
-  *fLog << "%) Evts survived the cosmics rejection!" << endl;
-  *fLog << endl;
-
-  return kTRUE;
-}
-
+    if (GetNumExecutions()==0)
+        return kTRUE;
+
+    *fLog << inf << endl;
+    *fLog << GetDescriptor() << " execution statistics:" << endl;
+    *fLog << dec << setfill(' ');
+
+    *fLog << " " << setw(7) << fCut[1] << " (" << setw(3) ;
+    *fLog << (int)(fCut[1]*100/GetNumExecutions()) ;
+    *fLog << "%) Evts skipped due to: Cosmics Rejection applied " ;
+    *fLog << " (with fMaxEmptyPixels = " << fMaxEmptyPixels << ")" << endl;
+
+    *fLog << " " << setw(7) << fCut[0] << " (" << setw(3) ;
+    *fLog << (int)(fCut[0]*100/GetNumExecutions()) ;
+    *fLog << "%) Evts survived the cosmics rejection!" << endl;
+    *fLog << endl;
+
+    return kTRUE;
+}
+

trunk/MagicSoft/Mars/mfilter/MFCosmics.h

@@ -1 +1 @@
 #ifndef MARS_MFCosmics
 #define MARS_MFCosmics
-
-/////////////////////////////////////////////////////////////////////////////
-//                                                                         //
-// MFCosmics                                                               //
-//                                                                         //
-// Filter against cosmics (used especially by the calibration              //
-//                                                                         //
-/////////////////////////////////////////////////////////////////////////////
 
 #ifndef MARS_MFilter
@@ -22 +14 @@
 {
 private:
+    MPedestalCam        *fPedestals; // Pedestals of all pixels in the camera
+    MExtractedSignalCam *fSignals;   // Calibration events of all pixels in the camera
 
-  MPedestalCam             *fPedestals;    // Pedestals of all pixels in the camera
-  MExtractedSignalCam      *fSignals;      // Calibration events of all pixels in the camera
+    MRawEvtData         *fRawEvt;    // raw event data (time slices)
 
-  MRawEvtData              *fRawEvt;       // raw event data (time slices)
+    Int_t   fCut[2];
+    Bool_t  fResult;
 
-  Int_t   fCut[2];   
-  Bool_t  fResult;
-  Int_t   fMaxEmptyPixels;                 // Maximum number of empty pixels before declaring event as cosmic
-  Float_t fSqrtHiGainSamples;              // Square root of the number of used Hi-Gain Samples
-  
-  Bool_t ReInit(MParList *pList); 
-  Int_t PreProcess(MParList *pList);
-  Int_t Process();
-  Int_t PostProcess();
+    Int_t   fMaxEmptyPixels;         // Maximum number of empty pixels before declaring event as cosmic
+    Float_t fSqrtHiGainSamples;      // Square root of the number of used Hi-Gain Samples
 
-  Bool_t CosmicsRejection();
-  
-  Bool_t IsExpressionTrue() const { return fResult; }
+    Bool_t ReInit(MParList *pList);
+    Int_t  PreProcess(MParList *pList);
+    Int_t  Process();
+    Int_t  PostProcess();
+
+    Bool_t CosmicsRejection() const;
+
+    Bool_t IsExpressionTrue() const { return fResult; }
   
 public:
+    MFCosmics(const char *name=NULL, const char *title=NULL);
 
-  MFCosmics(const char *name=NULL, const char *title=NULL);
+    void  SetMaxEmptyPixels(const Int_t n) { fMaxEmptyPixels = n; }
+    Int_t GetMaxEmptyPixels() const        { return fMaxEmptyPixels; }
 
-  void SetMaxEmptyPixels(const Int_t n)             { fMaxEmptyPixels = n; }
-  Int_t GetMaxEmptyPixels()          const    { return fMaxEmptyPixels; }
-
-  ClassDef(MFCosmics, 0)   // Filter to perform a cosmics rejection
+    ClassDef(MFCosmics, 0)   // Filter to perform a cosmics rejection
 };