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 yA, const double yB, const double z)
	Constructor. More...

double	getYA () const
	Get light yield due to muon itself. 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	getP (size_t M) const
	Get probability to observe a hit with given multiplicity. 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	__yA
	light yield due to minimum ionizing particle [npe] 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 27 of file JNPE.hh.

Member Typedef Documentation

◆ JL1_t

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

inherited

Type definition of expectation values of number of multiple hits.

Definition at line 34 of file JK40.hh.

Constructor & Destructor Documentation

◆ JNPE() [1/2]

JFIT::JNPE::JNPE ( )

inline

Default constructor.

Definition at line 37 of file JNPE.hh.

            :
       JK40(),
       __yA(0.0),
       __yB(0.0),
       __z (0.0)
     {}

◆ JNPE() [2/2]

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

inline

Constructor.

Parameters

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

Definition at line 53 of file JNPE.hh.

                          :
       JK40(k40),
       __yA(yA),
       __yB(yB),
       __z (z)
     {}

Member Function Documentation

◆ getYA()

double JFIT::JNPE::getYA ( ) const

inline

Get light yield due to muon itself.

Returns: light yield [npe]

Definition at line 69 of file JNPE.hh.

     {
       return __yA;
     }

◆ getYB()

double JFIT::JNPE::getYB ( ) const

inline

Get light yield due to bremsstrahlung.

Returns: light yield [npe/GeV]

Definition at line 80 of file JNPE.hh.

     {
       return __yB;
     }

◆ getZ()

double JFIT::JNPE::getZ ( ) const

inline

Get position along muon path.

Returns: position along muon path [m]

Definition at line 91 of file JNPE.hh.

     {
       return __z;
     }

◆ getH0() [1/2]

double JFIT::JNPE::getH0 ( ) const

inline

Expected number of photo-electrons for random background hypothesis.

Returns: light yield [npe]

Definition at line 102 of file JNPE.hh.

     {
       return this->getY0();
     }

◆ getH0() [2/2]

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 114 of file JNPE.hh.

     {
       return this->getY1(M);
     }

◆ getH1()

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 126 of file JNPE.hh.

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

◆ getP() [1/3]

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 147 of file JNPE.hh.

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

◆ getChi2()

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 160 of file JNPE.hh.

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

◆ getP() [2/3]

double JFIT::JK40::getP

inline

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 129 of file JK40.hh.

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

◆ getY0()

double JFIT::JK40::getY0 ( ) const

inlineinherited

Get expectation value for number of single hits.

Returns: expectation value

Definition at line 101 of file JK40.hh.

     {
       return __y0;
     }

◆ getY1()

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 113 of file JK40.hh.

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

◆ getP() [3/3]

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

◆ __yA

double JFIT::JNPE::__yA

protected

light yield due to minimum ionizing particle [npe]

Definition at line 166 of file JNPE.hh.

◆ __yB

double JFIT::JNPE::__yB

protected

light yield due to bremsstrahlung [npe/GeV]

Definition at line 167 of file JNPE.hh.

◆ __z

double JFIT::JNPE::__z

protected

position along muon path [m]

Definition at line 168 of file JNPE.hh.

◆ __y0

double JFIT::JK40::__y0

protectedinherited

expectation value of number of single hits

Definition at line 142 of file JK40.hh.

◆ __y1

JL1_t JFIT::JK40::__y1

protectedinherited

expectation values of number of multiple hits

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

◆ JL1_t

Constructor & Destructor Documentation

◆ JNPE() [1/2]

◆ JNPE() [2/2]

Member Function Documentation

◆ getYA()

◆ getYB()

◆ getZ()

◆ getH0() [1/2]

◆ getH0() [2/2]

◆ getH1()

◆ getP() [1/3]

◆ getChi2()

◆ getP() [2/3]

◆ getY0()

◆ getY1()

◆ getP() [3/3]

Member Data Documentation

◆ __yA

◆ __yB

◆ __z

◆ __y0

◆ __y1