Jpp
JSimbad.cc
Go to the documentation of this file.
1 
2 #include <string>
3 #include <iostream>
4 #include <iomanip>
5 
7 
8 #include "JDetector/JDetector.hh"
12 #include "JDetector/JK40Rates.hh"
19 
22 #include "JSupport/JSupport.hh"
23 #include "JSupport/JMeta.hh"
25 
26 #include "JTools/JRange.hh"
27 
28 #include "Jeep/JPrint.hh"
29 #include "Jeep/JParser.hh"
30 #include "Jeep/JMessage.hh"
31 
32 
33 /**
34  * \file
35  * Example program to simulate summary data.
36  *
37  * \author mdejong
38  */
39 int main(int argc, char **argv)
40 {
41  using namespace std;
42  using namespace JPP;
43  using namespace KM3NETDAQ;
44 
45  JSingleFileScanner<JDAQSummaryslice> inputFile;
46  JFileRecorder<JTYPELIST<JDAQSummaryslice, JMeta>::typelist> outputFile;
47  JLimit_t& numberOfEvents = inputFile.getLimit();
48  string detectorFile;
49  JPMTParametersMap pmtParameters;
50  JPMTRateRange R_Hz;
51  int debug;
52 
53  try {
54 
55  JParser<> zap("Example program to simulate summary data.");
56 
57  zap['f'] = make_field(inputFile, "input DAQ file");
58  zap['o'] = make_field(outputFile, "output file") = "simbad.root";
59  zap['n'] = make_field(numberOfEvents) = JLimit::max();
60  zap['a'] = make_field(detectorFile , "detector file");
61  zap['P'] = make_field(pmtParameters, "PMT simulation data (or corresponding file name)") = JPARSER::initialised();
62  zap['R'] = make_field(R_Hz) = JPMTRateRange();
63  zap['d'] = make_field(debug, "debug") = 0;
64 
65  zap(argc, argv);
66  }
67  catch(const exception &error) {
68  FATAL(error.what() << endl);
69  }
70 
71  cout.tie(&cerr);
72 
73  JDetector detector;
74 
75  try {
76  load(detectorFile, detector);
77  }
78  catch(const JException& error) {
79  FATAL(error);
80  }
81 
82  JPMTParametersMap::Throw(true);
83 
84  if (!pmtParameters.is_valid()) {
85  FATAL("Invalid PMT parameters " << pmtParameters << endl);
86  }
87 
88  const JModuleRouter moduleRouter(detector);
89  JDetectorSimulator simbad (detector);
90  JSummaryRouter summaryRouter;
91 
92  try {
93  simbad.reset(new JK40RunByRunSimulator(summaryRouter, JK40Rates()));
94  simbad.reset(new JPMTRunByRunSimulator(summaryRouter, pmtParameters, detector, R_Hz));
95  simbad.reset(new JCLBDefaultSimulator());
96  }
97  catch(const JException& error) {
98  FATAL(error.what() << endl);
99  }
100 
101  outputFile.open();
102 
103  if (!outputFile.is_open()) {
104  FATAL("Error opening file " << outputFile << endl);
105  }
106 
107  outputFile.put(JMeta(argc, argv));
108 
109  while (inputFile.hasNext()) {
110 
111  STATUS("event: " << setw(10) << inputFile.getCounter() << '\r'); DEBUG(endl);
112 
113  JDAQSummaryslice* summary = inputFile.next();
114 
115  summaryRouter.update(summary);
116 
117  summaryRouter.correct(dynamic_cast<const JPMTDefaultSimulatorInterface&>(simbad.getPMTSimulator()));
118 
119  outputFile.put(JSummaryslice(summary->getDAQChronometer(),simbad));
120  }
121  STATUS(endl);
122 
123  outputFile.close();
124 }
JPMTParametersMap.hh
JMeta.hh
JK40DefaultSimulator.hh
JSingleFileScanner.hh
JFileRecorder.hh
JMessage.hh
JPrint.hh
JK40RunByRunSimulator.hh
JPARSER::initialised
Empty structure for specification of parser element that is initialised (i.e.
Definition: JParser.hh:63
JDETECTOR::load
void load(const JString &file_name, JDetector &detector)
Load detector from input file.
Definition: JDetectorToolkit.hh:456
JSUPPORT::JLimit_t
JLimit JLimit_t
Type definition of limit.
Definition: JLimit.hh:215
JK40Rates.hh
JPARSER::JParser
Utility class to parse command line options.
Definition: JParser.hh:1493
JDAQSummaryslice.hh
JSupport.hh
JPP
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Definition: JAAnetToolkit.hh:37
JPMTRunByRunSimulator.hh
JRange.hh
JDetectorSimulator.hh
JCLBDefaultSimulator.hh
debug
int debug
debug level
Definition: JSirene.cc:59
KM3NETDAQ::JSummaryslice
Auxiliary class to create summary data.
Definition: JSummaryslice.hh:33
JSUPPORT::JLimit::getLimit
const JLimit & getLimit() const
Get limit.
Definition: JLimit.hh:73
KM3NETDAQ::JDAQSummaryslice
Data summary slice.
Definition: JDAQSummaryslice.hh:31
main
int main(int argc, char **argv)
Definition: JSimbad.cc:39
STATUS
#define STATUS(A)
Definition: JMessage.hh:63
JParser.hh
JModuleMapper.hh
JDetectorToolkit.hh
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition: JParser.hh:1954
JSummaryslice.hh
KM3NETDAQ::JDAQChronometer::getDAQChronometer
const JDAQChronometer & getDAQChronometer() const
Get DAQ chronometer.
Definition: JDAQChronometer.hh:88
JSUPPORT::JMeta
Auxiliary class for ROOT I/O of application specific meta data.
Definition: JMeta.hh:71
DEBUG
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62
std
Definition: jaanetDictionary.h:36
JPMTDefaultSimulator.hh
KM3NETDAQ
KM3NeT DAQ data structures and auxiliaries.
Definition: DataQueue.cc:39
JDetector.hh
FATAL
#define FATAL(A)
Definition: JMessage.hh:67
outputFile
string outputFile
Definition: JDAQTimesliceSelector.cc:37