MPulseIntegral.cc

Go to the documentation of this file.

#include "MPulseIntegral.hh"
#include "QEvent.hh"
#include "QEventList.hh"
#include "QBaselineData.hh"
#include "QPulse.hh"
#include "QBaseType.hh"
#include "QHeader.hh"
#include "QPulseInfo.hh"
#include "QRunDataHandle.hh"
 
 
REGISTER_MODULE(MPulseIntegral)
 
using namespace Diana;
 
void MPulseIntegral::Init(QEvent& ev)
{
    fPulseLabel = GetString("PulseLabel","DAQ@Pulse",false);
    fBaselineOwner = GetString("BaselineOwner","BaselineModule",false);
    ev.RequireByLabel<QPulse>(fPulseLabel);
    ev.Require<QBaselineData>(fBaselineOwner,"BaselineData");
    ev.Require<QHeader>("DAQ","Header");
    ev.Require<QPulseInfo>("DAQ","PulseInfo");
 
    ev.Add<QDouble>("Integral");
 
    fNormalizeAmplitude = GetBool("NormalizeAmplitude",false,false);
    if(fNormalizeAmplitude) {
        fAmplitudeLabel = GetString("AmplitudeLabel",Q_STRING_DEFAULT); 
        ev.RequireByLabel<QDouble>(fAmplitudeLabel);
        ev.Add<QDouble>("NormalizedAmplitude");
    }
 
    fTomV = GetBool("ConvertTomV",true,false);
}
 
void MPulseIntegral::Do(QEvent& ev)
{
    const QHeader& header = ev.Get<QHeader>("DAQ","Header");
    const int channel = ev.Get<QPulseInfo>("DAQ","PulseInfo").GetChannelId();
    int run = header.GetRun();
    QRunDataHandle rHandle(run);
    GlobalData().Get("",&rHandle,"");
    const QRunData& runData = rHandle.Get();
    const QChannelRunData& chanRunData = runData.GetChannelRunData(channel);
 
    double ADC2mV = 1;
    if(fTomV) ADC2mV = chanRunData.fADC2mV;
    const double samplingFreq = chanRunData.fSamplingFrequency;
 
    const QVector& samples = ev.GetByLabel<QPulse>(fPulseLabel).GetSamples()*ADC2mV;
    const double baseline = ev.Get<QBaselineData>(fBaselineOwner.c_str(),"BaselineData").GetBaseline();
    const size_t N = samples.Size();
    double integral = samples.Sum(N,0)-baseline*N;
    integral /= samplingFreq;
    ev.Get<QDouble>("Integral") = integral;
 
    if(fNormalizeAmplitude) {
        ev.Get<QDouble>("NormalizedAmplitude") = ev.GetByLabel<QDouble>(fAmplitudeLabel)/integral;
    }
 
 
}
 
void MPulseIntegral::Done()
{
}