63 int main(
int argc,
char **argv)
88 JParser<> zap(
"Main program to trigger events in acoustic data.");
90 zap[
'f'] =
make_field(inputFile,
"output of JConvertDB -q toashort");
91 zap[
'n'] =
make_field(numberOfEvents) = JLimit::max();
95 zap[
'a'] =
make_field(detectorFile,
"detector file");
100 zap[
'p'] =
make_field(precision,
"precision time-of-arrival") = 1.0e-6;
105 catch(
const exception &error) {
106 FATAL(error.what() << endl);
131 if (
i->getFloor() == 0) {
138 catch(
const exception&) {
144 i->getPosition() + pos,
145 i->getT0() * 1.0e-9);
149 for (tripods_container::const_iterator
i =
tripods.begin();
i !=
tripods.end(); ++
i) {
151 i->getUTMPosition() -
detector.getUTMPosition());
154 for (transmitters_container::const_iterator
i = transmitters.begin();
i != transmitters.end(); ++
i) {
159 catch(
const exception&) {
178 public map<int, JQuantile>
180 long long int sum()
const
182 long long int count = 0;
184 for (const_iterator
i = this->begin();
i != this->end(); ++
i) {
185 count +=
i->second.getCount();
192 map_type number_of_entries;
193 map_type number_of_aliens;
194 map_type number_of_triggers;
195 map_type number_of_events;
205 while (inputFile.hasNext()) {
207 if (inputFile.getCounter()%10000 == 0) {
208 STATUS(
"counter: " << setw(8) << inputFile.getCounter() <<
'\r' << flush);
DEBUG(endl);
214 oid = parameters->
DETID;
217 if (oid != parameters->
DETID) {
218 FATAL(
"Invalid detector identifier " << parameters->
DETID <<
" != " << oid << endl);
225 number_of_entries[parameters->
EMITTERID].put(1);
227 if (emitters.has(
id) && receivers.has(parameters->
DOMID)) {
234 catch(
const exception&) {
235 number_of_aliens[parameters->
EMITTERID].put(1);
240 for (map_type::const_iterator
i = number_of_entries.begin();
i != number_of_entries.end(); ++
i) {
241 STATUS(
"number of entries: " << setw(3) <<
i->first <<
' ' << setw(8) <<
i->second.getCount() << endl);
244 for (map_type::const_iterator
i = number_of_aliens.begin();
i != number_of_aliens.end(); ++
i) {
245 STATUS(
"number of aliens: " << setw(3) <<
i->first <<
' ' << setw(8) <<
i->second.getCount() << endl);
260 buffer_type::iterator __end = unique(receiver->second.begin(), receiver->second.end(),
JTransmission::equals(precision));
264 for (buffer_type::const_iterator p = receiver->second.begin(); p != __end; ++p) {
266 buffer.push_back(*p);
271 sort(buffer.begin(), buffer.end());
277 for (buffer_type::const_iterator p = buffer.begin(); p != buffer.end(); ++p) {
279 buffer_type::const_iterator q = p;
283 while (++q != buffer.end() && q->getToE() - p->getToE() <=
parameters.TMax_s) {}
289 JEvent event(oid, number_of_triggers.sum(),
i->first, p, q);
291 number_of_triggers[
i->first].put(event.size());
293 DEBUG(
"trigger: " << endl << event);
299 }
else if (match(out, event)) {
316 STATUS(
"number of toes: " << setw(3) <<
i->first <<
' ' << setw(8) << buffer.size());
317 if (Q.getCount() > 0) {
319 ' ' <<
FIXED(6,1) << Q.getXmin() <<
320 ' ' <<
FIXED(6,1) << Q.getXmax());
324 for (JTriggerOutput::const_iterator event = data.begin();
event != data.end(); ++event) {
326 number_of_events[
event->getID()].put(event->size());
333 for (map_type::const_iterator
i = number_of_triggers.begin();
i != number_of_triggers.end(); ++
i) {
334 STATUS(
"number of triggers: " << setw(3) <<
i->first <<
' ' << setw(8) <<
i->second.getCount() << endl);
337 for (map_type::const_iterator
i = number_of_events.begin();
i != number_of_events.end(); ++
i) {
338 STATUS(
"number of events: " << setw(3) <<
i->first <<
' ' << setw(8) <<
i->second.getCount() <<
' ' <<
FIXED(7,3) <<
i->second.getMean() << endl);
Utility class to parse command line options.
Q(UTCMax_s-UTCMin_s)-livetime_s
JPredicate< JResult_t T::*, JComparison::eq > make_predicate(JResult_t T::*member, const JResult_t value)
Helper method to create predicate for data member.
int main(int argc, char *argv[])
JContainer< std::vector< JTransmitter > > transmitters_container
void merge(const JEvent &event)
Merge event.
std::vector< T >::difference_type distance(typename std::vector< T >::const_iterator first, typename PhysicsEvent::const_iterator< T > second)
Specialisation of STL distance.
General purpose class for hash map of unique elements.
ROOT TTree parameter settings.
Recording of objects on file according a format that follows from the file name extension.
std::string DETID
constraint
*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
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Auxiliary data structure for floating point format specification.
V(JDAQEvent-JTriggerReprocessor)*1.0/(JDAQEvent+1.0e-10)
Data structure for detector geometry and calibration.
JMODEL::JString getString(const JFit &fit)
Get model parameters of string.
Data structure for hydrophone.
const JPolynome f1(1.0, 2.0, 3.0)
Function.
Auxiliary class for defining the range of iterations of objects.
static const JSoundVelocity getSoundVelocity(1541.0,-17.0e-3,-2000.0)
Function object for velocity of sound.
I/O formatting auxiliaries.
JContainer< std::vector< JHydrophone > > hydrophones_container
Data structure for transmitter.
Auxiliary wrapper for I/O of container with optional comment (see JComment).
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
int getID() const
Get identifier.
JPosition3D getPosition(const Vec &pos)
Get position.
JTransmission getTransmission(const JToAshort &data, const JAbstractSoundVelocity &V) const
Get transmission.
ROOT TTree parameter settings.
JContainer< std::vector< JTripod > > tripods_container
General purpose messaging.
Implementation for depth dependend velocity of sound.
Scanning of objects from multiple files according a format that follows from the extension of each fi...
void load(const std::string &file_name, JDetector &detector)
Load detector from input file.
Acoustic trigger parameters.
Auxiliary class to define a range between two values.
General purpose class for object reading from a list of file names.
Utility class to parse command line options.
static const int PIEZO_DISABLE
Enable (disable) use of piezo if this status bit is 0 (1);.
Auxiliary class to compare transmissions.
static const int HYDROPHONE_DISABLE
Enable (disable) use of hydrophone if this status bit is 0 (1);.
Data structure for position in three dimensions.
const JLimit & getLimit() const
Get limit.
bool equals(const JFirst_t &first, const JSecond_t &second, const double precision=std::numeric_limits< double >::min())
Check equality.
static JEmitterID getEmitterID
Function object for emitter identification.
do set_variable DETECTOR_TXT $WORKDIR detector
Data structure for tripod.
Match of two events considering overlap in time.
int EMITTERID
waveform identifier
#define DEBUG(A)
Message macros.