Jpp test-rotations-old
the software that should make you happy
Loading...
Searching...
No Matches
JSeabed.cc File Reference

Auxiliary application to apply tilt angles to seabed. More...

#include <string>
#include <iostream>
#include <iomanip>
#include <vector>
#include "JDetector/JTripod.hh"
#include "JDetector/JDetector.hh"
#include "JDetector/JDetectorToolkit.hh"
#include "JSupport/JMeta.hh"
#include "Jeep/JContainer.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 application to apply tilt angles to seabed.

Author
mdejong

Definition in file software/JAcoustics/JSeabed.cc.

Function Documentation

◆ main()

int main ( int argc,
char ** argv )

Definition at line 22 of file software/JAcoustics/JSeabed.cc.

23{
24 using namespace std;
25 using namespace JPP;
26
28
29 string detectorFile;
30 string tripodsFile;
31 double Tx;
32 double Ty;
33 int debug;
34
35 try {
36
37 JParser<> zap("Auxiliary application to apply tilt angles to seabed.");
38
39 zap['a'] = make_field(detectorFile, "detector file.");
40 zap['T'] = make_field(tripodsFile, "tripods file.");
41 zap['x'] = make_field(Tx, "tilt angle x-z [rad]");
42 zap['y'] = make_field(Ty, "tilt angle y-z [rad]");
43 zap['d'] = make_field(debug) = 1;
44
45 zap(argc, argv);
46 }
47 catch(const exception& error) {
48 FATAL(error.what() << endl);
49 }
50
51
52 tripods_container tripods;
53
54 try {
55 tripods.load(tripodsFile.c_str());
56 }
57 catch(const exception& error) {
58 FATAL(error.what() << endl);
59 }
60
61
63
64 try {
65 load(detectorFile, detector);
66 }
67 catch(const exception& error) {
68 FATAL(error.what() << endl);
69 }
70
71
72 JUTMPosition pos = detector.getUTMPosition();
73
74 if (!tripods.empty()) {
75
76 pos = JUTMPosition(0.0, 0.0, 0.0);
77
78 for (tripods_container::const_iterator i = tripods.begin(); i != tripods.end(); ++i) {
79 pos += i->getUTMPosition();
80 }
81
82 pos /= tripods.size();
83 }
84
85
86 detector.comment.add(JMeta(argc, argv));
87
88 for (JDetector::iterator i = detector.begin(); i != detector.end(); ++i) {
89
90 const double dz = Tx * (i->getX() + detector.getX() - pos.getX()) + Ty * (i->getY() + detector.getY() - pos.getY());
91
92 *i += JPosition3D(0.0, 0.0, dz);
93 }
94
95 store(detectorFile, detector);
96
97
98 tripods.comment.add(JMeta(argc, argv));
99
100 for (tripods_container::iterator i = tripods.begin(); i != tripods.end(); ++i) {
101
102 const double dz = Tx * (i->getX() - pos.getX()) + Ty * (i->getY() - pos.getY());
103
104 *i += JUTMPosition(0.0, 0.0, dz);
105 }
106
107 tripods.store(tripodsFile.c_str());
108
109 return 0;
110}
#define FATAL(A)
Definition JMessage.hh:67
int debug
debug level
Definition JSirene.cc:72
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition JParser.hh:2142
Detector data structure.
Definition JDetector.hh:96
Data structure for position in three dimensions.
Utility class to parse command line options.
Definition JParser.hh:1698
Data structure for UTM position.
double getY() const
Get y.
double getX() const
Get x.
const JUTMPosition & getUTMPosition() const
Get UTM position.
void load(const std::string &file_name, JDetector &detector)
Load detector from input file.
void store(const std::string &file_name, const JDetector &detector)
Store detector to output file.
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Detector file.
Definition JHead.hh:227
JComment & add(const std::string &comment)
Add comment.
Definition JComment.hh:100
Auxiliary wrapper for I/O of container with optional comment (see JComment).
Definition JContainer.hh:42
void store(const char *file_name) const
Store to output file.
void load(const char *file_name)
Load from input file.
Auxiliary class for ROOT I/O of application specific meta data.
Definition JMeta.hh:72