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
JShower3EZRegressor.hh
Go to the documentation of this file.
1 #ifndef __JFIT__JSHOWER3EZREGRESSOR__
2 #define __JFIT__JSHOWER3EZREGRESSOR__
3 
4 #include "JPhysics/JPDFTypes.hh"
5 #include "JPhysics/JPDFTable.hh"
6 #include "JPhysics/JPDFToolkit.hh"
7 #include "JPhysics/JNPETable.hh"
8 
9 #include "JTools/JFunction1D_t.hh"
10 #include "JTools/JFunctionalMap_t.hh"
11 #include "JPhysics/JConstants.hh"
12 #include "JTools/JRange.hh"
13 #include "JTools/JResult.hh"
14 #include "JTools/JFunction1D_t.hh"
15 #include "JTools/JFunctionalMap_t.hh"
16 #include "JTools/JAbstractHistogram.hh"
17 
18 #include "JGeometry3D/JVector3D.hh"
19 #include "JGeometry3D/JVersor3D.hh"
20 #include "JGeometry3D/JDirection3D.hh"
21 
22 #include "JMath/JZero.hh"
23 
24 #include "JFit/JTimeRange.hh"
25 #include "JFit/JPMTW0.hh"
26 #include "JFit/JSimplex.hh"
27 #include "JFit/JGandalf.hh"
28 #include "JFit/JMEstimator.hh"
29 #include "JFit/JRegressor.hh"
30 #include "JFit/JShower3EZ.hh"
31 #include "JFit/JFitToolkit.hh"
32 #include "JFit/JLine3Z.hh"
33 
34 #include "Jeep/JMessage.hh"
35 
36 /**
37  * \file
38  * Data regression method for JFIT::JShower3EZ.
39  * \author mdejong, vcarretero
40  */
41 
42 namespace JFIT {
43 
44  using JPHYSICS::JPDFType_t;
45  using JPHYSICS::DIRECT_LIGHT_FROM_EMSHOWER;
46  using JPHYSICS::SCATTERED_LIGHT_FROM_EMSHOWER;
47  using JTOOLS::get_value;
48 
49  /**
50  * Constrain PMT angle to [0,pi].
51  *
52  * \param angle angle [rad]
53  * \return angle [rad]
54  */
55  inline double getPMTAngle(const double angle)
56  {
57  const double epsilon = 1.0e-6;
58  const JTOOLS::JRange<double> range(epsilon, JMATH::PI - epsilon);
59 
60  return range.constrain(fabs(angle));
61  }
62 
63 
64  /**
65  * Function to constrain the versor and energy during the fit, to prevent unphysical values.
66  *
67  * \param value model (I/O)
68  */
69  void model(JShower3EZ& value)
70  {
71  using namespace std;
72 
73 
74  double Tx = value.getDX();
75  double Ty = value.getDY();
76  double E = max(0.0,value.getE());
77  const double u = hypot(Tx, Ty);
78 
79  if (u > 1.0) {
80  Tx /= u;
81  Ty /= u;
82  }
83 
84  value = JShower3EZ(static_cast<const JPoint4D&>(value), JVersor3Z(Tx,Ty), E, value.getBy());
85 
86  }
87 
88  /**
89  * Regressor function object for JShower3EZ fit using JSimplex minimiser.
90  */
91  template<>
92  struct JRegressor<JShower3EZ, JSimplex> :
93  public JAbstractRegressor<JShower3EZ, JSimplex>
94  {
95  using JAbstractRegressor<JShower3EZ, JSimplex>::operator();
96 
97  typedef JTOOLS::JSplineFunction1S_t JFunction1D_t;
98  typedef JTOOLS::JMapList<JTOOLS::JPolint0FunctionalMap,
99  JTOOLS::JMapList<JTOOLS::JPolint0FunctionalMap,
100  JTOOLS::JMapList<JTOOLS::JPolint0FunctionalGridMap,
101  JTOOLS::JMapList<JTOOLS::JPolint0FunctionalGridMap> > > > JPDFMaplist_t;
102  typedef JPHYSICS::JPDFTable<JFunction1D_t, JPDFMaplist_t> JPDF_t;
103 
104  typedef JTOOLS::JMapList<JTOOLS::JPolint0FunctionalMap,
105  JTOOLS::JMapList<JTOOLS::JPolint0FunctionalMap,
106  JTOOLS::JMapList<JTOOLS::JPolint0FunctionalGridMap,
107  JTOOLS::JMapList<JTOOLS::JPolint0FunctionalGridMap> > > > JNPEMaplist_t;
108  typedef JPHYSICS::JNPETable<double, double, JNPEMaplist_t> JNPE_t;
109 
110  /**
111  * Parameterized constructor
112  *
113  * The PDF file descriptor should contain the wild card character JPHYSICS::WILDCARD which
114  * will be replaced by the corresponding PDF types listed in JRegressor<JShower3Z, JGandalf>::pdf_t.
115  *
116  * \param fileDescriptor PDF file descriptor
117  */
118 
119  JRegressor(const std::string& fileDescriptor):
120  estimator(new JMEstimatorNull())
121  {
122  using namespace std;
123  using namespace JPP;
124 
125  const JPDF_t::JSupervisor supervisor(new JPDF_t::JDefaultResult(JMATH::zero));
126 
127  for (int i = 0; i != NUMBER_OF_PDFS; ++i) {
128 
129  try {
130 
131  JPDF_t pdf;
132 
133  const string file_name = getFilename(fileDescriptor, pdf_t[i]);
134 
135  NOTICE("loading PDF from file " << file_name << "... " << flush);
136 
137  pdf.load(file_name.c_str());
138 
139  NOTICE("OK" << endl);
140 
141  pdf.setExceptionHandler(supervisor);
142 
143  npe[ i ] = JNPE_t(pdf);
144  }
145  catch(const JException& error) {
146  FATAL(error.what() << endl);
147  }
148  }
149 
150  // Add PDFs
151  for (int i = 1; i < NUMBER_OF_PDFS; i += 2) {
152 
153  npe[ i ].add(npe[i-1]);
154 
155  JNPE_t buffer;
156 
157  npe[i-1].swap(buffer);
158  }
159  }
160 
161  /**
162  * Fit function.
163  * This method is used to determine the chi2 of given PMT with respect to shower hypothesis.
164  *
165  * \param shower shower
166  * \param pmt pmt
167  * \return chi2
168  */
169  double operator()(const JShower3EZ& shower, const JPMTW0& pmt) const
170  {
171  using namespace JPP;
172 
173  JPosition3D D(pmt.getPosition());
174  JDirection3D U(pmt.getDirection());
175 
176  D.sub(shower.getPosition());
177 
178  const double z = D.getDot(shower.getDirection());
179  const double x = D.getX();
180  const double y = D.getY();
181  const double cd = z/D.getLength(); // cosine angle between shower direction and PMT position
182 
183  U.rotate(JRotation3Z(-atan2(y,x))); // rotate PMT axis to x-z plane
184 
185  const double theta = getPMTAngle(U.getTheta());
186  const double phi = getPMTAngle(U.getPhi());
187 
188  JNPE_t::result_type H0 = getH0(pmt.getR()); // background hypothesis value for time integrated PDF.
189  JNPE_t::result_type H1 = getH1(D.getLength(), cd, theta, phi, shower.getE()); // signal hypothesis value for time integrated PDF.
190 
191  if (get_value(H1) >= Vmax_npe) {
192  H1 *= Vmax_npe / get_value(H1);
193  }
194 
195  H1 += H0; // now H1 is signal + background
196 
197  const bool hit = pmt.getN() != 0;
198  const double u = getChi2(get_value(H1), hit); // -log(lik)
199 
200  return estimator->getRho(u);
201  }
202 
203  /**
204  * Get background hypothesis value for time integrated PDF.
205  *
206  * \param R_Hz rate [Hz]
207  * \return hypothesis value
208  */
209  JNPE_t::result_type getH0(const double R_Hz) const
210  {
211  return JNPE_t::result_type(R_Hz * 1e-9 * T_ns.getLength());
212  }
213 
214  /**
215  * Get signal hypothesis value for time integrated PDF.
216  *
217  * \param D PMT distance from shower [m]
218  * \param cd cosine angle between shower direction and PMT position
219  * \param theta PMT zenith angle [deg]
220  * \param phi PMT azimuth angle [deg]
221  * \param E shower energy [GeV]
222  * \return hypothesis value
223  */
224  JNPE_t::result_type getH1(const double D,
225  const double cd,
226  const double theta,
227  const double phi,
228  const double E) const
229  {
230  JNPE_t::result_type h1 = JMATH::zero;
231 
232  for (int i = 0; i != NUMBER_OF_PDFS; ++i) {
233 
234  if (!npe[i].empty() && D <= npe[i].getXmax()) {
235 
236  try {
237 
238  JNPE_t::result_type y1 = E * npe[i](std::max(D, npe[i].getXmin()), cd, theta, phi);
239 
240  // safety measure
241 
242  if(y1 < 0){
243  y1 = 0;
244  }
245 
246  h1 += y1;
247 
248  }
249  catch(JLANG::JException& error) {
250  ERROR(error << std::endl);
251  }
252  }
253  }
254 
255  return h1;
256  }
257 
258 
259  static JTimeRange T_ns; //!< Time window with respect to Cherenkov hypothesis [ns]
260  static double Vmax_npe; //!< Maximal integral of PDF [npe]
261 
262  static const int NUMBER_OF_PDFS = 2;
263 
264  static const JPDFType_t pdf_t[NUMBER_OF_PDFS];
265 
266  JNPE_t npe[NUMBER_OF_PDFS]; //!< PDF
267 
268  JLANG::JSharedPointer<JMEstimator> estimator; //!< M-Estimator function
269  };
270 
271 
272  /**
273  * Regressor function object for JShower3EZ fit using JGandalf minimiser.
274  */
275  template<>
276  struct JRegressor<JShower3EZ, JGandalf> :
277  public JAbstractRegressor<JShower3EZ, JGandalf>
278  {
279  using JAbstractRegressor<JShower3EZ, JGandalf>::operator();
280 
281  typedef JTOOLS::JSplineFunction1S_t JFunction1D_t;
282  typedef JTOOLS::JMapList<JTOOLS::JPolint0FunctionalMap,
283  JTOOLS::JMapList<JTOOLS::JPolint0FunctionalMap,
284  JTOOLS::JMapList<JTOOLS::JPolint0FunctionalGridMap,
285  JTOOLS::JMapList<JTOOLS::JPolint0FunctionalGridMap> > > > JPDFMaplist_t;
286  typedef JPHYSICS::JPDFTable<JFunction1D_t, JPDFMaplist_t> JPDF_t;
287 
288  typedef JTOOLS::JMapList<JTOOLS::JPolint1FunctionalMap,
289  JTOOLS::JMapList<JTOOLS::JPolint1FunctionalMapH,
290  JTOOLS::JMapList<JTOOLS::JPolint1FunctionalGridMap,
291  JTOOLS::JMapList<JTOOLS::JPolint1FunctionalGridMap> > > > JNPEMaplist_t;
292  typedef JPHYSICS::JNPETable<double, double, JNPEMaplist_t> JNPE_t;
293 
294  /**
295  * Parameterized constructor
296  *
297  * The PDF file descriptor should contain the wild card character JPHYSICS::WILDCARD which
298  * will be replaced by the corresponding PDF types listed in JRegressor<JShower3Z, JGandalf>::pdf_t.
299  *
300  * \param fileDescriptor PDF file descriptor
301  */
302 
303  JRegressor(const std::string& fileDescriptor):
304  estimator(new JMEstimatorNull())
305  {
306  using namespace std;
307  using namespace JPP;
308 
309  const JPDF_t::JSupervisor supervisor(new JPDF_t::JDefaultResult(JMATH::zero));
310 
311  for (int i = 0; i != NUMBER_OF_PDFS; ++i) {
312 
313  try {
314 
315  JPDF_t pdf;
316 
317  const string file_name = getFilename(fileDescriptor, pdf_t[i]);
318 
319  NOTICE("loading PDF from file " << file_name << "... " << flush);
320 
321  pdf.load(file_name.c_str());
322 
323  NOTICE("OK" << endl);
324 
325  pdf.setExceptionHandler(supervisor);
326 
327  npe[i] = JNPE_t(pdf);
328  }
329  catch(const JException& error) {
330  FATAL(error.what() << endl);
331  }
332  }
333 
334  // Add PDFs
335  for (int i = 1; i < NUMBER_OF_PDFS; i += 2) {
336 
337  npe[ i ].add(npe[i-1]);
338 
339  JNPE_t buffer;
340 
341  npe[i-1].swap(buffer);
342  }
343  }
344 
345  /**
346  * Fit function.
347  * This method is used to determine the chi2 of given PMT with respect to shower hypothesis.
348  *
349  * \param shower shower
350  * \param pmt pmt
351  * \return chi2
352  */
353  result_type operator()(const JShower3EZ& shower, const JPMTW0& pmt) const
354  {
355  using namespace JPP;
356  using namespace std;
357 
358  JPosition3D D(pmt.getPosition());
359  JDirection3D U(pmt.getDirection());
360 
361  D.sub(shower.getPosition());
362 
363  const double x = D.getX();
364  const double y = D.getY();
365  const double d = D.getLength();
366  const double cd = D.getDot(shower.getDirection())/d; // cosine angle between shower direction and PMT position
367 
368  U.rotate(JRotation3Z(-atan2(y,x))); // rotate PMT axis to x-z plane
369 
370  const double theta = getPMTAngle(U.getTheta());
371  const double phi = getPMTAngle(U.getPhi());
372 
373  JNPE_t::result_type H0 = getH0(pmt.getR()); // background hypothesis value for time integrated PDF.
374  JNPE_t::result_type H1 = getH1(d, cd, theta, phi, shower.getE()); // signal hypothesis value for time integrated PDF.
375 
376  if (get_value(H1) >= Vmax_npe) {
377  H1 *= Vmax_npe / get_value(H1);
378  }
379 
380  double signal_npe = get_value(H1);
381 
382  H1 += H0; // now H1 is signal + background
383 
384  double expectation = get_value(H1);
385 
386  const bool hit = pmt.getN() != 0;
387  const double u = H1.getChi2(hit);
388 
389  result_type result;
390 
391  result.chi2 = estimator->getRho(u);
392 
393  double energy_gradient = signal_npe/shower.getE(); // dP/dE
394  if(hit) energy_gradient *= -exp(-expectation)/(1-exp(-expectation)); //dchi2/d(H1), if !hit is 1
395 
396  result.gradient = JShower3EZ(JPoint4D(JVector3D(0, // d(cos_th0)/d(x)
397  0, // d(cos_th0)/d(y)
398  0), // d(cos_th0)/d(z)
399  0.0), // d(cos_th0)/d(t)
400  JVersor3Z(x/d, // d(cos_th0)/d(dx)
401  y/d), // d(cos_th0)/d(dy)
402  energy_gradient); // d(chi2)/d(E)
403 
404  result.gradient.mul(estimator->getPsi(u));
405  static_cast<JShower3Z&>(result.gradient).mul(H1.getDerivativeOfChi2(hit)); // x d(chi2)/d(cos_th0)
406 
407  return result;
408  }
409 
410  /**
411  * Get background hypothesis value for time integrated PDF.
412  *
413  * \param R_Hz rate [Hz]
414  * \return hypothesis value
415  */
416  JNPE_t::result_type getH0(const double R_Hz) const
417  {
418  return JNPE_t::result_type(R_Hz * 1e-9 * T_ns.getLength(), 0.0);
419  }
420 
421  /**
422  * Get signal hypothesis value for time integrated PDF.
423  *
424  * \param D PMT distance from shower [m]
425  * \param cd cosine angle between shower direction and PMT position
426  * \param theta PMT zenith angle [deg]
427  * \param phi PMT azimuth angle [deg]
428  * \param E shower energy [GeV]
429  * \return hypothesis value
430  */
431  JNPE_t::result_type getH1(const double D,
432  const double cd,
433  const double theta,
434  const double phi,
435  const double E) const
436  {
437  JNPE_t::result_type h1 = JMATH::zero;
438 
439  for (int i = 0; i != NUMBER_OF_PDFS; ++i) {
440 
441  if (!npe[i].empty() && D <= npe[i].getXmax()) {
442 
443  try {
444 
445  JNPE_t::result_type y1 = E * npe[i](std::max(D, npe[i].getXmin()), cd, theta, phi);
446 
447  if (get_value(y1) > 0.0) {
448  h1 += y1;
449  }
450 
451  }
452  catch(JLANG::JException& error) {
453  ERROR(error << std::endl);
454  }
455  }
456  }
457 
458  return h1;
459  }
460 
461 
462  static JTimeRange T_ns; //!< Time window with respect to Cherenkov hypothesis [ns]
463  static double Vmax_npe; //!< Maximal integral of PDF [npe]
464 
465  static const int NUMBER_OF_PDFS = 2;
466 
467  static const JFIT::JPDFType_t pdf_t[NUMBER_OF_PDFS];
468 
469  JNPE_t npe[NUMBER_OF_PDFS]; //!< PDF
470 
471  JLANG::JSharedPointer<JMEstimator> estimator; //!< M-Estimator function
472  };
473 
474 /**
475  * Regressor function object for JShower3EZ fit using Abstract minimiser, that just computes the chi2 without a fit.
476  */
477  template<>
478  struct JRegressor<JShower3EZ, JAbstractMinimiser> :
479  public JAbstractRegressor<JShower3EZ, JAbstractMinimiser>
480  {
481  using JAbstractRegressor<JShower3EZ, JAbstractMinimiser>::operator();
482 
483  typedef JTOOLS::JSplineFunction1S_t JFunction1D_t;
484  typedef JTOOLS::JMapList<JTOOLS::JPolint0FunctionalMap,
485  JTOOLS::JMapList<JTOOLS::JPolint0FunctionalMap,
486  JTOOLS::JMapList<JTOOLS::JPolint0FunctionalGridMap,
487  JTOOLS::JMapList<JTOOLS::JPolint0FunctionalGridMap> > > > JPDFMaplist_t;
488  typedef JPHYSICS::JPDFTable<JFunction1D_t, JPDFMaplist_t> JPDF_t;
489 
490  typedef JTOOLS::JMapList<JTOOLS::JPolint0FunctionalMap,
491  JTOOLS::JMapList<JTOOLS::JPolint0FunctionalMap,
492  JTOOLS::JMapList<JTOOLS::JPolint0FunctionalGridMap,
493  JTOOLS::JMapList<JTOOLS::JPolint0FunctionalGridMap> > > > JNPEMaplist_t;
494  typedef JPHYSICS::JNPETable<double, double, JNPEMaplist_t> JNPE_t;
495 
496  /**
497  * Parameterized constructor
498  *
499  * The PDF file descriptor should contain the wild card character JPHYSICS::WILDCARD which
500  * will be replaced by the corresponding PDF types listed in JRegressor<JShower3Z, JGandalf>::pdf_t.
501  *
502  * \param fileDescriptor PDF file descriptor
503  */
504 
505  JRegressor(const std::string& fileDescriptor):
506  estimator(new JMEstimatorNull())
507  {
508  using namespace std;
509  using namespace JPP;
510 
511  const JPDF_t::JSupervisor supervisor(new JPDF_t::JDefaultResult(JMATH::zero));
512 
513  for (int i = 0; i != NUMBER_OF_PDFS; ++i) {
514 
515  try {
516 
517  JPDF_t pdf;
518 
519  const string file_name = getFilename(fileDescriptor, pdf_t[i]);
520 
521  NOTICE("loading PDF from file " << file_name << "... " << flush);
522 
523  pdf.load(file_name.c_str());
524 
525  NOTICE("OK" << endl);
526 
527  pdf.setExceptionHandler(supervisor);
528 
529  npe[ i ] = JNPE_t(pdf);
530  }
531  catch(const JException& error) {
532  FATAL(error.what() << endl);
533  }
534  }
535 
536  // Add PDFs
537  for (int i = 1; i < NUMBER_OF_PDFS; i += 2) {
538 
539  npe[ i ].add(npe[i-1]);
540 
541  JNPE_t buffer;
542 
543  npe[i-1].swap(buffer);
544  }
545  }
546 
547  /**
548  * Fit function.
549  * This method is used to determine the chi2 of given PMT with respect to shower hypothesis.
550  *
551  * \param shower shower
552  * \param pmt pmt
553  * \return chi2
554  */
555  double operator()(const JShower3EZ& shower, const JPMTW0& pmt) const
556  {
557  using namespace JPP;
558 
559  JPosition3D D(pmt.getPosition());
560  JDirection3D U(pmt.getDirection());
561 
562  D.sub(shower.getPosition());
563 
564  const double z = D.getDot(shower.getDirection());
565  const double x = D.getX();
566  const double y = D.getY();
567  const double cd = z/D.getLength(); // cosine angle between shower direction and PMT position
568 
569  U.rotate(JRotation3Z(-atan2(y,x))); // rotate PMT axis to x-z plane
570 
571  const double theta = getPMTAngle(U.getTheta());
572  const double phi = getPMTAngle(U.getPhi());
573 
574  JNPE_t::result_type H0 = getH0(pmt.getR()); // background hypothesis value for time integrated PDF.
575  JNPE_t::result_type H1 = getH1(D.getLength(), cd, theta, phi, shower.getE()); // signal hypothesis value for time integrated PDF.
576 
577  if (get_value(H1) >= Vmax_npe) {
578  H1 *= Vmax_npe / get_value(H1);
579  }
580 
581  H1 += H0; // now H1 is signal + background
582 
583  const bool hit = pmt.getN() != 0;
584  const double u = getChi2(get_value(H1), hit); // -log(lik)
585 
586  return estimator->getRho(u);
587  }
588 
589  /**
590  * Get background hypothesis value for time integrated PDF.
591  *
592  * \param R_Hz rate [Hz]
593  * \return hypothesis value
594  */
595  JNPE_t::result_type getH0(const double R_Hz) const
596  {
597  return JNPE_t::result_type(R_Hz * 1e-9 * T_ns.getLength());
598  }
599 
600  /**
601  * Get signal hypothesis value for time integrated PDF.
602  *
603  * \param D PMT distance from shower [m]
604  * \param cd cosine angle between shower direction and PMT position
605  * \param theta PMT zenith angle [deg]
606  * \param phi PMT azimuth angle [deg]
607  * \param E shower energy [GeV]
608  * \return hypothesis value
609  */
610  JNPE_t::result_type getH1(const double D,
611  const double cd,
612  const double theta,
613  const double phi,
614  const double E) const
615  {
616  JNPE_t::result_type h1 = JMATH::zero;
617 
618  for (int i = 0; i != NUMBER_OF_PDFS; ++i) {
619 
620  if (!npe[i].empty() && D <= npe[i].getXmax()) {
621 
622  try {
623 
624  JNPE_t::result_type y1 = E * npe[i](std::max(D, npe[i].getXmin()), cd, theta, phi);
625 
626  // safety measure
627 
628  if(y1 < 0){
629  y1 = 0;
630  }
631 
632  h1 += y1;
633 
634  }
635  catch(JLANG::JException& error) {
636  ERROR(error << std::endl);
637  }
638  }
639  }
640 
641  return h1;
642  }
643 
644 /**
645  * Get signal hypothesis value for time integrated PDF.
646  *
647  * \param shower shower
648  * \param pmt pmt
649  * \return hypothesis value
650  */
651  JNPE_t::result_type getH1(const JShower3EZ& shower, const JPMTW0& pmt) const
652  {
653  using namespace JPP;
654 
655  JPosition3D D(pmt.getPosition());
656  JDirection3D U(pmt.getDirection());
657 
658 
659  const double z = D.getDot(shower.getDirection());
660  const double x = D.getX();
661  const double y = D.getY();
662  const double cd = z/D.getLength(); // cosine angle between shower direction and PMT position
663 
664  U.rotate(JRotation3Z(-atan2(y,x))); // rotate PMT axis to x-z plane
665 
666  const double theta = getPMTAngle(U.getTheta());
667  const double phi = getPMTAngle(U.getPhi());
668 
669  JNPE_t::result_type H1 = getH1(D.getLength(), cd, theta, phi, 1.0); // signal hypothesis value for time integrated PDF. E=1 because it is linear with E.
670 
671  if (get_value(H1) >= Vmax_npe) {
672  H1 *= Vmax_npe / get_value(H1);
673  }
674 
675  return H1;
676  }
677  static JTimeRange T_ns; //!< Time window with respect to Cherenkov hypothesis [ns]
678  static double Vmax_npe; //!< Maximal integral of PDF [npe]
679 
680  static const int NUMBER_OF_PDFS = 2;
681 
682  static const JPDFType_t pdf_t[NUMBER_OF_PDFS];
683 
684  JNPE_t npe[NUMBER_OF_PDFS]; //!< PDF
685 
686  JLANG::JSharedPointer<JMEstimator> estimator; //!< M-Estimator function
687 
688  };
689 
690  /**
691  * PDF types.
692  */
693  const JPDFType_t JRegressor<JShower3EZ, JSimplex>::pdf_t[] = { DIRECT_LIGHT_FROM_EMSHOWER,
694  SCATTERED_LIGHT_FROM_EMSHOWER };
695 
696  const JPDFType_t JRegressor<JShower3EZ, JGandalf>::pdf_t[] = { DIRECT_LIGHT_FROM_EMSHOWER,
697  SCATTERED_LIGHT_FROM_EMSHOWER };
698  const JPDFType_t JRegressor<JShower3EZ, JAbstractMinimiser>::pdf_t[] = { DIRECT_LIGHT_FROM_EMSHOWER,
699  SCATTERED_LIGHT_FROM_EMSHOWER };
700  /**
701  * Default values.
702  */
703  JTimeRange JRegressor<JShower3EZ, JSimplex>::T_ns;
704  double JRegressor<JShower3EZ, JSimplex>::Vmax_npe = std::numeric_limits<double>::max();
705 
706  JTimeRange JRegressor<JShower3EZ, JGandalf>::T_ns;
707  double JRegressor<JShower3EZ, JGandalf>::Vmax_npe = std::numeric_limits<double>::max();
708 
709  JTimeRange JRegressor<JShower3EZ, JAbstractMinimiser>::T_ns;
710  double JRegressor<JShower3EZ, JAbstractMinimiser>::Vmax_npe = std::numeric_limits<double>::max();
711 
712 
713 
714 }
715 
716 #endif
JFitToolkit.hh
Auxiliary methods to evaluate Poisson probabilities and chi2.
JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::getH1
JNPE_t::result_type getH1(const double D, const double cd, const double theta, const double phi, const double E) const
Get signal hypothesis value for time integrated PDF.
Definition: JShower3EZRegressor.hh:610
JFIT::JRegressor
Template definition of a data regressor of given model.
Definition: JRegressor.hh:70
JFIT::JRegressor< JShower3EZ, JGandalf >::JPDF_t
JPHYSICS::JPDFTable< JFunction1D_t, JPDFMaplist_t > JPDF_t
Definition: JShower3EZRegressor.hh:286
JTOOLS::JSplineFunction1S_t
Type definition of a spline interpolation method based on a JCollection with JResultPDF result type...
Definition: JFunction1D_t.hh:51
JLANG::JException
General exception.
Definition: JException.hh:24
JIO::JObjectBinaryIO::load
void load(const char *file_name)
Load from input file.
Definition: JObjectBinaryIO.hh:29
JRegressor.hh
General purpose data regression method.
JFIT::JPMTW0::getR
double getR() const
Get rate.
Definition: JPMTW0.hh:56
JFIT::JMEstimatorNull
Null M-estimator.
Definition: JMEstimator.hh:92
JPDFToolkit.hh
Auxiliary methods for PDF calculations.
E
then usage $script< input file >[option[primary[working directory]]] nWhere option can be E
Definition: JMuonPostfit.sh:40
JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::JRegressor
JRegressor(const std::string &fileDescriptor)
Parameterized constructor.
Definition: JShower3EZRegressor.hh:505
JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::JPDFMaplist_t
JTOOLS::JMapList< JTOOLS::JPolint0FunctionalMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalGridMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalGridMap > > > > JPDFMaplist_t
Definition: JShower3EZRegressor.hh:487
JResult.hh
This include file containes various data structures that can be used as specific return types for the...
JFIT::JPMTW0
Auxiliary class for handling PMT geometry, rate and response.
Definition: JPMTW0.hh:22
JFIT::JPoint4D
Data structure for vertex fit.
Definition: JPoint4D.hh:22
JPHYSICS::SCATTERED_LIGHT_FROM_EMSHOWER
scattered light from EM shower
Definition: JPDFTypes.hh:38
JGEOMETRY3D::JDirection3D::getDirection
const JDirection3D & getDirection() const
Get direction.
Definition: JDirection3D.hh:107
JFIT::JRegressor< JShower3EZ, JGandalf >::JPDFMaplist_t
JTOOLS::JMapList< JTOOLS::JPolint0FunctionalMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalGridMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalGridMap > > > > JPDFMaplist_t
Definition: JShower3EZRegressor.hh:285
JMATH::zero
static const JZero zero
Function object to assign zero value.
Definition: JZero.hh:105
JTOOLS::JPolint1FunctionalMapH
Type definition of a 1st degree polynomial interpolation based on a JMap implementation.
Definition: JFunctionalMap_t.hh:145
JFIT::JRegressor< JShower3EZ, JSimplex >::JNPEMaplist_t
JTOOLS::JMapList< JTOOLS::JPolint0FunctionalMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalGridMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalGridMap > > > > JNPEMaplist_t
Definition: JShower3EZRegressor.hh:107
JZero.hh
Definition of zero value for any class.
JFIT::JRegressor< JShower3EZ, JGandalf >::pdf_t
static const JFIT::JPDFType_t pdf_t[NUMBER_OF_PDFS]
Definition: JShower3EZRegressor.hh:467
JGEOMETRY3D::JVersor3Z::getDY
double getDY() const
Get y direction.
Definition: JVersor3Z.hh:158
JFunctionalMap_t.hh
Various implementations of functional maps.
JFIT::getPMTAngle
double getPMTAngle(const double angle)
Constrain PMT angle to [0,pi].
Definition: JShower3EZRegressor.hh:55
JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::JNPEMaplist_t
JTOOLS::JMapList< JTOOLS::JPolint0FunctionalMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalGridMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalGridMap > > > > JNPEMaplist_t
Definition: JShower3EZRegressor.hh:493
JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::T_ns
static JTimeRange T_ns
Time window with respect to Cherenkov hypothesis [ns].
Definition: JShower3EZRegressor.hh:677
JPDFTypes.hh
Numbering scheme for PDF types.
JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::getH0
JNPE_t::result_type getH0(const double R_Hz) const
Get background hypothesis value for time integrated PDF.
Definition: JShower3EZRegressor.hh:595
JFIT::JRegressor< JShower3EZ, JSimplex >::operator()
double operator()(const JShower3EZ &shower, const JPMTW0 &pmt) const
Fit function.
Definition: JShower3EZRegressor.hh:169
JTOOLS::JPolint1FunctionalGridMap
Type definition of a 1st degree polynomial interpolation based on a JGridMap implementation.
Definition: JFunctionalMap_t.hh:91
JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::pdf_t
static const JPDFType_t pdf_t[NUMBER_OF_PDFS]
Definition: JShower3EZRegressor.hh:682
JFIT::JShower3EZ::getE
double getE() const
Get E.
Definition: JShower3EZ.hh:86
JFunction1D_t.hh
JFIT::JShower3EZ
Data structure for fit of straight line in positive z-direction with energy.
Definition: JShower3EZ.hh:27
JLANG::JSharedPointer
The template JSharedPointer class can be used to share a pointer to an object.
Definition: JSharedPointer.hh:28
JGEOMETRY3D::JVector3D
Data structure for vector in three dimensions.
Definition: JVector3D.hh:34
JTOOLS::JMapList
Map list.
Definition: JMapList.hh:24
JFIT::JRegressor< JShower3EZ, JGandalf >::JNPE_t
JPHYSICS::JNPETable< double, double, JNPEMaplist_t > JNPE_t
Definition: JShower3EZRegressor.hh:292
JTOOLS::JPolint1FunctionalMap
Type definition of a 1st degree polynomial interpolation based on a JMap implementation.
Definition: JFunctionalMap_t.hh:55
JFIT::JRegressor< JShower3EZ, JSimplex >::Vmax_npe
static double Vmax_npe
Maximal integral of PDF [npe].
Definition: JShower3EZRegressor.hh:260
JFIT::JRegressor< JShower3EZ, JGandalf >::getH1
JNPE_t::result_type getH1(const double D, const double cd, const double theta, const double phi, const double E) const
Get signal hypothesis value for time integrated PDF.
Definition: JShower3EZRegressor.hh:431
JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::operator()
double operator()(const JShower3EZ &shower, const JPMTW0 &pmt) const
Fit function.
Definition: JShower3EZRegressor.hh:555
NOTICE
#define NOTICE(A)
Definition: JMessage.hh:64
ERROR
#define ERROR(A)
Definition: JMessage.hh:66
JConstants.hh
Physics constants.
JPHYSICS::DIRECT_LIGHT_FROM_EMSHOWER
direct light from EM shower
Definition: JPDFTypes.hh:37
JMATH::PI
static const double PI
Mathematical constants.
Definition: JMath/JConstants.hh:20
JTOOLS::JRange::constrain
T constrain(argument_type x) const
Constrain value to range.
Definition: JRange.hh:350
JFIT::JRegressor< JShower3EZ, JSimplex >::JPDFMaplist_t
JTOOLS::JMapList< JTOOLS::JPolint0FunctionalMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalGridMap, JTOOLS::JMapList< JTOOLS::JPolint0FunctionalGridMap > > > > JPDFMaplist_t
Definition: JShower3EZRegressor.hh:101
JGEOMETRY3D::JPosition3D::getPosition
const JPosition3D & getPosition() const
Get position.
Definition: JPosition3D.hh:130
JPHYSICS::JPDFType_t
JPDFType_t
PDF types.
Definition: JPDFTypes.hh:24
JFIT::JRegressor< JShower3EZ, JSimplex >::estimator
JLANG::JSharedPointer< JMEstimator > estimator
M-Estimator function.
Definition: JShower3EZRegressor.hh:268
JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::getH1
JNPE_t::result_type getH1(const JShower3EZ &shower, const JPMTW0 &pmt) const
Get signal hypothesis value for time integrated PDF.
Definition: JShower3EZRegressor.hh:651
JFIT::JRegressor< JShower3EZ, JSimplex >::getH1
JNPE_t::result_type getH1(const double D, const double cd, const double theta, const double phi, const double E) const
Get signal hypothesis value for time integrated PDF.
Definition: JShower3EZRegressor.hh:224
JMessage.hh
General purpose messaging.
JFIT::JShower3Z
Data structure for cascade in positive z-direction.
Definition: JShower3Z.hh:32
JGEOMETRY3D::JVersor3Z::getDirection
const JVersor3Z & getDirection() const
Get direction.
Definition: JVersor3Z.hh:81
FATAL
#define FATAL(A)
Definition: JMessage.hh:67
JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::Vmax_npe
static double Vmax_npe
Maximal integral of PDF [npe].
Definition: JShower3EZRegressor.hh:678
JFIT::JGandalf
Fit method based on the Levenberg-Marquardt method.
Definition: JGandalf.hh:84
range
z range($ZMAX-$ZMIN)< $MINIMAL_DZ." fi fi typeset -Z 4 STRING typeset -Z 2 FLOOR JPlot1D -f $
Definition: module-Z:fit.sh:118
JFIT::JRegressor< JShower3EZ, JGandalf >::JRegressor
JRegressor(const std::string &fileDescriptor)
Parameterized constructor.
Definition: JShower3EZRegressor.hh:303
JTOOLS::JMultiFunction::setExceptionHandler
void setExceptionHandler(const typename function_type::supervisor_type &supervisor)
Set the supervisor for handling of exceptions.
Definition: JMultiFunction.hh:157
d
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonStart.sh:47
JLANG::JException::what
virtual const char * what() const override
Get error message.
Definition: JException.hh:64
JFIT::JPMTW0::getN
int getN() const
Get number of hits.
Definition: JPMTW0.hh:67
JRange.hh
Auxiliary class to define a range between two values.
JFIT::JSimplex
Simple fit method based on Powell&#39;s algorithm, see reference: Numerical Recipes in C++...
Definition: JRegressorHelper.hh:39
JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::JPDF_t
JPHYSICS::JPDFTable< JFunction1D_t, JPDFMaplist_t > JPDF_t
Definition: JShower3EZRegressor.hh:488
JGEOMETRY3D::JVersor3Z::getDX
double getDX() const
Get x direction.
Definition: JVersor3Z.hh:147
JTOOLS::JPolint0FunctionalGridMap
Type definition of a zero degree polynomial interpolation based on a JGridMap implementation.
Definition: JFunctionalMap_t.hh:82
JFIT::JShower3EZ::getBy
double getBy() const
Get bjorken y.
Definition: JShower3EZ.hh:76
JFIT::JRegressor< JShower3EZ, JGandalf >::operator()
result_type operator()(const JShower3EZ &shower, const JPMTW0 &pmt) const
Fit function.
Definition: JShower3EZRegressor.hh:353
JEEP::getFilename
std::string getFilename(const std::string &file_name)
Get file name part, i.e. part after last JEEP::PATHNAME_SEPARATOR if any.
Definition: JeepToolkit.hh:128
JFIT::JRegressor< JShower3EZ, JGandalf >::Vmax_npe
static double Vmax_npe
Maximal integral of PDF [npe].
Definition: JShower3EZRegressor.hh:463
JFIT::JRegressor< JShower3EZ, JGandalf >::estimator
JLANG::JSharedPointer< JMEstimator > estimator
M-Estimator function.
Definition: JShower3EZRegressor.hh:471
JFIT::JRegressor< JShower3EZ, JGandalf >::T_ns
static JTimeRange T_ns
Time window with respect to Cherenkov hypothesis [ns].
Definition: JShower3EZRegressor.hh:462
JFIT::JRegressor< JShower3EZ, JSimplex >::getH0
JNPE_t::result_type getH0(const double R_Hz) const
Get background hypothesis value for time integrated PDF.
Definition: JShower3EZRegressor.hh:209
JGEOMETRY3D::JPosition3D
Data structure for position in three dimensions.
Definition: JPosition3D.hh:36
JFIT::model
void model(JModel_t &value)
Auxiliary function to constrain model during fit.
Definition: JGandalf.hh:56
JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::JNPE_t
JPHYSICS::JNPETable< double, double, JNPEMaplist_t > JNPE_t
Definition: JShower3EZRegressor.hh:494
JTOOLS::u
double u[N+1]
Definition: JPolint.hh:865
JFIT::getChi2
double getChi2(const double P)
Get chi2 corresponding to given probability.
Definition: JFitToolkit.hh:56
JGEOMETRY3D::JVersor3Z
Data structure for normalised vector in positive z-direction.
Definition: JVersor3Z.hh:39
JFIT::JRegressor< JShower3EZ, JSimplex >::JPDF_t
JPHYSICS::JPDFTable< JFunction1D_t, JPDFMaplist_t > JPDF_t
Definition: JShower3EZRegressor.hh:102
JFIT::JAbstractRegressor
Abstract class for global fit method.
Definition: JRegressor.hh:77
D
do echo Generating $dir eval D
Definition: JDrawLED.sh:53
JTOOLS::get_value
JResultEvaluator< JResult_t >::result_type get_value(const JResult_t &value)
Helper method to recursively evaluate a to function value.
Definition: JResult.hh:998
JFIT::JRegressor< JShower3EZ, JSimplex >::T_ns
static JTimeRange T_ns
Time window with respect to Cherenkov hypothesis [ns].
Definition: JShower3EZRegressor.hh:259
exp
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable STRING $argv[2] set_array QUANTILES set_variable FORMULA *[0] exp(-0.5 *(x-[1])*(x-[1])/([2]*[2]))" set_variable MODULE `getModule -a $DETECTOR -L "$STRING 0"` source JAcousticsToolkit.sh typeset -A TRIPODS get_tripods $WORKDIR/tripod.txt TRIPODS XMEAN
Definition: toashort_monitor.sh:46
JFIT::JRegressor< JShower3EZ, JGandalf >::getH0
JNPE_t::result_type getH0(const double R_Hz) const
Get background hypothesis value for time integrated PDF.
Definition: JShower3EZRegressor.hh:416
JFIT::JAbstractMinimiser
Abstract minimiser.
Definition: JRegressor.hh:25
JMEstimator.hh
Maximum likelihood estimator (M-estimators).
JTOOLS::JPolint0FunctionalMap
Type definition of a zero degree polynomial interpolation based on a JMap implementation.
Definition: JFunctionalMap_t.hh:46
JFIT::JRegressor< JShower3EZ, JGandalf >::JNPEMaplist_t
JTOOLS::JMapList< JTOOLS::JPolint1FunctionalMap, JTOOLS::JMapList< JTOOLS::JPolint1FunctionalMapH, JTOOLS::JMapList< JTOOLS::JPolint1FunctionalGridMap, JTOOLS::JMapList< JTOOLS::JPolint1FunctionalGridMap > > > > JNPEMaplist_t
Definition: JShower3EZRegressor.hh:291
GAUSS_LEGENDRE::epsilon
const double epsilon
Definition: JQuadrature.cc:21
JACOUSTICS::JTimeRange
JTOOLS::JRange< double > JTimeRange
Type definition for time range (unit [s]).
Definition: JAcoustics/JTimeRange.hh:19
JFIT::JRegressor< JShower3EZ, JSimplex >::JNPE_t
JPHYSICS::JNPETable< double, double, JNPEMaplist_t > JNPE_t
Definition: JShower3EZRegressor.hh:108
JFIT::JRegressor< JShower3EZ, JAbstractMinimiser >::estimator
JLANG::JSharedPointer< JMEstimator > estimator
M-Estimator function.
Definition: JShower3EZRegressor.hh:686
JFIT::JRegressor< JShower3EZ, JSimplex >::JRegressor
JRegressor(const std::string &fileDescriptor)
Parameterized constructor.
Definition: JShower3EZRegressor.hh:119
Generated by   doxygen 1.8.5