Jpp 20.0.0-27-g39925593c-D
the software that should make you happy
Loading...
Searching...
No Matches
Functions
JTimesliceReprocessor.cc File Reference

Auxiliary program to re-process KM3NETDAQ::JDAQTimeslice data. More...

#include <string>
#include <iostream>
#include <iomanip>
#include "km3net-dataformat/offline/Head.hh"
#include "km3net-dataformat/offline/MultiHead.hh"
#include "km3net-dataformat/online/JDAQ.hh"
#include "km3net-dataformat/online/JDAQClock.hh"
#include "JDAQ/JDAQTimesliceIO.hh"
#include "JDAQ/JDAQEventIO.hh"
#include "JDAQ/JDAQSummarysliceIO.hh"
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JDetector/JModuleRouter.hh"
#include "JTrigger/JHit.hh"
#include "JTrigger/JHitToolkit.hh"
#include "JTrigger/JSuperFrame1D.hh"
#include "JTrigger/JSuperFrame2D.hh"
#include "JTrigger/JTimeslice.hh"
#include "JTrigger/JBuildL1.hh"
#include "JTrigger/JBuildL2.hh"
#include "JTrigger/JTimesliceL1.hh"
#include "JTrigger/JTriggerParameters.hh"
#include "JTrigger/JTriggerToolkit.hh"
#include "JLang/JObjectMultiplexer.hh"
#include "JROOT/JROOTClassSelector.hh"
#include "JSupport/JMonteCarloFileSupportkit.hh"
#include "JSupport/JSingleFileScanner.hh"
#include "JSupport/JFileRecorder.hh"
#include "JSupport/JSupport.hh"
#include "JSupport/JMeta.hh"
#include "JSupport/JTriggerParametersSupportkit.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 re-process KM3NETDAQ::JDAQTimeslice data.

Author
mdejong

Definition in file JTimesliceReprocessor.cc.

Function Documentation

◆ main()

int main ( int argc,
char ** argv )

Definition at line 49 of file JTimesliceReprocessor.cc.

