1 #ifndef __JPHYSICS_JPHYSICSTOOLKIT__
2 #define __JPHYSICS_JPHYSICSTOOLKIT__
12 namespace JPP {
using namespace JPHYSICS; }
42 const double p = sqrt((E + M) * (E - M));
43 const double beta = p / E;
45 return THETA_MCS * Q * sqrt(x/X0) * (1.0 + 0.038*log(x/X0)) / (beta*p);
61 using namespace JTOOLS;
double getThetaMCS(const double E, const double x, const double X0, const double M, const double Q)
Get multiple Coulomb scattering angle.