Jpp
JORCAShowerPositionFit.cc
Go to the documentation of this file.
1 #include <string>
2 
3 #include "evt/Head.hh"
4 #include "evt/Evt.hh"
5 
6 #include "JDAQ/JDAQEvent.hh"
7 #include "JDAQ/JDAQTimeslice.hh"
8 
11 
13 #include "JSupport/JSupport.hh"
15 #include "JSupport/JMeta.hh"
16 #include "JSupport/JTreeScanner.hh"
17 
18 #include "Jeep/JParser.hh"
19 #include "Jeep/JMessage.hh"
20 
21 #include "JFit/JEvt.hh"
25 
26 #include "JLang/JSharedPointer.hh"
27 
28 int main(int argc, char **argv){
29 
30  using namespace std;
31  using namespace JPP;
32  using namespace KM3NETDAQ;
33 
34  typedef JParallelFileScanner< JTypeList<JDAQEvent, JEvt> > JParallelFileScanner_t;
36 
37  JParallelFileScanner_t inputFile;
38  JFileRecorder<typelist> outputFile;
39  JLimit_t& numberOfEvents = inputFile.getLimit();
40  string detectorFile;
41  string pdfFile;
42  JShowerParameters showerParams;
43  int debug;
44 
45  try{
46  JParser<> zap;
47 
48  zap['f'] = make_field(inputFile);
49  zap['o'] = make_field(outputFile) = "PositionFit.root";
50  zap['a'] = make_field(detectorFile);
51  zap['F'] = make_field(pdfFile) = "PDF_ShowerPositionFit.root";
52  zap['n'] = make_field(numberOfEvents) = JLimit::max();
53  zap['@'] = make_field(showerParams) = JPARSER::initialised();
54  zap['d'] = make_field(debug) = 2;
55 
56  if (zap.read(argc, argv) != 0) return 1;
57  }
58  catch(const exception& error){
59  FATAL(error.what() << endl);
60  }
61 
62  using namespace JTRIGGER;
63  using namespace JDETECTOR;
64  using namespace JGEOMETRY3D;
65  using namespace JTOOLS;
66 
67  JDetector detector;
68 
69  try {
70  load(detectorFile, detector);
71  }
72  catch(const JException& error) {
73  FATAL(error);
74  }
75 
76  outputFile.open();
77  outputFile.put(JMeta(argc, argv));
78 
79  const JSharedPointer<const JModuleRouter> moduleRouter(new JModuleRouter(detector));
80 
81  const JORCAShowerPositionFit ShowerPosFit(moduleRouter, showerParams, pdfFile);
82 
83  while (inputFile.hasNext()) {
84 
85  STATUS("event: " << setw(10) << inputFile.getCounter() << '\r'); DEBUG(endl);
86 
87  JParallelFileScanner_t::multi_pointer_type ps = inputFile.next();
88 
89  JDAQEvent* tev = ps;
90  JEvt* in = ps;
91 
92  JEvt out;
93 
94  const JDAQTimeslice timeSliceBuildL0(*tev, true);
95  const JDAQTimeslice timeSliceBuildL2(*tev, false);
96 
97  ShowerPosFit.getJEvt(timeSliceBuildL0, timeSliceBuildL2, *in, out);
98 
99  outputFile.put(out);
100  }
101  STATUS(endl);
102 
103  JMultipleFileScanner<JRemove<typelist, JEvt>::typelist> io(inputFile);
104 
105  io >> outputFile;
106 
107  outputFile.close();
108 }
JMeta.hh
KM3NETDAQ::JDAQEvent
DAQ Event.
Definition: JDAQEvent.hh:34
JFileRecorder.hh
JMessage.hh
JPARSER::initialised
Empty structure for specification of parser element that is initialised (i.e.
Definition: JParser.hh:63
JShowerPositionFitParameters_t.hh
JDETECTOR::load
void load(const JString &file_name, JDetector &detector)
Load detector from input file.
Definition: JDetectorToolkit.hh:456
JSharedPointer.hh
JSUPPORT::JLimit_t
JLimit JLimit_t
Type definition of limit.
Definition: JLimit.hh:215
KM3NETDAQ::JDAQTimeslice
Data time slice.
Definition: JDAQTimeslice.hh:36
JPARSER::JParser
Utility class to parse command line options.
Definition: JParser.hh:1493
JDAQTimeslice.hh
JTreeScanner.hh
JSupport.hh
JShowerParameters.hh
JPP
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Definition: JAAnetToolkit.hh:37
JEvt.hh
debug
int debug
debug level
Definition: JSirene.cc:59
JSUPPORT::JLimit::getLimit
const JLimit & getLimit() const
Get limit.
Definition: JLimit.hh:73
JModuleRouter.hh
JParallelFileScanner.hh
JLANG::JTypeList
Type list.
Definition: JTypeList.hh:22
STATUS
#define STATUS(A)
Definition: JMessage.hh:63
JParser.hh
JDetectorToolkit.hh
JDETECTOR::JModuleRouter
Router for direct addressing of module data in detector data structure.
Definition: JModuleRouter.hh:34
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition: JParser.hh:1954
JDETECTOR::JDetector
Detector data structure.
Definition: JDetector.hh:80
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
KM3NETDAQ
KM3NeT DAQ data structures and auxiliaries.
Definition: DataQueue.cc:39
JORCAShowerPositionFit.hh
JDAQEvent.hh
JGEOMETRY3D
Auxiliary classes and methods for 3D geometrical objects and operations.
Definition: JAngle3D.hh:18
JTRIGGER
Checksum.
Definition: JSupport/JSupport.hh:35
JTOOLS
Auxiliary classes and methods for multi-dimensional interpolations and histograms.
Definition: JAbstractCollection.hh:9
FATAL
#define FATAL(A)
Definition: JMessage.hh:67
outputFile
string outputFile
Definition: JDAQTimesliceSelector.cc:37
JDETECTOR
Auxiliary classes and methods for detector calibration.
Definition: JAnchor.hh:12
main
int main(int argc, char **argv)
Definition: JORCAShowerPositionFit.cc:28
JPARSER::JParser::read
int read(const int argc, const char *const argv[])
Parse the program's command line options.
Definition: JParser.hh:1791