Class for defining the SPlot technique. More...

#include <LauSPlot.hh>

Inheritance diagram for LauSPlot:

Public Types
typedef std::set< TString >	NameSet
	Type to store names, e.g. of the discriminating/control variables. More...

typedef std::map< TString, Double_t >	NumbMap
	Type to associate a category name with a double precision number, e.g. a yield or PDF value for a given species. More...

typedef std::map< TString, TLeaf * >	LeafMap
	Type to associate a variable name with the leaf of the input tree. More...

typedef std::multimap< TString, std::pair< TString, TString > >	TwoDMap
	Type to associate the name of the species that have 2D PDFs with the names of the two variables involved in each such PDF. More...

Public Member Functions
	LauSPlot (const TString &fileName, const TString &treeName, Int_t firstExpt, Int_t nExpt, const NameSet &variableNames, const LauSPlot::NumbMap &freeSpecies, const LauSPlot::NumbMap &fixdSpecies, const TwoDMap &twodimPDFs, Bool_t sigSplit=kFALSE, Bool_t scfDPSmeared=kFALSE)
	Constructor. More...

virtual	~LauSPlot ()
	Destructor. More...

void	runCalculations (const TString &option="q")
	Method to calculate the sWeights and cN coeffs. More...

void	writeOutResults ()
	Save the sWeight results as a friend tree to the input tree (in the same file) More...

const std::vector < LauSPlot::NumbMap > &	totalPdf () const
	Access the per-event total PDF values for each species. More...

Protected Member Functions
Bool_t	readInput () const
	Check whether the input tree has been successfully read. More...

void	readInput (Bool_t ok)
	Set that the input tree has been successfully read. More...

Bool_t	signalSplit () const
	Check whether the signal is split into Truth Matched and Self Cross Feed. More...

Bool_t	scfDPSmear () const
	Check whether the Self Cross Feed is smeared in the DP. More...

Bool_t	definedCNBranches () const
	Check whether the cN branches have been already created. More...

void	definedCNBranches (Bool_t defined)
	Set that the cN branches have been already defined. More...

Bool_t	definedSWeightBranches () const
	Check whether the sWeights branches have been already created. More...

void	definedSWeightBranches (Bool_t defined)
	Set that the sWeights branches have been already defined. More...

void	openInputFileAndTree ()
	Method to open the file in "update" mode and grab the input tree for reading. More...

void	readInputInfo ()
	Read the leaf structure from the tree and check the status of the read (calls LauSPlot::readInputLeaves and LauSPlot::checkLeaves) More...

Bool_t	readInputLeaves ()
	Read the leaf structure from the tree and setup the leaf map. More...

Bool_t	checkLeaves () const
	Check whether the leaf structure makes sense given the PDFs we are expecting. More...

void	createCNTree ()
	Create (if not already done) the tree for storing the cN coeffs. More...

void	defineCNBranches ()
	Create the branches for each cN coefficient. More...

void	createSWeightTree ()
	Create (if not already done) the tree for storing the sWeights. More...

void	defineSWeightBranches ()
	Create the branches to store the sWeights. More...

void	setExperiment (Int_t iExpt)
	Set the event list to contain only events from the given experiment. More...

void	readExpt ()
	Reads the values of each PDF likelihood for every event in the experiment. More...

void	checkFitter () const
	Make sure that we're using Minuit. More...

void	initialiseFitter (const TString &opt)
	Initialise Minuit, set the verbosity. More...

void	setFitParameters () const
	Add the species yields as fit parameters and fix them as appropriate. More...

void	runFit ()
	Perform the minimisation wrt the yields alone. More...

void	retrieveFittedParameters (const TString &opt)
	Update the yields with the newly fitted values and print them (unless print option is "Q"). More...

void	printCovMatrixElements (const Double_t *covmat=0) const
	Print the supplied covariance matrix or, if pointer is null, the one previously calculated. More...

void	printSumOfWeights (const TString &exclName) const
	Print the sum of sWeights for all species. More...

Bool_t	calcCovMatrix ()

void	calcTotPDFValues (const TString &exclName)
	Calculate the total likelihood for each species by multiply together all the PDFs for that species. More...

void	calcCNCoeffs (const TString &exclName, const Double_t *covmat=0)
	Computes the cN for the extended sPlots from the covariance matrix. More...

void	calcSWeights (const TString &exclName, Double_t *covmat=0)
	Computes the sWeights from the PDFs and covariance matrix. More...

void	copyEventWeights (Int_t iEvent)
	Copy the sWeight of a given event into LauSPlot::sWeightsCurrent_, from which they can be stored in the output tree. More...

void	fillCNBranches ()
	Fill the cN branches. More...

void	fillSWeightBranches ()
	Fill the sWeights branches. More...

void	addFriendTree ()
	Add the sWeightTree as a friend tree of the input tree. More...

Private Attributes
TString	fileName_
	The name of the data file. More...

TString	inputTreeName_
	The name of the input tree (containing the per-event llhds) More...

TString	cnTreeName_
	The name of the cn tree (containing the cN coefficients) More...

TString	sweightTreeName_
	The name of the sweight tree (containing the sWeights) More...

TFile *	file_
	Pointer to the data file object. More...

TTree *	inputTree_
	Pointer to the input tree. More...

TTree *	cnTree_
	Pointer to the output tree containing the cN coefficients. More...

TTree *	sweightTree_
	Pointer to the output tree containing the sWeights. More...

TEventList *	eventList_
	Pointer to an event list, that is used to loop through the experiments. More...

LeafMap	leaves_
	Collection to hold pointers to the leaves of the input tree. More...

NameSet	variableNames_
	The names of the discriminating variables. More...

NumbMap	freeSpecies_
	The names and estimated yields of the free species. More...

NumbMap	fixdSpecies_
	The names and estimated yields of the fixed species. More...

NumbMap	origFreeSpecies_
	The names and estimated yields of the free species - need to keep the original values. More...

NumbMap	origFixdSpecies_
	The names and estimated yields of the fixed species - need to keep the original values. More...

TwoDMap	twodimPDFs_
	The names of the species that have 2D PDFs and the names of the variables involved. More...

const Bool_t	signalSplit_
	Is the signal split into TM and SCF? More...

const Bool_t	scfDPSmear_
	If so then is the SCF smeared in the DP? More...

Bool_t	readInput_
	Flag whether the input tree has been successfully read. More...

Bool_t	definedCNBranches_
	Flag whether the cN branches have already been created. More...

Bool_t	definedSWeightBranches_
	Flag whether the sWeights branches have already been created. More...

Int_t	firstExpt_
	First experiment. More...

Int_t	nExpt_
	Number of experiments. More...

Int_t	iExpt_
	The current experiment. More...

Int_t	nEvents_
	Number of events in current experiment. More...

Int_t	nDiscVars_
	Number of discriminating variables. More...

Int_t	nFreeSpecies_
	Number of species free to float in the fit. More...

Int_t	nFixdSpecies_
	Number of species fixed in the fit. More...

Int_t	nSpecies_
	Total number of species (free + fixed) More...

std::vector< NumbMap >	pdfTot_
	The per-event values of the total PDF for each species. More...

std::vector< std::map< TString, NumbMap > >	discPdf_
	The per-event values of the PDFs for each species for each disc variable. More...

std::vector< Double_t >	scfFrac_
	The per-event values of the SCF fraction. More...

TMatrixD	covMat_
	The calculated covariance matrix. More...

std::vector< std::map< TString, NumbMap > >	sWeights_
	The per-event values of the computed sWeights (for each species and for each combination of excluded vars) More...

std::map< TString, NumbMap >	sWeightsCurrent_
	The current-event values of the computed sWeights. More...

std::map< TString, NumbMap >	cN_
	The extended sPlot coefficients (for each species and for each combination of excluded vars) More...

Detailed Description

Class for defining the SPlot technique.

Class for defining the SPlot technique based on TSplot from ROOT by the following authors: Muriel Pivk, Anna Kreshuk (10/2005). (Original copyright notice below)

Code extended to deal with the following two extra scenarios:

Extended sPlots (see appendix B of sPlot paper)
2D PDFs, i.e. one PDF for 2 variables

When performing a multidimensional fit, the sWeights are calculated excluding each dimension in turn and excluding none of the variables. This allows sPlots of the fit variables to be made by using the weight calculated when that variable is excluded, while variables not in the fit can be plotted from the complete information.

Definition at line 55 of file LauSPlot.hh.

Member Typedef Documentation

typedef std::map<TString,TLeaf*> LauSPlot::LeafMap

Type to associate a variable name with the leaf of the input tree.

Definition at line 65 of file LauSPlot.hh.

typedef std::set<TString> LauSPlot::NameSet

Type to store names, e.g. of the discriminating/control variables.

Definition at line 59 of file LauSPlot.hh.

typedef std::map<TString,Double_t> LauSPlot::NumbMap

Type to associate a category name with a double precision number, e.g. a yield or PDF value for a given species.

Definition at line 62 of file LauSPlot.hh.

typedef std::multimap< TString, std::pair<TString,TString> > LauSPlot::TwoDMap

Type to associate the name of the species that have 2D PDFs with the names of the two variables involved in each such PDF.

Definition at line 68 of file LauSPlot.hh.

Constructor & Destructor Documentation

LauSPlot::LauSPlot	(	const TString &	fileName,
		const TString &	treeName,
		Int_t	firstExpt,
		Int_t	nExpt,
		const NameSet &	variableNames,
		const LauSPlot::NumbMap &	freeSpecies,
		const LauSPlot::NumbMap &	fixdSpecies,
		const TwoDMap &	twodimPDFs,
		Bool_t	sigSplit = `kFALSE`,
		Bool_t	scfDPSmeared = `kFALSE`
	)

Constructor.

Parameters

[in]	fileName	the name of the data file (containing the input tree)
[in]	treeName	the name of the input tree (containing the per-event llhds)
[in]	firstExpt	the first experiment
[in]	nExpt	the number of experiments
[in]	variableNames	the names of the discriminating variables
[in]	freeSpecies	the species that are free to float in the fit
[in]	fixdSpecies	the species that are fixed in the fit
[in]	twodimPDFs	the species that have 2D PDFs and the name of the variables involved in each such PDF
[in]	sigSplit	boolean flag to check whether the signal is split into Truth Matched and Self Cross Feed
[in]	scfDPSmeared	boolean flag to check whether the SCF is smeared in the DP

Definition at line 49 of file LauSPlot.cc.

LauSPlot::~LauSPlot ( )

virtual

Destructor.

Definition at line 95 of file LauSPlot.cc.

Member Function Documentation

void LauSPlot::addFriendTree ( )

protected

Add the sWeightTree as a friend tree of the input tree.

Definition at line 1246 of file LauSPlot.cc.

void LauSPlot::calcCNCoeffs	(	const TString &	exclName,
		const Double_t *	covmat = `0`
	)

protected

Computes the cN for the extended sPlots from the covariance matrix.

Parameters

[in]	exclName	the name of excluded variable (or "none")
[in]	covmat	2D array containing the covariace matrix elements

Definition at line 1122 of file LauSPlot.cc.

Bool_t LauSPlot::calcCovMatrix ( )

protected

Returns: true if calculation succeeded or false in case the covariance matrix can not be inverted

Definition at line 874 of file LauSPlot.cc.

void LauSPlot::calcSWeights	(	const TString &	exclName,
		Double_t *	covmat = `0`
	)

protected

Computes the sWeights from the PDFs and covariance matrix.

Parameters

[in]	exclName	the name of excluded variable (or "none")
[in]	covmat	2D array containing the covariace matrix elements

Definition at line 1148 of file LauSPlot.cc.

void LauSPlot::calcTotPDFValues ( const TString & exclName )

protected

Calculate the total likelihood for each species by multiply together all the PDFs for that species.

Parameters

[in] exclName the name of excluded variable (or "none")

Definition at line 928 of file LauSPlot.cc.

void LauSPlot::checkFitter ( ) const

protected

Make sure that we're using Minuit.

Definition at line 722 of file LauSPlot.cc.

Bool_t LauSPlot::checkLeaves ( ) const

protected

Check whether the leaf structure makes sense given the PDFs we are expecting.

Definition at line 189 of file LauSPlot.cc.

void LauSPlot::copyEventWeights ( Int_t iEvent )

protected

Copy the sWeight of a given event into LauSPlot::sWeightsCurrent_, from which they can be stored in the output tree.

Parameters

[in] iEvent the requested event

Definition at line 1212 of file LauSPlot.cc.

void LauSPlot::createCNTree ( )

protected

Create (if not already done) the tree for storing the cN coeffs.

Definition at line 358 of file LauSPlot.cc.

void LauSPlot::createSWeightTree ( )

protected

Create (if not already done) the tree for storing the sWeights.

Definition at line 372 of file LauSPlot.cc.

void LauSPlot::defineCNBranches ( )

protected

Create the branches for each cN coefficient.

Definition at line 386 of file LauSPlot.cc.

Bool_t LauSPlot::definedCNBranches ( ) const

inlineprotected

Check whether the cN branches have been already created.

Returns: true/false whether cN branches have been already created

Definition at line 143 of file LauSPlot.hh.

void LauSPlot::definedCNBranches ( Bool_t defined )

inlineprotected

Set that the cN branches have been already defined.

Parameters

[in] defined set that the cN branches have been already defined

Definition at line 149 of file LauSPlot.hh.

Bool_t LauSPlot::definedSWeightBranches ( ) const

inlineprotected

Check whether the sWeights branches have been already created.

Returns: true/false whether sWeights branches have been already created

Definition at line 155 of file LauSPlot.hh.

void LauSPlot::definedSWeightBranches ( Bool_t defined )

inlineprotected

Set that the sWeights branches have been already defined.

Parameters

[in] defined set that the sWeights branches have been already defined

Definition at line 161 of file LauSPlot.hh.

void LauSPlot::defineSWeightBranches ( )

protected

Create the branches to store the sWeights.

Definition at line 424 of file LauSPlot.cc.

void LauSPlot::fillCNBranches ( )

protected

Fill the cN branches.

Definition at line 1182 of file LauSPlot.cc.

void LauSPlot::fillSWeightBranches ( )

protected

Fill the sWeights branches.

Definition at line 1193 of file LauSPlot.cc.

void LauSPlot::initialiseFitter ( const TString & opt )

protected

Initialise Minuit, set the verbosity.

Parameters

[in] opt option to set the print level (opt = "Q", "V" or "W" == "VV")

Definition at line 735 of file LauSPlot.cc.

void LauSPlot::openInputFileAndTree ( )

protected

Method to open the file in "update" mode and grab the input tree for reading.

Definition at line 107 of file LauSPlot.cc.

void LauSPlot::printCovMatrixElements ( const Double_t * covmat = 0 ) const

protected

Print the supplied covariance matrix or, if pointer is null, the one previously calculated.

Parameters

[in] covmat 2D array containing the covariace matrix elements

Definition at line 814 of file LauSPlot.cc.

void LauSPlot::printSumOfWeights ( const TString & exclName ) const

protected

Print the sum of sWeights for all species.

Parameters

[in] exclName the name of variable excluded (or "none")

Definition at line 859 of file LauSPlot.cc.

void LauSPlot::readExpt ( )

protected

Reads the values of each PDF likelihood for every event in the experiment.

Definition at line 494 of file LauSPlot.cc.

Bool_t LauSPlot::readInput ( ) const

inlineprotected

Check whether the input tree has been successfully read.

Returns: true/false whether the task have been successfully performed

Definition at line 119 of file LauSPlot.hh.

void LauSPlot::readInput ( Bool_t ok )

inlineprotected

Set that the input tree has been successfully read.

Parameters

[in] ok set that the input tree has been succeessfully read

Definition at line 125 of file LauSPlot.hh.

void LauSPlot::readInputInfo ( )

protected

Read the leaf structure from the tree and check the status of the read (calls LauSPlot::readInputLeaves and LauSPlot::checkLeaves)

Definition at line 130 of file LauSPlot.cc.

Bool_t LauSPlot::readInputLeaves ( )

protected

Read the leaf structure from the tree and setup the leaf map.

Definition at line 146 of file LauSPlot.cc.

void LauSPlot::retrieveFittedParameters ( const TString & opt )

protected

Update the yields with the newly fitted values and print them (unless print option is "Q").

Parameters

[in] opt option to set the print level (opt = "Q", "V" or "W")

Definition at line 837 of file LauSPlot.cc.

void LauSPlot::runCalculations ( const TString & option = "q" )

Method to calculate the sWeights and cN coeffs.

Parameters

[in] option control print level (default set to no print out)
"Q" - no print out (default)
"V" - prints the estimated # of events in species
"VV" - as "V" + the MINUIT printing + sums of weights for control

Definition at line 626 of file LauSPlot.cc.

void LauSPlot::runFit ( )

protected

Perform the minimisation wrt the yields alone.

Definition at line 792 of file LauSPlot.cc.

Bool_t LauSPlot::scfDPSmear ( ) const

inlineprotected

Check whether the Self Cross Feed is smeared in the DP.

Returns: true/false whether the Self Cross Feed is smeared in the DP

Definition at line 137 of file LauSPlot.hh.

void LauSPlot::setExperiment ( Int_t iExpt )

protected

Set the event list to contain only events from the given experiment.

Parameters

[in] iExpt the required experiment number

Definition at line 459 of file LauSPlot.cc.

void LauSPlot::setFitParameters ( ) const

protected

Add the species yields as fit parameters and fix them as appropriate.

Definition at line 761 of file LauSPlot.cc.

Bool_t LauSPlot::signalSplit ( ) const

inlineprotected

Check whether the signal is split into Truth Matched and Self Cross Feed.

Returns: true/false whether the signal is split into TM and SCF

Definition at line 131 of file LauSPlot.hh.

const std::vector<LauSPlot::NumbMap>& LauSPlot::totalPdf ( ) const

inline

Access the per-event total PDF values for each species.

Returns: the per-event total PDF values

Definition at line 112 of file LauSPlot.hh.

void LauSPlot::writeOutResults ( )

Save the sWeight results as a friend tree to the input tree (in the same file)

Definition at line 1224 of file LauSPlot.cc.

Member Data Documentation

std::map<TString,NumbMap> LauSPlot::cN_

private

The extended sPlot coefficients (for each species and for each combination of excluded vars)

Definition at line 349 of file LauSPlot.hh.

TTree* LauSPlot::cnTree_

private

Pointer to the output tree containing the cN coefficients.

Definition at line 282 of file LauSPlot.hh.

TString LauSPlot::cnTreeName_

private

The name of the cn tree (containing the cN coefficients)

Definition at line 274 of file LauSPlot.hh.

TMatrixD LauSPlot::covMat_

private

The calculated covariance matrix.

Definition at line 342 of file LauSPlot.hh.

Bool_t LauSPlot::definedCNBranches_

private

Flag whether the cN branches have already been created.

Definition at line 313 of file LauSPlot.hh.

Bool_t LauSPlot::definedSWeightBranches_

private

Flag whether the sWeights branches have already been created.

Definition at line 315 of file LauSPlot.hh.

std::vector<std::map<TString,NumbMap> > LauSPlot::discPdf_

private

The per-event values of the PDFs for each species for each disc variable.

Definition at line 337 of file LauSPlot.hh.

TEventList* LauSPlot::eventList_

private

Pointer to an event list, that is used to loop through the experiments.

Definition at line 287 of file LauSPlot.hh.

TFile* LauSPlot::file_

private

Pointer to the data file object.

Definition at line 278 of file LauSPlot.hh.

TString LauSPlot::fileName_

private

The name of the data file.

