MBadPulse.cc

Go to the documentation of this file.

#include "MBadPulse.hh"
#include "QEvent.hh"
#include "QCountPulsesData.hh"
#include "QBaselineData.hh"
#include "QPulseParameters.hh"
 
using namespace Diana;
 
REGISTER_MODULE(MBadPulse)
 
 
void MBadPulse::GetCuts(const QEvent& ev)
{
        fMaxBaselineSlope = GetDouble("MaxBaselineSlope", 0);
    units = GetString("Units","Absolute",false);
    if(units.compare("Absolute")!=0 and units.compare("Relative")!=0)Panic("Units not supported, can choose between 'Absolute' (mV/sample) or 'Relative' (BaselineRMS/WindowLength)");
    // If MinBaselineSlope is specified, fMinBaselineSlope < BaselineSlope < fMaxBaselineSlope; 
    // else, -fMaxBaselineSlope < BaselineSlope < fMaxBaselineSlope
    fMinBaselineSlope = GetDouble("MinBaselineSlope", -fMaxBaselineSlope);
        fBaselineOwner = GetString("BaselineOwner","BaselineModule",false);
        if(fMaxBaselineSlope) {
                ev.Require<QBaselineData>(fBaselineOwner,"BaselineData");
        }
 
    fFullWindowBaselineSlope = GetDouble("MaxFullWindowBaselineSlope", 0, false);
    fFullWindowBaselineOwner = GetString("FullWindowBaselineOwner", "BaselineModule_FullWindow", false);
    if(fFullWindowBaselineSlope) {
        ev.Require<QBaselineData>(fFullWindowBaselineOwner, "BaselineData");
    }
 
        fnPeaks = GetInt("NumberOfPeaks", -1);
        fCountPulsesOwner = GetString("CountPulsesOwner","BCountPulses",false);
        if(fnPeaks >= 0) {
                ev.Require<QCountPulsesData>(fCountPulsesOwner,"CountPulsesData");
        }
 
        fPulseBasicParametersOwner = GetString("PulseBasicParametersOwner","PulseBasicParameters",false);
        ev.Require<QPulseParameters>(fPulseBasicParametersOwner,"Parameters");
        fCheckSaturation  = GetBool("CheckSaturation",true,false);
}
 
bool MBadPulse::Filter(const QEvent& ev)
{
        // perform selection: true if the event passes the cuts
        bool result = true;
        if(fCheckSaturation) { 
                //Note - May 2018:
                //As the variables IsSaturatedLow/IsSaturatedHigh are assigned up to now, this check only rejects events which have baselines which saturated the ADC range
                const QPulseParameters& bParams = ev.Get<QPulseParameters>(fPulseBasicParametersOwner.c_str(),"Parameters");
                if(bParams.fIsSaturatedLow || bParams.fIsSaturatedHigh) result = false;
        }
 
        if(fnPeaks >= 0)
        {
                if(ev.Get<QCountPulsesData>(fCountPulsesOwner.c_str(),"CountPulsesData").GetNumberOfPulses() != fnPeaks) 
                        result = false;
        }
        if(fMaxBaselineSlope){
                double bs1=0;
        if(units.compare("Absolute")==0){
            bs1= ev.Get<QBaselineData>(fBaselineOwner.c_str(),"BaselineData").GetBaselineSlope();
        }else if(units.compare("Relative")==0){
            bs1= ev.Get<QBaselineData>(fBaselineOwner.c_str(),"BaselineData").GetBaselineSlopeRMSWindow();
        }    
                if(bs1 > fMaxBaselineSlope || bs1 < fMinBaselineSlope)
                        result = false;
        }
    if(fFullWindowBaselineSlope)
    {
                double bs2=0;
        if(units.compare("Absolute")==0){
            bs2 = ev.Get<QBaselineData>(fFullWindowBaselineOwner.c_str(),"BaselineData").GetBaselineSlope();
        }else if(units.compare("Relative")==0){
            bs2= ev.Get<QBaselineData>(fFullWindowBaselineOwner.c_str(),"BaselineData").GetBaselineSlopeRMSWindow();
        }    
        if(bs2 > fFullWindowBaselineSlope || bs2 < (-1*fFullWindowBaselineSlope))
            result = false;
    }
    
        return result;
}