QCreateSidePulseByPulseSplitting.cc

Go to the documentation of this file.

#include "QCreateSidePulseByPulseSplitting.hh"
#include "QBool.hh"
#include "QReader.hh"
#include "QEvent.hh"
#include "QHeader.hh"
#include "QPulse.hh"
#include "QPulseFiller.hh"
#include "QPulseInfo.hh"
#include "QSidePulses.hh"
#include "QEventAssembler.hh"
#include "QCachedEventList.hh"
#include "QGlobalDataManager.hh"
#include "QMessage.hh"
#include <sstream>
 
Q_BEGIN_NAMESPACE
 
QCreateSidePulseByPulseSplitting::QCreateSidePulseByPulseSplitting()
{
}
 
void QCreateSidePulseByPulseSplitting::Init(const QReader& reader,QEventAssembler &eva)
{
    std::vector<std::string> def;
    def.push_back("DAQ@Pulse");
    fPulseLabels = reader.GetVectorString("PulseLabels",def,false);
    for(size_t p = 0; p < fPulseLabels.size(); p++) {
        std::string owner = "CreateSidePulseByPulseSplitting_";
        QEventLabel label(fPulseLabels[p]);
        owner += label.owner;
        eva.Add<QPulse>(owner.c_str(),label.name.c_str());
        eva.Get<QPulse>(owner.c_str(),label.name.c_str()).SetIsJustAdded();
        eva.GetEvent().RequireByLabel<QPulse>(label);
    }
 
}
 
void QCreateSidePulseByPulseSplitting::FillList(const QEvent& ev, QCachedEventList& list)
{
    const QHeader& header = ev.Get<QHeader>("DAQ","Header");
    const QPulseInfo& pulseInfo = ev.Get<QPulseInfo>("DAQ","PulseInfo");
 
    // the side event has the same structure of the main event
    QEventAssembler evaA,evaB;
    evaA.Assign(ev);
    evaB.Assign(ev);
    // same structure but everything is not valid
    evaB.Clear();
 
    // set the raw event
 
        // fill Master sample
    QPulseInfo& sPulseInfo = evaB.Get<QPulseInfo>("DAQ","PulseInfo");  
    sPulseInfo.SetChannelId(pulseInfo.GetChannelId());
 
    QSampleInfo& sMaster = sPulseInfo.GetMasterSample();
    sMaster = pulseInfo.GetMasterSample();
    sMaster.SetIsSide();
    sPulseInfo.Validate();
 
 
    int pulseLength = 0;
 
    for(size_t p = 0; p < fPulseLabels.size(); p++) {
        const QPulse& pulse = ev.GetByLabel<QPulse>(fPulseLabels[p].c_str());  
        const QVectorI & samples = pulse.GetSamplesADC();
        size_t newSize = samples.Size()/2;
        int offset = 0;
        if(newSize % 2 != 0) {
            newSize--;
            offset = 1;
        }
        QVectorI samplesA(newSize);
        for(size_t i = 0; i < newSize; i++) samplesA[i] = samples[i+offset]; 
        QVectorI samplesB(newSize);
        for(size_t i = 0; i < newSize; i++) samplesB[i] = samples[i+samples.Size()/2]; 
 
        QPulse pulseA,pulseB;
        pulseA.SetSamplesADC(samplesA);
        pulseB.SetSamplesADC(samplesB);
        pulseA.Validate();
        pulseB.Validate();
        std::string owner = "CreateSidePulseByPulseSplitting_";
        QEventLabel label(fPulseLabels[p]);
        owner += label.owner;
        evaA.Get<QPulse>(owner.c_str(),label.name.c_str()) = pulseA;
        evaB.Get<QPulse>(owner.c_str(),label.name.c_str()) = pulseB;
 
        if(p == 0) pulseLength = samples.Size()/2;
    }
 
    // fill header
    QHeader& sHeader = evaB.Get<QHeader>("DAQ","Header");  
    sHeader = header;
    Long64_t triggerTime = sHeader.GetTime().GetFromStartRunNs()+pulseLength;
    sHeader.GetTime().SetFromStartRunNs(triggerTime);
 
 
 
    // Specify that this event has been generated in an event generator
    evaB.Get<QBool>("Reader","GeneratedEvent")=true;
 
    // insert side event in the list
    list.Push(evaB.GetEvent());
    // insert this event at the end of the list
    list.Push(evaA.GetEvent());
}
 
Q_END_NAMESPACE