Jpp  18.3.0-rc.1
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Functions
JStingray.cc File Reference

Auxiliary program to manipulate MUPAGE data. More...

#include <string>
#include <iostream>
#include <iomanip>
#include <map>
#include "TRandom3.h"
#include "km3net-dataformat/offline/Head.hh"
#include "km3net-dataformat/offline/MultiHead.hh"
#include "km3net-dataformat/offline/Evt.hh"
#include "JSupport/JMultipleFileScanner.hh"
#include "JSupport/JMonteCarloFileSupportkit.hh"
#include "JSupport/JFileRecorder.hh"
#include "JSupport/JSupport.hh"
#include "JSupport/JMeta.hh"
#include "JPhysics/JGeane.hh"
#include "JPhysics/JGeanz.hh"
#include "JPhysics/JConstants.hh"
#include "Jeep/JPrint.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"

Go to the source code of this file.

Functions

int main (int argc, char **argv)
 

Detailed Description

Auxiliary program to manipulate MUPAGE data.

Author
mdejong

Definition in file JStingray.cc.

Function Documentation

int main ( int  argc,
char **  argv 
)

Definition at line 34 of file JStingray.cc.

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 }
JSUPPORT::JMeta
Auxiliary class for ROOT I/O of application specific meta data.
Definition: JMeta.hh:70
JSUPPORT::JFileRecorder
Object writing to file.
Definition: JFileRecorder.hh:42
JPARSER::JParser
Utility class to parse command line options.
Definition: JParser.hh:1514
E
then usage $script< input file >[option[primary[working directory]]] nWhere option can be E
Definition: JMuonPostfit.sh:40
STATUS
#define STATUS(A)
Definition: JMessage.hh:63
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
outputFile
string outputFile
Definition: JDAQTimesliceSelector.cc:37
JSUPPORT::JLimit
Auxiliary class for defining the range of iterations of objects.
Definition: JLimit.hh:41
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition: JParser.hh:1989
JMATH::pow
T pow(const T &x, const double y)
Power .
Definition: JMath.hh:97
JSUPPORT::JMultipleFileScanner< Head >
Template specialisation of JMultipleFileScanner for Monte Carlo header.
Definition: JMonteCarloFileSupportkit.hh:406
FATAL
#define FATAL(A)
Definition: JMessage.hh:67
JPHYSICS::getSpeedOfLight
const double getSpeedOfLight()
Get speed of light.
Definition: JPhysics/JConstants.hh:112
JSUPPORT::JMultipleFileScanner
General purpose class for object reading from a list of file names.
Definition: JMultipleFileScanner.hh:202
JSUPPORT::JLimit::getLimit
const JLimit & getLimit() const
Get limit.
Definition: JLimit.hh:84
JTOOLS::u
double u[N+1]
Definition: JPolint.hh:865
debug
int debug
debug level
Definition: archive-put-wiki-detectors.sh:92
Evt
The Evt class respresent a Monte Carlo (MC) event as well as an offline event.
Definition: Evt.hh:20
P
then $DIR JPlotNPE PDG P
Definition: JPlotNPE-PDG.sh:62
DEBUG
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62
Generated by   doxygen 1.8.5