Gamma Combination Author: Maximilian Schlupp, maxsc.nosp@m.hlup.nosp@m.p@gma.nosp@m.il.c.nosp@m.om Author: Konstantin Schubert, schub.nosp@m.ert..nosp@m.konst.nosp@m.anti.nosp@m.n@gma.nosp@m.il.c.nosp@m.om Date: October 2016. More...

#include <PDF_Datasets.h>

Inheritance diagram for PDF_Datasets:

[legend]

Collaboration diagram for PDF_Datasets:

[legend]

Public Member Functions
void	addFitObs (TString name)

bool	areObservglobalablesSet ()

bool	areParametersSet ()

void	deleteNLL ()

virtual RooFitResult *	fit (RooDataSet *dataToFit)

virtual RooFitResult *	fitBkg (RooDataSet *dataToFit)

virtual void	generateBkgToys (int SeedShift=0)

virtual void	generateBkgToysGlobalObservables (int SeedShift=0)

virtual void	generateToys (int SeedShift=0)

virtual void	generateToysGlobalObservables (int SeedShift=0)

OptParser *	getArg ()

TString	getBkgPdfName ()

RooDataSet *	getBkgToyObservables ()

TString	getConstraintName ()

RooDataSet *	getData ()

TString	getDataName ()

std::vector< TString >	getFitObs ()

int	getFitStatus ()

int	getFitStrategy ()

TString	getGlobalObsName ()

float	getMinNll ()

float	getMinNllBkg ()

float	getMinNllFree ()

float	getMinNllScan ()

TString	getObsName ()

TString	getParName ()

TString	getPdfName ()

RooDataSet *	getToyObservables ()

std::map< TString, TString >	getUnblindRegions ()

RooWorkspace *	getWorkspace ()

void	initBkgPDF (const TString &name)

void	initConstraints (const TString &setName)

void	initData (const TString &name)

void	initGlobalObservables (const TString &setName)

void	initObservables (const TString &setName)

virtual void	initObservables ()

void	initParameters (const TString &setName)

void	initParameters (const vector< TString > &parNames)

virtual void	initParameters ()

void	initPDF (const TString &name)

bool	isDataInitialized ()

bool	isPdfInitialized ()

bool	notSetupToFit (bool fitToys)

	PDF_Datasets (RooWorkspace w, int nObs, OptParser opt)
	Gamma Combination Author: Maximilian Schlupp, maxsc.nosp@m.hlup.nosp@m.p@gma.nosp@m.il.c.nosp@m.om Author: Konstantin Schubert, schub.nosp@m.ert..nosp@m.konst.nosp@m.anti.nosp@m.n@gma.nosp@m.il.c.nosp@m.om Date: October 2016. More...

	PDF_Datasets (RooWorkspace *w)

void	print ()

void	printParameters ()

void	setFitStatus (int stat=0)

void	setFitStrategy (int strat=0)

void	setMinNll (float mnll)

void	setMinNllFree (float mnll)

void	setMinNllScan (float mnll)

void	setNCPU (int n)

void	setToyData (RooDataSet *ds)

void	setVarRange (const TString &varName, const TString &rangeName, const double &rangeMin, const double &rangeMax)

void	unblind (TString var, TString unblindRegs)

	~PDF_Datasets ()

Public Member Functions inherited from PDF_Abs
bool	bkgpdfset ()

virtual void	build ()
	Shortcut for buildCov() and buildPdf(). More...

void	buildCov ()
	Build both the covariance and the correlation matrix (members covMatrix and corMatrix) from the stat and syst correlation matrics and the respective errors. More...

virtual void	buildPdf ()

virtual bool	checkConsistency ()
	Perform a couple of consistency checks to make it easier to find bugs: More...

void	deleteToys ()

TString	getBaseName ()
	Return the base name, which is the name without any unique ID. More...

RooAbsPdf *	getBkgPdf ()

TString	getCorrelationSourceString ()

TString	getErrorSourceString ()

int	getGcId ()

vector< TString >	getLatexObservables ()

TString	getName ()

int	getNobs ()

RooArgList *	getObservables ()

TString	getObservableSourceString ()

float	getObservableValue (TString obsname)
	Return the numerical value of an observable of a given name. More...

RooArgList *	getParameters ()

RooAbsPdf *	getPdf ()

