MSampleInfoFilter.cc

Go to the documentation of this file.

#include "MSampleInfoFilter.hh"
 
#include "QEvent.hh"
#include "QSampleInfo.hh"
#include "QRawEvent.hh"
#include "QChannelRunData.hh"
#include "QRunDataHandle.hh"
 
 
REGISTER_MODULE(MSampleInfoFilter)
 
 
void MSampleInfoFilter::GetCuts(const Diana::QEvent& ev)
{
    // Get parameters from cfg
    fType = GetString("Type","Noise");
    fKeepType = GetBool("KeepType",true); 
    fTimeThreshold = GetDouble("PulserTimeThreshold", 0.04);
    ev.Require<QPulseInfo>("DAQ","PulseInfo");
    ev.Require<QHeader>("DAQ","Header");
    return;
}
 
bool MSampleInfoFilter::Filter(const Diana::QEvent& ev)
{
    const QHeader& header = ev.Get<QHeader>("DAQ","Header");
    const QPulseInfo& pulseInfo = ev.Get<QPulseInfo>("DAQ","PulseInfo");
    const int chan = pulseInfo.GetChannelId();
    QRunDataHandle rHandle(header.GetRun());
    GlobalData().Get("",&rHandle,"");
    const QRunData& runData = rHandle.Get();
    const QChannelRunData& channelRunData = runData.GetChannelRunData(chan);
 
    double sampFreq = channelRunData.fSamplingFrequency;
 
    // Within this threshold, pulser triggers and corresponding derivative triggers are considered as one
    double sampleThreshold = fTimeThreshold * sampFreq;
    
    const std::vector<QSampleInfo>& samplesInfo = pulseInfo.GetSamplesInfo();
    if(!samplesInfo.empty()) {
        std::vector<QSampleInfo>::const_iterator sampleInfo = samplesInfo.begin();
        while(sampleInfo != samplesInfo.end()) {
            const int trigpos = pulseInfo.GetMasterSample().GetSampleIndex();
            bool same = (abs(trigpos-(int)sampleInfo->GetSampleIndex())<sampleThreshold);
            if(fKeepType){
                if(fType == "Noise" && !sampleInfo->GetIsNoise()) return false;
                else if(fType == "Signal" && !sampleInfo->GetIsSignal() && same) return false;
                else if(fType == "Heater" && !sampleInfo->GetIsPulser() && same) return false;
                else if(fType == "LED" && !sampleInfo->GetIsLED() && same) return false;
            }else{
                if(fType == "Noise" && sampleInfo->GetIsNoise()) return false;
                else if(fType == "Signal" && sampleInfo->GetIsSignal()) return false;
                else if(fType == "Heater" && sampleInfo->GetIsPulser()) return false;
                else if(fType == "LED" && sampleInfo->GetIsLED()) return false;
            }
            sampleInfo++;
        }
    }
    return true;
}