34 int main(
int argc,
char **argv)
41 JLimit_t& numberOfEvents = inputFile.getLimit();
50 JParser<> zap(
"Auxiliary program to manipulate MUPAGE data.");
54 zap[
'n'] =
make_field(numberOfEvents) = JLimit::max();
55 zap[
'P'] =
make_field(P,
"Keep muons in event with given probability (<= 1, 1 == all)") = 1.0;
56 zap[
'u'] =
make_field(
u,
"Reject events as a function of multiplicity (>= 0, 0 == none)") = 0.0;
57 zap[
'E'] =
make_field(E,
"Energy scaling factor") = 1.0;
63 catch(
const exception &error) {
64 FATAL(error.what() << endl);
67 gRandom->SetSeed(seed);
93 if (gRandom->Rndm() < P)
96 i =
event->mc_trks.erase(i);
109 i =
event->mc_trks.erase(i);
114 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.
Template specialisation of JMultipleFileScanner for Monte Carlo header.
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 double MASS_MUON
muon mass [GeV]
static const JGeaneWater gWater
Function object for energy loss of muon in sea water.
const double getSpeedOfLight()
Get speed of light.
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.
Auxiliary class for defining the range of iterations of objects.