Jpp test-rotations-old
the software that should make you happy
Loading...
Searching...
No Matches
JMuonGandalf.cc
Go to the documentation of this file.
1#include <string>
2#include <iostream>
3#include <iomanip>
4#include <memory>
5
9
10#include "JDAQ/JDAQEventIO.hh"
13
15
19
21
22#include "JSupport/JSupport.hh"
26#include "JSupport/JMeta.hh"
27
32
33#include "Jeep/JParser.hh"
34#include "Jeep/JMessage.hh"
35
36
37/**
38 * \file
39 *
40 * Program to perform JFIT::JRegressor<JFIT::JLine3Z,JFIT::JGandalf> fit with I/O of JFIT::JEvt data.
41 *
42 * \author mdejong, gmaggi
43 */
44int main(int argc, char **argv)
45{
46 using namespace std;
47 using namespace JPP;
48 using namespace KM3NETDAQ;
49
51 typedef JParallelFileScanner< JTypeList<JDAQEvent, JFIT::JEvt> > JParallelFileScanner_t;
52 typedef JParallelFileScanner_t::multi_pointer_type multi_pointer_type;
54 JACOUSTICS::JEvt>::typelist calibration_types;
55 typedef JMultipleFileScanner<calibration_types> JCalibration_t;
56
57 JParallelFileScanner_t inputFile;
59 JLimit_t& numberOfEvents = inputFile.getLimit();
60 string detectorFile;
61 JCalibration_t calibrationFile;
62 double Tmax_s;
63 string pdfFile;
64 JMuonGandalfParameters_t parameters;
65 size_t threads; // number of parallel threads
66 int debug;
67
68 try {
69
70 JParser<> zap("Program to perform PDF fit of muon trajectory to data.");
71
72 zap['f'] = make_field(inputFile);
73 zap['o'] = make_field(outputFile) = "gandalf.root";
74 zap['a'] = make_field(detectorFile);
75 zap['+'] = make_field(calibrationFile) = JPARSER::initialised();
76 zap['T'] = make_field(Tmax_s) = 100.0;
77 zap['n'] = make_field(numberOfEvents) = JLimit::max();
78 zap['P'] = make_field(pdfFile);
79 zap['@'] = make_field(parameters) = JPARSER::initialised();
80 zap['N'] = make_field(threads, "number of threads") = 0;
81 zap['d'] = make_field(debug) = 1;
82
83 zap(argc, argv);
84 }
85 catch(const exception& error) {
86 FATAL(error.what() << endl);
87 }
88
89
91
93
94 try {
95 load(detectorFile, detector);
96 }
97 catch(const JException& error) {
98 FATAL(error);
99 }
100
101 unique_ptr<JDynamics> dynamics;
102
103 if (!calibrationFile.empty()) {
104
105 try {
106
107 dynamics.reset(new JDynamics(detector, Tmax_s));
108
109 dynamics->load(calibrationFile);
110 }
111 catch(const exception& error) {
112 FATAL(error.what());
113 }
114 }
115
116 outputFile.open();
117 outputFile.put(JMeta(argc, argv));
118
119 {
120 const JModuleRouter router(dynamics ? dynamics->getDetector() : detector);
121
122 JSummaryFileRouter summary(inputFile);
123
124 const JMuonGandalf::storage_type storage(pdfFile, JTimeRange(parameters.TMin_ns, parameters.TMax_ns), parameters.TTS_ns);
125
126 JMuonGandalf fit(parameters, storage, debug);
127
128 JMultiThreadedReconstruction<JMuonGandalf> fits(fit, outputFile, threads, 2 * threads);
129
130 while (inputFile.hasNext()) {
131
132 STATUS("event: " << setw(10) << inputFile.getCounter() << '\r'); DEBUG(endl);
133
134 multi_pointer_type ps = inputFile.next();
135
136 const JDAQEvent* tev = ps;
137 const JFIT::JEvt* in = ps;
138
139 summary.update(*tev);
140
141 if (dynamics) {
142 dynamics->update(*tev);
143 }
144
145 auto input = fit.getInput(router, summary, *tev, *in, dynamics ? dynamics->getCoverage() : coverage_type());
146
147 if (threads > 0)
148 fits.enqueue(input);
149 else
150 outputFile.put(fit(input));
151 }
152 STATUS(endl);
153 }
154
156
157 io >> outputFile;
158
159 outputFile.close();
160}
161
string outputFile
Data structure for detector geometry and calibration.
Dynamic detector calibration.
Recording of objects on file according a format that follows from the file name extension.
General purpose messaging.
#define DEBUG(A)
Message macros.
Definition JMessage.hh:62
#define STATUS(A)
Definition JMessage.hh:63
#define FATAL(A)
Definition JMessage.hh:67
int debug
debug level
Definition JSirene.cc:72
ROOT I/O of application specific meta data.
Direct access to module in detector data structure.
int main(int argc, char **argv)
Parallel scanning of objects from a single file or multiple files according a format that follows fro...
Utility class to parse command line options.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition JParser.hh:2142
ROOT TTree parameter settings of various packages.
Detector data structure.
Definition JDetector.hh:96
Router for direct addressing of module data in detector data structure.
Data structure for set of track fit results.
General exception.
Definition JException.hh:24
Utility class to parse command line options.
Definition JParser.hh:1698
Thread pool for event-by-event reconstruction.
void setSplitLevel(int value)
Set split level.
Object writing to file.
General purpose class for object reading from a list of file names.
General purpose class for parallel reading of objects from a single file or multiple files.
Object reading from a list of files.
File router for fast addressing of summary data.
void update(const JDAQHeader &header)
Update router.
JTOOLS::JRange< double > JTimeRange
Type definition for time range (unit [s]).
void load(const std::string &file_name, JDetector &detector)
Load detector from input file.
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
JTreeParameters & getTreeParameters()
Template definition for method returning TTree parameters.
KM3NeT DAQ data structures and auxiliaries.
Definition DataQueue.cc:39
Detector file.
Definition JHead.hh:227
Acoustic event fit.
Orientation of module.
Dynamic detector calibration.
Definition JDynamics.hh:86
Data structure for coverage of detector by dynamical calibrations.
Definition JCoverage.hh:19
Auxiliary class for recursive type list generation.
Definition JTypeList.hh:351
Type list.
Definition JTypeList.hh:23
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Definition JParser.hh:68
double TMin_ns
minimal time w.r.t. Cherenkov hypothesis [ns]
double TMax_ns
maximal time w.r.t. Cherenkov hypothesis [ns]
Wrapper class to make final fit of muon trajectory.
input_type getInput(const JModuleRouter &router, const JSummaryRouter &summary, const JDAQEvent &event, const JEvt &in, const coverage_type &coverage) const
Get input data.
Auxiliary class for defining the range of iterations of objects.
Definition JLimit.hh:45
const JLimit & getLimit() const
Get limit.
Definition JLimit.hh:84
static counter_type max()
Get maximum counter value.
Definition JLimit.hh:128
Auxiliary class for ROOT I/O of application specific meta data.
Definition JMeta.hh:72