Class for defining a 2D histogram PDF. More...

#include <Lau2DHistPdf.hh>

Inheritance diagram for Lau2DHistPdf:

Public Member Functions
	Lau2DHistPdf (const std::vector< TString > &theVarNames, const TH2 *hist, const LauFitData &minVals, const LauFitData &maxVals, Bool_t useInterpolation=kTRUE, Bool_t fluctuateBins=kFALSE)
	Constructor. More...

virtual	~Lau2DHistPdf ()
	Destructor. More...

	Lau2DHistPdf (const Lau2DHistPdf &other)
	Copy constructor. More...

Double_t	getMinX () const
	Get the minimum value of x-axis abscissa. More...

Double_t	getMaxX () const
	Get the maximum value of x-axis abscissa. More...

Double_t	getRangeX () const
	Get the range of x-axis abscissa values. More...

Double_t	getMinY () const
	Get the minimum value of y-axis abscissa. More...

Double_t	getMaxY () const
	Get the maximum value of y-axis abscissa. More...

Double_t	getRangeY () const
	Get the range of y-axis abscissa values. More...

virtual void	cacheInfo (const LauFitDataTree &inputData)
	Cache information from data. More...

virtual void	calcLikelihoodInfo (const LauAbscissas &abscissas)
	Calculate the likelihood (and intermediate info) for a given value of the abscissas. More...

virtual void	calcLikelihoodInfo (UInt_t iEvt)
	Calculate the likelihood (and intermediate info) for a given event number. More...

virtual Double_t	getLikelihood (const TString &theVarName) const
	Get likelihood only for the given variable. More...

virtual void	checkPositiveness ()
	Check that PDF is positive. More...

virtual void	calcNorm ()
	Calculate the normalisation. More...

virtual void	calcPDFHeight (const LauKinematics *kinematics)
	Calculate the PDF height. More...

virtual LauFitData	generate (const LauKinematics *kinematics)
	Generate an event from the PDF. More...

Public Member Functions inherited from LauAbsPdf
	LauAbsPdf (const TString &theVarName, const std::vector< LauParameter * > &params, Double_t minAbscissa, Double_t maxAbscissa)
	Constructor for a 1D PDF. More...

	LauAbsPdf (const std::vector< TString > &theVarNames, const std::vector< LauParameter * > &params, const LauFitData &minAbscissas, const LauFitData &maxAbscissas)
	Constructor for a multidimensional PDF. More...

virtual	~LauAbsPdf ()
	Destructor. More...

virtual const TString &	varName () const
	Retrieve the name of the abscissa. More...

virtual std::vector< TString >	varNames () const
	Retrieve the names of the abscissas. More...

virtual UInt_t	nParameters () const
	Retrieve the number of PDF parameters. More...

virtual UInt_t	nFixedParameters () const
	Retrieve the number of fixed PDF parameters. More...

virtual UInt_t	nInputVars () const
	Retrieve the number of abscissas. More...

virtual Bool_t	isDPDependent () const
	Specifies whether or not the PDF is DP dependent. More...

virtual Double_t	getParValue (const TString &parName) const
	Retrieve the value of the specified parameter. More...

virtual Double_t	getParMin (const TString &parName) const
	Retrieve the minimum value of the specified parameter. More...

virtual Double_t	getParMax (const TString &parName) const
	Retrieve the maximum value of the specified parameter. More...

virtual Double_t	getParRange (const TString &parName) const
	Retrieve the range of the specified parameter. More...

virtual Bool_t	parFixed (const TString &parName) const
	Retrieve whether the specified parameter is fixed. More...

virtual Bool_t	parClone (const TString &parName) const
	Retrieve whether the specified parameter is a clone. More...

virtual Double_t	getMinAbscissa () const
	Retrieve the minimum value of the (primary) abscissa. More...

virtual Double_t	getMaxAbscissa () const
	Retrieve the maximum value of the (primary) abscissa. More...

virtual Double_t	getRange () const
	Retrieve the range of the (primary) abscissa. More...

virtual Double_t	getMinAbscissa (const TString &theVarName) const
	Retrieve the minimum value of the specified abscissa. More...

virtual Double_t	getMaxAbscissa (const TString &theVarName) const
	Retrieve the maximum value of the specified abscissa. More...

virtual Double_t	getRange (const TString &theVarName) const
	Retrieve the range of the specified abscissa. More...

