SignalPDF Class Reference

Parametrization of signal PDF in a single pixel. More...

#include <SignalPDF.h>

Classes
struct	PDFExtra
	Extra information about single PDF peak. More...
struct	PDFPeak
	Single PDF peak. More...

Public Types
enum	EPeakType {   c_Unknown = 0 ,   c_Direct = 1 ,   c_Reflected = 2 }
	Enumerator for single PDF peak types. More...

Public Member Functions
	SignalPDF (int pixelID, const TOPNominalTTS &tts)
	Class constructor.
void	append (const PDFPeak &peak)
	Appends single PDF peak.
void	append (const PDFExtra &extra)
	Appends extra information of single PDF peak.
void	normalize (double a)
	Normalize PDF peaks by dividing nph with given constant.
double	getSum () const
	Returns a sum of nph of all peaks.
int	getPixelID () const
	Returns pixel ID.
const TOPNominalTTS *	getTTS () const
	Returns TTS.
const std::vector< PDFPeak > &	getPDFPeaks () const
	Returns PDF peaks.
const std::vector< PDFExtra > &	getPDFExtraInfo () const
	Returns PDF extra info Note: available only if stored on demand.
bool	isExtraInfoAvailable () const
	Checks if extra info is available.
double	getPDFValue (double time, double timeErr, double sigt=0) const
	Returns PDF value at given time.
double	getIntegral (double minTime, double maxTime) const
	Returns integral of PDF from minTime to maxTime.

Private Attributes
int	m_pixelID = 0
	pixel ID
const TOPNominalTTS *	m_tts = 0
	nominal TTS
std::vector< PDFPeak >	m_peaks
	PDF peaks.
std::vector< PDFExtra >	m_extras
	extra information (stored only on request)

Detailed Description

Parametrization of signal PDF in a single pixel.

PDF is parametrized as a sum of Gaussian distributions.

Definition at line 25 of file SignalPDF.h.

Member Enumeration Documentation

EPeakType

enum EPeakType

Enumerator for single PDF peak types.

Enumerator
c_Unknown	unknown
c_Direct	direct photon
c_Reflected	reflected photon

Definition at line 32 of file SignalPDF.h.

                     {
        c_Unknown = 0, 
        c_Direct = 1,   
        c_Reflected = 2 
      };

Constructor & Destructor Documentation

SignalPDF()

SignalPDF	(	int	pixelID,
		const TOPNominalTTS &	tts
	)

Class constructor.

Parameters

pixelID	pixel ID
tts	TTS

Definition at line 21 of file SignalPDF.cc.

                                                             :
      m_pixelID(pixelID), m_tts(&tts)
    {}

Member Function Documentation

append() [1/2]

void append ( const PDFExtra &  extra )

inline

Appends extra information of single PDF peak.

Parameters

extra

extra information

Definition at line 91 of file SignalPDF.h.

{ m_extras.push_back(extra);}

append() [2/2]

void append ( const PDFPeak &  peak )

inline

Appends single PDF peak.

Parameters

peak	single PDF peak

Definition at line 85 of file SignalPDF.h.

{ m_peaks.push_back(peak);}

getIntegral()

double getIntegral	(	double	minTime,
		double	maxTime
	)		const

Returns integral of PDF from minTime to maxTime.

Parameters

minTime	integral lower limit
maxTime	integral upper limit

Returns

integral of PDF

Definition at line 41 of file SignalPDF.cc.

    {
      double s = 0;
      for (const auto& peak : m_peaks) {
        if (peak.t0 < minTime or peak.t0 > maxTime) continue;
        s += peak.nph;
      }
      return s;
    }

getPDFExtraInfo()

const std::vector< PDFExtra > & getPDFExtraInfo ( ) const

inline

Returns PDF extra info Note: available only if stored on demand.

Returns

PDF extra info

Definition at line 127 of file SignalPDF.h.

{return m_extras;}

getPDFPeaks()

const std::vector< PDFPeak > & getPDFPeaks ( ) const

inline

Returns PDF peaks.

Returns

PDF peaks

Definition at line 120 of file SignalPDF.h.

{return m_peaks;}

getPDFValue()

double getPDFValue	(	double	time,
		double	timeErr,
		double	sigt = 0
	)		const

Returns PDF value at given time.

Parameters

time	photon hit time
timeErr	uncertainty of hit time (sigma)
sigt	additional time smearing

Returns

PDF value

Definition at line 25 of file SignalPDF.cc.

    {
      double wid0 = timeErr * timeErr + sigt * sigt;
      double f = 0;
      for (const auto& peak : m_peaks) {
        for (const auto& gaus : m_tts->getTTS()) {
          double sigma = peak.wid + gaus.sigma * gaus.sigma + wid0; // sigma squared!
          double x = pow(time - peak.t0 - gaus.position, 2) / sigma;
          if (x > 100) continue;
          sigma = sqrt(2 * M_PI * sigma);
          f += peak.nph * gaus.fraction / sigma * exp(-x / 2);
        }
      }
      return f;
    }

getPixelID()

int getPixelID ( ) const

inline

Returns pixel ID.

Returns

pixel ID

Definition at line 108 of file SignalPDF.h.

{return m_pixelID;}

getSum()

double getSum ( ) const

inline

Returns a sum of nph of all peaks.

Definition at line 168 of file SignalPDF.h.

    {
      double s = 0;
      for (const auto& peak : m_peaks) s += peak.nph;
      return s;
    }

getTTS()

const TOPNominalTTS * getTTS ( ) const

inline

Returns TTS.

Returns

TTS

Definition at line 114 of file SignalPDF.h.

{return m_tts;}

isExtraInfoAvailable()

bool isExtraInfoAvailable ( ) const

inline

Checks if extra info is available.

Returns

true if available

Definition at line 133 of file SignalPDF.h.

{return m_extras.size() == m_peaks.size();}

normalize()

void normalize ( double  a )

inline

Normalize PDF peaks by dividing nph with given constant.

Parameters

a	normalization constant

Definition at line 163 of file SignalPDF.h.

    {
      for (auto& peak : m_peaks) peak.nph /= a;
    }

Member Data Documentation

m_extras

std::vector<PDFExtra> m_extras

private

extra information (stored only on request)

Definition at line 157 of file SignalPDF.h.

m_peaks

std::vector<PDFPeak> m_peaks

private

PDF peaks.

Definition at line 156 of file SignalPDF.h.

m_pixelID

int m_pixelID = 0

private

pixel ID

Definition at line 154 of file SignalPDF.h.

m_tts

const TOPNominalTTS* m_tts = 0

private

nominal TTS

Definition at line 155 of file SignalPDF.h.

---

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

• top/reconstruction_cpp/include/SignalPDF.h
• top/reconstruction_cpp/src/SignalPDF.cc