18 int main(
int argc,
char **argv)
35 catch(
const exception &error) {
36 FATAL(error.what() << endl);
42 TH2D h2(
"h2", NULL, 1000, 0.0, 8.0, 1000, 0.0, 1000.0);
44 for (
int ix = 1; ix <= h2.GetXaxis()->GetNbins(); ++ix) {
45 for (
int iy = 1; iy <= h2.GetYaxis()->GetNbins(); ++iy) {
47 const double x = h2.GetXaxis()->GetBinCenter(ix);
48 const double y = h2.GetYaxis()->GetBinCenter(iy);
50 const double E =
pow(10.0, x);
56 h2.SetBinContent(ix, iy, z);
Utility class to parse command line options.
int main(int argc, char *argv[])
static const JGeaneWater gWater
Function object for energy loss of muon in sea water.
double getThetaMCS(const double E, const double x, const double X0, const double M, const double Q)
Get multiple Coulomb scattering angle.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
T pow(const T &x, const double y)
Power .
General purpose messaging.
Utility class to parse command line options.