60int main(
int argc,
char **argv)
73 JLimit_t& numberOfEvents = inputFile.getLimit();
79 disable_container disable;
88 JParser<> zap(
"Application to fit position calibration model to acoustic data.");
90 zap[
'f'] =
make_field(inputFile,
"output of JAcousticEventBuilder[.sh]");
102 "Precede name of data structure by a '+' or '-' "
103 "to add or remove data types in the output, respectively."
106 zap[
'q'] =
make_field(squash,
"squash meta data");
112 catch(
const exception &error) {
113 FATAL(error.what() << endl);
131 for (JDetector::const_iterator i =
detector.begin(); i !=
detector.end(); ++i) {
132 receivers[i->getID()] = i->getLocation();
135 for (tripods_container::const_iterator i = tripods.begin(); i != tripods.end(); ++i) {
136 emitters[i->getID()] =
JEmitter(i->getID(), i->getUTMPosition() -
detector.getUTMPosition());
139 for (transmitters_container::const_iterator i = transmitters.begin(); i != transmitters.end(); ++i) {
141 emitters[i->getID()] =
JEmitter(i->getID(), i->getPosition() +
detector.getModule(i->getLocation()).getPosition());
143 catch(
const exception&) {}
158 TH1D h0(
"chi2/NDF", NULL, 50000, 0.0, 1000.0);
159 TH1D h1(
"h1", NULL, 51, -0.5, 50.5);
160 TH1D hn(
"hn", NULL, 100, 0.0, 6.0);
170 for (Int_t i = 1; i <= HA->GetXaxis()->GetNbins(); ++i) {
171 HA->GetXaxis()->SetBinLabel(i,
MAKE_CSTRING(geometry.at(i-1).first));
172 HB->GetXaxis()->SetBinLabel(i,
MAKE_CSTRING(geometry.at(i-1).first));
175 if (inputFile.size() > 1u) {
179 for (
const string& file_name : inputFile) {
185 const JEvent* evt = in.next();
192 zmap[evt->begin()->getToE()] = file_name;
201 inputFile.push_back(i->second);
214 catch(
const exception&) {}
216 int counter[] = { 0, 0 };
230 if (!evt->empty() && emitters.
has(evt->
getID())) {
231 if (utc(evt->begin()->getToA()) && utc(evt->rbegin()->getToA())) {
232 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()) {
247 zbuf.erase(zbuf.begin(), p);
263 for (buffer_type::iterator evt = p; evt != q; ++evt) {
269 const JEmitter& emitter = emitters [evt->getID()];
270 const double weight =
getWeight(evt->getID());
272 for (JEvent::const_iterator i = evt->begin(); i != __end; ++i) {
277 if (receivers.
has(i->getID()) && geometry.
hasLocation(receivers[i->getID()]) && i->getQ() >= parameters.
Qmin * (parameters.
Qmin <= 1.0 ? i->getW() : 1.0)) {
279 data.push_back(
JHit(emitter,
281 receivers[i->getID()],
292 for (data_type::const_iterator hit = data.begin(); hit != data.end(); ++hit) {
293 HA[hit->getID()]->Fill(geometry.
getIndex(hit->getString()), hit->getFloor(), 1.0);
298 const auto result = katoomba(data.begin(), data.end());
300 for (data_type::const_iterator hit =
result.begin; hit !=
result.end; ++hit) {
301 HB[hit->getID()]->Fill(geometry.
getIndex(hit->getString()), hit->getFloor(), 1.0);
306 cout <<
"result:" <<
' '
310 for (data_type::const_iterator hit =
result.begin; hit !=
result.end; ++hit) {
337 for (buffer_type::iterator i = p; i != q; ++i) {
343 for (JEvent::iterator hit = out.begin(); hit != out.end(); ++hit) {
373 STATUS(
"Number of events written / rejected: " << counter[0] <<
" / " << counter[1] << endl);
ROOT TTree parameter settings.
Data structure for detector geometry and calibration.
Recording of objects on file according a format that follows from the file name extension.
Fit function of acoustic model.
General purpose class for hash map of unique elements.
Data structure for hydrophone.
Dynamic ROOT object management.
General purpose messaging.
#define DEBUG(A)
Message macros.
Data structure for optical module.
Scanning of objects from multiple files according a format that follows from the extension of each fi...
Utility class to parse command line options.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
#define MAKE_CSTRING(A)
Make C-string.
Auxiliary class to define a range between two values.
Scanning of objects from a single file according a format that follows from the extension of each fil...
Acoustic transmission identifier.
Data structure for transmitter.
Data structure for tripod.
std::vector< T >::difference_type distance(typename std::vector< T >::const_iterator first, typename PhysicsEvent::const_iterator< T > second)
Specialisation of STL distance.
Utility class to parse command line options.
Auxiliary class to manage set of compatible ROOT objects (e.g. histograms) using unique keys.
JTreeWriter object output.
JTreeWriter< T > & getTreeWriter()
Get TreeWriter.
General purpose class for object reading from a list of file names.
virtual bool hasNext() override
Check availability of next element.
virtual const pointer_type & next() override
Get next element.
const std::string & getFilename() const
Get current file name.
Object reading from a list of files.
virtual bool hasNext() override
Check availability of next element.
size_t getMinimumNumberOfEmitters(T __begin, T __end)
Get minimum number of emitters for any string in data.
double getWeight(T __begin, T __end)
Get total weight of data points.
JContainer< std::vector< JTripod > > tripods_container
JContainer< std::vector< JTransmitter > > transmitters_container
JContainer< std::vector< JHydrophone > > hydrophones_container
size_t getNumberOfEmitters(T __begin, T __end)
Get number of emitters.
static const JSoundVelocity getSoundVelocity(1541.0, -17.0e-3, -2000.0)
Function object for velocity of sound.
static JDetectorMechanics getMechanics
Function object to get string mechanics.
JEvt getEvt(const JHead &header, const JModel &model)
Get event.
floor_range getRangeOfFloors(const JDetector &detector)
Get range of floors.
void load(const std::string &file_name, JDetector &detector)
Load detector from input file.
static const JStringCounter getNumberOfStrings
Function object to count unique strings.
JComparator< JResult_t T::*, JComparison::lt > make_comparator(JResult_t T::*member)
Helper method to create comparator between values of data member.
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
std::vector< JHitW0 > buffer_type
hits
std::vector< event_type > data_type
int main(int argc, char **argv)
Auxiliary data structure for floating point format specification.
int getID() const
Get emitter identifier.
const int getDetectorID() const
Get detector identifier.
double Qmin
minimal quality transmission
double deadTime_s
dead time between events [s]
size_t Nmin
minimum number of emitters
double sigma_s
time-of-arrival resolution [s]
double Tmax_s
time window to combine events [s]
double chi2perNDF
maximal chi2/NDF to store event
bool hasLocation(const JLocation &location) const
Check if this detector has given location.
Global fit of prameterised detector geometry to acoustics data.
JKatoomba< JGandalf > gandalf
JKatoomba< JAbstractMinimiser > evaluator
Template definition of fit function of acoustic model.
Implementation for depth dependend velocity of sound.
JSoundVelocity & set(const double z0)
Set depth.
Acoustic transmission identifier.
int getID() const
Get identifier.
Auxiliary data structure to unify weights of acoustics data according to the number of pings per emit...
Data structure for measured coincidence rates of all pairs of PMTs in optical module.
Auxiliary wrapper for I/O of container with optional comment (see JComment).
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Auxiliary class for ROOT class selection.
bool is_valid() const
Get status of given data type.
Auxiliary class for defining the range of iterations of objects.
static counter_type max()
Get maximum counter value.