MTimeCutFilter.cc

Go to the documentation of this file.

#include "MTimeCutFilter.hh"
#include "QEvent.hh"
#include "QEventList.hh"
#include "QRawEvent.hh"
 
//MinTime and MaxTime must be given in seconds to the cfg file
 
REGISTER_MODULE(MTimeCutFilter)
 
using namespace Diana;
 
 
void MTimeCutFilter::GetCuts(const QEvent& ev)
{
 
  fMinTime = GetDouble("MinTime",(double)Q_DOUBLE_DEFAULT);
  fMaxTime = GetDouble("MaxTime",(double)Q_DOUBLE_DEFAULT);
 
  fIsWarned = false;
 
}
 
 
bool MTimeCutFilter::Filter(const QEvent& ev)
{
 
  bool passed = false;
 
  QHeader header = ev.Get<QHeader>("DAQ","Header");
  double time = header.GetTime().GetFromStartRunNs()/1.e09;
 
 
  if(fMinTime==Q_DOUBLE_DEFAULT && fMaxTime!=Q_DOUBLE_DEFAULT)  
    {
      // applies only upper limit
      if (time <= fMaxTime) passed = true;
    }
 
 if(fMinTime!=Q_DOUBLE_DEFAULT && fMaxTime==Q_DOUBLE_DEFAULT)  
   {
     // applies only lower limit
     if (time >= fMinTime) passed = true;
   }
 
  // no range specified
 if(fMaxTime==Q_DOUBLE_DEFAULT && fMinTime==Q_DOUBLE_DEFAULT && !fIsWarned) 
   {
     Warn("No range specified. The variable will not be filtered");
     fIsWarned= true;
   }
 
 // wrong range (max < min)
 if(fMaxTime!=Q_DOUBLE_DEFAULT && fMinTime!=Q_DOUBLE_DEFAULT && fMaxTime <= fMinTime && !fIsWarned) 
   {
     Warn("Wrong Range (Max <= Min). The variable will not be filtered");
     fIsWarned = true;
   }
 
 else
   {
     if(time >= fMinTime && time <= fMaxTime) passed = true;
   }
 
 
  return passed;
}