MFilterInIntervalI.cc

Go to the documentation of this file.

#include "MFilterInIntervalI.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(MFilterInIntervalI)
 
void MFilterInIntervalI::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 = GetInt("MinValue", Q_INT_DEFAULT);
    fMaxValue = GetInt("MaxValue", Q_INT_DEFAULT);
 
    fError = false;
    if (fMaxValue == Q_INT_DEFAULT && fMinValue == Q_INT_DEFAULT) {
        Error("No range specified. The variable (%s) will not be filtered", fVariableLabel.GetStringLabel().c_str());
        fError = true;
    }
    else if (fMaxValue != Q_INT_DEFAULT && fMinValue != Q_INT_DEFAULT && fMaxValue < fMinValue) {
        Error("Bad range (max < min). The variable (%s) will not be filtered", fVariableLabel.GetStringLabel().c_str());
        fError = true;
    }
 
    ev.RequireByLabel<QInt>(fVariableLabel);
}
 
bool MFilterInIntervalI::Filter(const Diana::QEvent& ev)
{
    // get variable from event
    const int variable = ev.GetByLabel<QInt>(fVariableLabel);
 
    if (fError)
        return false;
    if (fMinValue == Q_INT_DEFAULT && fMaxValue != Q_INT_DEFAULT)
        return variable <= fMaxValue; // applies only upper limit
    if (fMaxValue == Q_INT_DEFAULT && fMinValue != Q_INT_DEFAULT)
        return variable >= fMinValue; // applies only lower limit
 
    return (variable >= fMinValue && variable <= fMaxValue);
}