20int main(
int argc,
char* argv[])
29 JParser<> zap(
"Example program to histogram photon emission profile from showers.");
36 catch(
const exception &error) {
37 FATAL(error.what() << endl);
45 TH1D h1(
"Probability", NULL, 1000, -1.0, +1.0);
46 TH1D h2(
"Integral", NULL, 1000, -1.0, +1.0);
48 for(
int i = 1; i <= h1.GetNbinsX(); ++i) {
50 const double x = h1.GetBinCenter(i);
52 h1.SetBinContent(i, geanx(x));
53 h2.SetBinContent(i, geanx(-1.0,x));