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.
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	pdfA
	PDF for minimum ionising particle.
JPDF_t	pdfB
	PDF for average energy losses.
JPDF_t	pdfC
	PDF for delta-rays.

Detailed Description

Auxiliary data structure for muon PDF.

Definition at line 135 of file JPDF_t.hh.

Member Typedef Documentation

JFunction1D_t

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

Definition at line 139 of file JPDF_t.hh.

JPDFMaplist_t

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

Definition at line 143 of file JPDF_t.hh.

JPDF_t

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

Definition at line 144 of file JPDF_t.hh.

result_type

JFunction1D_t::result_type JMuonPDF_t::result_type

Definition at line 145 of file JPDF_t.hh.

Constructor & Destructor Documentation

JMuonPDF_t()

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

calculate()

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 235 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) * JDeltaRays::getEnergyLossFromMuon(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 JMuonPDF_t::pdfA

PDF for minimum ionising particle.

Definition at line 261 of file JPDF_t.hh.

pdfB

JPDF_t JMuonPDF_t::pdfB

PDF for average energy losses.

Definition at line 262 of file JPDF_t.hh.

pdfC

JPDF_t JMuonPDF_t::pdfC

PDF for delta-rays.

Definition at line 263 of file JPDF_t.hh.

---

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

• JPDF_t.hh