Jpp 20.0.0-rc.2
the software that should make you happy
Loading...
Searching...
No Matches
Functions
JPrintDetector.cc File Reference

Auxiliary program to print detector file in human friendly format. More...

#include <string>
#include <iostream>
#include <iomanip>
#include <set>
#include <iterator>
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JDetector/JDetectorAddressMap.hh"
#include "JDetector/JDetectorSupportkit.hh"
#include "JDetector/JModule.hh"
#include "JMath/JConstants.hh"
#include "JTools/JRange.hh"
#include "JSystem/JDateAndTime.hh"
#include "JGeometry3D/JGeometry3DToolkit.hh"
#include "JGeometry3D/JCylinder3D.hh"
#include "JGeometry3D/JQuaternion3D.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 print detector file in human friendly format.

Author
mdejong

Definition in file JPrintDetector.cc.

Function Documentation

◆ main()

int main ( int argc,
char ** argv )

Definition at line 70 of file JPrintDetector.cc.

71{
72 using namespace std;
73 using namespace JPP;
74
75 string detectorFile;
76 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) =
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;
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 == default_t) {
129
130 cout << detector << endl;
131
132 } else if (option == modules_t) {
133
134 const JDetectorBuilder& demo = getDetectorBuilder(detector.getID());
135
136 for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {
137
138 cout << "Module";
139 cout << ' ' << noshowpos << setw(8) << right << module->getID();
140 cout << ' ' << noshowpos << setw(3) << right << module->getString();
141 cout << ' ' << noshowpos << setw(2) << right << module->getFloor();
142 cout << ' ' << FIXED(7,2) << module->getX();
143 cout << ' ' << FIXED(7,2) << module->getY();
144 cout << ' ' << FIXED(7,2) << module->getZ();
145 cout << ' ' << FIXED(8,2) << module->getT0();
146
147 if (module->getFloor() != 0) {
148 {
149 JModule buffer = demo.getModule(module->getID(), module->getLocation());
150
151 const JQuaternion3D Q = getRotation(buffer, *module);
152 const JQuaternion3D::decomposition q1(Q, JVector3Z_t);
153
154 const double phi = (JVector3Z_t.getDot(q1.twist) >= 0.0 ? +1.0 : -1.0) * q1.twist.getAngle();
155
156 cout << ' ' << FIXED(7,2) << phi;
157 }
158 {
159 const JQuaternion3D Q = module->getQuaternion();
160 const JQuaternion3D::decomposition q1(Q, JVector3Z_t);
161
162 const double phi = (JVector3Z_t.getDot(q1.twist) >= 0.0 ? +1.0 : -1.0) * q1.twist.getAngle();
163
164 cout << ' ' << FIXED(7,2) << phi;
165 }
166 {
167 JModule buffer = demo.getModule(module->getID(), module->getLocation());
168
169 const JRotation3D R = getRotation(buffer, *module);
170
171 buffer.rotate(R);
172
173 cout << ' ' << (JModule::compare(buffer, *module, precision) ? "==" : "!=")
174 << ' ' << (demo.get(module->getID()) == demo.getDefaultModuleAddressMap() ? "default" : "custom");
175 }
176 }
177
178 cout << endl;
179 }
180
181 } else if (option == pmts_t) {
182
183 const JDetectorAddressMap& demo = getDetectorAddressMap(detector.getID());
184
185 for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {
186
187 const JModuleAddressMap& memo = demo.get(module->getID());
188
189 for (size_t i = 0; i != module->size(); ++i) {
190
191 const JPMT& pmt = module->getPMT(i);
192 const JPMTPhysicalAddress& address = memo.getPMTPhysicalAddress(i);
193
194 cout << "PMT[" << FILL(2,'0') << i << FILL() << "]";
195 cout << ' ' << address;
196 cout << ' ' << noshowpos << setw(8) << right << pmt.getID();
197 cout << ' ' << FIXED(7,2) << pmt.getX();
198 cout << ' ' << FIXED(7,2) << pmt.getY();
199 cout << ' ' << FIXED(7,2) << pmt.getZ();
200 cout << ' ' << FIXED(6,3) << pmt.getDX();
201 cout << ' ' << FIXED(6,3) << pmt.getDY();
202 cout << ' ' << FIXED(6,3) << pmt.getDZ();
203 cout << ' ' << FIXED(8,2) << pmt.getT0();
204 cout << ' ' << noshowpos << setw(8) << right << pmt.getStatus();
205 cout << endl;
206 }
207 }
208
209 } else if (option == geometry_t ||
210 option == GEOMETRY_t) {
211
212 if (option == geometry_t) {
213
214 cout << "Maximal distance [m] = " << FIXED(6,1) << getMaximalDistance(detector) << endl;
215 cout << "Minimal distance [m] = " << FIXED(6,1) << getMinimalDistance(detector) << endl;
216 cout << "Maximal time [ns] = " << FIXED(6,1) << getMaximalTime (detector) << endl;
217
218 } else {
219
220 cout << "set_variable MAXIMAL_DISTANCE_M " << FIXED(6,1) << getMaximalDistance(detector) << ";" << endl;
221 cout << "set_variable MINIMAL_DISTANCE_M " << FIXED(6,1) << getMinimalDistance(detector) << ";" << endl;
222 }
223
224 } else if (option == comment_t) {
225
226 cout << detector.comment << endl;
227
228 } else if (option == header_t ||
229 option == HEADER_t) {
230
231 if (option == header_t) {
232
233 cout << detector.getProperties() << endl;
234 cout << "name = " << (isARCADetector(detector) ? "ARCA" :
235 isORCADetector(detector) ? "ORCA" :
236 "unknown") << endl;
237 cout << "validity = "
238 << JDateAndTime(detector.getLowerLimit(), true).toString() << ' '
239 << JDateAndTime(detector.getUpperLimit(), true).toString() << endl;
240
241 } else {
242
243 cout << "set_variable UTM_EAST " << FIXED(12,2) << detector.getUTMEast() << ";" << endl;
244 cout << "set_variable UTM_NORTH " << FIXED(12,2) << detector.getUTMNorth() << ";" << endl;
245 cout << "set_variable UTM_Z " << FIXED(12,2) << detector.getUTMZ() << ";" << endl;
246 cout << "set_variable UTM_ZONE " << FIXED(12,2) << detector.getUTMZone() << ";" << endl;
247 cout << "set_variable UTM_WGS " << FIXED(12,2) << detector.getWGS() << ";" << endl;
248 }
249
250 } else if (option == version_t ||
251 option == VERSION_t) {
252
253 if (option == version_t) {
254
255 cout << detector.getVersion() << endl;
256
257 } else {
258
259 cout << "set_variable DETECTOR_VERSION " << detector.getVersion() << ";" << endl;
260 }
261
262 } else if (option == identifier_t ||
263 option == IDENTIFIER_t) {
264
265 if (option == identifier_t) {
266
267 cout << "Detector " << detector.getID() << endl;
268
269 } else {
270
271 cout << "set_variable DETECTOR_ID " << detector.getID() << ";" << endl;
272 }
273
274 } else if (option == can_t ||
275 option == CAN_t) {
276
277 const JCylinder3D cylinder(detector.begin(), detector.end());
278
279 const double V = (cylinder.getZmax() - cylinder.getZmin()) * PI * cylinder.getRadius() * cylinder.getRadius();
280 const double D = getMaximalDistance(detector);
281
282 if (option == can_t) {
283
284 cout << "X = " << FIXED(7,1) << cylinder.getX() << endl;
285 cout << "Y = " << FIXED(7,1) << cylinder.getY() << endl;
286 cout << "Zmin = " << FIXED(7,1) << cylinder.getZmin() << endl;
287 cout << "Zmax = " << FIXED(7,1) << cylinder.getZmax() << endl;
288 cout << "Radius = " << FIXED(7,1) << cylinder.getRadius() << endl;
289 cout << "Depth = " << FIXED(7,1) << detector.getUTMZ() << endl;
290 cout << "Volume = " << SCIENTIFIC(12,3) << V << endl;
291 cout << "Distance = " << FIXED(9,3) << D << endl;
292
293 } else {
294
295 cout << "set_variable CAN_X_M " << FIXED(7,1) << cylinder.getX() << ";" << endl;
296 cout << "set_variable CAN_Y_M " << FIXED(7,1) << cylinder.getY() << ";" << endl;
297 cout << "set_variable CAN_ZMIN_M " << FIXED(7,1) << cylinder.getZmin() << ";" << endl;
298 cout << "set_variable CAN_ZMAX_M " << FIXED(7,1) << cylinder.getZmax() << ";" << endl;
299 cout << "set_variable CAN_RADIUS_M " << FIXED(7,1) << cylinder.getRadius() << ";" << endl;
300 cout << "set_variable CAN_DEPTH_M " << FIXED(7,1) << detector.getUTMZ() << ";" << endl;
301 cout << "set_variable CAN_VOLUME_M3 " << SCIENTIFIC(12,3) << V << ";" << endl;
302 cout << "set_variable CAN_DISTANCE_M " << FIXED(9,3) << D << ";" << endl;
303 }
304
305 } else if (option == center_t ||
306 option == CENTER_t) {
307
308 const JCenter3D center(detector.begin(), detector.end());
309
310 if (option == center_t) {
311
312 cout << "center = ";
313 cout << showpos << FIXED(8,3) << center.getX() << ' ';
314 cout << showpos << FIXED(8,3) << center.getY() << ' ';
315 cout << showpos << FIXED(8,3) << center.getZ() << endl;
316
317 } else {
318
319 cout << "set_variable CENTER_X_M " << FIXED(7,1) << center.getX() << ";" << endl;
320 cout << "set_variable CENTER_Y_M " << FIXED(7,1) << center.getY() << ";" << endl;
321 cout << "set_variable CENTER_Z_M " << FIXED(7,1) << center.getZ() << ";" << endl;
322 }
323
324 } else if (option == summary_t ||
325 option == SUMMARY_t) {
326
327 const int numberOfStrings = getNumberOfStrings(detector);
328 const int numberOfFloors = getNumberOfFloors (detector);
329 const int numberOfModules = getNumberOfModules(detector);
330 const int numberOfPMTs = getNumberOfPMTs (detector);
331
332 typedef JTOOLS::JRange<int> JRange_t;
333
334 const JRange_t string_t(make_array(detector.begin(), detector.end(), &JModule::getString));
335 const JRange_t floor_t (make_array(detector.begin(), detector.end(), &JModule::getFloor));
336
337 set<int> strings;
338 set<int> modules;
339
340 for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {
341 strings.insert(module->getString());
342 modules.insert(module->getID());
343 }
344
345 if (option == summary_t) {
346
347 cout << "Number of strings = " << setw(4) << numberOfStrings << endl;
348 cout << "Number of floors = " << setw(4) << numberOfFloors << endl;
349 cout << "Number of modules = " << setw(4) << numberOfModules << endl;
350 cout << "Number of PMTs = " << setw(4) << numberOfPMTs << endl;
351 cout << "First string = " << setw(4) << string_t.first << endl;
352 cout << "Last string = " << setw(4) << string_t.second << endl;
353 cout << "First floor = " << setw(4) << floor_t .first << endl;
354 cout << "Last floor = " << setw(4) << floor_t .second << endl;
355
356 } else {
357
358 cout << "let \"NUMBER_OF_STRINGS = " << setw(5) << numberOfStrings << "\";" << endl;
359 cout << "let \"NUMBER_OF_FLOORS = " << setw(5) << numberOfFloors << "\";" << endl;
360 cout << "let \"NUMBER_OF_MODULES = " << setw(5) << numberOfModules << "\";" << endl;
361 cout << "let \"NUMBER_OF_PMTS = " << setw(5) << numberOfPMTs << "\";" << endl;
362 cout << "let \"FIRST_STRING = " << setw(5) << string_t.first << "\";" << endl;
363 cout << "let \"LAST_STRING = " << setw(5) << string_t.second << "\";" << endl;
364 cout << "let \"FIRST_FLOOR = " << setw(5) << floor_t .first << "\";" << endl;
365 cout << "let \"LAST_FLOOR = " << setw(5) << floor_t .second << "\";" << endl;
366 cout << "STRINGS=(";
367 copy(strings.begin(), strings.end(), ostream_iterator<int>(cout, " "));
368 cout << ");" << endl;
369 cout << "MODULES=(";
370 copy(modules.begin(), modules.end(), ostream_iterator<int>(cout, " "));
371 cout << ");" << endl;
372
373 // detector map
374
375 for (const auto i : strings) {
376 cout << "typeset -A " << getString(i) << ";" << endl;
377 }
378 for (const auto& i : detector) {
379 cout << getString(i.getString()) << "[" << i.getFloor() << "]=" << i.getID() << ";" << endl;
380 }
381 }
382 }
383}
FATAL
#define FATAL(A)
Definition JMessage.hh:67
debug
int debug
debug level
Definition JSirene.cc:72
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition JParser.hh:2142
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::JDetectorAddressMap::get
const JModuleAddressMap & get(const int id) const
Get module address map.
Definition JDetectorAddressMap.hh:69
JDETECTOR::JDetectorAddressMap::getDefaultModuleAddressMap
virtual const JModuleAddressMap & getDefaultModuleAddressMap() const =0
Get default module address map.
JDETECTOR::JDetector
Detector data structure.
Definition JDetector.hh:96
JDETECTOR::JModuleAddressMap
Lookup table for PMT addresses in optical module.
Definition JModuleAddressMap.hh:88
JDETECTOR::JModuleAddressMap::getPMTPhysicalAddress
const JPMTPhysicalAddress & getPMTPhysicalAddress(const int tdc) const
Get PMT physical address.
Definition JModuleAddressMap.hh:269
JDETECTOR::JModule
Data structure for a composite optical module.
Definition JModule.hh:75
JDETECTOR::JModule::rotate
void rotate(const JRotation3D &R)
Rotate module.
Definition JModule.hh:314
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
JPARSER::JParser
Utility class to parse command line options.
Definition JParser.hh:1698
JTOOLS::JRange
Range of values.
Definition JRange.hh:42
JAANET::copy
void copy(const Head &from, JHead &to)
Copy header from from to to.
Definition JHead.cc:163
JACOUSTICS::getString
JMODEL::JString getString(const JFit &fit)
Get model parameters of string.
Definition JAcoustics/JEvtToolkit.hh:102
JCALIBRATE::string_t
const char *const string_t
routing by string
Definition JCalibrateMuon.hh:173
JDETECTOR::getNumberOfPMTs
int getNumberOfPMTs(const JModule &module)
Get number of PMTs.
Definition JDetectorToolkit.hh:202
JDETECTOR::getNumberOfFloors
int getNumberOfFloors(const JDetector &detector)
Get number of floors.
Definition JDetectorToolkit.hh:271
JDETECTOR::getDetectorBuilder
JDetectorBuilder & getDetectorBuilder()
Get detector builder.
Definition JDetectorSupportkit.hh:528
JDETECTOR::load
void load(const std::string &file_name, JDetector &detector)
Load detector from input file.
Definition JDetectorToolkit.hh:536
JDETECTOR::isORCADetector
bool isORCADetector(const JDetectorHeader &header)
Check if given detector header is compatible with that of ORCA.
Definition JDetectorHeader.hh:309
JDETECTOR::isARCADetector
bool isARCADetector(const JDetectorHeader &header)
Check if given detector header is compatible with tat of ARCA.
Definition JDetectorHeader.hh:297
JDETECTOR::getDetectorAddressMap
JDetectorAddressMap & getDetectorAddressMap()
Get detector address map.
Definition JDetectorSupportkit.hh:542
JDETECTOR::getNumberOfModules
int getNumberOfModules(const JDetector &detector, const bool option=false)
Get number of modules.
Definition JDetectorToolkit.hh:314
JDETECTOR::getMinimalDistance
double getMinimalDistance(const JDetector &detector)
Get minimal distance between modules in detector.
Definition JDetectorToolkit.hh:108
JDETECTOR::getMaximalDistance
double getMaximalDistance(const JDetector &detector, const bool option=false)
Get maximal distance between modules in detector.
Definition JDetectorToolkit.hh:78
JPP
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Definition JAAnetToolkit.hh:45
KM3NETDAQ::getMaximalTime
double getMaximalTime(const double R_Hz)
Get maximal time for given rate.
Definition JHighRateVeto.hh:35
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
JDETECTOR::JDetectorBuilder::getModule
const JModule & getModule(const int id=-1, const JLocation &location=JLocation()) const
Get module.
Definition JDetectorSupportkit.hh:85
JEEP::JStatus::getStatus
int getStatus() const
Get status.
Definition JStatus.hh:63
JGEOMETRY3D::JQuaternion3D::decomposition
Auxiliary data structure for decomposition of quaternion in twist and swing quaternions.
Definition JQuaternion3D.hh:1127
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.12.0