virtual LauFitData	getMinAbscissas () const
	Retrieve the minimum values of all the abscissas. More...

virtual LauFitData	getMaxAbscissas () const
	Retrieve the maximum values of all the abscissas. More...

virtual LauFitData	getRanges () const
	Retrieve the ranges of all the abscissas. More...

virtual void	setParValue (const TString &parName, Double_t value)
	Change the value of the specified parameter. More...

virtual void	setParMin (const TString &parName, Double_t minValue)
	Change the minimum value of the specified parameter. More...

virtual void	setParMax (const TString &parName, Double_t maxValue)
	Change the maximum value of the specified parameter. More...

virtual void	setParRange (const TString &parName, Double_t minValue, Double_t maxValue)
	Change the range of the specified parameter. More...

virtual void	fixPar (const TString &parName)
	Fix the specified parameter. More...

virtual void	floatPar (const TString &parName)
	Float the specified parameter. More...

virtual void	updatePulls ()
	Update the pulls for all parameters. More...

virtual Double_t	getUnNormLikelihood () const
	Retrieve the unnormalised likelihood value. More...

virtual Double_t	getNorm () const
	Retrieve the normalisation factor. More...

virtual Double_t	getLikelihood () const
	Retrieve the normalised likelihood value. More...

virtual Double_t	getMaxHeight () const
	Retrieve the maximum height. More...

virtual void	setRandomFun (TRandom *randomFun)
	Set the random function used for toy MC generation. More...

virtual const std::vector < LauParameter * > &	getParameters () const
	Retrieve the parameters of the PDF, e.g. so that they can be loaded into a fit. More...

virtual std::vector < LauParameter * > &	getParameters ()
	Retrieve the parameters of the PDF, e.g. so that they can be loaded into a fit. More...

virtual Bool_t	heightUpToDate () const
	Check if the maximum height of the PDF is up to date. More...

virtual Bool_t	cachePDF () const
	Check if the PDF is to be cached. More...

virtual Int_t	nNormPoints () const
	Retrieve the number of points to integrate over when normalising. More...

virtual void	nNormPoints (Int_t nPoints)
	Set the number of points to integrate over when normalising. More...

virtual IntMethod	integMethod () const
	Retrieve the integration method used to normalise the PDF. More...

virtual void	integMethod (IntMethod method)
	Set the integration method used to normalise the PDF. More...

Protected Member Functions
void	doBinFluctuation ()
	Fluctuate the histogram bin contents in accordance with their errors. More...

void	checkNormalisation ()
	Check the normalisation calculation. More...

Double_t	getBinHistValue (Int_t i, Int_t j) const
	Get the bin content from the histogram. More...

Double_t	interpolateXY (Double_t x, Double_t y) const
	Perform the interpolation (unnormalised) More...

Double_t	interpolateXYNorm (Double_t x, Double_t y) const
	Perform the interpolation and divide by the normalisation. More...

Protected Member Functions inherited from LauAbsPdf
virtual void	cachePDF (Bool_t doCachePDF)
	Set whether the PDF is to be cached. More...

virtual void	heightUpToDate (Bool_t hutd)
	Set whether the height is up to date. More...

virtual Double_t	integrGaussLegendre ()
	Integrate the PDF using the Gauss-Legendre method. More...

virtual Double_t	integTrapezoid ()
	Integrate the PDF using the simple trapezoid method. More...

virtual void	setNorm (Double_t norm)
	Set the normalisation factor. More...

virtual void	setMaxHeight (Double_t maxHeight)
	Set the maximum height. More...

virtual void	setMinAbscissa (const TString &theVarName, Double_t minAbscissa)
	Set the minimum value of the specified abscissa. More...

virtual void	setMaxAbscissa (const TString &theVarName, Double_t maxAbscissa)
	Set the maximum value of the specified abscissa. More...

virtual void	setRange (const TString &theVarName, Double_t minAbscissa, Double_t maxAbscissa)
	Set the range of the specified abscissa. More...

virtual Bool_t	checkRange (const LauAbscissas &abscissas) const
	Check that all abscissas are within their allowed ranges. More...

virtual void	setUnNormPDFVal (Double_t unNormPDFVal)
	Set the unnormalised likelihood. More...

virtual LauParameter *	findParameter (const TString &parName)
	Retrieve the specified parameter. More...

virtual const LauParameter *	findParameter (const TString &parName) const
	Retrieve the specified parameter. More...

