JShowerNPE_t Struct Reference

Auxiliary data structure for shower PDF. More...

#include <JNPE_t.hh>

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

Public Member Functions
	JShowerNPE_t (const std::string &fileDescriptor, const int numberOfPoints=0)
	Constructor.
double	calculate (const double E, const double D, const double cd, const double theta, const double phi) const
	Get PDF.

Private Attributes
int	numberOfPoints
JNPE_t	npe
	PDF for shower.
JNPE_t	F [2]
	PDF for shower.

Detailed Description

Auxiliary data structure for shower PDF.

Definition at line 175 of file JNPE_t.hh.

Member Typedef Documentation

JNPEMaplist_t

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

Definition at line 180 of file JNPE_t.hh.

JNPE_t

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

Definition at line 181 of file JNPE_t.hh.

Constructor & Destructor Documentation

JShowerNPE_t()

JShowerNPE_t::JShowerNPE_t ( const std::string & fileDescriptor, const int numberOfPoints = 0 )

inline

Constructor.

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

Parameters

fileDescriptor	PDF file descriptor
numberOfPoints	number of points for shower elongation

Definition at line 192 of file JNPE_t.hh.

                                                      :
    numberOfPoints(numberOfPoints)
  {
    using namespace std;
    using namespace JPP;
    
    const JPDFType_t pdf_t[] = { SCATTERED_LIGHT_FROM_EMSHOWER,
                                 DIRECT_LIGHT_FROM_EMSHOWER };
 
    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) {
 
      const string file_name = getFilename(fileDescriptor, pdf_t[i]);
 
      cout << "loading input from file " << file_name << "... " << flush;
 
      JPDF_t pdf;
 
      pdf.load(file_name.c_str());
 
      pdf.setExceptionHandler(supervisor);
 
      if (npe.empty())
        npe = JNPE_t(pdf);
      else
        npe.add(JNPE_t(pdf));
 
      F[i] = JNPE_t(pdf);
 
      cout << "OK" << endl;
    }
  }

Member Function Documentation

calculate()

double JShowerNPE_t::calculate ( const double E, const double D, const double cd, 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

.

Parameters

E	shower energy at minimum distance of approach [GeV]
D	distance [m]
cd	cosine emission angle
theta	PMT zenith angle [rad]
phi	PMT azimuth angle [rad]

Returns

hypothesis value

Definition at line 249 of file JNPE_t.hh.

  {
    using namespace std;
    using namespace JPP;
 
    double Y = 0.0;
  
    if (numberOfPoints > 0) {
 
      const double W = 1.0 / (double) numberOfPoints;
 
      for (int i = 0; i != numberOfPoints; ++i) {
 
        const double z = geanz.getLength(E, (i + 0.5) / (double) numberOfPoints);
 
        const double __D  = sqrt(D*D - 2.0*(D*cd)*z + z*z);
        const double __cd = (D * cd - z) / __D;
 
        try {
          Y += W * npe (__D, __cd, theta, phi);
        }
        catch(const exception& error) {
          //cerr << error.what() << endl;
        }
      }
 
    } else {
 
      try {
        Y = npe(D, cd, theta, phi);
      }
      catch(const exception& error) {
        //cerr << error.what() << endl;
      }
    }
 
    return E * Y;
  }

Member Data Documentation

numberOfPoints

int JShowerNPE_t::numberOfPoints

private

Definition at line 293 of file JNPE_t.hh.

npe

JNPE_t JShowerNPE_t::npe

private

PDF for shower.

Definition at line 294 of file JNPE_t.hh.

F

JNPE_t JShowerNPE_t::F[2]

private

PDF for shower.

Definition at line 295 of file JNPE_t.hh.

---

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

• JNPE_t.hh