MOverlappedWindowFilter.cc

Go to the documentation of this file.

#include "MOverlappedWindowFilter.hh"
#include "QEvent.hh"
#include "QRawEvent.hh"
#include "QRunDataHandle.hh"
#include <cmath>
 
using namespace Diana;
 
REGISTER_MODULE(MOverlappedWindowFilter)
/*
MOverlappedWindowFilter::MOverlappedWindowFilter(QSequence* s) : QFilter("OverlappedWindowFilter",s)
{
    // called once in a sequence
}
 
MOverlappedWindowFilter::~MOverlappedWindowFilter()
{
    // called once in a sequence
}*/
 
 
void MOverlappedWindowFilter::GetCuts(const QEvent& ev)
{
    // Get parameters from cfg
    fPulseLabel = GetString("PulseLabel","DAQ@Pulse",false);
 
    fRun = -1;
}
 
bool MOverlappedWindowFilter::Filter(const QEvent& ev)
{
 
    const QHeader& hdr = ev.Get<QHeader>("DAQ","Header");
 
    QRunDataHandle rHandle(hdr.GetRun());
    GlobalData().Get("",&rHandle,"");
    const QRunData& runData = rHandle.Get();
 
    int thisRun = runData.fNumber;
    if(thisRun != fRun) 
    {
        fChannels.clear();
        fRun = thisRun;
    }
 
    // check that this window is not overlapped with a previous one
    //    const QRawEvent& raw = ev.GetRawEvent();
    //    const int channel = raw.GetPulse().GetChannelId();
    const QPulseInfo& pulseInfo = ev.Get<QPulseInfo>("DAQ","PulseInfo");
    int channel = pulseInfo.GetChannelId();
 
    //    const unsigned long long time = ev.GetTime().GetFromStartRunNs();
    const unsigned long long time = hdr.GetTime().GetFromStartRunNs();
    const double samplFreq = runData.GetChannelRunData(channel).fSamplingFrequency;
    //    const size_t N = raw.GetPulse().GetSamples().Size();
    //const QPulse& pulse = ev.Get<QPulse>("DAQ","Pulse");
    const QPulse &pulse = ev.GetByLabel<QPulse>(fPulseLabel);
    const size_t N = pulse.GetSamples().Size();
    //    const int timedSample = pulse.GetTimedSample();
    const int timedSample = pulseInfo.GetMasterSample().GetSampleIndex();
    std::map<int, unsigned long long>::iterator chanIter = fChannels.find(channel);
 
    const unsigned long long firstTimeInWindow = time - (unsigned long long)(double(timedSample)*1.e09/samplFreq);
    const unsigned long long lastTimeInWindow = firstTimeInWindow + (unsigned long long)(double(N)*1.e09/samplFreq);
    if(chanIter != fChannels.end())  {
        if(firstTimeInWindow <= chanIter->second) {
            return false;
        } else {
            chanIter->second = lastTimeInWindow;
            return true;
        }
 
    } else {
        fChannels[channel] = lastTimeInWindow;
        return true;
    }
 
 
    return true;
}