virtual TRandom *	getRandomFun () const
	Retrieve the random function used for MC generation. More...

virtual std::vector < LauAbscissas > &	getAbscissas ()
	Retrieve the abscissa(s) More...

virtual const std::vector < LauAbscissas > &	getAbscissas () const
	Retrieve the abscissa(s) More...

virtual std::vector< Double_t > &	getUnNormPDFValues ()
	Retrieve the cached unnormalised likelihood values. More...

virtual const std::vector < Double_t > &	getUnNormPDFValues () const
	Retrieve the cached unnormalised likelihood values. More...

virtual void	addParameters (std::vector< LauParameter * > &params)
	Add parameters to the PDF. More...

virtual Bool_t	withinNormCalc () const
	Check whether the calcNorm method is running. More...

virtual void	withinNormCalc (Bool_t yorn)
	Set flag to track whether the calcNorm method is running. More...

virtual Bool_t	withinGeneration () const
	Check whether the generate method is running. More...

virtual void	withinGeneration (Bool_t yorn)
	Set flag to track whether the generate method is running. More...

virtual Bool_t	normWeightsDone () const
	Check whether the normalisation weights have been calculated. More...

virtual void	normWeightsDone (Bool_t yorn)
	Set whether the normalisation weights have been calculated. More...

virtual void	getNormWeights ()
	Calculate the weights and abscissas used for normalisation. More...

virtual const std::vector < LauAbscissas > &	normAbscissas () const
	Retrieve the abscissa points used for normalisation. More...

virtual const std::vector < Double_t > &	normWeights () const
	Retrieve the weights used for normalisation. More...

Private Attributes
TH2 *	hist_
	The underlying histogram. More...

TH1 *	xProj_
	Projection of histogram x-axis. More...

TH1 *	yProj_
	Projection of histogram y-axis. More...

Lau1DHistPdf *	xVarPdf_
	1D PDF for x variable More...

Lau1DHistPdf *	yVarPdf_
	1D PDF for y variable More...

TString	xName_
	x variable name More...

TString	yName_
	y variable name More...

Int_t	nBinsX_
	The number of bins on the x-axis of the histogram. More...

Int_t	nBinsY_
	The number of bins on the y-axis of the histogram. More...

Double_t	minX_
	The histogram x-axis minimum. More...

Double_t	maxX_
	The histogram x-axis maximum. More...

Double_t	minY_
	The histogram y-axis minimum. More...

Double_t	maxY_
	The histogram y-axis maximum. More...

Double_t	rangeX_
	The histogram x-axis range. More...

Double_t	rangeY_
	The histogram y-axis range. More...

Double_t	binXWidth_
	The histogram x-axis bin width. More...

Double_t	binYWidth_
	The histogram y-axis bin width. More...

Double_t	invBinXWidth_
	The histogram x-axis inverse bin width. More...

Double_t	invBinYWidth_
	The histogram y-axis inverse bin width. More...

Bool_t	useInterpolation_
	Control boolean for using the linear interpolation. More...

Bool_t	fluctuateBins_
	Control boolean for performing the fluctuation of the histogram bin contents. More...

Additional Inherited Members
Public Types inherited from LauAbsPdf
enum	IntMethod { GaussLegendre, Trapezoid }
	The possible numerical intergration methods. More...

typedef std::vector< Double_t >	LauAbscissas
	The type used for containing multiple abscissa values. More...

Detailed Description

Class for defining a 2D histogram PDF.

Class for defining a 2D histogram PDF. Employs linear interpolation to get the histogram value based on how far away a point is to nearby bin centres. The returned values are normalised to the total area.

Definition at line 33 of file Lau2DHistPdf.hh.

Constructor & Destructor Documentation

Lau2DHistPdf::Lau2DHistPdf	(	const std::vector< TString > &	theVarNames,
		const TH2 *	hist,
		const LauFitData &	minVals,
		const LauFitData &	maxVals,
		Bool_t	useInterpolation = `kTRUE`,
		Bool_t	fluctuateBins = `kFALSE`
	)

Constructor.

Parameters

[in]	theVarNames	the names of the abscissa variables
[in]	hist	the 2D histogram from which the PDF should be constructed
[in]	minVals	the minimum values of the abscissas
[in]	maxVals	the maximum values of the abscissas
[in]	useInterpolation	boolean flag to determine whether linear interpolation between bins should be used or simply the raw bin values
[in]	fluctuateBins	boolean flag to determine whether the bin contents should be fluctuated in accordance with their errors

Definition at line 33 of file Lau2DHistPdf.cc.

Lau2DHistPdf::~Lau2DHistPdf ( )

virtual

Destructor.

Definition at line 151 of file Lau2DHistPdf.cc.

Lau2DHistPdf::Lau2DHistPdf ( const Lau2DHistPdf & other )

Copy constructor.

Definition at line 161 of file Lau2DHistPdf.cc.

Member Function Documentation

void Lau2DHistPdf::cacheInfo ( const LauFitDataTree & inputData )

virtual

Cache information from data.

Parameters

[in] inputData the input data

Reimplemented from LauAbsPdf.

Definition at line 431 of file Lau2DHistPdf.cc.

void Lau2DHistPdf::calcLikelihoodInfo ( const LauAbscissas & abscissas )

virtual

Calculate the likelihood (and intermediate info) for a given value of the abscissas.

Parameters

[in] abscissas the abscissa values

Implements LauAbsPdf.

Definition at line 390 of file Lau2DHistPdf.cc.

void Lau2DHistPdf::calcLikelihoodInfo ( UInt_t iEvt )

virtual

Calculate the likelihood (and intermediate info) for a given event number.

Parameters

[in] iEvt event number

Reimplemented from LauAbsPdf.

Definition at line 409 of file Lau2DHistPdf.cc.

void Lau2DHistPdf::calcNorm ( )

virtual

Calculate the normalisation.

Reimplemented from LauAbsPdf.

Definition at line 202 of file Lau2DHistPdf.cc.

void Lau2DHistPdf::calcPDFHeight ( const LauKinematics * kinematics )

virtual

Calculate the PDF height.

Parameters

[in] kinematics the current DP kinematics (not strictly required in this case since PDF has no DP dependence)

Implements LauAbsPdf.

Definition at line 190 of file Lau2DHistPdf.cc.

void Lau2DHistPdf::checkNormalisation ( )

protected

Check the normalisation calculation.

Definition at line 227 of file Lau2DHistPdf.cc.

virtual void Lau2DHistPdf::checkPositiveness ( )

inlinevirtual

Check that PDF is positive.

Check that PDF is positive. Dealt with in getBinHistValue, so nothing to do here.

Implements LauAbsPdf.

Definition at line 105 of file Lau2DHistPdf.hh.

void Lau2DHistPdf::doBinFluctuation ( )

protected

Fluctuate the histogram bin contents in accordance with their errors.

Definition at line 441 of file Lau2DHistPdf.cc.

LauFitData Lau2DHistPdf::generate ( const LauKinematics * kinematics )

virtual

Generate an event from the PDF.

Parameters

[in] kinematics the current DP kinematics (used to determine the DP position by some PDFs that have dependence on it)

Returns: generated data

Reimplemented from LauAbsPdf.

Definition at line 458 of file Lau2DHistPdf.cc.

Double_t Lau2DHistPdf::getBinHistValue	(	Int_t	i,
		Int_t	j
	)		const

protected

Get the bin content from the histogram.

Parameters

[in]	i	the x-axis bin number
[in]	j	the y-axis bin number

Returns: the bin conent

Definition at line 257 of file Lau2DHistPdf.cc.

Double_t Lau2DHistPdf::getLikelihood ( const TString & theVarName ) const

virtual

Get likelihood only for the given variable.

Parameters

[in] theVarName the variable name (must be one of the two dimensions of this histogram!)

Returns: likelihood value

Reimplemented from LauAbsPdf.

Definition at line 419 of file Lau2DHistPdf.cc.

Double_t Lau2DHistPdf::getMaxX ( ) const

inline

Get the maximum value of x-axis abscissa.

Definition at line 59 of file Lau2DHistPdf.hh.

Double_t Lau2DHistPdf::getMaxY ( ) const

inline

Get the maximum value of y-axis abscissa.

Definition at line 68 of file Lau2DHistPdf.hh.

Double_t Lau2DHistPdf::getMinX ( ) const

inline

Get the minimum value of x-axis abscissa.

Definition at line 56 of file Lau2DHistPdf.hh.

Double_t Lau2DHistPdf::getMinY ( ) const

inline

Get the minimum value of y-axis abscissa.

Definition at line 65 of file Lau2DHistPdf.hh.

Double_t Lau2DHistPdf::getRangeX ( ) const

inline

