40{
44
47 bool logx;
49 std::string option;
51
52
53 try {
54
55 JParser<> zap(
"Example program to histogram neutrino effective mass for triggered events.");
56
59 zap[
'X'] =
make_field(logx,
"Use logarithm of energy");
61 zap[
'O'] =
make_field(option,
"Result option") = Mass_t, Volume_t;
63
64 zap(argc, argv);
65 }
66 catch(const exception &error) {
67 FATAL(error.what() << endl);
68 }
69
70 try {
71
72 double Xmin = numeric_limits<double>::max();
73 double Xmax = numeric_limits<double>::lowest();
74
76
78
79 const JVolume volume(scanner->getHeader(), logx);
80
81 if (volume.getXmin() < Xmin) { Xmin = volume.getXmin(); }
82 if (volume.getXmax() > Xmax) { Xmax = volume.getXmax(); }
83 }
84
86
88
90
91 NOTICE(
"Scanning file type " << scanner->getName() << endl);
92 DEBUG (header << endl);
93
94 const JEvtWeight& weighter = scanner->getEvtWeighter();
95
97 FATAL(
"Mising normalisation in header." << endl);
98 }
99
101
104 else
106
107 while (scanner->hasNext()) {
108
109 STATUS(
"event: " << setw(10) << scanner->getCounter() <<
'\r');
DEBUG(endl);
110
111 const Evt*
event = scanner->next();
113
116 double w = (
getVolume(weighter, *event) /
117 (logx ? log(10.0) *
primary.E * hm->GetBinWidth(1) : hm->GetBinWidth(1)) /
118 Wnorm);
119
120 if (option == Mass_t) {
122 } else if (option == Volume_t) {
124 }
125
127 }
129 }
130
132
133 out << hm;
134
135 out.Close();
136 }
138 FATAL(error << endl);
139 }
140}
#define DEBUG(A)
Message macros.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
int numberOfBins
number of bins for average CDF integral of optical module
const JHead & getHeader() const
Get header.
bool is_valid(T JHead::*pd) const
Check validity of given data member in JHead.
Utility class to parse command line options.
Auxiliary class to manage set of compatible ROOT objects (e.g. histograms) using unique keys.
double getVolume(const JType< JEvtWeightGSeaGen > &type, const Evt &evt)
Get volume of given event according given weighter.
const Trk & get_primary(const Evt &evt)
Auxiliary function to retrieve the primary track of an event.
static const double PI
Mathematical constants.
static const double DENSITY_SEA_WATER
Fixed environment values.
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
KM3NeT DAQ data structures and auxiliaries.
The Evt class respresent a Monte Carlo (MC) event as well as an offline event.
Abstract base class for event weighing.
Auxiliary class for histogramming of effective volume.
JRange_t cosT
Cosine zenith angle range
double numberOfEvents
Number of events.
Auxiliary class for organising Monte Carlo file scanners associated with event weighters.
std::vector< filescanner_type >::const_iterator const_iterator
std::vector< filescanner_type >::iterator iterator
Auxiliary base class for list of file names.
The Trk class represents a Monte Carlo (MC) particle as well as a reconstructed track/shower.