StatTools interface for making plots and datacards
cmsrel CMSSW_7_4_7
cd CMSSW_7_4_7/src/
cmsenv
git cms-init
git clone https://github.com/cms-analysis/HiggsAnalysis-CombinedLimit.git HiggsAnalysis/CombinedLimit
git clone https://github.com/cms-analysis/CombineHarvester.git CombineHarvester #combine harvester integration is underway
mkdir UWAnalysis
git clone [email protected]:lmdodd/StatTools.git UWAnalysis/StatTools
Combine hard to set up. Always make a clean build, as scram doesn't always see updates to src/LinkDef.h. See https://twiki.cern.ch/twiki/bin/viewauth/CMS/SWGuideHiggsAnalysisCombinedLimit for more
cd $CMSSW_BASE/src/HiggsAnalysis/CombinedLimit
git fetch origin
git checkout v6.2.1
scramv1 b clean; scramv1 b
To build the package
cd $CMSSW_BASE/src/UWAnalysis/StatTools/
export USER_CXXFLAGS="-Wno-error=unused-but-set-variable -Wno-error=sign-compare"
source setup.sh
scram b -j8
to make plots make a folder in data
cd UWAnalysis/StatTools/data
cd sm_higgs_ltau
./makeMuTauPlots
root -l -b -q makeTemplatePlotsMuTau.C #check the output directory here
to run limit example: remember to change 'aux_shapes' in cpp file. works for XTT <-> HTT
(remember to change the aux shapes directory in HTT)
HTT
cd output/htt_cards
combineTool.py -M T2W -i {cmb,et,mt}/* -o workspace.root --parallel 4
nohup combineTool.py -M Asymptotic -d */*/workspace.root --there -n .limit --parallel 4 &
combineTool.py -M CollectLimits */*/*.limit.* --use-dirs -o limits.json
plotLimits.py limits_{et,mt,cmb}.json:exp0 --auto-style --logy # will plot the expected limit comparisons of the channels