CMSDAS-MonoZ-Tutorial

Quick start with combine:

Setting up CMSSW:

export SCRAM_ARCH=slc7_amd64_gcc700
cmsrel CMSSW_10_2_13 
cd CMSSW_10_2_13/src 
cmsenv

fetching combine:

git clone https://github.com/cms-analysis/HiggsAnalysis-CombinedLimit.git HiggsAnalysis/CombinedLimit 
cd $CMSSW_BASE/src/HiggsAnalysis/CombinedLimit 
git fetch origin 
git checkout v8.0.1

# and combinetool
cd $CMSSW_BASE/src 
bash <(curl -s https://raw.githubusercontent.com/cms-analysis/CombineHarvester/master/CombineTools/scripts/sparse-checkout-https.sh) 

# compile everything
cd $CMSSW_BASE/src/ 
scramv1 b clean; scramv1 b -j 10

Fetching a Mono-Z datacard:

cp -r /eos/user/c/cmsdas/long-exercises/MonoZ/datacards/cards-DMSimp_MonoZLL_NLO_Axial_1000_MXd-1 .

# remove the old combine output
rm cards-DMSimp_MonoZLL_NLO_Axial_1000_MXd-1/combined.root

Transform the datacard into a worspace format (this step can be skipped if you run combine on combined.dat directly, see below)

text2workspace.py cards-DMSimp_MonoZLL_NLO_Axial_1000_MXd-1/combined.dat -o cards-DMSimp_MonoZLL_NLO_Axial_1000_MXd-1/combined.root

run combine

combine  -M AsymptoticLimits --datacard cards-DMSimp_MonoZLL_NLO_Axial_1000_MXd-1/combined.root

Making impacts plot

export PARAM="--rMin=-1 --rMax=4 --cminFallbackAlgo Minuit2,Migrad,0:0.05 --X-rtd MINIMIZER_analytic --X-rtd FAST_VERTICAL_MORPH"

combineTool.py -M Impacts -d cards-DMSimp_MonoZLL_NLO_Axial_1000_MXd-1/combined.root -m 125 --robustFit 1 --doInitialFit $PARAM
combineTool.py -M Impacts -d cards-DMSimp_MonoZLL_NLO_Axial_1000_MXd-1/combined.root -m 125 --robustFit 1 --doFits $PARAM
combineTool.py -M Impacts -d cards-DMSimp_MonoZLL_NLO_Axial_1000_MXd-1/combined.root -m 125 -o impactsTEST.json $PARAM

# plot everything now:
plotImpacts.py -i impactsTEST.json -o impactsTEST