59 int main(
int argc,
char **argv)
75 transmitters_container transmitters;
76 hydrophones_container hydrophones;
79 disable_container disable;
85 JParser<> zap(
"Application to fit position calibration model to acoustic data.");
87 zap[
'f'] =
make_field(inputFile,
"output of JAcousticEventBuilder[.sh]");
89 zap[
'n'] =
make_field(numberOfEvents) = JLimit::max();
97 zap[
'u'] =
make_field(unify,
"unify weighing of pings");
100 "Precede name of data structure by a '+' or '-' "
101 "to add or remove data types in the output, respectively."
107 catch(
const exception &error) {
108 FATAL(error.what() << endl);
128 for (JDetector::const_iterator i =
detector.begin(); i !=
detector.end(); ++i) {
129 receivers[i->getID()] = i->getLocation();
132 for (tripods_container::const_iterator i =
tripods.begin(); i !=
tripods.end(); ++i) {
133 emitters[i->getID()] =
JEmitter(i->getID(),
134 i->getUTMPosition() -
detector.getUTMPosition());
137 for (transmitters_container::const_iterator i = transmitters.begin(); i != transmitters.end(); ++i) {
139 emitters[i->getID()] =
JEmitter(i->getID(),
142 catch(
const exception&) {
163 TH1D h0(
"chi2/NDF", NULL, 50000, 0.0, 1000.0);
164 TH1D h1(
"h1", NULL, 51, -0.5, 50.5);
165 TH1D hn(
"hn", NULL, 100, 0.0, 6.0);
169 range.getLength() + 1,
range.getLowerLimit() - 0.5,
range.getUpperLimit() + 0.5));
173 range.getLength() + 1,
range.getLowerLimit() - 0.5,
range.getUpperLimit() + 0.5));
175 for (Int_t i = 1; i <= HA->GetXaxis()->GetNbins(); ++i) {
176 HA->GetXaxis()->SetBinLabel(i,
MAKE_CSTRING(geometry.at(i-1).first));
177 HB->GetXaxis()->SetBinLabel(i,
MAKE_CSTRING(geometry.at(i-1).first));
185 for (
const string& file_name : inputFile) {
195 else if (oid != evt->
getOID())
196 FATAL(
"Invalid detector identifier " << evt->
getOID() <<
" != " << oid << endl);
199 zmap[evt->begin()->getToE()] = file_name;
208 inputFile.push_back(i->second);
218 int counter[] = { 0, 0 };
220 typedef deque<JEvent> buffer_type;
222 for (buffer_type zbuf; inputFile.hasNext(); ) {
224 STATUS(inputFile.getFilename() <<
'\r');
DEBUG(endl);
228 for (
const string file_name = inputFile.getFilename(); inputFile.hasNext() && file_name == inputFile.getFilename(); ) {
230 const JEvent* evt = inputFile.next();
232 if (emitters.has(evt->
getID())) {
233 zbuf.push_back(*evt);
237 sort(zbuf.begin(), zbuf.end());
239 for (buffer_type::iterator p = zbuf.begin(), q; p != zbuf.end(); p = q) {
241 for (q = p; ++q != zbuf.end() && q->begin()->getToE() <= p->rbegin()->getToE() +
parameters.Tmax_s; ) {}
243 if (q == zbuf.end()) {
245 if (inputFile.hasNext()) {
247 zbuf.erase(zbuf.begin(), p);
259 for (buffer_type::const_iterator i = p; i != q; ++i) {
260 numberOfPings[i->getID()] += 1;
264 DEBUG(
"Number of pings " << setw(2) << i->first <<
' ' << setw(3) << i->second << endl);
267 int minimum_number_of_pings = numeric_limits<int>::max();
270 minimum_number_of_pings = min(minimum_number_of_pings, i->second);
275 for (buffer_type::iterator evt = p; evt != q; ++evt) {
277 sort(evt->begin(), evt->end(),
compare);
281 const JEmitter& emitter = emitters[evt->getID()];
282 const double signal = (unify ? (double) minimum_number_of_pings / (
double) numberOfPings[evt->getID()] : 1.0);
284 for (JEvent::const_iterator i = evt->begin(); i != __end; ++i) {
289 if (receivers.has(i->getID()) && geometry.
hasLocation(receivers[i->getID()]) && i->getQ() >=
parameters.Qmin) {
293 receivers[i->getID()],
300 for (data_type::const_iterator hit = data.begin(); hit != data.end(); ++hit) {
301 HA[hit->getID()]->Fill(geometry.
getIndex(hit->getString()), hit->getFloor(), 1.0);
306 const auto result = katoomba(data.begin(), data.end());
308 for (data_type::const_iterator hit =
result.begin; hit !=
result.end; ++hit) {
309 HB[hit->getID()]->Fill(geometry.
getIndex(hit->getString()), hit->getFloor(), 1.0);
314 cout <<
"result:" <<
' '
318 for (data_type::const_iterator hit =
result.begin; hit !=
result.end; ++hit) {
343 if (selection.is_valid<
JEvent>()) {
345 for (buffer_type::iterator i = p; i != q; ++i) {
351 for (JEvent::iterator hit = out.begin(); hit != out.end(); ++hit) {
372 STATUS(
"Number of events written / rejected: " << counter[0] <<
" / " << counter[1] << endl);
Worker class for fit function of acoustic model.
Utility class to parse command line options.
double getValue(const JScale_t scale)
Get numerical value corresponding to scale.
int main(int argc, char *argv[])
JComparator< JResult_t T::*, JComparison::lt > make_comparator(JResult_t T::*member)
Helper method to create comparator between values of data member.
static JDetectorMechanics getMechanics
Function object to get string mechanics.
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.
JKatoomba< JAbstractMinimiser > evaluator
Recording of objects on file according a format that follows from the file name extension.
Auxiliary class for ROOT class selection.
JEvt getEvt(const JHead &header, const JModel &model)
Get event.
*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
#define MAKE_CSTRING(A)
Make C-string.
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Dynamic ROOT object management.
Auxiliary data structure for floating point format specification.
V(JDAQEvent-JTriggerReprocessor)*1.0/(JDAQEvent+1.0e-10)
floor_range getRangeOfFloors(const JDetector &detector)
Get range of floors.
Data structure for detector geometry and calibration.
Data structure for hydrophone.
size_t getNumberOfEmitters(T __begin, T __end)
Get number of emitters.
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.
Auxiliary class to manage set of compatible ROOT objects (e.g. histograms) using unique keys...
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
set_variable E_E log10(E_{fit}/E_{#mu})"
const std::string & getOID() const
Get detector identifier.
const array_type< JValue_t > & make_array(const JValue_t(&array)[N])
Method to create array of values.
JPosition3D getPosition(const Vec &pos)
Get position.
static struct JACOUSTICS::@4 compare
Auxiliary data structure to sort transmissions.
General purpose messaging.
static const JStringCounter getNumberOfStrings
Function object to count unique strings.
Implementation for depth dependend velocity of sound.
Scanning of objects from multiple files according a format that follows from the extension of each fi...
z range($ZMAX-$ZMIN)< $MINIMAL_DZ." fi fi typeset -Z 4 STRING typeset -Z 2 FLOOR JPlot1D -f $
void load(const std::string &file_name, JDetector &detector)
Load detector from input file.
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.
Acoustic transmission identifier.
bool hasLocation(const JLocation &location) const
Check if this detector has given location.
Custom probability density function of time-of-arrival.
const JLimit & getLimit() const
Get limit.
Fit functions of acoustic model.
do set_variable DETECTOR_TXT $WORKDIR detector
int getID() const
Get identifier.
Data structure for tripod.
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 JAcoustics sh $DETECTOR_ID source JAcousticsToolkit sh CHECK_EXIT_CODE typeset A EMITTERS get_tripods $WORKDIR tripod txt EMITTERS get_transmitters $WORKDIR transmitter txt EMITTERS for EMITTER in
Acoustic transmission identifier.
Template definition of fit function of acoustic model.
#define DEBUG(A)
Message macros.
JKatoomba< JGandalf > gandalf
Data structure for optical module.