Jpp  18.2.1
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::JNPE Struct Reference

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

#include <JNPE.hh>

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

Public Member Functions

 JNPE ()
 Default constructor. More...
 
 JNPE (const double y0, const double y1, const double yA, const double yB, const double z)
 Constructor. 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 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 light yield due to random background. More...
 
double getH0 () const
 Expected number of photo-electrons for random background hypothesis. More...
 
double getP (const bool hit) const
 Get probability for observing a hit or not. More...
 
double getChi2 (const bool hit) const
 Get chi2 for observing a hit or not. 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
 light yield due to random background [npe] 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 30 of file JNPE.hh.

Constructor & Destructor Documentation

JFIT::JNPE::JNPE ( )
inline

Default constructor.

Definition at line 36 of file JNPE.hh.

36  :
37  JK40(),
38  y1(0.0),
39  yA(0.0),
40  yB(0.0),
41  z (0.0)
42  {}
double z
position along muon path [m]
Definition: JNPE.hh:161
double yB
light yield due to bremsstrahlung [npe/GeV]
Definition: JNPE.hh:160
JK40()
Default constructor.
Definition: JK40.hh:29
double y1
light yield due to minimum ionizing particle [npe]
Definition: JNPE.hh:158
double yA
light yield due to delta-rays [npe*m/GeV]
Definition: JNPE.hh:159
JFIT::JNPE::JNPE ( const double  y0,
const double  y1,
const double  yA,
const double  yB,
const double  z 
)
inline

Constructor.

Parameters
y0light yield due to random background [npe]
y1light yield due to minimum ionizing particle [npe]
yAlight yield due to delta-rays [npe*m/GeV]
yBlight yield due to bremsstrahlung [npe/GeV]
zposition along muon path [m]

Definition at line 54 of file JNPE.hh.

58  :
59  JK40(y0),
60  y1(y1),
61  yA(yA),
62  yB(yB),
63  z (z)
64  {}
double z
position along muon path [m]
Definition: JNPE.hh:161
double yB
light yield due to bremsstrahlung [npe/GeV]
Definition: JNPE.hh:160
JK40()
Default constructor.
Definition: JK40.hh:29
double y0
light yield due to random background [npe]
Definition: JK40.hh:90
double y1
light yield due to minimum ionizing particle [npe]
Definition: JNPE.hh:158
double yA
light yield due to delta-rays [npe*m/GeV]
Definition: JNPE.hh:159

Member Function Documentation

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

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  }
double y1
light yield due to minimum ionizing particle [npe]
Definition: JNPE.hh:158
double JFIT::JNPE::getYA ( ) const
inline

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  }
double yA
light yield due to delta-rays [npe*m/GeV]
Definition: JNPE.hh:159
double JFIT::JNPE::getYB ( ) const
inline

Get light yield due to bremsstrahlung.

Returns
light yield [npe/GeV]

Definition at line 94 of file JNPE.hh.

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

Get position along muon path.

Returns
position along muon path [m]

Definition at line 105 of file JNPE.hh.

106  {
107  return z;
108  }
double z
position along muon path [m]
Definition: JNPE.hh:161
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_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 
124  return (this->getY1() +
125  this->getYA() * getDeltaRaysFromMuon(E) +
126  this->getYB() * E);
127  else
128  return 0.0;
129  }
then usage $script< input file >[option[primary[working directory]]] nWhere option can be E
Definition: JMuonPostfit.sh:40
double getYB() const
Get light yield due to bremsstrahlung.
Definition: JNPE.hh:94
static const double MASS_MUON
muon mass [GeV]
static const double INDEX_OF_REFRACTION_WATER
Average index of refraction of water corresponding to the group velocity.
static const JGeaneWater gWater
Function object for energy loss of muon in sea water.
Definition: JGeane.hh:381
double getYA() const
Get light yield due to delta-rays.
Definition: JNPE.hh:83
double getY1() const
Get light yield due to minimum ionizing particle.
Definition: JNPE.hh:72
virtual double getE(const double E, const double dx) const override
Get energy of muon after specified distance.
Definition: JGeane.hh:257
double getDeltaRaysFromMuon(const double E)
Equivalent EM-shower energy due to delta-rays per unit muon track length.
Definition: JPDFToolkit.hh:67
double getZ() const
Get position along muon path.
Definition: JNPE.hh:105
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_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  }
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
double getH0() const
Expected number of photo-electrons for random background hypothesis.
Definition: JK40.hh:60
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
inline

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  }
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
double getH0() const
Expected number of photo-electrons for random background hypothesis.
Definition: JK40.hh:60
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  }
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  }
double getY0() const
Get light yield due to random background.
Definition: JK40.hh:49
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 getH0() const
Expected number of photo-electrons for random background hypothesis.
Definition: JK40.hh:60
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::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 getH0() const
Expected number of photo-electrons for random background hypothesis.
Definition: JK40.hh:60
double getChi2(const double P)
Get chi2 corresponding to given probability.
Definition: JFitToolkit.hh:56

Member Data Documentation

double JFIT::JNPE::y1
protected

light yield due to minimum ionizing particle [npe]

Definition at line 158 of file JNPE.hh.

double JFIT::JNPE::yA
protected

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

Definition at line 159 of file JNPE.hh.

double JFIT::JNPE::yB
protected

light yield due to bremsstrahlung [npe/GeV]

Definition at line 160 of file JNPE.hh.

double JFIT::JNPE::z
protected

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: