36int main(
int argc,
char **argv)
43 JLimit_t& numberOfEvents = inputFile.getLimit();
52 JParser<> zap(
"Auxiliary program to manipulate MUPAGE data.");
57 zap[
'P'] =
make_field(P,
"Keep muons in event with given probability (<= 1, 1 == all)") = 1.0;
58 zap[
'u'] =
make_field(u,
"Reject events as a function of multiplicity (>= 0, 0 == none)") = 0.0;
59 zap[
'E'] =
make_field(E,
"Energy scaling factor") = 1.0;
65 catch(
const exception &error) {
66 FATAL(error.what() << endl);
94 if (gRandom->Rndm() < P)
97 i =
event->mc_trks.erase(i);
104 i->E =
gWater(i->E, i->t * getSpeedOfLight());
107 if (i->E > MASS_MUON)
110 i =
event->mc_trks.erase(i);
115 if (gRandom->Rndm() <
pow(event->mc_trks.size(),-u)) {
Recording of objects on file according a format that follows from the file name extension.
Longitudinal emission profile EM-shower.
General purpose messaging.
#define DEBUG(A)
Message macros.
Scanning of objects from multiple files according a format that follows from the extension of each fi...
Utility class to parse command line options.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
I/O formatting auxiliaries.
int main(int argc, char **argv)
ROOT TTree parameter settings of various packages.
Utility class to parse command line options.
General purpose class for object reading from a list of file names.
virtual bool hasNext() override
Check availability of next element.
counter_type getCounter() const
Get counter.
virtual const pointer_type & next() override
Get next element.
T pow(const T &x, const double y)
Power .
static const JGeaneWater gWater
Function object for energy loss of muon in sea water.
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
The Evt class respresent a Monte Carlo (MC) event as well as an offline event.
Template definition of random value generator.
Auxiliary class for defining the range of iterations of objects.
static counter_type max()
Get maximum counter value.