TConstants.h

Go to the documentation of this file.

#ifndef RAD_TConstants
#define RAD_TConstants
 
#include <cmath>
#include <iostream>
#define M_ALPHA    0.00729735307964482 //  = 1/137.0359895
#define M_HC     197.327053
#define M_ME       0.51099907
#define M_MMU    105.658357
#define M_MTAU  1776.99
#define M_MPI    139.57018
#define M_MK0    497.672
#define M_MKC    493.677
#define M_META   547.75
class TConstants{
protected:
  double fAlphaScale;
  double fAlpha;
  double fAlpha2;
  double fAlpha3;
  double fAlphaPi;
public:
  TConstants(){
    fAlphaScale = 1;
    Init();
  };
  ~TConstants(){};
  void Init(){
    fAlpha = M_ALPHA*fAlphaScale;
    fAlpha2 = fAlpha*fAlpha;
    fAlpha3 = fAlpha2*fAlpha;
    fAlphaPi = fAlpha/M_PI;
  }
  void SetAlphaScale(const double &x){fAlphaScale = x; Init();}
  inline double GetAlphaScale(){return fAlphaScale;}
  inline double Alpha(){return fAlpha;}
  inline double Alpha2(){return fAlpha2;}
  inline double Alpha3(){return fAlpha3;}
  static double Pi(){return M_PI;}
  static double Pi2(){return M_PI*M_PI;}
  inline double AlphaPi(){return fAlphaPi;}
  static double HC(){return M_HC;}
  static double HC2(){return M_HC*M_HC;}
  static double Me(){return M_ME;}
  static double Me2(){return M_ME*M_ME;}
  static double Mmu(){return M_MMU;}
  static double Mmu2(){return M_MMU*M_MMU;}
  static double Mtau(){return M_MTAU;}
  static double Mtau2(){return M_MTAU*M_MTAU;}
  static double Mpi(){return M_MPI;}
  static double Mpi2(){return M_MPI*M_MPI;}
  static double Meta(){return M_META;}
  static double Meta2(){return M_META*M_META;}
  static double MK0(){return M_MK0;}
  static double MK02(){return M_MK0*M_MK0;}
  static double MKC(){return M_MKC;}
  static double MKC2(){return M_MKC*M_MKC;}
  void Print(){
    std::cout<<"Alpha scale      = "<<fAlphaScale<<std::endl;
  }
};
 
extern TConstants *gConst;
 
#endif //#ifndef RAD_TConstants