RdStationRiseTimeCalculator.h

Go to the documentation of this file.

#ifndef _RdStationRiseTimeCalculator_h_
#define _RdStationRiseTimeCalculator_h_

#include <fwk/VModule.h>
#include <io/IoCodes.h>  // needed for "Status"
#include <revt/Channel.h>
#include <revt/Station.h>
#include <vector>


namespace evt {
  class Event;
}

namespace revt {
  class REvent;
  class Station;
}

namespace RdStationRiseTimeCalculator {

  class RdStationRiseTimeCalculator : public fwk::VModule {

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

    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;

    // Construction of the cumulative function in the rise time window
    void MakeCumulative(const revt::StationTimeSeries& Stationtrace,
                        const unsigned int startRiseTimeWindow, const unsigned int stopRiseTimeWindow,
                        std::vector<double>& cumulative, double& cumulativeSum, unsigned int chwithmaxSNR) const;
    // Rise time calculation
    void GetRiseTime(const unsigned int startCumulative, const unsigned int stopCumulative,
                     std::vector<double>& normalizedCumulative, unsigned int& minCumulative, unsigned int& maxCumulative) const;

    REGISTER_MODULE("RdStationRiseTimeCalculator", RdStationRiseTimeCalculator);

  };

}


#endif