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_t >	JPDF_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. More...
result_type	calculate (const double E, const double D, const double cd, const double theta, const double phi, const double t1) const
	Get PDF. More...

Public Attributes
JPDF_t	pdfA
	PDF for shower. More...

Detailed Description

Auxiliary data structure for shower PDF.

Definition at line 268 of file JPDF_t.hh.

Member Typedef Documentation

JFunction1D_t

typedef JPP::JSplineFunction1D<JPP::JSplineElement2S<double, double>, JPP::JCollection, JPP::JResultPDF<double> > JShowerPDF_t::JFunction1D_t

Definition at line 272 of file JPDF_t.hh.

JPDFMaplist_t

typedef JPP::JMAPLIST<JPP::JPolint1FunctionalMap, JPP::JPolint1FunctionalMap, JPP::JPolint0FunctionalGridMap, JPP::JPolint0FunctionalGridMap>::maplist JShowerPDF_t::JPDFMaplist_t

Definition at line 277 of file JPDF_t.hh.

JPDF_t

typedef JPP::JPDFTable<JFunction1D_t, JPDFMaplist_t> JShowerPDF_t::JPDF_t

Definition at line 278 of file JPDF_t.hh.

result_type

typedef JFunction1D_t::result_type JShowerPDF_t::result_type

Definition at line 279 of file JPDF_t.hh.

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

fileDescriptor	PDF file descriptor
TTS	TTS [ns]
numberOfPoints	number of points for Gauss-Hermite integration of TTS
epsilon	precision for Gauss-Hermite integration of TTS

Definition at line 293 of file JPDF_t.hh.

  {
    using namespace std;
    using namespace JPP;
 
    const JPDFType_t pdf_t[] = { SCATTERED_LIGHT_FROM_EMSHOWER,
                                 DIRECT_LIGHT_FROM_EMSHOWER };
 
    const  int N = sizeof(pdf_t) / sizeof(pdf_t[0]);
 
    const JPDF_t::JSupervisor supervisor(new JPDF_t::JDefaultResult(zero));
 
    for (int i = 0; i != N; ++i) {
 
      const string file_name = getFilename(fileDescriptor, pdf_t[i]);
 
      cout << "loading input from file " << file_name << "... " << flush;
 
      JPDF_t pdf;
 
      pdf.load(file_name.c_str());
 
      pdf.setExceptionHandler(supervisor);
 
      if (pdfA.empty())
        pdfA = pdf;
      else
        pdfA.add(pdf);
 
      cout << "OK" << endl;
    }
 
    if        (TTS > 0.0) {
 
      cout << "bluring PDFs... " << flush;
 
      pdfA.blur(TTS, numberOfPoints, epsilon);
 
      cout << "OK" << endl;
 
    } else if (TTS < 0.0) {
 
      THROW(JValueOutOfRange, "Illegal value of TTS [ns]: " << TTS);
    }
  }

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

.

Parameters

E	shower energy [GeV]
D	distance [m]
cd	cosine emission angle
theta	PMT zenith angle [rad]
phi	PMT azimuth angle [rad]
t1	arrival time relative to Cherenkov hypothesis [ns]

Returns

hypothesis value

Definition at line 357 of file JPDF_t.hh.

  {
    using namespace JPP;
 
    result_type h1 = pdfA(D, cd, theta, phi, t1) * E;
 
    // safety measures
 
    if (h1.f <= 0.0) {
      h1.f  = 0.0;
      h1.fp = 0.0;
    }
 
    if (h1.v <= 0.0) {
      h1.v  = 0.0;
    }
 
    return h1;
  }

Member Data Documentation

pdfA

JPDF_t JShowerPDF_t::pdfA

PDF for shower.

Definition at line 382 of file JPDF_t.hh.

---

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

• JPDF_t.hh