RdAntennaChannelToStationConverter.h

Go to the documentation of this file.

#ifndef _RdAntennaChannelToStationConverter_RdAntennaChannelToStationConverter_h_
#define _RdAntennaChannelToStationConverter_RdAntennaChannelToStationConverter_h_

#include <fwk/VModule.h>
#include <utl/Trace.h>
#include <utl/Vector.h>
#include <string>


namespace RdAntennaChannelToStationConverter {

  class RdAntennaChannelToStationConverter : public fwk::VModule {

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

  private:
    enum UsedDirection {
      eRdReconstruction,
      eSdReconstruction,
      eFdReconstruction,
      eMcTruth,
      eIndividualPerStation,
      eReference,
      eLineOfSightStationShowerMaximum
    };

    bool GetDirection(const evt::Event& event, utl::Vector& showerAxis, double& azimuth, double& zenith) const;

    double GetXmaxEstimator(const evt::Event& event) const;

    UsedDirection fUsedDirection = eReference;

    std::string fDirectionOfXmax;
    std::string fXmaxEstimator;

    bool fCropFreq = true;
    double fLowSmear = 0;
    double fUpSmear = 0;
    double fZeroLimit = 1e-5;

    bool fUseUserDefinedWeights = false;
    double fWeight12 = 1;
    double fWeight13 = 0;
    double fWeight23 = 0;

    bool fCalculateWeightsFromExpectedEField = true;
    bool fCalculateWeightsFromAntennaPattern = false;

    bool fRejectSaturatedStations = true;

    unsigned int fWhichEye = 6;

    // specify which implementation for the interpolation of the antenna pattern should be used
    std::string fInterpolationMode = "LinearInterpolation";

    utl::Vector3C rotateYaxis(utl::Vector3C vecC, double angle);

    utl::Vector3C rotateZaxis(utl::Vector3C vecC, double angle);

    REGISTER_MODULE("RdAntennaChannelToStationConverter", RdAntennaChannelToStationConverter);

  };

}


#endif