Jpp  15.0.5
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JShowerDirectionPrefit.cc
Go to the documentation of this file.
1 #include <string>
2 #include <limits>
3 
6 
7 #include "JDAQ/JDAQEventIO.hh"
10 
13 
14 #include "JDynamics/JDynamics.hh"
15 
17 #include "JSupport/JSupport.hh"
19 #include "JSupport/JMeta.hh"
20 #include "JSupport/JTreeScanner.hh"
22 
23 #include "Jeep/JParser.hh"
24 #include "Jeep/JMessage.hh"
25 
26 #include "JReconstruction/JEvt.hh"
30 
31 #include "JLang/JSharedPointer.hh"
32 
33 
34 /**
35  * \file
36  *
37  * Program to perform EM Shower Fit for ORCA with I/O of JFIT::JEvt data.
38  * The reconstruction is made at the PMT level.
39  * \author adomi
40  */
41 
42 int main(int argc, char** argv){
43 
44  using namespace std;
45  using namespace JPP;
46  using namespace KM3NETDAQ;
47 
49  typedef JParallelFileScanner< JTypeList<JDAQEvent, JFIT::JEvt> > JParallelFileScanner_t;
50  typedef JParallelFileScanner_t::multi_pointer_type multi_pointer_type;
52  JACOUSTICS::JEvt>::typelist calibration_types;
53  typedef JMultipleFileScanner<calibration_types> JCalibration_t;
54 
55  JParallelFileScanner_t inputFile;
58  JLimit_t& numberOfEvents = inputFile.getLimit();
59  string detectorFile;
60  JCalibration_t calibrationFile;
61  double Tmax_s;
62  string pdfFile;
63  int debug;
64 
65  try {
66 
67  JParser<> zap;
68 
69  zap['f'] = make_field(inputFile) ;
70  zap['o'] = make_field(outputFile) = "JORCAShowerDirectionPrefit_Output.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  getMechanics.load(detector.getID());
95 
96  unique_ptr<JDynamics> dynamics;
97 
98  try {
99 
100  dynamics.reset(new JDynamics(detector, Tmax_s));
101 
102  dynamics->load(calibrationFile);
103  }
104  catch(const exception& error) {
105  if (!calibrationFile.empty()) {
106  FATAL(error.what());
107  }
108  }
109 
110  const JModuleRouter router(dynamics ? dynamics->getDetector() : detector);
111 
112  JSummaryFileRouter summary(inputFile, parameters.R_Hz);
113 
114  JShowerDirectionPrefit fit(parameters, router, summary, pdfFile, debug);
115 
116  outputFile.open();
117  outputFile.put(JMeta(argc, argv));
118 
119  while (inputFile.hasNext()) {
120 
121  STATUS("event: " << setw(10) << inputFile.getCounter() << '\r'); DEBUG(endl);
122 
123  multi_pointer_type ps = inputFile.next();
124 
125  JDAQEvent* tev = ps;
126  JFIT::JEvt* in = ps;
127 
128  summary.update(*tev);
129 
130  if (dynamics) {
131  dynamics->update(*tev);
132  }
133 
134  // select start values
135  JFIT::JEvt cp = *in;
136 
137  cp.select(parameters.numberOfPrefits, qualitySorter);
138 
139  if (!cp.empty()) {
140  cp.select(JHistory::is_event(cp.begin()->getHistory()));
141  }
142 
143  // fit
144 
145  JFIT::JEvt out = fit(*tev, cp);
146 
147  if (dynamics) {
148 
149  const JDynamics::coverage_type coverage = dynamics->getCoverage();
150 
151  for (JFIT::JEvt::iterator i = out.begin(); i != out.end(); ++i) {
152  i->setW(JPP_COVERAGE_ORIENTATION, coverage.orientation);
153  i->setW(JPP_COVERAGE_POSITION, coverage.position);
154  }
155  }
156 
157  // apply default sorter
158 
159  sort(out.begin(), out.end(), qualitySorter);
160 
161  // reduce output fits
162 
163  JFIT::JEvt::iterator __end = out.end();
164 
165  advance(__end = out.begin(), min(parameters.numberOfOutfits, (size_t) distance(out.begin(), out.end())));
166 
167  JFIT::JEvt reduced_out;
168 
169  copy(out.begin(), __end, back_inserter(reduced_out));
170 
171  copy(in->begin(), in->end(), back_inserter(reduced_out));
172 
173  outputFile.put(reduced_out);
174 
175  }
176  STATUS(endl);
177 
179 
180  io >> outputFile;
181 
182  outputFile.close();
183 }
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:1500
General exception.
Definition: JException.hh:23
int main(int argc, char *argv[])
Definition: Main.cc:15
ROOT TTree parameter settings of various packages.
class to handle the direction fit of the shower reconstruction, mainly dedicated for ORCA ...
double position
coverage of detector by available position calibration [0,1]
Definition: JDynamics.hh:536
static JDetectorMechanics getMechanics
Function object to get string mechanics.
Definition: JMechanics.hh:243
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
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:66
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:535
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:196
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:1961
Auxiliary class to test history.
Definition: JHistory.hh:104
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.
void load(const std::string &file_name)
Load mechanical model parameters from file.
Definition: JMechanics.hh:142
int debug
debug level
Definition: JSirene.cc:63
General purpose messaging.
Data structure for coverage of dynamic calibrations.
Definition: JDynamics.hh:534
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:61
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.
Utility class to parse command line options.
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
void copy(const Head &from, JHead &to)
Copy header from from to to.
Definition: JHead.cc:139
then cp
Orientation of module.
Object reading from a list of files.
const JLimit & getLimit() const
Get limit.
Definition: JLimit.hh:73
do set_variable DETECTOR_TXT $WORKDIR detector
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 source JAcoustics sh $DETECTOR_ID CHECK_EXIT_CODE typeset A TRIPODS get_tripods $WORKDIR tripod txt TRIPODS for EMITTER in
Definition: JCanberra.sh:41
void select(const JSelector_t &selector)
Select fits.
then usage $script< detector file >< tripodfile >< stage > input file nInput files correspond to the output of JAcousticsEventBuilder[.sh] nFirst stage eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY eval JPrintDetector a $DETECTOR O CAN source JAcoustics sh $DETECTOR_ID typeset A CONFIGURATION for key in Tmax_s
static const int JPP_COVERAGE_POSITION
coverage of dynamic position calibration from any Jpp application