#include <MethodBergerBoosScan.h>

Inheritance diagram for MethodBergerBoosScan:

Collaboration diagram for MethodBergerBoosScan:

Public Member Functions
TH2F *	calcPValues (TH2F better, TH2F all, TH2F bg)
	Calculates 2D p-Values in (scanpoint,BergerBoos_id) Histogram. More...

void	drawBBPoints (TString varX, TString varY, int runMin=1, int runMax=1, bool save=true)
	Draws 2D Histogram showing the BB points in varX-varY space. More...

void	getBestPValue (TH1F hCL, TH2F pValues)
	Returns a 1D Histogram containing the 1-CL curve Input is a 2D Histogram with the p-Values in the dimensions (scanpoint, BergerBoos_id) More...

int	getNBergerBoosPointsPerScanpoint ()

	MethodBergerBoosScan (MethodProbScan *s, TString d="XX")
	Initialize from a previous Prob scan, setting the profile likelihood. More...

void	readScan1dTrees (int runMin=1, int runMax=1)
	Return number of BB points per scan point. More...

int	scan1d (int nRun=1)
	Perform the 1d BergerBoos scan. More...

void	setNBergerBoosPointsPerScanpoint (int n)

void	setNewBergerBoosPoint (int m)
	Set number of BB points per scan point. More...

	~MethodBergerBoosScan ()

Public Member Functions inherited from MethodPluginScan
int	getNtoys ()

double	getPvalue1d (RooSlimFitResult plhScan, double chi2minGlobal, ToyTree t=0, int id=0)

	MethodPluginScan (MethodProbScan *s)
	Initialize from a previous Prob scan, setting the profile likelihood. More...

	MethodPluginScan (MethodProbScan s, PDF_Datasets pdf, OptParser *opt)
	Constructor, mainly to ensure compatibility with MethodDatasetsPluginScan. More...

	MethodPluginScan (Combiner *comb)
	Initialize from a Combiner object. More...

virtual void	readScan1dTrees (int runMin=1, int runMax=1, TString fName="default")
	Read in the TTrees that were produced by scan1d(). More...

void	readScan2dTrees (int runMin=1, int runMax=1)
	Read in the TTrees that were produced by scan2d(). More...

virtual void	scan2d (int nRun=1)
	Perform the 2d Plugin scan. More...

void	setNtoysPerPoint (int n)

void	setParevolPLH (MethodProbScan *s)
	Set the parameter evolution over the profile likelihood that was previously computed by a MethodProbScan scanner. More...

Public Member Functions inherited from MethodAbsScan
virtual void	calcCLintervals (int CLsType=0)
	Calculate the CL intervals from the CL curve. More...

void	calcCLintervalsSimple (int CLsType=0)

virtual void	checkCLs ()

void	confirmSolutions ()
	Refit all possible solutions with the scan parameter left free to confirm the solutions. More...

void	doInitialFit (bool force=false)
	Try to find global mininum of the PDF. More...

const vector< RooSlimFitResult * > &	getAllResults ()

OptParser *	getArg ()

const std::pair< double, double >	getBorders (const TGraph &graph, const double confidence_level, bool qubic=false)
	determines the borders of the confidence interval by linear or qubic interpolation. More...

const std::pair< double, double >	getBorders_CLs (const TGraph &graph, const double confidence_level, bool qubic=false)

float	getChi2minBkg ()

float	getChi2minGlobal ()

float	getCL (double val)

CLInterval	getCLinterval (int iSol=0, int sigma=1)
	Get the CL interval that includes the best-fit value. More...

CLInterval	getCLintervalCentral (int sigma=1)
	Get the CL interval that includes the best-fit value. More...

Combiner *	getCombiner () const

const vector< RooSlimFitResult * > &	getCurveResults ()

int	getDrawSolution ()

bool	getFilled ()

int	getFillStyle ()

TH1F *	getHchisq ()

TH2F *	getHchisq2d ()

TH1F *	getHCL ()

TH2F *	getHCL2d ()

TH1F *	getHCLs ()

TH2F *	getHCLs2d ()

TH1F *	getHCLsErr1Dn ()

TH1F *	getHCLsErr1Up ()

TH1F *	getHCLsErr2Dn ()

TH1F *	getHCLsErr2Up ()

