JMuonPDF_t Struct Reference

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
	JMuonPDF_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 R, const double theta, const double phi, const double t1) const
	Get PDF. More...

Public Attributes
JPDF_t	pdfA
	PDF for minimum ionising particle. More...
JPDF_t	pdfB
	PDF for average energy losses. More...
JPDF_t	pdfC
	PDF for delta-rays. More...

Detailed Description

Auxiliary data structure for muon PDF.

Definition at line 135 of file JPDF_t.hh.

Member Typedef Documentation

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

Definition at line 139 of file JPDF_t.hh.

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

Definition at line 143 of file JPDF_t.hh.

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

Definition at line 144 of file JPDF_t.hh.

typedef JFunction1D_t::result_type JMuonPDF_t::result_type

Definition at line 145 of file JPDF_t.hh.

Constructor & Destructor Documentation

JMuonPDF_t::JMuonPDF_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 159 of file JPDF_t.hh.

  {
    using namespace std;
    using namespace JPP;

    const JPDFType_t pdf_t[] = { DIRECT_LIGHT_FROM_MUON,
                                 SCATTERED_LIGHT_FROM_MUON,
                                 DIRECT_LIGHT_FROM_EMSHOWERS,
                                 SCATTERED_LIGHT_FROM_EMSHOWERS,
                                 DIRECT_LIGHT_FROM_DELTARAYS,
                                 SCATTERED_LIGHT_FROM_DELTARAYS };

    const  int N = sizeof(pdf_t) / sizeof(pdf_t[0]);
    
    JPDF_t pdf[N];

    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;

      pdf[i].load(file_name.c_str());

      pdf[i].setExceptionHandler(supervisor);

      cout << "OK" << endl;
    }

    // Add PDFs

    cout << "adding PDFs... " << flush;

    pdfA = pdf[1];  pdfA.add(pdf[0]);
    pdfB = pdf[3];  pdfB.add(pdf[2]);
    pdfC = pdf[5];  pdfC.add(pdf[4]);

    cout << "OK" << endl;

    if        (TTS > 0.0) {

      cout << "bluring PDFs... " << flush;

      pdfA.blur(TTS, numberOfPoints, epsilon);
      pdfB.blur(TTS, numberOfPoints, epsilon);
      pdfC.blur(TTS, numberOfPoints, epsilon);

      cout << "OK" << endl;

    } else if (TTS < 0.0) {

      THROW(JValueOutOfRange, "Illegal value of TTS [ns]: " << TTS); 
    }
  }

Member Function Documentation

result_type JMuonPDF_t::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

.

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 233 of file JPDF_t.hh.

  {
    using namespace JPP;

    result_type h1 = (pdfA(R, theta, phi, t1)                               +
                      pdfB(R, theta, phi, t1) * E                           +
                      pdfC(R, theta, phi, t1) * 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

JPDF_t JMuonPDF_t::pdfA

PDF for minimum ionising particle.

Definition at line 259 of file JPDF_t.hh.

JPDF_t JMuonPDF_t::pdfB

PDF for average energy losses.

Definition at line 260 of file JPDF_t.hh.

JPDF_t JMuonPDF_t::pdfC

PDF for delta-rays.

Definition at line 261 of file JPDF_t.hh.

---

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

• JPDF_t.hh