Jpp 20.0.0-195-g190c9e876
the software that should make you happy
Loading...
Searching...
No Matches
Functions
JTriggerEfficiency.cc File Reference

Auxiliary program to trigger Monte Carlo events. More...

#include <string>
#include <iostream>
#include <iomanip>
#include <vector>
#include <set>
#include "TH1D.h"
#include "TRandom3.h"
#include "km3net-dataformat/offline/Head.hh"
#include "km3net-dataformat/offline/MultiHead.hh"
#include "km3net-dataformat/offline/Evt.hh"
#include "km3net-dataformat/offline/Hit.hh"
#include "km3net-dataformat/definitions/pmt_status.hh"
#include "JDAQ/JDAQTimesliceIO.hh"
#include "JDAQ/JDAQEventIO.hh"
#include "JDAQ/JDAQSummarysliceIO.hh"
#include "JDAQ/JDAQToolkit.hh"
#include "JTimeslice/JEventTimeslice.hh"
#include "JSummaryslice/JSummaryslice.hh"
#include "JAAnet/JHead.hh"
#include "JAAnet/JHeadToolkit.hh"
#include "JAAnet/JAAnetToolkit.hh"
#include "JPhysics/JK40Rates.hh"
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JDetector/JDetectorSimulator.hh"
#include "JDetector/JModuleMapper.hh"
#include "JDetector/JPMTRouter.hh"
#include "JDetector/JTimeRange.hh"
#include "JDetector/JPMTParametersMap.hh"
#include "JDetector/JK40DefaultSimulator.hh"
#include "JDetector/JPMTDefaultSimulator.hh"
#include "JDetector/JCLBDefaultSimulator.hh"
#include "JTrigger/JHit.hh"
#include "JTrigger/JHitToolkit.hh"
#include "JTrigger/JTimeslice.hh"
#include "JTrigger/JSuperFrame1D.hh"
#include "JTrigger/JSuperFrame2D.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/JTriggerBits.hh"
#include "JTrigger/JEventOverlap.hh"
#include "JTrigger/JTimesliceRouter.hh"
#include "JTrigger/JTriggeredEvent.hh"
#include "JTrigger/JTimesliceL1.hh"
#include "JTrigger/JTriggerParameters.hh"
#include "JTrigger/JEventToolkit.hh"
#include "JTrigger/JSummaryRouter.hh"
#include "JTrigger/JTriggerToolkit.hh"
#include "JTrigger/JK40RunByRunSimulator.hh"
#include "JTrigger/JPMTRunByRunSimulator.hh"
#include "JTrigger/JCLBRunByRunSimulator.hh"
#include "JSupport/JMultipleFileScanner.hh"
#include "JSupport/JTreeScanner.hh"
#include "JSupport/JFileRecorder.hh"
#include "JSupport/JMonteCarloFileSupportkit.hh"
#include "JSupport/JTriggerParametersSupportkit.hh"
#include "JSupport/JSupportToolkit.hh"
#include "JSupport/JSupport.hh"
#include "JSupport/JMeta.hh"
#include "JROOT/JRandom.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 trigger Monte Carlo events.

The options

  • -a refers to the detector configuration that is used to convert Monte Carlo true information to raw data; and
  • -b refers to the detector configuration that is used to convert raw data to calibrated data.

The event counter of the triggered events, as returned by KM3NETDAQ::JDAQEvent::getCounter(), is set to the corresponding index of the Monte Carlo event in the output.
This allows for correlating the DAQ event to the Monte Carlo event.
The time of the DAQ hits can be correlated to the time of the Monte Carlo hits using the frame index and the time of the Monte Carlo event (Evt::mc_t), respectively (see also time_converter).
The universal time of the DAQ event is set to the universal time of the Monte Carlo event with a correction for the time of the event (read hits) within the time slice.
In run-by-run mode, the trigger parameters as well as the singles rates and status of the PMTs are read from the given raw data file.
The two-fold (and higher) coincidence rates should still be provided on the command line.

