Jpp  17.1.0
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JToAshort.cc
Go to the documentation of this file.
1 #include <iostream>
2 #include <iomanip>
3 #include <vector>
4 #include <map>
5 #include <algorithm>
6 
7 #include "TROOT.h"
8 #include "TFile.h"
9 
10 #include "JDB/JToAshort.hh"
11 #include "JDB/JSupport.hh"
12 
13 #include "JLang/JPredicate.hh"
15 
17 
18 #include "Jeep/JPrint.hh"
19 #include "Jeep/JParser.hh"
20 #include "Jeep/JMessage.hh"
21 
22 
23 namespace {
24 
25  using namespace JPP;
26 
27  /**
28  * Auxiliary class to extend JDATABASE::JToAshort.
29  */
30  struct JToAshort_t :
31  public JToAshort
32  {
33  /**
34  * Constructor.
35  *
36  * \param data data
37  * \param counter counter
38  */
39  JToAshort_t( const JToAshort& data, const int counter) :
40  JToAshort(data),
41  counter(counter)
42  {}
43 
44 
45  /**
46  * Write data to output stream.
47  *
48  * \param out output stream
49  * \param object object
50  * \return output stream
51  */
52  friend inline std::ostream& operator<<(std::ostream& out, const JToAshort_t& object)
53  {
54  using namespace std;
55  using namespace JPP;
56 
57  out << setw(6) << object.RUN << ' '
58  << setw(8) << object.counter << ' '
59  << FIXED(20,5) << object.UNIXTIMEBASE << ' '
60  << setw(10) << object.DOMID << ' '
61  << FIXED(9,6) << object.TOA_S << ' '
62  << FIXED(5,0) << object.QUALITYFACTOR;
63 
64  return out;
65  }
66 
67  int counter;
68 
69 
70  /**
71  * Compare times of arrival.
72  *
73  * \param first first time-of-arrival
74  * \param second second time-of-arrival
75  * \return true if first time-of-arrival less than second; else false
76  */
77  friend inline bool operator<(const JToAshort_t& first, const JToAshort_t& second)
78  {
79  return (first.UNIXTIMEBASE + first.TOA_S < second.UNIXTIMEBASE + second.TOA_S);
80  }
81  };
82 }
83 
84 
85 /**
86  * \file
87  *
88  * Example program to test acoustic data.
89  * \author mdejong
90  */
91 int main(int argc, char **argv)
92 {
93  using namespace std;
94  using namespace JPP;
95 
97  JLimit_t& numberOfEvents = inputFile.getLimit();
98  double precision;
99  bool allowed;
100  int debug;
101 
102  try {
103 
104  JParser<> zap("Example program to test acoustic data.");
105 
106  zap['f'] = make_field(inputFile);
107  zap['n'] = make_field(numberOfEvents) = JLimit::max();
108  zap['A'] = make_field(allowed);
109  zap['e'] = make_field(precision) = 1.0e-7;
110  zap['d'] = make_field(debug) = 2;
111 
112  zap(argc, argv);
113  }
114  catch(const exception &error) {
115  FATAL(error.what() << endl);
116  }
117 
118 
119  typedef vector<JToAshort_t> buffer_type; // acoustic data type
120 
121  map<int, map< int, buffer_type > > data; // emitter -> receiver -> data
122 
123  for (int counter = 0; inputFile.hasNext(); ++counter) {
124 
125  STATUS("counter: " << setw(8) << counter << '\r' << flush); DEBUG(endl);
126 
127  const JToAshort* parameters = inputFile.next();
128 
129  data[parameters->EMITTERID][parameters->DOMID].push_back(JToAshort_t(*parameters, counter));
130  }
131  STATUS(endl);
132 
133 
134  for (map<int, map< int, buffer_type> >::iterator i = data.begin(); i != data.end(); ++i) {
135 
136  for (map< int, buffer_type>::iterator module = i->second.begin(); module != i->second.end(); ++module) {
137 
138  buffer_type& buffer = module->second;
139 
140  sort(buffer.begin(), buffer.end());
141 
142  for (buffer_type::const_iterator p = buffer.begin(); p != buffer.end(); ++p) {
143  for (buffer_type::const_iterator q = buffer.begin(); q != p; ++q) {
144 
145  if (labs(p->DOMID - q->DOMID) == 0 &&
146  fabs(p->QUALITYFACTOR - q->QUALITYFACTOR) <= precision &&
147  fabs((p->UNIXTIMEBASE + p->TOA_S) -
148  (q->UNIXTIMEBASE + q->TOA_S)) <= precision) {
149 
150  if (!allowed || (p->TOA_S < TOAMAX_S &&
151  q->TOA_S < TOAMAX_S)) {
152 
153  cout << *p << endl;
154  cout << *q << endl;
155 
156  cout << "Difference between absolute times "
157  << SCIENTIFIC(12,3) << ((p->UNIXTIMEBASE + p->TOA_S) -
158  (q->UNIXTIMEBASE + q->TOA_S)) << endl;
159  }
160  }
161  }
162  }
163  }
164  }
165 }
Utility class to parse command line options.
Definition: JParser.hh:1500
int main(int argc, char *argv[])
Definition: Main.cc:15
bool operator<(const Head &first, const Head &second)
Less than operator.
Definition: JHead.hh:1722
#define STATUS(A)
Definition: JMessage.hh:63
*fatal Wrong number of arguments esac JCookie sh typeset Z DETECTOR typeset Z SOURCE_RUN typeset Z TARGET_RUN set_variable PARAMETERS_FILE $WORKDIR parameters
Definition: diff-Tuna.sh:38
static counter_type max()
Get maximum counter value.
Definition: JLimit.hh:117
Auxiliary data structure for floating point format specification.
Definition: JManip.hh:446
then echo The file $DIR KM3NeT_00000001_00000000 root already please rename or remove it first
Auxiliary class for defining the range of iterations of objects.
Definition: JLimit.hh:41
static const double TOAMAX_S
Maximal allowed time-of-arrival [s].
I/O formatting auxiliaries.
Acoustics support kit.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Definition: JParser.hh:1961
ROOT TTree parameter settings.
int debug
debug level
Definition: JSirene.cc:67
General purpose messaging.
#define FATAL(A)
Definition: JMessage.hh:67
Scanning of objects from multiple files according a format that follows from the extension of each fi...
General purpose class for object reading from a list of file names.
Utility class to parse command line options.
std::ostream & operator<<(std::ostream &stream, const CLBCommonHeader &header)
const JLimit & getLimit() const
Get limit.
Definition: JLimit.hh:73
Auxiliary data structure for floating point format specification.
Definition: JManip.hh:484
int EMITTERID
waveform identifier
Definition: JToAshort.hh:27
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62