50{
51 using namespace std;
52 using namespace JPP;
53 using namespace KM3NETDAQ;
54
55 typedef JAllTypes_t typelist;
56
57 JSingleFileScanner<JDAQTimesliceTypes_t> inputFile;
58 JFileRecorder<typelist> outputFile;
59 JLimit_t& numberOfEvents = inputFile.getLimit();
60 string detectorFile;
61 JTriggerParameters parameters;
62 bool reuse_parameters;
63 JROOTClassSelector selector;
64 int debug;
65
66 try {
67
68 JParser<> zap("Auxiliary program to re-process time slice data.");
69
70 zap['f'] = make_field(inputFile);
71 zap['o'] = make_field(outputFile) = "timeslice.root";
72 zap['n'] = make_field(numberOfEvents) = JLimit::max();
73 zap['a'] = make_field(detectorFile);
74 zap['@'] = make_field(parameters) = JPARSER::initialised();
75 zap['U'] = make_field(reuse_parameters);
76 zap['C'] = make_field(selector, "timeslice selector, e.g. JDAQTimesliceL1.") = getROOTClassSelection<JDAQTimesliceTypes_t>();
77 zap['d'] = make_field(debug) = 0;
78
79 zap(argc, argv);
80 }
81 catch(const exception &error) {
82 FATAL(error.what() << endl);
83 }
84
85
86 setDAQLongprint(debug >= JEEP::debug_t);
87
88 DEBUG("Frame time [ms] " << getFrameTime() * 1e-6 << endl);
89 DEBUG("Reset time [ms] " << getRTS() * 1e-6 << endl);
90 DEBUG("Trigger" << endl << parameters << endl);
91
92 JDetector detector;
93
94 try {
95 load(detectorFile, detector);
96 }
97 catch(const JException& error) {
98 FATAL(error);
99 }
100
101 if (reuse_parameters) {
102
103 try {
104
105 parameters = getTriggerParameters(inputFile);
106
107 NOTICE("Set trigger parameters from input." << endl);
108 }
109 catch(const JException& error) {
110 FATAL("No trigger parameters from input." << endl);
111 }
112 }
113 // detector
114
115 if (parameters.disableHighRateVeto) {
116
117 NOTICE("Disabling high-rate veto of all PMTs." << endl);
118
119 detector.setPMTStatus(HIGH_RATE_VETO_DISABLE);
120 }
121
122 parameters.set(getMaximalDistance(detector));
123
124 const JModuleRouter moduleRouter(detector);
125
126 if (parameters.writeSummary()) { WARNING("Discard writeSummary option during reprocesing of data." << endl); }
127
128 //typedef JHit hit_type;
129 //typedef int hit_type;
130 typedef double hit_type;
131
132 typedef JSuperFrame1D<hit_type> JSuperFrame1D_t;
133 typedef JSuperFrame2D<hit_type> JSuperFrame2D_t;
134 typedef JTimeslice <hit_type> JTimeslice_t;
135 typedef JBuildL1 <hit_type> JBuildL1_t;
136 typedef JBuildL2 <hit_type> JBuildL2_t;
137
138 const JBuildL1_t buildL1(parameters);
139 const JBuildL2_t buildL2(parameters.L2);
140 const JBuildL2_t buildSN(parameters.SN);
141
142 JTimesliceRouter timesliceRouter(parameters.numberOfBins);
143
144 outputFile.open();
145
146 if (!outputFile.is_open()) {
147 FATAL("Error opening file " << outputFile << endl);
148 }
149
150 outputFile.put(JMeta(argc, argv));
151 outputFile.put(parameters);
152
153 struct counts_type {
154
155 void add(const size_t value)
156 {
157 this->value += value;
158 this->count += 1;
159 }
160
161 size_t value = 0; // number of entries
162 size_t count = 0; // number of optical modules
163 };
164
165 map<string, counts_type> ls;
166
167 JObjectMultiplexer<JDAQTimesliceTypes_t, JDAQTimeslice> in(inputFile, selector);
168
169 counter_type counter = 0;
170
171 for ( ; in.hasNext() && counter != inputFile.getLimit(); ++counter) {
172
173 STATUS("event: " << setw(10) << counter << '\r'); DEBUG(endl);
174
175 const JDAQTimeslice* timeslice = in.next();
176
177 DEBUG(*timeslice << endl);
178
179 timesliceRouter.configure(*timeslice);
180
181 JTimeslice_t timesliceL0(timeslice->getDAQChronometer());
182 JTimeslice_t timesliceL1(timeslice->getDAQChronometer());
183 JTimeslice_t timesliceL2(timeslice->getDAQChronometer());
184 JTimeslice_t timesliceSN(timeslice->getDAQChronometer());
185
186 for (JDAQTimeslice::const_iterator super_frame = timeslice->begin(); super_frame != timeslice->end(); ++super_frame) {
187
188 if (moduleRouter.hasModule(super_frame->getModuleID())) {
189
190 // calibration
191
192 const JModule& module = moduleRouter.getModule(super_frame->getModuleID());
193 JSuperFrame2D_t& buffer = JSuperFrame2D_t::demultiplex(*super_frame, module);
194
195 // Apply high-rate veto
196
197 buffer.applyHighRateVeto(parameters.highRateVeto_Hz);
198
199 // L0
200
201 timesliceL0.push_back(JSuperFrame1D_t(buffer));
202
203 // L1
204
205 timesliceL1.push_back(JSuperFrame1D_t(super_frame->getDAQChronometer(),
206 super_frame->getModuleIdentifier(),
207 module.getPosition()));
208
209 buildL1(*timesliceL0.rbegin(), back_inserter(*timesliceL1.rbegin()));
210
211 // L2
212
213 timesliceL2.push_back(JSuperFrame1D_t(super_frame->getDAQChronometer(),
214 super_frame->getModuleIdentifier(),
215 module.getPosition()));
216
217 buildL2(buffer, *timesliceL1.rbegin(), back_inserter(*timesliceL2.rbegin()));
218
219 // SN
220 {
221 JTimeslice_t::value_type tv(super_frame->getDAQChronometer(),
222 super_frame->getModuleIdentifier(),
223 module.getPosition());
224
225 buildSN(buffer, *timesliceL1.rbegin(), back_inserter(tv));
226
227 if (!tv.empty()) {
228 timesliceSN.push_back(tv);
229 }
230 }
231 }
232 }
233
234 for (const auto& i : timesliceL1) { ls["L1"].add(i.size()); }
235 for (const auto& i : timesliceL2) { ls["L2"].add(i.size()); }
236 for (const auto& i : timesliceSN) { ls["SN"].add(i.size()); }
237
238 if (parameters.writeL0()) {
239 outputFile.put(*timeslice);
240 }
241
242 if (parameters.writeL1()) {
243 outputFile.put(JTimesliceL1<JDAQTimesliceL1>(timesliceL1, timesliceRouter, moduleRouter, parameters.TMaxLocal_ns));
244 }
245
246 if (parameters.writeL2()) {
247 outputFile.put(JTimesliceL1<JDAQTimesliceL2>(timesliceL2, timesliceRouter, moduleRouter, parameters.L2.TMaxLocal_ns));
248 }
249
250 if (parameters.writeSN()) {
251 outputFile.put(JTimesliceL1<JDAQTimesliceSN>(timesliceSN, timesliceRouter, moduleRouter, parameters.SN.TMaxLocal_ns));
252 }
253 }
254 STATUS(endl);
255
256 for (const auto& i : ls) {
257 NOTICE(left << setw(4) << i.first << right << FIXED(7,2) << (i.second.count != 0 ? (double) i.second.value / ((double) i.second.count * 1.0e-9 * getFrameTime()) : 0.0) << " Hz" << endl);
258 }
259
260 JSingleFileScanner<JRemove<typelist, JDAQTimesliceTypes_t>::typelist> io(inputFile);
261
262 io >> outputFile;
263
264 outputFile.close();
265}
outputFile
string outputFile
Definition JDAQTimesliceSelector.cc:37
DEBUG
#define DEBUG(A)
Message macros.
Definition JMessage.hh:62
STATUS
#define STATUS(A)
Definition JMessage.hh:63
NOTICE
#define NOTICE(A)
Definition JMessage.hh:64
FATAL
#define FATAL(A)
Definition JMessage.hh:67
debug
int debug
debug level
Definition JSirene.cc:72
WARNING
#define WARNING(A)
Definition JMessage.hh:65
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition JParser.hh:2142
JDETECTOR::JDetector
Detector data structure.
Definition JDetector.hh:96
JDETECTOR::JModuleRouter
Router for direct addressing of module data in detector data structure.
Definition JModuleRouter.hh:38
JDETECTOR::JModule
Data structure for a composite optical module.
Definition JModule.hh:76
JLANG::JException
General exception.
Definition JException.hh:24
JOSCPROB::JOscProbInterpolator
Template definition of a multi-dimensional oscillation probability interpolation table.
Definition JOscProbInterpolator.hh:71
KM3NETDAQ::JDAQChronometer::getDAQChronometer
const JDAQChronometer & getDAQChronometer() const
Get DAQ chronometer.
Definition JDAQChronometer.hh:88
JDETECTOR::load
void load(const std::string &file_name, JDetector &detector)
Load detector from input file.
Definition JDetectorToolkit.hh:520
JDETECTOR::getMaximalDistance
double getMaximalDistance(const JDetector &detector, const bool option=false)
Get maximal distance between modules in detector.
Definition JDetectorToolkit.hh:78
JEEP::debug_t
@ debug_t
debug
Definition JMessage.hh:29
JPP
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Definition JAAnetToolkit.hh:45
JSUPPORT::getTriggerParameters
JTriggerParameters getTriggerParameters(const JMultipleFileScanner_t &file_list)
Get trigger parameters.
Definition JTriggerParametersSupportkit.hh:148
KM3NETDAQ
KM3NeT DAQ data structures and auxiliaries.
Definition DataQueue.cc:39
KM3NETDAQ::getFrameTime
double getFrameTime()
Get frame time duration.
Definition JDAQClock.hh:162
KM3NETDAQ::setDAQLongprint
void setDAQLongprint(const bool option)
Set DAQ print option.
Definition JDAQPrint.hh:28
KM3NETDAQ::getRTS
double getRTS()
Get TDC dynamic range.
Definition JDAQClock.hh:173
HIGH_RATE_VETO_DISABLE
static const int HIGH_RATE_VETO_DISABLE
Enable (disable) use of high-rate veto test if this status bit is 0 (1);.
Definition pmt_status.hh:14
FIXED
Auxiliary data structure for floating point format specification.
Definition JManip.hh:448
JAANET::detector
Detector file.
Definition JHead.hh:227
JACOUSTICS::hit_type
Transmission with position.
Definition JAcousticsTriggerProcessor.cc:72
JPARSER::initialised
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Definition JParser.hh:68
JROOT::JROOTClassSelector
Auxiliary class to select ROOT class based on class name.
Definition JROOTClassSelector.hh:38
JSUPPORT::JLimit
Auxiliary class for defining the range of iterations of objects.
Definition JLimit.hh:45
JSUPPORT::JLimit::max
static counter_type max()
Get maximum counter value.
Definition JLimit.hh:128
JSUPPORT::JMeta
Auxiliary class for ROOT I/O of application specific meta data.
Definition JMeta.hh:72
JSYSTEM::ls
Auxiliary data structure to list files in directory.
Definition JFilesystem.hh:20
JTRIGGER::JTimesliceRouter
Router for fast addressing of hits in KM3NETDAQ::JDAQTimeslice data structure as a function of the op...
Definition JTimesliceRouter.hh:94
Generated by doxygen 1.10.0