#ifndef physics_constants #define physics_constants const double alphaQED=1/137.036; const double hbarcSqr=389.37966; // in GeV*GeV*ub const double mElectron=0.51099907e-3; const double mMuon=0.105658389; const complex I_(0,1); const double PI_=2*atan2(1,0);