MFilterInInterval.cc

Go to the documentation of this file.

#include "MFilterInInterval.hh"
#include "QEvent.hh"
#include "QBaseType.hh"
#include <map>
#include <set>
#include <string>
 
using std::map;
using std::set;
using std::string;
 
using namespace Diana;
 
REGISTER_MODULE(MFilterInInterval)
 
void MFilterInInterval::GetCuts(const Diana::QEvent& ev)
{
  // Get parameters from cfg
  
  // define variable label (variable on which the cut is applied)
  fVariableLabel = GetString("VariableLabel","");
 
  // define minimum and maximum value for the variable
  fMinValue = GetDouble("MinValue",Q_DOUBLE_DEFAULT);
  fMaxValue = GetDouble("MaxValue",Q_DOUBLE_DEFAULT);
 
  fIsWarned = false;
 
  ev.RequireByLabel<QDouble>(fVariableLabel);
}
 
bool MFilterInInterval::Filter(const Diana::QEvent& ev)
{
  // get variable from event
  const double variable = ev.GetByLabel<QDouble>(fVariableLabel);
  
  bool passed = false;
  
  // perform selection: true if the event passes the cuts  
  if(fMinValue==Q_DOUBLE_DEFAULT && fMaxValue!=Q_DOUBLE_DEFAULT)  
    {
      // applies only upper limit
      if (variable<=fMaxValue) passed = true;
    }
  if(fMaxValue==Q_DOUBLE_DEFAULT && fMinValue!=Q_DOUBLE_DEFAULT)  
    {
      // applies only lower limit
      if (variable>=fMinValue) passed = true;
    }
  // no range specified
  if(fMaxValue==Q_DOUBLE_DEFAULT && fMinValue==Q_DOUBLE_DEFAULT && !fIsWarned) 
    {
      Warn("No range specified. The variable will not be filtered");
      fIsWarned= true;
    }
  // wrong range (max < min)
  if(fMaxValue!=Q_DOUBLE_DEFAULT && fMinValue!=Q_DOUBLE_DEFAULT && fMaxValue <= fMinValue && !fIsWarned) 
    {
      Warn("Wrong Range (Max %f <= Min %f). The variable will not be filtered",fMaxValue,fMinValue);
      fIsWarned= true;
    }
  else 
    {
      if (variable>fMinValue && variable<fMaxValue) passed = true;
    }
 
  return passed;
 
}