Jpp 19.3.0-rc.1
the software that should make you happy
Loading...
Searching...
No Matches
Functions
JDetectorDB.cc File Reference

Auxiliary program to download or update detector from database. More...

#include <string>
#include <iostream>
#include <iomanip>
#include <vector>
#include <set>
#include <map>
#include <algorithm>
#include <cmath>
#include "JDB/JDB.hh"
#include "JDB/JSelector.hh"
#include "JDB/JSelectorSupportkit.hh"
#include "JDB/JCLBMap.hh"
#include "JDB/JDetectorIntegration_t.hh"
#include "JDB/JPBSSequence.hh"
#include "JDB/JProductRouter.hh"
#include "JDB/JDBToolkit.hh"
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JDetector/JDetectorSupportkit.hh"
#include "JDetector/JDetectorCalibration.hh"
#include "JLang/JEquation_t.hh"
#include "JLang/JComparator.hh"
#include "JLang/JPredicate.hh"
#include "JLang/JVectorize.hh"
#include "JSupport/JMeta.hh"
#include "JSon/JSupport.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 download or update detector from database.

Author
mdejong

Definition in file JDetectorDB.cc.

Function Documentation

◆ main()

int main ( int argc,
char ** argv )

Definition at line 74 of file JDetectorDB.cc.

