MFilterSideAmplitude.cc

Go to the documentation of this file.

#include "MFilterSideAmplitude.hh"
#include "QEvent.hh"
#include "QEventList.hh"
#include "QBaseType.hh"
#include <cmath>
 
 
REGISTER_MODULE(MFilterSideAmplitude)
 
using namespace Diana;
 
 
void MFilterSideAmplitude::GetCuts(const QEvent& ev)
{
  fNSigma = GetDouble("NumberOfSigma",3,false);
  fSideAmplitudeOwner = GetString("SideAmplitudeOwner","SideAmplitude",false);
  
  ev.RequireByLabel<QDouble>(QEventLabel(fSideAmplitudeOwner.c_str(),"Amplitude"));
  ev.RequireByLabel<QDouble>(QEventLabel(fSideAmplitudeOwner.c_str(),"NSigma"));
}
 
bool MFilterSideAmplitude::Filter(const QEvent& ev, const QEventList& neighbours)
{
  
  bool passed = true;
 
  for(size_t i = 0; i < neighbours.Size(); i++) {
    const QDouble& n_amplitude = neighbours[i].GetByLabel<QDouble>(QEventLabel(fSideAmplitudeOwner.c_str(),"Amplitude"));
    const QDouble& n_nsigma = neighbours[i].GetByLabel<QDouble>(QEventLabel(fSideAmplitudeOwner.c_str(),"NSigma"));
    if (n_amplitude>0. && n_nsigma>fNSigma) passed = false;
  }
 
  return passed;
}