39 static const int WILDCARD = -1;
41 static const std::string reset_t =
"reset";
42 static const std::string set_t =
"set";
43 static const std::string add_t =
"add";
44 static const std::string sub_t =
"sub";
45 static const std::string rot_t =
"rot";
46 static const std::string mul_t =
"mul";
47 static const std::string tilt_t =
"tilt";
49 static const std::string assign_t =
"assign";
50 static const std::string string_t =
"string";
52 static const std::string rand_t =
"rand";
53 static const std::string randset_t = rand_t + set_t;
54 static const std::string randadd_t = rand_t + add_t;
55 static const std::string randsub_t = rand_t + sub_t;
56 static const std::string randrot_t = rand_t + rot_t;
57 static const std::string randmul_t = rand_t + mul_t;
58 static const std::string randtilt_t = rand_t + tilt_t;
60 static const std::string RESET_t =
"RESET";
61 static const std::string SET_t =
"SET";
62 static const std::string ADD_t =
"ADD";
63 static const std::string SUB_t =
"SUB";
64 static const std::string ROT_t =
"ROT";
74 inline bool apply(
JModule& module,
const std::string& action)
76 if (action == reset_t) {
82 for (JModule::iterator pmt = module.begin(); pmt != module.end(); ++pmt) {
87 }
else if (action == RESET_t) {
92 }
else if (action == SET_t) {
118 inline bool apply(
JModule& module,
const std::string& action,
const double value)
120 if (action == set_t) {
124 }
else if (action == add_t) {
128 }
else if (action == sub_t) {
132 }
else if (action == rot_t) {
142 }
else if (action == mul_t) {
150 }
else if (action == SET_t) {
154 }
else if (action == ADD_t) {
158 }
else if (action == SUB_t) {
162 }
else if (action == ROT_t) {
166 }
else if (action == randadd_t) {
168 module.
add(gRandom->Gaus(0.0, value));
170 }
else if (action == randsub_t) {
172 module.
sub(gRandom->Gaus(0.0, value));
174 }
else if (action == randrot_t){
184 }
else if (action == randmul_t) {
195 }
else if (action == assign_t) {
197 module.
setID((
int) value);
199 }
else if (action == string_t) {
221 inline bool apply(
JModule& module,
const std::string& action,
const double first,
const double second)
223 if (action == tilt_t) {
225 const double Tx =
first;
226 const double Ty = second;
227 const double Tz = sqrt(1.0 - Tx*Tx - Ty*Ty);
229 const double x = Tx * module.
getZ() + module.
getX();
230 const double y = Ty * module.
getZ() + module.
getY();
231 const double z = Tz * module.
getZ();
252 inline bool apply(
JModule& module,
const std::string& action,
const JVector3D& pos)
255 gRandom->Gaus(0.0, pos.
getY()),
256 gRandom->Gaus(0.0, pos.
getZ()));
260 else if (action == add_t)
262 else if (action == sub_t)
264 else if (action == randset_t)
266 else if (action == randadd_t)
268 else if (action == randsub_t)
289 else if (action == ADD_t)
291 else if (action == SUB_t)
308 inline bool apply(
JPMT& pmt,
const std::string& action,
const std::string& value)
314 else if (action == reset_t)
321 catch(
const std::exception&) {
345 return (action !=
"" && !data.empty());
355 bool apply(
JModule& module)
const
357 switch (data.size()) {
360 return ::apply(module, action);
363 return ::apply(module, action, data[0]);
366 return ::apply(module, action, data[0], data[1]);
369 return ::apply(module, action,
JVector3D(data[0], data[1], data[2]));
372 return ::apply(module, action,
JQuaternion3D(data[0], data[1], data[2], data[3]));
387 friend inline std::istream&
operator>>(std::istream&
in, JModifier& modifier)
389 if (in >> modifier.action) {
391 modifier.data.clear();
393 for (
double x; in >> x; ) {
394 modifier.data.push_back(x);
397 in.clear(std::ios_base::eofbit);
411 friend inline std::ostream&
operator<<(std::ostream& out,
const JModifier& modifier)
413 out << modifier.action;
438 inline const JModifier& getModifier(
const int id,
const JModifier& modifier)
444 const string::size_type pos = modifier.action.find(rand_t);
446 if (pos != string::npos) {
448 JModifier&
result = buffer[id][modifier.action][modifier.data.size()];
450 if (!result.is_valid()) {
452 result.action = modifier.action.substr(pos + rand_t.length());
454 for (
size_t i = 0; i != modifier.data.size(); ++i) {
455 result.data.push_back(gRandom->Gaus(0.0, modifier.data[i]));
471 struct JPMTModifier {
485 bool apply(
JPMT& pmt)
const
487 return ::apply(pmt, action, value);
498 friend inline std::istream&
operator>>(std::istream& in, JPMTModifier& modifier)
500 return in >> modifier.action >> modifier.value;
511 friend inline std::ostream&
operator<<(std::ostream& out,
const JPMTModifier& modifier)
513 out << modifier.action;
515 out << modifier.value;
537 static const char SEPARATOR =
'-';
554 friend inline std::istream&
operator>>(std::istream& in, JRange_t&
range)
556 if (in >> range.first) {
558 range.second = range.first;
560 if (in.peek() == (int) JRange_t::SEPARATOR) {
583 friend inline std::ostream&
operator<<(std::ostream& out,
const JRange_t& range)
585 return out << range.first << JRange_t::SEPARATOR << range.second;
597 inline void print(std::ostream& out,
const JModule& module,
const JModifier& modifier)
602 out <<
"Modifier" <<
' '
604 << setw(10) << module.
getID() <<
' '
605 <<
"action" <<
' ' << modifier.action <<
JEEPZ() << modifier.data << endl;
616 inline void print(std::ostream& out,
const JPMTIdentifier& pmt,
const JPMTModifier& modifier)
621 out <<
"PMT modifier" <<
' '
622 <<
"(" << setw(10) << pmt.
getID() <<
"," << setw(2) << pmt.
getPMTAddress() <<
")" <<
' '
623 <<
"action" <<
' ' << modifier.action <<
' '
624 <<
"value" <<
' ' << modifier.value << endl;
724 int main(
int argc,
char **argv)
749 JParser<> zap(
"Auxiliary program to modify detector.");
767 catch(
const exception &error) {
768 FATAL(error.what() << endl);
784 if ((keep.empty() ? 0 : 1 +
785 rm .empty() ? 0 : 1 +
786 id .empty() ? 0 : 1) > 1) {
787 FATAL(
"Use either option -k, -r or -m." << endl);
793 if (
detector.setToLatestVersion()) {
794 NOTICE(
"Set detector version to " <<
detector.getVersion() << endl);
808 istringstream
is(*i);
819 for (JDetector::iterator module =
detector.begin(); module !=
detector.end(); ) {
821 bool __rm__ = !keep.empty() && rm.empty();
836 __rm__ = find(
id.begin(),
id.end(), module->
getID()) ==
id.end();
839 const auto range = del.equal_range(module->
getString());
841 for (
auto i = range.first; i != range.second; ++i) {
842 if (i->second == module->
getFloor()) {
856 for (JDetector::iterator module =
detector.begin(); module !=
detector.end(); ++module) {
858 if (module->
getID() == i->first || i->first == WILDCARD ){
861 print(cout, *module, i->second);
864 if (!i->second.apply(*module)) {
865 ERROR(
"No valid action: " << i->first <<
' ' << i->second << endl);
874 for (JDetector::iterator module =
detector.begin(); module !=
detector.end(); ++module) {
876 if (module->
getString() == i->first || i->first == WILDCARD) {
878 const JModifier modifier = getModifier(module->
getString(), i->second);
881 print(cout, *module, i->second);
884 if (!modifier.apply(*module)) {
885 ERROR(
"No valid action: " << i->first <<
' ' << i->second << endl);
894 for (JDetector::iterator module =
detector.begin(); module !=
detector.end(); ++module) {
898 if (module->
getString() == i->first || i->first == WILDCARD) {
900 const JModifier modifier = getModifier(module->
getString(), i->second);
903 print(cout, *module, i->second);
906 if (!modifier.apply(*module)) {
907 ERROR(
"No valid action: " << i->first <<
' ' << i->second << endl);
917 for (JDetector::iterator module =
detector.begin(); module !=
detector.end(); ++module) {
919 if (module->
getID() == i->first.getModuleID() || i->first.getModuleID() == WILDCARD) {
922 print(cout, i->first, i->second);
925 if (i->first.getPMTAddress() == WILDCARD) {
928 if (!i->second.apply(module->
getPMT(pmt))) {
929 ERROR(
"No valid action: " << i->first <<
' ' << i->second << endl);
933 }
else if (i->first.getPMTAddress() < 0 ||
935 !i->second.apply(module->
getPMT(i->first.getPMTAddress()))) {
936 ERROR(
"No valid action: " << i->first <<
' ' << i->second << endl);
946 FATAL(
"Invalid detector identifier " <<
detector.getID() << endl);
951 for (JDetector::iterator module =
detector.begin(); module !=
detector.end(); ++module) {
960 print(cout,
id, i->second);
963 if (!i->second.apply(module->
getPMT(
id.getPMTAddress()))) {
964 ERROR(
"No valid action: " << i->first <<
' ' << i->second << endl);
void set(const int bit)
Set PMT status.
Utility class to parse command line options.
static const JGetPMTStatusBit getPMTStatusBit
Function object to map key to PMT status bit.
Q(UTCMax_s-UTCMin_s)-livetime_s
int main(int argc, char *argv[])
int getFloor() const
Get floor number.
Data structure for a composite optical module.
static const double HYDROPHONE_DELAYTIME_US
Hydrophone delay time [us].
JQuaternion3D getConjugate() const
Get conjugate of this quaternion.
void setLocation(const JLocation &location)
Set location.
bool hasDetectorAddressMap(const int id)
Check if detector address map is available.
Utility class to parse parameter values.
std::iterator_traits< T >::value_type getAverage(T __begin, T __end)
Get average.
int getNumberOfPMTs(const JModule &module)
Get number of PMTs.
void subT0(const double t0)
Subtract time offset.
Auxiliary class for controlling PMT status.
Lookup table for PMT addresses in detector.
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Data structure for time calibration.
void reset(const int bit)
Reset PMT status.
Data structure for detector geometry and calibration.
const JModuleAddressMap & get(const int id) const
Get module address map.
then echo The file $DIR KM3NeT_00000001_00000000 root already please rename or remove it first
const JQuaternion3D & getQuaternion() const
Get quaternion.
JModule & sub(const JVector3D &pos)
Subtract position.
Lookup table for PMT addresses in optical module.
I/O formatting auxiliaries.
static const double PIEZO_DELAYTIME_US
Piezo delay time [us].
Data structure for vector in three dimensions.
void setQuaternion(const JQuaternion3D &quaternion)
Set quaternion.
Logical location of module.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
const array_type< JValue_t > & make_array(const JValue_t(&array)[N])
Method to create array of values.
int getID() const
Get identifier.
This class represents a rotation around the z-axis.
bool is_valid(const json &js)
Check validity of JSon data.
Data structure for PMT geometry, calibration and status.
void store(const std::string &file_name, const JDetector &detector)
Store detector to output file.
JPosition3D getPosition(const Vec &pos)
Get position.
int getPMTAddress() const
Get PMT address (= TDC).
double getY() const
Get y position.
Auxiliary data structure for streaming of STL containers.
const JPosition3D & getPosition() const
Get position.
JDetectorAddressMap & getDetectorAddressMap()
Get detector address map.
const JPMT & getPMT(const int index) const
Get PMT.
void rotate(const JRotation3D &R)
Rotate module.
General purpose messaging.
Auxiliary data structure for sequence of same character.
Data structure for unit quaternion in three dimensions.
Direct access to module in detector data structure.
z range($ZMAX-$ZMIN)< $MINIMAL_DZ." fi fi typeset -Z 4 STRING typeset -Z 2 FLOOR JPlot1D -f $
int getString() const
Get string number.
std::istream & operator>>(std::istream &in, JAANET::JHead &header)
Read header from input.
const JPMTAddressTranslator & getAddressTranslator(const int tdc) const
Get PMT address translator.
void load(const std::string &file_name, JDetector &detector)
Load detector from input file.
void setT0(const double t0)
Set time offset.
Auxiliary class to define a range between two values.
Utility class to parse command line options.
Wrapper class around string.
void setCalibration(const JCalibration &cal)
Set calibration.
void setID(const int id)
Set identifier.
double getX() const
Get x position.
Base class for data structures with artithmetic capabilities.
void addT0(const double t0)
Add time offset.
Data structure for position in three dimensions.
Data structure for PMT physical address.
do set_variable DETECTOR_TXT $WORKDIR detector
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2] eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY source JAcoustics sh $DETECTOR_ID CHECK_EXIT_CODE typeset A TRIPODS get_tripods $WORKDIR tripod txt TRIPODS for EMITTER in
JModule & set(const JVector3D &pos)
Set position.
double getZ() const
Get z position.
JModule & add(const JVector3D &pos)
Add position.
double getT0() const
Get time offset.