RdDirectionConvergenceChecker.h

Go to the documentation of this file.

#ifndef _RdDirectionConvergenceChecker_RdDirectionConvergenceChecker_h_
#define _RdDirectionConvergenceChecker_RdDirectionConvergenceChecker_h_

#include <fwk/VModule.h>

#include <utl/AugerUnits.h>

namespace RdDirectionConvergenceChecker {

  class RdDirectionConvergenceChecker : public fwk::VModule {

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

  private:
    double fAngleEpsilon = 0.5 * utl::degree;

    unsigned long int fMaxIterations = 10;

    bool fCheckMirrorZenith = false;

    // vectors to store the azimuth and zenith angles of all iterations for an event
    std::vector<double> fAllAzimuths;
    std::vector<double> fAllZeniths;

    double GetAngularDifference(const double zen1, const double azim1, const double zen2, const double azim2);

    REGISTER_MODULE("RdDirectionConvergenceChecker", RdDirectionConvergenceChecker);
  };

}


#endif