Author
mdejong

Definition in file JTriggerEfficiency.cc.

Function Documentation

◆ main()

int main ( int argc,
char ** argv )

Definition at line 134 of file JTriggerEfficiency.cc.

135{
136 using namespace std;
137 using namespace JPP;
138 using namespace KM3NETDAQ;
139
140 typedef JTYPELIST<JAllTypes_t, JRootTypes_t>::typelist typelist;
141
142 JMultipleFileScanner<Evt> inputFile;
143 JFileRecorder<typelist> outputFile;
144 JLimit_t& numberOfEvents = inputFile.getLimit();
145 string detectorFileA;
146 string detectorFileB;
147 int run;
148 JTriggerParameters parameters;
149 bool triggeredEventsOnly;
150 JPMTParametersMap pmtParameters;
151 JK40Rates rates_Hz;
152 JMixedK40Rates mixed_Hz;
153 string runbyrun;
154 double sigma_ns;
155 JRandom seed;
156 int debug;
157
158 try {
159
160 JParser<> zap("Auxiliary program to trigger Monte Carlo events.");
161
162 zap['f'] = make_field(inputFile, "input file (output of detector simulation)");
163 zap['o'] = make_field(outputFile, "output file") = "trigger_efficieny.root";
164 zap['n'] = make_field(numberOfEvents) = JLimit::max();
165 zap['a'] = make_field(detectorFileA, "detector used for conversion from Monte Carlo truth to raw data.");
166 zap['b'] = make_field(detectorFileB, "detector used for conversion of raw data to calibrated data.") = "";
167 zap['R'] = make_field(run, "run number") = -1;
168 zap['r'] = make_field(runbyrun, "DAQ file for run-by-run mode") = JPARSER::initialised();
169 zap['@'] = make_field(parameters, "Trigger parameters (or corresponding file name)") = JPARSER::initialised();
170 zap['O'] = make_field(triggeredEventsOnly, "optionally write only triggered events.");
171 zap['P'] = make_field(pmtParameters, "PMT simulation data (or corresponding file name)") = JPARSER::initialised();
172 zap['B'] = make_field(rates_Hz, "background rates [Hz]") = JPARSER::initialised();
173 zap['X'] = make_field(mixed_Hz, "mixed-L1/L0 background rates [Hz]") = JPARSER::initialised();
174 zap['s'] = make_field(sigma_ns, "intrinsic time smearing of K40 coincidences [ns]") = JK40DefaultSimulatorInterface::getSigma();
175 zap['S'] = make_field(seed, "seed") = 0;
176 zap['d'] = make_field(debug, "debug") = 0;
177
178 zap(argc, argv);
179 }
180 catch(const exception &error) {
181 FATAL(error.what() << endl);
182 }
183
184 seed.set(gRandom);
185
186 JK40DefaultSimulatorInterface::setSigma(sigma_ns);
187
188 setDAQLongprint(debug >= JEEP::debug_t);
189
190 if (detectorFileB == "") {
191 detectorFileB = detectorFileA;
192 }
193
194
195 JDetector detectorA;
196 JDetector detectorB;
197
198 try {
199 load(detectorFileA, detectorA);
200 load(detectorFileB, detectorB);
201 }
202 catch(const JException& error) {
203 FATAL(error);
204 }
205
206 const double TA_ns = getT0(detectorA);
207 const double TB_ns = getT0(detectorB);
208
209 JPMTParametersMap::Throw(true);
210
211 if (!pmtParameters.is_valid()) {
212 FATAL("Invalid PMT parameters " << pmtParameters << endl);
213 }
214
215 if (pmtParameters.getQE() != 1.0) {
216
217 WARNING("Correct background rates with global efficiency " << pmtParameters.getQE() << endl);
218
219 rates_Hz.correct(pmtParameters.getQE());
220 }
221
222 Head header;
223
224 try {
225 header = getHeader(inputFile);
226 }
227 catch(const JException& error) {
228 FATAL(error);
229 }
230
231 const JModuleRouter moduleRouter(detectorB);
232 JDetectorSimulator simbad (detectorA);
233 JSummaryRouter summaryRouter;
234
235 JTreeScanner<JDAQSummaryslice, JDAQEvaluator> scanner; // run-by-run
236
237 set<int> fs; // frame indices of summary data already written
238
239 JDAQUTCExtended UTC = JDAQUTCExtended::getInstance(); // start of run
240
241 if (runbyrun != "") {
242
243 NOTICE("Using run-by-run: " << runbyrun << endl);
244
245 scanner.configure(runbyrun);
246
247 if (!scanner.empty()) {
248 UTC = scanner.begin()->getTimesliceStart();
249 } else {
250 FATAL("Run-by-run simulation misses summary data." << endl);
251 }
252
253 try {
254 simbad.reset(new JK40RunByRunSimulator(summaryRouter, rates_Hz, mixed_Hz));
255 simbad.reset(new JPMTRunByRunSimulator(summaryRouter, pmtParameters, detectorA));
256 simbad.reset(new JCLBRunByRunSimulator(summaryRouter));
257 }
258 catch(const JException& error) {
259 FATAL(error.what() << endl);
260 }
261
262 try {
263
264 parameters = getTriggerParameters(runbyrun);
265
266 NOTICE("Set trigger parameters from run-by-run input." << endl);
267 }
268 catch(const JException& error) {
269 WARNING("No trigger parameters from run-by-run input;\nrun with default/user input." << endl);
270 }
271
272 // set live time
273
274 JHead buffer(header);
275
276 buffer.DAQ.livetime_s = getTimeDuration(getUTCTimeRange(scanner));
277 buffer.push(&JHead::DAQ);
278
279 copy(buffer, header);
280
281 } else {
282
283 NOTICE("Using fixed rates [Hz]: " << rates_Hz << endl);
284
285 try {
286 simbad.reset(new JK40DefaultSimulator(rates_Hz, mixed_Hz));
287 simbad.reset(new JPMTDefaultSimulator(pmtParameters, detectorA));
288 simbad.reset(new JCLBDefaultSimulator());
289 }
290 catch(const JException& error) {
291 FATAL(error.what() << endl);
292 }
293 }
294
295 // detector
296
297 if (parameters.disableHighRateVeto) {
298
299 NOTICE("Disabling high-rate veto of all PMTs." << endl);
300
301 detectorB.setPMTStatus(HIGH_RATE_VETO_DISABLE);
302 }
303
304 parameters.set(getMaximalDistance(detectorB));
305
306 DEBUG("Trigger:" << endl << parameters << endl);
307 DEBUG("PMT parameters:" << endl << pmtParameters << endl);
308
309 const double Tmax = max(getMaximalTime(detectorA),
310 getMaximalTime(detectorB));
311
312 const JTimeRange period(-(Tmax + parameters.TMaxLocal_ns),
313 +(Tmax + parameters.TMaxLocal_ns));
314
315 const JTimeRange frame_time(0.0, getFrameTime());
316
317 // hit type should be identical to that in JDataFilter.cc
318
319 typedef double hit_type;
320
321 typedef JSuperFrame1D<hit_type> JSuperFrame1D_t;
322 typedef JSuperFrame2D<hit_type> JSuperFrame2D_t;
323 typedef JTimeslice <hit_type> JTimeslice_t;
324 typedef JBuildL1 <hit_type> JBuildL1_t;
325 typedef JBuildL2 <hit_type> JBuildL2_t;
326
327 const JBuildL1_t buildL1(parameters);
328 const JBuildL2_t buildL2(parameters.L2);
329 const JBuildL2_t buildSN(parameters.SN);
330
331 JTimesliceRouter timesliceRouter(parameters.numberOfBins);
332
333 const JTrigger3DMuon trigger3DMuon (parameters);
334 const JTrigger3DShower trigger3DShower(parameters);
335 const JTriggerMXShower triggerMXShower(parameters, detectorB);
336
337
338 TH1D h1("Trigger bits", NULL, NUMBER_OF_TRIGGER_BITS, -0.5, NUMBER_OF_TRIGGER_BITS - 0.5);
339
340
341 outputFile.open();
342
343 if (!outputFile.is_open()) {
344 FATAL("Error opening file " << outputFile << endl);
345 }
346
347 outputFile.put(*gRandom);
348 outputFile.put(JMeta(argc, argv));
349 outputFile.put(header);
350 outputFile.put(parameters);
351
352 JHead head;
353 int trigger_counter = 0;
354
355 for (string file_name; inputFile.hasNext(); ) {
356
357 STATUS("event: " << setw(10) << inputFile.getCounter() << '\r'); DEBUG(endl);
358
359 Evt* event = inputFile.next();
360
361 DEBUG(*event << endl);
362
363 if (file_name != inputFile.getFilename()) {
364
365 file_name = inputFile.getFilename();
366 head = getHeader(file_name);
367 }
368
369 event->mc_run_id = head.start_run.run_id;
370
371 if (run == -1) {
372 run = event->mc_run_id;
373 }
374
375 bool trigger = false;
376 int frame_index = inputFile.getCounter() + 1; // event / slice
377 double t1 = gRandom->Rndm() * getFrameTime(); // time of event within time slice
378
379 JDAQUTCExtended utc(UTC.getTimeNanoSecond() + getTimeOfFrame(frame_index));
380
381 if (runbyrun != "") {
382
383 if (event->mc_event_time == TTimeStamp(0)) {
384 FATAL("Monte Carlo event time undefined." << endl);
385 }
386
387 utc = getDAQUTCExtended(event->mc_event_time);
388 t1 = getTimeDifference(summaryRouter.getTimesliceStart(), utc) * 1.0e9;
389
390 if (!frame_time(t1)) {
391
392 const Long64_t index = scanner.find(utc);
393 JDAQSummaryslice* p = scanner.getEntry(index);
394
395 if (getTimeDifference(p->getTimesliceStart(), utc) < 0.0 && index != 0) {
396 p = scanner.getEntry(index - 1);
397 }
398
399 summaryRouter.update(p);
400 summaryRouter.correct(dynamic_cast<const JPMTDefaultSimulatorInterface&>(simbad.getPMTSimulator()));
401 }
402
403 t1 = getTimeDifference(summaryRouter.getTimesliceStart(), utc) * 1.0e9;
404
405 run = summaryRouter.getRunNumber();
406 frame_index = summaryRouter.getFrameIndex();
407 utc = summaryRouter.getTimesliceStart();
408
409 DEBUG("event: "
410 << FIXED(15,3) << event->mc_event_time.AsDouble() << " [s] "
411 << utc << ' '
412 << FIXED(12,0) << t1 << " [ns] "
413 << frame_time(t1) << endl);
414 }
415
416 if (!event->mc_hits.empty() && frame_time(t1)) {
417
418 JTimeRange timeRange = getTimeRange(*event, period);
419
420 double t0 = 0.5 * (timeRange.getLowerLimit() + timeRange.getUpperLimit());
421
422 DEBUG("Start time: " << FIXED(12,2) << t0 << ' ' << FIXED(12,2) << t1 << ' ' << FIXED(9,2) << TA_ns << endl);
423
424 t1 += TA_ns; // add time offset of detector
425
426 event->mc_t = getTimeOfFrame(frame_index) + t1 - t0; // set to time since start of data taking run
427
428 timeRange.add(event->mc_t);
429 timeRange.add(period);
430
431 const JDAQChronometer chronometer(detectorB.getID(), run, frame_index, utc);
432
433 const JEventTimeslice timeslice(chronometer, simbad, *event, period);
434
435 DEBUG(timeslice << endl);
436
437
438 timesliceRouter.configure(timeslice);
439
440 JTimeslice_t timesliceL0(timeslice.getDAQChronometer());
441 JTimeslice_t timesliceL1(timeslice.getDAQChronometer());
442 JTimeslice_t timesliceL2(timeslice.getDAQChronometer());
443 JTimeslice_t timesliceSN(timeslice.getDAQChronometer());
444
445 for (JDAQTimeslice::const_iterator super_frame = timeslice.begin(); super_frame != timeslice.end(); ++super_frame) {
446
447 if (moduleRouter.hasModule(super_frame->getModuleID())) {
448
449 // calibration
450
451 const JModule& module = moduleRouter.getModule(super_frame->getModuleID());
452 const JSuperFrame2D_t& buffer = JSuperFrame2D_t::demultiplex(*super_frame, module);
453
454 // L0
455
456 timesliceL0.push_back(JSuperFrame1D_t(buffer));
457
458 DEBUG("L0 " << setw(8) << timesliceL0.rbegin()->getModuleID() << ' ' << setw(8) << timesliceL0.rbegin()->size() << LAMBDA([ps = timesliceL0.rbegin()](ostream& out) { for (const auto& i : *ps) { out << " " << i; } }) << endl);
459
460 // L1
461
462 timesliceL1.push_back(JSuperFrame1D_t(super_frame->getDAQChronometer(),
463 super_frame->getModuleIdentifier(),
464 module.getPosition()));
465
466 buildL1(*timesliceL0.rbegin(), back_inserter(*timesliceL1.rbegin()));
467
468 DEBUG("L1 " << setw(8) << timesliceL1.rbegin()->getModuleID() << ' ' << setw(8) << timesliceL1.rbegin()->size() << LAMBDA([ps = timesliceL1.rbegin()](ostream& out) { for (const auto& i : *ps) { out << " " << i; } }) << endl);
469
470 // L2
471
472 timesliceL2.push_back(JSuperFrame1D_t(super_frame->getDAQChronometer(),
473 super_frame->getModuleIdentifier(),
474 module.getPosition()));
475
476 buildL2(buffer, *timesliceL1.rbegin(), back_inserter(*timesliceL2.rbegin()));
477
478 DEBUG("L2 " << setw(8) << timesliceL2.rbegin()->getModuleID() << ' ' << setw(8) << timesliceL2.rbegin()->size() << LAMBDA([ps = timesliceL2.rbegin()](ostream& out) { for (const auto& i : *ps) { out << " " << i; } }) << endl);
479
480 // SN
481 {
482 JTimeslice_t::value_type tv(super_frame->getDAQChronometer(),
483 super_frame->getModuleIdentifier(),
484 module.getPosition());
485
486 buildSN(buffer, *timesliceL1.rbegin(), back_inserter(tv));
487
488 if (!tv.empty()) {
489
490 timesliceSN.push_back(tv);
491
492 DEBUG("SN " << setw(8) << timesliceSN.rbegin()->getModuleID() << ' ' << setw(8) << timesliceSN.rbegin()->size() << endl);
493 DEBUG("L2 " << setw(8) << timesliceL2.rbegin()->getModuleID() << ' ' << setw(8) << timesliceL2.rbegin()->size() << LAMBDA([ps = timesliceL2.rbegin()](ostream& out) { for (const auto& i : *ps) { out << " " << i; } }) << endl);
494 }
495 }
496 }
497 }
498
499
500 // Trigger
501
502 JTriggerInput trigger_input(timesliceL2);
503 JTriggerOutput trigger_output;
504
505 trigger3DMuon (trigger_input, back_inserter(trigger_output));
506 trigger3DShower(trigger_input, back_inserter(trigger_output));
507 triggerMXShower(trigger_input, timesliceL0, back_inserter(trigger_output));
508
509 trigger_output.merge(JEventOverlap(parameters.TMaxEvent_ns));
510
511 for (JTriggerOutput::const_iterator to = trigger_output.begin(); to != trigger_output.end(); ++to) {
512
513 for (int i = 0; i != h1.GetNbinsX(); ++i) {
514 if (to->hasTriggerBit(i)) {
515 h1.Fill((double) i);
516 }
517 }
518
519 JTimeRange eventTime = getTimeRange(*to).add(getTimeOfRTS(*to));
520
521 eventTime.add(TA_ns);
522 eventTime.sub(TB_ns);
523
524 if (debug >= status_t) {
525 cout << "Event time: "
526 << to->getFrameIndex() << ' '
527 << eventTime << ' '
528 << timeRange << ' '
529 << (timeRange.overlap(eventTime) ? "Y" : "N") << endl;
530 }
531
532 if (timeRange.overlap(eventTime)) {
533
534 JTriggeredEvent tev(*to,
535 timesliceRouter,
536 moduleRouter,
537 parameters.TMaxLocal_ns,
538 getTimeRange(parameters));
539
540 tev.setCounter(trigger_counter); // set the event counter to the index of the Monte Carlo event in the output.
541
542 outputFile.put(tev);
543
544 trigger = true;
545 }
546 }
547
548
549 if (!triggeredEventsOnly || trigger) {
550
551 if (parameters.writeL0()) {
552 outputFile.put(timeslice);
553 }
554
555 if (parameters.writeL1()) {
556 outputFile.put(JTimesliceL1<JDAQTimesliceL1>(timesliceL1, timesliceRouter, moduleRouter, parameters.TMaxLocal_ns));
557 }
558
559 if (parameters.writeL2()) {
560 outputFile.put(JTimesliceL1<JDAQTimesliceL2>(timesliceL2, timesliceRouter, moduleRouter, parameters.L2.TMaxLocal_ns));
561 }
562
563 if (parameters.writeSN()) {
564 outputFile.put(JTimesliceL1<JDAQTimesliceSN>(timesliceSN, timesliceRouter, moduleRouter, parameters.SN.TMaxLocal_ns));
565 }
566
567 if (parameters.writeSummary()) {
568
569 if (fs.count(chronometer.getFrameIndex()) == 0) {
570
571 outputFile.put(JSummaryslice(chronometer,simbad));
572
573 fs.insert(chronometer.getFrameIndex());
574 }
575 }
576 }
577 }
578
579 if (!triggeredEventsOnly || trigger) {
580
581 outputFile.put(*event);
582
583 ++trigger_counter;
584 }
585 }
586 STATUS(endl);
587
588 {
589 JMultipleFileScanner<JTYPELIST<JMetaTypes_t, MultiHead>::typelist> io(inputFile);
590
591 io >> outputFile;
592 }
593
594 outputFile.put(h1);
595 outputFile.put(*gRandom);
596 outputFile.close();
597}
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:74
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
JAANET::JHead
Monte Carlo run header.
Definition JHead.hh:1236
JAANET::JHead::start_run
JAANET::start_run start_run
Definition JHead.hh:1601
JAANET::JHead::DAQ
JAANET::DAQ DAQ
Definition JHead.hh:1625
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::JDetector::setPMTStatus
void setPMTStatus(const int bit)
Set status of all PMTs.
Definition JDetector.hh:207
JDETECTOR::JK40DefaultSimulatorInterface::getSigma
static double getSigma()
Get intrinsic time smearing of K40 coincidences.
Definition JK40DefaultSimulatorInterface.hh:378
JDETECTOR::JK40DefaultSimulatorInterface::setSigma
static void setSigma(const double sigma)
Set intrinsic time smearing of K40 coincidences.
Definition JK40DefaultSimulatorInterface.hh:389
JDETECTOR::JK40DefaultSimulator
Default implementation of the simulation of K40 background.
Definition JK40DefaultSimulator.hh:35
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
JDETECTOR::JPMTParametersMap
Auxiliary class for map of PMT parameters.
Definition JPMTParametersMap.hh:93
JDETECTOR::JPMTParametersMap::getQE
double getQE(const JPMTIdentifier &id) const
Get QE of given PMT.
Definition JPMTParametersMap.hh:164
JDETECTOR::JPMTParametersMap::is_valid
bool is_valid() const
Check validity of PMT parameters.
Definition JPMTParametersMap.hh:255
JGEOMETRY3D::JPosition3D::getPosition
const JPosition3D & getPosition() const
Get position.
Definition JPosition3D.hh:130
JLANG::JException
General exception.
Definition JException.hh:24
JLANG::JException::what
virtual const char * what() const override
Get error message.
Definition JException.hh:64
JLANG::JObjectID::getID
int getID() const
Get identifier.
Definition JObjectID.hh:50
JLANG::JThrow::Throw
static void Throw(const bool option)
Enable/disable throw option.
Definition JThrow.hh:37
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
JSUPPORT::JMultipleFileScanner::hasNext
virtual bool hasNext() override
Check availability of next element.
Definition JMultipleFileScanner.hh:409
JSUPPORT::JMultipleFileScanner::getCounter
counter_type getCounter() const
Get counter.
Definition JMultipleFileScanner.hh:382
JSUPPORT::JMultipleFileScanner::next
virtual const pointer_type & next() override
Get next element.
Definition JMultipleFileScanner.hh:457
JSUPPORT::JMultipleFileScanner::getFilename
const std::string & getFilename() const
Get current file name.
Definition JMultipleFileScanner.hh:371
JSUPPORT::JSummaryRouter
Router for fast addressing of summary data in KM3NETDAQ::JDAQSummaryslice data structure as a functio...
Definition JSupport/JSummaryRouter.hh:38
JSUPPORT::JSummaryRouter::update
void update(const JDAQSummaryslice *ps)
Update router.
Definition JSupport/JSummaryRouter.hh:53
JSUPPORT::JTreeScanner
Template definition for direct access of elements in ROOT TChain.
Definition JTreeScanner.hh:698
JTOOLS::JRange::overlap
bool overlap(const range_type &range) const
Test overlap with given range.
Definition JRange.hh:382
JTOOLS::JRange::getLowerLimit
T getLowerLimit() const
Get lower limit.
Definition JRange.hh:202
JTOOLS::JRange::sub
range_type & sub(argument_type x)
Subtract offset.
Definition JRange.hh:460
JTOOLS::JRange::add
range_type & add(argument_type x)
Add offset.
Definition JRange.hh:446
JTOOLS::JRange::getUpperLimit
T getUpperLimit() const
Get upper limit.
Definition JRange.hh:213
JTRIGGER::JCLBRunByRunSimulator
CLB simulation based on run-by-run information.
Definition JCLBRunByRunSimulator.hh:34
JTRIGGER::JK40RunByRunSimulator
K40 simulation based on run-by-run information.
Definition JK40RunByRunSimulator.hh:35
JTRIGGER::JPMTRunByRunSimulator
PMT simulation based on run-by-run information.
Definition JPMTRunByRunSimulator.hh:37
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::JTriggeredEvent
Auxiliary class to build KM3NETDAQ::JDAQEvent for a triggered event.
Definition JTriggeredEvent.hh:50
KM3NETDAQ::JDAQChronometer::getRunNumber
int getRunNumber() const
Get run number.
Definition JDAQChronometer.hh:121
KM3NETDAQ::JDAQChronometer::getTimesliceStart
JDAQUTCExtended getTimesliceStart() const
Get start of timeslice.
Definition JDAQChronometer.hh:144
KM3NETDAQ::JDAQChronometer::getFrameIndex
int getFrameIndex() const
Get frame index.
Definition JDAQChronometer.hh:132
KM3NETDAQ::JDAQUTCExtended
Data structure for UTC time.
Definition JDAQUTCExtended.hh:23
KM3NETDAQ::JDAQUTCExtended::getInstance
static const JDAQUTCExtended & getInstance()
Get arbitrary offset (e.g.
Definition JDAQUTCExtended.hh:192
KM3NETDAQ::JDAQUTCExtended::getTimeNanoSecond
double getTimeNanoSecond() const
Get time (limited to 16 ns cycles).
Definition JDAQUTCExtended.hh:109
JAANET::getTimeRange
JTimeRange getTimeRange(const Evt &event)
Get time range (i.e. time between earliest and latest hit) of Monte Carlo event.
Definition JAAnetToolkit.hh:140
JAANET::copy
void copy(const Head &from, JHead &to)
Copy header from from to to.
Definition JHead.cc:163
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::status_t
@ status_t
status
Definition JMessage.hh:30
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:48
JSUPPORT::getTimeDuration
double getTimeDuration(const JDAQUTCTimeRange &utc)
Get time duration of given UTC time range.
Definition JSupportToolkit.hh:60
JSUPPORT::getUTCTimeRange
JDAQUTCTimeRange getUTCTimeRange()
Get UTC time range.
Definition JSupportToolkit.hh:145
JSUPPORT::getTriggerParameters
JTriggerParameters getTriggerParameters(const JMultipleFileScanner_t &file_list)
Get trigger parameters.
Definition JTriggerParametersSupportkit.hh:148
JSUPPORT::getHeader
Head getHeader(const JMultipleFileScanner_t &file_list)
Get Monte Carlo header.
Definition JMonteCarloFileSupportkit.hh:547
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::getTimeOfFrame
double getTimeOfFrame(const int frame_index)
Get start time of frame in ns since start of run for a given frame index.
Definition JDAQClock.hh:185
KM3NETDAQ::getMaximalTime
double getMaximalTime(const double R_Hz)
Get maximal time for given rate.
Definition JHighRateVeto.hh:35
KM3NETDAQ::getTimeDifference
double getTimeDifference(const JDAQChronometer &first, const JDAQChronometer &second)
Get time difference between two chronometers.
Definition JDAQChronometer.hh:267
KM3NETDAQ::getDAQUTCExtended
JDAQUTCExtended getDAQUTCExtended(const TTimeStamp &t0, const double t1=0.0)
Get DAQ UTC time.
Definition JDAQToolkit.hh:26
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
Evt
The Evt class respresent a Monte Carlo (MC) event as well as an offline event.
Definition Evt.hh:21
FIXED
Auxiliary data structure for floating point format specification.
Definition JManip.hh:448
Head
The Head class reflects the header of Monte-Carlo event files, which consists of keys (also referred ...
Definition Head.hh:65
JAANET::start_run::run_id
int run_id
MC run number.
Definition JHead.hh:143
JACOUSTICS::JEventOverlap
Match of two events considering overlap in time and position.
Definition JAcoustics/JEventOverlap.hh:20
JACOUSTICS::hit_type
Transmission with position.
Definition JBillabong.cc:70
JLANG::JTypeList
Type list.
Definition JTypeList.hh:23
JMATH::JRandom
Template definition of random value generator.
Definition JMath/JRandom.hh:26
JPARSER::initialised
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Definition JParser.hh:68
JPHYSICS::JK40Rates
Auxiliary class for K40 rates.
Definition JK40Rates.hh:41
JPHYSICS::JK40Rates::correct
void correct(const double QE)
Correct rates for global efficiency,.
Definition JK40Rates.hh:130
JPHYSICS::JMixedK40Rates
Auxiliary class for mixed-L1/L0 K40 rates.
Definition JK40Rates.hh:229
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
KM3NETDAQ::JEventTimeslice
Timeslice with Monte Carlo event.
Definition JEventTimeslice.hh:39
KM3NETDAQ::JSummaryslice
Auxiliary class to create summary data.
Definition JSummaryslice.hh:41
LAMBDA
Auxiliary data structure to convert (lambda) function to printable object.
Definition JManip.hh:726
Generated by doxygen 1.12.0