MBaselineFilter.cc

Go to the documentation of this file.

#include "MBaselineFilter.hh"
#include "QEvent.hh"
#include "QBaseType.hh"
#include "QBaselineData.hh"
#include <map>
#include <set>
#include <string>
 
using std::map;
using std::set;
using std::string;
using std::cerr;
using std::endl;
 
using namespace Diana;
 
REGISTER_MODULE(MBaselineFilter)
 
void MBaselineFilter::GetCuts(const Diana::QEvent& ev)
{
  // Get parameters from cfg
  
  // define minimum and maximum value for the variable
  fMinValue = GetDouble("MinValue",Q_DOUBLE_DEFAULT);
  fMaxValue = GetDouble("MaxValue",Q_DOUBLE_DEFAULT);
  
  fBaselineOwner = GetString("BaselineOwner","BaselineModule",false);
  
  fVariable = GetString("Variable","Baseline",false);
 
  fIsWarned = false;
 
  ev.Require<QBaselineData>((string)fBaselineOwner,"BaselineData");
}
 
bool MBaselineFilter::Filter(const Diana::QEvent& ev)
{
 
  const QBaselineData& bData = ev.Get<QBaselineData>(((string)fBaselineOwner).c_str(),"BaselineData");
 
  // get variable from event
  
  Double_t variable;
  
  if(fVariable.CompareTo("Baseline")==0){
        variable = bData.GetBaseline();
  
  }else if(fVariable.CompareTo("BaselineFlatRMS")==0){
        variable = bData.GetBaselineFlatRMS();
  
  }else if(fVariable.CompareTo("BaselineIntercept")==0){
        variable = bData.GetBaselineIntercept();
  
  }else if(fVariable.CompareTo("BaselineSlope")==0){
        variable = bData.GetBaselineSlope();
  
  }else if(fVariable.CompareTo("BaselineSlopeRMSWindow")==0){
        variable = bData.GetBaselineSlopeRMSWindow();
  
  }else if(fVariable.CompareTo("BaselineRMS")==0){
        variable = bData.GetBaselineRMS();
  
  }else if(fVariable.CompareTo("RightBaseline")==0){
        variable = bData.GetRightBaseline();
  
  }else if(fVariable.CompareTo("RightBaselineRMS")==0){
        variable = bData.GetRightBaselineInRMS();
  }else if(fVariable.CompareTo("RightLeftBaseline")==0){
        variable = bData.GetRightLeftBaseline();
  }else if(fVariable.CompareTo("RightLeftBaselineInRMS")==0){
        variable = bData.GetRightLeftBaselineInRMS();
  
  }else{
        Panic(TString("'")+fVariable+TString("' if not a valid choice for the 'Variable' field."));
  }
  
  
  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;
 
}