XeTest Class Reference

Test and Demonstration package. More...

#include <XeTest.h>

Inheritance diagram for XeTest:

Public Member Functions
void	computeSensitivityBands (bool s1=true, bool publishedBkg=true, int nSimulations=100, int runNumber=10)
	Compute Sensitivity bands of sigma.
void	PLexclusionComputedBackground (bool withCLs=true, double sigToE=1.)
	Compute the CLs limits making poisson a maximum likekihood.
void	plotAcceptance ()
	Plot cut acceptances, for both smeared and unsmeared cuts.
void	plotFFSD ()
	Thorough comparison of FF and Structure functions for SD parametrizations.
void	plotLimitsSD (int runNumber=10)
	Plot theoretical cross section limits for various SD parametrizations.
void	plotLimitsSI (int runNumber)
	Plots the expected number of events for Xenon100(r) Run 10 and the theoretical exclusion limit in case of no candidates.
void	plotLogLL (bool s1=true, bool publishedBkg=true, bool doFitSystBkTvalue=false, double specialMass=0., int nMasses=0, int runNumber=10)
	Plot log likelihood as a function of number of events.
void	plotS1SD (double mass=50.*GeV, int runNumber=10)
	Plots S1 and Er distributions for various models of SD.
void	plotS1SI (double mass=50.*GeV, int runNumber=10)
	Plots S1 and Er distributions for SI, various LEff t-values.
void	plotSigmaSD (double mass=50.GeV, double sig=1.E-45 cm2)
	Plot Energy recoil for various SD parametrizations.
void	plotSigmaSI (double mass=50.*GeV, double sigma=2.E-45, int A=NATURAL)
	Plot rate for Xenon target, spin independent.
bool	printTheSetup ()
	Print the result of the setup.
void	runAnalysis (bool s1=true, bool publishedBkg=true, bool doFitSystBkTvalue=false, double specialMass=0., int runNumber=10)
	Compute limit of sigma SI for runs 8 and 10 Profile likelihood.
void	setupTheAnalysis (int mode=PL_ANALYSIS, bool forSigmaPositive=true, bool doFitSystBkTvalue=true, bool doFitLeffTValue=true, bool withS1Shape=true)
	Technical method for interactive operation: setup an analysis.
void	setupTheRun (bool publishedBackground=true, double mass=DEFAULT_WIMP_MASS, int runNumber=10)
	technical method for interactive operation: setup a run
void	showRunData (bool publishedEr=true, double mass=DEFAULT_WIMP_MASS, int runNumber=10, int analysisMode=PL_ANALYSIS, int displayMode=FLATTENED_S2_VS_S1)
	Show data, S1S2 bands for signal and background.
void	studyPoissonCoverage ()
	Compute the coverage of Poisson confidence levels.
void	studyQValuesForBackground (int nm=0, int nSimulations=100000, bool s1=true, bool publishedBkg=true, int runNumber=10)
	Show how far we are from Wilk's theorem.
void	tableCLs (bool withCLs=true)
	Print a table of simple CLs exclusoin.
void	tableCombinedCLs (bool withCLs=true)
	Print a table of combined CLs thru Fisher method.
void	tableCombinedPL ()
	Print a table of combined PL.
void	testExclusionPL (int sMin=0, int sMax=10)
	shows how the Exclusion class runs on simple Signal+background PL
void	testLK ()
	ML fit of the parameters of a distribution made of an exponential distribution plus a flat background.
void	testYellin (int s=5, int b=0, int n=100)
	Demonstrates Yellin's method.
Public Attributes
S1S2Bands *	bands
XeSetOfCuts *	cuts
Exclusion *	exclusion
XeGraph *	gr
SDInteraction *	interSD
SIInteraction *	interSI
XeMultiGraph *	mg
S1S2PL *	pl
PValue *	pv
XeRun *	run
OneDimSimulatedDataSet *	simData
Target *	target
Wimp *	wimp

Detailed Description

Test and Demonstration package.

Member Function Documentation

void XeTest::computeSensitivityBands	(	bool	s1 = `true`,
		bool	publishedBkg = `true`,
		int	nSimulations = `100`,
		int	runNumber = `10`
	)

Compute Sensitivity bands of sigma.

Parameters:

	s1	with s1 profile likelihood
	publishedBkg	use PRL background instead of computed one
	nSimulations	number of runs per mass points
	runNumber	Where to take data from

void XeTest::PLexclusionComputedBackground	(	bool	withCLs = `true`,
		double	sigToE = `1.`
	)

Compute the CLs limits making poisson a maximum likekihood.

Parameters:

	withCls	Apply CLs correction?
	sigToE	Conversion from cm2 to number of events

void XeTest::plotLimitsSD ( int runNumber = 10 )

Plot theoretical cross section limits for various SD parametrizations.

Parameters:

runNumber

Run number to be mimicked

void XeTest::plotLimitsSI ( int runNumber )

Plots the expected number of events for Xenon100(r) Run 10 and the theoretical exclusion limit in case of no candidates.

Parameters:

runNumber

Run number to be mimicked

void XeTest::plotLogLL	(	bool	s1 = `true`,
		bool	publishedBkg = `true`,
		bool	doFitSystBkTvalue = `false`,
		double	specialMass = `0.`,
		int	nMasses = `0`,
		int	runNumber = `10`
	)

Plot log likelihood as a function of number of events.

Parameters:

	s1	with s1 profile likelihood
	publishedBkg	use PRL background instead of computed one
	doFitSystBkTvalue	do we also fit the syst unc. on background?
	specialMass	if >0., do it for one mass only nMasses number of masses to be drawns
	runNumber	Where to take data from

void XeTest::plotS1SD	(	double	mass = `50.*GeV`,
		int	runNumber = `10`
	)

Plots S1 and Er distributions for various models of SD.

Parameters:

	mass	Wimp mass
	runNumber	Run number to be mimicked

void XeTest::plotS1SI	(	double	mass = `50.*GeV`,
		int	runNumber = `10`
	)

Plots S1 and Er distributions for SI, various LEff t-values.

Parameters:

	mass	Wimp mass
	runNumber	Run number to be mimicked

void XeTest::plotSigmaSD	(	double	mass = `50.*GeV`,
		double	sig = `1.E-45 *cm2`
	)

Plot Energy recoil for various SD parametrizations.

Parameters:

	mass	Wimp mass
	sigma	wimp-nucleon cross section

void XeTest::plotSigmaSI	(	double	mass = `50.*GeV`,
		double	sigma = `2.E-45`,
		int	A = `NATURAL`
	)

Plot rate for Xenon target, spin independent.

Parameters:

	mass	Wimp mass
	sigma	wimp-nucleon cross section
	A	NATURAL, DEPLETED, or nuclear number

bool XeTest::printTheSetup ( )

Print the result of the setup.

Returns:: ok or not

void XeTest::runAnalysis	(	bool	s1 = `true`,
		bool	publishedBkg = `true`,
		bool	doFitSystBkTvalue = `false`,
		double	specialMass = `0.`,
		int	runNumber = `10`
	)

Compute limit of sigma SI for runs 8 and 10 Profile likelihood.

Parameters:

	s1	with s1 profile likelihood
	publishedBkg	use PRL background instead of computed one
	doFitSystBkTvalue	do we also fit the syst unc. on background?
	specialMass	if >0., print detailed computation for this mass
	runNumber	Where to take data from

void XeTest::setupTheAnalysis	(	int	mode = `PL_ANALYSIS`,
		bool	forSigmaPositive = `true`,
		bool	doFitSystBkTvalue = `true`,
		bool	doFitLeffTValue = `true`,
		bool	withS1Shape = `true`
	)

Technical method for interactive operation: setup an analysis.

Parameters:

	mode	either PL_ANALYSIS or CUTS_ANALYSIS
	forSigmaPositive	Do we force sigmaHat to be positive?
	doFitSystBkTvalue	do we also fit the syst unc. on background?
	doFitLeffTValue	Do we want to fit LEFF t-value?
	withS1Shape	consider S1 shape?

void XeTest::setupTheRun	(	bool	publishedBackground = `true`,
		double	mass = `DEFAULT_WIMP_MASS`,
		int	runNumber = `10`
	)

technical method for interactive operation: setup a run

Parameters:

	publishedBackground	use PRL background instead of computed one
	mass	Wimp Mass for initial print-out
	runNumber	Where to take data from

void XeTest::showRunData	(	bool	publishedEr = `true`,
		double	mass = `DEFAULT_WIMP_MASS`,
		int	runNumber = `10`,
		int	analysisMode = `PL_ANALYSIS`,
		int	displayMode = `FLATTENED_S2_VS_S1`
	)

Show data, S1S2 bands for signal and background.

Parameters:

	mass	Wimp Mass for initial print-out
	runNumber	Where to take data from
	analysisMode	PL_ANALYSIS or CUTS_ANALYSIS
	displayMode	S2_VS_S1, _OVER_S1_VS_S1, or FLATTENED_S2_VS_S1

void XeTest::studyQValuesForBackground	(	int	nm = `0`,
		int	nSimulations = `100000`,
		bool	s1 = `true`,
		bool	publishedBkg = `true`,
		int	runNumber = `10`
	)

Show how far we are from Wilk's theorem.

Parameters:

	number	of masses to study (default: all of them)
	nSimulations	number of runs per mass points
	s1	with s1 profile likelihood
	publishedBkg	use PRL background instead of computed one
	runNumber	Where to take data from