#include <JNPE_t.hh>

Public Types
typedef JPP::JMAPLIST< JPP::JPolint1FunctionalMap, JPP::JPolint1FunctionalGridMap, JPP::JPolint1FunctionalGridMap >::maplist	JNPEMaplist_t

typedef JPP::JNPETable< double, double, JNPEMaplist_t >	JNPE_t

Public Member Functions
	JMuonNPE_t (const std::string &fileDescriptor)
	Constructor.

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

Static Private Member Functions
static double	getNPE (const std::vector< JNPE_t > &NPE, const double R, const double theta, const double phi)
	Get number of photo-electrons.

Private Attributes
std::vector< JNPE_t >	Y1
	light from muon

std::vector< JNPE_t >	YA
	light from delta-rays

std::vector< JNPE_t >	YB
	light from EM showers

Detailed Description

Definition at line 24 of file JNPE_t.hh.

Member Typedef Documentation

◆ JNPEMaplist_t

JPP::JMAPLIST<JPP::JPolint1FunctionalMap,JPP::JPolint1FunctionalGridMap,JPP::JPolint1FunctionalGridMap>::maplist JMuonNPE_t::JNPEMaplist_t

Definition at line 28 of file JNPE_t.hh.

◆ JNPE_t

JPP::JNPETable<double, double, JNPEMaplist_t> JMuonNPE_t::JNPE_t

Definition at line 29 of file JNPE_t.hh.

Constructor & Destructor Documentation

◆ JMuonNPE_t()

JMuonNPE_t::JMuonNPE_t ( const std::string & fileDescriptor )

inline

Constructor.

The PDF file descriptor should contain the wild card character JPHYSICS::WILDCARD.

Parameters

fileDescriptor PDF file descriptor

Definition at line 39 of file JNPE_t.hh.

  {
    using namespace std;
    using namespace JPP;
 
    const JPDFType_t pdf_t[] = { DIRECT_LIGHT_FROM_MUON,
                                 SCATTERED_LIGHT_FROM_MUON,
                                 DIRECT_LIGHT_FROM_DELTARAYS,
                                 SCATTERED_LIGHT_FROM_DELTARAYS,
                                 DIRECT_LIGHT_FROM_EMSHOWERS,
                                 SCATTERED_LIGHT_FROM_EMSHOWERS };
 
    const  int N = sizeof(pdf_t) / sizeof(pdf_t[0]);
    
    typedef JPP::JSplineFunction1D<JSplineElement2D<double, double>,
                                   JCollection,
                                   double>                          JFunction1D_t;
    typedef JPDFTable<JFunction1D_t, JNPEMaplist_t>                 JPDF_t;
 
 
    const JNPE_t::JSupervisor supervisor(new JNPE_t::JDefaultResult(zero));
 
    for (int i = 0; i != N; ++i) {
 
      JPDF_t pdf;
 
      const JPDFType_t type      = pdf_t[i];
      const string     file_name = getFilename(fileDescriptor, type);
 
      cout << "loading PDF from file " << file_name << "... " << flush;
 
      pdf.load(file_name.c_str());
 
      cout << "OK" << endl;
 
      pdf.setExceptionHandler(supervisor);
 
      if      (is_bremsstrahlung(type))
        YB.push_back(JNPE_t(pdf));
      else if (is_deltarays(type))
        YA.push_back(JNPE_t(pdf));
      else
        Y1.push_back(JNPE_t(pdf));
    }
 
    // Add PDFs
 
    cout << "adding PDFs... " << flush;
 
    Y1[1].add(Y1[0]); Y1.erase(Y1.begin());
    YA[1].add(YA[0]); YA.erase(YA.begin());
    YB[1].add(YB[0]); YB.erase(YB.begin());
  
    cout << "OK" << endl;
  }

Member Function Documentation

◆ calculate()

double JMuonNPE_t::calculate	(	const double	E,
		const double	R,
		const double	theta,
		const double	phi ) 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]

Returns: number of photo-electrons

Definition at line 108 of file JNPE_t.hh.

  {
    using namespace JPP;
 
    const double y1 = getNPE(Y1, R, theta, phi);
    const double yA = getNPE(YA, R, theta, phi);
    const double yB = getNPE(YB, R, theta, phi);
 
    if (E >= MASS_MUON * INDEX_OF_REFRACTION_WATER)
      return y1  +  getDeltaRaysFromMuon(E) * yA  +  E * yB;
    else
      return 0.0;
  }

◆ getNPE()

static double JMuonNPE_t::getNPE	(	const std::vector< JNPE_t > &	NPE,
		const double	R,
		const double	theta,
		const double	phi )

inlinestaticprivate

Get number of photo-electrons.

Parameters

NPE	NPE tables
R	distance between muon and PMT [m]
theta	zenith angle orientation PMT [rad]
phi	azimuth angle orientation PMT [rad]

Returns: number of photo-electrons

Definition at line 139 of file JNPE_t.hh.

  {
    using namespace std;
    using namespace JPP;
 
    double npe = 0.0;
 
    for (vector<JNPE_t>::const_iterator i = NPE.begin(); i != NPE.end(); ++i) {
 
      if (R <= i->getXmax()) {
 
        try {
 
          const double y = get_value((*i)(std::max(R, i->getXmin()), theta, phi));
 
          if (y > 0.0) {
            npe += y;
          }
        }
        catch(const exception& error) {
          cerr << error.what() << endl;
        }
      }
    }
 
    return npe;
  }