Probability Density Functions of the time response of a PMT (C-like interface) More...

#include <JLED.hh>

Inheritance diagram for JPHYSICS::JLED_C:

Public Member Functions
	JLED_C (const double Area, const JAbstractLED LED, double(QE)(const double), double(Pmt)(const double), const double L_abs, const double L_s, double(Km3)(const double), const double P_atm, const double wavelength, const double Tmin_ns, const double Tmax_ns, const JQuadrature &engine=JCotangent(20), const int numberOfPoints=20, const double epsilon=1e-12)
	Constructor. More...

virtual double	getLightFromLED (const double ct, const double phi, const double dt) const
	Light yield from LED (number of p.e. More...

virtual double	getPhotocathodeArea () const
	Photo-cathode area of PMT. More...

virtual double	getQE (const double lambda) const
	Quantum efficiency of PMT (incl. More...

virtual double	getAngularAcceptance (const double ct) const
	Angular acceptence of PMT (normalised to one at cos() = -1). More...

virtual double	getAbsorptionLength (const double lambda) const
	Absorption length. More...

virtual double	getScatteringLength (const double lambda) const
	Scattering length. More...

virtual double	getScatteringProbability (const double ct) const
	Model specific function to describe light scattering in water (integral over full solid angle normalised to one). 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. More...

virtual double	getIndexOfRefractionPhase (const double lambda) const
	Index of refraction (phase velocity). More...

virtual double	getDispersionPhase (const double lambda) const
	Dispersion of light for phase velocity. More...

virtual double	getDispersionGroup (const double lambda) const
	Dispersion of light for group velocity. More...

Public Attributes
const double	P
	Dispersion parameters (x = 1/lambda) More...

const double	a0
	offset More...

const double	a1
	dn/dP More...

const double	a2
	d^1n/(dx)^1 More...

const double	a3
	d^2n/(dx)^2 More...

const double	a4
	d^3n/(dx)^3 More...

Protected Attributes
const double	A
	photo-cathode area [m2] More...

const JAbstractLED *	led
	Pointer to interface for emission profile from LED. More...

double(*	qe )(const double lambda)
	Quantum efficiency of PMT (incl. More...

double	l_abs
	Absorption length. More...

double	ls
	Scattering length. More...

double(*	pmt )(const double ct)
	Angular acceptence of PMT (normalised to one at cos() = -1) More...

double(*	km3 )(const double ct)
	Model specific function to describe light scattering in water. More...

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 (C-like interface)

Definition at line 278 of file JLED.hh.

Constructor & Destructor Documentation

JPHYSICS::JLED_C::JLED_C	(	const double	Area,
		const JAbstractLED *	LED,
		double(*)(const double)	QE,
		double(*)(const double)	Pmt,
		const double	L_abs,
		const double	L_s,
		double(*)(const double)	Km3,
		const double	P_atm,
		const double	wavelength,
		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

Area	photo-cathode area [m^2]
LED	pointer to interface for emission profile from LED
QE	pointer to function for quantum efficiency of PMT
Pmt	pointer to function for angular acceptence of PMT
L_abs	absorption length [m]
L_s	scattering length [m]
Km3	pointer to model specific function to describe light scattering in water
P_atm	ambient pressure [atm]
wavelength	wavelength of light [nm]
Tmin_ns	minimal time of emmision [ns]
Tmax_ns	minimal time of emmision [ns]
engine	scattering angle integrator
numberOfPoints	number of points for integration
epsilon	precision of points for integration

Definition at line 301 of file JLED.hh.

                                                 :
       JLED(wavelength, Tmin_ns, Tmax_ns, engine, numberOfPoints, epsilon),
       JDispersion(P_atm),
       A    (Area ),
       led  (LED  ),
       qe   (QE   ),
       l_abs(L_abs),
       ls   (L_s  ),
       pmt  (Pmt  ),
       km3  (Km3  )
     {}

Member Function Documentation

virtual double JPHYSICS::JLED_C::getLightFromLED	(	const double	ct,
		const double	phi,
		const double	dt
	)		const

inlinevirtual

Light yield from LED (number of p.e.

per unit solid angle per unit time).

Parameters

ct	zenith angle of emission
phi	azimuth angle of emission
dt	time of emission [ns]

Returns: d^2P / dOmega dt [npe/ns/sr]

Implements JPHYSICS::JAbstractLED.

Definition at line 341 of file JLED.hh.

     { 
       return led->getLightFromLED(ct, phi, dt);
     }

virtual double JPHYSICS::JLED_C::getPhotocathodeArea ( ) const

inlinevirtual

Photo-cathode area of PMT.

Returns: photo-cathode area [m^2]

Implements JPHYSICS::JAbstractPMT.

Definition at line 354 of file JLED.hh.

     { 
       return A;
     }

virtual double JPHYSICS::JLED_C::getQE ( const double lambda ) const

inlinevirtual

Quantum efficiency of PMT (incl.

absorption in glass, gel, etc.).

Parameters

lambda wavelenth [nm]

Returns: QE

Implements JPHYSICS::JAbstractPMT.

Definition at line 366 of file JLED.hh.

     { 
       return qe(lambda);
     }

virtual double JPHYSICS::JLED_C::getAngularAcceptance ( const double ct ) const

inlinevirtual

Angular acceptence of PMT (normalised to one at cos() = -1).

Parameters

ct	cosine angle of incidence

Returns: relative efficiency

Implements JPHYSICS::JAbstractPMT.

Definition at line 378 of file JLED.hh.

     { 
       return pmt(ct);
     }

virtual double JPHYSICS::JLED_C::getAbsorptionLength ( const double lambda ) const

inlinevirtual

Absorption length.

Parameters

lambda wavelenth [nm]

Returns: absorption length [m]

Implements JPHYSICS::JAbstractMedium.

Definition at line 390 of file JLED.hh.

     { 
       return l_abs;
     }

virtual double JPHYSICS::JLED_C::getScatteringLength ( const double lambda ) const

inlinevirtual

Scattering length.

Parameters

lambda wavelenth [nm]

Returns: scattering length [m]

Implements JPHYSICS::JAbstractMedium.

Definition at line 402 of file JLED.hh.

     { 
       return ls;
     }

virtual double JPHYSICS::JLED_C::getScatteringProbability ( const double ct ) const

inlinevirtual

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

Parameters

ct	cosine scattering angle

Returns: probability

Implements JPHYSICS::JAbstractMedium.

Definition at line 415 of file JLED.hh.

     { 
       return km3(ct);
     }

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

inlineinherited

Probability density function for direct light from LED.

Parameters

D_m	distance between LED and PMT [m]
cd	cosine angle LED orientation and LED - PMT position
theta	zenith angle orientation PMT
phi	azimuth angle orientation PMT
t_ns	time difference relative to direct light [ns]

Returns: dP/dt [npe/ns]

Definition at line 109 of file JLED.hh.

     {
       const double sd =  sqrt((1.0 + cd)*(1.0 - cd));
       const double d  =  D_m;                                  // photon path [m]
 
       const double A  =  getPhotocathodeArea();
 
       const double px =  sin(theta)*cos(phi);
       //const double py =  sin(theta)*sin(phi);
       const double pz =  cos(theta);
 
       const double ct = sd*px + cd*pz;                         // cosine angle of incidence on PMT
 
       const double qe    = getQE(wavelength);
       const double l_abs = getAbsorptionLength(wavelength);
       const double ls    = getScatteringLength(wavelength);
 
       const double U  = getAngularAcceptance(ct);
       const double V  = exp(-d/l_abs) * exp(-d/ls);            // absorption & scattering
       const double W  = A/(d*d);                               // solid angle
   
       return qe * getLightFromLED(cd, 0.0, t_ns) * U * V * W;
     }

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

inlineinherited

Probability density function for scattered light from LED.

Parameters

D_m	distance between LED and PMT [m]
cd	cosine angle LED orientation and LED - PMT position
theta	zenith angle orientation PMT
phi	azimuth angle orientation PMT
t_ns	time difference relative to direct light [ns]

Returns: dP/dt [npe/ns]

Definition at line 148 of file JLED.hh.

     {
       using namespace std;
       using JTOOLS::PI;
       using JTOOLS::C;
 
       double value = 0;
 
       if (t_ns >= tmin) {
 
         const double ng =  getIndexOfRefractionGroup(wavelength);
 
         const double sd =  sqrt((1.0 + cd)*(1.0 - cd));
         const double l  =  D_m;                                    // distance [m]
 
         const double A  =  getPhotocathodeArea();
         //const double sr =  sqrt(A/PI) / D_m;
         //const double cr =  sqrt((1.0 + sr)*(1.0 - sr));
 
         const double px =  sin(theta)*cos(phi);
         const double py =  sin(theta)*sin(phi);
         const double pz =  cos(theta);
 
         const double qx = cd*px +  0  - sd*pz;
         const double qy =         1*py;
         const double qz = sd*px +  0  + cd*pz;
 
         const double qe    = getQE(wavelength);
         const double l_abs = getAbsorptionLength(wavelength);
         const double ls    = getScatteringLength(wavelength);
 
         const double Jc = 1.0 / ls;                                // dN/dx
 
         const double xmin = tmin;
         const double xmax = std::min(tmax, t_ns);
 
         JQuadrature buffer;
         
         if (xmax > xmin) {
 
           for (JQuadrature::const_iterator i = main_engine.begin(); i != main_engine.end(); ++i) {
 
             const double t  = 0.5 * (xmax + xmin)  +  i->getX() * 0.5 * (xmax - xmin);
             const double dt = i->getY() * 0.5 * (xmax - xmin);
 
             buffer[t] = dt;
           }
 
         } else {
 
           buffer[tmin] = 1.0;
         }
 
         for (JQuadrature::const_iterator i = buffer.begin(); i != buffer.end(); ++i) {
 
           const double t  = i->getX();
           const double dt = i->getY();
 
           const double d  = l + C*(t_ns - t)/ng;                   // photon path
           const double V  = exp(-d/l_abs);                         // absorption
 
           for (JQuadrature::const_iterator j = beta_engine.begin(); j != beta_engine.end(); ++j) {
 
             const double cb = j->getX();
             const double dc = j->getY();
             
             const double sb = sqrt((1.0 + cb)*(1.0 - cb));
             
             const double v  = 0.5 * (d + l) * (d - l) / (d - l*cb);
             const double u  = d - v;
             
             if (u <= 0) continue;
             if (v <= 0) continue;
             
             const double cts = (l*cb - v) / u;                   // cosine scattering angle
             
             const double W  = min(A/(v*v), 2.0*PI);              // solid angle
             const double Ja = getScatteringProbability(cts);     // P(cos(),phi)
             const double Jd = ng * (1.0 - cts) / C;              // dt/du
             
             const double ca = (l - v*cb) / u;
             const double sa =      v*sb  / u;
             
             for (std::vector<JElement3D_t>::const_iterator k = phi_engine.begin(); k != phi_engine.end(); ++k) {
               
               const double cp = k->getX();
               const double sp = k->getY();
               const double dp = k->getZ();
 
               const double dom = dc * dp * v*v / (u*u);
               
               const double ct0 = cd*ca - sd*sa*cp;
               const double ph0 = atan2(sa*sp, cd*sa*cp + sd*ca);
               
               const double vx  = -sb*cp * qx;
               const double vy  = -sb*sp * qy;
               const double vz  =  cb    * qz;
               
               const double U  =
                 getAngularAcceptance(vx + vy + vz) * getLightFromLED(ct0, +ph0, t) +
                 getAngularAcceptance(vx - vy + vz) * getLightFromLED(ct0, -ph0, t);
               
               value += qe * dt * dom * Ja * Jc * U * V * W / Jd;
             }
           }
         }
       }
 
       return value; 
     }                            

virtual double JPHYSICS::JDispersionInterface::getIndexOfRefractionPhase ( const double lambda ) const

pure virtualinherited

Index of refraction for phase velocity.

Parameters

lambda wavelenth [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

lambda wavelength 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

lambda wavelenth [nm]

Returns: index of refraction

Definition at line 52 of file JDispersionInterface.hh.

     {
       const double n = getIndexOfRefractionPhase(lambda);
       const double y = getDispersionPhase(lambda);
       
       return n / (1.0 + y*lambda/n);
     }

virtual double JPHYSICS::JDispersionInterface::getDispersionGroup ( const double lambda ) const

pure virtualinherited

Dispersion of light for group velocity.

Parameters

lambda wavelength 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

lambda wavelength of light [nm]

Returns: index of refraction

Definition at line 76 of file JDispersionInterface.hh.

     {
       const double n  = getIndexOfRefractionPhase(lambda);
       const double ng = getIndexOfRefractionGroup(lambda);
 
       return (ng * n - 1.0) / sqrt(n*n - 1.0);
     }

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

lambda wavelength of light [nm]

Returns: index of refraction

Definition at line 91 of file JDispersionInterface.hh.

     {
       const double ng = getIndexOfRefractionGroup(lambda);
 
       return sqrt(ng*ng - 1.0);
     }

virtual double JPHYSICS::JDispersion::getIndexOfRefractionPhase ( const double lambda ) const

inlinevirtualinherited

Index of refraction (phase velocity).

Parameters

lambda wavelenth [nm]

Returns: index of refraction

Implements JPHYSICS::JDispersionInterface.

Definition at line 51 of file JDispersion.hh.

     {
       const double x = 1.0 / lambda;
 
       return a0  +  a1*P  +  x*(a2 + x*(a3 + x*a4));
     }

virtual double JPHYSICS::JDispersion::getDispersionPhase ( const double lambda ) const

inlinevirtualinherited

Dispersion of light for phase velocity.

Parameters

lambda wavelength of light [nm]

Returns: dn/dlambda

Implements JPHYSICS::JDispersionInterface.

Definition at line 65 of file JDispersion.hh.

     {
       const double x = 1.0 / lambda;
       
       return -x*x*(a2 + x*(2.0*a3 + x*3.0*a4));
     }

virtual double JPHYSICS::JDispersion::getDispersionGroup ( const double lambda ) const

inlinevirtualinherited

Dispersion of light for group velocity.

Parameters

lambda wavelength of light [nm]

Returns: dn/dlambda

Implements JPHYSICS::JDispersionInterface.

Definition at line 79 of file JDispersion.hh.

     {
       const double x = 1.0 / lambda;
 
       const double n   = getIndexOfRefractionPhase(lambda);
       const double np  = getDispersionPhase(lambda);
       const double npp = x*x*x*(2.0*a2 + x*(6.0*a3 + x*12.0*a4));
       const double ng  = n / (1.0 + np*lambda/n);
       
       return ng*ng * (2*np*np - n*npp) * lambda / (n*n*n);
     }

Member Data Documentation

const double JPHYSICS::JLED_C::A

protected

photo-cathode area [m2]

Definition at line 424 of file JLED.hh.

const JAbstractLED* JPHYSICS::JLED_C::led

protected

Pointer to interface for emission profile from LED.

Definition at line 429 of file JLED.hh.

double(* JPHYSICS::JLED_C::qe)(const double lambda)

protected

Quantum efficiency of PMT (incl.

absorption in glass, gel, etc.)

Parameters

lambda wavelenth [nm]

Returns: QE

Definition at line 437 of file JLED.hh.

double JPHYSICS::JLED_C::l_abs

protected

Absorption length.

Definition at line 442 of file JLED.hh.

double JPHYSICS::JLED_C::ls

protected

Scattering length.

Definition at line 447 of file JLED.hh.

double(* JPHYSICS::JLED_C::pmt)(const double ct)

protected

Angular acceptence of PMT (normalised to one at cos() = -1)

Parameters

ct	cosine angle of incidence

Returns: relative efficiency

Definition at line 455 of file JLED.hh.

double(* JPHYSICS::JLED_C::km3)(const double ct)

protected

Model specific function to describe light scattering in water.

Parameters

ct	cosine scattering angle

Returns: probability

Definition at line 463 of file JLED.hh.

double JPHYSICS::JLED::wavelength

protectedinherited

Definition at line 264 of file JLED.hh.

double JPHYSICS::JLED::tmin

protectedinherited

Definition at line 266 of file JLED.hh.

double JPHYSICS::JLED::tmax

protectedinherited

Definition at line 267 of file JLED.hh.

JQuadrature JPHYSICS::JLED::main_engine

protectedinherited

Definition at line 269 of file JLED.hh.

JQuadrature JPHYSICS::JLED::beta_engine

protectedinherited

Definition at line 270 of file JLED.hh.

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

protectedinherited

Definition at line 271 of file JLED.hh.

const double JPHYSICS::JDispersion::P

inherited

Dispersion parameters (x = 1/lambda)

ambient pressure [atm]

Definition at line 95 of file JDispersion.hh.

const double JPHYSICS::JDispersion::a0

inherited

offset

Definition at line 96 of file JDispersion.hh.

const double JPHYSICS::JDispersion::a1

inherited

dn/dP

Definition at line 97 of file JDispersion.hh.

const double JPHYSICS::JDispersion::a2

inherited

d^1n/(dx)^1

Definition at line 98 of file JDispersion.hh.

const double JPHYSICS::JDispersion::a3

inherited

d^2n/(dx)^2

Definition at line 99 of file JDispersion.hh.

const double JPHYSICS::JDispersion::a4

inherited

d^3n/(dx)^3

Definition at line 100 of file JDispersion.hh.

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

JLED.hh

Public Member Functions

Public Attributes

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation