Jpp master_rocky-44-g75b7c4f75
the software that should make you happy
Loading...
Searching...
No Matches
Functions
JTriggerProcessor.cc File Reference

Auxiliary program to trigger KM3NETDAQ::JDAQTimeslice data. More...

#include <string>
#include <iostream>
#include <iomanip>
#include <set>
#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/JHitL0.hh"
#include "JTrigger/JHitL1.hh"
#include "JTrigger/JBuildL1.hh"
#include "JTrigger/JBuildL2.hh"
#include "JTrigger/JTrigger3DShower.hh"
#include "JTrigger/JTriggerMXShower.hh"
#include "JTrigger/JTrigger3DMuon.hh"
#include "JTrigger/JTriggerNB.hh"
#include "JTrigger/JTriggerBits.hh"
#include "JTrigger/JEventOverlap.hh"
#include "JTrigger/JTimesliceRouter.hh"
#include "JTrigger/JTriggeredEvent.hh"
#include "JTrigger/JTimesliceL1.hh"
#include "JTrigger/JTriggerParameters.hh"
#include "JTrigger/JTriggerToolkit.hh"
#include "JROOT/JROOTClassSelector.hh"
#include "JLang/JObjectMultiplexer.hh"
#include "JSupport/JMonteCarloFileSupportkit.hh"
#include "JSupport/JMultipleFileScanner.hh"
#include "JSupport/JFileRecorder.hh"
#include "JSupport/JSupport.hh"
#include "JSupport/JMeta.hh"
#include "Jeep/JTimer.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 trigger KM3NETDAQ::JDAQTimeslice data.

Random KM3NETDAQ::JDAQTimeslice data can be produced with JRandomTimesliceWriter.cc application.

Author
mdejong

Definition in file JTriggerProcessor.cc.

Function Documentation

◆ main()

int main ( int argc,
char ** argv )

Definition at line 62 of file JTriggerProcessor.cc.

