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

Auxiliary program to compose detector from separate calibrations. More...

#include <string>
#include <iostream>
#include <iomanip>
#include <fstream>
#include <vector>
#include <map>
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorCalibration.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JDetector/JDetectorAddressMapToolkit.hh"
#include "JDetector/JModuleRouter.hh"
#include "JDetector/JPMTRouter.hh"
#include "JSupport/JMeta.hh"
#include "JSon/JSon.hh"
#include "JSystem/JDateAndTime.hh"
#include "Jeep/JeepToolkit.hh"
#include "Jeep/JPrint.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 compose detector from separate calibrations.


Note that the validity range as specified in the header of the detector data structure is determined by the maximal date and time value at JSON::ValidFrom_t and the minimal date and time value at JSON::ValidThrough_t, respectively.

Author
mdejong

Definition in file JConstructDetector.cc.

Function Documentation

int main ( int  argc,
char **  argv 
)

Definition at line 61 of file JConstructDetector.cc.

62 {
63  using namespace std;
64  using namespace JPP;
65 
66  string detectorFile;
67  vector<string> inputFile;
68  string outputFile;
69  int debug;
70 
71  try {
72 
73  JParser<> zap("Auxiliary program to compose detector from separate calibrations.");
74 
75  zap['a'] = make_field(detectorFile, "detector file");
76  zap['f'] = make_field(inputFile, "detector calibration files in JSON format "\
77  "(wild card \'" << FILENAME_WILD_CARD << "\' will be replaced by corresponding calibration set)");
78  zap['o'] = make_field(outputFile, "detector file") = "";
79  zap['d'] = make_field(debug) = 1;
80 
81  zap(argc, argv);
82  }
83  catch(const exception &error) {
84  FATAL(error.what() << endl);
85  }
86 
87 
88  JDetector detector;
89 
90  try {
91  load(detectorFile, detector);
92  }
93  catch(const JException& error) {
94  FATAL(error);
95  }
96 
97  if (!hasDetectorAddressMap(detector.getID())) {
98  FATAL("No detector address map for detector identier " << detector.getID() << endl);
99  }
100 
101  const JDetectorAddressMap& demo = getDetectorAddressMap(detector.getID());
102 
103  const JModuleRouter moduleRouter(detector);
104  const JPMTRouter pmtRouter (detector);
105 
106  detector.setUTCTimeRange(JUTCTimeRange());
107 
108 
109  enum {
110  UNDEFINED = 0,
111  DEFINED = 1
112  };
113 
114  struct status_type {
115 
116  status_type() :
117  value(UNDEFINED)
118  {}
119 
120  operator int() const { return value; }
121 
122  status_type& operator=(const int value)
123  {
124  this->value = value;
125 
126  return *this;
127  }
128 
129  int value;
130  };
131 
132  typedef map<int, status_type> map_type;
133 
134  map_type tcal;
135  map_type pcal;
136  map_type rcal;
137  map_type acal;
138  map_type ccal;
139  map_type scal[2];
140 
141 
142  for (const auto& file_name : inputFile) {
143 
144  vector<string> buffer;
145 
146  if (hasWildCard(file_name))
147  buffer = { setWildCard(file_name, TCAL),
148  setWildCard(file_name, PCAL),
149  setWildCard(file_name, RCAL),
150  setWildCard(file_name, ACAL),
151  setWildCard(file_name, CCAL),
152  setWildCard(file_name, SCAL) };
153  else
154  buffer = { file_name };
155 
156  for (const auto& file_name : buffer) {
157 
158  const JSon js(file_name);
159 
160  json::const_iterator data;
161 
162  if (is_valid(js) && (data = js.find(Data_t)) != js.end()) {
163 
164  detector.comment.add(MAKE_STRING("calibration=" << (js.contains(Comment_t) ? js[Comment_t] : "?")));
165 
166  for (size_t i = 0; i != data->size(); ++i) {
167 
168  from_json(data->at(i), detector.getUTCTimeRange());
169 
170  if (data->at(i).contains(DetID_t)) {
171  if (data->at(i)[DetID_t].get<int>() != detector.getID()) {
172  FATAL("Detector identifier mismatch " << data->at(i)[DetID_t].get<int>() << " != " << detector.getID() << endl);
173  }
174  }
175 
176  if (data->at(i).contains(PMTT0s_t)) {
177 
178  for (const auto& element : data->at(i)[PMTT0s_t].get<JPMTCalibration>()) {
179 
180  if (pmtRouter.hasPMT(element.getID())) {
181 
182  JPMT& pmt = detector.getPMT(pmtRouter.getAddress(element.getID()));
183 
184  pmt.setCalibration(element.getCalibration());
185 
186  tcal[pmt.getID()] = DEFINED;
187  }
188  }
189  }
190 
191  if (data->at(i).contains(DOMPositions_t)) {
192 
193  for (const auto& element : data->at(i)[DOMPositions_t].get<JModulePosition>()) {
194 
195  if (moduleRouter.hasModule(element.getID())) {
196 
197  JModule& module = detector.getModule(moduleRouter.getAddress(element.getID()));
198 
199  module.set(element.getPosition());
200 
201  pcal[module.getID()] = DEFINED;
202  }
203  }
204  }
205 
206  if (data->at(i).contains(BasePositions_t)) {
207 
208  for (const auto& element : data->at(i)[BasePositions_t].get<JModulePosition>()) {
209 
210  if (moduleRouter.hasModule(element.getID())) {
211 
212  JModule& module = detector.getModule(moduleRouter.getAddress(element.getID()));
213 
214  module.set(element.getPosition());
215 
216  pcal[module.getID()] = DEFINED;
217  }
218  }
219  }
220 
221  if (data->at(i).contains(DOMRotations_t)) {
222 
223  for (const auto& element : data->at(i)[DOMRotations_t].get<JModuleRotation>()) {
224 
225  if (moduleRouter.hasModule(element.getID())) {
226 
227  JModule& module = detector.getModule(moduleRouter.getAddress(element.getID()));
228 
229  const JModule buffer = getModule(demo.get(module.getID()), module.getID());
230 
231  if (module.size() != buffer.size()) {
232  FATAL("Module size " << module.size() << " != " << buffer.size() << endl);
233  }
234 
235  const JPosition3D center = module.getCenter();
236 
237  for (size_t i = 0; i != module.size(); ++i) {
238  module.getPMT(i).setAxis(buffer.getPMT(i).getAxis());
239  }
240 
241  module.rotate(element.getQuaternion());
242 
243  module.setPosition(module.getCenter());
244  module.set(center);
245 
246  rcal[module.getID()] = DEFINED;
247  }
248  }
249  }
250 
251  if (data->at(i).contains(DOMAcousticT0_t)) {
252 
253  for (const auto& element : data->at(i)[DOMAcousticT0_t].get<JModuleCalibration>()) {
254 
255  if (moduleRouter.hasModule(element.getID())) {
256 
257  JModule& module = detector.getModule(moduleRouter.getAddress(element.getID()));
258 
259  module.setCalibration(element.getCalibration());
260 
261  acal[module.getID()] = DEFINED;
262  }
263  }
264  }
265 
266  if (data->at(i).contains(BaseAcousticT0_t)) {
267 
268  for (const auto& element : data->at(i)[BaseAcousticT0_t].get<JModuleCalibration>()) {
269 
270  if (moduleRouter.hasModule(element.getID())) {
271 
272  JModule& module = detector.getModule(moduleRouter.getAddress(element.getID()));
273 
274  module.setCalibration(element.getCalibration());
275 
276  acal[module.getID()] = DEFINED;
277  }
278  }
279  }
280 
281  if (data->at(i).contains(DOMCompassRotations_t)) {
282 
283  for (const auto& element : data->at(i)[DOMCompassRotations_t].get<JCompassRotation>()) {
284 
285  if (moduleRouter.hasModule(element.getID())) {
286 
287  JModule& module = detector.getModule(moduleRouter.getAddress(element.getID()));
288 
289  module.setQuaternion(element.getQuaternion());
290 
291  ccal[module.getID()] = DEFINED;
292  }
293  }
294  }
295 
296  if (data->at(i).contains(BaseCompassRotations_t)) {
297 
298  for (const auto& element : data->at(i)[BaseCompassRotations_t].get<JCompassRotation>()) {
299 
300  if (moduleRouter.hasModule(element.getID())) {
301 
302  JModule& module = detector.getModule(moduleRouter.getAddress(element.getID()));
303 
304  module.setQuaternion(element.getQuaternion());
305 
306  ccal[module.getID()] = DEFINED;
307  }
308  }
309  }
310 
311  if (data->at(i).contains(PMTStatusInfo_t)) {
312 
313  for (const auto& element : data->at(i)[PMTStatusInfo_t].get<JPMTStatus>()) {
314 
315  if (pmtRouter.hasPMT(element.getID())) {
316 
317  JPMT& pmt = detector.getPMT(pmtRouter.getAddress(element.getID()));
318 
319  pmt.setStatus(element.getStatus());
320 
321  scal[0][pmt.getID()] = DEFINED;
322  }
323  }
324  }
325 
326  if (data->at(i).contains(DOMStatusInfo_t)) {
327 
328  for (const auto& element : data->at(i)[DOMStatusInfo_t].get<JModuleStatus>()) {
329 
330  if (moduleRouter.hasModule(element.getID())) {
331 
332  JModule& module = detector.getModule(moduleRouter.getAddress(element.getID()));
333 
334  module.setStatus(element.getStatus());
335 
336  scal[1][module.getID()] = DEFINED;
337  }
338  }
339  }
340 
341  if (data->at(i).contains(BaseStatusInfo_t)) {
342 
343  for (const auto& element : data->at(i)[BaseStatusInfo_t].get<JModuleStatus>()) {
344 
345  if (moduleRouter.hasModule(element.getID())) {
346 
347  JModule& module = detector.getModule(moduleRouter.getAddress(element.getID()));
348 
349  module.setStatus(element.getStatus());
350 
351  scal[1][module.getID()] = DEFINED;
352  }
353  }
354  }
355  }
356  }
357  }
358  }
359 
360  for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {
361 
362  if (pcal [module->getID()] != DEFINED) { ERROR("Module " << setw(10) << module->getID() << ' ' << getLabel(module->getLocation()) << " no position calibration." << endl); }
363  if (rcal [module->getID()] != DEFINED &&
364  module->getFloor() != 0) { ERROR("Module " << setw(10) << module->getID() << ' ' << getLabel(module->getLocation()) << " no rotation calibration." << endl); }
365  if (acal [module->getID()] != DEFINED) { ERROR("Module " << setw(10) << module->getID() << ' ' << getLabel(module->getLocation()) << " no acoustics calibration." << endl); }
366  if (ccal [module->getID()] != DEFINED) { ERROR("Module " << setw(10) << module->getID() << ' ' << getLabel(module->getLocation()) << " no compass calibration." << endl); }
367  if (scal[1][module->getID()] != DEFINED) { ERROR("Module " << setw(10) << module->getID() << ' ' << getLabel(module->getLocation()) << " no status calibration." << endl); }
368 
369  for (JModule::const_iterator pmt = module->begin(); pmt != module->end(); ++pmt) {
370  if (tcal [pmt->getID()] != DEFINED) { ERROR("PMT " << setw(8) << pmt->getID() << " no time calibration." << endl); }
371  if (scal[0][pmt->getID()] != DEFINED) { ERROR("PMT " << setw(8) << pmt->getID() << " no status calibration." << endl); }
372  }
373  }
374 
375 
376  detector.comment.add(JMeta(argc,argv));
377 
378  try {
379  if (outputFile != "")
380  store(outputFile, detector);
381  else
382  cout << detector << endl;
383  }
384  catch(const JException& error) {
385  FATAL(error);
386  }
387 
388  return 0;
389 }
JDETECTOR::JPMTRouter
Router for direct addressing of PMT data in detector data structure.
Definition: JPMTRouter.hh:35
JSUPPORT::JMeta
Auxiliary class for ROOT I/O of application specific meta data.
Definition: JMeta.hh:70
JPARSER::JParser
Utility class to parse command line options.
Definition: JParser.hh:1514
JLANG::JException
General exception.
Definition: JException.hh:23
JSON::DOMStatusInfo_t
static const std::string DOMStatusInfo_t
Definition: JSon/JSupport.hh:199
JDETECTOR::JLocation::getFloor
int getFloor() const
Get floor number.
Definition: JLocation.hh:145
JDETECTOR::JModule
Data structure for a composite optical module.
Definition: JModule.hh:68
JSON::ACAL
static const std::string ACAL
acoustic time offsets (piezo sensor or hydrophone)
Definition: JSon/JSupport.hh:35
JDETECTOR::getLabel
std::string getLabel(const JLocation &location)
Get module label for monitoring and other applications.
Definition: JLocation.hh:246
JSON::CCAL
static const std::string CCAL
compass alignment (a.k.a. quaternion calibration)
Definition: JSon/JSupport.hh:36
JSON::is_valid
bool is_valid(const json &js)
Check validity of JSon data.
Definition: JSon/JSupport.hh:246
JDETECTOR::JDetector
Detector data structure.
Definition: JDetector.hh:89
JSON::PMTT0s_t
static const std::string PMTT0s_t
Definition: JSon/JSupport.hh:194
JDETECTOR::JModuleRouter
Router for direct addressing of module data in detector data structure.
Definition: JModuleRouter.hh:36
JDETECTOR::hasDetectorAddressMap
bool hasDetectorAddressMap(const int id)
Check if detector address map is available.
Definition: JDetectorAddressMapToolkit.hh:524
JDATABASE::from_json
void from_json(const json &js, JDBString &object)
Convert JSon to database string.
Definition: JSonSupportkit.hh:58
JGEOMETRY3D::JAxis3D::setAxis
void setAxis(const JAxis3D &axis)
Set axis.
Definition: JAxis3D.hh:109
JDETECTOR::JDetectorAddressMap
Lookup table for PMT addresses in detector.
Definition: JDetectorAddressMap.hh:25
outputFile
string outputFile
Definition: JDAQTimesliceSelector.cc:37
JSON::BaseStatusInfo_t
static const std::string BaseStatusInfo_t
Definition: JSon/JSupport.hh:203
MAKE_STRING
#define MAKE_STRING(A)
Make string.
Definition: JPrint.hh:127
JSON::TCAL
static const std::string TCAL
PMT time offsets.
Definition: JSon/JSupport.hh:32
JSON::Comment_t
static const std::string Comment_t
Definition: JSon/JSupport.hh:157
JAANET::detector
Detector file.
Definition: JHead.hh:226
JGEOMETRY3D::JQuaternion3D::setQuaternion
void setQuaternion(const JQuaternion3D &quaternion)
Set quaternion.
Definition: JQuaternion3D.hh:498
JDETECTOR::JLocation::getLocation
const JLocation & getLocation() const
Get location.
Definition: JLocation.hh:69
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition: JParser.hh:1989
JLANG::JObjectID::getID
int getID() const
Get identifier.
Definition: JObjectID.hh:50
JDETECTOR::JPMT
Data structure for PMT geometry, calibration and status.
Definition: JPMT.hh:43
JDETECTOR::store
void store(const std::string &file_name, const JDetector &detector)
Store detector to output file.
Definition: JDetectorToolkit.hh:578
ERROR
#define ERROR(A)
Definition: JMessage.hh:66
JSON::DOMPositions_t
static const std::string DOMPositions_t
Definition: JSon/JSupport.hh:195
JDETECTOR::JModule::getCenter
JVector3D getCenter() const
Get center of module based on crossing point of PMT axes.
Definition: JModule.hh:215
JDETECTOR::getDetectorAddressMap
JDetectorAddressMap & getDetectorAddressMap()
Get detector address map.
Definition: JDetectorAddressMapToolkit.hh:388
JSON::RCAL
static const std::string RCAL
optical module orientations
Definition: JSon/JSupport.hh:34
JDETECTOR::JModule::getPMT
const JPMT & getPMT(const int index) const
Get PMT.
Definition: JModule.hh:173
JDETECTOR::JModule::rotate
void rotate(const JRotation3D &R)
Rotate module.
Definition: JModule.hh:315
JSON::DOMCompassRotations_t
static const std::string DOMCompassRotations_t
Definition: JSon/JSupport.hh:198
JEEP::JStatus::setStatus
void setStatus(const JStatus &status)
Set status.
Definition: JStatus.hh:97
FATAL
#define FATAL(A)
Definition: JMessage.hh:67
JEEP::setWildCard
std::string setWildCard(const std::string &file_name, const std::string &value)
Get file name by setting wild card to given value.
Definition: JeepToolkit.hh:66
JUTC::JUTCTimeRange
UTC time range [s].
Definition: JUTCTimeRange.hh:25
JSON::Data_t
static const std::string Data_t
Definition: JSon/JSupport.hh:185
JDETECTOR::load
void load(const std::string &file_name, JDetector &detector)
Load detector from input file.
Definition: JDetectorToolkit.hh:485
JSON::SCAL
static const std::string SCAL
(module|PMT) status
Definition: JSon/JSupport.hh:37
JDETECTOR::JCalibration::setCalibration
void setCalibration(const JCalibration &cal)
Set calibration.
Definition: JDetector/JCalibration.hh:89
JSON::BasePositions_t
static const std::string BasePositions_t
Definition: JSon/JSupport.hh:200
JSON::DOMRotations_t
static const std::string DOMRotations_t
Definition: JSon/JSupport.hh:196
JGEOMETRY3D::JAxis3D::getAxis
const JAxis3D & getAxis() const
Get axis.
Definition: JAxis3D.hh:98
JEEP::FILENAME_WILD_CARD
static const char FILENAME_WILD_CARD
wild card character for file name substitution
Definition: JeepToolkit.hh:44
JGEOMETRY3D::JPosition3D
Data structure for position in three dimensions.
Definition: JPosition3D.hh:36
JSON::BaseCompassRotations_t
static const std::string BaseCompassRotations_t
Definition: JSon/JSupport.hh:202
JSON::JSon
Auxiliary class to load json data from file.
Definition: JSon.hh:61
detector
do set_variable DETECTOR_TXT $WORKDIR detector
Definition: detector-XY:fit1d.sh:38
JSON::PMTStatusInfo_t
static const std::string PMTStatusInfo_t
Definition: JSon/JSupport.hh:210
JSON::PCAL
static const std::string PCAL
(optical|base) module positions
Definition: JSon/JSupport.hh:33
JSON::DetID_t
static const std::string DetID_t
Definition: JSon/JSupport.hh:219
JDETECTOR::JModule::set
JModule & set(const JVector3D &pos)
Set position.
Definition: JModule.hh:408
JSON::DOMAcousticT0_t
static const std::string DOMAcousticT0_t
Definition: JSon/JSupport.hh:197
getModule
const JModule & getModule(const JDetector &detector, const JModuleLocation &location)
find module with a given string and floor number
Definition: JBeaconSimulator.cc:928
debug
int debug
debug level
Definition: archive-put-wiki-detectors.sh:92
JSON::BaseAcousticT0_t
static const std::string BaseAcousticT0_t
Definition: JSon/JSupport.hh:201
JGEOMETRY3D::JPosition3D::setPosition
void setPosition(const JVector3D &pos)
Set position.
Definition: JPosition3D.hh:152
JEEP::hasWildCard
bool hasWildCard(const std::string &file_name)
Check presence of wild card.
Definition: JeepToolkit.hh:53
Generated by   doxygen 1.8.5