Changeset 3112 for trunk/MagicSoft/Mars/mfilter/MFCosmics.cc

Timestamp:

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

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mfilter/MFCosmics.cc (modified) (7 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;
+}
+