MPulseBasicParameters.cc

Go to the documentation of this file.

#include "MPulseBasicParameters.hh"
#include "QPulseBasicParams.hh"
#include "QEvent.hh"
#include "QRawEvent.hh"
#include "QPulseParameters.hh"
#include "QBaselineData.hh"
#include "QBaseType.hh"
#include "QRunDataHandle.hh"
 
using namespace Diana;
 
REGISTER_MODULE(MPulseBasicParameters)
 
 
void MPulseBasicParameters::Init(QEvent& ev)
{
    fIsOnline = GetBool("IsOnline",false,false);
    ev.Add<QPulseParameters>("Parameters");
    ev.Add<QDouble>("MaxBaseline");
    ev.Add<QDouble>("MaxBaselineInRMS");
    ev.Add<QDouble>("MaxBaselineInWindow");
    ev.Add<QDouble>("MaxMinInWindow");
    ev.Add<QDouble>("MaxMinInWindowInRMS");
    ev.Add<QDouble>("RMS");
    ev.Add<QDouble>("Max");
    ev.Add<QDouble>("Min");
    ev.Add<QInt>("MaxPosInWindow");
    fPulseLabel = GetString("PulseLabel","DAQ@Pulse",false);
    fBaselineOwner = GetString("BaselineOwner","BaselineModule",false);
    fTomV = GetBool("ConvertTomV",true,false);
}
 
void MPulseBasicParameters::Do(QEvent& ev)
{
    double  MaxBaselineInWindow =0.;
    double  MaxMinInWindow =0.;
    int MaxPosInWindow = 0;
 
    const QVector &pulse = ev.GetByLabel<QPulse>(fPulseLabel).GetSamples();
    const int run = ev.Get<QHeader>("DAQ","Header").GetRun();
    const int channel = ev.Get<QPulseInfo>("DAQ","PulseInfo").GetChannelId();
 
    int triggerPosition = ev.Get<QPulseInfo>("DAQ","PulseInfo").GetMasterSample().GetSampleIndex();
    QPulseParameters &prep = ev.Get<QPulseParameters>("Parameters");
    QDouble& maxBaseline = ev.Get<QDouble>("MaxBaseline");
    const QBaselineData& bData = ev.Get<QBaselineData>(fBaselineOwner.c_str(),"BaselineData");
 
    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 int ADCMax = chanRunData.fADCMax;
    const int ADCMin = chanRunData.fADCMin;
    double ADC2Time = 1. / chanRunData.fSamplingFrequency;
    if(fIsOnline) {
        QPulseBasicParams bParams(triggerPosition,
                bData.GetBaseline()/ADC2mV, 
                                  bData.GetBaselineRMS()/ADC2mV, chanRunData.fSamplingFrequency);
        bParams.FindMaximumPosition(pulse);
        bParams.ComputeHeight(pulse);
        bParams.CheckSaturation(pulse,ADCMax,ADCMin);
        prep.fMaxPosition=bParams.GetMaxPosition();
        prep.fIsSaturatedHigh=bParams.GetIsSaturatedHigh();
        prep.fIsSaturatedLow=bParams.GetIsSaturatedLow();
 
        maxBaseline = bParams.GetMaxBaseline()*ADC2mV;
 
 
 
    } else {
        QPulseBasicParams bParams(pulse,triggerPosition,
                bData.GetBaseline()/ADC2mV, 
                                  bData.GetBaselineRMS()/ADC2mV, ADCMax, ADCMin, chanRunData.fSamplingFrequency);
 
        if(bParams.GetMaxBaseline() != Q_DOUBLE_DEFAULT)
            maxBaseline = bParams.GetMaxBaseline()*ADC2mV;
        prep.fMaxPosition = bParams.GetMaxPosition();
        prep.fIsSaturatedHigh = bParams.GetIsSaturatedHigh();
        prep.fIsSaturatedLow = bParams.GetIsSaturatedLow();
        prep.fStartPosition = bParams.GetStart();
        prep.fEndPosition = bParams.GetEnd();
        prep.fDecayTime = bParams.GetDecay()*ADC2Time;
        prep.fSlowDecayTime = bParams.GetSlowDecay()*ADC2Time;
        prep.fRiseTime = bParams.GetRaise()*ADC2Time;
 
    }
    double maxInWindow = ADC2mV*(pulse.GetMax());
    double minInWindow = ADC2mV*(pulse.GetMin());
    MaxPosInWindow = pulse.GetMaxIndex();
 
    MaxBaselineInWindow = maxInWindow-bData.GetBaseline();
    MaxMinInWindow = maxInWindow-minInWindow;
 
 
    if(maxBaseline != Q_DOUBLE_DEFAULT)
        ev.Get<QDouble>("MaxBaselineInRMS") = maxBaseline/bData.GetBaselineRMS();
    ev.Get<QDouble>("MaxBaselineInWindow") = MaxBaselineInWindow;
    ev.Get<QDouble>("MaxMinInWindow") = MaxMinInWindow;
    ev.Get<QDouble>("RMS") = ADC2mV*(pulse.GetRMS());
    ev.Get<QDouble>("Min") = (ADC2mV*pulse).GetMin();
    ev.Get<QDouble>("Max") =  (ADC2mV*pulse).GetMax();
    ev.Get<QDouble>("MaxMinInWindowInRMS") = MaxMinInWindow/bData.GetBaselineRMS();
    ev.Get<QInt>("MaxPosInWindow") = MaxPosInWindow;
 
}
 
void MPulseBasicParameters::Done()
{
        // called at the end of the event loop
}