RdChannelRiseTimeCalculator.h

Go to the documentation of this file.

#ifndef _RdChannelRiseTimeCalculator_h_
#define _RdChannelRiseTimeCalculator_h_

#include <fwk/VModule.h>
#include <revt/Channel.h>
#include <vector>


namespace evt {
  class Event;
}

namespace revt {
  class REvent;
  class Channel;
}

namespace RdChannelRiseTimeCalculator {

  class RdChannelRiseTimeCalculator : public fwk::VModule {

  public:
    RdChannelRiseTimeCalculator() { }
    virtual ~RdChannelRiseTimeCalculator() { }

    fwk::VModule::ResultFlag Init() override;
    fwk::VModule::ResultFlag Run(evt::Event& event) override;
    fwk::VModule::ResultFlag Finish() override { return eSuccess; }

  private:
    unsigned int fRiseTimeWindow = 0;

    // calculates Mean and RMS in selected noise window
    void Noisefinder(const revt::ChannelTimeSeries& channeltrace,
                     double& noiseRMS, const double noiseWindowStart, const double noiseWindowStop) const;

    // calculates the amplitude and the position of the maximum in the signal search window
    void Pulsefinder(const revt::ChannelTimeSeries& channeltrace,
                     double& peakAmplitude, double& peakTime, double& peakTimeError,
                     const double signalSearchWindowStart, const double signalSearchWindowStop,
                     unsigned int& sample) const;

    // Construction of the cumulative function in the rise time window
    void MakeCumulative(const revt::ChannelTimeSeries& channeltrace,
                        const unsigned int startRiseTimeWindow, const unsigned int stopRiseTimeWindow,
                        std::vector<double>& cumulative, double& cumulativeSum) const;

    // Rise time calculation
    void GetRiseTime(const unsigned int startCumulative, const unsigned int stopCumulative,
                     std::vector<double>& normalizedCumulative, double &riseTime, const double samplingFrequency) const;

    REGISTER_MODULE("RdChannelRiseTimeCalculator", RdChannelRiseTimeCalculator);

  };

}


#endif