#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. More...

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

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. More...

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

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

Detailed Description

Definition at line 24 of file JNPE_t.hh.

Member Typedef Documentation

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

Definition at line 28 of file JNPE_t.hh.

typedef 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 ( const std::string & fileDescriptor )

inline

Constructor.

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

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_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 string file_name = getFilename(fileDescriptor, pdf_t[i]);
 
       cout << "loading PDF from file " << file_name << "... " << flush;
 
       pdf.load(file_name.c_str());
 
       cout << "OK" << endl;
 
       pdf.setExceptionHandler(supervisor);
 
       if (!is_bremsstrahlung(pdf_t[i]))
         Y1.push_back(JNPE_t(pdf));
       else
         YB.push_back(JNPE_t(pdf));
     }
 
     // Add PDFs
 
     cout << "adding PDFs... " << flush;
 
     Y1[1].add(Y1[0]); Y1.erase(Y1.begin());
     YB[1].add(YB[0]); YB.erase(YB.begin());
 
     cout << "OK" << endl;
   }

Member Function Documentation

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 102 of file JNPE_t.hh.

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

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 131 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;
   }