Get the range of x-axis abscissa values.

Definition at line 62 of file Lau2DHistPdf.hh.

Double_t Lau2DHistPdf::getRangeY ( ) const

inline

Get the range of y-axis abscissa values.

Definition at line 71 of file Lau2DHistPdf.hh.

Double_t Lau2DHistPdf::interpolateXY	(	Double_t	x,
		Double_t	y
	)		const

protected

Perform the interpolation (unnormalised)

Parameters

[in]	x	the x-axis abscissa value
[in]	y	the y-aixs abscissa value

Returns: the unnormalised PDF value

Definition at line 284 of file Lau2DHistPdf.cc.

Double_t Lau2DHistPdf::interpolateXYNorm	(	Double_t	x,
		Double_t	y
	)		const

protected

Perform the interpolation and divide by the normalisation.

Parameters

[in]	x	the x-axis abscissa value
[in]	y	the y-axis abscissa value

Returns: the normalised PDF value

Definition at line 276 of file Lau2DHistPdf.cc.

Member Data Documentation

Double_t Lau2DHistPdf::binXWidth_

private

The histogram x-axis bin width.

Definition at line 190 of file Lau2DHistPdf.hh.

Double_t Lau2DHistPdf::binYWidth_

private

The histogram y-axis bin width.

Definition at line 192 of file Lau2DHistPdf.hh.

Bool_t Lau2DHistPdf::fluctuateBins_

private

Control boolean for performing the fluctuation of the histogram bin contents.

Definition at line 202 of file Lau2DHistPdf.hh.

TH2* Lau2DHistPdf::hist_

private

The underlying histogram.

Definition at line 156 of file Lau2DHistPdf.hh.

Double_t Lau2DHistPdf::invBinXWidth_

private

The histogram x-axis inverse bin width.

Definition at line 195 of file Lau2DHistPdf.hh.

Double_t Lau2DHistPdf::invBinYWidth_

private

The histogram y-axis inverse bin width.

Definition at line 197 of file Lau2DHistPdf.hh.

Double_t Lau2DHistPdf::maxX_

private

The histogram x-axis maximum.

Definition at line 179 of file Lau2DHistPdf.hh.

Double_t Lau2DHistPdf::maxY_

private

The histogram y-axis maximum.

Definition at line 183 of file Lau2DHistPdf.hh.

Double_t Lau2DHistPdf::minX_

private

The histogram x-axis minimum.

Definition at line 177 of file Lau2DHistPdf.hh.

Double_t Lau2DHistPdf::minY_

private

The histogram y-axis minimum.

Definition at line 181 of file Lau2DHistPdf.hh.

Int_t Lau2DHistPdf::nBinsX_

private

The number of bins on the x-axis of the histogram.

Definition at line 173 of file Lau2DHistPdf.hh.

Int_t Lau2DHistPdf::nBinsY_

private

The number of bins on the y-axis of the histogram.

Definition at line 175 of file Lau2DHistPdf.hh.

Double_t Lau2DHistPdf::rangeX_

private

The histogram x-axis range.

Definition at line 185 of file Lau2DHistPdf.hh.

Double_t Lau2DHistPdf::rangeY_

private

The histogram y-axis range.

Definition at line 187 of file Lau2DHistPdf.hh.

Bool_t Lau2DHistPdf::useInterpolation_

private

Control boolean for using the linear interpolation.

Definition at line 200 of file Lau2DHistPdf.hh.

TString Lau2DHistPdf::xName_

private

x variable name

Definition at line 168 of file Lau2DHistPdf.hh.

TH1* Lau2DHistPdf::xProj_

private

Projection of histogram x-axis.

Definition at line 158 of file Lau2DHistPdf.hh.

Lau1DHistPdf* Lau2DHistPdf::xVarPdf_

private

1D PDF for x variable

Definition at line 163 of file Lau2DHistPdf.hh.

TString Lau2DHistPdf::yName_

private

y variable name

Definition at line 170 of file Lau2DHistPdf.hh.

TH1* Lau2DHistPdf::yProj_

private

Projection of histogram y-axis.

Definition at line 160 of file Lau2DHistPdf.hh.

Lau1DHistPdf* Lau2DHistPdf::yVarPdf_

private

1D PDF for y variable

Definition at line 165 of file Lau2DHistPdf.hh.

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

inc/Lau2DHistPdf.hh
src/Lau2DHistPdf.cc

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation