Jpp 20.0.0-rc.2
the software that should make you happy
Loading...
Searching...
No Matches
JFIT::JNPEHit Struct Reference

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

#include <JNPEHit.hh>

Inheritance diagram for JFIT::JNPEHit:
JFIT::JNPE JFIT::JK40

Public Member Functions

 JNPEHit ()
 Default constructor.
 
 JNPEHit (const JNPE &npe, const size_t ns, const double ps)
 Constructor.
 
size_t getN () const
 Get number of hits.
 
double getPS () const
 Get survival probability of one photo-electron.
 
double getP (const double E_GeV) const
 Get muon energy dependent probability.
 
double getChi2 (const double E_GeV) const
 Get chi2.
 
double getY1 () const
 Get light yield due to minimum ionizing particle.
 
double getYA () const
 Get light yield due to delta-rays.
 
double getYB () const
 Get light yield due to bremsstrahlung.
 
double getZ () const
 Get position along muon path.
 
double getH1 (const double E_GeV) const
 Expected number of photo-electrons for muon hypothesis as a function of muon energy.
 
double getY0 () const
 Get light yield due to random background.
 
double getH0 () const
 Expected number of photo-electrons for random background hypothesis.
 
double getP (const bool hit) const
 Get probability for observing a hit or not.
 
double getChi2 (const bool hit) const
 Get chi2 for observing a hit or not.
 

Protected Attributes

size_t ns
 number of hits
 
double ps
 survival probability of one photo-electron
 
double y1
 light yield due to minimum ionizing particle [npe]
 
double yA
 light yield due to delta-rays [npe*m/GeV]
 
double yB
 light yield due to bremsstrahlung [npe/GeV]
 
double z
 position along muon path [m]
 
double y0
 light yield due to random background [npe]
 

Detailed Description

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

Definition at line 19 of file JNPEHit.hh.

Constructor & Destructor Documentation

◆ JNPEHit() [1/2]

JFIT::JNPEHit::JNPEHit ( )
inline

Default constructor.

Definition at line 25 of file JNPEHit.hh.

25 :
26 JNPE(),
27 ns (0),
28 ps (1.0)
29 {}
size_t ns
number of hits
Definition JNPEHit.hh:102
double ps
survival probability of one photo-electron
Definition JNPEHit.hh:103
JNPE()
Default constructor.
Definition JNPE.hh:35

◆ JNPEHit() [2/2]

JFIT::JNPEHit::JNPEHit ( const JNPE & npe,
const size_t ns,
const double ps )
inline

Constructor.

Parameters
npenpe
nsnumber of hits
pssurvival probability of one photo-electron

Definition at line 39 of file JNPEHit.hh.

41 :
42 JNPE(npe),
43 ns (ns),
44 ps (ps)
45 {}

Member Function Documentation

◆ getN()

size_t JFIT::JNPEHit::getN ( ) const
inline

Get number of hits.

Returns
number of hits

Definition at line 53 of file JNPEHit.hh.

54 {
55 return ns;
56 }

◆ getPS()

double JFIT::JNPEHit::getPS ( ) const
inline

Get survival probability of one photo-electron.

Returns
survival probability of one photo-electron

Definition at line 64 of file JNPEHit.hh.

65 {
66 return ps;
67 }

◆ getP() [1/2]

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

Get muon energy dependent probability.

Parameters
E_GeVenergy [GeV]
Returns
probability

Definition at line 76 of file JNPEHit.hh.

77 {
78 const double x = this->getH1(E_GeV) + this->getH0(); // expectation value for total number of photo-electrons
79
80 const double p0 = exp(-x) * (1.0 + // Poisson probability of 0 photo-electrons OR
81 x * (1.0 - ps)); // Poisson probability of 1 photo-electron AND NOT survive
82
83 if (this->getN() > 0)
84 return 1.0 - p0;
85 else
86 return p0;
87 }
double getH0() const
Expected number of photo-electrons for random background hypothesis.
Definition JK40.hh:60
size_t getN() const
Get number of hits.
Definition JNPEHit.hh:53
double getH1(const double E_GeV) const
Expected number of photo-electrons for muon hypothesis as a function of muon energy.
Definition JNPE.hh:116

◆ getChi2() [1/2]

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

Get chi2.

Parameters
E_GeVenergy [GeV]
Returns
chi2

Definition at line 96 of file JNPEHit.hh.

97 {
98 return -log(this->getP(E_GeV));
99 }
double getP(const double E_GeV) const
Get muon energy dependent probability.
Definition JNPEHit.hh:76

◆ getY1()

double JFIT::JNPE::getY1 ( ) const
inlineinherited

Get light yield due to minimum ionizing particle.

Returns
light yield [npe]

