42 TH2D h2(
"h2", NULL, 1000, 0.0, 8.0, 1000, 0.0, 1000.0);
int main(int argc, char *argv[])
then usage $script< input file >[option[primary[working directory]]] nWhere option can be E
double getThetaMCS(const double E, const double x, const double X0, const double M, const double Q)
Get multiple Coulomb scattering angle.