MGoodWaveform.cc

Go to the documentation of this file.

#include "MGoodWaveform.hh"
#include "QEvent.hh"
#include "QCountPulsesData.hh"
#include "QBaselineData.hh"
#include "QPulseParameters.hh"
#include "QPulseInfo.hh"
 
using namespace Diana;
 
REGISTER_MODULE(MGoodWaveform)
 
 
void MGoodWaveform::GetCuts(const QEvent& ev)
{
    fnPeaks = GetInt("NumberOfPeaks", -1);
    fBaselineSlope = GetDouble("MaxBaselineSlope", 0, false);
    ev.Require<QCountPulsesData>("BCountPulses","CountPulsesData");
    ev.Require<QBaselineData>("BaselineModule","BaselineData");
    ev.Require<QPulseParameters>("PulseBasicParameters","Parameters");
    ev.Require<QPulseInfo>("DAQ","PulseInfo");
}
 
bool MGoodWaveform::Filter(const QEvent& ev)
{
    // perform selection: true if the event passes the cuts
        bool result = true;
    const QPulseParameters& bParams = ev.Get<QPulseParameters>("PulseBasicParameters","Parameters");
    if(bParams.fIsSaturatedLow || bParams.fIsSaturatedHigh) return false;
        if(fBaselineSlope)
        {
                double bs = ev.Get<QBaselineData>("BaselineModule","BaselineData").GetBaselineSlope();
                if(bs > fBaselineSlope || bs < (-1*fBaselineSlope))
                        result = false;
        }
 
 
    const QPulseInfo& pi = ev.Get<QPulseInfo>("DAQ","PulseInfo");
    int numberOfPulses = ev.Get<QCountPulsesData>("BCountPulses","CountPulsesData").GetNumberOfPulses();
 
    if(pi.GetIsNoise()) {
        if(numberOfPulses > 0) result = false;
    } else if(pi.GetIsSignal() || pi.GetIsPulser()) {
        if(fnPeaks >= 0  && numberOfPulses != fnPeaks) result = false;
    } else {
        result = false;
    }
 
    return result;
}