RdStationEFieldVectorCalculator.h

Go to the documentation of this file.

#ifndef _RdStationEFieldVectorCalculator_RdStationEFieldVectorCalculator_h_
#define _RdStationEFieldVectorCalculator_RdStationEFieldVectorCalculator_h_

#include <fwk/VModule.h>
#include <revt/Station.h>
#include <TVector3.h>


namespace RdStationEFieldVectorCalculator {

  class RdStationEFieldVectorCalculator : public fwk::VModule {

  public:
    fwk::VModule::ResultFlag Init() override;
    fwk::VModule::ResultFlag Run(evt::Event& event) override;
    fwk::VModule::ResultFlag Finish() override;

  private:
    // Peakposistion from what trace? 0=EW, 1=NS, 2=V, 3=Trace with max. amp
    unsigned int fTrace = 0;

    // where to search for the maximum peak: -1/+1 = before/after max. of evenlope of mag, 0 = everywhere
    int fPeakSearch = 0;

    // search window for maximum puls (only has affect if fPeakSearch != 0)
    double fPeakSearchWindowWidth = 0;

    unsigned int fPeakSearchWindow = 0;
    bool fMeanEFieldVector = false;
    double fChargeExcessStrength = 0;

    double GetAngleToLorentzForceError(double snr);
    std::pair<double, double> GetFWHM(const revt::StationTimeSeries& stationEnvelope, const double maxSignal,
      const double tPeakPos);
    TVector3 GetMeanEfieldInRange(const revt::StationTimeSeries& stationTrace,
      const std::pair<double, double>& range, const TVector3& firstGuess);
    double NormalizeAngle(const double x);

    REGISTER_MODULE("RdStationEFieldVectorCalculator", RdStationEFieldVectorCalculator);
  };

}

#endif