void	getSubCorrelationStat (TMatrixDSym &target, vector< int > &indices)
	Return a submatrix of the statistical correlation matrix, defined by the rows and columns provided. More...

void	getSubCorrelationSyst (TMatrixDSym &target, vector< int > &indices)
	Return a submatrix of the systematic correlation matrix, defined by the rows and columns provided. More...

RooArgList *	getTheory ()

TString	getTitle ()

unsigned long long	getUniqueGlobalID ()

TString	getUniqueID ()

bool	hasObservable (TString obsname)
	Check if this PDF has an observable of the given name. More...

virtual void	initRelations ()

bool	isCrossCorPdf ()

void	loadExtParameters (RooFitResult *r)
	Set all parameters to values found in a provided fit result. More...

	PDF_Abs (int nObs)
	Gamma Combination Author: Till Moritz Karbach, morit.nosp@m.z.ka.nosp@m.rbach.nosp@m.@cer.nosp@m.n.ch Date: August 2012. More...

	PDF_Abs (int nObs, ParametersAbs &pars)

void	print () const
	Print this PDF in a verbose way: More...

void	printObservables ()

void	printParameters ()

void	resetCorrelations ()
	Set all correlations to zero. More...

bool	ScaleError (TString obsname, float scale)
	Scale the error of a given observable. More...

virtual void	setCorrelations (TString c)

void	setErrorSourceString (TString source)

void	setGcId (int id)

void	setName (TString myName)

void	setNObs (int val)

void	setObservable (TString name, float value)
	Set the observed central value of an observable. More...

virtual void	setObservables (TString c)

void	setObservableSourceString (TString source)

void	setObservablesToy ()
	Set all observables to 'toy' values drawn from the PDF using the current parameter values. More...

void	setObservablesTruth ()
	Set all observables to 'truth' values computed from the current parameters. More...

void	setSystCorrelation (TMatrixDSym &corSystMatrix)
	Set an external systematic correlation matrix. More...

void	setTitle (TString t)

virtual void	setUncertainties (TString c)

void	setUncertainty (TString obsName, float stat, float syst)
	Set the uncertainties of an observable. More...

void	storeErrorsInObs ()
	Store the errors as RooFit error into the observables to have them easily available for the pull computation. More...

virtual bool	test ()
	Test PDF implementation. More...

void	uniquify (int uID)
	used to uniquify all names when added More...

virtual	~PDF_Abs ()
	Clean off all objects in the trash bin. More...

Public Attributes
const TString	globalObsDataSnapshotName = "globalObsDataSnapshotName"

const TString	globalObsToySnapshotName = "globalObsToySnapshotName"

float	minNll

int	NCPU

Public Attributes inherited from PDF_Abs
TMatrixDSym	corMatrix

TString	corSource

TMatrixDSym	corStatMatrix

TMatrixDSym	corSystMatrix

TMatrixDSym	covMatrix

TString	obsErrSource

TString	obsValSource

vector< double >	StatErr

vector< double >	SystErr

Protected Member Functions
void	initializeRandomGenerator (int seedShift)
	Initializes the random generator. More...

Protected Member Functions inherited from PDF_Abs
void	addToTrash (TObject *)
	Add an object to the trash bin which gets emptied when this PDF object gets deleted. More...

void	getSubMatrix (TMatrixDSym &target, TMatrixDSym &source, vector< int > &indices)
	Return a submatrix of a given input matrix, defined by the rows and columns provided. More...

Protected Attributes
RooAbsPdf *	_constraintPdf

RooAbsReal *	_NLL

bool	areObsSet

bool	areParsSet

bool	areRangesSet

OptParser *	arg

TString	constraintName

RooDataSet *	data

TString	dataName

std::vector< TString >	fitObs

int	fitStatus

int	fitStrategy

TString	globalObsName

TString	globalParsName

bool	isBkgPdfSet

bool	isDataSet

bool	isPdfSet

bool	isToyDataSet

float	minNllBkg

float	minNllFree

float	minNllScan

TString	obsName

TString	parName

TString	pdfBkgName

TString	pdfName

std::map< TString, TString >	unblindRegions

RooWorkspace *	wspc

Protected Attributes inherited from PDF_Abs
int	gcId

bool	isBkgPdfSet

