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

Timestamp:

02/23/04 00:10:35 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mfilter/MFCosmics.cc (modified) (7 diffs)

trunk/MagicSoft/Mars/mfilter/MFCosmics.cc

@@ -28 +28 @@
 //
 //   Filter to reject cosmics by the criterion that at least 
-//   fMaxEmptyPixels pixels have values of lower than 3 Pedestal RMS. 
-//   fMaxEmptyPixels is set to 230 by default which is slightly higher 
+//   fMaxEmptyPixels of the pixels have values of lower than 3 Pedestal RMS. 
+//   fMaxEmptyPixels is set to 0.4 by default which is slightly higher 
 //   than the number of outer pixels in MAGIC (for the case that 
 //   the outer pixels have some defect).
@@ -67 +67 @@
 MFCosmics::MFCosmics(const char *name, const char *title)
     : fPedestals(NULL), fSignals(NULL),
-      fRawEvt(NULL), fMaxEmptyPixels(230)
+      fRawEvt(NULL), fMaxEmptyPixels(0.4)
 {
     fName  = name  ? name  : "MFCosmics";
@@ -102 +102 @@
         return kFALSE;
     }
+
 
     memset(fCut, 0, sizeof(fCut));
@@ -150 +151 @@
 
     Int_t cosmicpix = 0;
+    Int_t allpix    = 0;
 
     //
@@ -164 +166 @@
         const Float_t sumhi  = sig.GetExtractedSignalHiGain();
 
+        allpix++;
         //
         // We consider a pixel as presumably due to cosmics
@@ -176 +179 @@
     // presumed pixels due to cosmics, then the event is discarted.
     //
-    return cosmicpix > fMaxEmptyPixels;
+    return cosmicpix > fMaxEmptyPixels*allpix;
 }
 
@@ -188 +191 @@
     *fLog << dec << setfill(' ');
 
+    *fLog << " " << setw(7) << fCut[0] << " (" << setw(3) ;
+    *fLog << (int)(fCut[0]*100/GetNumExecutions()) ;
+    *fLog << "%) Evts skipped due to: Cosmics Rejection applied " ;
+    *fLog << " (with fMaxEmptyPixels = " << fMaxEmptyPixels*100 << "%)" << endl;
+
     *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;