Changeset 2282 for trunk/MagicSoft/Mars/mfilter

Timestamp:

07/17/03 08:06:40 (21 years ago)

Author:

wittek

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mfilter

Files:

• MFCT1SelBasic.cc (modified) (7 diffs)
• MFCT1SelBasic.h (modified) (3 diffs)

trunk/MagicSoft/Mars/mfilter/MFCT1SelBasic.cc

@@ -30 +30 @@
 //  This is a class to evaluate basic cuts
 //
-//  WHAT ARE THE BASIC CUTS?
-//
 //  to be called after the calibration (when the number of photons is
 //               available for all pixels)
 //
+//  The basic cuts are :
+//
+//      remove bad runs
+//      thetamin < theta < thetamax
+//      software trigger fullfilled (with minimum no.of photons = minphotons)
+//
+//
 /////////////////////////////////////////////////////////////////////////////
 
@@ -44 +49 @@
 
 #include "MCerPhotEvt.h"
-//nclude "MRawRunHeader.h"
+#include "MRawRunHeader.h"
 
 #include "MGeomPix.h"
@@ -95 +100 @@
 Int_t MFCT1SelBasic::PreProcess(MParList *pList)
 {
-    /*
     fRawRun = (MRawRunHeader*)pList->FindObject("MRawRunHeader");
     if (!fRawRun)
@@ -102 +106 @@
         return kFALSE;
     }
-    */
 
     fMcEvt = (MMcEvt*)pList->FindObject("MMcEvt");
@@ -148 +151 @@
 // Evaluate basic cuts
 // 
-// if cuts are fulfilled      : fResult = kTRUE;
-// if they are not fullfilled : fResult = kFALSE;
+//     bad events    : fResult = kTRUE;
+//     good events   : fResult = kFALSE;
 //
 Int_t MFCT1SelBasic::Process()
@@ -157 +160 @@
     fResult  = kFALSE;
 
+    // remove bad runs for MC gammas
+    if (fMcEvt->GetEnergy() == 0.0  &&  fMcEvt->GetImpact() == 0.0)
+    {
+      if (fRawRun->GetRunNumber() == 601  ||
+          fRawRun->GetRunNumber() == 613  ||
+          fRawRun->GetRunNumber() == 614    )
+        return Set(1);
+    }
+
     if (theta<fThetaMin)
-        return Set(1);
+        return Set(2);
 
     if (theta>fThetaMax)
-        return Set(2);
+        return Set(3);
 
     if (!SwTrigger())
-        return Set(3);
+        return Set(4);
 
     fCut[0]++;
@@ -234 +246 @@
     *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: Zenith angle < " << fThetaMin << endl;
+    *fLog << "%) Evts skipped due to: bad run " << endl;
 
     *fLog << " " << setw(7) << fCut[2] << " (" << setw(3) ;
     *fLog << (int)(fCut[2]*100/GetNumExecutions()) ;
-    *fLog << "%) Evts skipped due to: Zenith angle > " << fThetaMax << endl;
+    *fLog << "%) Evts skipped due to: Zenith angle < " << fThetaMin << endl;
 
     *fLog << " " << setw(7) << fCut[3] << " (" << setw(3) ;
     *fLog << (int)(fCut[3]*100/GetNumExecutions()) ;
+    *fLog << "%) Evts skipped due to: Zenith angle > " << fThetaMax << endl;
+
+    *fLog << " " << setw(7) << fCut[4] << " (" << setw(3) ;
+    *fLog << (int)(fCut[4]*100/GetNumExecutions()) ;
     *fLog << "%) Evts skipped due to: Software trigger not fullfilled" ;
     *fLog << " (with fMinPhotons = " << fMinPhotons << ")" << endl;

trunk/MagicSoft/Mars/mfilter/MFCT1SelBasic.h

@@ -18 +18 @@
 class MCerPhotEvt;
 class MPedestalCam;
-//class MRawRunHeader;
+class MRawRunHeader;
 
 class MFCT1SelBasic : public MFilter
@@ -27 +27 @@
     const MCerPhotEvt   *fEvt;      // Cerenkov Photon Event 
     //const MPedestalCam  *fPed;      // Pedestal information
-    //const MRawRunHeader *fRawRun;
+    const MRawRunHeader *fRawRun;
 
     Float_t     fMinPhotons;
@@ -33 +33 @@
     Float_t     fThetaMax;
 
-    Int_t       fCut[4];
+    Int_t       fCut[5];
 
     Bool_t      fResult;