CherenkovModel.h

Go to the documentation of this file.

#ifndef _atm_CherenkovModel_h_
#define _atm_CherenkovModel_h_


namespace atm {

  class CherenkovModel {

  public:
    //Global Constants.
    static const int kN_WAVE_BIN = 20;
    static const double kWAVE_MAX = 449.3;
    static const double kWAVE_MIN = 303.1;
    static const double kWAVE_STEP = (kWAVE_MAX - kWAVE_MIN) / kN_WAVE_BIN;

  public:
    CherenkovModel();

    double GaisserHillas(double chi,double chi0,double chimax, double Nmax,double lambda) const;
    double Energy(double beta) const;
    double Beta(double energy) const;
    double CKVThreshold(double n_index) const;
    double Distance(double x1[],double x2[]) const;
    double GDistance(double x1[],double x2[]) const;
    double TrackLengthH(double energy,double chi,double chimax) const;
    double TrackLengthN(double energy,double chi,double chimax) const;
    double CKVPhot(double trl,double n_index,double beta,double dslant,double l1,double l2) const;
    void CKVIntegral(double& xlo,double& xhi,double& xb,double& integral) const;
    void CKVDirect(double bin,double theta,double height,double rp,double chimax,std::vector<double>& DirCh)const;
    void CKVBeam(double h1,double h2,double theta,double chimax,std::vector<double>& BeamCh) const;
    void CKVRayleigh(double h1,double h2,double theta,double rr,std::vector<double>& RayCh) const;
    void CKVAerosol(double h1,double h2,double theta,double rr,std::vector<double>& AerCh) const;

    virtual ~CherenkovModel() { }

  private:
    double Fu(double x,double a) const;

  };

}


#endif