Jpp  18.6.0-rc.1
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JPlatypus_t.hh
Go to the documentation of this file.
1 #ifndef __JACOUSTICS__JPLATYPUS_T__
2 #define __JACOUSTICS__JPLATYPUS_T__
3 
4 #include <string>
5 #include <type_traits>
6 #include <functional>
7 #include <future>
8 #include <mutex>
9 #include <thread>
10 #include <vector>
11 
12 #include "JLang/JObjectIterator.hh"
13 
14 #include "JMath/JQuantile_t.hh"
15 
16 #include "JTools/JHashMap.hh"
17 
18 #include "JAcoustics/JHit.hh"
19 #include "JAcoustics/JEmitter.hh"
22 #include "JAcoustics/JSuperEvt.hh"
24 
25 
26 namespace JACOUSTICS {}
27 namespace JPP { using namespace JACOUSTICS; }
28 
29 namespace JACOUSTICS {
30 
31  /**
32  * Thread pool for global fits using super events.
33  */
34  class JPlatypus {
35  public:
36 
39 
40  /**
41  * Constructor.
42  *
43  * \param geometry detector geometry
44  * \param emitters emitters geometry
45  * \param velocity sound velocity
46  * \param parameters parameters
47  * \param input input
48  * \param ns number of threads
49  */
50  JPlatypus(const JGeometry& geometry,
51  const JEmitters& emitters,
52  const JSoundVelocity& velocity,
55  const size_t ns) :
56  input(input),
57  stop(false)
58  {
59  using namespace std;
60  using namespace JPP;
61 
62  Q.reset();
63 
64  for (size_t i = 0; i < ns; ++i) {
65 
66  thread worker([this, geometry, emitters, velocity, parameters]() {
67 
69 
70  for (JKatoomba<JGandalf> katoomba(geometry, velocity, parameters.option); ; ) {
71 
72  katoomba.estimator.reset(getMEstimator(parameters.mestimator));
73 
74  {
75  unique_lock<mutex> lock(in);
76 
77  cv.wait(lock, [this]() { return stop || this->input.hasNext(); });
78 
79  if (stop && !this->input.hasNext()) {
80  return;
81  }
82 
83  const JSuperEvt* evt = this->input.next();
84 
85  data.clear();
86 
87  for (JSuperEvt::rx_type::const_iterator hit = evt->rx.begin(); hit != evt->rx.end(); ++hit) {
88 
89  data.push_back(JHit(emitters[hit->id],
90  hit->counter,
91  JLocation(hit->string, hit->floor),
92  hit->toa,
93  parameters.sigma_s,
94  hit->weight));
95  }
96 
97  katoomba.value = getModel(*evt);
98  }
99 
100  const double chi2 = katoomba (data.begin(), data.end()) / katoomba.estimator->getRho(1.0);
101  const double ndf = getWeight(data.begin(), data.end()) - katoomba.value.getN();
102 
103  {
104  unique_lock<mutex> lock(out);
105 
106  Q.put(chi2 / ndf);
107  }
108  }
109  });
110 
111  workers.emplace_back(std::move(worker));
112  }
113  }
114 
115 
116  /**
117  * Destructor.
118  */
120  {
121  using namespace std;
122 
123  {
124  unique_lock<mutex> lock(in);
125 
126  stop = true;
127  }
128 
129  cv.notify_all();
130 
131  for (auto& worker : workers) {
132  worker.join();
133  }
134  }
135 
136  static JMATH::JQuantile_t Q; //!< chi2/NDF
137 
138  private:
141  std::mutex in;
142  std::mutex out;
143  std::condition_variable cv;
144  bool stop;
145  };
146 }
147 
148 #endif
JModel getModel(const JEvt &evt)
Get model.
Acoustic hit.
double getWeight(T __begin, T __end)
Get total weight of data points.
Definition: JKatoomba_t.hh:61
General purpose class for hash map of unique keys.
Definition: JHashMap.hh:72
General purpose class for hash map of unique elements.
static JMATH::JQuantile_t Q
chi2/NDF
Definition: JPlatypus_t.hh:136
std::vector< size_t > ns
void put(const double x)
Put value.
Definition: JQuantile_t.hh:46
*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
Interface of object iteration for a single data type.
Acoustic emitter.
JPlatypus(const JGeometry &geometry, const JEmitters &emitters, const JSoundVelocity &velocity, const JFitParameters &parameters, input_type &input, const size_t ns)
Constructor.
Definition: JPlatypus_t.hh:50
Acoustics hit.
Acoustic super event fit toolkit.
Acoustic fit parameters.
Template specialisation of fit function of acoustic model based on JGandalf minimiser.
Definition: JKatoomba_t.hh:678
virtual const pointer_type & next()=0
Get next element.
Logical location of module.
Definition: JLocation.hh:37
Acoustic event fit.
JTOOLS::JHashMap< int, JEmitter > JEmitters
Definition: JPlatypus_t.hh:37
virtual bool hasNext()=0
Check availability of next element.
std::vector< std::thread > workers
Definition: JPlatypus_t.hh:139
Implementation for depth dependend velocity of sound.
Acoustic super event fit.
Definition: JSuperEvt.hh:36
~JPlatypus()
Destructor.
Definition: JPlatypus_t.hh:119
Thread pool for global fits using super events.
Definition: JPlatypus_t.hh:34
Fit functions of acoustic model.
then if[[!-f $DETECTOR]] then JDetector sh $DETECTOR fi cat $WORKDIR trigger_parameters txt<< EOFtrigger3DMuon.enabled=1;trigger3DMuon.numberOfHits=5;trigger3DMuon.gridAngle_deg=1;ctMin=0.0;TMaxLocal_ns=15.0;EOF set_variable TRIGGEREFFICIENCY_TRIGGERED_EVENTS_ONLY INPUT_FILES=() for((i=1;$i<=$NUMBER_OF_RUNS;++i));do JSirene.sh $DETECTOR $JPP_DATA/genhen.km3net_wpd_V2_0.evt.gz $WORKDIR/sirene_ ${i}.root JTriggerEfficiency.sh $DETECTOR $DETECTOR $WORKDIR/sirene_ ${i}.root $WORKDIR/trigger_efficiency_ ${i}.root $WORKDIR/trigger_parameters.txt $JPP_DATA/PMT_parameters.txt INPUT_FILES+=($WORKDIR/trigger_efficiency_ ${i}.root) done for ANGLE_DEG in $ANGLES_DEG[*];do set_variable SIGMA_NS 3.0 set_variable OUTLIERS 3 set_variable OUTPUT_FILE $WORKDIR/matrix\[${ANGLE_DEG}\deg\].root $JPP_DIR/examples/JReconstruction-f"$INPUT_FILES[*]"-o $OUTPUT_FILE-S ${SIGMA_NS}-A ${ANGLE_DEG}-O ${OUTLIERS}-d ${DEBUG}--!fiif[[$OPTION=="plot"]];then if((0));then for H1 in h0 h1;do JPlot1D-f"$WORKDIR/matrix["${^ANGLES_DEG}" deg].root:${H1}"-y"1 2e3"-Y-L TR-T""-\^"number of events [a.u.]"-> o chi2
Definition: JMatrixNZ.sh:106
double sigma_s
time-of-arrival resolution [s]
int mestimator
M-estimator.
JMEstimator * getMEstimator(const int type)
Get M-Estimator.
Definition: JMEstimator.hh:203
Auxiliary data structure for average.
Definition: JQuantile_t.hh:21
JLANG::JObjectIterator< JSuperEvt > input_type
Definition: JPlatypus_t.hh:38
void reset()
Reset.
Definition: JQuantile_t.hh:34
std::condition_variable cv
Definition: JPlatypus_t.hh:143