Auxiliary data structure for muon 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::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
	JPDF (const std::string &file_name, const double TTS, const int numberOfPoints=25, const double epsilon=1.0e-10)
	Constructor.

result_type	calculate (const double E, const double R, const double theta, const double phi, const double t1) const
	Get PDF.

Public Attributes
JPDF_t	pdf
	PDF.

int	type
	PDF type.

Detailed Description

Auxiliary data structure for muon PDF.

Definition at line 26 of file JPDF_t.hh.

Member Typedef Documentation

◆ JFunction1D_t

JPP::JSplineFunction1D<JPP::JSplineElement2S<double, double>, JPP::JCollection, JPP::JResultPDF<double> > JPDF::JFunction1D_t

Definition at line 30 of file JPDF_t.hh.

◆ JPDFMaplist_t

JPP::JMAPLIST<JPP::JPolint1FunctionalMap,JPP::JPolint0FunctionalGridMap,JPP::JPolint0FunctionalGridMap>::maplist JPDF::JPDFMaplist_t

Definition at line 34 of file JPDF_t.hh.

◆ JPDF_t

JPP::JPDFTable<JFunction1D_t, JPDFMaplist_t> JPDF::JPDF_t

Definition at line 35 of file JPDF_t.hh.

◆ result_type

JFunction1D_t::result_type JPDF::result_type

Definition at line 36 of file JPDF_t.hh.

Constructor & Destructor Documentation

◆ JPDF()

JPDF::JPDF	(	const std::string &	file_name,
		const double	TTS,
		const int	numberOfPoints = 25,
		const double	epsilon = 1.0e-10 )

inline

Constructor.

The TTS corresponds to the additional time smearing applied to the PDFs.

Parameters

file_name	file name
TTS	TTS [ns]
numberOfPoints	number of points for Gauss-Hermite integration of TTS
epsilon	precision for Gauss-Hermite integration of TTS

Definition at line 49 of file JPDF_t.hh.

  {
    using namespace std;
    using namespace JPP;
 
    const JPDF_t::JSupervisor supervisor(new JPDF_t::JDefaultResult(zero));
 
    cout << "loading input from file " << file_name << "... " << flush;
 
    pdf.load(file_name.c_str());
 
    pdf.setExceptionHandler(supervisor);
 
    cout << "OK" << endl;
 
    type = getPDFType(file_name);
 
    if        (TTS > 0.0) {
 
      cout << "bluring PDFs... " << flush;
 
      pdf.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 JPDF::calculate	(	const double	E,
		const double	R,
		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

E	muon energy at minimum distance of approach [GeV]
R	minimum distance of approach [m]
theta	PMT zenith angle [rad]
phi	PMT azimuth angle [rad]
t1	arrival time relative to Cherenkov hypothesis [ns]

Returns: hypothesis value

Definition at line 97 of file JPDF_t.hh.

  {
    using namespace JPP;
 
    result_type h1 = pdf(R, theta, phi, t1);
 
    if        (is_bremsstrahlung(type)) {
      h1 *= E;
    } else if (is_deltarays(type)) {
      h1 *= getDeltaRaysFromMuon(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

◆ pdf

JPDF_t JPDF::pdf

PDF.

Definition at line 127 of file JPDF_t.hh.

◆ type

int JPDF::type

PDF type.

Definition at line 128 of file JPDF_t.hh.

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

JPDF_t.hh

Public Types

Public Member Functions

Public Attributes