laura is hosted by Hepforge, IPPP Durham
Laura++  v3r0p1
A maximum likelihood fitting package for performing Dalitz-plot analysis.

Class to allow having different PDFs in different regions of the DP. More...

#include <LauDPDepMapPdf.hh>

Inheritance diagram for LauDPDepMapPdf:
LauAbsPdf

Public Types

enum  DPAxis {
  M12, M13, M23, CentreDist,
  MMIN, MMAX
}
 Define possibilties for the DP axes. More...
 
- 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...
 

Public Member Functions

 LauDPDepMapPdf (const std::vector< LauAbsPdf * > &pdfs, const LauDaughters *daughters, const TH2 *dpHisto, Bool_t upperHalf=kFALSE)
 Constructor - map described by 2D histogram. More...
 
 LauDPDepMapPdf (const std::vector< LauAbsPdf * > &pdfs, const LauDaughters *daughters, const TH1 *dpAxisHisto, DPAxis dpAxis)
 Constructor - map described by 1D histogram of a DP "axis". More...
 
virtual ~LauDPDepMapPdf ()
 Destructor. More...
 
virtual Bool_t isDPDependent () const
 Specifies whether or not the PDF is DP dependent. 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 abscissa. More...
 
virtual void calcLikelihoodInfo (UInt_t iEvt)
 Calculate the likelihood (and intermediate info) for a given event number. More...
 
virtual void calcNorm ()
 Calculate the normalisation. More...
 
virtual void calcPDFHeight (const LauKinematics *kinematics)
 Calculate the PDF height. More...
 