75{
76 using namespace std;
77 using namespace JPP;
78
79 JServer server;
80 string usr;
81 string pwd;
82 string cookie;
83 string id;
84 vector<JEquation_t> calset;
85 string detectorFile;
86 string outputFile;
87 int run;
88 string variant;
89 JCounter wizzard;
90 double precision;
91 int debug;
92
93 try {
94
95 JParser<> zap("Auxiliary program to download or update detector from database.");
96
97 zap['s'] = make_field(server) = getServernames();
98 zap['u'] = make_field(usr) = "";
99 zap['!'] = make_field(pwd) = "";
100 zap['C'] = make_field(cookie) = "";
101 zap['D'] = make_field(id, "detector identifier") = "";
102 zap['r'] = make_field(run, "run number") = -1;
103 zap['@'] = make_field(calset, "calibration sets \"<key> = <value>[; <key> = <value>]\"" << endl
104 << "possible keys: " << JEEPZ() << make_array(getCalibrationType.begin(),
105 getCalibrationType.end(),
106 &JCalibrationType::nick_name) << endl
107 << "possible values: <calibration identifier> or \"?\"") = JPARSER::initialised();
108 zap['a'] = make_field(detectorFile) = "";
109 zap['o'] = make_field(outputFile) = "";
110 zap['V'] = make_field(variant, "detector version") = getDetectorVersions<string>(), "";
111 zap['W'] = make_field(wizzard, "-W add base modules; -WW set time offsets of acoustic sensors");
112 zap['p'] = make_field(precision) = 1.0e-5;
113 zap['d'] = make_field(debug) = 2;
114
115 zap(argc, argv);
116 }
117 catch(const exception &error) {
118 FATAL(error.what() << endl);
119 }
120
121
122 map<string, string> calibration;
123
124 for (vector<JEquation_t>::const_iterator i = calset.begin(); i != calset.end(); ++i) {
125
126 if (getCalibrationType.hasNickname(i->getKey()))
127 FATAL("Invalid calibration set \"" << i->getKey() << "\"" << endl);
128 else
129 calibration[i->getKey()] = i->getValue();
130 }
131
132 if (calibration.size() != 0u &&
133 calibration.size() != 3u &&
134 calibration.size() != 6u) {
135 FATAL("Inconsistent specificaton of calibration sets, either {}, {t,p,r}cal or {t,p,r,a,c,s}cal.");
136 }
137
138 const bool file = (detectorFile != "");
139 const bool tpr = (calibration.count(TCAL) != 0 &&
140 calibration.count(PCAL) != 0 &&
141 calibration.count(RCAL) != 0 &&
142 calibration.count(ACAL) == 0 &&
143 calibration.count(CCAL) == 0 &&
144 calibration.count(SCAL) == 0);
145 const bool tpracs = (calibration.count(TCAL) != 0 &&
146 calibration.count(PCAL) != 0 &&
147 calibration.count(RCAL) != 0 &&
148 calibration.count(ACAL) != 0 &&
149 calibration.count(CCAL) != 0 &&
150 calibration.count(SCAL) != 0);
151 const bool daq = (run != -1);
152
153 if (file && id != "") { FATAL("Inconsistent specificaton of detector file and detector identifier."); }
154 if (file && tpr) { FATAL("Inconsistent specificaton of detector file and {t,p,r}cal calibration sets."); }
155 if (file && tpracs) { FATAL("Inconsistent specificaton of detector file and {t,p,r,a,c,s}cal calibration sets."); }
156 if (file && daq) { FATAL("Inconsistent specificaton of detector file and run number."); }
157
158 if (tpr && id == "") { FATAL("Inconsistent specificaton of {t,p,r}cal calibration sets and no detector identifier."); }
159 if (tpr && variant != "") { FATAL("Inconsistent specificaton of {t,p,r}cal calibration sets and variant."); }
160 if (tpr && daq) { FATAL("Inconsistent specificaton of {t,p,r}cal calibration sets and run number."); }
161
162 if (tpracs && id == "") { FATAL("Inconsistent specificaton of {t,p,r,a,c,s}cal calibration sets and no detector identifier."); }
163 if (tpracs && variant == "") { FATAL("Inconsistent specificaton of {t,p,r,a,c,s}cal calibration sets and no variant."); }
164 if (tpracs && daq) { FATAL("Inconsistent specificaton of {t,p,r,a,c,s}cal calibration sets and run number."); }
165
166 if (daq && id == "") { FATAL("Inconsistent specificaton of run number and no detector identifier."); }
167
168 if (wizzard &&
169 variant != "" &&
170 getDetectorVersion(variant) < JDetectorVersion::V4) {
171 FATAL("Inconsistent detector variant and option -W.");
172 }
173
174 try {
175 JDB::reset(usr, pwd, cookie);
176 }
177 catch(const exception& error) {
178 FATAL(error.what() << endl);
179 }
180
181 JDetector detector;
182
183 try {
184
185 if (file)
186 load(detectorFile, detector);
187 else if (tpr)
188 *(JDB::get()->DetX)(id.c_str(),
189 calibration[TCAL] != "?" ? calibration[TCAL].c_str() : NULL,
190 calibration[PCAL] != "?" ? calibration[PCAL].c_str() : NULL,
191 calibration[RCAL] != "?" ? calibration[RCAL].c_str() : NULL) >> detector;
192 else if (tpracs)
193 *(JDB::get()->DetX)(id.c_str(),
194 calibration[TCAL] != "?" ? calibration[TCAL].c_str() : NULL,
195 calibration[PCAL] != "?" ? calibration[PCAL].c_str() : NULL,
196 calibration[RCAL] != "?" ? calibration[RCAL].c_str() : NULL,
197 calibration[ACAL] != "?" ? calibration[ACAL].c_str() : NULL,
198 calibration[CCAL] != "?" ? calibration[CCAL].c_str() : NULL,
199 calibration[SCAL] != "?" ? calibration[SCAL].c_str() : NULL,
200 getDetectorVersion(variant)) >> detector;
201 else if (daq && variant == "")
202 *(JDB::get()->DetX)(id.c_str(), run) >> detector;
203 else if (daq && variant != "")
204 *(JDB::get()->DetX)(id.c_str(), run, getDetectorVersion(variant)) >> detector;
205 else
206 FATAL("Invalid options.");
207 }
208 catch(const exception& error) {
209 FATAL(error.what() << endl);
210 }
211
212
213 detector.comment.add(JMeta(argc,argv));
214
215 if (detector.setToLatestVersion()) {
216 NOTICE("Set detector version to " << detector.getVersion() << endl);
217 }
218
219
220 if (wizzard) {
221
222 vector<JCLBMap> clbmap;
223
224 try {
225
226 JSelector selector = getSelector<JCLBMap>(getDetector(detector.getID()));
227
228 selector.add(&JCLBMap::FLOORID, 0);
229
230 ResultSet& rs = getResultSet(getTable<JCLBMap>(), selector);
231
232 rs >> clbmap;
233
234 rs.Close();
235 }
236 catch(const exception& error) {
237 FATAL(error.what() << endl);
238 }
239
240 map<int, map<int, JPosition3D> > zmap;
241
242 for (JDetector::const_iterator module = detector.begin(); module != detector.end(); ++module) {
243 zmap[module->getString()][module->getFloor()] = module->getPosition();
244 }
245
246 for (vector<JCLBMap>::const_iterator i = clbmap.begin(); i != clbmap.end(); ++i) {
247
248 if (i->FLOORID == 0) {
249
250 JDetector::iterator p = find_if(detector.begin(), detector.end(), make_predicate(&JModule::getID, i->DOMID, JComparison::eq()));
251
252 if (p == detector.end()) {
253
254 if (zmap[i->DUID].count(1) != 0) {
255
256 JModule module(i->DOMID, JLocation(i->DUID, i->FLOORID));
257
258 module.setPosition(JPosition3D(zmap[i->DUID][1].getX(),
259 zmap[i->DUID][1].getY(),
260 isARCADetector(detector) ? ARCA_TBARZ_M :
261 isORCADetector(detector) ? ORCA_TBARZ_M :
262 0.0));
263
264 NOTICE("Adding module " << module << endl);
265
266 detector.push_back(module);
267
268 } else {
269
270 ERROR("Missing floor in string " << FILL(4,'0') << i->DUID << FILL() << endl);
271 }
272
273 } else if (zmap[i->DUID].count(1) != 0) {
274
275 if (fabs(p->getX() - zmap[i->DUID][1].getX()) > precision ||
276 fabs(p->getY() - zmap[i->DUID][1].getY()) > precision) {
277
278 NOTICE("Offset module " << getLabel(p->getLocation()) << " (x,y) position "
279 << SCIENTIFIC(12,3) << (p->getX() - zmap[i->DUID][1].getX()) << ' '
280 << SCIENTIFIC(12,3) << (p->getY() - zmap[i->DUID][1].getY()) << endl);
281
282 p->setPosition(JPosition3D(zmap[i->DUID][1].getX(),
283 zmap[i->DUID][1].getY(),
284 p->getZ()));
285 }
286
287 } else {
288
289 FATAL("Missing floor in string " << setw(4) << i->DUID << endl);
290 }
291 }
292 }
293 }
294
295 if (wizzard > 1) {
296
297 try {
298
299 JDetectorIntegration_t detector_t;
300
301 ResultSet& rs = getResultSet(getTable<JDetectorIntegration_t>());
302
303 if (! (rs >> detector_t)) {
304 THROW(JDatabaseException, "Error reading " << getTable<JDetectorIntegration_t>());
305 }
306
307 detector_t.configure(getDetector(detector.getID()));
308
309 for (const JPBS_t& pbs : { PBS::ACOUSTIC_SENSOR, PBS::HYDROPHONE }) {
310
311 const JProductRouter router(detector_t, getPBSSequences(pbs));
312
313 JDetectorIntegration_t::range_type range = detector_t.find(pbs);
314
315 for (JDetectorIntegration_t::range_const_iterator i = range.first; i != range.second; ++i) {
316
317 const JLocation_t location = router.getLocation(i->first);
318
319 JDetector::iterator p = find_if(detector.begin(),
320 detector.end(),
321 make_predicate(&JModule::getLocation,
322 JLocation(location.string, (pbs != PBS::HYDROPHONE ? location.floor : 0)),
323 JComparison::eq()));
324
325 if (p != detector.end()) {
326
327 try {
328
329 const double t0 = getAverage(make_array(p->begin(), p->end(), &JPMT::getT0), 0.0) - getDelayTime(i->first);
330
331 if (t0 != p->getT0()) {
332
333 NOTICE("Set module time calibration " << getLabel(p->getLocation()) << ' ' << setw(10) << p->getID() << ' ' << setw(28) << left << i->first << right << ' ' << FIXED(12,3) << t0 << endl);
334
335 p->setCalibration(t0);
336 }
337 }
338 catch(const exception& error) {
339 if (i->first.getVariant() != NOPIEZO) {
340 ERROR(location << ' ' << error.what() << endl);
341 }
342 }
343
344 } else {
345
346 ERROR("No module for UPI " << i->first << " at " << location << endl);
347 }
348 }
349 }
350 }
351 catch(const exception& error) {
352 FATAL(error.what() << endl);
353 }
354 }
355
356 if (outputFile != "")
357 store(outputFile, detector);
358 else
359 cout << detector;
360}
outputFile
string outputFile
Definition JDAQTimesliceSelector.cc:37
THROW
#define THROW(JException_t, A)
Marco for throwing exception with std::ostream compatible message.
Definition JException.hh:712
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
make_field
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition JParser.hh:2142
JDATABASE::JSelector
Auxiliary class for specifying selection of database data.
Definition JDB/JSelector.hh:71
JDATABASE::JSelector::add
JSelector & add(const JSelector &selection)
Add selection.
Definition JDB/JSelector.hh:134
JDETECTOR::JCalibration::getT0
double getT0() const
Get time offset.
Definition JDetector/JCalibration.hh:100
JDETECTOR::JDetector
Detector data structure.
Definition JDetector.hh:96
JDETECTOR::JLocation
Logical location of module.
Definition JLocation.hh:40
JDETECTOR::JLocation::getLocation
const JLocation & getLocation() const
Get location.
Definition JLocation.hh:70
JDETECTOR::JModule
Data structure for a composite optical module.
Definition JModule.hh:75
JGEOMETRY3D::JPosition3D
Data structure for position in three dimensions.
Definition JPosition3D.hh:38
JLANG::JDatabaseException
Database exception.
Definition JException.hh:684
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
JDATABASE::PBS::HYDROPHONE
static const JPBS_t HYDROPHONE(4, 5)
PBS of hydrophone
JDATABASE::PBS::ACOUSTIC_SENSOR
static const JPBS_t ACOUSTIC_SENSOR(3, 4, 3, 6, 2)
PBS of piezo sensor
JDATABASE::getResultSet
ResultSet & getResultSet(const std::string &query)
Get result set.
Definition JDB.hh:438
JDATABASE::getPBSSequences
static JGetPBSSequences getPBSSequences
Function object to get PBS sequences as a function of PBS.
Definition JPBSSequence.hh:276
JDATABASE::getTable
const char * getTable()
Get table name.
Definition JDB.hh:373
JDATABASE::getDetector
JDetectorsHelper & getDetector()
Auxiliary function for helper object initialisation.
Definition JDBToolkit.hh:378
JDATABASE::getServernames
std::vector< JServer > getServernames()
Get list of names of available database servers.
Definition JDB.hh:108
JDETECTOR::getLabel
std::string getLabel(const JLocation &location)
Get module label for monitoring and other applications.
Definition JLocation.hh:247
JDETECTOR::getDetectorVersion
static const JGetDetectorVersion getDetectorVersion
Function object to map detector version to numerical value.
Definition JDetectorVersion.hh:311
JDETECTOR::load
void load(const std::string &file_name, JDetector &detector)
Load detector from input file.
Definition JDetectorToolkit.hh:508
JDETECTOR::store
void store(const std::string &file_name, const JDetector &detector)
Store detector to output file.
Definition JDetectorToolkit.hh:605
JDETECTOR::getDetectorVersions
std::vector< T > getDetectorVersions()
Get detector versions.
Definition JDetectorVersion.hh:354
JLANG::make_array
const array_type< JValue_t > & make_array(const JValue_t(&array)[N])
Method to create array of values.
Definition JVectorize.hh:54
JMATH::getAverage
std::iterator_traits< T >::value_type getAverage(T __begin, T __end)
Get average.
Definition JMath.hh:494
JPP
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Definition JAAnetToolkit.hh:43
JSON::TCAL
static const std::string TCAL
PMT time offsets.
Definition JSon/JSupport.hh:32
JSON::PCAL
static const std::string PCAL
(optical|base) module positions
Definition JSon/JSupport.hh:33
JSON::SCAL
static const std::string SCAL
(module|PMT) status
Definition JSon/JSupport.hh:37
JSON::RCAL
static const std::string RCAL
optical module orientations
Definition JSon/JSupport.hh:34
JSON::ACAL
static const std::string ACAL
acoustic time offsets (piezo sensor or hydrophone)
Definition JSon/JSupport.hh:35
JSON::CCAL
static const std::string CCAL
compass alignment (a.k.a. quaternion calibration)
Definition JSon/JSupport.hh:36
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::calibration
Calibration.
Definition JHead.hh:330
JAANET::detector
Detector file.
Definition JHead.hh:227
JDATABASE::JDetectorIntegration_t
Detector integration.
Definition JDetectorIntegration_t.hh:189
JDATABASE::JDetectorIntegration_t::range_const_iterator
map_type::const_iterator range_const_iterator
Definition JDetectorIntegration_t.hh:195
JDATABASE::JDetectorIntegration_t::configure
void configure(const std::string &detid, const bool option=false)
Configure detector integration for given detector identifier.
Definition JDetectorIntegration_t.hh:219
JDATABASE::JDetectorIntegration_t::find
range_type find(const JUPI_t &upi) const
Find range of products with given UPI.
Definition JDetectorIntegration_t.hh:492
JDATABASE::JLocation_t
Auxiliary data structure for location of product in detector.
Definition JLocation_t.hh:26
JDATABASE::JLocation_t::string
int string
position in detector
Definition JLocation_t.hh:105
JDATABASE::JLocation_t::floor
int floor
position in string
Definition JLocation_t.hh:106
JDATABASE::JPBS_t
Product breakdown structure (PBS).
Definition JPBS_t.hh:29
JDATABASE::JProductRouter
Auxiliary class to map UPI to location in detector.
Definition JProductRouter.hh:27
JDATABASE::JServer
Wrapper class for server name.
Definition JDB.hh:54
JDATABASE::getSelector
Template definition for getting table specific selector.
Definition JSelectorSupportkit.hh:35
JDETECTOR::JDetectorVersion::V4
@ V4
Version with quaternion and time offset per module.
Definition JDetectorVersion.hh:52
JEEP::JEEPZ
Auxiliary data structure for streaming of STL containers.
Definition JStreamToolkit.hh:606
JPARSER::initialised
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Definition JParser.hh:68
JSON::JCalibrationType::nick_name
std::string nick_name
nick name (e.g. "xcal")
Definition JSon/JSupport.hh:51
JSON::getCalibrationType
Auxiliary data structure for correspondence between nick and full name of calibration types.
Definition JSon/JSupport.hh:61
JSON::getCalibrationType::hasNickname
bool hasNickname(const std::string &type) const
Has calibration type.
Definition JSon/JSupport.hh:82
JSUPPORT::JMeta
Auxiliary class for ROOT I/O of application specific meta data.
Definition JMeta.hh:72
SCIENTIFIC
Auxiliary data structure for floating point format specification.
Definition JManip.hh:488
Generated by doxygen 1.12.0