Auxiliary class for simultaneously handling light yields and response of PMT. More...

#include <JNPEHit.hh>

Inheritance diagram for JFIT::JNPEHit:

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

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

	JNPEHit (const JNPE &npe, const size_t n)
	Constructor. More...

size_t	getN () const
	Get number of hits. More...

double	getP () const
	Get probability. More...

double	getP (const double E_GeV) const
	Get muon energy dependent probability. More...

double	getChi2 (const double E_GeV) const
	Get chi2. 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	getP (size_t M) const
	Get probability to observe a hit with given multiplicity. 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	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	getP (size_t M) const
	Get probability to observe a hit with given multiplicity. 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...

Protected Attributes
size_t	__n

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 simultaneously handling light yields and response of PMT.

Definition at line 19 of file JNPEHit.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

◆ JNPEHit() [1/2]

JFIT::JNPEHit::JNPEHit ( )

inline

Default constructor.

Definition at line 28 of file JNPEHit.hh.

               :
       JNPE(),
       __n (0)
     {}

◆ JNPEHit() [2/2]

JFIT::JNPEHit::JNPEHit	(	const JNPE &	npe,
		const size_t	n
	)

inline

Constructor.

Parameters

npe	npe
n	number of hits

Definition at line 40 of file JNPEHit.hh.

                             :
       JNPE(npe),
       __n (n)
     {}

Member Function Documentation

◆ getN()

size_t JFIT::JNPEHit::getN ( ) const

inline

Get number of hits.

Returns: number of hits

Definition at line 52 of file JNPEHit.hh.

     {
       return __n;
     }

◆ getP() [1/6]

double JFIT::JNPEHit::getP ( ) const

inline

Get probability.

Returns: probability

Definition at line 63 of file JNPEHit.hh.

     {
       return this->getP(this->getN());
     }

◆ getP() [2/6]

double JFIT::JNPEHit::getP ( const double E_GeV ) const

inline

Get muon energy dependent probability.

Parameters

E_GeV energy [GeV]

Returns: probability

Definition at line 75 of file JNPEHit.hh.

     {
       return this->getP(E_GeV, this->getN() > 0);
     }

◆ getChi2() [1/3]

double JFIT::JNPEHit::getChi2 ( const double E_GeV ) const

inline

Get chi2.

Parameters

E_GeV energy [GeV]

Returns: chi2

Definition at line 87 of file JNPEHit.hh.

     {
       return this->getChi2(E_GeV, this->getN() > 0);
     }

◆ getP() [3/6]

double JFIT::JNPE::getP

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

◆ getP() [4/6]

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

◆ getChi2() [2/3]

double JFIT::JNPE::getChi2

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

◆ getYA()

double JFIT::JNPE::getYA ( ) const

inlineinherited

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

inlineinherited

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

inlineinherited

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

inlineinherited

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

inlineinherited

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

inlineinherited

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() [5/6]

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

inlineinherited

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

◆ getP() [6/6]

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

◆ getChi2() [3/3]

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

inlineinherited

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

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

Member Data Documentation

◆ __n

size_t JFIT::JNPEHit::__n

protected

Definition at line 93 of file JNPEHit.hh.

◆ __yA

double JFIT::JNPE::__yA

protectedinherited

light yield due to minimum ionizing particle [npe]

Definition at line 166 of file JNPE.hh.

◆ __yB

double JFIT::JNPE::__yB

protectedinherited

light yield due to bremsstrahlung [npe/GeV]

Definition at line 167 of file JNPE.hh.

◆ __z

double JFIT::JNPE::__z

protectedinherited

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:

JNPEHit.hh

Public Types

Public Member Functions

Protected Attributes

Detailed Description

Member Typedef Documentation

◆ JL1_t

Constructor & Destructor Documentation

◆ JNPEHit() [1/2]

◆ JNPEHit() [2/2]

Member Function Documentation

◆ getN()

◆ getP() [1/6]

◆ getP() [2/6]

◆ getChi2() [1/3]

◆ getP() [3/6]

◆ getP() [4/6]

◆ getChi2() [2/3]

◆ getYA()

◆ getYB()

◆ getZ()

◆ getH0() [1/2]

◆ getH0() [2/2]

◆ getH1()

◆ getP() [5/6]

◆ getP() [6/6]

◆ getChi2() [3/3]

◆ getY0()

◆ getY1()

Member Data Documentation

◆ __n

◆ __yA

◆ __yB

◆ __z

◆ __y0

◆ __y1