Auxiliary class for handling various light yields. More...

#include <JNPE.hh>

Inheritance diagram for JFIT::JNPE:

Public Types
typedef std::vector< double >	JL1_t
	Type definition of expectation values of number of multiple hits. More...

Public Member Functions
	JNPE ()
	Default constructor. More...

	JNPE (const JK40 &k40, const double y1, const double yA, const double yB, const double z)
	Constructor. More...

double	getY1 () const
	Get light yield due to muon itself. More...

double	getYA () const
	Get light yield due to delta-rays. More...

double	getYB () const
	Get light yield due to bremsstrahlung. More...

double	getZ () const
	Get position along muon path. More...

double	getH0 () const
	Expected number of photo-electrons for random background hypothesis. More...

double	getH0 (const size_t M) const
	Expected number of photo-electrons for random background hypothesis. More...

double	getH1 (const double E_GeV) const
	Expected number of photo-electrons for muon hypothesis as a function of muon energy. More...

double	getP (const double E_GeV, const bool hit) const
	Get probability for observing a hit or not as a function of muon energy. More...

double	getChi2 (const double E_GeV, const bool hit) const
	Get chi2 for observing a hit or not as a function of muon energy. More...

double	getY0 () const
	Get expectation value for number of single hits. More...

double	getY1 (const size_t M) const
	Get expectation value for number of multiple hits. More...

double	getP (size_t M) const
	Get probability to observe a hit with given multiplicity. More...

Protected Attributes
double	y1
	light yield due to minimum ionizing particle [npe] More...

double	yA
	light yield due to delta-rays [npe*m/GeV] More...

double	yB
	light yield due to bremsstrahlung [npe/GeV] More...

double	z
	position along muon path [m] More...

double	__y0
	expectation value of number of single hits More...

JL1_t	__y1
	expectation values of number of multiple hits More...

Detailed Description

Auxiliary class for handling various light yields.

Note that the effective light yield due to bremsstrahlung is proportional to the muon energy. The position along the muon path is used to correct for the energy loss of the muon between the reference position (z = 0) and the point of emission of the light.

Definition at line 28 of file JNPE.hh.

Member Typedef Documentation

typedef std::vector<double> JFIT::JK40::JL1_t

inherited

Type definition of expectation values of number of multiple hits.

Definition at line 33 of file JK40.hh.

Constructor & Destructor Documentation

JFIT::JNPE::JNPE ( )

inline

Default constructor.

Definition at line 34 of file JNPE.hh.

            :
       JK40(),
       y1(0.0),
       yA(0.0),
       yB(0.0),
       z (0.0)
     {}

JFIT::JNPE::JNPE	(	const JK40 &	k40,
		const double	y1,
		const double	yA,
		const double	yB,
		const double	z
	)

inline

Constructor.

Parameters

k40	K40 light yields [npe]
y1	light yield due to minimum ionizing particle [npe]
yA	light yield due to delta-rays [npe*m/GeV]
yB	light yield due to bremsstrahlung [npe/GeV]
z	position along muon path [m]

Definition at line 52 of file JNPE.hh.

                          :
       JK40(k40),
       y1(y1),
       yA(yA),
       yB(yB),
       z (z)
     {}

Member Function Documentation

double JFIT::JNPE::getY1 ( ) const

inline

Get light yield due to muon itself.

Returns: light yield [npe]

Definition at line 70 of file JNPE.hh.

     {
       return y1;
     }

double JFIT::JNPE::getYA ( ) const

inline

Get light yield due to delta-rays.

Returns: light yield [npe*m/GeV]

Definition at line 81 of file JNPE.hh.

     {
       return yA;
     }

double JFIT::JNPE::getYB ( ) const

inline

Get light yield due to bremsstrahlung.

Returns: light yield [npe/GeV]

Definition at line 92 of file JNPE.hh.

     {
       return yB;
     }

double JFIT::JNPE::getZ ( ) const

inline

Get position along muon path.

Returns: position along muon path [m]

