Binned one dimensional PDF (a projection of PDF to time axis) More...

#include <PDF1Dim.h>

Collaboration diagram for PDF1Dim:

Public Member Functions
	PDF1Dim (const PDFConstructor &pdfConstructor, double binSize, double timeWindow)
	Full constructor. More...

double	getBinSize () const
	Returns actual bin size used. More...

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

TH1F	getHistogram (std::string name, std::string title) const
	Returns binned one dimensional PDF (projection to time axis) More...

const std::vector< double > &	getPhotonTimes () const
	Returns photon times of selected photons in a given slot. More...

int	getNumOfPhotons () const
	Returns number of photons. More...

double	getExpectedSignal () const
	Returns expected number of signal photons. More...

double	getExpectedDeltaPhotons () const
	Returns expected number of delta-ray photons. More...

double	getExpectedBG () const
	Returns expected number of background photons. More...

double	getTminPDF () const
	Returns minimal time of signal PDF. More...

double	getTmaxPDF () const
	Returns maximal time of signal PDF. More...

double	getTminFot () const
	Returns minimal time of selected photons. More...

double	getTmaxFot () const
	Returns maximal time of selected photons. More...

double	getMinT0 () const
	Returns lower edge of the T0 search region. More...

double	getMaxT0 () const
	Returns upper edge of the T0 search region. More...

int	getNumBinsT0 () const
	Returns number of bins the T0 search region. More...

double	getLogL (double timeShift) const
	Returns log likelihood. More...

Private Attributes
int	m_moduleID = 0
	slot number

int	m_numBins = 0
	number of bins for signal PDF

double	m_minTime = 0
	lower edge of the first bin

double	m_maxTime = 0
	upper edge of the last bin

std::vector< double >	m_times
	photon times, from a given slot

double	m_tminPDF = 0
	minimal time of signal PDF

double	m_tmaxPDF = 0
	maximal time of signal PDF

double	m_tminFot = 0
	minimal time of photons

double	m_tmaxFot = 0
	maximal time of photons

double	m_minT0 = 0
	minimal T0

double	m_maxT0 = 0
	maximal T0

int	m_numBinsT0 = 0
	number of bins for T0 finder w/ same bin size as PDF

double	m_expectedBG = 0
	expected number of background photons

double	m_expectedSignal = 0
	expected number of signal photons

double	m_expectedDelta = 0
	expected number of delta-ray photons

double	m_bkg = 0
	background [photons/bin]

double	m_binSize = 0
	bin size

std::vector< double >	m_logF
	log(PDF) values

double	m_logBkg = 0
	log(m_bkg)

Detailed Description

Binned one dimensional PDF (a projection of PDF to time axis)

Definition at line 25 of file PDF1Dim.h.

Constructor & Destructor Documentation

◆ PDF1Dim()

PDF1Dim	(	const PDFConstructor &	pdfConstructor,
		double	binSize,
		double	timeWindow
	)

Full constructor.

Parameters

pdfConstructor	reconstruction object
binSize	approximate bin size (use getBinSize() to get the actual one)
timeWindow	full time window size

Definition at line 22 of file PDF1Dim.cc.

                                                                                            :
       m_moduleID(pdfConstructor.getModuleID()), m_minTime(TOPRecoManager::getMinTime()),
       m_maxTime(TOPRecoManager::getMaxTime())
     {
  
       if (binSize <= 0) B2FATAL("TOP::PDF1Dim: bin size must be positive");
       if (not pdfConstructor.isValid()) {
         B2ERROR("TOP::PDF1Dim: PDFConstructor is not valid, cannot continue");
         return;
       }
  
       // binning
       m_numBins = (m_maxTime - m_minTime) / binSize + 1;
       m_binSize = (m_maxTime - m_minTime) / m_numBins;
  
       // photon times
       for (const auto& hit : pdfConstructor.getSelectedHits()) {
         m_times.push_back(hit.time);
       }
       if (not m_times.empty()) {
         std::sort(m_times.begin(), m_times.end());
         m_tminFot = m_times.front();
         m_tmaxFot = m_times.back();
       }
  
       m_expectedSignal = pdfConstructor.getExpectedSignalPhotons();
       m_expectedDelta = pdfConstructor.getExpectedDeltaPhotons();
       m_expectedBG = std::max(pdfConstructor.getBkgRate() * timeWindow, 0.1);
  
       // temporary histogram to make 1D projection of PDF
       TH1F pdfHisto("pdf_temporary", "", m_numBins, m_minTime, m_maxTime);
  
       // first add background and delta-ray PDF's
       m_bkg = m_expectedBG / timeWindow * m_binSize;
       const auto& deltaRayPDF = pdfConstructor.getDeltaRayPDF();
       for (int i = 0; i < m_numBins; i++) {
         double t = pdfHisto.GetBinCenter(i + 1);
         pdfHisto.SetBinContent(i + 1, m_bkg + m_expectedDelta * deltaRayPDF.getPDFValue(t) * m_binSize);
       }
  
       // then fill signal PDF
       bool start = true;
       for (const auto& pixelPDF : pdfConstructor.getSignalPDF()) {
         for (const auto& peak : pixelPDF.getPDFPeaks()) {
           pdfHisto.Fill(peak.t0, m_expectedSignal * peak.nph);
           if (start) {
             start = false;
             m_tminPDF = peak.t0;
             m_tmaxPDF = peak.t0;
           }
           m_tminPDF = std::min(m_tminPDF, peak.t0);
           m_tmaxPDF = std::max(m_tmaxPDF, peak.t0);
         }
       }
  
       // finally make look-up table for log(PDF)
       for (int i = 0; i < m_numBins; i++) {
         m_logF.push_back(log(pdfHisto.GetBinContent(i + 1)));
       }
       m_logBkg = log(m_bkg);
  
       // determine T0 search region
       m_minT0 = m_tminFot - m_tmaxPDF - 2 * m_binSize;
       double maxT0 = m_tmaxFot - m_tminPDF;
       m_numBinsT0 = std::max(int((maxT0 - m_minT0) / m_binSize), 3) + 1;
       m_maxT0 = m_minT0 + m_binSize * m_numBinsT0;
  
     }