Definition at line 71 of file JNPE.hh.

72 {
73 return y1;
74 }
double y1
light yield due to minimum ionizing particle [npe]
Definition JNPE.hh:131

◆ getYA()

double JFIT::JNPE::getYA ( ) const
inlineinherited

Get light yield due to delta-rays.

Returns
light yield [npe*m/GeV]

Definition at line 82 of file JNPE.hh.

83 {
84 return yA;
85 }
double yA
light yield due to delta-rays [npe*m/GeV]
Definition JNPE.hh:132

◆ getYB()

double JFIT::JNPE::getYB ( ) const
inlineinherited

Get light yield due to bremsstrahlung.

Returns
light yield [npe/GeV]

Definition at line 93 of file JNPE.hh.

94 {
95 return yB;
96 }
double yB
light yield due to bremsstrahlung [npe/GeV]
Definition JNPE.hh:133

◆ getZ()

double JFIT::JNPE::getZ ( ) const
inlineinherited

Get position along muon path.

Returns
position along muon path [m]

Definition at line 104 of file JNPE.hh.

105 {
106 return z;
107 }
double z
position along muon path [m]
Definition JNPE.hh:134

◆ 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_GeVenergy [GeV]
Returns
light yield [npe]

Definition at line 116 of file JNPE.hh.

117 {
118 using namespace JPP;
119
120 const double E = gWater.getE(E_GeV, this->getZ());
121
122 if (E >= MASS_MUON * INDEX_OF_REFRACTION_WATER)
123 return (this->getY1() +
124 this->getYA() * getDeltaRaysFromMuon(E) +
125 this->getYB() * E);
126 else
127 return 0.0;
128 }
double getDeltaRaysFromMuon(const double E, const JRange< double > T_GeV=JRange< double >(DELTARAY_TMIN, DELTARAY_TMAX))
Equivalent EM-shower energy due to delta-rays per unit muon track length.
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
double getYB() const
Get light yield due to bremsstrahlung.
Definition JNPE.hh:93
double getY1() const
Get light yield due to minimum ionizing particle.
Definition JNPE.hh:71
double getZ() const
Get position along muon path.
Definition JNPE.hh:104
double getYA() const
Get light yield due to delta-rays.
Definition JNPE.hh:82

◆ getY0()

double JFIT::JK40::getY0 ( ) const
inlineinherited

Get light yield due to random background.

Returns
light yield [npe]

Definition at line 49 of file JK40.hh.

50 {
51 return y0;
52 }
double y0
light yield due to random background [npe]
Definition JK40.hh:90

◆ getH0()

double JFIT::JK40::getH0 ( ) const
inlineinherited

Expected number of photo-electrons for random background hypothesis.

Returns
light yield [npe]

Definition at line 60 of file JK40.hh.

61 {
62 return this->getY0();
63 }
double getY0() const
Get light yield due to random background.
Definition JK40.hh:49

◆ getP() [2/2]

double JFIT::JK40::getP ( const bool hit) const
inlineinherited

Get probability for observing a hit or not.

Parameters
hithit
Returns
probability

Definition at line 72 of file JK40.hh.

73 {
74 return JFIT::getP(this->getH0(), hit);
75 }
double getP(const double expval, bool hit)
Get Poisson probability to observe a hit or not for given expectation value for the number of hits.

◆ getChi2() [2/2]

double JFIT::JK40::getChi2 ( const bool hit) const
inlineinherited

Get chi2 for observing a hit or not.

Parameters
hithit
Returns
probability

Definition at line 84 of file JK40.hh.

85 {
86 return JFIT::getChi2(this->getH0(), hit);
87 }
double getChi2(const double P)
Get chi2 corresponding to given probability.

Member Data Documentation

◆ ns

size_t JFIT::JNPEHit::ns
protected

number of hits

Definition at line 102 of file JNPEHit.hh.

◆ ps

double JFIT::JNPEHit::ps
protected

survival probability of one photo-electron

Definition at line 103 of file JNPEHit.hh.

◆ y1

double JFIT::JNPE::y1
protectedinherited

light yield due to minimum ionizing particle [npe]

Definition at line 131 of file JNPE.hh.

◆ yA

double JFIT::JNPE::yA
protectedinherited

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

Definition at line 132 of file JNPE.hh.

◆ yB

double JFIT::JNPE::yB
protectedinherited

light yield due to bremsstrahlung [npe/GeV]

Definition at line 133 of file JNPE.hh.

◆ z

double JFIT::JNPE::z
protectedinherited

position along muon path [m]

Definition at line 134 of file JNPE.hh.

◆ y0

double JFIT::JK40::y0
protectedinherited

light yield due to random background [npe]

Definition at line 90 of file JK40.hh.


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