Definition at line 103 of file JNPE.hh.

     {
       return z;
     }

double JFIT::JNPE::getH0 ( ) const

inline

Expected number of photo-electrons for random background hypothesis.

Returns: light yield [npe]

Definition at line 114 of file JNPE.hh.

     {
       return this->JK40::getY0();
     }

double JFIT::JNPE::getH0 ( const size_t M ) const

inline

Expected number of photo-electrons for random background hypothesis.

Parameters

M	multiplicity

Returns: light yield [npe]

Definition at line 126 of file JNPE.hh.

     {
       return this->JK40::getY1(M);
     }

double JFIT::JNPE::getH1 ( const double E_GeV ) const

inline

Expected number of photo-electrons for muon hypothesis as a function of muon energy.

Parameters

E_GeV energy [GeV]

Returns: light yield [npe]

Definition at line 138 of file JNPE.hh.

     {
       using namespace JPP;
 
       const double E = gWater.getE(E_GeV, this->getZ());
 
       if (E >= MASS_MUON * INDEX_OF_REFRACTION_WATER)
         return (this->getY1()                            + 
                 this->getYA() * getDeltaRaysFromMuon(E)  +
                 this->getYB() * E);
       else
         return 0.0;
     }

double JFIT::JNPE::getP	(	const double	E_GeV,
		const bool	hit
	)		const

inline

Get probability for observing a hit or not as a function of muon energy.

Parameters

E_GeV	energy [GeV]
hit	hit

Returns: probability

Definition at line 160 of file JNPE.hh.

     {
       return JFIT::getP(this->getH1(E_GeV) + this->getH0(), hit);
     }

double JFIT::JNPE::getChi2	(	const double	E_GeV,
		const bool	hit
	)		const

inline

Get chi2 for observing a hit or not as a function of muon energy.

Parameters

E_GeV	energy [GeV]
hit	hit

Returns: probability

Definition at line 173 of file JNPE.hh.

     {
       return JFIT::getChi2(this->getH1(E_GeV) + this->getH0(), hit);
     }

double JFIT::JK40::getY0 ( ) const

inlineinherited

Get expectation value for number of single hits.

Returns: expectation value

Definition at line 100 of file JK40.hh.

     {
       return __y0;
     }

double JFIT::JK40::getY1 ( const size_t M ) const

inlineinherited

Get expectation value for number of multiple hits.

Parameters

M	multiplicity (M >= 2)

Returns: expectation value

Definition at line 112 of file JK40.hh.

     {
       if (M >= 2 && M - 2 < __y1.size()) 
         return __y1[M-2];
       else
         return 0.0;
     }

double JFIT::JK40::getP ( size_t M ) const

inlineinherited

Get probability to observe a hit with given multiplicity.

Note that if M = 0, the probability corresponds to observing no hits.

Parameters

M	multiplicity (M >= 0)

Returns: probability

Definition at line 128 of file JK40.hh.

     {
       double expval = 0.0;
       
       if (M < 2)
         expval = getY0();
       else
         expval = getY1(M);
 
       return JFIT::getP(expval, M != 0);
     }

Member Data Documentation

double JFIT::JNPE::y1

protected

light yield due to minimum ionizing particle [npe]

Definition at line 179 of file JNPE.hh.

double JFIT::JNPE::yA

protected

light yield due to delta-rays [npe*m/GeV]

Definition at line 180 of file JNPE.hh.

double JFIT::JNPE::yB

protected

light yield due to bremsstrahlung [npe/GeV]

Definition at line 181 of file JNPE.hh.

double JFIT::JNPE::z

protected

position along muon path [m]

Definition at line 182 of file JNPE.hh.

double JFIT::JK40::__y0

protectedinherited

expectation value of number of single hits

Definition at line 141 of file JK40.hh.

JL1_t JFIT::JK40::__y1

protectedinherited

expectation values of number of multiple hits

Definition at line 142 of file JK40.hh.

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

JNPE.hh

Public Types

Public Member Functions

Protected Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation