RdHASLDFFitter.h

Go to the documentation of this file.

#ifndef _RdHASLDFFitter_RdHASLDFFitter_h_
#define _RdHASLDFFitter_RdHASLDFFitter_h_

#include <fwk/VModule.h>
#include <utl/PhysicalFunctions.h>

#include <string>


namespace RdHASLDFFitter {

  class RdHASLDFFitter : public fwk::VModule {

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

  private:
    // Parameter for energy calibration Srad -> Eem
    double fS19 = 0;
    double fGamma = 0;

    // Parameter for density correction Erad -> Srad
    double fAvgDensity = 0;
    double fp0 = 0;
    double fp1 = 0;

    // correction factor to S19
    double fSTEPFCCorrectionToSimulatedRadioEnergyScale = 0;

    unsigned int fMinimumNumberOfStations = 0;

    bool fFitDistanceToXmax = true;
    bool fFitCore = true;

    bool fAddLowestSignalAsError = false;
    double fAddRelativeMaxSignalError = 0;

    std::string fXmaxType;
    utl::HadronicInteractionModel fHadronicInteractionModel =
      utl::HadronicInteractionModel::eQGSJETII04;
    std::string fFitDirection;

    bool fUseSaturatedStations = false;
    bool fUseParametrizationToDisentanglePolarisation = true;
    bool fUseSoftExponent = true;

    bool fApplyParametricCorrectionToElectromagneticEnergyUncertainty = true;

    int fSuccessfulFits = 0;
    int fFailedFits = 0;

    bool GetXmaxEstimator(const evt::Event& event, double& xmax) const;

    void ElectromagneticEnergy(const double sRad, const double sRadErr,
      double& electromagneticEnergyRec, double& electromagneticEnergyRecErr, const double zenith) const;

    void RadiationDensityCorrection(const double densityMax, const double densityMaxErr,
      double& radiationDensityCorrection, double& radiationDensityCorrectionErr) const;

    REGISTER_MODULE("RdHASLDFFitter", RdHASLDFFitter);

  };

}


#endif