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

#include <TOP1Dpdf.h>

Collaboration diagram for TOP1Dpdf:

Public Member Functions
	TOP1Dpdf (TOPreco &reco, const StoreArray< TOPDigit > &digits, int moduleID, double binSize)
	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	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_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 39 of file TOP1Dpdf.h.

Constructor & Destructor Documentation

◆ TOP1Dpdf()

TOP1Dpdf	(	TOPreco &	reco,
		const StoreArray< TOPDigit > &	digits,
		int	moduleID,
		double	binSize
	)

Full constructor.

Parameters

reco	reconstruction object
moduleID	slot number
binSize	approximate bin size (use getBinSize() to get the actual one)
digits	collection of TOP digits

Definition at line 30 of file TOP1Dpdf.cc.

                              : digits) {
         if (digit.getHitQuality() == TOPDigit::c_Good) {
           double t = digit.getTime();
           m_tminFot = std::min(m_tminFot, t);
           m_tmaxFot = std::max(m_tmaxFot, t);
           if (digit.getModuleID() == m_moduleID) {
             m_times.push_back(t);
           } else {
             m_expectedBG += 1;
           }
         }
       }
       m_expectedBG /= 15; // per module
       int nbins = (m_tmaxFot - m_tminFot) / m_binSize;
       m_bkg = std::max(m_expectedBG, 1.0) / std::max(nbins, m_numBins);
  
       // temporary histogram to make 1D projection of PDF
       TH1F pdfHisto("pdf_temporary", "", m_numBins, m_minTime, m_maxTime);
  
       // first add background
       for (int i = 0; i < m_numBins; i++) pdfHisto.SetBinContent(i + 1, m_bkg);
  
       // then fill signal PDF
       bool start = true;
       for (int pix = 1; pix <= 512; pix++) {
         for (int k = 0; k < reco.getNumofPDFPeaks(pix); k++) {
           float pos = 0;
           float wid = 0;
           float nph = 0;
           reco.getPDFPeak(pix, k, pos, wid, nph);
           m_expectedSignal += nph;
           double time = pos;
           pdfHisto.Fill(time, nph);
           if (start) {
             start = false;
             m_tminPDF = time;
             m_tmaxPDF = time;
           }
           m_tminPDF = std::min(m_tminPDF, time);
           m_tmaxPDF = std::max(m_tmaxPDF, time);
         }
       }
  
       // 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;
  
     }
  
  
     double TOP1Dpdf::getLogL(double timeShift) const
     {
  
       double logL = 0;
       for (const auto& time : m_times) {
         double f = m_logBkg;

Member Function Documentation

◆ getBinSize()

double getBinSize ( ) const

inline

Returns actual bin size used.

Returns: bin size

Definition at line 57 of file TOP1Dpdf.h.

◆ getExpectedBG()

double getExpectedBG ( ) const

inline

Returns expected number of background photons.

Returns: number of background photons

Definition at line 95 of file TOP1Dpdf.h.

◆ getExpectedSignal()

double getExpectedSignal ( ) const

inline

Returns expected number of signal photons.

Returns: number of signal photons

Definition at line 89 of file TOP1Dpdf.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 119 of file TOP1Dpdf.cc.

◆ getLogL()

double getLogL ( double timeShift ) const

Returns log likelihood.

Parameters

timeShift time shift of PDF

Returns: log likelihood

Definition at line 102 of file TOP1Dpdf.cc.

◆ getMaxT0()

double getMaxT0 ( ) const

inline

Returns upper edge of the T0 search region.

Returns: upper edge

Definition at line 131 of file TOP1Dpdf.h.

◆ getMinT0()

double getMinT0 ( ) const

inline

Returns lower edge of the T0 search region.

Returns: lower edge

Definition at line 125 of file TOP1Dpdf.h.

◆ getModuleID()

int getModuleID ( ) const

inline

Returns slot number.

Returns: slot number

Definition at line 63 of file TOP1Dpdf.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 137 of file TOP1Dpdf.h.

◆ getNumOfPhotons()

int getNumOfPhotons ( ) const

inline

Returns number of photons.

Returns: number of photons

Definition at line 83 of file TOP1Dpdf.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 77 of file TOP1Dpdf.h.

◆ getTmaxFot()

double getTmaxFot ( ) const

inline

Returns maximal time of selected photons.

Returns: maximal time

Definition at line 119 of file TOP1Dpdf.h.

◆ getTmaxPDF()

double getTmaxPDF ( ) const

inline

Returns maximal time of signal PDF.

Returns: maximal time

Definition at line 107 of file TOP1Dpdf.h.

◆ getTminFot()

double getTminFot ( ) const

inline

Returns minimal time of selected photons.

Returns: minimal time

Definition at line 113 of file TOP1Dpdf.h.

◆ getTminPDF()

double getTminPDF ( ) const

inline

Returns minimal time of signal PDF.

Returns: minimal time

Definition at line 101 of file TOP1Dpdf.h.

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

top/reconstruction/include/TOP1Dpdf.h
top/reconstruction/src/TOP1Dpdf.cc

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ TOP1Dpdf()

Member Function Documentation

◆ getBinSize()

◆ getExpectedBG()

◆ getExpectedSignal()

◆ getHistogram()

◆ getLogL()

◆ getMaxT0()

◆ getMinT0()

◆ getModuleID()

◆ getNumBinsT0()

◆ getNumOfPhotons()

◆ getPhotonTimes()

◆ getTmaxFot()

◆ getTmaxPDF()

◆ getTminFot()

◆ getTminPDF()