int	iToyObs

vector< TString >	latexObservables

bool	m_isCrossCorPdf

TString	name

int	nObs

int	nToyObs

RooArgList *	observables

RooArgList *	parameters

RooAbsPdf *	pdf

RooAbsPdf *	pdfBkg

RooArgList *	theory

TString	title

RooDataSet *	toyBkgObservables

RooDataSet *	toyObservables

map< string, TObject * >	trash

Detailed Description

Gamma Combination Author: Maximilian Schlupp, maxsc.nosp@m.hlup.nosp@m.p@gma.nosp@m.il.c.nosp@m.om Author: Konstantin Schubert, schub.nosp@m.ert..nosp@m.konst.nosp@m.anti.nosp@m.n@gma.nosp@m.il.c.nosp@m.om Date: October 2016.

Abstract class for handling generic PDFs The RooArgLists for Observables and Parameters must be provided by an external workspace. Either there are named sets within the workspace or define the variables for each set manually. The user should inherit from this class. The user has to implement the PDFs methods: fit(), generateToys()

Definition at line 22 of file PDF_Datasets.h.

Constructor & Destructor Documentation

§ PDF_Datasets() [1/2]

PDF_Datasets::PDF_Datasets	(	RooWorkspace *	w,
		int	nObs,
		OptParser *	opt
	)

Gamma Combination Author: Maximilian Schlupp, maxsc.nosp@m.hlup.nosp@m.p@gma.nosp@m.il.c.nosp@m.om Author: Konstantin Schubert, schub.nosp@m.ert..nosp@m.konst.nosp@m.anti.nosp@m.n@gma.nosp@m.il.c.nosp@m.om Date: October 2016.

Definition at line 13 of file PDF_Datasets.cpp.

§ PDF_Datasets() [2/2]

PDF_Datasets::PDF_Datasets ( RooWorkspace * w )

Definition at line 34 of file PDF_Datasets.cpp.

§ ~PDF_Datasets()

PDF_Datasets::~PDF_Datasets ( )

Definition at line 41 of file PDF_Datasets.cpp.

Member Function Documentation

§ addFitObs()

void PDF_Datasets::addFitObs ( TString name )

inline

Definition at line 47 of file PDF_Datasets.h.

§ areObservglobalablesSet()

bool PDF_Datasets::areObservglobalablesSet ( )

inline

Definition at line 83 of file PDF_Datasets.h.

§ areParametersSet()

bool PDF_Datasets::areParametersSet ( )

inline

Definition at line 84 of file PDF_Datasets.h.

§ deleteNLL()

void PDF_Datasets::deleteNLL ( )

inline

Definition at line 28 of file PDF_Datasets.h.

§ fit()

RooFitResult * PDF_Datasets::fit ( RooDataSet * dataToFit )

virtual

Reimplemented in PDF_DatasetTutorial.

Definition at line 272 of file PDF_Datasets.cpp.

§ fitBkg()

RooFitResult * PDF_Datasets::fitBkg ( RooDataSet * dataToFit )

virtual

Definition at line 303 of file PDF_Datasets.cpp.

§ generateBkgToys()

void PDF_Datasets::generateBkgToys ( int SeedShift = 0 )

virtual

Definition at line 353 of file PDF_Datasets.cpp.

§ generateBkgToysGlobalObservables()

void PDF_Datasets::generateBkgToysGlobalObservables ( int SeedShift = 0 )

virtual

Definition at line 236 of file PDF_Datasets.cpp.

§ generateToys()

void PDF_Datasets::generateToys ( int SeedShift = 0 )

virtual

Reimplemented in PDF_DatasetTutorial.

Definition at line 340 of file PDF_Datasets.cpp.

§ generateToysGlobalObservables()

void PDF_Datasets::generateToysGlobalObservables ( int SeedShift = 0 )

virtual

Definition at line 254 of file PDF_Datasets.cpp.

§ getArg()

OptParser * PDF_Datasets::getArg ( )

Definition at line 230 of file PDF_Datasets.cpp.

§ getBkgPdfName()

TString PDF_Datasets::getBkgPdfName ( )

inline

Definition at line 65 of file PDF_Datasets.h.

§ getBkgToyObservables()

RooDataSet* PDF_Datasets::getBkgToyObservables ( )

inline

