QTimeNeighbour.cc

Go to the documentation of this file.

#include "QTimeNeighbour.hh"
#include "QReader.hh"
#include "QEvent.hh"
#include "QRawEvent.hh"
#include "QBool.hh"
 
Q_BEGIN_NAMESPACE
 
void QTimeNeighbour::Init(const QReader& reader)
{
    fWindowSizeMs = reader.GetDouble("WindowSizeMs",200);
    fNeighbourDistanceMs = reader.GetDouble("NeighbourDistanceMs",100);
    fNeighbourOffsetMs = reader.GetDouble("NeighbourOffsetNs",0);
    fFiltVar = reader.GetString("MasterFilteringVariable","",false);
    fOffsetNs = fNeighbourOffsetMs*1e6;
}
 
QReaderNeighbour::NeighbourInfo QTimeNeighbour::GetNeighbourInfo(const QEvent& thisEvent, const QEvent& otherEvent)
{
 
    QReaderNeighbour::NeighbourInfo info;
    const QHeader& thisHeader = thisEvent.Get<QHeader>("DAQ","Header");
    const QHeader& otherHeader = otherEvent.Get<QHeader>("DAQ","Header");
 
 
    unsigned long long stime = thisHeader.GetTime().GetStartRunUnix();
    unsigned long long sothertime = otherHeader.GetTime().GetStartRunUnix();
 
    //If different start runs then they are not neighbours (debatable)
    if(stime != sothertime) {
        info.fIsInWindow = false;
        info.fIsNeighbour = false;
        return info;
    }
 
 
 
    long long time = thisHeader.GetTime().GetFromStartRunNs()+fOffsetNs;
    long long othertime = otherHeader.GetTime().GetFromStartRunNs();
 
    //da mettere qui l'offset!
    double diff = llabs(time-othertime);
    diff /= 1.e06;
    int channel =  thisEvent.Get<QPulseInfo>("DAQ","PulseInfo").GetChannelId();     
    int ochannel =  otherEvent.Get<QPulseInfo>("DAQ","PulseInfo").GetChannelId();     
 
    info.fIsInWindow = false;
    info.fIsNeighbour = false;
 
    if(diff < fWindowSizeMs) {
        info.fIsInWindow = true;
        if(diff < fNeighbourDistanceMs && channel != ochannel) info.fIsNeighbour = true;
    }
    return info;
 
 
}
 
bool QTimeNeighbour::GetWantsNeighbours(const QEvent& thisEvent){
    return (bool)thisEvent.GetByLabel<Diana::QBool>(fFiltVar);
}
 
Q_END_NAMESPACE