quickstats 0.8.2.0

A tool for quick statistical analysis for HEP experiments

Setup

Clone the repository:

git clone ssh://[e-mail]:7999/clcheng/quickstats.git

1. CERN User

To set up from lxplus, just do

source setup.sh

2. Genearl User

To set up locally, make sure you have pyROOT 6.24+ installed (using conda is recommended), and do

pip install quickstats

Installing pyROOT

Simplest way to install pyROOT is via conda

conda install -c conda-forge ROOT

Important: First-time compilation

To compile c++ dependencies, do this for first time use

quickstats compile

Command Line Tools

Run Nuisance Parameter Pulls and Ranking

quickstats run_pulls -i <input_ws_path> -d <dataset_name> -p <np_name/pattern> --poi <poi_name> --parallel -1 -o <output_dir>

The following options are available

Option	Description	Default
-i/--input_file	Path to the input workspace file	-
-w/--workspace	Name of workspace. Auto-detect by default.	None
-m/--model_config	Name of model config. Auto-detect by default.	None
-d/--data	Name of dataset	"combData"
-p/--parameter	Nuisance parameter(s) to run pulls on. Multiple parameters are separated by commas. Wildcards are accepted. All NPs will be run over by default	""
-x/--poi	POIs to measure. If empty, impact on POI will not be calculated.	""
-r/--profile	Parameters to profile	""
-f/--fix	Parameters to fix	""
-s/--snapshot	Name of initial snapshot	"nominalNuis"
-o/--outdir	Output directory	"pulls"
-t/--minimizer_type	Minimizer type	"Minuit2"
-a/--minimizer_algo	Minimizer algorithm	"Migrad"
-c/--num_cpu	Number of CPUs to use per parameter	1
--binned/--unbinned	Whether to use binned likelihood	True
-q/--precision	Precision for scan	0.001
-e/--eps	Tolerance	1.0
-l/--log_level	Log level	"INFO"
--eigen/--no-eigen	Compute eigenvalues and vectors	False
--strategy	Default fit strategy	0
--fix-cache/--no-fix-cache	Fix StarMomentMorph cache	True
--fix-multi/--no-fix-multi	Fix MultiPdf level 2	True
--offset/--no-offset	Offset likelihood	True
--optimize/--no-optimize	Optimize constant terms	True
--max_calls	Maximum number of function calls	-1
--max_iters	Maximum number of Minuit iterations	-1
--parallel	Parallelize job across different nuisanceparameters using N workers. Use -1 for N_CPU workers.	0
--cache/--no-cache	Cache existing result	True
--exclude	Exclude NPs (wildcard is accepted)	""

Plot Nuisance Parameter Pulls and Ranking

quickstats plot_pulls --help

Likelihood Fit (Best-fit)

quickstats likelihood_fit --help

Run Likelihood Scan

quickstats likelihood_scan --help

Asymptotic CLs Limit

quickstats cls_limit --help

CLs Limit Scan

quickstats limit_scan --help

Generate Asimov dataset

quickstats generate_standard_asimov --help

Inspect Workspace

quickstats inspect_workspace --help

Create Workspace from XML Cards

quickstats build_xml_ws --help

Modify Workspace from XML Cards or Json Config

quickstats modify_ws --help

Combine Workspace from XML Cards or Json Config

quickstats combine_ws --help

Compare Workspaces

quickstats compare_ws --help

Run Event Loop from Custom Config File

quickstats process_rfile --help