Jpp 20.0.0-27-g39925593c-D
the software that should make you happy
Loading...
Searching...
No Matches
JDetector.cc
Go to the documentation of this file.
1
2#include <stdlib.h>
3#include <string>
4#include <iostream>
5#include <sstream>
6#include <fstream>
7#include <deque>
8
9#include "JDetector/JDetector.hh"
10#include "JDetector/JDetectorToolkit.hh"
11#include "JDetector/JDetectorSupportkit.hh"
12#include "JDetector/JDetectorParameters.hh"
13#include "JDetector/JAnchor.hh"
14#include "JSupport/JMeta.hh"
15
16#include "Jeep/JeepToolkit.hh"
17#include "Jeep/JParser.hh"
18#include "Jeep/JMessage.hh"
19
20namespace {
21
22 static const char* const ARCA_t = "ARCA"; // KM3NeT/ARCA detector
23 static const char* const ORCA_t = "ORCA"; // KM3NeT/ORCA detector
24 static const char* const Antares_t = "Antares"; // Antares detector
25 static const char* const NOPMTs_t = "NOPMTs"; // no PMTs detector
26}
27
28
29/**
30 * \file
31 *
32 * Auxiliary program to generate detector files.
33 *
34 * The footprint of the detector can be provided through an input file (option -f).
35 *
36 * When no input file for the footprint is specified,
37 * the footprint is according the so-called WPD reference detector which is subsequently scaled
38 * according the specified line distance (JDETECTOR::JDetectorParameters::lineDistance_m).\n
39 * The specified number of lines (JDETECTOR::JDetectorParameters::numberOfLines) is then selected from
40 * this footprint according the specified detector type (JDETECTOR::JDetectorParameters::detectorType).\n
41 * The list of options includes:
42 * 1. remove outer positions (massive);
43 * 2. remove inner positions (hollow);
44 * \author mdejong
45 */
46int main(int argc, char **argv)
47{
48 using namespace std;
49 using namespace JPP;
50
51 string inputFile;
52 string outputFile;
53 int detectorID;
54 string variant;
55 JDetectorParameters parameters;
56 double Zmin_m;
57 string option;
58 int debug;
59
60 // default values
61
62 parameters.detectorType = 1; //
63 parameters.floorDistance_m = 0.0; // [m]
64 parameters.numberOfFloors = 18; //
65 parameters.lineDistance_m = 0.0; // [m]
66 parameters.numberOfLines = 0; //
67
68 try {
69
70 JParser<> zap("Auxiliary program to generate detector files.");
71
72 JDetectorParametersHelper helper(parameters);
73
74 zap['f'] = make_field(inputFile, "anchor positions") = "";
75 zap['o'] = make_field(outputFile, "detector file");
76 zap['D'] = make_field(detectorID, "detector identifier");
77 zap['V'] = make_field(variant, "detector version") = getDetectorVersions<string>();
78 zap['@'] = make_field(helper, "detector parameters") = JPARSER::initialised();
79 zap['z'] = make_field(Zmin_m, "distance from seabed [m]") = 100.;
80 zap['O'] = make_field(option, "detector type") = ARCA_t, ORCA_t, Antares_t, NOPMTs_t;
81 zap['d'] = make_field(debug) = 0;
82
83 zap(argc, argv);
84 }
85 catch(const exception &error) {
86 FATAL(error.what() << endl);
87 }
88
89
90 JDetectorHeader header;
91
92 if (option == ARCA_t)
93 header = getARCADetectorHeader();
94 else if (option == ORCA_t)
95 header = getORCADetectorHeader();
96 else if (option == Antares_t)
97 ;
98 else if (option == NOPMTs_t)
99 ;
100 else
101 FATAL("Invalid option " << option << endl);
102
103
104 JDetector detector(detectorID, variant, header);
105
106 if (!hasDetectorAddressMap(detector.getID())) {
107 FATAL("No detector address map for detector identifier " << detector.getID() << endl);
108 }
109
110 detector.comment.add(JMeta(argc,argv));
111
112 const JDetectorAddressMap& demo = getDetectorAddressMap(detector.getID());
113
114
115 deque<JAnchor> footprint;
116
117 if (inputFile != "") {
118
119 const string file_name = getFullFilename(LD_LIBRARY_PATH, inputFile);
120
121 NOTICE("Reading file " << file_name << ".. " << flush);
122
123 ifstream in(file_name.c_str());
124
125 if (in) {
126
127 for (JAnchor anchor; in >> anchor && footprint.size() < parameters.numberOfLines; ) {
128 footprint.push_back(anchor);
129 }
130
131 in.close();
132
133 NOTICE("OK" << endl);
134
135 } else {
136
137 FATAL("Error opening file " << inputFile);
138 }
139
140 } else {
141
142 // WPD footprint
143
144 footprint.push_back(JAnchor( 1, 6, -7));
145 footprint.push_back(JAnchor( 2, 188, 9));
146 footprint.push_back(JAnchor( 3, 131, 153));
147 footprint.push_back(JAnchor( 4, -85, 120));
148 footprint.push_back(JAnchor( 5, -175, 39));
149 footprint.push_back(JAnchor( 6, -64, -181));
150 footprint.push_back(JAnchor( 7, 90, -178));
151 footprint.push_back(JAnchor( 8, 368, 20));
152 footprint.push_back(JAnchor( 9, 251, 159));
153 footprint.push_back(JAnchor( 10, 264, -146));
154 footprint.push_back(JAnchor( 11, 189, 293));
155 footprint.push_back(JAnchor( 12, 12, 309));
156 footprint.push_back(JAnchor( 13, -176, 338));
157 footprint.push_back(JAnchor( 14, -266, 149));
158 footprint.push_back(JAnchor( 15, -336, -1));
159 footprint.push_back(JAnchor( 16, -248, -174));
160 footprint.push_back(JAnchor( 17, -196, -348));
161 footprint.push_back(JAnchor( 18, -18, -350));
162 footprint.push_back(JAnchor( 19, 211, -339));
163 footprint.push_back(JAnchor( 20, 541, -12));
164 footprint.push_back(JAnchor( 21, 442, 183));
165 footprint.push_back(JAnchor( 22, 422, -125));
166 footprint.push_back(JAnchor( 23, 364, 287));
167 footprint.push_back(JAnchor( 24, 342, -323));
168 footprint.push_back(JAnchor( 25, 297, 499));
169 footprint.push_back(JAnchor( 26, 51, 452));
170 footprint.push_back(JAnchor( 27, -78, 488));
171 footprint.push_back(JAnchor( 28, -281, 508));
172 footprint.push_back(JAnchor( 29, -325, 304));
173 footprint.push_back(JAnchor( 30, -455, 128));
174 footprint.push_back(JAnchor( 31, -541, 15));
175 footprint.push_back(JAnchor( 32, -429, -142));
176 footprint.push_back(JAnchor( 33, -327, -316));
177 footprint.push_back(JAnchor( 34, -255, -488));
178 footprint.push_back(JAnchor( 35, -103, -508));
179 footprint.push_back(JAnchor( 36, 74, -469));
180 footprint.push_back(JAnchor( 37, 289, -452));
181 footprint.push_back(JAnchor( 38, 715, 0));
182 footprint.push_back(JAnchor( 39, 653, 154));
183 footprint.push_back(JAnchor( 40, 614, -141));
184 footprint.push_back(JAnchor( 41, 530, 269));
185 footprint.push_back(JAnchor( 42, 512, -323));
186 footprint.push_back(JAnchor( 43, 473, 472));
187 footprint.push_back(JAnchor( 44, 458, -473));
188 footprint.push_back(JAnchor( 45, 395, 632));
189 footprint.push_back(JAnchor( 46, 205, 630));
190 footprint.push_back(JAnchor( 47, 1, 667));
191 footprint.push_back(JAnchor( 48, -160, 644));
192 footprint.push_back(JAnchor( 49, -352, 639));
193 footprint.push_back(JAnchor( 50, -413, 460));
194 footprint.push_back(JAnchor( 51, -564, 286));
195 footprint.push_back(JAnchor( 52, -636, 125));
196 footprint.push_back(JAnchor( 53, -717, 23));
197 footprint.push_back(JAnchor( 54, -605, -139));
198 footprint.push_back(JAnchor( 55, -518, -279));
199 footprint.push_back(JAnchor( 56, -482, -448));
200 footprint.push_back(JAnchor( 57, -379, -595));
201 footprint.push_back(JAnchor( 58, -211, -619));
202 footprint.push_back(JAnchor( 59, 8, -652));
203 footprint.push_back(JAnchor( 60, 214, -619));
204 footprint.push_back(JAnchor( 61, 395, -617));
205 footprint.push_back(JAnchor( 62, 886, -12));
206 footprint.push_back(JAnchor( 63, 802, 155));
207 footprint.push_back(JAnchor( 64, 844, -128));
208 footprint.push_back(JAnchor( 65, 697, 285));
209 footprint.push_back(JAnchor( 66, 731, -275));
210 footprint.push_back(JAnchor( 67, 621, 450));
211 footprint.push_back(JAnchor( 68, 600, -497));
212 footprint.push_back(JAnchor( 69, 568, 593));
213 footprint.push_back(JAnchor( 70, 517, -586));
214 footprint.push_back(JAnchor( 71, 430, 776));
215 footprint.push_back(JAnchor( 72, 240, 758));
216 footprint.push_back(JAnchor( 73, 70, 783));
217 footprint.push_back(JAnchor( 74, -75, 803));
218 footprint.push_back(JAnchor( 75, -252, 809));
219 footprint.push_back(JAnchor( 76, -478, 771));
220 footprint.push_back(JAnchor( 77, -510, 590));
221 footprint.push_back(JAnchor( 78, -643, 468));
222 footprint.push_back(JAnchor( 79, -741, 284));
223 footprint.push_back(JAnchor( 80, -845, 175));
224 footprint.push_back(JAnchor( 81, -907, -6));
225 footprint.push_back(JAnchor( 82, -792, -142));
226 footprint.push_back(JAnchor( 83, -737, -340));
227 footprint.push_back(JAnchor( 84, -639, -492));
228 footprint.push_back(JAnchor( 85, -518, -596));
229 footprint.push_back(JAnchor( 86, -455, -751));
230 footprint.push_back(JAnchor( 87, -260, -814));
231 footprint.push_back(JAnchor( 88, -103, -819));
232 footprint.push_back(JAnchor( 89, 106, -764));
233 footprint.push_back(JAnchor( 90, 274, -754));
234 footprint.push_back(JAnchor( 91, 417, -762));
235 footprint.push_back(JAnchor( 92, 1094, -4));
236 footprint.push_back(JAnchor( 93, 1010, 120));
237 footprint.push_back(JAnchor( 94, 992, -171));
238 footprint.push_back(JAnchor( 95, 898, 324));
239 footprint.push_back(JAnchor( 96, 901, -298));
240 footprint.push_back(JAnchor( 97, 768, 471));
241 footprint.push_back(JAnchor( 98, 818, -477));
242 footprint.push_back(JAnchor( 99, 731, 632));
243 footprint.push_back(JAnchor(100, 713, -633));
244 footprint.push_back(JAnchor(101, 625, 772));
245 footprint.push_back(JAnchor(102, 661, -759));
246 footprint.push_back(JAnchor(103, 530, 949));
247 footprint.push_back(JAnchor(104, 345, 960));
248 footprint.push_back(JAnchor(105, 221, 931));
249 footprint.push_back(JAnchor(106, -29, 962));
250 footprint.push_back(JAnchor(107, -155, 967));
251 footprint.push_back(JAnchor(108, -382, 917));
252 footprint.push_back(JAnchor(109, -545, 965));
253 footprint.push_back(JAnchor(110, -634, 779));
254 footprint.push_back(JAnchor(111, -731, 600));
255 footprint.push_back(JAnchor(112, -852, 471));
256 footprint.push_back(JAnchor(113, -903, 322));
257 footprint.push_back(JAnchor(114, -976, 166));
258 footprint.push_back(JAnchor(115, -1100, 11));
259 footprint.push_back(JAnchor(116, -1006, -131));
260 footprint.push_back(JAnchor(117, -867, -332));
261 footprint.push_back(JAnchor(118, -800, -504));
262 footprint.push_back(JAnchor(119, -706, -593));
263 footprint.push_back(JAnchor(120, -604, -800));
264 footprint.push_back(JAnchor(121, -503, -933));
265 footprint.push_back(JAnchor(122, -318, -927));
266 footprint.push_back(JAnchor(123, -177, -926));
267 footprint.push_back(JAnchor(124, -10, -894));
268 footprint.push_back(JAnchor(125, 155, -919));
269 footprint.push_back(JAnchor(126, 381, -967));
270 footprint.push_back(JAnchor(127, 531, -900));
271 footprint.push_back(JAnchor(128, 1091, 311));
272 footprint.push_back(JAnchor(129, 1061, -285));
273 footprint.push_back(JAnchor(130, 966, 481));
274 footprint.push_back(JAnchor(131, 958, -457));
275 footprint.push_back(JAnchor(132, 905, 658));
276 footprint.push_back(JAnchor(133, 863, -630));
277 footprint.push_back(JAnchor(134, 788, 742));
278 footprint.push_back(JAnchor(135, 686, 911));
279 footprint.push_back(JAnchor(136, 277, 1113));
280 footprint.push_back(JAnchor(137, 92, 1094));
281 footprint.push_back(JAnchor(138, -84, 1127));
282 footprint.push_back(JAnchor(139, -302, 1062));
283 footprint.push_back(JAnchor(140, -732, 892));
284 footprint.push_back(JAnchor(141, -793, 761));
285 footprint.push_back(JAnchor(142, -901, 661));
286 footprint.push_back(JAnchor(143, -1012, 456));
287 footprint.push_back(JAnchor(144, -1101, 280));
288 footprint.push_back(JAnchor(145, -1077, -326));
289 footprint.push_back(JAnchor(146, -973, -457));
290 footprint.push_back(JAnchor(147, -883, -602));
291 footprint.push_back(JAnchor(148, -822, -771));
292 footprint.push_back(JAnchor(149, -689, -910));
293 footprint.push_back(JAnchor(150, -457, -1066));
294 footprint.push_back(JAnchor(151, -268, -1067));
295 footprint.push_back(JAnchor(152, -58, -1080));
296 footprint.push_back(JAnchor(153, 72, -1100));
297 footprint.push_back(JAnchor(154, 280, -1081));
298
299 // scale positions to specified distance
300
301 for (deque<JAnchor>::iterator i = footprint.begin(); i != footprint.end(); ++i) {
302 i->mul(parameters.lineDistance_m * 0.5 / 90.0);
303 }
304
305 if (parameters.numberOfLines > footprint.size()) {
306 FATAL("Available footprint inconsistent with specified number of lines "
307 << parameters.numberOfLines << " > " << footprint.size() << endl);
308 }
309
310 // select specified number of lines
311
312 sort(footprint.begin(), footprint.end());
313
314 if (parameters.detectorType == 2) { // hollow detector
315 while (parameters.numberOfLines < footprint.size()) {
316 footprint.pop_front();
317 }
318 } else { // massive detector
319 while (parameters.numberOfLines < footprint.size()) {
320 footprint.pop_back();
321 }
322 }
323
324 // reset object identifiers
325
326 for (deque<JAnchor>::iterator i = footprint.begin(); i != footprint.end(); ++i) {
327 static_cast<JObjectID&>(*i) = JObjectID(distance(footprint.begin(),i) + 1);
328 }
329 }
330
331
332 for (deque<JAnchor>::iterator i = footprint.begin(); i != footprint.end(); ++i) {
333 DEBUG(i->getID() << ' ' << i->getX() << ' ' << i->getY() << endl);
334 }
335
336
337 // build detector
338
339 int PMT = 1;
340
341 for (deque<JAnchor>::const_iterator anchor = footprint.begin(); anchor != footprint.end(); ++anchor) {
342
343 {
344 const JVector3D position(anchor->getX(),
345 anchor->getY(),
346 option == ARCA_t ? ARCA_TBARZ_M :
347 option == ORCA_t ? ORCA_TBARZ_M :
348 0.0);
349
350 const JLocation location(anchor->getID(), 0);
351
352 const int id = demo.getModuleID(location);
353
354 JModule module(id, location);
355
356 module.set(position);
357
358 module.compile();
359
360 detector.push_back(module);
361 }
362
363 for (unsigned int floor = 1; floor <= parameters.numberOfFloors; ++floor) {
364
365 const JVector3D position(anchor->getX(),
366 anchor->getY(),
367 Zmin_m + (floor - 1) * parameters.floorDistance_m);
368
369 const JLocation location(anchor->getID(), floor);
370
371 const int id = demo.getModuleID(location);
372
373 JModule module;
374
375 if (option == ARCA_t || option == ORCA_t)
376 module = getModule<JKM3NeT_t> (id, location);
377 else if (option == Antares_t)
378 module = getModule<JAntares_t>(id, location);
379 else if (option == NOPMTs_t)
380 module = JModule(id, location);
381 else
382 FATAL("Invalid option " << option << endl);
383
384 module.add(position);
385
386 for (JModule::iterator pmt = module.begin(); pmt != module.end(); ++pmt, ++PMT) {
387 pmt->setID(PMT++);
388 //pmt->setCalibration(JCalibration((R + position.getZ())/C));
389 }
390
391 detector.push_back(module);
392 }
393 }
394
395 // Store detector.
396
397 try {
398 store(outputFile, detector);
399 }
400 catch(const JException& error) {
401 FATAL(error);
402 }
403
404 return 0;
405}
PMT
JDAQPMTIdentifier PMT
Command line options.
Definition JDAQTimesliceSelector.cc:35
outputFile
string outputFile
Definition JDAQTimesliceSelector.cc:37
JDetectorSupportkit.hh
Detector support kit.
main
int main(int argc, char **argv)
Definition JDetector.cc:46
JDetector.hh
Data structure for detector geometry and calibration.
JMessage.hh
General purpose messaging.
DEBUG
#define DEBUG(A)
Message macros.
Definition JMessage.hh:62
NOTICE
#define NOTICE(A)
Definition JMessage.hh:64
FATAL
#define FATAL(A)
Definition JMessage.hh:67
debug
int debug
debug level
Definition JSirene.cc:72
JMeta.hh
ROOT I/O of application specific meta data.
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
JeepToolkit.hh
Auxiliary methods for handling file names, type names and environment.
distance
std::vector< T >::difference_type distance(typename std::vector< T >::const_iterator first, typename PhysicsEvent::const_iterator< T > second)
Specialisation of STL distance.
Definition build/antares-dataformat/PhysicsEvent.hh:434
JDETECTOR::JAnchor
Data structure for anchor position on sea bed.
Definition JAnchor.hh:27
JDETECTOR::JDetectorAddressMap
Lookup table for PMT addresses in detector.
Definition JDetectorAddressMap.hh:26
JDETECTOR::JDetectorHeader
Data structure for detector header.
Definition JDetectorHeader.hh:43
JDETECTOR::JDetectorParametersHelper
Auxiliary class for I/O of JDetectorParameters data structure.
Definition JDetectorParameters.hh:102
JDETECTOR::JDetectorParameters
Data structure for parameters for detector geometry.
Definition JDetectorParameters.hh:29
JDETECTOR::JDetectorParameters::numberOfLines
JParameter< unsigned int > numberOfLines
Definition JDetectorParameters.hh:47
JDETECTOR::JDetectorParameters::numberOfFloors
JParameter< unsigned int > numberOfFloors
Definition JDetectorParameters.hh:45
JDETECTOR::JDetector
Detector data structure.
Definition JDetector.hh:96
JDETECTOR::JLocation
Logical location of module.
Definition JLocation.hh:40
JDETECTOR::JModule
Data structure for a composite optical module.
Definition JModule.hh:76
JGEOMETRY3D::JVector3D
Data structure for vector in three dimensions.
Definition JVector3D.hh:36
JLANG::JException
General exception.
Definition JException.hh:24
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
JPP
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Definition JAAnetToolkit.hh:45
JAANET::detector
Detector file.
Definition JHead.hh:227
JPARSER::initialised
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Definition JParser.hh:68
JSUPPORT::JMeta
Auxiliary class for ROOT I/O of application specific meta data.
Definition JMeta.hh:72
Generated by doxygen 1.10.0