Jpp test-rotations-old-533-g2bdbdb559
the software that should make you happy
Loading...
Searching...
No Matches
JShowerPDF_t Struct Reference

Auxiliary data structure for shower PDF. More...

#include <JPDF_t.hh>

Public Types

typedef JPP::JSplineFunction1D< JPP::JSplineElement2S< double, double >, JPP::JCollection, JPP::JResultPDF< double > > JFunction1D_t
 
typedef JPP::JMAPLIST< JPP::JPolint1FunctionalMap, JPP::JPolint1FunctionalMap, JPP::JPolint0FunctionalGridMap, JPP::JPolint0FunctionalGridMap >::maplist JPDFMaplist_t
 
typedef JPP::JPDFTable< JFunction1D_t, JPDFMaplist_tJPDF_t
 
typedef JFunction1D_t::result_type result_type
 

Public Member Functions

 JShowerPDF_t (const std::string &fileDescriptor, const double TTS, const int numberOfPoints=25, const double epsilon=1.0e-10)
 Constructor.
 
result_type calculate (const double E, const double D, const double cd, const double theta, const double phi, const double t1) const
 Get PDF.
 

Public Attributes

JPDF_t pdfA
 PDF for shower.
 

Detailed Description

Auxiliary data structure for shower PDF.

Definition at line 270 of file JPDF_t.hh.

Member Typedef Documentation

◆ JFunction1D_t

◆ JPDFMaplist_t

◆ JPDF_t

◆ result_type

Constructor & Destructor Documentation

◆ JShowerPDF_t()

JShowerPDF_t::JShowerPDF_t ( const std::string & fileDescriptor,
const double TTS,
const int numberOfPoints = 25,
const double epsilon = 1.0e-10 )
inline

Constructor.

The PDF file descriptor should contain the wild card character JPHYSICS::WILDCARD.
The TTS corresponds to the additional time smearing applied to the PDFs.

Parameters
fileDescriptorPDF file descriptor
TTSTTS [ns]
numberOfPointsnumber of points for Gauss-Hermite integration of TTS
epsilonprecision for Gauss-Hermite integration of TTS

Definition at line 295 of file JPDF_t.hh.

299 {
300 using namespace std;
301 using namespace JPP;
302
303 const JPDFType_t pdf_t[] = {
304 SCATTERED_LIGHT_FROM_EMSHOWER,
305 DIRECT_LIGHT_FROM_EMSHOWER
306 };
307
308 const int N = sizeof(pdf_t) / sizeof(pdf_t[0]);
309
310 const JPDF_t::JSupervisor supervisor(new JPDF_t::JDefaultResult(zero));
311
312 for (int i = 0; i != N; ++i) {
313
314 const string file_name = getFilename(fileDescriptor, pdf_t[i]);
315
316 cout << "loading input from file " << file_name << "... " << flush;
317
318 JPDF_t pdf;
319
320 pdf.load(file_name.c_str());
321
322 pdf.setExceptionHandler(supervisor);
323
324 if (pdfA.empty())
325 pdfA = pdf;
326 else
327 pdfA.add(pdf);
328
329 cout << "OK" << endl;
330 }
331
332 if (TTS > 0.0) {
333
334 cout << "bluring PDFs... " << flush;
335
336 pdfA.blur(TTS, numberOfPoints, epsilon);
337
338 cout << "OK" << endl;
339
340 } else if (TTS < 0.0) {
341
342 THROW(JValueOutOfRange, "Illegal value of TTS [ns]: " << TTS);
343 }
344 }
#define THROW(JException_t, A)
Marco for throwing exception with std::ostream compatible message.
int numberOfPoints
Definition JResultPDF.cc:22
Exception for accessing a value in a collection that is outside of its range.
void blur(const double TTS, const int numberOfPoints=25, const double epsilon=1.0e-10, const double quantile=0.99)
Blur PDF.
Definition JPDFTable.hh:111
void add(const JMultiFunction_t &input)
Add function.
JPDFType_t
PDF types.
Definition JPDFTypes.hh:24
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
void load(const char *file_name)
Load from input file.
JPDF_t pdfA
PDF for shower.
Definition JPDF_t.hh:386
JPP::JPDFTable< JFunction1D_t, JPDFMaplist_t > JPDF_t
Definition JPDF_t.hh:280

Member Function Documentation

◆ calculate()

result_type JShowerPDF_t::calculate ( const double E,
const double D,
const double cd,
const double theta,
const double phi,
const double t1 ) const
inline

Get PDF.

The orientation of the PMT should be defined according this documentation.
In this, the zenith and azimuth angles are limited to

\[\left[0, \pi\right]\]

.

Parameters
Eshower energy [GeV]
Ddistance [m]
cdcosine emission angle
thetaPMT zenith angle [rad]
phiPMT azimuth angle [rad]
t1arrival time relative to Cherenkov hypothesis [ns]
Returns
hypothesis value

Definition at line 361 of file JPDF_t.hh.

367 {
368 using namespace JPP;
369
370 result_type h1 = pdfA(D, cd, theta, phi, t1) * E;
371
372 // safety measures
373
374 if (h1.f <= 0.0) {
375 h1.f = 0.0;
376 h1.fp = 0.0;
377 }
378
379 if (h1.v <= 0.0) {
380 h1.v = 0.0;
381 }
382
383 return h1;
384 }
JFunction1D_t::result_type result_type
Definition JPDF_t.hh:281

Member Data Documentation

◆ pdfA

JPDF_t JShowerPDF_t::pdfA

PDF for shower.

Definition at line 386 of file JPDF_t.hh.


The documentation for this struct was generated from the following file: