Parametrization of delta-ray PDF in pixels of single module. More...

#include <DeltaRayPDF.h>

Collaboration diagram for DeltaRayPDF:

[legend]

Classes
struct	GausXY
	Normal (Gaussian) distribution: an entry for the table. More...

Public Member Functions
	DeltaRayPDF (int moduleID)
	Class constructor. More...

void	prepare (const TOPTrack &track, const Const::ChargedStable &hypothesis)
	Prepare the object. More...

int	getModuleID () const
	Returns slot ID. More...

double	getNumPhotons () const
	Returns number of photons. More...

double	getFraction () const
	Returns fraction of delta-ray photons in the time window. More...

double	getPDFValue (int pixelID, double time) const
	Returns PDF value at given time and pixel. More...

double	getPDFValue (double time, double dt0=0, double acc=1) const
	Returns PDF value at given time and integrated over all pixels. More...

double	getIntegral (double minTime, double maxTime) const
	Returns integral of PDF from minTime to maxTime. More...

Private Member Functions
double	angularDistr (double kz) const
	Angular distribution of photons from delta rays w/ total reflection requirement. More...

double	timeDistr (double t, double t0) const
	Time distribution of photons from delta rays (normalized). More...

double	smearedTimeDistr (double t, double t0) const
	Smeared time distribution of photons from delta rays (normalized). More...

double	peakFraction (double tmin, double tmax, double t0) const
	Fraction of delta-ray photons within given propagation time interval for single peak at t0. More...

double	totalFraction (double tmin, double tmax) const
	Total fraction of delta-ray photons within given propagation time interval. More...

double	directFraction (double z) const
	Fraction of direct photons from delta-rays, e.g direct/(direct+reflected) More...

double	photonYield (double beta, int PDGCode) const
	Photon yield from delta-rays per track length in quartz for nominal photon detection efficiency. More...

Private Attributes
int	m_moduleID
	slot ID

const BackgroundPDF *	m_background = 0
	background PDF

const PixelPositions *	m_pixelPositions = 0
	pixel positions

double	m_zD = 0
	detector (photo-cathode) position in z

double	m_zM = 0
	spherical mirror position in z

double	m_phaseIndex = 0
	phase refractive index

double	m_groupIndex = 0
	group refractive index

double	m_dispersion = 0
	dispersion coefficient

double	m_angularNorm = 0
	angular distribution normalization constant

std::vector< double >	m_norms
	relative angular distribution normalization constants (cumulative)

std::vector< GausXY >	m_tableGaus
	table of normal (Gaussian) distribution

double	m_xE = 0
	average photon emission position in x

double	m_yE = 0
	average photon emission position in y

double	m_zE = 0
	average photon emission position in z

double	m_dirFrac = 0
	fraction of direct photons

double	m_dirT0 = 0
	minimal propagation time of direct photons

double	m_reflT0 = 0
	minimal propagation time of reflected photons

double	m_TOF = 0
	time-of-flight of particle

double	m_fraction = 0
	fraction of delta-ray photons within time window

double	m_numPhotons = 0
	number of photons

std::vector< double >	m_pixelAcceptances
	pixel angular acceptances for direct peak (index = pixelID - 1)

Detailed Description

Parametrization of delta-ray PDF in pixels of single module.

Definition at line 27 of file DeltaRayPDF.h.

Constructor & Destructor Documentation

◆ DeltaRayPDF()

DeltaRayPDF ( int moduleID )

explicit

Class constructor.

Parameters

moduleID slot ID

Definition at line 24 of file DeltaRayPDF.cc.

                                         :
       m_moduleID(moduleID), m_background(TOPRecoManager::getBackgroundPDF(moduleID))
     {
       if (not m_background) {
         B2ERROR("TOP::DeltaRayPDF: background PDF not found");
         return;
       }
  
       const auto* yScanner = TOPRecoManager::getYScanner(moduleID);
       if (not yScanner) {
         B2ERROR("TOP::DeltaRayPDF: YScanner not found");
         return;
       }
  
       m_pixelPositions = &(yScanner->getPixelPositions());
  
       m_zD = yScanner->getPrism().zD;
       m_zM = yScanner->getBars().back().zR;
       double meanE = yScanner->getMeanEnergyBeta1();
       double sigE = yScanner->getRMSEnergyBeta1();
       m_phaseIndex = TOPGeometryPar::Instance()->getPhaseIndex(meanE);
       m_groupIndex = TOPGeometryPar::Instance()->getGroupIndex(meanE);
       double dng_de = TOPGeometryPar::Instance()->getGroupIndexDerivative(meanE);
       m_dispersion = dng_de / m_groupIndex * sigE;
  
       const int N = 1000;
       double dkz = 1.0 / N;
       m_norms.push_back(0);
       for (int i = 0; i < N; i++) {
         double kz = (i + 0.5) * dkz;
         m_angularNorm += angularDistr(kz) * dkz;
         m_norms.push_back(m_angularNorm);
       }
       for (auto& norm : m_norms) norm /= m_angularNorm;
  
       const int Np = 101;
       double dx = 6.0 / (Np - 1);
       for (int i = 0; i < Np; i++) {
         double x = i * dx - 3.0;
         m_tableGaus.push_back(GausXY(x));
       }
       double sum = 0;
       for (const auto& gaus : m_tableGaus) sum += gaus.y;
       for (auto& gaus : m_tableGaus) gaus.y /= sum;
     }

