Arbeletche2021CherenkovAngularModel.h

Go to the documentation of this file.

#ifndef _atm_Arbeletche2021CherenkovAngularModel_h_
#define _atm_Arbeletche2021CherenkovAngularModel_h_

namespace atm {

  class Arbeletche2021CherenkovAngularModel {
  private:
    // fields computed once on construction
    double fOmegaEnergyDep, fEpsilon, fMaxTheta;
    // fields computed by UpdateParameters(s, n)
    mutable double fNu, fOmega1, fOmega2, fNormalization;
    mutable double fThetaEm, fCosThetaEm, fSinThetaEm;
  
    
    void UpdateParameters(const double showerAge, const double refractiveIndex) const;
    double FunctionK(const double theta) const;
    double UnnormalizedIntegral(const double lowAngle, const double highAngle) const;
    double UnnormalizedIntegralLeft(const double lowAngle, const double highAngle) const;
    double UnnormalizedIntegralRight(const double lowAngle, const double highAngle) const;
    double UnnormalizedIntegralSineLeft(const double lowAngle, const double highAngle) const;
    double UnnormalizedIntegralSineRight(const double lowAngle, const double highAngle) const;
  
  public:
    Arbeletche2021CherenkovAngularModel(const double showerEnergyTeV, const double maxTheta);

    double PDF(const double theta, const double showerAge, const double refractiveIndex) const;

    double CDF(const double theta, const double showerAge, const double refractiveIndex) const;

    double Integral(const double lowAngle, const double highAngle, const double showerAge, const double refractiveIndex) const;
  };

} // namespace atm

#endif /* _atm_Arbeletche2021CherenkovAngularModel_h_ */