63{
64 using namespace std;
65 using namespace KM3NETDAQ;
66 using namespace JPP;
67
68 typedef JTYPELIST<JAllTypes_t>::typelist typelist;
69
70 JMultipleFileScanner<JDAQTimesliceTypes_t> inputFile;
71 JLimit_t& numberOfEvents = inputFile.getLimit();
72 JFileRecorder<typelist> outputFile;
73 string detectorFile;
74 JTriggerParameters parameters;
75 JDAQClock clock;
76 JROOTClassSelector selector;
77 bool profile;
78 int debug;
79
80 try {
81
82 JParser<> zap("Auxiliary program to trigger time slice data.");
83
84 zap['f'] = make_field(inputFile, "input file.");
85 zap['o'] = make_field(outputFile, "output file.") = "trigger_processor.dat";
86 zap['n'] = make_field(numberOfEvents) = JLimit::max();
87 zap['a'] = make_field(detectorFile, "detector file.");
88 zap['@'] = make_field(parameters, "trigger parameters") = JPARSER::initialised();
89 zap['c'] = make_field(clock, "DAQ clock") = JDAQClock::KM3NET, JDAQClock::PPM_DU, JDAQClock::ANTARES;
90 zap['C'] = make_field(selector, "timeslice selector, e.g. JDAQTimesliceL1.") = getROOTClassSelection<JDAQTimesliceTypes_t>();
91 zap['P'] = make_field(profile, "profile CPU");
92 zap['d'] = make_field(debug, "debug flag.") = 1;
93
94 zap(argc, argv);
95 }
96 catch(const exception &error) {
97 FATAL(error.what() << endl);
98 }
99
100
101 clock.set();
102
103 setDAQLongprint(debug >= JEEP::debug_t);
104
105 if (!profile) {
106 JTimer::disable();
107 }
108
109 DEBUG("Frame time [ms] " << getFrameTime() * 1e-6 << endl);
110 DEBUG("Reset time [ms] " << getRTS() * 1e-6 << endl);
111 DEBUG("Trigger" << endl << parameters << endl);
112
113 JDetector detector;
114
115 try {
116 load(detectorFile, detector);
117 }
118 catch(const JException& error) {
119 FATAL(error);
120 }
121
122 // detector
123
124 if (parameters.disableHighRateVeto) {
125
126 NOTICE("Disabling high-rate veto of all PMTs." << endl);
127
128 detector.setPMTStatus(HIGH_RATE_VETO_DISABLE);
129 }
130
131 parameters.set(getMaximalDistance(detector));
132
133 const JModuleRouter moduleRouter(detector);
134
135 //typedef JHit hit_type;
136 //typedef int hit_type;
137 typedef double hit_type;
138
139 typedef JSuperFrame1D<hit_type> JSuperFrame1D_t;
140 typedef JSuperFrame2D<hit_type> JSuperFrame2D_t;
141 typedef JTimeslice <hit_type> JTimeslice_t;
142 typedef JBuildL1 <hit_type> JBuildL1_t;
143 typedef JBuildL2 <hit_type> JBuildL2_t;
144
145 const JBuildL1_t buildL1(parameters);
146 const JBuildL2_t buildL2(parameters.L2);
147 const JBuildL2_t buildSN(parameters.SN);
148 const JBuildL2_t buildNB(parameters.NB);
149
150 JTimer timerCC("Calibration");
151 JTimer timerL0("L0");
152 JTimer timerL1("L1");
153 JTimer timerL2("L2");
154 JTimer timerSN("SN");
155 JTimer timerRX("Timeslice router");
156 JTimer timerTR("Trigger");
157 JTimer timerTX("Trigger router");
158 JTimer timerTW("Timeslice writer");
159 JTimer timerSW("Summary writer");
160
161 JTimesliceRouter timesliceRouter(parameters.numberOfBins);
162
163 const JTriggerNB triggerNB (parameters);
164 const JTrigger3DMuon trigger3DMuon (parameters);
165 const JTrigger3DShower trigger3DShower(parameters);
166 const JTriggerMXShower triggerMXShower(parameters, detector);
167
168
169 outputFile.open();
170
171 if (!outputFile.is_open()) FATAL("Error opening file " << outputFile << endl);
172
173 outputFile.put(JMeta(argc, argv));
174 outputFile.put(parameters);
175
176 unsigned int numberOfTriggers = 0;
177
178 JObjectMultiplexer<JDAQTimesliceTypes_t, JDAQTimeslice> in(inputFile, selector);
179
180 counter_type counter = 0;
181
182 for ( ; in.hasNext() && counter != inputFile.getLimit(); ++counter) {
183
184 STATUS("event: " << setw(10) << counter << '\r'); DEBUG(endl);
185
186 const JDAQTimeslice* timeslice = in.next();
187
188 DEBUG(*timeslice << endl);
189
190 timerRX.start();
191
192 timesliceRouter.configure(*timeslice);
193
194 timerRX.stop();
195
196
197 JTimeslice_t timesliceL0(timeslice->getDAQChronometer());
198 JTimeslice_t timesliceL1(timeslice->getDAQChronometer());
199 JTimeslice_t timesliceL2(timeslice->getDAQChronometer());
200 JTimeslice_t timesliceSN(timeslice->getDAQChronometer());
201 JTimeslice_t timesliceNB(timeslice->getDAQChronometer());
202
203 for (JDAQTimeslice::const_iterator super_frame = timeslice->begin(); super_frame != timeslice->end(); ++super_frame) {
204
205 if (moduleRouter.hasModule(super_frame->getModuleID())) {
206
207 // calibration
208
209 timerCC.start();
210
211 const JModule& module = moduleRouter.getModule(super_frame->getModuleID());
212 JSuperFrame2D_t& buffer = JSuperFrame2D_t::demultiplex(*super_frame, module);
213
214 timerCC.stop();
215
216 // Apply high-rate veto
217
218 buffer.applyHighRateVeto(parameters.highRateVeto_Hz);
219
220 // L0
221
222 timerL0.start();
223
224 timesliceL0.push_back(JSuperFrame1D_t(buffer));
225
226 timerL0.stop();
227
228 // Nano-beacon trigger
229
230 if (parameters.triggerNB.enabled) {
231
232 JSuperFrame2D_t::iterator __end = partition(buffer.begin(), buffer.end(), parameters.triggerNB.pmts);
233
234 if (buffer.begin() != __end) {
235
236 timesliceNB.push_back(JSuperFrame1D_t(super_frame->getDAQChronometer(),
237 super_frame->getModuleIdentifier(),
238 module.getPosition()));
239
240 JSuperFrame1D_t zbuf;
241
242 buildL1(buffer.begin(), __end , back_inserter(zbuf));
243
244 buildNB(buffer.begin() , __end, zbuf, back_inserter(*timesliceNB.rbegin()));
245 }
246 }
247
248 // L1
249
250 timerL1.start();
251
252 timesliceL1.push_back(JSuperFrame1D_t(super_frame->getDAQChronometer(),
253 super_frame->getModuleIdentifier(),
254 module.getPosition()));
255
256 buildL1(*timesliceL0.rbegin(), back_inserter(*timesliceL1.rbegin()));
257
258 timerL1.stop();
259
260 // L2
261
262 timerL2.start();
263
264 timesliceL2.push_back(JSuperFrame1D_t(super_frame->getDAQChronometer(),
265 super_frame->getModuleIdentifier(),
266 module.getPosition()));
267
268 buildL2(buffer, *timesliceL1.rbegin(), back_inserter(*timesliceL2.rbegin()));
269
270 timerL2.stop();
271
272 // SN
273
274 timerSN.start();
275
276 timesliceSN.push_back(JSuperFrame1D_t(super_frame->getDAQChronometer(),
277 super_frame->getModuleIdentifier(),
278 module.getPosition()));
279
280 buildSN(buffer, *timesliceL1.rbegin(), back_inserter(*timesliceSN.rbegin()));
281
282 timerSN.stop();
283 }
284 }
285
286 // Trigger
287
288 if (parameters.triggerNB.enabled) {
289
290 timerTX.start();
291
292 const JTriggerInput trigger_input(timesliceNB);
293
294 for (JTriggerInput::const_iterator hit = trigger_input.begin(); hit != trigger_input.end(); ++hit) {
295
296 if (parameters.triggerNB.write()) {
297
298 JTriggeredEvent tev(timesliceNB.getDAQChronometer(),
299 getTriggerMask(triggerNB.getTriggerBit()),
300 *hit,
301 timesliceRouter,
302 moduleRouter,
303 parameters.TMaxLocal_ns,
304 parameters.triggerNB.DMax_m,
305 getTimeRange(parameters.triggerNB));
306
307 outputFile.put(tev);
308 }
309 }
310
311 timerTX.stop();
312 }
313
314 timerTR.start();
315
316 JTriggerInput trigger_input(timesliceL2);
317 JTriggerOutput trigger_output;
318
319 trigger3DMuon (trigger_input, back_inserter(trigger_output));
320 trigger3DShower(trigger_input, back_inserter(trigger_output));
321 triggerMXShower(trigger_input, timesliceL0, back_inserter(trigger_output));
322
323 trigger_output.merge(JEventOverlap(parameters.TMaxEvent_ns));
324
325 numberOfTriggers += trigger_output.size();
326
327 timerTR.stop();
328
329 DEBUG("Number of triggers: " << trigger_output.size() << endl);
330
331 for (JTriggerOutput::const_iterator event = trigger_output.begin(); event != trigger_output.end(); ++event) {
332
333 timerTX.start();
334
335 JTriggeredEvent tev(*event,
336 timesliceRouter,
337 moduleRouter,
338 parameters.TMaxLocal_ns,
339 getTimeRange(parameters));
340
341 outputFile.put(tev);
342
343 timerTX.stop();
344 }
345
346 if (parameters.writeL0()) {
347
348 timerTW.start();
349
350 outputFile.put(*timeslice);
351
352 timerTW.stop();
353 }
354
355 if (parameters.writeL1()) {
356
357 timerTW.start();
358
359 outputFile.put(JTimesliceL1<JDAQTimesliceL1>(timesliceL1, timesliceRouter, moduleRouter, parameters.TMaxLocal_ns));
360
361 timerTW.stop();
362 }
363
364 if (parameters.writeL2()) {
365
366 timerTW.start();
367
368 outputFile.put(JTimesliceL1<JDAQTimesliceL2>(timesliceL2, timesliceRouter, moduleRouter, parameters.L2.TMaxLocal_ns));
369
370 timerTW.stop();
371 }
372
373 if (parameters.writeSN()) {
374
375 timerTW.start();
376
377 outputFile.put(JTimesliceL1<JDAQTimesliceSN>(timesliceSN, timesliceRouter, moduleRouter, parameters.SN.TMaxLocal_ns));
378
379 timerTW.stop();
380 }
381
382 if (parameters.writeSummary()) {
383
384 timerSW.start();
385
386 outputFile.put(JDAQSummaryslice(*timeslice));
387
388 timerSW.stop();
389 }
390 }
391 STATUS(endl);
392
393 if (debug >= JEEP::notice_t&& counter != 0) {
394
395 const double factor = 1.0 / (double) counter;
396
397 for (const JTimer* p : { &timerCC, &timerL0, &timerL1, &timerL2, &timerSN, &timerRX, &timerTR, &timerTX, &timerTW, &timerSW }) {
398 p->print(cout, factor);
399 }
400
401 NOTICE("Number of trigger/slices " << numberOfTriggers << "/" << counter << endl);
402 NOTICE("Trigger rate [Hz] " << numberOfTriggers * 1.0e9 * factor / getFrameTime() << endl);
403 }
404
405 JMultipleFileScanner<JTYPELIST<JAAnetTypes_t, JMetaTypes_t>::typelist> io(inputFile);
406
407 io >> outputFile;
408
409 outputFile.close();
410}
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:69
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition JParser.hh:2142
JACOUSTICS::JTriggerOutput::merge
void merge(const JMatch_t &match)
Merge events.
Definition JAcoustics/JTriggerOutput.hh:40
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:75
JEEP::JTimer
Auxiliary class for CPU timing and usage.
Definition JTimer.hh:33
JEEP::JTimer::disable
static void disable()
Disable timers.
Definition JTimer.hh:84
JGEOMETRY3D::JPosition3D::getPosition
const JPosition3D & getPosition() const
Get position.
Definition JPosition3D.hh:130
JLANG::JException
General exception.
Definition JException.hh:24
JLANG::JObjectMultiplexer
Auxiliary class for multiplexing object iterators.
Definition JObjectMultiplexer.hh:35
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
General purpose class for object reading from a list of file names.
Definition JMultipleFileScanner.hh:306
JTRIGGER::JSuperFrame1D
1-dimensional frame with time calibrated data from one optical module.
Definition JSuperFrame1D.hh:40
JTRIGGER::JSuperFrame2D
2-dimensional frame with time calibrated data from one optical module.
Definition JSuperFrame2D.hh:46
JTRIGGER::JTimesliceL1
Auxiliary class to build JDAQTimeslice for L1 timeslice.
Definition JTimesliceL1.hh:38
JTRIGGER::JTriggerInput
Data structure for input to trigger algorithm.
Definition JTriggerInput.hh:36
JTRIGGER::JTriggerNB
Nano-beacon trigger.
Definition JTriggerNB.hh:21
JTRIGGER::JTriggeredEvent
Auxiliary class to build JDAQEvent for a triggered event.
Definition JTriggeredEvent.hh:50
KM3NETDAQ::JDAQChronometer::getDAQChronometer
const JDAQChronometer & getDAQChronometer() const
Get DAQ chronometer.
Definition JDAQChronometer.hh:88
KM3NETDAQ::JDAQClock
Auxiliary class to set DAQ system clock parameters.
Definition JDAQClock.hh:27
KM3NETDAQ::JDAQClock::set
void set(const JDAQClock_t clock)
Set clock type.
Definition JDAQClock.hh:64
JDETECTOR::load
void load(const std::string &file_name, JDetector &detector)
Load detector from input file.
Definition JDetectorToolkit.hh:505
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
JEEP::notice_t
@ notice_t
notice
Definition JMessage.hh:32
JPP
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Definition JAAnetToolkit.hh:43
JROOT::getROOTClassSelection
std::set< JROOTClassSelector > getROOTClassSelection(const bool option=false)
Get ROOT class selection.
Definition JROOTClassSelector.hh:318
JROOT::counter_type
Long64_t counter_type
Type definition for counter.
Definition JROOT/JCounter.hh:24
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:13
JAANET::detector
Detector file.
Definition JHead.hh:227
JACOUSTICS::JEventOverlap
Match of two events considering overlap in time.
Definition JAcoustics/JEventOverlap.hh:20
JACOUSTICS::hit_type
Acoustic hit.
Definition JBillabong.cc:70
JLANG::JTypeList
Type list.
Definition JTypeList.hh:23
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
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.12.0