Member Function Documentation

◆ getBinSize()

double getBinSize ( ) const

inline

Returns actual bin size used.

Returns: bin size

Definition at line 41 of file PDF1Dim.h.

41 {return m_binSize;}

◆ getExpectedBG()

double getExpectedBG ( ) const

inline

Returns expected number of background photons.

Returns: number of background photons

Definition at line 85 of file PDF1Dim.h.

◆ getExpectedDeltaPhotons()

double getExpectedDeltaPhotons ( ) const

inline

Returns expected number of delta-ray photons.

Returns: number of delta-ray photons

Definition at line 79 of file PDF1Dim.h.

◆ getExpectedSignal()

double getExpectedSignal ( ) const

inline

Returns expected number of signal photons.

Returns: number of signal photons

Definition at line 73 of file PDF1Dim.h.

◆ getHistogram()

TH1F getHistogram	(	std::string	name,
		std::string	title
	)		const

Returns binned one dimensional PDF (projection to time axis)

Parameters

name	histogram name
title	histogram title

Returns: PDF (signal + background)

Definition at line 108 of file PDF1Dim.cc.

◆ getLogL()

double getLogL ( double timeShift ) const

Returns log likelihood.

Parameters

timeShift time shift of PDF

Returns: log likelihood

Definition at line 92 of file PDF1Dim.cc.

◆ getMaxT0()

double getMaxT0 ( ) const

inline

Returns upper edge of the T0 search region.

Returns: upper edge

Definition at line 121 of file PDF1Dim.h.

◆ getMinT0()

double getMinT0 ( ) const

inline

Returns lower edge of the T0 search region.

Returns: lower edge

Definition at line 115 of file PDF1Dim.h.

◆ getModuleID()

int getModuleID ( ) const

inline

Returns slot number.

Returns: slot number

Definition at line 47 of file PDF1Dim.h.

◆ getNumBinsT0()

int getNumBinsT0 ( ) const

inline

Returns number of bins the T0 search region.

Bin size is the same as for PDF.

Returns: number of bins

Definition at line 127 of file PDF1Dim.h.

◆ getNumOfPhotons()

int getNumOfPhotons ( ) const

inline

Returns number of photons.

Returns: number of photons

Definition at line 67 of file PDF1Dim.h.

◆ getPhotonTimes()

const std::vector<double>& getPhotonTimes ( ) const

inline

Returns photon times of selected photons in a given slot.

Returns: times of TOPDigit::c_Good digits

Definition at line 61 of file PDF1Dim.h.

◆ getTmaxFot()

double getTmaxFot ( ) const

inline

Returns maximal time of selected photons.

Returns: maximal time

Definition at line 109 of file PDF1Dim.h.

◆ getTmaxPDF()

double getTmaxPDF ( ) const

inline

Returns maximal time of signal PDF.

Returns: maximal time

Definition at line 97 of file PDF1Dim.h.

◆ getTminFot()

double getTminFot ( ) const

inline

Returns minimal time of selected photons.

Returns: minimal time

Definition at line 103 of file PDF1Dim.h.

◆ getTminPDF()

double getTminPDF ( ) const

inline

Returns minimal time of signal PDF.

Returns: minimal time

Definition at line 91 of file PDF1Dim.h.

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

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

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ PDF1Dim()

Member Function Documentation

◆ getBinSize()

◆ getExpectedBG()

◆ getExpectedDeltaPhotons()

◆ getExpectedSignal()

◆ getHistogram()

◆ getLogL()

◆ getMaxT0()

◆ getMinT0()

◆ getModuleID()

◆ getNumBinsT0()

◆ getNumOfPhotons()

◆ getPhotonTimes()

◆ getTmaxFot()

◆ getTmaxPDF()

◆ getTminFot()

◆ getTminPDF()