TH1F *	getHCLsExp ()

TH1F *	getHCLsFreq ()

int	getLineColor ()

int	getLineStyle ()

TString	getMethodName () const

TString	getName () const

int	getNObservables ()

int	getNPoints1d ()

int	getNPoints2dx ()

int	getNPoints2dy ()

const RooArgSet *	getObservables ()

TString	getObsName ()

TString	getParsName ()

RooRealVar *	getScanVar1 ()

TString	getScanVar1Name ()

float	getScanVar1Solution (int i=0)
	Get value of scan parameter 1 a certain solution. More...

RooRealVar *	getScanVar2 ()

TString	getScanVar2Name ()

float	getScanVar2Solution (int i=0)
	Get value of scan parameter 2 a certain solution (only meaningful for 2d scan). More...

float	getScanVarSolution (int iVar, int iSol)
	Get value of scan parameter at a certain solution. More...

RooSlimFitResult *	getSolution (int i=0)
	Return a solution corresponding to a minimum of the profile likelihoood. More...

vector< RooSlimFitResult * >	getSolutions ()

int	getTextColor ()

const RooArgSet *	getTheory ()

TString	getTitle ()

RooWorkspace *	getWorkspace ()

virtual void	initScan ()

void	loadParameters (RooSlimFitResult *r)
	Load the values given by an (external) fit result. More...

virtual bool	loadScanner (TString fName="")
	Save a scanner from plots/scanner. More...

bool	loadSolution (int i=0)
	Load the values at a specific minimum into the workspace. More...

	MethodAbsScan ()

	MethodAbsScan (Combiner *c)

	MethodAbsScan (OptParser *opt)

void	plot1d (TString var)
	Make a 1d plot of the NLL in var. More...

void	plot2d (TString varx, TString vary)
	Make a 2d plot of the PDF in varx and vary. More...

void	plotOn (OneMinusClPlotAbs *plot, int CLsType=0)

void	plotPulls (int nSolution=0)
	Make a pull plot of observables corresponding to the given solution. More...

virtual void	print ()

void	printCLintervals (int CLsType)
	Print CL intervals. More...

void	printLocalMinima ()
	Print local minima solutions. More...

void	saveLocalMinima (TString fName="")
	Save local minima solutions. More...

void	saveScanner (TString fName="")
	Save this scanner to a root file placed into plots/scanner. More...

virtual int	scan1d ()

virtual int	scan2d ()

void	setChi2minGlobal (double x)
	Set the global minimum manually. More...

void	setDrawSolution (int code=0)

void	setFilled (bool filled)

void	setFillStyle (int c)

void	setHchisq (TH1F *h)

void	setHCL (TH1F *h)

void	setLineColor (int c)

void	setLineStyle (int c)

void	setNPoints1d (int n)

void	setNPoints2dx (int n)

void	setNPoints2dy (int n)

void	setPValueCorrector (PValueCorrection *pvalCor)

void	setScanVar1 (TString var)

void	setScanVar2 (TString var)

void	setSolutions (vector< RooSlimFitResult *> s)
	Helper function to copy over solutions from another scanner. More...

void	setTextColor (int c)

void	setTitle (TString s)

void	setVerbose (bool yesNo=true)

void	setXscanRange (float min, float max)

void	setYscanRange (float min, float max)

	~MethodAbsScan ()

Public Attributes
TTree *	BBtree

TString	dir

TFile *	file

Public Attributes inherited from MethodAbsScan
vector< RooSlimFitResult * >	allResults
	All fit results we encounter along the scan. More...

vector< CLInterval >	clintervals1sigma
	all 1 sigma intervals that were found by calcCLintervals() More...

vector< CLInterval >	clintervals2sigma
	all 2 sigma intervals that were found by calcCLintervals() More...

vector< CLInterval >	clintervals3sigma
	all 3 sigma intervals that were found by calcCLintervals() More...

vector< RooSlimFitResult * >	curveResults
	All fit results of the the points that make it into the 1-CL curve. More...

vector< vector< RooSlimFitResult * > >	curveResults2d
	All fit results of the the points that make it into the 1-CL curve. More...

vector< RooSlimFitResult * >	solutions
	Local minima filled by saveSolutions() and saveSolutions2d(). More...

