Jpp
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
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/Evt.hh"
11 
12 #include "JSupport/JMultipleFileScanner.hh"
13 #include "JSupport/JMonteCarloFileSupportkit.hh"
14 #include "JSupport/JFileRecorder.hh"
15 #include "JSupport/JSupport.hh"
16 #include "JSupport/JMeta.hh"
17 
18 #include "JPhysics/JGeane.hh"
19 #include "JPhysics/JGeanz.hh"
20 #include "JTools/JConstants.hh"
21 
22 #include "Jeep/JPrint.hh"
23 #include "Jeep/JParser.hh"
24 #include "Jeep/JMessage.hh"
25 
26 
27 /**
28  * \file
29  *
30  * Auxiliary program to manipulate MUPAGE data.
31  * \author mdejong
32  */
33 int main(int argc, char **argv)
34 {
35  using namespace std;
36  using namespace JPP;
37 
38  JMultipleFileScanner<Evt> inputFile;
39  JFileRecorder<JTYPELIST<JAAnetTypes_t, JMeta, JRootTypes_t>::typelist> outputFile;
40  JLimit_t& numberOfEvents = inputFile.getLimit();
41  double P;
42  double u;
43  double E;
44  UInt_t seed;
45  int debug;
46 
47  try {
48 
49  JParser<> zap("Auxiliary program to manipulate MUPAGE data.");
50 
51  zap['f'] = make_field(inputFile);
52  zap['o'] = make_field(outputFile) = "mupage.root";
53  zap['n'] = make_field(numberOfEvents) = JLimit::max();
54  zap['P'] = make_field(P, "Keep muons in event with given probability (<= 1, 1 == all)") = 1.0;
55  zap['u'] = make_field(u, "Reject events as a function of multiplicity (>= 0, 0 == none)") = 0.0;
56  zap['E'] = make_field(E, "Energy scaling factor") = 1.0;
57  zap['S'] = make_field(seed) = 0;
58  zap['d'] = make_field(debug) = 2;
59 
60  zap(argc, argv);
61  }
62  catch(const exception &error) {
63  FATAL(error.what() << endl);
64  }
65 
66  gRandom->SetSeed(seed);
67 
68 
69  outputFile.open();
70 
71  if (!outputFile.is_open()) {
72  FATAL("Error opening file " << outputFile << endl);
73  }
74 
75  outputFile.put(JMeta(argc, argv));
76  outputFile.put(*gRandom);
77 
78  JMultipleFileScanner<Head> io(inputFile);
79 
80  io >> outputFile;
81 
82  while (inputFile.hasNext()) {
83 
84  STATUS("event: " << setw(10) << inputFile.getCounter() << '\r'); DEBUG(endl);
85 
86  Evt* event = inputFile.next();
87 
88  // randomly reject muons in event
89 
90  for (vector<Trk>::iterator i = event->mc_trks.begin(); i != event->mc_trks.end(); ) {
91 
92  if (gRandom->Rndm() < P)
93  ++i;
94  else
95  i = event->mc_trks.erase(i);
96  }
97 
98  // apply energy loss
99 
100  for (vector<Trk>::iterator i = event->mc_trks.begin(); i != event->mc_trks.end(); ) {
101 
102  i->E = gWater(i->E, i->t * getSpeedOfLight());
103  i->E *= E;
104 
105  if (i->E > MASS_MUON)
106  ++i;
107  else
108  i = event->mc_trks.erase(i);
109  }
110 
111  // randonly reject events as a function of muon multiplicity
112 
113  if (gRandom->Rndm() < pow(event->mc_trks.size(),-u)) {
114  outputFile.put(*event);
115  }
116  }
117  STATUS(endl);
118 
119  outputFile.put(*gRandom);
120  outputFile.close();
121 }
JSUPPORT::JMeta
Auxiliary class for ROOT I/O of application specific meta data.
Definition: JMeta.hh:71
JSUPPORT::JFileRecorder
Object writing to file.
Definition: JFileRecorder.hh:41
JPARSER::JParser
Utility class to parse command line options.
Definition: JParser.hh:1493
JSupport.hh
ROOT TTree parameter settings.
JTOOLS::MASS_MUON
static const double MASS_MUON
muon mass [GeV]
Definition: JConstants.hh:59
JTOOLS::getSpeedOfLight
const double getSpeedOfLight()
Number of bytes in a gigabyte.
Definition: JConstants.hh:89
JGeane.hh
Energy loss of muon.
STATUS
#define STATUS(A)
Definition: JMessage.hh:63
JFileRecorder.hh
Recording of objects on file according a format that follows from the file name extension.
JPHYSICS::gWater
static const JGeaneWater gWater
Function object for energy loss of muon in sea water.
Definition: JGeane.hh:328
outputFile
string outputFile
Definition: JDAQTimesliceSelector.cc:37
JSUPPORT::JLimit
Auxiliary class for defining the range of iterations of objects.
Definition: JLimit.hh:41
JConstants.hh
Constants.
JPrint.hh
I/O formatting auxiliaries.
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition: JParser.hh:1954
JMeta.hh
ROOT I/O of application specific meta data.
debug
int debug
debug level
Definition: JSirene.cc:61
JMessage.hh
General purpose messaging.
JSUPPORT::JMultipleFileScanner< Head >
Template specialisation of JMultipleFileScanner for Monte Carlo header.
Definition: JMonteCarloFileSupportkit.hh:324
FATAL
#define FATAL(A)
Definition: JMessage.hh:67
JMultipleFileScanner.hh
Scanning of objects from multiple files according a format that follows from the extension of each fi...
JSUPPORT::JMultipleFileScanner
General purpose class for object reading from a list of file names.
Definition: JMultipleFileScanner.hh:167
JParser.hh
Utility class to parse command line options.
JSUPPORT::JLimit::getLimit
const JLimit & getLimit() const
Get limit.
Definition: JLimit.hh:73
JGeanz.hh
Longitudinal emission profile EM-shower.
JTOOLS::u
double u[N+1]
Definition: JPolint.hh:706
Evt::mc_trks
std::vector< Trk > mc_trks
MC: list of MC truth tracks.
Definition: Evt.hh:45
Evt
The Evt class respresent a Monte Carlo (MC) event as well as an offline event.
Definition: Evt.hh:19
E
then usage $script[input file[working directory[option]]] nWhere option can be E
Definition: JMuonPostfit.sh:37
P
then $DIR JPlotNPE PDG P
Definition: JPlotNPE-PDG.sh:60
DEBUG
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62
main
int main(int argc, char *argv[])
Definition: Main.cpp:15
Generated by   doxygen 1.8.5