RdStationSimulationRejector.h

Go to the documentation of this file.

#ifndef _RdStationSimulationRejector_RdStationSimulationRejector_h_
#define _RdStationSimulationRejector_RdStationSimulationRejector_h_


#include <fwk/VModule.h>

#include <cmath>


namespace RdStationSimulationRejector {

  class RdStationSimulationRejector : public fwk::VModule {

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

  private:
    // parameterisation at which radial axis distance thinning might become a problem.
    // depends on your thinning level, cf. PhD thesis of Felix Schlüter
    double thinning_distance(const double zenith)
    { return (1.93259388e-02 * exp(zenith * 7.89363304e+00) + 6.37926253e+02) * 0.85; }

    bool fRejectStationAffectedByThinning = false;
    bool fRemoveNonTriggeredRdStationsFromEvent = true;

    REGISTER_MODULE("RdStationSimulationRejector", RdStationSimulationRejector);
  };

}

#endif