Definition at line 270 of file LauSPlot.hh.

Int_t LauSPlot::firstExpt_

private

First experiment.

Definition at line 318 of file LauSPlot.hh.

NumbMap LauSPlot::fixdSpecies_

private

The names and estimated yields of the fixed species.

Definition at line 297 of file LauSPlot.hh.

NumbMap LauSPlot::freeSpecies_

private

The names and estimated yields of the free species.

Definition at line 295 of file LauSPlot.hh.

Int_t LauSPlot::iExpt_

private

The current experiment.

Definition at line 322 of file LauSPlot.hh.

TTree* LauSPlot::inputTree_

private

Pointer to the input tree.

Definition at line 280 of file LauSPlot.hh.

TString LauSPlot::inputTreeName_

private

The name of the input tree (containing the per-event llhds)

Definition at line 272 of file LauSPlot.hh.

LeafMap LauSPlot::leaves_

private

Collection to hold pointers to the leaves of the input tree.

Definition at line 290 of file LauSPlot.hh.

Int_t LauSPlot::nDiscVars_

private

Number of discriminating variables.

Definition at line 326 of file LauSPlot.hh.

Int_t LauSPlot::nEvents_

private

Number of events in current experiment.

Definition at line 324 of file LauSPlot.hh.

Int_t LauSPlot::nExpt_

private

Number of experiments.

Definition at line 320 of file LauSPlot.hh.

Int_t LauSPlot::nFixdSpecies_

private

Number of species fixed in the fit.

Definition at line 330 of file LauSPlot.hh.

Int_t LauSPlot::nFreeSpecies_

private

Number of species free to float in the fit.

Definition at line 328 of file LauSPlot.hh.

Int_t LauSPlot::nSpecies_

private

Total number of species (free + fixed)

Definition at line 332 of file LauSPlot.hh.

NumbMap LauSPlot::origFixdSpecies_

private

The names and estimated yields of the fixed species - need to keep the original values.

Definition at line 301 of file LauSPlot.hh.

NumbMap LauSPlot::origFreeSpecies_

private

The names and estimated yields of the free species - need to keep the original values.

Definition at line 299 of file LauSPlot.hh.

std::vector<NumbMap> LauSPlot::pdfTot_

private

The per-event values of the total PDF for each species.

Definition at line 335 of file LauSPlot.hh.

Bool_t LauSPlot::readInput_

private

Flag whether the input tree has been successfully read.

Definition at line 311 of file LauSPlot.hh.

const Bool_t LauSPlot::scfDPSmear_

private

If so then is the SCF smeared in the DP?

Definition at line 308 of file LauSPlot.hh.

std::vector<Double_t> LauSPlot::scfFrac_

private

The per-event values of the SCF fraction.

Definition at line 339 of file LauSPlot.hh.

const Bool_t LauSPlot::signalSplit_

private

Is the signal split into TM and SCF?

Definition at line 306 of file LauSPlot.hh.

std::vector< std::map<TString,NumbMap> > LauSPlot::sWeights_

private

The per-event values of the computed sWeights (for each species and for each combination of excluded vars)

Definition at line 345 of file LauSPlot.hh.

std::map<TString,NumbMap> LauSPlot::sWeightsCurrent_

private

The current-event values of the computed sWeights.

Definition at line 347 of file LauSPlot.hh.

TTree* LauSPlot::sweightTree_

private

Pointer to the output tree containing the sWeights.

Definition at line 284 of file LauSPlot.hh.

TString LauSPlot::sweightTreeName_

private

The name of the sweight tree (containing the sWeights)

Definition at line 276 of file LauSPlot.hh.

TwoDMap LauSPlot::twodimPDFs_

private

The names of the species that have 2D PDFs and the names of the variables involved.

Definition at line 303 of file LauSPlot.hh.

NameSet LauSPlot::variableNames_

private

The names of the discriminating variables.

Definition at line 293 of file LauSPlot.hh.

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

inc/LauSPlot.hh
src/LauSPlot.cc

Public Types

Public Member Functions

Protected Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation