Jpp  18.3.0
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JShowerPositionFit.cc
Go to the documentation of this file.
1 #include <string>
2 #include <limits>
3 
7 
8 #include "JDAQ/JDAQEventIO.hh"
11 
14 
15 #include "JDynamics/JDynamics.hh"
16 
18 #include "JSupport/JSupport.hh"
20 #include "JSupport/JMeta.hh"
21 #include "JSupport/JTreeScanner.hh"
23 
24 #include "Jeep/JParser.hh"
25 #include "Jeep/JMessage.hh"
26 
27 #include "JReconstruction/JEvt.hh"
30 
31 #include "JLang/JSharedPointer.hh"
32 
33 
34 /**
35  * \file
36  *
37  * Program to perform EM Shower Second Position Fit for ORCA with I/O of JFIT::JEvt data.
38  * The reconstruction is made at the HIT level.
39  * \author adomi
40  */
41 
42 int main(int argc, char** argv){
43 
44  using namespace std;
45  using namespace JPP;
46 
48  typedef JParallelFileScanner< JTypeList<JDAQEvent, JFIT::JEvt> > JParallelFileScanner_t;
49  typedef JParallelFileScanner_t::multi_pointer_type multi_pointer_type;
51  JACOUSTICS::JEvt>::typelist calibration_types;
52  typedef JMultipleFileScanner<calibration_types> JCalibration_t;
53 
54 
55  JParallelFileScanner_t inputFile;
58  JLimit_t& numberOfEvents = inputFile.getLimit();
59  string detectorFile;
60  string pdfFile;
61  JCalibration_t calibrationFile;
62  double Tmax_s;
63  int debug;
64 
65  try {
66 
67  JParser<> zap;
68 
69  zap['f'] = make_field(inputFile) ;
70  zap['o'] = make_field(outputFile) = "JShowerPositionFit.root";
71  zap['a'] = make_field(detectorFile) ;
72  zap['+'] = make_field(calibrationFile) = JPARSER::initialised();
73  zap['T'] = make_field(Tmax_s) = 100.0;
74  zap['F'] = make_field(pdfFile) ;
76  zap['n'] = make_field(numberOfEvents) = JLimit::max();
77  zap['d'] = make_field(debug) = 2;
78 
79  zap(argc, argv);
80  }
81  catch(const exception& error) {
82  FATAL(error.what() << endl);
83  }
84 
86 
87  try {
88  load(detectorFile, detector);
89  }
90  catch(const JException& error) {
91  FATAL(error);
92  }
93 
94  unique_ptr<JDynamics> dynamics;
95 
96  try {
97 
98  dynamics.reset(new JDynamics(detector, Tmax_s));
99 
100  dynamics->load(calibrationFile);
101  }
102  catch(const exception& error) {
103  if (!calibrationFile.empty()) {
104  FATAL(error.what());
105  }
106  }
107 
108  const JModuleRouter router(dynamics ? dynamics->getDetector() : detector);
109 
110  outputFile.open();
111  outputFile.put(JMeta(argc, argv));
112 
113  JSummaryFileRouter summary(inputFile, parameters.R_Hz);
114 
115  JShowerPositionFit fit(parameters, router, summary, pdfFile, debug);
116 
117  while (inputFile.hasNext()) {
118 
119  STATUS("event: " << setw(10) << inputFile.getCounter() << '\r'); DEBUG(endl);
120 
121  multi_pointer_type ps = inputFile.next();
122 
123  JDAQEvent* tev = ps;
124  JFIT::JEvt* in = ps;
125 
126  summary.update(*tev);
127 
128  if (dynamics) {
129  dynamics->update(*tev);
130  }
131 
132  // select start values
133  JFIT::JEvt cp = *in;
134 
135  cp.select(parameters.numberOfPrefits, qualitySorter);
136 
137  if (!cp.empty()) {
138  cp.select(JHistory::is_event(cp.begin()->getHistory()));
139  }
140 
141  // fit
142  JFIT::JEvt out = fit(*tev, cp);
143 
144  if (dynamics) {
145 
146  const JDynamics::coverage_type coverage = dynamics->getCoverage();
147 
148  for (JFIT::JEvt::iterator i = out.begin(); i != out.end(); ++i) {
149  i->setW(JPP_COVERAGE_ORIENTATION, coverage.orientation);
150  i->setW(JPP_COVERAGE_POSITION, coverage.position);
151  }
152  }
153 
154  // apply default sorter
155 
156  sort(out.begin(), out.end(), qualitySorter);
157 
158  // reduce output fits
159 
160  JFIT::JEvt::iterator __end = out.end();
161 
162  advance(__end = out.begin(), min(parameters.numberOfOutfits, (size_t) distance(out.begin(), out.end())));
163 
164  JFIT::JEvt reduced_out;
165 
166  copy(out.begin(), __end, back_inserter(reduced_out));
167 
168  copy(in->begin(), in->end(), back_inserter(reduced_out));
169 
170  outputFile.put(reduced_out);
171 
172  }
173  STATUS(endl);
174 
176 
177  io >> outputFile;
178 
179  outputFile.close();
180 }
Auxiliary class for ROOT I/O of application specific meta data.
Definition: JMeta.hh:70
Object writing to file.
Utility class to parse command line options.
Definition: JParser.hh:1514
General exception.
Definition: JException.hh:24
int main(int argc, char *argv[])
Definition: Main.cc:15
ROOT TTree parameter settings of various packages.
double position
coverage of detector by available position calibration [0,1]
Definition: JDynamics.hh:575
std::vector< T >::difference_type distance(typename std::vector< T >::const_iterator first, typename PhysicsEvent::const_iterator< T > second)
Specialisation of STL distance.
#define STATUS(A)
Definition: JMessage.hh:63
void update(const JDAQHeader &header)
Update router.
Detector data structure.
Definition: JDetector.hh:89
Recording of objects on file according a format that follows from the file name extension.
Router for direct addressing of module data in detector data structure.
*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
General purpose class for parallel reading of objects from a single file or multiple files...
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Definition: JParser.hh:83
class to handle the second position fit of the shower reconstruction, mainly dedicated for ORCA ...
string outputFile
bool qualitySorter(const JFit &first, const JFit &second)
Comparison of fit results.
double orientation
coverage of detector by available orientation calibration [0,1]
Definition: JDynamics.hh:574
Type list.
Definition: JTypeList.hh:22
Auxiliary class for defining the range of iterations of objects.
Definition: JLimit.hh:41
Detector file.
Definition: JHead.hh:226
Acoustic event fit.
Auxiliary class for recursive type list generation.
Definition: JTypeList.hh:351
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition: JParser.hh:1989
Auxiliary class to test history.
Definition: JHistory.hh:105
Parallel scanning of objects from a single file or multiple files according a format that follows fro...
ROOT I/O of application specific meta data.
Dynamic detector calibration.
File router for fast addressing of summary data.
General purpose messaging.
Data structure for coverage of dynamic calibrations.
Definition: JDynamics.hh:573
counter_type advance(counter_type &counter, const counter_type value, const counter_type limit=std::numeric_limits< counter_type >::max())
Advance counter.
#define FATAL(A)
Definition: JMessage.hh:67
Direct access to module in detector data structure.
static const int JPP_COVERAGE_ORIENTATION
coverage of dynamic orientation calibration from any Jpp application
Dynamic detector calibration.
Definition: JDynamics.hh:81
void load(const std::string &file_name, JDetector &detector)
Load detector from input file.
Data structure for set of track fit results.
General purpose class for object reading from a list of file names.
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2] eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY JAcoustics sh $DETECTOR_ID source JAcousticsToolkit sh CHECK_EXIT_CODE typeset A EMITTERS get_tripods $WORKDIR tripod txt EMITTERS get_transmitters $WORKDIR transmitter txt EMITTERS for EMITTER in
Definition: JCanberra.sh:48
Utility class to parse command line options.
void copy(const Head &from, JHead &to)
Copy header from from to to.
Definition: JHead.cc:162
then cp
Orientation of module.
Object reading from a list of files.
const JLimit & getLimit() const
Get limit.
Definition: JLimit.hh:84
do set_variable DETECTOR_TXT $WORKDIR detector
void select(const JSelector_t &selector)
Select fits.
int debug
debug level
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62
static const int JPP_COVERAGE_POSITION
coverage of dynamic position calibration from any Jpp application