Definition at line 67 of file PDF_Datasets.h.

§ getConstraintName()

TString PDF_Datasets::getConstraintName ( )

inline

Definition at line 50 of file PDF_Datasets.h.

§ getData()

RooDataSet* PDF_Datasets::getData ( )

inline

Definition at line 52 of file PDF_Datasets.h.

§ getDataName()

TString PDF_Datasets::getDataName ( )

inline

Definition at line 51 of file PDF_Datasets.h.

§ getFitObs()

std::vector<TString> PDF_Datasets::getFitObs ( )

inline

Definition at line 55 of file PDF_Datasets.h.

§ getFitStatus()

int PDF_Datasets::getFitStatus ( )

inline

Definition at line 53 of file PDF_Datasets.h.

§ getFitStrategy()

int PDF_Datasets::getFitStrategy ( )

inline

Definition at line 54 of file PDF_Datasets.h.

§ getGlobalObsName()

TString PDF_Datasets::getGlobalObsName ( )

inline

Definition at line 57 of file PDF_Datasets.h.

§ getMinNll()

float PDF_Datasets::getMinNll ( )

inline

Definition at line 58 of file PDF_Datasets.h.

§ getMinNllBkg()

float PDF_Datasets::getMinNllBkg ( )

inline

Definition at line 60 of file PDF_Datasets.h.

§ getMinNllFree()

float PDF_Datasets::getMinNllFree ( )

inline

Definition at line 59 of file PDF_Datasets.h.

§ getMinNllScan()

float PDF_Datasets::getMinNllScan ( )

inline

Definition at line 61 of file PDF_Datasets.h.

§ getObsName()

TString PDF_Datasets::getObsName ( )

inline

Definition at line 62 of file PDF_Datasets.h.

§ getParName()

TString PDF_Datasets::getParName ( )

inline

Definition at line 63 of file PDF_Datasets.h.

§ getPdfName()

TString PDF_Datasets::getPdfName ( )

inline

Definition at line 64 of file PDF_Datasets.h.

§ getToyObservables()

RooDataSet* PDF_Datasets::getToyObservables ( )

inline

Definition at line 66 of file PDF_Datasets.h.

§ getUnblindRegions()

std::map<TString,TString> PDF_Datasets::getUnblindRegions ( )

inline

Definition at line 56 of file PDF_Datasets.h.

§ getWorkspace()

RooWorkspace* PDF_Datasets::getWorkspace ( )

inline

Definition at line 68 of file PDF_Datasets.h.

§ initBkgPDF()

void PDF_Datasets::initBkgPDF ( const TString & name )

Definition at line 162 of file PDF_Datasets.cpp.

§ initConstraints()

void PDF_Datasets::initConstraints ( const TString & setName )

Definition at line 46 of file PDF_Datasets.cpp.

§ initData()

void PDF_Datasets::initData ( const TString & name )

Definition at line 56 of file PDF_Datasets.cpp.

§ initGlobalObservables()

void PDF_Datasets::initGlobalObservables ( const TString & setName )

Definition at line 89 of file PDF_Datasets.cpp.

§ initializeRandomGenerator()

void PDF_Datasets::initializeRandomGenerator ( int seedShift )

protected

Initializes the random generator.

If seedShift is set to zero, the machine environment is used to generate a hopefully unique random seed. If seedShift is nonzero, a deterministic seed is calculated from the seedShift several command line call parameters.

Definition at line 374 of file PDF_Datasets.cpp.

§ initObservables() [1/2]

void PDF_Datasets::initObservables ( const TString & setName )

Definition at line 78 of file PDF_Datasets.cpp.

§ initObservables() [2/2]

void PDF_Datasets::initObservables ( )

virtual

Reimplemented from PDF_Abs.

Definition at line 96 of file PDF_Datasets.cpp.

§ initParameters() [1/3]

void PDF_Datasets::initParameters ( const TString & setName )

Definition at line 105 of file PDF_Datasets.cpp.

§ initParameters() [2/3]

void PDF_Datasets::initParameters ( const vector< TString > & parNames )

Definition at line 119 of file PDF_Datasets.cpp.

§ initParameters() [3/3]

void PDF_Datasets::initParameters ( )

virtual

Reimplemented from PDF_Abs.

Definition at line 135 of file PDF_Datasets.cpp.

§ initPDF()

void PDF_Datasets::initPDF ( const TString & name )

Definition at line 144 of file PDF_Datasets.cpp.

§ isDataInitialized()

bool PDF_Datasets::isDataInitialized ( )

inline

Definition at line 86 of file PDF_Datasets.h.

§ isPdfInitialized()

bool PDF_Datasets::isPdfInitialized ( )

inline

Definition at line 85 of file PDF_Datasets.h.

§ notSetupToFit()

bool PDF_Datasets::notSetupToFit ( bool fitToys )

inline

Definition at line 87 of file PDF_Datasets.h.

§ print()

void PDF_Datasets::print ( )

Definition at line 205 of file PDF_Datasets.cpp.

§ printParameters()

void PDF_Datasets::printParameters ( )

Definition at line 216 of file PDF_Datasets.cpp.

§ setFitStatus()

void PDF_Datasets::setFitStatus ( int stat = 0 )

inline

Definition at line 70 of file PDF_Datasets.h.

§ setFitStrategy()

void PDF_Datasets::setFitStrategy ( int strat = 0 )

inline

Definition at line 71 of file PDF_Datasets.h.

§ setMinNll()

void PDF_Datasets::setMinNll ( float mnll )

inline

Definition at line 72 of file PDF_Datasets.h.

§ setMinNllFree()

void PDF_Datasets::setMinNllFree ( float mnll )

inline

Definition at line 73 of file PDF_Datasets.h.

§ setMinNllScan()

void PDF_Datasets::setMinNllScan ( float mnll )

inline

Definition at line 74 of file PDF_Datasets.h.

§ setNCPU()

void PDF_Datasets::setNCPU ( int n )

inline

Definition at line 75 of file PDF_Datasets.h.

§ setToyData()

void PDF_Datasets::setToyData ( RooDataSet * ds )

Definition at line 199 of file PDF_Datasets.cpp.

§ setVarRange()

void PDF_Datasets::setVarRange	(	const TString &	varName,
		const TString &	rangeName,
		const double &	rangeMin,
		const double &	rangeMax
	)

Definition at line 178 of file PDF_Datasets.cpp.

§ unblind()

void PDF_Datasets::unblind	(	TString	var,
		TString	unblindRegs
	)

Definition at line 392 of file PDF_Datasets.cpp.

Member Data Documentation

§ _constraintPdf

RooAbsPdf* PDF_Datasets::_constraintPdf

protected

Definition at line 103 of file PDF_Datasets.h.

§ _NLL

RooAbsReal* PDF_Datasets::_NLL

protected

Definition at line 102 of file PDF_Datasets.h.

§ areObsSet

bool PDF_Datasets::areObsSet

protected

Definition at line 119 of file PDF_Datasets.h.

§ areParsSet

bool PDF_Datasets::areParsSet

protected

Definition at line 120 of file PDF_Datasets.h.

§ areRangesSet

bool PDF_Datasets::areRangesSet

protected

Definition at line 121 of file PDF_Datasets.h.

§ arg

OptParser* PDF_Datasets::arg

protected

Definition at line 112 of file PDF_Datasets.h.

§ constraintName

TString PDF_Datasets::constraintName

protected

Definition at line 109 of file PDF_Datasets.h.

§ data

RooDataSet* PDF_Datasets::data

protected

Definition at line 101 of file PDF_Datasets.h.

§ dataName

TString PDF_Datasets::dataName

protected

Definition at line 108 of file PDF_Datasets.h.

§ fitObs

std::vector<TString> PDF_Datasets::fitObs

protected

Definition at line 126 of file PDF_Datasets.h.

§ fitStatus

int PDF_Datasets::fitStatus

protected

Definition at line 114 of file PDF_Datasets.h.

§ fitStrategy

int PDF_Datasets::fitStrategy

protected

Definition at line 113 of file PDF_Datasets.h.

§ globalObsDataSnapshotName

const TString PDF_Datasets::globalObsDataSnapshotName = "globalObsDataSnapshotName"

Definition at line 93 of file PDF_Datasets.h.

§ globalObsName

TString PDF_Datasets::globalObsName

protected

Definition at line 111 of file PDF_Datasets.h.

§ globalObsToySnapshotName

const TString PDF_Datasets::globalObsToySnapshotName = "globalObsToySnapshotName"

Definition at line 95 of file PDF_Datasets.h.

§ globalParsName

TString PDF_Datasets::globalParsName

protected

Definition at line 110 of file PDF_Datasets.h.

§ isBkgPdfSet

bool PDF_Datasets::isBkgPdfSet

protected

Definition at line 123 of file PDF_Datasets.h.

§ isDataSet

bool PDF_Datasets::isDataSet

protected

Definition at line 124 of file PDF_Datasets.h.

§ isPdfSet

bool PDF_Datasets::isPdfSet

protected

Definition at line 122 of file PDF_Datasets.h.

§ isToyDataSet

bool PDF_Datasets::isToyDataSet

protected

Definition at line 125 of file PDF_Datasets.h.

§ minNll

float PDF_Datasets::minNll

Definition at line 91 of file PDF_Datasets.h.

§ minNllBkg

float PDF_Datasets::minNllBkg

protected

Definition at line 117 of file PDF_Datasets.h.

§ minNllFree

float PDF_Datasets::minNllFree

protected

Definition at line 115 of file PDF_Datasets.h.

§ minNllScan

float PDF_Datasets::minNllScan

protected

Definition at line 118 of file PDF_Datasets.h.

§ NCPU

int PDF_Datasets::NCPU

Definition at line 87 of file PDF_Datasets.h.

§ obsName

TString PDF_Datasets::obsName

protected

Definition at line 106 of file PDF_Datasets.h.

§ parName

TString PDF_Datasets::parName

protected

Definition at line 107 of file PDF_Datasets.h.

§ pdfBkgName

TString PDF_Datasets::pdfBkgName

protected

Definition at line 105 of file PDF_Datasets.h.

§ pdfName

TString PDF_Datasets::pdfName

protected

Definition at line 104 of file PDF_Datasets.h.

§ unblindRegions

std::map<TString,TString> PDF_Datasets::unblindRegions

protected

Definition at line 127 of file PDF_Datasets.h.

§ wspc

RooWorkspace* PDF_Datasets::wspc

protected

Definition at line 100 of file PDF_Datasets.h.

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

/Users/matt/Scratch/lhcb/CombinationCode/gammacombo/core/include/PDF_Datasets.h
/Users/matt/Scratch/lhcb/CombinationCode/gammacombo/core/src/PDF_Datasets.cpp

Public Member Functions

Public Attributes

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

§ PDF_Datasets() [1/2]

§ PDF_Datasets() [2/2]

§ ~PDF_Datasets()

Member Function Documentation

§ addFitObs()

§ areObservglobalablesSet()

§ areParametersSet()

§ deleteNLL()

§ fit()

§ fitBkg()

§ generateBkgToys()

§ generateBkgToysGlobalObservables()

§ generateToys()

§ generateToysGlobalObservables()

§ getArg()

§ getBkgPdfName()

§ getBkgToyObservables()

§ getConstraintName()

§ getData()

§ getDataName()

§ getFitObs()

§ getFitStatus()

§ getFitStrategy()

§ getGlobalObsName()

§ getMinNll()

§ getMinNllBkg()

§ getMinNllFree()

§ getMinNllScan()

§ getObsName()

§ getParName()

§ getPdfName()

§ getToyObservables()

§ getUnblindRegions()

§ getWorkspace()

§ initBkgPDF()

§ initConstraints()

§ initData()

§ initGlobalObservables()

§ initializeRandomGenerator()

§ initObservables() [1/2]

§ initObservables() [2/2]

§ initParameters() [1/3]

§ initParameters() [2/3]

§ initParameters() [3/3]

§ initPDF()

§ isDataInitialized()

§ isPdfInitialized()

§ notSetupToFit()

§ print()

§ printParameters()

§ setFitStatus()

§ setFitStrategy()

§ setMinNll()

§ setMinNllFree()

§ setMinNllScan()

§ setNCPU()

§ setToyData()

§ setVarRange()

§ unblind()

Member Data Documentation

§ _constraintPdf

§ _NLL

§ areObsSet

§ areParsSet

§ areRangesSet

§ arg

§ constraintName

§ data

§ dataName

§ fitObs

§ fitStatus

§ fitStrategy

§ globalObsDataSnapshotName

§ globalObsName