Jpp 20.0.0-27-g39925593c-D
the software that should make you happy
Loading...
Searching...
No Matches
JPrintDetector.cc
Go to the documentation of this file.
1#include <string>
2#include <iostream>
3#include <iomanip>
4#include <set>
5#include <iterator>
6
7#include "JDetector/JDetector.hh"
8#include "JDetector/JDetectorToolkit.hh"
9#include "JDetector/JDetectorAddressMap.hh"
10#include "JDetector/JDetectorSupportkit.hh"
11#include "JDetector/JModule.hh"
12
13#include "JMath/JConstants.hh"
14#include "JTools/JRange.hh"
15#include "JSystem/JDateAndTime.hh"
16#include "JGeometry3D/JGeometry3DToolkit.hh"
17#include "JGeometry3D/JCylinder3D.hh"
18#include "JGeometry3D/JQuaternion3D.hh"
19
20#include "Jeep/JPrint.hh"
21#include "Jeep/JParser.hh"
22#include "Jeep/JMessage.hh"
23
24
25namespace {
26 /**
27 * Print options.
28 */
29 const char* const default_t = "default"; //!< as-is
30
31 const char* const modules_t = "modules"; //!< modules
32 const char* const pmts_t = "pmts"; //!< PMTs
33 const char* const geometry_t = "geometry"; //!< geometry
34
35 const char* const comment_t = "comment"; //!< comment
36 const char* const header_t = "header"; //!< header
37 const char* const identifier_t = "identifier"; //!< identifier
38 const char* const version_t = "version"; //!< version
39 const char* const can_t = "can"; //!< can
40 const char* const center_t = "center"; //!< center
41 const char* const summary_t = "summary"; //!< summary
42
43 const char* const GEOMETRY_t = "GEOMETRY"; //!< geometry in zsh format
44 const char* const HEADER_t = "HEADER"; //!< header in zsh format
45 const char* const IDENTIFIER_t = "IDENTIFIER"; //!< identifier in zsh format
46 const char* const VERSION_t = "VERSION"; //!< version in zsh format
47 const char* const CAN_t = "CAN"; //!< can in zsh format
48 const char* const CENTER_t = "CENTER"; //!< center in zsh format
49 const char* const SUMMARY_t = "SUMMARY"; //!< summary in zsh format
50
51 /**
52 * Get zsh array name for given string.
53 *
54 * \param id string number
55 * \return zsh array name
56 */
57 inline std::string getString(const int id)
58 {
59 return MAKE_STRING("STRING_" << FILL(4,'0') << id);
60 }
61}
62
63
64/**
65 * \file
66 *
67 * Auxiliary program to print detector file in human friendly format.
68 * \author mdejong
69 */
70int main(int argc, char **argv)
71{
72 using namespace std;
73 using namespace JPP;
74
75 string detectorFile;
76 set<string> option;
77 double precision;
78 int debug;
79
80 try {
81
82 JParser<> zap("Auxiliary program to print detector file in human friendly format.");
83
84 zap['a'] = make_field(detectorFile);
85 zap['O'] = make_field(option, "possible options "
86
87 << default_t << ' '
88
89 << pmts_t << ' '
90 << modules_t << ' '
91 << geometry_t << ' '
92
93 << comment_t << ' '
94 << header_t << ' '
95 << version_t << ' '
96 << identifier_t << ' '
97 << can_t << ' '
98 << center_t << ' '
99 << summary_t << ' '
100
101 << GEOMETRY_t << ' '
102 << HEADER_t << ' '
103 << VERSION_t << ' '
104 << IDENTIFIER_t << ' '
105 << CAN_t << ' '
106 << CENTER_t << ' '
107 << SUMMARY_t) = JPARSER::initialised();
108
109 zap['p'] = make_field(precision, "precision for match with reference module") = 1.0e-5;
110 zap['d'] = make_field(debug) = 1;
111
112 zap(argc, argv);
113 }
114 catch(const exception &error) {
115 FATAL(error.what() << endl);
116 }
117
118
119 JDetector detector;
120
121 try {
122 load(detectorFile, detector);
123 }
124 catch(const JException& error) {
125 FATAL(error);
126 }
127
128 if (option.empty() || option.count(default_t)) {
129
130 cout << detector << endl;
131
132 }
133
134 if (option.count(modules_t)) {
135
136 const JDetectorBuilder& demo = getDetectorBuilder(detector.getID());
137
138 for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {
139
140 cout << "Module";
141 cout << ' ' << noshowpos << setw(8) << right << module->getID();
142 cout << ' ' << noshowpos << setw(3) << right << module->getString();
143 cout << ' ' << noshowpos << setw(2) << right << module->getFloor();
144 cout << ' ' << FIXED(7,2) << module->getX();
145 cout << ' ' << FIXED(7,2) << module->getY();
146 cout << ' ' << FIXED(7,2) << module->getZ();
147 cout << ' ' << FIXED(8,2) << module->getT0();
148
149 if (module->getFloor() != 0) {
150 {
151 JModule buffer = demo.getModule(module->getID(), module->getLocation());
152
153 const JQuaternion3D Q = getRotation(buffer, *module);
154 const JQuaternion3D::decomposition q1(Q, JVector3Z_t);
155
156 const double phi = (JVector3Z_t.getDot(q1.twist) >= 0.0 ? +1.0 : -1.0) * q1.twist.getAngle();
157
158 cout << ' ' << FIXED(7,2) << phi;
159 }
160 {
161 const JQuaternion3D Q = module->getQuaternion();
162 const JQuaternion3D::decomposition q1(Q, JVector3Z_t);
163
164 const double phi = (JVector3Z_t.getDot(q1.twist) >= 0.0 ? +1.0 : -1.0) * q1.twist.getAngle();
165
166 cout << ' ' << FIXED(7,2) << phi;
167 }
168 {
169 JModule buffer = demo.getModule(module->getID(), module->getLocation());
170
171 const JRotation3D R = getRotation(buffer, *module);
172
173 buffer.rotate(R);
174
175 cout << ' ' << (JModule::compare(buffer, *module, precision) ? "==" : "!=")
176 << ' ' << (demo.get(module->getID()) == demo.getDefaultModuleAddressMap() ? "default" : "custom");
177 }
178 }
179
180 cout << endl;
181 }
182
183 }
184
185 if (option.count(pmts_t)) {
186
187 const JDetectorAddressMap& demo = getDetectorAddressMap(detector.getID());
188
189 for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {
190
191 const JModuleAddressMap& memo = demo.get(module->getID());
192
193 for (size_t i = 0; i != module->size(); ++i) {
194
195 const JPMT& pmt = module->getPMT(i);
196 const JPMTPhysicalAddress& address = memo.getPMTPhysicalAddress(i);
197
198 cout << "PMT[" << FILL(2,'0') << i << FILL() << "]";
199 cout << ' ' << address;
200 cout << ' ' << noshowpos << setw(8) << right << pmt.getID();
201 cout << ' ' << FIXED(7,2) << pmt.getX();
202 cout << ' ' << FIXED(7,2) << pmt.getY();
203 cout << ' ' << FIXED(7,2) << pmt.getZ();
204 cout << ' ' << FIXED(6,3) << pmt.getDX();
205 cout << ' ' << FIXED(6,3) << pmt.getDY();
206 cout << ' ' << FIXED(6,3) << pmt.getDZ();
207 cout << ' ' << FIXED(8,2) << pmt.getT0();
208 cout << ' ' << noshowpos << setw(8) << right << pmt.getStatus();
209 cout << endl;
210 }
211 }
212
213 }
214
215 if (option.count(geometry_t) ||
216 option.count(GEOMETRY_t)) {
217
218 if (option.count(geometry_t)) {
219
220 cout << "Maximal distance [m] = " << FIXED(6,1) << getMaximalDistance(detector) << endl;
221 cout << "Minimal distance [m] = " << FIXED(6,1) << getMinimalDistance(detector) << endl;
222 cout << "Maximal time [ns] = " << FIXED(6,1) << getMaximalTime (detector) << endl;
223
224 } else {
225
226 cout << "set_variable MAXIMAL_DISTANCE_M " << FIXED(6,1) << getMaximalDistance(detector) << ";" << endl;
227 cout << "set_variable MINIMAL_DISTANCE_M " << FIXED(6,1) << getMinimalDistance(detector) << ";" << endl;
228 }
229
230 }
231
232 if (option.count(comment_t)) {
233
234 cout << detector.comment << endl;
235
236 }
237
238 if (option.count(header_t) ||
239 option.count(HEADER_t)) {
240
241 if (option.count(header_t)) {
242
243 cout << detector.getProperties() << endl;
244 cout << "name = " << (isARCADetector(detector) ? "ARCA" :
245 isORCADetector(detector) ? "ORCA" :
246 "unknown") << endl;
247 cout << "validity = "
248 << JDateAndTime(detector.getLowerLimit(), true).toString() << ' '
249 << JDateAndTime(detector.getUpperLimit(), true).toString() << endl;
250
251 } else {
252
253 cout << "set_variable UTM_EAST " << FIXED(12,2) << detector.getUTMEast() << ";" << endl;
254 cout << "set_variable UTM_NORTH " << FIXED(12,2) << detector.getUTMNorth() << ";" << endl;
255 cout << "set_variable UTM_Z " << FIXED(12,2) << detector.getUTMZ() << ";" << endl;
256 cout << "set_variable UTM_ZONE " << FIXED(12,2) << detector.getUTMZone() << ";" << endl;
257 cout << "set_variable UTM_WGS " << FIXED(12,2) << detector.getWGS() << ";" << endl;
258 }
259
260 }
261
262 if (option.count(version_t) ||
263 option.count(VERSION_t)) {
264
265 if (option.count(version_t)) {
266
267 cout << detector.getVersion() << endl;
268
269 } else {
270
271 cout << "set_variable DETECTOR_VERSION " << detector.getVersion() << ";" << endl;
272 }
273
274 }
275
276 if (option.count(identifier_t) ||
277 option.count(IDENTIFIER_t)) {
278
279 if (option.count(identifier_t)) {
280
281 cout << "Detector " << detector.getID() << endl;
282
283 } else {
284
285 cout << "set_variable DETECTOR_ID " << detector.getID() << ";" << endl;
286 }
287
288 }
289
290 if (option.count(can_t) ||
291 option.count(CAN_t)) {
292
293 const JCylinder3D cylinder(detector.begin(), detector.end());
294
295 const double V = (cylinder.getZmax() - cylinder.getZmin()) * PI * cylinder.getRadius() * cylinder.getRadius();
296 const double D = getMaximalDistance(detector);
297
298 if (option.count(can_t)) {
299
300 cout << "X = " << FIXED(7,1) << cylinder.getX() << endl;
301 cout << "Y = " << FIXED(7,1) << cylinder.getY() << endl;
302 cout << "Zmin = " << FIXED(7,1) << cylinder.getZmin() << endl;
303 cout << "Zmax = " << FIXED(7,1) << cylinder.getZmax() << endl;
304 cout << "Radius = " << FIXED(7,1) << cylinder.getRadius() << endl;
305 cout << "Depth = " << FIXED(7,1) << detector.getUTMZ() << endl;
306 cout << "Volume = " << SCIENTIFIC(12,3) << V << endl;
307 cout << "Distance = " << FIXED(9,3) << D << endl;
308
309 } else {
310
311 cout << "set_variable CAN_X_M " << FIXED(7,1) << cylinder.getX() << ";" << endl;
312 cout << "set_variable CAN_Y_M " << FIXED(7,1) << cylinder.getY() << ";" << endl;
313 cout << "set_variable CAN_ZMIN_M " << FIXED(7,1) << cylinder.getZmin() << ";" << endl;
314 cout << "set_variable CAN_ZMAX_M " << FIXED(7,1) << cylinder.getZmax() << ";" << endl;
315 cout << "set_variable CAN_RADIUS_M " << FIXED(7,1) << cylinder.getRadius() << ";" << endl;
316 cout << "set_variable CAN_DEPTH_M " << FIXED(7,1) << detector.getUTMZ() << ";" << endl;
317 cout << "set_variable CAN_VOLUME_M3 " << SCIENTIFIC(12,3) << V << ";" << endl;
318 cout << "set_variable CAN_DISTANCE_M " << FIXED(9,3) << D << ";" << endl;
319 }
320
321 }
322
323 if (option.count(center_t) ||
324 option.count(CENTER_t)) {
325
326 const JCenter3D center(detector.begin(), detector.end());
327
328 if (option.count(center_t)) {
329
330 cout << "center = ";
331 cout << showpos << FIXED(8,3) << center.getX() << ' ';
332 cout << showpos << FIXED(8,3) << center.getY() << ' ';
333 cout << showpos << FIXED(8,3) << center.getZ() << endl;
334
335 } else {
336
337 cout << "set_variable CENTER_X_M " << FIXED(7,1) << center.getX() << ";" << endl;
338 cout << "set_variable CENTER_Y_M " << FIXED(7,1) << center.getY() << ";" << endl;
339 cout << "set_variable CENTER_Z_M " << FIXED(7,1) << center.getZ() << ";" << endl;
340 }
341
342 }
343
344 if (option.count(summary_t) ||
345 option.count(SUMMARY_t)) {
346
347 const int numberOfStrings = getNumberOfStrings(detector);
348 const int numberOfFloors = getNumberOfFloors (detector);
349 const int numberOfModules = getNumberOfModules(detector);
350 const int numberOfPMTs = getNumberOfPMTs (detector);
351
352 typedef JTOOLS::JRange<int> JRange_t;
353
354 const JRange_t string_t(make_array(detector.begin(), detector.end(), &JModule::getString));
355 const JRange_t floor_t (make_array(detector.begin(), detector.end(), &JModule::getFloor));
356
357 set<int> strings;
358 set<int> modules;
359
360 for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {
361 strings.insert(module->getString());
362 modules.insert(module->getID());
363 }
364
365 if (option.count(summary_t)) {
366
367 cout << "Number of strings = " << setw(4) << numberOfStrings << endl;
368 cout << "Number of floors = " << setw(4) << numberOfFloors << endl;
369 cout << "Number of modules = " << setw(4) << numberOfModules << endl;
370 cout << "Number of PMTs = " << setw(4) << numberOfPMTs << endl;
371 cout << "First string = " << setw(4) << string_t.first << endl;
372 cout << "Last string = " << setw(4) << string_t.second << endl;
373 cout << "First floor = " << setw(4) << floor_t .first << endl;
374 cout << "Last floor = " << setw(4) << floor_t .second << endl;
375
376 } else {
377
378 cout << "let \"NUMBER_OF_STRINGS = " << setw(5) << numberOfStrings << "\";" << endl;
379 cout << "let \"NUMBER_OF_FLOORS = " << setw(5) << numberOfFloors << "\";" << endl;
380 cout << "let \"NUMBER_OF_MODULES = " << setw(5) << numberOfModules << "\";" << endl;
381 cout << "let \"NUMBER_OF_PMTS = " << setw(5) << numberOfPMTs << "\";" << endl;
382 cout << "let \"FIRST_STRING = " << setw(5) << string_t.first << "\";" << endl;
383 cout << "let \"LAST_STRING = " << setw(5) << string_t.second << "\";" << endl;
384 cout << "let \"FIRST_FLOOR = " << setw(5) << floor_t .first << "\";" << endl;
385 cout << "let \"LAST_FLOOR = " << setw(5) << floor_t .second << "\";" << endl;
386 cout << "STRINGS=(";
387 copy(strings.begin(), strings.end(), ostream_iterator<int>(cout, " "));
388 cout << ");" << endl;
389 cout << "MODULES=(";
390 copy(modules.begin(), modules.end(), ostream_iterator<int>(cout, " "));
391 cout << ");" << endl;
392
393 // detector map
394
395 for (const auto i : strings) {
396 cout << "typeset -A " << getString(i) << ";" << endl;
397 }
398 for (const auto& i : detector) {
399 cout << getString(i.getString()) << "[" << i.getFloor() << "]=" << i.getID() << ";" << endl;
400 }
401 }
402 }
403}
JDateAndTime.hh
Date and time functions.
JDetectorSupportkit.hh
Detector support kit.
JDetector.hh
Data structure for detector geometry and calibration.
JConstants.hh
Mathematical constants.
JMessage.hh
General purpose messaging.
FATAL
#define FATAL(A)
Definition JMessage.hh:67
debug
int debug
debug level
Definition JSirene.cc:72
JModule.hh
Data structure for optical module.
JParser.hh
Utility class to parse command line options.
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition JParser.hh:2142
main
int main(int argc, char **argv)
Definition JPrintDetector.cc:70
JPrint.hh
I/O formatting auxiliaries.
MAKE_STRING
#define MAKE_STRING(A)
Make string.
Definition JPrint.hh:63
JRange.hh
Auxiliary class to define a range between two values.
JDETECTOR::JCalibration::getT0
double getT0() const
Get time offset.
Definition JDetector/JCalibration.hh:100
JDETECTOR::JDetectorAddressMap
Lookup table for PMT addresses in detector.
Definition JDetectorAddressMap.hh:26
JDETECTOR::JDetector
Detector data structure.
Definition JDetector.hh:96
JDETECTOR::JModuleAddressMap
Lookup table for PMT addresses in optical module.
Definition JModuleAddressMap.hh:88
JDETECTOR::JModule
Data structure for a composite optical module.
Definition JModule.hh:76
JDETECTOR::JModule::rotate
void rotate(const JRotation3D &R)
Rotate module.
Definition JModule.hh:315
JDETECTOR::JPMTPhysicalAddress
Data structure for PMT physical address.
Definition JPMTPhysicalAddress.hh:31
JDETECTOR::JPMT
Data structure for PMT geometry, calibration and status.
Definition JPMT.hh:49
JGEOMETRY3D::JQuaternion3D
Data structure for unit quaternion in three dimensions.
Definition JQuaternion3D.hh:309
JGEOMETRY3D::JVector3D::getY
double getY() const
Get y position.
Definition JVector3D.hh:104
JGEOMETRY3D::JVector3D::getZ
double getZ() const
Get z position.
Definition JVector3D.hh:115
JGEOMETRY3D::JVector3D::getX
double getX() const
Get x position.
Definition JVector3D.hh:94
JGEOMETRY3D::JVersor3D::getDY
double getDY() const
Get y direction.
Definition JVersor3D.hh:106
JGEOMETRY3D::JVersor3D::getDX
double getDX() const
Get x direction.
Definition JVersor3D.hh:95
JGEOMETRY3D::JVersor3D::getDZ
double getDZ() const
Get z direction.
Definition JVersor3D.hh:117
JLANG::JException
General exception.
Definition JException.hh:24
JLANG::JObjectID::getID
int getID() const
Get identifier.
Definition JObjectID.hh:50
JOSCPROB::JOscProbInterpolator
Template definition of a multi-dimensional oscillation probability interpolation table.
Definition JOscProbInterpolator.hh:71
JPARSER::JParser
Utility class to parse command line options.
Definition JParser.hh:1698
JTOOLS::JRange
Range of values.
Definition JRange.hh:42
JACOUSTICS::getString
JMODEL::JString getString(const JFit &fit)
Get model parameters of string.
Definition JAcoustics/JEvtToolkit.hh:102
JPP
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Definition JAAnetToolkit.hh:45
FILL
Auxiliary data structure for sequence of same character.
Definition JManip.hh:330
FIXED
Auxiliary data structure for floating point format specification.
Definition JManip.hh:448
JAANET::detector
Detector file.
Definition JHead.hh:227
JDETECTOR::JDetectorBuilder
Auxiliary interface for building detector.
Definition JDetectorSupportkit.hh:75
JGEOMETRY3D::JQuaternion3D::decomposition
Auxiliary data structure for decomposition of quaternion in twist and swing quaternions.
Definition JQuaternion3D.hh:1127
JLANG::JStatus::getStatus
status_type getStatus(const JType< status_type > &type) const
Get status.
Definition JStatus.hh:60
JPARSER::initialised
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Definition JParser.hh:68
JSYSTEM::JDateAndTime
Auxiliary class for date and time.
Definition JDateAndTime.hh:81
JSYSTEM::JDateAndTime::toString
std::string toString() const
Get ASCII formatted date and time.
Definition JDateAndTime.hh:253
SCIENTIFIC
Auxiliary data structure for floating point format specification.
Definition JManip.hh:488
Generated by doxygen 1.10.0