Jpp  18.6.0-rc.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
JPHYSICS::JLED Class Referenceabstract

Probability Density Functions of the time response of a PMT. More...

#include <JLED.hh>

Inheritance diagram for JPHYSICS::JLED:
JPHYSICS::JDispersionInterface JPHYSICS::JAbstractPMT JPHYSICS::JAbstractLED JPHYSICS::JAbstractMedium JPHYSICS::JLED_C

Public Member Functions

 JLED (const double lambda, const double Tmin_ns, const double Tmax_ns, const JQuadrature &engine=JCotangent(20), const int numberOfPoints=20, const double epsilon=1e-12)
 Constructor. More...
 
double getDirectLightFromLED (const double D_m, const double cd, const double theta, const double phi, const double t_ns) const
 Probability density function for direct light from LED. More...
 
double getScatteredLightFromLED (const double D_m, const double cd, const double theta, const double phi, const double t_ns) const
 Probability density function for scattered light from LED. More...
 
virtual double getIndexOfRefractionPhase (const double lambda) const =0
 Index of refraction for phase velocity. More...
 
virtual double getDispersionPhase (const double lambda) const =0
 Dispersion of light for phase velocity. More...
 
virtual double getIndexOfRefractionGroup (const double lambda) const
 Index of refraction for group velocity. More...
 
virtual double getDispersionGroup (const double lambda) const =0
 Dispersion of light for group velocity. More...
 
double getKappa (const double lambda) const
 Get effective index of refraction for muon light. More...
 
double getKmin (const double lambda) const
 Get smallest index of refraction for Bremsstrahlung light (i.e. point at which dt/dz = 0). More...
 
virtual double getPhotocathodeArea () const =0
 Photo-cathode area of PMT. More...
 
