MPulseTimeConstants.cc

Go to the documentation of this file.

#include "MPulseTimeConstants.hh"
#include "QEvent.hh"
#include "QEventList.hh"
#include "QBaseType.hh"
#include "QInterval.hh"
#include "QRawEvent.hh"
#include "QBaselineData.hh"
#include "QPulseParameters.hh"
#include "QRunDataHandle.hh"
 
 
REGISTER_MODULE(MPulseTimeConstants)
 
using namespace Diana;
using namespace std;
 
void MPulseTimeConstants::Init(QEvent& ev)
{
    fTimeConstants = 0;
    fRiseVars.clear(); fDecayVars.clear();
 
    int nrise = GetInt("NumberOfRiseTimes", 1);
    if(nrise < 0) Panic("NumberOfRiseTimes must be >= 0");
    vector<QInterval> rises;
    for(int i = 1; i < nrise+1; i++) {
        char buf[256];
        char buf2[256];
        snprintf(buf,256,"RiseTime%d_Interval",i);
        vector<int> in; in.push_back(10); in.push_back(90);
        in = GetVectorInt(buf,in);
        if(in.size() != 2) Panic("%s needs 2 parameters instead of %d",buf,in.size());
        else if(in[0] > in[1]) Panic("%s: first element (%d) must be smaller than second element (%d)",buf,in[0],in[1]);
        else if(in[1] > 100) Panic("%s: elements cannot be greater than 100",buf);
        else if(in[0] < 0 ) Panic("%s: elements cannot be negative",buf);
 
        QInterval rise;
        rise.SetMin(double(in[0])/100); rise.SetMax(double(in[1])/100);
        rises.push_back(rise);
 
        snprintf(buf,256,"RiseTime%02d%02d",in[0],in[1]);
        snprintf(buf2,256,"RiseTime%02d%02dWindow",in[0],in[1]);
        ev.Add<QDouble>(buf);
        ev.Add<QDouble>(buf2);
        fRiseVars.push_back(buf);
    }
 
    int ndecay = GetInt("NumberOfDecayTimes", 1);
    if(ndecay < 0) Panic("NumberOfDecayTimes must be >= 0");
    vector<QInterval> decays;
    for(int i = 1; i < ndecay+1; i++) {
        char buf[256];
        char buf2[256];
        snprintf(buf,256,"DecayTime%d_Interval",i);
        vector<int> in; in.push_back(90); in.push_back(30);
        in = GetVectorInt(buf,in);
        if(in.size() != 2) Panic("%s needs 2 parameters instead of %d",buf,in.size());
        else if(in[1] > in[0]) Panic("%s: first element (%d) must be greater than second element (%d)",buf,in[0],in[1]);
        else if(in[0] > 100) Panic("%s: elements cannot be greater than 100",buf);
        else if(in[1] < 0 ) Panic("%s: elements cannot be negative",buf);
 
        QInterval decay;
        decay.SetMin(double(in[1])/100); decay.SetMax(double(in[0])/100);
        decays.push_back(decay);
 
        snprintf(buf,256,"DecayTime%02d%02d",in[0],in[1]);
        snprintf(buf2,256,"DecayTime%02d%02dWindow",in[0],in[1]);
        ev.Add<QDouble>(buf);
        ev.Add<QDouble>(buf2);
        fDecayVars.push_back(buf);
    }
 
    fTimeConstants = new QPulseTimeConstants(rises,decays);
 
    fPulseLabel = GetString("PulseLabel","DAQ@Pulse",false);
    fInputExtraLabel = GetString("InputExtraLabel","",false);
    fBaselineOwner = GetString("BaselineOwner","BaselineModule",false); 
    fPulseBasicParamsOwner = GetString("PulseBasicParametersOwner","PulseBasicParameters",false);
    if(!fInputExtraLabel.empty()) {
        fInputExtraLabel.insert(0,"_");
        fBaselineOwner += fInputExtraLabel;
        fPulseBasicParamsOwner += fInputExtraLabel;
    }
    fTomV = GetBool("ConvertTomV",true,false);
 
    ev.RequireByLabel<QPulse>(fPulseLabel);
    ev.Require<QBaselineData>(fBaselineOwner,"BaselineData");
    ev.Require<QPulseParameters>(fPulseBasicParamsOwner,"Parameters");
    ev.Require<QHeader>("DAQ","Header");
    ev.Require<QPulseInfo>("DAQ","PulseInfo");
 
}
 
void MPulseTimeConstants::Do(QEvent& ev)
{
    const QHeader& header = ev.Get<QHeader>("DAQ","Header");
    const int channel = ev.Get<QPulseInfo>("DAQ","PulseInfo").GetChannelId();
    const 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();
    const double baseline = ev.Get<QBaselineData>(fBaselineOwner.c_str(),"BaselineData").GetBaseline()/ADC2mV;
    int maxPos = ev.Get<QPulseParameters>(fPulseBasicParamsOwner.c_str(),"Parameters").fMaxPosition;
 
 
    vector<double> outputRise,outputDecay;
    QError err = fTimeConstants->Compute(samples,maxPos,baseline,outputRise,outputDecay);
    if(err != QERR_SUCCESS) return;
    
    if(outputRise.size() != fRiseVars.size()) Panic("Output rise time variables differ in number from expected");
    for(size_t i = 0; i < outputRise.size(); i++) {
        ev.Get<QDouble>(fRiseVars[i].c_str()) = outputRise[i]/samplingFreq;
        ev.Get<QDouble>((fRiseVars[i]+"Window").c_str()) = outputRise[i]/(double)samples.Size();
    }
 
    if(outputDecay.size() != fDecayVars.size()) Panic("Output decay time variables differ in number from expected");
    for(size_t i = 0; i < outputDecay.size(); i++) {
        ev.Get<QDouble>(fDecayVars[i].c_str()) = outputDecay[i]/samplingFreq;
        ev.Get<QDouble>((fDecayVars[i]+"Window").c_str()) = outputDecay[i]/(double)samples.Size();
    }
}
 
void MPulseTimeConstants::Done()
{
    if(fTimeConstants) delete fTimeConstants;
}