- Public Member Functions inherited from LauAbsPdf
 LauAbsPdf (const TString &theVarName, const std::vector< LauAbsRValue * > &params, Double_t minAbscissa, Double_t maxAbscissa)
 Constructor for a 1D PDF. More...
 
 LauAbsPdf (const std::vector< TString > &theVarNames, const std::vector< LauAbsRValue * > &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 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 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 getLikelihood (const TString &theVarName) const
 For multidimentional PDFs, retrieve the normalised likelihood value of a named variable. More...
 
virtual Double_t getMaxHeight () const
 Retrieve the maximum height. More...
 
virtual LauFitData generate (const LauKinematics *kinematics)
 Generate an event from the PDF. More...
 
virtual void setRandomFun (TRandom *randomFun)
 Set the random function used for toy MC generation. More...
 
virtual const std::vector
< LauAbsRValue * > & 
getParameters () const
 Retrieve the parameters of the PDF, e.g. so that they can be loaded into a fit. More...
 
virtual std::vector
< LauAbsRValue * > & 
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 void heightUpToDate (Bool_t hutd)
 Set whether the height 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

UInt_t determineDPRegion (Double_t m13Sq, Double_t m23Sq) const
 Determine the DP region. More...
 
- Protected Member Functions inherited from LauAbsPdf
virtual void cachePDF (Bool_t doCachePDF)
 Set whether the PDF is to be cached. 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 LauAbsRValuefindParameter (const TString &parName)
 Retrieve the specified parameter. More...
 
virtual const LauAbsRValuefindParameter (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< LauAbsRValue * > &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 Member Functions

 LauDPDepMapPdf (const LauDPDepMapPdf &other)
 Copy constructor - not implemented. More...
 
LauDPDepMapPdfoperator= (const LauDPDepMapPdf &other)
 Copy assignment operator - not implemented. More...
 

Private Attributes

LauDaughtersdaughters_
 Daughter particles. More...
 
std::vector< LauAbsPdf * > pdfs_
 The PDFs. More...
 
Lau2DAbsDPdpDependence_
 2D histogram - DP More...
 
TH1 * dpAxisDependence_
 1D histogram - DP axis More...
 
DPAxis dpAxis_
 The DP axis we depend on. More...
 
std::vector< UInt_t > indices_
 Cached indices values. More...
 

Detailed Description

Class to allow having different PDFs in different regions of the DP.

Class to allow having different PDFs in different regions of the DP. The DP regions are specified in a 2D histogram or in a 1D histogram of some projection of the DP. Each PDF is then associated with that region.

Definition at line 40 of file LauDPDepMapPdf.hh.

Member Enumeration Documentation

Define possibilties for the DP axes.

Enumerator
M12 

dependence is on m^2_12

M13 

dependence is on m^2_13

M23 

dependence is on m^2_23

CentreDist 

dependence is on the distance from the DP centre

MMIN 

dependence is on the minimum of m^2_13 and m^2_23

MMAX 

dependence is on the maximum of m^2_13 and m^2_23

Definition at line 44 of file LauDPDepMapPdf.hh.

Constructor & Destructor Documentation

LauDPDepMapPdf::LauDPDepMapPdf ( const std::vector< LauAbsPdf * > &  pdfs,
const LauDaughters daughters,
const TH2 *  dpHisto,
Bool_t  upperHalf = kFALSE 
)

Constructor - map described by 2D histogram.

Parameters
[in]pdfsthe PDFs
[in]daughtersthe daughter particles
[in]dpHistothe 2D histogram
[in]upperHalfspecifies whether only upper half of DP is supplied in the case of a symmetric DP

Definition at line 34 of file LauDPDepMapPdf.cc.

LauDPDepMapPdf::LauDPDepMapPdf ( const std::vector< LauAbsPdf * > &  pdfs,
const LauDaughters daughters,
const TH1 *  dpAxisHisto,
DPAxis  dpAxis 
)

Constructor - map described by 1D histogram of a DP "axis".

Parameters
[in]pdfsthe PDFs
[in]daughtersthe daughter particles
[in]dpAxisHistothe 1D histogram
[in]dpAxisthe DP axis

Definition at line 94 of file LauDPDepMapPdf.cc.

LauDPDepMapPdf::~LauDPDepMapPdf ( )
virtual

Destructor.

Definition at line 159 of file LauDPDepMapPdf.cc.

LauDPDepMapPdf::LauDPDepMapPdf ( const LauDPDepMapPdf other)
private

Copy constructor - not implemented.

Member Function Documentation

void LauDPDepMapPdf::cacheInfo ( const LauFitDataTree inputData)
virtual

Cache information from data.

Parameters
[in]inputDatathe input data

Reimplemented from LauAbsPdf.

Definition at line 297 of file LauDPDepMapPdf.cc.

void LauDPDepMapPdf::calcLikelihoodInfo ( const LauAbscissas abscissas)
virtual

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

Parameters
[in]abscissasthe values of the abscissa(s)

Implements LauAbsPdf.

Definition at line 208 of file LauDPDepMapPdf.cc.

void LauDPDepMapPdf::calcLikelihoodInfo ( UInt_t  iEvt)
virtual

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

Parameters
[in]iEvtevent number

Reimplemented from LauAbsPdf.

Definition at line 340 of file LauDPDepMapPdf.cc.

void LauDPDepMapPdf::calcNorm ( )
virtual

Calculate the normalisation.

Reimplemented from LauAbsPdf.

Definition at line 258 of file LauDPDepMapPdf.cc.

void LauDPDepMapPdf::calcPDFHeight ( const LauKinematics kinematics)
virtual

Calculate the PDF height.

Parameters
[in]kinematicsthe current DP kinematics

Implements LauAbsPdf.

Definition at line 264 of file LauDPDepMapPdf.cc.

UInt_t LauDPDepMapPdf::determineDPRegion ( Double_t  m13Sq,
Double_t  m23Sq 
) const
protected

Determine the DP region.

Parameters
[in]m13Sqthe invariant mass squared of daughters 1 and 3
[in]m23Sqthe invariant mass squared of daughters 2 and 3

Definition at line 167 of file LauDPDepMapPdf.cc.

virtual Bool_t LauDPDepMapPdf::isDPDependent ( ) const
inlinevirtual

Specifies whether or not the PDF is DP dependent.

Returns
true if the PDF is DP-dependent (the default)

Reimplemented from LauAbsPdf.

Definition at line 83 of file LauDPDepMapPdf.hh.

LauDPDepMapPdf& LauDPDepMapPdf::operator= ( const LauDPDepMapPdf other)
private

Copy assignment operator - not implemented.

Member Data Documentation

LauDaughters* LauDPDepMapPdf::daughters_
private

Daughter particles.

Definition at line 128 of file LauDPDepMapPdf.hh.

DPAxis LauDPDepMapPdf::dpAxis_
private

The DP axis we depend on.

Definition at line 139 of file LauDPDepMapPdf.hh.

TH1* LauDPDepMapPdf::dpAxisDependence_
private

1D histogram - DP axis

Definition at line 136 of file LauDPDepMapPdf.hh.

Lau2DAbsDP* LauDPDepMapPdf::dpDependence_
private

2D histogram - DP

Definition at line 134 of file LauDPDepMapPdf.hh.

std::vector<UInt_t> LauDPDepMapPdf::indices_
private

Cached indices values.

Definition at line 142 of file LauDPDepMapPdf.hh.

std::vector<LauAbsPdf*> LauDPDepMapPdf::pdfs_
private

The PDFs.

Definition at line 131 of file LauDPDepMapPdf.hh.


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