virtual double getQE (const double lambda) const =0
 Quantum efficiency of PMT (incl. More...
 
virtual double getAngularAcceptance (const double ct) const =0
 Angular acceptence of PMT. More...
 
virtual double getLightFromLED (const double ct, const double phi, const double dt) const =0
 Light yield from LED (number of p.e. More...
 
virtual double getAbsorptionLength (const double lambda) const =0
 Absorption length. More...
 
virtual double getScatteringLength (const double lambda) const =0
 Scattering length. More...
 
virtual double getScatteringProbability (const double ct) const =0
 Model specific function to describe light scattering in water (integral over full solid angle normalised to unity). More...
 

Protected Attributes

double wavelength
 
double tmin
 
double tmax
 
JQuadrature main_engine
 
JQuadrature beta_engine
 
std::vector< JElement3D_t > phi_engine
 

Detailed Description

Probability Density Functions of the time response of a PMT.

Definition at line 61 of file JLED.hh.

Constructor & Destructor Documentation

JPHYSICS::JLED::JLED ( const double  lambda,
const double  Tmin_ns,
const double  Tmax_ns,
const JQuadrature &  engine = JCotangent(20),
const int  numberOfPoints = 20,
const double  epsilon = 1e-12 
)
inline

Constructor.

Parameters
lambdawavelength photon [nm]
Tmin_nsminimal time of emmision [ns]
Tmax_nsminimal time of emmision [ns]
enginescattering angle integrator
numberOfPointsnumber of points for integration
epsilonprecision of points for integration

Definition at line 78 of file JLED.hh.

83  :
84  wavelength(lambda),
85  tmin(Tmin_ns),
86  tmax(Tmax_ns),
87  main_engine(JGaussLegendre(numberOfPoints,epsilon)),
88  beta_engine(engine)
89  {
90  const double dp = PI / numberOfPoints;
91 
92  for (double x = 0.5*dp; x < PI; x += dp) {
93  phi_engine.push_back(JElement3D_t(cos(x),sin(x), dp));
94  }
95  }
JPHYSICS::JElement3D_t
JTOOLS::JElement3D< double, double > JElement3D_t
Definition: JLED.hh:28
JPHYSICS::JLED::tmin
double tmin
Definition: JLED.hh:263
JPHYSICS::JLED::main_engine
JQuadrature main_engine
Definition: JLED.hh:266
JPHYSICS::JLED::tmax
double tmax
Definition: JLED.hh:264
sin
then usage set_variable ACOUSTICS_WORKDIR $WORKDIR set_variable FORMULA sin([0]+2 *$PI *([1]+[2]*x)*x)" set_variable DY 1.0e-8 mkdir $WORKDIR for DETECTOR in $DETECTORS[*]
JPHYSICS::JLED::phi_engine
std::vector< JElement3D_t > phi_engine
Definition: JLED.hh:268
JPHYSICS::JLED::beta_engine
JQuadrature beta_engine
Definition: JLED.hh:267
JMATH::PI
static const double PI
Mathematical constants.
Definition: JMath/JConstants.hh:20
numberOfPoints
int numberOfPoints
Definition: JResultPDF.cc:22
GAUSS_LEGENDRE::epsilon
const double epsilon
Definition: JQuadrature.cc:21
JPHYSICS::JLED::wavelength
double wavelength
Definition: JLED.hh:261

Member Function Documentation

double JPHYSICS::JLED::getDirectLightFromLED ( const double  D_m,
const double  cd,
const double  theta,
const double  phi,
const double  t_ns 
) const
inline

Probability density function for direct light from LED.

Parameters
D_mdistance between LED and PMT [m]
cdcosine angle LED orientation and LED - PMT position
thetazenith angle orientation PMT
phiazimuth angle orientation PMT
t_nstime difference relative to direct light [ns]
Returns
dP/dt [npe/ns]

Definition at line 108 of file JLED.hh.

113  {
114  const double sd = sqrt((1.0 + cd)*(1.0 - cd));
115  const double d = D_m; // photon path [m]
116 
117  const double A = getPhotocathodeArea();
118 
119  const double px = sin(theta)*cos(phi);
120  //const double py = sin(theta)*sin(phi);
121  const double pz = cos(theta);
122 
123  const double ct = sd*px + cd*pz; // cosine angle of incidence on PMT
124 
125  const double qe = getQE(wavelength);
126  const double l_abs = getAbsorptionLength(wavelength);
127  const double ls = getScatteringLength(wavelength);
128 
129  const double U = getAngularAcceptance(ct);
130  const double V = exp(-d/l_abs) * exp(-d/ls); // absorption & scattering
131  const double W = A/(d*d); // solid angle
132 
133  return qe * getLightFromLED(cd, 0.0, t_ns) * U * V * W;
134  }
V
V(JDAQEvent-JTriggerReprocessor)*1.0/(JDAQEvent+1.0e-10)
sin
then usage set_variable ACOUSTICS_WORKDIR $WORKDIR set_variable FORMULA sin([0]+2 *$PI *([1]+[2]*x)*x)" set_variable DY 1.0e-8 mkdir $WORKDIR for DETECTOR in $DETECTORS[*]
JPHYSICS::JAbstractMedium::getScatteringLength
virtual double getScatteringLength(const double lambda) const =0
Scattering length.
JPHYSICS::JAbstractMedium::getAbsorptionLength
virtual double getAbsorptionLength(const double lambda) const =0
Absorption length.
JPHYSICS::JAbstractPMT::getPhotocathodeArea
virtual double getPhotocathodeArea() const =0
Photo-cathode area of PMT.
d
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonStart.sh:47
JPHYSICS::JAbstractPMT::getAngularAcceptance
virtual double getAngularAcceptance(const double ct) const =0
Angular acceptence of PMT.
JPHYSICS::JAbstractLED::getLightFromLED
virtual double getLightFromLED(const double ct, const double phi, const double dt) const =0
Light yield from LED (number of p.e.
JPHYSICS::JAbstractPMT::getQE
virtual double getQE(const double lambda) const =0
Quantum efficiency of PMT (incl.
A
source $JPP_DIR setenv csh $JPP_DIR &dev null eval JShellParser o a A
Definition: JShellParser.csh:15
exp
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable STRING $argv[2] set_array QUANTILES set_variable FORMULA *[0] exp(-0.5 *(x-[1])*(x-[1])/([2]*[2]))" set_variable MODULE `getModule -a $DETECTOR -L "$STRING 0"` source JAcousticsToolkit.sh typeset -A TRIPODS get_tripods $WORKDIR/tripod.txt TRIPODS XMEAN
Definition: toashort_monitor.sh:46
JPHYSICS::JLED::wavelength
double wavelength
Definition: JLED.hh:261
double JPHYSICS::JLED::getScatteredLightFromLED ( const double  D_m,
const double  cd,
const double  theta,
const double  phi,
const double  t_ns 
) const
inline

Probability density function for scattered light from LED.

Parameters
D_mdistance between LED and PMT [m]
cdcosine angle LED orientation and LED - PMT position
thetazenith angle orientation PMT
phiazimuth angle orientation PMT
t_nstime difference relative to direct light [ns]
Returns
dP/dt [npe/ns]

Definition at line 147 of file JLED.hh.

152  {
153  using namespace std;
154 
155  double value = 0;
156 
157  if (t_ns >= tmin) {
158 
159  const double ng = getIndexOfRefractionGroup(wavelength);
160 
161  const double sd = sqrt((1.0 + cd)*(1.0 - cd));
162  const double l = D_m; // distance [m]
163 
164  const double A = getPhotocathodeArea();
165  //const double sr = sqrt(A/PI) / D_m;
166  //const double cr = sqrt((1.0 + sr)*(1.0 - sr));
167 
168  const double px = sin(theta)*cos(phi);
169  const double py = sin(theta)*sin(phi);
170  const double pz = cos(theta);
171 
172  const double qx = cd*px + 0 - sd*pz;
173  const double qy = 1*py;
174  const double qz = sd*px + 0 + cd*pz;
175 
176  const double qe = getQE(wavelength);
177  const double l_abs = getAbsorptionLength(wavelength);
178  const double ls = getScatteringLength(wavelength);
179 
180  const double Jc = 1.0 / ls; // dN/dx
181 
182  const double xmin = tmin;
183  const double xmax = std::min(tmax, t_ns);
184 
185  JQuadrature buffer;
186 
187  if (xmax > xmin) {
188 
189  for (JQuadrature::const_iterator i = main_engine.begin(); i != main_engine.end(); ++i) {
190 
191  const double t = 0.5 * (xmax + xmin) + i->getX() * 0.5 * (xmax - xmin);
192  const double dt = i->getY() * 0.5 * (xmax - xmin);
193 
194  buffer[t] = dt;
195  }
196 
197  } else {
198 
199  buffer[tmin] = 1.0;
200  }
201 
202  for (JQuadrature::const_iterator i = buffer.begin(); i != buffer.end(); ++i) {
203 
204  const double t = i->getX();
205  const double dt = i->getY();
206 
207  const double d = l + C*(t_ns - t)/ng; // photon path
208  const double V = exp(-d/l_abs); // absorption
209 
210  for (JQuadrature::const_iterator j = beta_engine.begin(); j != beta_engine.end(); ++j) {
211 
212  const double cb = j->getX();
213  const double dc = j->getY();
214 
215  const double sb = sqrt((1.0 + cb)*(1.0 - cb));
216 
217  const double v = 0.5 * (d + l) * (d - l) / (d - l*cb);
218  const double u = d - v;
219 
220  if (u <= 0) continue;
221  if (v <= 0) continue;
222 
223  const double cts = (l*cb - v) / u; // cosine scattering angle
224 
225  const double W = min(A/(v*v), 2.0*PI); // solid angle
226  const double Ja = getScatteringProbability(cts); // P(cos(),phi)
227  const double Jd = ng * (1.0 - cts) / C; // dt/du
228 
229  const double ca = (l - v*cb) / u;
230  const double sa = v*sb / u;
231 
232  for (std::vector<JElement3D_t>::const_iterator k = phi_engine.begin(); k != phi_engine.end(); ++k) {
233 
234  const double cp = k->getX();
235  const double sp = k->getY();
236  const double dp = k->getZ();
237 
238  const double dom = dc * dp * v*v / (u*u);
239 
240  const double ct0 = cd*ca - sd*sa*cp;
241  const double ph0 = atan2(sa*sp, cd*sa*cp + sd*ca);
242 
243  const double vx = -sb*cp * qx;
244  const double vy = -sb*sp * qy;
245  const double vz = cb * qz;
246 
247  const double U =
248  getAngularAcceptance(vx + vy + vz) * getLightFromLED(ct0, +ph0, t) +
249  getAngularAcceptance(vx - vy + vz) * getLightFromLED(ct0, -ph0, t);
250 
251  value += qe * dt * dom * Ja * Jc * U * V * W / Jd;
252  }
253  }
254  }
255  }
256 
257  return value;
258  }
GAUSS_LEGENDRE::xmax
const double xmax
Definition: JQuadrature.cc:24
JPHYSICS::JLED::tmin
double tmin
Definition: JLED.hh:263
k
then fatal No hydrophone data file $HYDROPHONE_TXT fi sort gr k
Definition: hydrophone-phi:fit.sh:35
JPHYSICS::JLED::main_engine
JQuadrature main_engine
Definition: JLED.hh:266
JPHYSICS::JLED::tmax
double tmax
Definition: JLED.hh:264
V
V(JDAQEvent-JTriggerReprocessor)*1.0/(JDAQEvent+1.0e-10)
sin
then usage set_variable ACOUSTICS_WORKDIR $WORKDIR set_variable FORMULA sin([0]+2 *$PI *([1]+[2]*x)*x)" set_variable DY 1.0e-8 mkdir $WORKDIR for DETECTOR in $DETECTORS[*]
JPHYSICS::C
static const double C
Physics constants.
Definition: JPhysics/JConstants.hh:25
JPHYSICS::JLED::phi_engine
std::vector< JElement3D_t > phi_engine
Definition: JLED.hh:268
JPHYSICS::JAbstractMedium::getScatteringLength
virtual double getScatteringLength(const double lambda) const =0
Scattering length.
JPHYSICS::JLED::beta_engine
JQuadrature beta_engine
Definition: JLED.hh:267
JPHYSICS::JAbstractMedium::getAbsorptionLength
virtual double getAbsorptionLength(const double lambda) const =0
Absorption length.
JMATH::PI
static const double PI
Mathematical constants.
Definition: JMath/JConstants.hh:20
JPHYSICS::JAbstractPMT::getPhotocathodeArea
virtual double getPhotocathodeArea() const =0
Photo-cathode area of PMT.
GAUSS_LEGENDRE::xmin
const double xmin
Definition: JQuadrature.cc:23
d
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonStart.sh:47
JPHYSICS::JDispersionInterface::getIndexOfRefractionGroup
virtual double getIndexOfRefractionGroup(const double lambda) const
Index of refraction for group velocity.
Definition: JDispersionInterface.hh:52
JPHYSICS::JAbstractPMT::getAngularAcceptance
virtual double getAngularAcceptance(const double ct) const =0
Angular acceptence of PMT.
JPHYSICS::JAbstractLED::getLightFromLED
virtual double getLightFromLED(const double ct, const double phi, const double dt) const =0
Light yield from LED (number of p.e.
cp
then cp
Definition: darkroomCalibration.sh:41
JTOOLS::JCollection< JElement2D_t >::const_iterator
container_type::const_iterator const_iterator
Definition: JCollection.hh:91
JTOOLS::j
int j
Definition: JPolint.hh:792
JPHYSICS::JAbstractPMT::getQE
virtual double getQE(const double lambda) const =0
Quantum efficiency of PMT (incl.
JTOOLS::v
data_type v[N+1][M+1]
Definition: JPolint.hh:866
JTOOLS::u
double u[N+1]
Definition: JPolint.hh:865
A
source $JPP_DIR setenv csh $JPP_DIR &dev null eval JShellParser o a A
Definition: JShellParser.csh:15
exp
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable STRING $argv[2] set_array QUANTILES set_variable FORMULA *[0] exp(-0.5 *(x-[1])*(x-[1])/([2]*[2]))" set_variable MODULE `getModule -a $DETECTOR -L "$STRING 0"` source JAcousticsToolkit.sh typeset -A TRIPODS get_tripods $WORKDIR/tripod.txt TRIPODS XMEAN
Definition: toashort_monitor.sh:46
JPHYSICS::JAbstractMedium::getScatteringProbability
virtual double getScatteringProbability(const double ct) const =0
Model specific function to describe light scattering in water (integral over full solid angle normali...
JPHYSICS::JLED::wavelength
double wavelength
Definition: JLED.hh:261
virtual double JPHYSICS::JDispersionInterface::getIndexOfRefractionPhase ( const double  lambda) const
pure virtualinherited

Index of refraction for phase velocity.

Parameters
lambdawavelenth [nm]
Returns
index of refraction

Implemented in JPHYSICS::JDispersion.

virtual double JPHYSICS::JDispersionInterface::getDispersionPhase ( const double  lambda) const
pure virtualinherited

Dispersion of light for phase velocity.

Parameters
lambdawavelength of light [nm]
Returns
dn/dlambda

Implemented in JPHYSICS::JDispersion.

virtual double JPHYSICS::JDispersionInterface::getIndexOfRefractionGroup ( const double  lambda) const
inlinevirtualinherited

Index of refraction for group velocity.

Parameters
lambdawavelenth [nm]
Returns
index of refraction

Definition at line 52 of file JDispersionInterface.hh.

53  {
54  const double n = getIndexOfRefractionPhase(lambda);
55  const double y = getDispersionPhase(lambda);
56 
57  return n / (1.0 + y*lambda/n);
58  }
JTOOLS::n
const int n
Definition: JPolint.hh:786
JPHYSICS::JDispersionInterface::getIndexOfRefractionPhase
virtual double getIndexOfRefractionPhase(const double lambda) const =0
Index of refraction for phase velocity.
JPHYSICS::JDispersionInterface::getDispersionPhase
virtual double getDispersionPhase(const double lambda) const =0
Dispersion of light for phase velocity.
virtual double JPHYSICS::JDispersionInterface::getDispersionGroup ( const double  lambda) const
pure virtualinherited

Dispersion of light for group velocity.

Parameters
lambdawavelength of light [nm]
Returns
dn/dlambda

Implemented in JPHYSICS::JDispersion.

double JPHYSICS::JDispersionInterface::getKappa ( const double  lambda) const
inlineinherited

Get effective index of refraction for muon light.

Parameters
lambdawavelength of light [nm]
Returns
index of refraction

Definition at line 76 of file JDispersionInterface.hh.

77  {
78  const double n = getIndexOfRefractionPhase(lambda);
79  const double ng = getIndexOfRefractionGroup(lambda);
80 
81  return (ng * n - 1.0) / sqrt(n*n - 1.0);
82  }
JTOOLS::n
const int n
Definition: JPolint.hh:786
JPHYSICS::JDispersionInterface::getIndexOfRefractionPhase
virtual double getIndexOfRefractionPhase(const double lambda) const =0
Index of refraction for phase velocity.
JPHYSICS::JDispersionInterface::getIndexOfRefractionGroup
virtual double getIndexOfRefractionGroup(const double lambda) const
Index of refraction for group velocity.
Definition: JDispersionInterface.hh:52
double JPHYSICS::JDispersionInterface::getKmin ( const double  lambda) const
inlineinherited

Get smallest index of refraction for Bremsstrahlung light (i.e. point at which dt/dz = 0).

Parameters
lambdawavelength of light [nm]
Returns
index of refraction

Definition at line 91 of file JDispersionInterface.hh.

92  {
93  const double ng = getIndexOfRefractionGroup(lambda);
94 
95  return sqrt(ng*ng - 1.0);
96  }
JPHYSICS::JDispersionInterface::getIndexOfRefractionGroup
virtual double getIndexOfRefractionGroup(const double lambda) const
Index of refraction for group velocity.
Definition: JDispersionInterface.hh:52
virtual double JPHYSICS::JAbstractPMT::getPhotocathodeArea ( ) const
pure virtualinherited

Photo-cathode area of PMT.

Returns
photo-cathode area [m^2]

Implemented in JPHYSICS::JPDF_C, and JPHYSICS::JLED_C.

virtual double JPHYSICS::JAbstractPMT::getQE ( const double  lambda) const
pure virtualinherited

Quantum efficiency of PMT (incl.

absorption in glass, gel, etc.).

Parameters
lambdawavelenth [nm]
Returns
QE

Implemented in JPHYSICS::JPDF_C, and JPHYSICS::JLED_C.

virtual double JPHYSICS::JAbstractPMT::getAngularAcceptance ( const double  ct) const
pure virtualinherited

Angular acceptence of PMT.

Parameters
ctcosine angle of incidence
Returns
acceptance

Implemented in JPHYSICS::JPDF_C, and JPHYSICS::JLED_C.

virtual double JPHYSICS::JAbstractLED::getLightFromLED ( const double  ct,
const double  phi,
const double  dt 
) const
pure virtualinherited

Light yield from LED (number of p.e.

per unit solid angle per unit time).

Parameters
ctzenith angle of emission
phiazimuth angle of emission
dttime of emission [ns]
Returns
d^2P / dOmega dt [npe/ns/sr]

Implemented in JPHYSICS::JLED_C, and LED.

virtual double JPHYSICS::JAbstractMedium::getAbsorptionLength ( const double  lambda) const
pure virtualinherited

Absorption length.

Parameters
lambdawavelenth [nm]
Returns
absorption length [m]

Implemented in JPHYSICS::JPDF_C, and JPHYSICS::JLED_C.

virtual double JPHYSICS::JAbstractMedium::getScatteringLength ( const double  lambda) const
pure virtualinherited

Scattering length.

Parameters
lambdawavelenth [nm]
Returns
scattering length [m]

Implemented in JPHYSICS::JPDF_C, and JPHYSICS::JLED_C.

virtual double JPHYSICS::JAbstractMedium::getScatteringProbability ( const double  ct) const
pure virtualinherited

Model specific function to describe light scattering in water (integral over full solid angle normalised to unity).

Parameters
ctcosine scattering angle
Returns
probability

Implemented in JPHYSICS::JPDF_C, and JPHYSICS::JLED_C.

Member Data Documentation

double JPHYSICS::JLED::wavelength
protected

Definition at line 261 of file JLED.hh.

double JPHYSICS::JLED::tmin
protected

Definition at line 263 of file JLED.hh.

double JPHYSICS::JLED::tmax
protected

Definition at line 264 of file JLED.hh.

JQuadrature JPHYSICS::JLED::main_engine
protected

Definition at line 266 of file JLED.hh.

JQuadrature JPHYSICS::JLED::beta_engine
protected

Definition at line 267 of file JLED.hh.

std::vector<JElement3D_t> JPHYSICS::JLED::phi_engine
protected

Definition at line 268 of file JLED.hh.

The documentation for this class was generated from the following file:
Generated by   doxygen 1.8.5