Jpp  19.0.0
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Public Member Functions | Protected Attributes | List of all members
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. 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 getY1 () const
 Get light yield due to minimum ionizing particle. 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 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 (const bool hit) const
 Get probability for observing a hit or not. 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 getChi2 (const bool hit) const
 Get chi2 for observing a hit or not. More...
 
double getY0 () const
 Get light yield due to random background. More...
 
double getH0 () const
 Expected number of photo-electrons for random background hypothesis. More...
 

Protected Attributes

size_t __n
 
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
 light yield due to random background [npe] More...
 

Detailed Description

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

Definition at line 19 of file JNPEHit.hh.

Constructor & Destructor Documentation

JFIT::JNPEHit::JNPEHit ( )
inline

Default constructor.

Definition at line 28 of file JNPEHit.hh.

28  :
29  JNPE(),
30  __n (0)
31  {}
JFIT::JNPEHit::__n
size_t __n
Definition: JNPEHit.hh:93
JFIT::JNPE::JNPE
JNPE()
Default constructor.
Definition: JNPE.hh:36
JFIT::JNPEHit::JNPEHit ( const JNPE &  npe,
const size_t  n 
)
inline

Constructor.

Parameters
npenpe
nnumber of hits

Definition at line 40 of file JNPEHit.hh.

41  :
42  JNPE(npe),
43  __n (n)
44  {}
JTOOLS::n
const int n
Definition: JPolint.hh:786
JFIT::JNPEHit::__n
size_t __n
Definition: JNPEHit.hh:93
JFIT::JNPE::JNPE
JNPE()
Default constructor.
Definition: JNPE.hh:36

Member Function Documentation

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

Get number of hits.

Returns
number of hits

Definition at line 52 of file JNPEHit.hh.

53  {
54  return __n;
55  }
JFIT::JNPEHit::__n
size_t __n
Definition: JNPEHit.hh:93
double JFIT::JNPEHit::getP ( ) const
inline

Get probability.

Returns
probability

Definition at line 63 of file JNPEHit.hh.

64  {
65  return this->getP(this->getN());
66  }
JFIT::JNPEHit::getN
size_t getN() const
Get number of hits.
Definition: JNPEHit.hh:52
JFIT::JNPEHit::getP
double getP() const
Get probability.
Definition: JNPEHit.hh:63
double JFIT::JNPEHit::getP ( const double  E_GeV) const
inline

Get muon energy dependent probability.

Parameters
E_GeVenergy [GeV]
Returns
probability

Definition at line 75 of file JNPEHit.hh.

76  {
77  return this->getP(E_GeV, this->getN() > 0);
78  }
JFIT::JNPEHit::getN
size_t getN() const
Get number of hits.
Definition: JNPEHit.hh:52
JFIT::JNPEHit::getP
double getP() const
Get probability.
Definition: JNPEHit.hh:63
double JFIT::JNPEHit::getChi2 ( const double  E_GeV) const
inline

Get chi2.

Parameters
E_GeVenergy [GeV]
Returns
chi2

Definition at line 87 of file JNPEHit.hh.

88  {
89  return this->getChi2(E_GeV, this->getN() > 0);
90  }
JFIT::JNPEHit::getN
size_t getN() const
Get number of hits.
Definition: JNPEHit.hh:52
JFIT::JNPEHit::getChi2
double getChi2(const double E_GeV) const
Get chi2.
Definition: JNPEHit.hh:87
double JFIT::JNPE::getY1 ( ) const
inlineinherited

Get light yield due to minimum ionizing particle.

Returns
light yield [npe]

Definition at line 72 of file JNPE.hh.

73  {
74  return y1;
75  }
JFIT::JNPE::y1
double y1
light yield due to minimum ionizing particle [npe]
Definition: JNPE.hh:158
double JFIT::JNPE::getYA ( ) const
inlineinherited

Get light yield due to delta-rays.

Returns
light yield [npe*m/GeV]

Definition at line 83 of file JNPE.hh.

84  {
85  return yA;
86  }
JFIT::JNPE::yA
double yA
light yield due to delta-rays [npe*m/GeV]
Definition: JNPE.hh:159
double JFIT::JNPE::getYB ( ) const
inlineinherited

Get light yield due to bremsstrahlung.

Returns
light yield [npe/GeV]

Definition at line 94 of file JNPE.hh.

95  {
96  return yB;
97  }
JFIT::JNPE::yB
double yB
light yield due to bremsstrahlung [npe/GeV]
Definition: JNPE.hh:160
double JFIT::JNPE::getZ ( ) const
inlineinherited

Get position along muon path.

Returns
position along muon path [m]

Definition at line 105 of file JNPE.hh.

106  {
107  return z;
108  }
JFIT::JNPE::z
double z
position along muon path [m]
Definition: JNPE.hh:161
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 117 of file JNPE.hh.

118  {
119  using namespace JPP;
120 
121  const double E = gWater.getE(E_GeV, this->getZ());
122 
123  if (E >= MASS_MUON * INDEX_OF_REFRACTION_WATER)
124  return (this->getY1() +
125  this->getYA() * getDeltaRaysFromMuon(E) +
126  this->getYB() * E);
127  else
128  return 0.0;
129  }
E
then usage $script< input file >[option[primary[working directory]]] nWhere option can be E
Definition: JMuonPostfit.sh:40
JFIT::JNPE::getYB
double getYB() const
Get light yield due to bremsstrahlung.
Definition: JNPE.hh:94
JPHYSICS::MASS_MUON
static const double MASS_MUON
muon mass [GeV]
Definition: JPhysics/JConstants.hh:69
JPHYSICS::INDEX_OF_REFRACTION_WATER
static const double INDEX_OF_REFRACTION_WATER
Average index of refraction of water corresponding to the group velocity.
Definition: JPhysics/JConstants.hh:47
JPHYSICS::gWater
static const JGeaneWater gWater
Function object for energy loss of muon in sea water.
Definition: JGeane.hh:381
JPHYSICS::getDeltaRaysFromMuon
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.
Definition: JPDFToolkit.hh:260
JFIT::JNPE::getYA
double getYA() const
Get light yield due to delta-rays.
Definition: JNPE.hh:83
JFIT::JNPE::getY1
double getY1() const
Get light yield due to minimum ionizing particle.
Definition: JNPE.hh:72
JPHYSICS::JGeaneWater::getE
virtual double getE(const double E, const double dx) const override
Get energy of muon after specified distance.
Definition: JGeane.hh:257
JFIT::JNPE::getZ
double getZ() const
Get position along muon path.
Definition: JNPE.hh:105
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_GeVenergy [GeV]
hithit
Returns
probability

Definition at line 139 of file JNPE.hh.

140  {
141  return JFIT::getP(this->getH1(E_GeV) + this->getH0(), hit);
142  }
JFIT::JNPE::getH1
double getH1(const double E_GeV) const
Expected number of photo-electrons for muon hypothesis as a function of muon energy.
Definition: JNPE.hh:117
JFIT::JK40::getH0
double getH0() const
Expected number of photo-electrons for random background hypothesis.
Definition: JK40.hh:60
JFIT::getP
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...
Definition: JFitToolkit.hh:41
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  }
JFIT::JK40::getH0
double getH0() const
Expected number of photo-electrons for random background hypothesis.
Definition: JK40.hh:60
JFIT::getP
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...
Definition: JFitToolkit.hh:41
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_GeVenergy [GeV]
hithit
Returns
probability

Definition at line 152 of file JNPE.hh.

153  {
154  return JFIT::getChi2(this->getH1(E_GeV) + this->getH0(), hit);
155  }
JFIT::JNPE::getH1
double getH1(const double E_GeV) const
Expected number of photo-electrons for muon hypothesis as a function of muon energy.
Definition: JNPE.hh:117
JFIT::JK40::getH0
double getH0() const
Expected number of photo-electrons for random background hypothesis.
Definition: JK40.hh:60
JFIT::getChi2
double getChi2(const double P)
Get chi2 corresponding to given probability.
Definition: JFitToolkit.hh:56
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  }
JFIT::JK40::getH0
double getH0() const
Expected number of photo-electrons for random background hypothesis.
Definition: JK40.hh:60
JFIT::getChi2
double getChi2(const double P)
Get chi2 corresponding to given probability.
Definition: JFitToolkit.hh:56
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  }
JFIT::JK40::y0
double y0
light yield due to random background [npe]
Definition: JK40.hh:90
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  }
JFIT::JK40::getY0
double getY0() const
Get light yield due to random background.
Definition: JK40.hh:49

Member Data Documentation

size_t JFIT::JNPEHit::__n
protected

Definition at line 93 of file JNPEHit.hh.

double JFIT::JNPE::y1
protectedinherited

light yield due to minimum ionizing particle [npe]

Definition at line 158 of file JNPE.hh.

double JFIT::JNPE::yA
protectedinherited

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

Definition at line 159 of file JNPE.hh.

double JFIT::JNPE::yB
protectedinherited

light yield due to bremsstrahlung [npe/GeV]

Definition at line 160 of file JNPE.hh.

double JFIT::JNPE::z
protectedinherited

position along muon path [m]

Definition at line 161 of file JNPE.hh.

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:
Generated by   doxygen 1.8.5