Statistical model for CMS-SMP-24-002

Introduction

This repository provides the statistical model for the CMS Zgamma cross section measurement and nTGC constraints at 13.6 TeV analysis on 2022 data.

The instructions given below detail the commands needed to reproduce the main results of the paper. For additional documentation on Combine, refer to the CMS Combine Documentation.

Analysis specific instructions

For the commands below, you may require running ulimit -s unlimited; ulimit -u unlimited to avoid memory issues. All commands are run with Combine v10, within CMSSW_11_3_4. All systematic naming conventions are given in `input/systematics.yml` . 

Datacards and commands

The `input` directory contains the datacard for this analysis. Below we provide the commands use to extract the best fit values and the confidence intervals, splitting the total uncertainty into statistical and systematic components. 

cd input/
text2workspace.py combinedcard.txt -m 125
combine  combinedcard.root -M MultiDimFit --saveWorkspace -n combinedcard.postfit --expectSignal=1 -t -1
combine higgsCombinecombinedcard.postfit.MultiDimFit.mH120.root -M MultiDimFit -n combinedcard.total --algo grid --snapshotName MultiDimFit --setParameterRanges r=0,4 --expectSignal=1 -t -1
combine higgsCombinecombinedcard.postfit.MultiDimFit.mH120.root -M MultiDimFit --algo grid --snapshotName MultiDimFit --setParameterRanges r=0,4  --freezeNuisanceGroups theory -n combinedcard.freeze_theory --expectSignal=1 -t -1
combine higgsCombinecombinedcard.postfit.MultiDimFit.mH120.root -M MultiDimFit --algo grid --snapshotName MultiDimFit --setParameterRanges r=0,4  --freezeNuisanceGroups theory,sys -n combinedcard.freeze_theory_sys --expectSignal=1 -t -1
plot1DScan.py higgsCombinecombinedcard.total.MultiDimFit.mH120.root --main-label "Total Uncert."  --others higgsCombinecombinedcard.freeze_theory.MultiDimFit.mH120.root:"freeze theory":8 higgsCombinecombinedcard.freeze_theory_sys.MultiDimFit.mH120.root:"stat only":6  --output breakdown_combine --y-max 10 --y-cut 40 --breakdown "theory,sys,stat"