AtmCorr.h

Go to the documentation of this file.

#ifndef __AtmCorr_h_
#define __AtmCorr_h_

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <fstream>

//----------------------------------------------------------------------------------------------
//   3 types of atmospheric effects at the same (r,DX) :
//          Type I --> Showers with the same (theta,psi) but different rho_ground
//          Type II--> Showers with the same (psi,rho_ground)  show theta depedences.
//          Type III--> Showers with the same (theta,rho_ground) show psi depedences due to the different atmospheric profiles
//----------------------------------------------------------------------------------------------

// Remarks :
//    -A wrong attenuation parameter can also induce azimuthal asymmetries that can be confused with atmospheric effects
//    -Type III corrections for AMIGA are due to imperfections in the Geometrical asymmetries

// Corrections calculated using FeQGSJetII


namespace AtmCorrNS {

  const double parAtmTypeI[2][4] = { {  0.0000 , 0.0735 , 0.0395 , 0.0227 } , { 1.0000 , 0.4797 , 0.7292 , 1.3356 }};

  const double f_mod1_par_WCD[4][3] = { { 0.000, -0.013, 0.089 } , { 0.000, 0.278, 0.028 } , { 0.000, 0.080, 0.052 } , { 0.000, -0.027, 0.083 }};
  const double f_mod2_par_WCD[4][3] = { { 0.000, -0.108, -0.040 } , { 0.000, -0.436, 0.038 } , { 0.000, -0.182, -0.007 } , { 0.000, -0.120, -0.041 }};

  const double f_mod1_par_Scin[4][3] = { { 0.000, -0.021, 0.064 } , { 0.000, 0.127, 0.078 } , { 0.000, -0.014, 0.067 } , { 0.000, -0.053, 0.069 }};
  const double f_mod2_par_Scin[4][3] = { { 0.000, -0.085, -0.018 } , { 0.000, -0.379, 0.047 } , { 0.000, -0.118, -0.019 } , { 0.000, -0.167, 0.015 }};

  const double f_mod1_par_MD[4][3] = { { 0.000, -0.083, 0.178 } , {0.000, 0.000, 0.000}, {0.000, 0.000, 0.000}, {0.000, 0.000, 0.000} };
  const double f_mod2_par_MD[4][3] = { { 0.000, 0.000, -0.072 } , {0.000, 0.000, 0.000}, {0.000, 0.000, 0.000}, {0.000, 0.000, 0.000} };

  class AtmCorr {

  private:
    int fDetectorType;
  public:
    AtmCorr(int DetectorType);
    ~AtmCorr() {};

    //Type I atmospheric effect
    double GetAlphaRho(double r, int icomp);
    double GetCorrTypeI(double r, int icomp, double rhoGround, double hGround);

    //Type II atmospheric effects
    double GetCorrTypeII(double r, double theta, int icomp, double hground);

    //Type III atmospheric effect
    double GetParCorrTypeIII(double r,  int icomp, int ipar);
    double GetCorrTypeIII(double r, double theta, double psi, int icomp);

  };
}
#endif // __AtmCorr_h_