#ifndef PHYSICSCONSTANTS_H 
#define  PHYSICSCONSTANTS_H 

#include "TROOT.h"
#include "TMath.h"

// -*- C++ -*-

// 09/23/2008 Kei Moriya

// constants taken from 2006 PDG //
// 0. constants_______________________________
const Double_t CLIGHT = 29.9792458;

// 1. mass____________________________________

// we don't need this, but just for consistency
const double m_photon = 0;

// baryons
const Double_t m_p = 0.93827203;
const Double_t m_n = 0.93956536;
// hyperons
const Double_t m_Lambda     = 1.115683;
const Double_t m_SigmaPlus  = 1.18937;
const Double_t m_SigmaZero  = 1.192642;
const Double_t m_SigmaMinus = 1.197449;
const Double_t m_Sigma1385Plus  = 1.3828;
const Double_t m_Sigma1385Zero  = 1.3837;
const Double_t m_Sigma1385Minus = 1.3872;
const Double_t m_Lambda1405 = 1.406;
const Double_t m_Lambda1520 = 1.5195;

const Double_t m_XiZero = 1.31486;
const Double_t m_XiMinus = 1.32171;

const Double_t m_XiMinus1530 = 1.5350;
const Double_t m_XiZero1530 = 1.53180;

// mesons
const Double_t m_piZero = 0.1349766;
const Double_t m_piPlus = 0.13957018;
const Double_t m_KZero  = 0.497648;
const Double_t m_KPlus  = 0.493677;
const Double_t m_KstarPlus = 0.89166;
const Double_t m_KstarZero = 0.89600;

const Double_t m_eta = 0.547853;
const Double_t m_omega = 0.78265;

const Double_t m_phi = 1.019455;

// leptons
const Double_t m_electron = 0.00051099892;
const Double_t m_muon     = 0.105658369;

// 2. width___________________________________

// hyperons
const Double_t Gamma_Sigma1385Plus  = 0.0358;
const Double_t Gamma_Sigma1385Zero  = 0.036;
const Double_t Gamma_Sigma1385Minus = 0.0394;
const Double_t Gamma_Lambda1405 = 0.050;
const Double_t Gamma_Lambda1520 = 0.0156;

// mesons
const Double_t Gamma_KstarPlus = 0.0508;
const Double_t Gamma_KstarZero = 0.0503;

const Double_t Gamma_phi = 0.00426;

// 3. hyperon decay parameters_______________
const Double_t alpha_Lambda     = 0.642;
const Double_t alphaZero_SigmaPlus  = -0.980;
const Double_t alphaPlus_SigmaPlus  = 0.068;
const Double_t alpha_SigmaMinus = -0.068;

//_________________________________________________________________

// BES analysis stuff
const Double_t m_Jpsi     = 3.096916;
const Double_t m_psiPrime = 3.68609;

const Double_t m_chic0 = 3.41475;
const Double_t m_chic1 = 3.51066;
const Double_t m_chic2 = 3.55620;

const Double_t m_gamma = 0.;
const Double_t m_K0    = 0.497614;

const Double_t Gamma_Jpsi     = 92.9 * 1.e-6; // 92.9 keV
const Double_t Gamma_psiPrime = 304. * 1.e-6; // 304  keV

const Double_t Gamma_chic0 = 0.0104;  // 10.4 MeV
const Double_t Gamma_chic1 = 0.00086; // 0.86 MeV
const Double_t Gamma_chic2 = 0.00197; // 1.97 MeV

#endif