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. More...
double	calculate (const double E, const double D, const double cd, const double theta, const double phi) const
	Get PDF. More...

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

Detailed Description

Auxiliary data structure for shower PDF.

Definition at line 167 of file JNPE_t.hh.

Member Typedef Documentation

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

Definition at line 172 of file JNPE_t.hh.

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

Definition at line 173 of file JNPE_t.hh.

Constructor & Destructor Documentation

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::WILD_CARD.

Parameters

fileDescriptor	PDF file descriptor
numberOfPoints	number of points for shower elongation

Definition at line 184 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

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 241 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

int JShowerNPE_t::numberOfPoints

private

Definition at line 285 of file JNPE_t.hh.

JNPE_t JShowerNPE_t::npe

private

PDF for shower.

Definition at line 286 of file JNPE_t.hh.

JNPE_t JShowerNPE_t::F[2]

private

PDF for shower.

Definition at line 287 of file JNPE_t.hh.

---

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

• JNPE_t.hh