MTutorialPulseShape.cc

Go to the documentation of this file.

 
#include "MTutorialPulseShape.hh"
#include "QEvent.hh"
#include "QRawEvent.hh"
#include "QVector.hh"
#include "QBaseType.hh"
#include "QBaselineData.hh"
 
REGISTER_MODULE(MTutorialPulseShape)
 
using namespace Diana;
 
void MTutorialPulseShape::Init(QEvent& ev)
{
    fAveragePulseInput = GetString("AveragePulseInput","NODEFAULT");
    
    // require presence in the events of quantities read by this module
    ev.Require<QPulseInfo>("DAQ","PulseInfo");
    ev.Require<QPulse>("DAQ","Pulse");
    ev.Require<QBaselineData>("BaselineModule","BaselineData");
 
    // add to event quantities written by this module.
    ev.Add<QDouble>("ChiSquare");
}
 
void MTutorialPulseShape::Do(QEvent& ev)
{
    // get channel number and samples from raw event
    const QPulseInfo& pi = ev.Get<QPulseInfo>("DAQ","PulseInfo");    
    const int chan = pi.GetChannelId();
    const QPulse& pulse = ev.Get<QPulse>("DAQ","Pulse");
    const QVector& samples = pulse.GetSamples();
 
    // get baseline value and rms computed by MBaselineModule in a previous sequence (diana_tutorial.cfg).
    const QBaselineData& baseData = ev.Get<QBaselineData>("BaselineModule","BaselineData");
    double baseline = baseData.GetBaseline();
    double baselineRMS = baseData.GetBaselineRMS();
 
    // get amplitude (Maximum minus baseline) computed in a previous sequence (diana_tutorial.cfg).
    const QDouble& amplitude = ev.Get<QDouble>("PulseBasicParameters","MaxBaseline");
 
    // get average pulse of this channel
    GlobalHandle<QVector> apHandle("Samples");
    apHandle.SetChannel(chan);
    GlobalData().Get("TutorialAveragePulse",&apHandle,fAveragePulseInput);
 
    if(!apHandle.IsValid()) {
        if(fChannelsWithoutAP.count(chan) == 0) {
            fChannelsWithoutAP.insert(chan);
            Warn("Average pulse not found for channel %04d",chan);
        }
        return;
    }
 
    const QVector& averagePulse = apHandle.Get();
 
    // compute amplitude and baseline of the average pulse
    double apBaseline = averagePulse[0];
    double apAmplitude = averagePulse.GetMax()-apBaseline;
 
    // estimate chi square
    double chisq = 0.;
    for(size_t i =0;  i < samples.Size(); i++) {
        double res = (samples[i]-baseline) - (averagePulse[i]-apBaseline)*amplitude/apAmplitude;
        chisq += res*res;
 
    }
    chisq /= (samples.Size() - 1)*baselineRMS*baselineRMS;
 
    // save chi square value of this event 
    ev.Get<QDouble>("ChiSquare") = chisq;
 
}
 
void MTutorialPulseShape::Done()
{
}