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

Auxiliary program to re-trigger KM3NETDAQ::JDAQEvent data. More...

#include <string>
#include <iostream>
#include <iomanip>
#include "km3net-dataformat/offline/Head.hh"
#include "km3net-dataformat/offline/Evt.hh"
#include "km3net-dataformat/offline/Hit.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 "JDAQ/JDAQEvaluator.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 "JLang/JPipe.hh"
#include "JROOT/JROOTClassSelector.hh"
#include "JSupport/JSingleFileScanner.hh"
#include "JSupport/JTreeScanner.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-trigger KM3NETDAQ::JDAQEvent data.

Author
mdejong

Definition in file JTriggerReprocessor.cc.

Function Documentation

int main ( int  argc,
char **  argv 
)

Definition at line 60 of file JTriggerReprocessor.cc.

61 {
62  using namespace std;
63  using namespace JPP;
64  using namespace KM3NETDAQ;
65 
66  typedef JAllTypes_t typelist;
67 
68  JSingleFileScanner<JDAQEvent> inputFile;
69  JFileRecorder<typelist> outputFile;
70  JLimit_t& numberOfEvents = inputFile.getLimit();
71  string detectorFile;
72  JTriggerParameters parameters;
73  bool reuse_parameters;
74  bool snapshot;
75  JROOTClassSelection selection = getROOTClassSelection<typelist>();
76  int debug;
77 
78  try {
79 
80  JParser<> zap("Auxiliary program to re-trigger event data.");
81 
82  zap['f'] = make_field(inputFile, "input file.");
83  zap['o'] = make_field(outputFile, "output file.") = "trigger_reprocessor.root";
84  zap['n'] = make_field(numberOfEvents) = JLimit::max();
85  zap['a'] = make_field(detectorFile, "detector file.");
86  zap['@'] = make_field(parameters, "trigger parameters") = JPARSER::initialised();
87  zap['U'] = make_field(reuse_parameters, "reuse trigger parameters from input file.");
88  zap['S'] = make_field(snapshot, "use snapshot hits instead of triggered hits.");
89  zap['C'] = make_field(selection, "selection of data types for output.") = JPARSER::initialised();
90  zap['d'] = make_field(debug, "debug flag.") = 1;
91 
92  zap(argc, argv);
93  }
94  catch(const exception &error) {
95  FATAL(error.what() << endl);
96  }
97 
98 
99  cout.tie(&cerr);
100 
101  setDAQLongprint(debug >= JEEP::debug_t);
102 
103  JDetector detector;
104 
105  try {
106  load(detectorFile, detector);
107  }
108  catch(const JException& error) {
109  FATAL(error);
110  }
111 
112  if (reuse_parameters) {
113 
114  try {
115 
116  parameters = getTriggerParameters(inputFile);
117 
118  NOTICE("Set trigger parameters from input." << endl);
119  }
120  catch(const JException& error) {
121  FATAL("No trigger parameters from input." << endl);
122  }
123  }
124 
125  parameters.set(getMaximalDistance(detector));
126 
127  parameters.triggerNB.write.prescale = 1;
128 
129  DEBUG("Frame time [ms] " << getFrameTime() * 1e-6 << endl);
130  DEBUG("Reset time [ms] " << getRTS() * 1e-6 << endl);
131  DEBUG("Trigger" << endl << parameters << endl);
132 
133  const JModuleRouter moduleRouter(detector);
134 
135  if (parameters.writeSummary()) { WARNING("Discard writeSummary option during reprocesing of data." << endl); }
136  if (parameters.writeL1()) { WARNING("Discard writeL1 option during reprocesing of data." << endl); }
137  if (parameters.writeL2()) { WARNING("Discard writeL2 option during reprocesing of data." << endl); }
138  if (parameters.writeSN()) { WARNING("Discard writeSN option during reprocesing of data." << endl); }
139 
140  //typedef JHit hit_type;
141  //typedef int hit_type;
142  typedef double hit_type;
143 
144  typedef JSuperFrame1D<hit_type> JSuperFrame1D_t;
145  typedef JSuperFrame2D<hit_type> JSuperFrame2D_t;
146  typedef JTimeslice <hit_type> JTimeslice_t;
147  typedef JBuildL1 <hit_type> JBuildL1_t;
148  typedef JBuildL2 <hit_type> JBuildL2_t;
149 
150  const JBuildL1_t buildL1(parameters);
151  const JBuildL2_t buildL2(parameters.L2);
152  const JBuildL2_t buildSN(parameters.SN);
153  const JBuildL2_t buildNB(parameters.NB);
154 
155  JTimesliceRouter timesliceRouter(parameters.numberOfBins);
156 
157  const JTriggerNB triggerNB (parameters);
158  const JTrigger3DMuon trigger3DMuon (parameters);
159  const JTrigger3DShower trigger3DShower(parameters);
160  const JTriggerMXShower triggerMXShower(parameters, detector);
161 
162 
163  outputFile.open();
164 
165  if (!outputFile.is_open()) {
166  FATAL("Error opening file " << outputFile << endl);
167  }
168 
169  outputFile.put(JMeta(argc, argv));
170  outputFile.put(parameters);
171 
172 
173  JTreeScanner<JDAQSummaryslice, JDAQEvaluator>::debug = debug;
174 
175  JTreeScanner<JDAQSummaryslice, JDAQEvaluator> scan(inputFile);
176 
177  if (scan.getEntries() == 0) {
178  FATAL("No summary data." << endl);
179  }
180 
181  while (inputFile.hasNext()) {
182 
183  STATUS("event: " << setw(10) << inputFile.getCounter() << '\r'); DEBUG(endl);
184 
185  JDAQEvent* evt = inputFile.next();
186  Long64_t index = scan.find(*evt);
187  JDAQSummaryslice* summary = scan.getEntry(index);
188 
189  if (evt->getFrameIndex() != summary->getFrameIndex()) {
190  DEBUG("Frame indices " << evt->getFrameIndex() << " != " << summary->getFrameIndex() << endl);
191  }
192 
193  JDAQTimeslice timeslice(*evt, *summary, snapshot);
194 
195  DEBUG(timeslice << endl);
196 
197  timesliceRouter.configure(timeslice);
198 
199 
200  JTimeslice_t timesliceL0(timeslice.getDAQChronometer());
201  JTimeslice_t timesliceL1(timeslice.getDAQChronometer());
202  JTimeslice_t timesliceL2(timeslice.getDAQChronometer());
203  JTimeslice_t timesliceSN(timeslice.getDAQChronometer());
204  JTimeslice_t timesliceNB(timeslice.getDAQChronometer());
205 
206  for (JDAQTimeslice::const_iterator super_frame = timeslice.begin(); super_frame != timeslice.end(); ++super_frame) {
207 
208  if (moduleRouter.hasModule(super_frame->getModuleID())) {
209 
210  // calibration
211 
212  const JModule& module = moduleRouter.getModule(super_frame->getModuleID());
213  JSuperFrame2D_t& buffer = JSuperFrame2D_t::demultiplex(*super_frame, module);
214 
215  // Apply high-rate veto
216 
217  buffer.applyHighRateVeto(parameters.highRateVeto_Hz);
218 
219  // L0
220 
221  timesliceL0.push_back(JSuperFrame1D_t(buffer));
222 
223  // Nano-beacon trigger
224 
225  if (parameters.triggerNB.enabled) {
226 
227  JSuperFrame2D_t::iterator __end = partition(buffer.begin(), buffer.end(), parameters.triggerNB.pmts);
228 
229  if (buffer.begin() != __end) {
230 
231  timesliceNB.push_back(JSuperFrame1D_t(super_frame->getDAQChronometer(),
232  super_frame->getModuleIdentifier(),
233  module.getPosition()));
234 
235 
236  JSuperFrame1D_t zbuf;
237 
238  buildL1(buffer.begin(), __end , back_inserter(zbuf));
239 
240 
241  buildNB(buffer.begin() , __end, zbuf, back_inserter(*timesliceNB.rbegin()));
242  }
243  }
244 
245  // L1
246 
247  timesliceL1.push_back(JSuperFrame1D_t(super_frame->getDAQChronometer(),
248  super_frame->getModuleIdentifier(),
249  module.getPosition()));
250 
251  buildL1(*timesliceL0.rbegin(), back_inserter(*timesliceL1.rbegin()));
252 
253  // L2
254 
255  timesliceL2.push_back(JSuperFrame1D_t(super_frame->getDAQChronometer(),
256  super_frame->getModuleIdentifier(),
257  module.getPosition()));
258 
259  buildL2(buffer, *timesliceL1.rbegin(), back_inserter(*timesliceL2.rbegin()));
260 
261  // SN
262 
263  timesliceSN.push_back(JSuperFrame1D_t(super_frame->getDAQChronometer(),
264  super_frame->getModuleIdentifier(),
265  module.getPosition()));
266 
267  buildSN(buffer, *timesliceL1.rbegin(), back_inserter(*timesliceSN.rbegin()));
268  }
269  }
270 
271 
272  // Trigger
273 
274  if (parameters.triggerNB.enabled) {
275 
276  const JTriggerInput trigger_input(timesliceNB);
277 
278  for (JTriggerInput::const_iterator hit = trigger_input.begin(); hit != trigger_input.end(); ++hit) {
279 
280  if (parameters.triggerNB.write()) {
281 
282  JTriggeredEvent tev(timesliceNB.getDAQChronometer(),
283  getTriggerMask(triggerNB.getTriggerBit()),
284  *hit,
285  timesliceRouter,
286  moduleRouter,
287  parameters.TMaxLocal_ns,
288  parameters.triggerNB.DMax_m,
289  getTimeRange(parameters.triggerNB));
290 
291  outputFile.put(tev);
292  }
293  }
294  }
295 
296  JTriggerInput trigger_input(timesliceL2);
297  JTriggerOutput trigger_output;
298 
299  trigger3DMuon (trigger_input, back_inserter(trigger_output));
300  trigger3DShower(trigger_input, back_inserter(trigger_output));
301  triggerMXShower(trigger_input, timesliceL0, back_inserter(trigger_output));
302 
303  trigger_output.merge(JEventOverlap(parameters.TMaxEvent_ns));
304 
305  DEBUG("Number of triggers: " << trigger_output.size() << endl);
306 
307  for (JTriggerOutput::const_iterator event = trigger_output.begin(); event != trigger_output.end(); ++event) {
308 
309  JTriggeredEvent tev(*event,
310  timesliceRouter,
311  moduleRouter,
312  parameters.TMaxLocal_ns,
313  getTimeRange(parameters));
314 
315  tev.setCounter(evt->getCounter());
316 
317  outputFile.put(tev);
318  }
319  }
320  STATUS(endl);
321 
322  JSingleFileScanner<typelist> io(inputFile);
323 
324  selection.remove<JDAQEvent>();
325 
326  io | JValve<typelist>(selection) | outputFile;
327 
328  outputFile.close();
329 }
JACOUSTICS::hit_type
Acoustic hit.
Definition: JAcousticsTriggerProcessor.cc:59
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:41
JPARSER::JParser
Utility class to parse command line options.
Definition: JParser.hh:1500
JLANG::JException
General exception.
Definition: JException.hh:23
WARNING
#define WARNING(A)
Definition: JMessage.hh:65
KM3NETDAQ::getTriggerMask
JTriggerMask_t getTriggerMask(const unsigned int bit)
Convert trigger bit to trigger mask.
Definition: JDAQTriggerMask.hh:36
JEEP::debug_t
debug
Definition: JMessage.hh:29
JDETECTOR::JModule
Data structure for a composite optical module.
Definition: JModule.hh:68
JTRIGGER::JTimesliceRouter
Router for fast addressing of hits in KM3NETDAQ::JDAQTimeslice data structure as a function of the op...
Definition: JTimesliceRouter.hh:91
STATUS
#define STATUS(A)
Definition: JMessage.hh:63
JDETECTOR::JDetector
Detector data structure.
Definition: JDetector.hh:89
JDETECTOR::JModuleRouter
Router for direct addressing of module data in detector data structure.
Definition: JModuleRouter.hh:36
JROOT::JROOTClassSelection
Auxiliary class for ROOT class selection.
Definition: JROOTClassSelector.hh:91
parameters
*fatal Wrong number of arguments esac JCookie sh typeset Z DETECTOR typeset Z SOURCE_RUN typeset Z TARGET_RUN set_variable PARAMETERS_FILE $WORKDIR parameters
Definition: diff-Tuna.sh:38
JPARSER::initialised
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Definition: JParser.hh:66
outputFile
string outputFile
Definition: JDAQTimesliceSelector.cc:37
JSUPPORT::JTreeScanner
Template definition for direct access of elements in ROOT TChain.
Definition: JTreeScanner.hh:91
JLANG::JValve
Auxiliary class for selection of data type.
Definition: JObjectIterator.hh:34
KM3NETDAQ::JDAQChronometer::getFrameIndex
int getFrameIndex() const
Get frame index.
Definition: JDAQChronometer.hh:132
JTRIGGER::JSuperFrame1D
1-dimensional frame with time calibrated data from one optical module.
Definition: JSuperFrame1D.hh:37
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:121
JLANG::JTypeList
Type list.
Definition: JTypeList.hh:22
KM3NETDAQ::setDAQLongprint
void setDAQLongprint(const bool option)
Set DAQ print option.
Definition: JDAQPrint.hh:28
JSUPPORT::JLimit
Auxiliary class for defining the range of iterations of objects.
Definition: JLimit.hh:41
JTRIGGER::JBuildL2
Template L2 builder.
Definition: JBuildL2.hh:45
JDETECTOR::getMaximalDistance
double getMaximalDistance(const JDetector &detector)
Get maximal distance between modules in detector.
Definition: JDetectorToolkit.hh:78
JAANET::detector
Detector file.
Definition: JHead.hh:224
JACOUSTICS::JTriggerOutput::merge
void merge(const JMatch_t &match)
Merge events.
Definition: JAcoustics/JTriggerOutput.hh:40
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition: JParser.hh:1961
KM3NETDAQ::getFrameTime
double getFrameTime()
Get frame time duration.
Definition: JDAQClock.hh:162
NOTICE
#define NOTICE(A)
Definition: JMessage.hh:64
KM3NETDAQ::JDAQTimeslice
Data time slice.
Definition: JDAQTimeslice.hh:30
JTRIGGER::JTriggeredEvent
Auxiliary class to build JDAQEvent for a triggered event.
Definition: JTriggeredEvent.hh:46
debug
int debug
debug level
Definition: JSirene.cc:63
JGEOMETRY3D::JPosition3D::getPosition
const JPosition3D & getPosition() const
Get position.
Definition: JPosition3D.hh:130
JTRIGGER::JTriggerNB
Nano-beacon trigger.
Definition: JTriggerNB.hh:19
JTRIGGER::JBuildL1
Template L1 hit builder.
Definition: JBuildL1.hh:85
FATAL
#define FATAL(A)
Definition: JMessage.hh:67
JDETECTOR::load
void load(const std::string &file_name, JDetector &detector)
Load detector from input file.
Definition: JDetectorToolkit.hh:478
KM3NETDAQ::getRTS
double getRTS()
Get TDC dynamic range.
Definition: JDAQClock.hh:173
DEBUG
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
Definition: JCalibratePMT.sh:21
JTRIGGER::JSuperFrame2D
2-dimensional frame with time calibrated data from one optical module.
Definition: JSuperFrame2D.hh:43
JTRIGGER::JTriggerInput
Data structure for input to trigger algorithm.
Definition: JTriggerInput.hh:32
JSUPPORT::JSingleFileScanner
Object reading from a list of files.
Definition: JSingleFileScanner.hh:75
JSUPPORT::JLimit::getLimit
const JLimit & getLimit() const
Get limit.
Definition: JLimit.hh:73
detector
do set_variable DETECTOR_TXT $WORKDIR detector
Definition: detector-XY:fit1d.sh:38
JACOUSTICS::JEventOverlap
Match of two events considering overlap in time.
Definition: JAcoustics/JEventOverlap.hh:20
KM3NETDAQ::JDAQTriggerCounter::getCounter
JTriggerCounter_t getCounter() const
Get trigger counter.
Definition: JDAQTriggerCounter.hh:99
JSUPPORT::getTriggerParameters
JTriggerParameters getTriggerParameters(const JMultipleFileScanner_t &file_list)
Get trigger parameters.
Definition: JTriggerParametersSupportkit.hh:148
KM3NETDAQ::JDAQTriggerCounter::next
JTriggerCounter_t next()
Increment trigger counter.
Definition: JDAQTriggerCounter.hh:121
JTRIGGER::JTimeslice
Time slice with calibrated data.
Definition: JTimeslice.hh:26
Generated by   doxygen 1.8.5