Member Function Documentation

◆ angularDistr()

double angularDistr ( double kz ) const

inlineprivate

Angular distribution of photons from delta rays w/ total reflection requirement.

Distribution is not normalized.

Parameters

kz	z-component of photon direction

Returns: distribution value

Definition at line 190 of file DeltaRayPDF.h.

     {
       if (std::abs(kz) < 1 / m_phaseIndex) {
         double x = sqrt((1 - 1 / pow(m_phaseIndex, 2)) / (1 - pow(kz, 2)));
         return 1 - acos(x) * 4 / M_PI;
       }
       return 1;
     }

◆ directFraction()

double directFraction ( double z ) const

private

Fraction of direct photons from delta-rays, e.g direct/(direct+reflected)

Parameters

z	local z position of track at TOP

Returns: fraction of direct photons

Definition at line 155 of file DeltaRayPDF.cc.

◆ getFraction()

double getFraction ( ) const

inline

Returns fraction of delta-ray photons in the time window.

Returns: fraction of delta-ray photons in the time window

Definition at line 60 of file DeltaRayPDF.h.

◆ getIntegral()

double getIntegral	(	double	minTime,
		double	maxTime
	)		const

inline

Returns integral of PDF from minTime to maxTime.

Parameters

minTime	integral lower limit
maxTime	integral upper limit

Returns: integral of PDF

Definition at line 213 of file DeltaRayPDF.h.

◆ getModuleID()

int getModuleID ( ) const

inline

Returns slot ID.

Returns: slot ID

Definition at line 48 of file DeltaRayPDF.h.

◆ getNumPhotons()

double getNumPhotons ( ) const

inline

Returns number of photons.

Returns: number of photons

Definition at line 54 of file DeltaRayPDF.h.

◆ getPDFValue() [1/2]

double getPDFValue	(	double	time,
		double	dt0 = `0`,
		double	acc = `1`
	)		const

inline

Returns PDF value at given time and integrated over all pixels.

Parameters

time	photon hit time
dt0	direct peak position correction
acc	acceptance correction factor for direct peak

Returns: PDF value (projection to time axis)

Definition at line 205 of file DeltaRayPDF.h.

◆ getPDFValue() [2/2]

double getPDFValue	(	int	pixelID,
		double	time
	)		const

Returns PDF value at given time and pixel.

Parameters

pixelID	pixel ID
time	photon hit time

Returns: PDF value

Definition at line 113 of file DeltaRayPDF.cc.

◆ peakFraction()

double peakFraction	(	double	tmin,
		double	tmax,
		double	t0
	)		const

private

Fraction of delta-ray photons within given propagation time interval for single peak at t0.

Parameters

tmin	time interval lower edge
tmax	time interval upper edge
t0	minimal possible propagation time

Returns: fraction of photons within time interval for single peak

Definition at line 138 of file DeltaRayPDF.cc.

◆ photonYield()

double photonYield	(	double	beta,
		int	PDGCode
	)		const

private

Photon yield from delta-rays per track length in quartz for nominal photon detection efficiency.

Parameters

beta	particle velocity
PDGCode	PDG code

Returns: photon yield per cm

Definition at line 169 of file DeltaRayPDF.cc.

◆ prepare()

void prepare	(	const TOPTrack &	track,
		const Const::ChargedStable &	hypothesis
	)

Prepare the object.

Parameters

track	track at TOP
hypothesis	particle hypothesis

Definition at line 70 of file DeltaRayPDF.cc.

◆ smearedTimeDistr()

double smearedTimeDistr	(	double	t,
		double	t0
	)		const

private

Smeared time distribution of photons from delta rays (normalized).

Parameters

t	photon propagation time
t0	minimal possible propagation time

Returns: normalized time distribution value at given propagaton time

Definition at line 126 of file DeltaRayPDF.cc.

◆ timeDistr()

double timeDistr	(	double	t,
		double	t0
	)		const

inlineprivate

Time distribution of photons from delta rays (normalized).

Parameters

t	photon propagation time
t0	minimal possible propagation time

Returns: normalized time distribution value at given propagaton time

Definition at line 199 of file DeltaRayPDF.h.

◆ totalFraction()

double totalFraction	(	double	tmin,
		double	tmax
	)		const

private

Total fraction of delta-ray photons within given propagation time interval.

Parameters

tmin	time interval lower edge
tmax	time interval upper edge

Returns: fraction of photons within time interval including direct and reflected peaks

Definition at line 148 of file DeltaRayPDF.cc.

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

top/reconstruction_cpp/include/DeltaRayPDF.h
top/reconstruction_cpp/src/DeltaRayPDF.cc

Classes

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ DeltaRayPDF()

Member Function Documentation

◆ angularDistr()

◆ directFraction()

◆ getFraction()

◆ getIntegral()

◆ getModuleID()

◆ getNumPhotons()

◆ getPDFValue() [1/2]

◆ getPDFValue() [2/2]

◆ peakFraction()

◆ photonYield()

◆ prepare()

◆ smearedTimeDistr()

◆ timeDistr()

◆ totalFraction()