Protected Attributes
int	nBBPoints
	number of sampled Berger Boos points per scan point More...

Protected Attributes inherited from MethodPluginScan
int	nToys
	number of toys to be generated at each scan point More...

MethodProbScan *	parevolPLH
	external scanner defining the parameter evolution: set to profileLH unless for the Hybrid Plugin More...

MethodProbScan *	profileLH
	external scanner holding the profile likelihood: DeltaChi2 of the scan PDF on data More...

Protected Attributes inherited from MethodAbsScan
OptParser *	arg
	command line options More...

double	chi2minBkg
	chi2 value at global minimum More...

double	chi2minGlobal
	chi2 value at global minimum More...

bool	chi2minGlobalFound
	flag to avoid finding minimum twice More...

Combiner *	combiner
	the combination More...

bool	drawFilled
	choose if Histogram is drawn filled or not More...

int	drawSolution
	Configure how to draw solutions on the plots. More...

int	fillStyle

RooFitResult *	globalMin
	parameter values at a global minimum More...

TH1F *	hChi2min
	histogram for the chi2min values before Prob() More...

TH2F *	hChi2min2d
	histogram for the chi2min values before Prob() More...

TH1F *	hCL
	1-CL curve More...

TH2F *	hCL2d
	1-CL curve More...

TH1F *	hCLs
	1-CL curve More...

TH2F *	hCLs2d
	1-CL curve More...

TH1F *	hCLsErr1Dn
	1-CL curve More...

TH1F *	hCLsErr1Up
	1-CL curve More...

TH1F *	hCLsErr2Dn
	1-CL curve More...

TH1F *	hCLsErr2Up
	1-CL curve More...

TH1F *	hCLsExp
	1-CL curve More...

TH1F *	hCLsFreq
	1-CL curve More...

int	lineColor

int	lineStyle

bool	m_initialized
	true if initScan() was called More...

bool	m_xrangeset
	true if the x range was set manually (setXscanRange()) More...

bool	m_yrangeset
	true if the y range was set manually (setYscanRange()) More...

TString	methodName
	Prob, ... More...

TString	name
	basename, e.g. ggsz More...

int	nPoints1d
	number of scan points used by 1d scan More...

int	nPoints2dx
	number of scan points used by 2d scan, x axis More...

int	nPoints2dy
	number of scan points used by 2d scan, y axis More...

int	nWarnings
	number of warnings printed in getScanVarSolution() More...

RooDataSet *	obsDataset
	save the nominal observables so we can restore them after we have fitted toys More...

TString	obsName
	dataset name of observables, derived from name More...

TString	parsName
	set name of physics parameters, derived from name More...

TString	pdfName
	PDF name in workspace, derived from name. More...

PValueCorrection *	pvalueCorrector

bool	pvalueCorrectorSet

TRandom3	rndm

TString	scanVar1
	scan parameter More...

TString	scanVar2
	second scan parameter if we're scanning 2d More...

RooDataSet *	startPars
	save the start parameter values before any scan More...

int	textColor
	color used for plotted central values More...

TString	thName
	set name of theory parameters, derived from name More...

TString	title
	nice string for the legends More...

TString	toysName
	set name of parameters to vary in toys More...

bool	verbose

RooWorkspace *	w

Additional Inherited Members
Protected Member Functions inherited from MethodPluginScan
TH1F *	analyseToys (ToyTree *t, int id=-1)
	Analyse toys that were written either by a scan or a by getPvalue(). More...

void	computePvalue1d (RooSlimFitResult plhScan, double chi2minGlobal, ToyTree t, int id, Fitter f, ProgressBar pb)
	Compute the p-value at a certain point in parameter space using the plugin method. More...

RooDataSet *	generateToys (int nToys)
	Generate toys. More...

RooSlimFitResult *	getParevolPoint (float scanpoint)
	Helper function for scan1d(). More...

double	importance (double pvalue)
	Importance sampling for low p-values: Returns a value between 0.05 and 1 which can be used to scale down the number of toys to be generated at each scan step. More...

Protected Member Functions inherited from MethodAbsScan
void	sortSolutions ()
	Sort solutions in order of increasing chi2. More...

Detailed Description

Definition at line 13 of file MethodBergerBoosScan.h.

