Jpp  master_rocky-40-g5f0272dcd
the software that should make you happy
JStingray.cc
Go to the documentation of this file.
1 
2 #include <string>
3 #include <iostream>
4 #include <iomanip>
5 #include <map>
6 
7 #include "TRandom3.h"
8 
9 #include "km3net-dataformat/offline/Head.hh"
10 #include "km3net-dataformat/offline/MultiHead.hh"
11 #include "km3net-dataformat/offline/Evt.hh"
12 
13 #include "JSupport/JMultipleFileScanner.hh"
14 #include "JSupport/JMonteCarloFileSupportkit.hh"
15 #include "JSupport/JFileRecorder.hh"
16 #include "JSupport/JSupport.hh"
17 #include "JSupport/JMeta.hh"
18 
19 #include "JPhysics/JGeane.hh"
20 #include "JPhysics/JGeanz.hh"
21 #include "JPhysics/JConstants.hh"
22 
23 #include "Jeep/JPrint.hh"
24 #include "Jeep/JParser.hh"
25 #include "Jeep/JMessage.hh"
26 
27 
28 /**
29  * \file
30  *
31  * Auxiliary program to manipulate MUPAGE data.
32  * \author mdejong
33  */
34 int main(int argc, char **argv)
35 {
36  using namespace std;
37  using namespace JPP;
38 
39  JMultipleFileScanner<Evt> inputFile;
40  JFileRecorder<JTYPELIST<JAAnetTypes_t, JMeta, JRootTypes_t>::typelist> outputFile;
41  JLimit_t& numberOfEvents = inputFile.getLimit();
42  double P;
43  double u;
44  double E;
45  UInt_t seed;
46  int debug;
47 
48  try {
49 
50  JParser<> zap("Auxiliary program to manipulate MUPAGE data.");
51 
52  zap['f'] = make_field(inputFile);
53  zap['o'] = make_field(outputFile) = "mupage.root";
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;
58  zap['S'] = make_field(seed) = 0;
59  zap['d'] = make_field(debug) = 2;
60 
61  zap(argc, argv);
62  }
63  catch(const exception &error) {
64  FATAL(error.what() << endl);
65  }
66 
67  gRandom->SetSeed(seed);
68 
69 
70  outputFile.open();
71 
72  if (!outputFile.is_open()) {
73  FATAL("Error opening file " << outputFile << endl);
74  }
75 
76  outputFile.put(JMeta(argc, argv));
77  outputFile.put(*gRandom);
78 
79  JMultipleFileScanner<Head> io(inputFile);
80 
81  io >> outputFile;
82 
83  while (inputFile.hasNext()) {
84 
85  STATUS("event: " << setw(10) << inputFile.getCounter() << '\r'); DEBUG(endl);
86 
87  Evt* event = inputFile.next();
88 
89  // randomly reject muons in event
90 
91  for (vector<Trk>::iterator i = event->mc_trks.begin(); i != event->mc_trks.end(); ) {
92 
93  if (gRandom->Rndm() < P)
94  ++i;
95  else
96  i = event->mc_trks.erase(i);
97  }
98 
99  // apply energy loss
100 
101  for (vector<Trk>::iterator i = event->mc_trks.begin(); i != event->mc_trks.end(); ) {
102 
103  i->E = gWater(i->E, i->t * getSpeedOfLight());
104  i->E *= E;
105 
106  if (i->E > MASS_MUON)
107  ++i;
108  else
109  i = event->mc_trks.erase(i);
110  }
111 
112  // randonly reject events as a function of muon multiplicity
113 
114  if (gRandom->Rndm() < pow(event->mc_trks.size(),-u)) {
115  outputFile.put(*event);
116  }
117  }
118  STATUS(endl);
119 
120  outputFile.put(*gRandom);
121  outputFile.close();
122 }
outputFile
string outputFile
Definition: JDAQTimesliceSelector.cc:37
JFileRecorder.hh
Recording of objects on file according a format that follows from the file name extension.
JGeane.hh
Energy loss of muon.
JGeanz.hh
Longitudinal emission profile EM-shower.
JMessage.hh
General purpose messaging.
DEBUG
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62
STATUS
#define STATUS(A)
Definition: JMessage.hh:63
FATAL
#define FATAL(A)
Definition: JMessage.hh:67
debug
int debug
debug level
Definition: JSirene.cc:69
JMeta.hh
ROOT I/O of application specific meta data.
JMultipleFileScanner.hh
Scanning of objects from multiple files according a format that follows from the extension of each fi...
JParser.hh
Utility class to parse command line options.
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition: JParser.hh:2142
JConstants.hh
Physics constants.
JPrint.hh
I/O formatting auxiliaries.
main
int main(int argc, char **argv)
Definition: JStingray.cc:34
JSupport.hh
ROOT TTree parameter settings of various packages.
JPARSER::JParser
Utility class to parse command line options.
Definition: JParser.hh:1698
JSUPPORT::JFileRecorder
Object writing to file.
Definition: JFileRecorder.hh:44
JSUPPORT::JMultipleFileScanner< Head >
Template specialisation of JMultipleFileScanner for Monte Carlo header.
Definition: JMonteCarloFileSupportkit.hh:416
JSUPPORT::JMultipleFileScanner::hasNext
virtual bool hasNext() override
Check availability of next element.
Definition: JMultipleFileScanner.hh:409
JSUPPORT::JMultipleFileScanner::getCounter
counter_type getCounter() const
Get counter.
Definition: JMultipleFileScanner.hh:382
JSUPPORT::JMultipleFileScanner::next
virtual const pointer_type & next() override
Get next element.
Definition: JMultipleFileScanner.hh:457
JMATH::pow
T pow(const T &x, const double y)
Power .
Definition: JMath.hh:97
JPHYSICS::MASS_MUON
static const double MASS_MUON
muon mass [GeV]
Definition: JPhysics/JConstants.hh:69
JPHYSICS::gWater
static const JGeaneWater gWater
Function object for energy loss of muon in sea water.
Definition: JGeane.hh:381
JPHYSICS::getSpeedOfLight
const double getSpeedOfLight()
Get speed of light.
Definition: JPhysics/JConstants.hh:112
JPP
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Definition: JAAnetToolkit.hh:43
JTOOLS::u
double u[N+1]
Definition: JPolint.hh:865
std
Definition: JSTDTypes.hh:14
Evt
The Evt class respresent a Monte Carlo (MC) event as well as an offline event.
Definition: Evt.hh:21
JSUPPORT::JLimit
Auxiliary class for defining the range of iterations of objects.
Definition: JLimit.hh:45
JSUPPORT::JMeta
Auxiliary class for ROOT I/O of application specific meta data.
Definition: JMeta.hh:72
Generated by doxygen 1.9.1