#ifndef MUON_H #define MUON_H #include /* for size_t */ #define EULER_CONSTANT 0.57721 double muon_get_range(double T, double rho); double muon_get_dEdx_rad(double T, double rho); double muon_get_dEdx(double T, double rho); #endif