Constructor & Destructor Documentation

§ MethodBergerBoosScan()

MethodBergerBoosScan::MethodBergerBoosScan	(	MethodProbScan *	s,
		TString	d = `"XX"`
	)

Initialize from a previous Prob scan, setting the profile likelihood.

This should be the default.

Definition at line 14 of file MethodBergerBoosScan.cpp.

§ ~MethodBergerBoosScan()

MethodBergerBoosScan::~MethodBergerBoosScan ( )

Definition at line 42 of file MethodBergerBoosScan.cpp.

Member Function Documentation

§ calcPValues()

TH2F * MethodBergerBoosScan::calcPValues	(	TH2F	better,
		TH2F	all,
		TH2F	bg
	)

Calculates 2D p-Values in (scanpoint,BergerBoos_id) Histogram.

Definition at line 50 of file MethodBergerBoosScan.cpp.

§ drawBBPoints()

void MethodBergerBoosScan::drawBBPoints	(	TString	varX,
		TString	varY,
		int	runMin = `1`,
		int	runMax = `1`,
		bool	save = `true`
	)

Draws 2D Histogram showing the BB points in varX-varY space.

Draws a 2D Histogram in var1:var2 space.

The boolian 'save' specifies if a copy of the plot will be saved in the plots folder.

Definition at line 78 of file MethodBergerBoosScan.cpp.

§ getBestPValue()

void MethodBergerBoosScan::getBestPValue	(	TH1F *	hCL,
		TH2F *	pValues
	)

Returns a 1D Histogram containing the 1-CL curve Input is a 2D Histogram with the p-Values in the dimensions (scanpoint, BergerBoos_id)

Definition at line 154 of file MethodBergerBoosScan.cpp.

§ getNBergerBoosPointsPerScanpoint()

int MethodBergerBoosScan::getNBergerBoosPointsPerScanpoint ( )

inline

Definition at line 20 of file MethodBergerBoosScan.h.

§ readScan1dTrees()

void MethodBergerBoosScan::readScan1dTrees	(	int	runMin = `1`,
		int	runMax = `1`
	)

Return number of BB points per scan point.

Reads in Scan1DTree result from scan1D() This is a slightly modified version of MethodPluginScan::readScan1dTrees()'s function to account for BergerBoos specifications.

Definition at line 179 of file MethodBergerBoosScan.cpp.

§ scan1d()

int MethodBergerBoosScan::scan1d ( int nRun = 1 )

virtual

Perform the 1d BergerBoos scan.

Saves chi2 values in a root tree, together with the full fit result for each toy. If a combined PDF for the toy generation is given by setToyGenCombo(), this will be used to generate the toys.

Parameters

nRun	Part of the root tree file name to facilitate parallel production.

Reimplemented from MethodPluginScan.

Definition at line 329 of file MethodBergerBoosScan.cpp.

§ setNBergerBoosPointsPerScanpoint()

void MethodBergerBoosScan::setNBergerBoosPointsPerScanpoint ( int n )

inline

Definition at line 23 of file MethodBergerBoosScan.h.

§ setNewBergerBoosPoint()

void MethodBergerBoosScan::setNewBergerBoosPoint ( int m )

Set number of BB points per scan point.

Set number of Berger Boos points drawn from the BergerBoos CL intervals defined in the workspace Samples one Berger Boos point Draws ALL parameters randomly according to their BergerBoos ranges defined in the workspace The parameter of interest (scan parameter) has to be set constant (and to its scan point value) after the function call.

Definition at line 496 of file MethodBergerBoosScan.cpp.

Member Data Documentation

§ BBtree

TTree* MethodBergerBoosScan::BBtree

Definition at line 35 of file MethodBergerBoosScan.h.

§ dir

TString MethodBergerBoosScan::dir

Definition at line 36 of file MethodBergerBoosScan.h.

§ file

TFile* MethodBergerBoosScan::file

Definition at line 34 of file MethodBergerBoosScan.h.

§ nBBPoints

int MethodBergerBoosScan::nBBPoints

protected

number of sampled Berger Boos points per scan point

Definition at line 40 of file MethodBergerBoosScan.h.

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

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

Public Member Functions

Public Attributes

Protected Attributes

Additional Inherited Members