The geometric acceptance multiplied by the efficiency of the $\it{p}_{T}^\text{miss} >$ 200 GeV selection as a function of the proper decay length c$\tau$ for a scalar particle S with a mass of 40 GeV. The acceptance region for DT is defined by requiring the LLP to decay in the region with |z| < 661 cm and 380 cm < r < 736 cm. The acceptance region for CSC is defined by requiring the LLP decay in the region with $|\eta| < 2.4$, r < 695.5 cm, and 661 cm < |z| < 1100 cm or in the region with $|\eta| < 2.4$, r < 270 cm, and 500 cm < |z| < 661 cm. Single CSC cluster requires exactly one LLP to decay in CSC; Single DT cluster requires exactly one LLP to decay in DT; Double cluster requires both LLP to decay in CSC or DT. The denominator in this plot includes all generated events. The nominator includes events that pass the acceptance requirements above and $\it{p}_{T}^\text{miss} >$ 200 GeV.

Distributions of the cluster time for signal, where S decaying to $d\bar{d}$ for a proper decay length c$\tau$ of 1 m and mass of 40 GeV, and for a background-enriched sample in data selected by inverting the $N_\text{hits}$ requirement.

The distributions of $N_\text{hits}$ for single CSC clusters are shown for signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m).

The distributions of $N_\text{hits}$ for single CSC clusters are shown for compared to the OOT background ($t_\text{clusters} < 12.5$ ns). The OOT background is representative of the overall background shape, because the background passing all the selections described above is dominated by pileup and underlying events.

The distributions of $\Delta\phi(\it{p}_{T}^\text{miss} \text{,cluster)}$ for single CSC clusters are shown for signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m), compared to the shape of background in a selection in which the cluster is not matched to any RPC hit and passes all other selections. The background is dominated by clusters from noise and low-pT particles.

The distributions of $N_\text{hits}$ for single CSC clusters are shown for signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m), compared to the shape of background in a selection in which the cluster is not matched to any RPC hit and passes all other selections. The background is dominated by clusters from noise and low-pT particles.

The distributions of $\Delta\phi(\it{p}_{T}^\text{miss} \text{,cluster)}$ for single CSC clusters are shown for signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m), compared to the shape of background in a selection in which the cluster is not matched to any RPC hit and passes all other selections. The background is dominated by clusters from noise and low-pT particles.

The signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m), background, and data distributions of $N_\text{clusters}$ passing the $N_\text{hits}$ selection in the search region for CSC-CSC category.

The signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m), background, and data distributions of $N_\text{clusters}$ passing the $N_\text{hits}$ selection in the search region for DT-DT category.

The signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m), background, and data distributions of $N_\text{clusters}$ passing the $N_\text{hits}$ selection in the search region for DT-CSC category.

The signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m), background, and data distributions of $N_\text{hits}$ in the search region of the single-CSC cluster category are shown. The $N_\text{hits}$ distribution includes only events in bins A and D. The right-hand bin in the Nhits distribution includes overflow events.

The signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m), background, and data distributions of $\Delta\phi(\it{p}_{T}^\text{miss} \text{,cluster)}$ in the search region of the single-CSC cluster category are shown. The $\Delta\phi(\it{p}_{T}^\text{miss} \text{,cluster)}$ distribution includes only events in bins A and B.

The signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m), background, and data distributions of $N_\text{hits}$ in the search region of the single-DT cluster category are shown. The $N_\text{hits}$ distribution includes only events in bins A and D. The right-hand bin in the Nhits distribution includes overflow events.

The signal (assuming B(H $\rightarrow$ SS) = 1%, S $\rightarrow d\bar{d}$, and c$\tau$ = 1 m), background, and data distributions of $\Delta\phi(\it{p}_{T}^\text{miss} \text{,cluster)}$ in the search region of the single-DT cluster category are shown. The $\Delta\phi(\it{p}_{T}^\text{miss} \text{,cluster)}$ distribution includes only events in bins A and B.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 3 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 7 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 0.4 GeV mass and $ S \rightarrow \pi^{0} \pi^{0}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 1 GeV mass and $ S \rightarrow \pi^{0} \pi^{0}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 7 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 0.4 GeV mass and $ S \rightarrow \pi^{+} \pi^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 1 GeV mass and $ S \rightarrow \pi^{+} \pi^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 1.5 GeV mass and $ S \rightarrow K^{+}K^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 1.5 GeV mass and $ S \rightarrow K^{0}K^{0}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 0.4 GeV mass and $ S \rightarrow \gamma\gamma$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 0.4 GeV mass and $ S \rightarrow \e^{+} \e^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) as functions of mass and c$\tau$ assuming S inherits all couplings from the Higgs boson evaluated at the LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) as functions of mass and c$\tau$ assuming S inherits all couplings from the Higgs boson evaluated at the LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) as functions of mass and c$\tau$ assuming S inherits all couplings from the Higgs boson evaluated at the LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) as functions of mass and c$\tau$ assuming S inherits all couplings from the Higgs boson evaluated at the LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) as functions of mass and c$\tau$ assuming S inherits all couplings from the Higgs boson evaluated at the LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) as functions of mass and c$\tau$ assuming S inherits all couplings from the Higgs boson evaluated at the LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) as functions of mass and c$\tau$ assuming S inherits all couplings from the Higgs boson evaluated at the LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) as functions of mass and c$\tau$ assuming S inherits all couplings from the Higgs boson evaluated at the LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) as functions of mass and c$\tau$ assuming S inherits all couplings from the Higgs boson evaluated at the LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) as functions of mass and c$\tau$ assuming S inherits all couplings from the Higgs boson evaluated at the LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for vector portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 2 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for vector portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for vector portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for vector portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for vector portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for gluon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 2 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 2 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 2 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for photon portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 4 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 4 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 4 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for higgs portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 1)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (2.5, 2.5)$ and 20 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 5 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 10 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 15 GeV LLP mass.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow \Psi\Psi$) for darkphoton portal, assuming $(\xi_{\omega}$, $\xi_{\Lambda}) = (1, 1)$ and 20 GeV LLP mass.

Measured min-$d_{xy}$ distribution in the 2-match category, after requiring the $d_{xy}$ muon to pass the isolation requirement $I_{\mathrm{rel}}^{\mathrm{PF}} <0.25$ (i. e., the B and D bins of the ABCD plane). Overlaid with a red histogram is the background predicted from the region of the ABCD plane failing the same requirement (the A and C bins), as well as three signal benchmark hypotheses (as defined in the legends), assuming $\alpha_D = \alpha_{\mathrm{EM}}$ (the fine-structure constant). The red hatched bands correspond to the background prediction uncertainty. The last bin includes the overflow.

Two-dimensional exclusion surfaces for $\Delta = 0.1 \, m_1$ as a function of the DM mass $m_1$ and the signal strength $y$, with $m_{A'} = 3 \, m_1$. Filled histograms denote observed limits on $\sigma(\mathrm{pp} \to A' \to \chi_2 \chi_1) \, \mathcal{B}(\chi_2 \to \chi_1 \mu^+ \mu^-)$. Solid (dashed) curves denote the observed (expected) exclusion limits at 95% CL, with 68% CL uncertainty bands around the expectation. Regions above the curves are excluded, depending on the $\alpha_D$ hypothesis: $\alpha_{\mathrm{D}} = \alpha_{\mathrm{EM}}$ (dark blue) or 0.1 (light magenta).

Two-dimensional exclusion surfaces for $\Delta = 0.4 \, m_1$ as a function of the DM mass $m_1$ and the signal strength $y$, with $m_{A'} = 3 \, m_1$. Filled histograms denote observed limits on $\sigma(\mathrm{pp} \to A' \to \chi_2 \chi_1) \, \mathcal{B}(\chi_2 \to \chi_1 \mu^+ \mu^-)$. Solid (dashed) curves denote the observed (expected) exclusion limits at 95% CL, with 68% CL uncertainty bands around the expectation. Regions above the curves are excluded, depending on the $\alpha_D$ hypothesis: $\alpha_{\mathrm{D}} = \alpha_{\mathrm{EM}}$ (dark blue) or 0.1 (light magenta).

Measured ratio of $\mathcal{B}_{4\mu}/\mathcal{B}_{2\mu}$

Measured branching fraction $\mathcal{B}_{4\mu}$

Measured branching fraction $\mathcal{B}_{4\mu}$

Observed and expected 95% CL upper limits on the product of the production cross section ($\sigma$) and the branching fraction, obtained after combining all categories with 138 $\mathrm{fb}^{−1}$ of data at $\sqrt{s}$ = 13 TeV for $\mathrm{V'}$ to $VV$ + $VH$ signal in HVT model C.

Total signal efficiency including all event categories as a function of the resonance mass $M_X$ for all DY/ggF signal models under consideration. The estimates from MC simulation (markers) are interpolated with a functional form (lines).

Total signal efficiency including all event categories as a function of the resonance mass $M_X$ for all VBF signal models under consideration. The estimates from MC simulation (markers) are interpolated with a functional form (lines).

Observed (points with error bars) and expected background (solid histograms) distributions for $E_{T}^{miss}$ in the signal region (a) SRL, (b) SRM and (c) SRH after the background-only fit applied to the CRs. The predicted signal distributions for the two models with a gluino mass of 2000 GeV and neutralino mass of 250 GeV (SRL), 1050 GeV (SRM) or 1950 GeV (SRH) are also shown for comparison. The uncertainties in the SM background are only statistical.

Observed and expected exclusion limit in the gluino-neutralino mass plane at 95% CL combined using the signal region with the best expected sensitivity at each point, for the full Run-2 dataset corresponding to an integrated luminosity of $139~\mathrm{fb}^{-1}$, for $\gamma/Z$ (a) and $\gamma/h$ (b) signal models. The black solid line corresponds to the expected limits at 95% CL, with the light (yellow) bands indicating the 1$\sigma$ exclusions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves, the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties. For each point in the higgsino-bino parameter space, the labels indicate the best-expected signal region, where L, M and H mean SRL, SRM and SRH, respectively.

Observed and expected exclusion limit in the gluino-neutralino mass plane at 95% CL combined using the signal region with the best expected sensitivity at each point, for the full Run-2 dataset corresponding to an integrated luminosity of $139~\mathrm{fb}^{-1}$, for $\gamma/Z$ (a) and $\gamma/h$ (b) signal models. The black solid line corresponds to the expected limits at 95% CL, with the light (yellow) bands indicating the 1$\sigma$ exclusions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves, the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties. For each point in the higgsino-bino parameter space, the labels indicate the best-expected signal region, where L, M and H mean SRL, SRM and SRH, respectively.

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Cutflow for the SRL selection, for two relevant signal points for both $\gamma/Z$ and $\gamma/h$ models, where the gluinos have mass of 2000 GeV and the neutralinos have a mass of 250 GeV (10000 generated events). The numbers are normalized to a luminosity of 139 $fb^{-1}$.

Cutflow for the SRM selection, for two relevant signal points for both $\gamma/Z$ and $\gamma/h$ models, where the gluinos have mass of 2000 GeV and the neutralinos have a mass of 1050 GeV (10000 generated events). The numbers are normalized to a luminosity of 139 $fb^{-1}$.

Cutflow for the SRH selection, for two relevant signal points for both $\gamma/Z$ and $\gamma/h$ models, where the gluinos have mass of 2000 GeV and the neutralinos have a mass of 1950 GeV (10000 generated events). The numbers are normalized to a luminosity of 139 $fb^{-1}$.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Trigger efficiency for simulated signal events as a function of the LLP decay position in the x&ndash;y plane (L_xy) for LLPs decaying in the barrel (|&eta;|&lt;1.4) for one of the low-E_T signal samples for HLT CalRatio triggers seeded by the high-E_T L1 triggers with E_T thresholds of 60 GeV and 100 GeV and by the two versions of the low-E_T L1 triggers. Only statistical uncertainties are shown.

Trigger efficiency for simulated signal events as a function of the LLP decay position in the x&ndash;y plane (L_xy) for LLPs decaying in the barrel (|&eta;|&lt;1.4) for one of the high-E_T signal samples for HLT CalRatio triggers seeded by the high-E_T L1 triggers with E_T thresholds of 60 GeV and 100 GeV and by the two versions of the low-E_T L1 triggers. Only statistical uncertainties are shown.

The NN output scores in the dijet control region for the low-E_T training with no adversary network. Statistical uncertainties are shown in all plots. In cases where training with adversary networks is considered, the systematic uncertainty related to modelling discrepancies is included as well.

The NN output scores in the dijet control region for the low-E_T training with an adversary network included. Statistical uncertainties are shown in all plots. In cases where training with adversary networks is considered, the systematic uncertainty related to modelling discrepancies is included as well.

The NN output scores in the dijet control region for the high-E_T training with no adversary network. Statistical uncertainties are shown in all plots. In cases where training with adversary networks is considered, the systematic uncertainty related to modelling discrepancies is included as well.

The NN output scores in the dijet control region for the high-E_T training with an adversary network included. Statistical uncertainties are shown in all plots. In cases where training with adversary networks is considered, the systematic uncertainty related to modelling discrepancies is included as well.

Distribution of the low-E_T per-event BDT in main data, BIB data and some of the benchmark signal samples after preselection. Only statistical uncertainties are shown.

Distribution of the high-E_T per-event BDT outputs in main data, BIB data and some of the benchmark signal samples after preselection. Only statistical uncertainties are shown.

Sequential impact of each requirement on the number of events passing the selection for the high-E_T selections. The signal columns represent the cumulative fraction of events passing the selection than the number of events.

Sequential impact of each requirement on the number of events passing the selection for the low-E_T selections. The signal columns represent the cumulative fraction of events passing the selection than the number of events.

Application of the modified ABCD method to the final high-E_T selections. The a priori estimate refers to the "pre-unblinding" case, where the data in region A are ignored by removing the Poisson constraint in that region and the signal strength is fixed to zero. This matches the simple N^bkg_A=(N^bkg_B&middot; N^bkg_C)/N^bkg_D relation. The a posteriori estimate refers to the "post-unblinding" case, including the observed data in region A in the background-only global fit, obtained by fixing the signal strength to 0 (background-only fit) or allowing it to float (signal-plus-background fit). The table also shows one set of representative signal yields in each selection for the signal-plus-background fit. Only statistical uncertainties are included in the quoted error of the background, while the uncertainties in the signal include those from both statistical and experimental sources.

Application of the modified ABCD method to the final low-E_T selections. The a priori estimate refers to the "pre-unblinding" case, where the data in region A are ignored by removing the Poisson constraint in that region and the signal strength is fixed to zero. This matches the simple N^bkg_A=(N^bkg_B&middot; N^bkg_C)/N^bkg_D relation. The a posteriori estimate refers to the "post-unblinding" case, including the observed data in region A in the background-only global fit, obtained by fixing the signal strength to 0 (background-only fit) or allowing it to float (signal-plus-background fit). The table also shows one set of representative signal yields in each selection for the signal-plus-background fit. Only statistical uncertainties are included in the quoted error of the background, while the uncertainties in the signal include those from both statistical and experimental sources.

95% CL expected and observed limits on the BR of SM Higgs bosons to pairs of neutral LLPs (B_{H&rarr; ss}), showing the &plusmn; 1 &sigma; (green) and &plusmn; 2 &sigma; (yellow) expected limit bands, as well as a comparison with the results from previous ATLAS searches [36,78]. The cross-section for SM Higgs boson gluon--gluon fusion production is assumed to be 48.6 pb.

Trigger efficiency of simulated signal events as a function of the LLP decay position in the z direction for LLPs decaying in the calorimeter endcaps (1.4 &le; |&eta;| &lt; 2.5) for one of the low-E_T signal samples for HLT CalRatio triggers seeded by the 60 GeV-high-E_T trigger and by the two versions of the low-E_T triggers.

Trigger efficiency of simulated signal events as a function of the LLP decay position in the z direction for LLPs decaying in the calorimeter endcaps (1.4 &le; |&eta;| &lt; 2.5) for one of the high-E_T signal samples for HLT CalRatio triggers seeded by the 60 GeV-high-E_T trigger and by the two versions of the low-E_T triggers.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The simulated efficiencies as a function of c times the mean proper lifetime (c&tau;) of s for several different MC samples. A weight-based extrapolation procedure is used to determine the efficiency at a given mean proper lifetime.

The event BDT in the dijet control region for the low-E_T training training. The uncertainties are combined statistical and ML modelling systematic uncertainty.

The event BDT in the dijet control region for the high-E_T training. The uncertainties are combined statistical and ML modelling systematic uncertainty.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; mass of 60 GeV compared to the results from the 2016-data analysis and subsequent combinations, where available.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; mass of 60 GeV compared to the results from the 2016-data analysis and subsequent combinations, where available.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; mass 125 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; mass 125 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available. The 125 GeV mediator is assumed to be the SM Higgs boson.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; mass 125 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available. The 125 GeV mediator is assumed to be the SM Higgs boson.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; mass of 200 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; mass of 400 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; mass 600 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; mass of 600 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; masses of of 1000 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; masses of of 1000 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; masses of of 1000 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available.

The 95% CL observed limits, expected limits and &plusmn; 1 &sigma; and 2 &sigma; bands for the &Phi; masses of of 1000 GeV, compared to the results from the 2016-data analysis and subsequent combinations, where available.

Efficiency for an event to enter Region A of the high-E_T selection, as a function of the two LLP kinematics, decay type and decay position. LLP p_T is binned in the ranges of [0, 50 , 100, 200 , 400, 1600] GeV (5 bins), the LLP decay position is binned in decay position in L_xy in [0, 1.5, 2, 2.5, 3, 3.5, 3.9, &infin;] m for LLPs with |&eta;| &lt; 1.5 and L_z in [0, 3.6, 4.2, 4.8, 5.5, 6, &infin; ] m for LLPs with |&eta;|geq1.5, 13 bins in all. Finally, four decay types are considered: LLPs decaying to pairs of c, b, t, or &tau; in bins 0, 1, 2, 3 respectively. The efficiency is presented as a function of "Bin Index", which is calculated as follows: Bin Index = (decay position bin index ) &times; (number of p_T bins &times; number of decay type bins) + p_T bin index * (number of decay type bin) + decay type bin index. The efficiency for a given pair-produced LLP sample can be obtained by summing the efficiency values for each event as extracted from this map, and dividing it by the total number of events in the sample. The efficiency map is symmetric between the LLPs, so the choice of LLP1 and LLP2 is arbitrary. For the high-E_T selections, for overall efficiencies above 0.5%, the results are typically accurate to around 25%, but below this the efficiency can be overestimated and therefore should not be used for re-interpration.

Efficiency for an event to enter Region A of the low-E_T selection, as a function of the two LLP kinematics, decay type and decay position. LLP p_T is binned in the ranges of [0, 50 , 100, 200 , 400, 1600] GeV (5 bins), the LLP decay position is binned in L_xy in [0, 1.5, 2, 2.5, 3, 3.5, 3.9, &infin;] m for LLPs with |&eta;| &lt; 1.5 and L_z in [0, 3.6, 4.2, 4.8, 5.5, 6, &infin; ] m for LLPs with |&eta;|geq1.5, 13 bins in all. Finally, four decay types are considered: LLPs decaying to pairs of c, b, t, or &tau; in bins 0, 1, 2, 3 respectively. The efficiency is presented as a function of "Bin Index", which is calculated as follows: Bin Index = (decay position bin index ) &times; (number of p_T bins &times; number of decay type bins) + p_T bin index * (number of decay type bin) + decay type bin index. The efficiency for a given pair-produced LLP sample can be obtained by summing the efficiency values for each event as extracted from this map, and dividing it by the total number of events in the sample. The efficiency map is symmetric between the LLPs, so the choice of LLP1 and LLP2 is arbitrary. For the low-E_T selections, for overall efficiencies above 0.15%, the results are typically accurate to around 33%, and below this the efficiency is typically accurate up to a factor of 3.

Expected event yields in each $m_{jj}$ bin for the different background processes in the SR of the VTR category, in the 2017 and 2018 samples. The background yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR. The expected signal contributions for a Higgs boson, produced in the non-VBF and VBF modes, decaying to invisible particles with a branching fraction of $\mathcal{B}(\text{H} \to \text{inv}) = 1$, and the observed event yields are also reported.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for MTR, 2017, signal region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for MTR, 2017, dimuon region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for MTR, 2018, signal region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for MTR, 2018, dimuon region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for VTR, 2017, signal region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for VTR, 2017, dimuon region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for VTR, 2018, signal region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for VTR, 2018, single electron region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for VTR, 2018, dielectron region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for VTR, 2018, dimuon region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for MTR, 2017, photon region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for MTR, 2018, photon region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for MTR, 2017, single electron region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for MTR, 2017, single muon region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for MTR, 2017, dielectron region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for MTR, 2018, single electron region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for MTR, 2018, single muon region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for MTR, 2018, dielectron region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for VTR, 2017, single electron region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for VTR, 2017, single muon region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for VTR, 2017, dielectron region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Event yields for the postfit background model from the background-only and signal-plus-background fits, as well as the observed event yields for VTR, 2018, single muon region. All signal and control regions from the 2017 and 2018 data sets are included in the fit.

Trigger efficiency measured as a function of $p_{T}^{miss}$ in data for the 2017 and 2018 samples for the MTR category. The efficiency is measured in a phase space with at least two jets in which the leading (trailing) jet has $p_{T}>80~\text{GeV}$ ($>40~\text{GeV}$). The two jets are required to form a system with $m_{jj}>200~\text{GeV}$ and $\Delta\Phi_{jj}<1.5$.

Trigger efficiency measured as a function of $p_{T}^{miss}$ in data for the 2017 and 2018 samples for the MTR category. The efficiency is measured in a phase space with at least two jets in which the leading (trailing) jet has $p_{T}>140~\text{GeV}$ ($>70~\text{GeV}$). The two jets are required to form a system with $m_{jj}>900~\text{GeV}$ and $\Delta\Phi_{jj}<1.8$.

Theory cross sections vs m($\widetilde{\chi}^0_1$) for SMS model TChiHH-G, $\widetilde{\chi}^0_1, \widetilde{\chi}^0_2$ NLSPs only.

The 95% CL observed and expected upper limits on the production cross sections of the TChiHH signal model as a function of the $\widetilde{\chi}^0_2$ and $\widetilde{\chi}^0_1$ masses. Relative to Fig. 13, the binning of the table is adjusted to align with the model points that were sampled. In addition, a small correction has been applied that affects the cross section upper limits for points with $m(\widetilde{\chi}^0_1) = 1$ GeV. For that row of points in the figure, the small region of model phase space reached by both the resolved and boosted selections was inadvertently double counted in computing the signal contribution to the predicted yield.

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Exclusion limits assuming the approximate-NNLO+NNLL cross sections

Cross section, 95% CL upper limit plus theory vs m($\widetilde{g}$) for SMS model T5HH.

Pre-fit background covariance matrix $\sigma_{xy}$ for the 22 analysis bins, ordered as in Fig. 10.

Pre-fit background correlation matrix $\rho_{xy}$ for the 22 analysis bins, ordered as in Fig. 10.

Significance vs $(m(\widetilde{\chi}^0_2), m(\widetilde{\chi}^0_1))$ for SMS model TChiHH.

Significance vs $m(\widetilde{g})$ for SMS model T5HH.

Likelihood profile at $(m(\widetilde{\chi}^0_2) = $275, $m(\widetilde{\chi}^0_1) = $1) for SMS model TChiHH.

Likelihood profile at $(m(\widetilde{\chi}^0_2) = $300, $m(\widetilde{\chi}^0_1) = $1) for SMS model TChiHH.

Likelihood profile at $(m(\widetilde{\chi}^0_2) = $300, $m(\widetilde{\chi}^0_1) = $50) for SMS model TChiHH.

Likelihood profile at $(m(\widetilde{\chi}^0_2) = $450, $m(\widetilde{\chi}^0_1) = $1) for SMS model TChiHH.

Likelihood profile at $(m(\widetilde{\chi}^0_2) = $750, $m(\widetilde{\chi}^0_1) = $1) for SMS model TChiHH.

Likelihood profile at $(m(\widetilde{\chi}^0_2) = $750, $m(\widetilde{\chi}^0_1) = $200) for SMS model TChiHH.

Efficiency vs $m(\widetilde{\chi}^0_1)$ for SMS model TChiHH-G. The denominator includes all 22 signal regions, and assumes $\mathcal{B}(\mathrm{H}$-->$\mathrm{b}\overline{\mathrm{b}})=100\%$.

Efficiency vs $m(\widetilde{g})$ for SMS model T5HH. The denominator includes all 22 signal regions, and no constraint on $\mathcal{B}(\mathrm{H}$-->$\mathrm{b}\overline{\mathrm{b}})$.

Efficiency vs $(m(\widetilde{\chi}^0_2), m(\widetilde{\chi}^0_1))$ for SMS model TChiHH. The denominator includes all 22 signal regions, and assumes $\mathcal{B}(\mathrm{H}$-->$\mathrm{b}\overline{\mathrm{b}})=100\%$.

Signal and control region populations for the resolved signature.

Signal and control region populations for the boosted signature.

Cut flow table for the resolved topology showing expected event yields at 137 fb$^{-1}$ for selected signal models as successive selections are applied.

Cut flow table for the boosted topology showing expected event yields at 137 fb$^{-1}$ for selected signal model points as successive selections are applied. Here the hadronic baseline is defined by $N_{\mathrm{jet}}\geq$2, \ptmiss$>$300 GeV, $\Delta\phi$ cuts, and the isolated lepton and track vetoes.

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Comparison of predicted and observed event yields in each of the control and signal regions in the dilepton channel after the fit to the data. The uncertainty band includes all uncertainties and their correlations.

Comparison of predicted and observed event yields in each of the control and signal regions in the single-lepton channels after the fit to the data. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $120\le p_T^H<200$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $200\le p_T^H<300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $120\le p_T^H<200$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $200\le p_T^H<300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 450$ GeV (yield only). The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton boosted SRs after the inclusive fit to the data for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton boosted SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for ${\Delta R^{{avg}}_{bb}}$ after the inclusive fit to the data in the single-lepton $CR^{5j}_{{\geq}4b\ lo}$ control region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for ${\Delta R^{{avg}}_{bb}}$ after the inclusive fit to the data in the single-lepton $CR^{5j}_{{\geq}4b\ hi}$ control region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit yields of signal ($S$) and total background ($B$) as a function of $\log (S/B)$, compared with data. Final-discriminant bins in all dilepton and single-lepton analysis regions are combined into bins of $\log (S/B)$, with the signal normalised to the SM prediction used for the computation of $\log (S/B)$. The signal is then shown normalised to the best-fit value and the SM prediction. The lower frame reports the ratio of data to background, and this is compared with the expected ${t\bar {t}H}$-signal-plus-background yield divided by the background-only yield for the best-fit signal strength (solid red line) and the SM prediction (dashed orange line).

Comparison between data and prediction for the reconstruction BDT score for the Higgs boson candidate identified using Higgs boson information, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the average $\Delta \eta $ between $b$-tagged jets, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the likelihood discriminant, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the average $\Delta R$ for all possible combinations of $b$-tagged jet pairs, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the DNN $P(H)$ output for the Higgs boson candidate after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the DNN $P(H)$ output for the Higgs boson candidate after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the reconstructed Higgs boson candidate mass for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit distribution of the reconstructed Higgs boson candidate mass for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit distribution of the reconstructed Higgs boson candidate mass for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Fitted values of the ${t\bar {t}H}$ signal strength parameter in the individual channels and in the inclusive signal-strength measurement.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the fit. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Pre-fit distribution of the number of jets in the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Pre-fit distribution of the number of jets in the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Pre-fit distribution of the number of jets in the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Post-fit distribution of the number of jets in the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the number of jets in the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the number of jets in the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Signal-strength measurements in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive signal strength.

95% CL simplified template cross-section upper limits in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive limit. The observed limits are shown (solid black lines), together with the expected limits both in the background-only hypothesis (dotted black lines) and in the SM hypothesis (dotted red lines). In the case of the expected limits in the background-only hypothesis, one- and two-standard-deviation uncertainty bands are also shown. The hatched uncertainty bands correspond to the theory uncertainty in the fiducial cross-section prediction in each bin.

The ratios $S/B$ (black solid line, referring to the vertical axis on the left) and $S/\sqrt{B}$ (red dashed line, referring to the vertical axis on the right) for each category in the inclusive analysis in the dilepton channel (left) and in the single-lepton channels (right), where $S$ ($B$) is the number of selected signal (background) events predicted by the simulation and normalised to a luminosity of 139 fb$^{-1}$ .

Comparison between data and prediction for the $\Delta R$ between the Higgs candidate and the ${t\bar {t}}$ candidate system, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the number of $b$-tagged jet pairs with an invariant mass within 30 GeV of 125 GeV, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the reconstruction BDT score for the Higgs boson candidate identified using Higgs boson information, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the $\Delta R$ between the two highest ${p_{{T}}}$ $b$-tagged jets, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the sum of $b$-tagging discriminants of jets from Higgs, hadronic top and leptonic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the sum of $b$-tagging discriminants of jets from Higgs, hadronic top and leptonic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the hadronic top candidate invariant mass, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the hadronic top candidate invariant mass, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the fraction of the sum of $b$-tagging discriminants of all jets not associated to the Higgs or hadronic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the fraction of the sum of $b$-tagging discriminants of all jets not associated to the Higgs or hadronic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $0\le {\hat{p}_{{T}}^{H}}<120$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $120\le {\hat{p}_{{T}}^{H}}<200$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $200\le {\hat{p}_{{T}}^{H}}<300$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $300\le {\hat{p}_{{T}}^{H}}<450$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for ${\hat{p}_{{T}}^{H}}\ge 450$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

95% confidence level upper limits on signal-strength measurements in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive signal-strength limit, after the fit used to extract multiple signal-strength parameters. The observed limits are shown (solid black lines), together with the expected limits both in the background-only hypothesis (dotted black lines) and in the SM hypothesis (dotted red lines). In the case of the expected limits in the background-only hypothesis, one- and two-standard-deviation uncertainty bands are also shown.

Post-fit correlation matrix (in percentages) between the $\mu $ values obtained in the STXS bins.

Performance of the Higgs boson reconstruction algorithms. For each row of `truth' ${\hat{p}_{{T}}^{H}}$, the matrix shows (in percentages) the fraction of Higgs boson candidates which are truth-matched to ${b\bar {b}}$ decays, with reconstructed $p_T^H$ in the various bins of the dilepton (left), single lepton resolved (middle) and boosted (right) channels.

Pre-fit event yields in the dilepton signal regions and control regions. All uncertainties are included except the $k({t\bar {t}+{\geq }1b})$ uncertainty that is not defined pre-fit. For the ${t\bar {t}H}$ signal, the pre-fit yield values correspond to the theoretical prediction and corresponding uncertainties. `Other sources' refers to s-channel, t-channel, $tW$, $tWZ$, $tZq$, $Z+$ jets and diboson events.

Post-fit event yields in the dilepton signal regions and control regions, after the inclusive fit in all channels. All uncertainties are included, taking into account correlations. For the ${t\bar {t}H}$ signal, the post-fit yield and uncertainties correspond to those in the inclusive signal-strength measurement. `Other sources' refers to s-channel, t-channel, $tW$, $tWZ$, $tZq$, $Z+$ jets and diboson events.

Pre-fit event yields in the single-lepton resolved and boosted signal regions and control regions. All uncertainties are included except the $k({t\bar {t}+{\geq }1b})$ uncertainty that is not defined pre-fit. For the ${t\bar {t}H}$ signal, the pre-fit yield values correspond to the theoretical prediction and corresponding uncertainties. `Other top sources' refers to s-channel, t-channel, $tWZ$ and $tZq$ events.

Post-fit event yields in the single-lepton resolved and boosted signal regions and control regions, after the inclusive fit in all channels. All uncertainties are included, taking into account correlations. For the ${t\bar {t}H}$ signal, the post-fit yield and uncertainties correspond to those in the inclusive signal-strength measurement. `Other top sources' refers to s-channel, t-channel, $tWZ$ and $tZq$ events.

Breakdown of the contributions to the uncertainties in $\mu$. The contributions from the different sources of uncertainty are evaluated after the fit. The $\Delta \mu $ values are obtained by repeating the fit after having fixed a certain set of nuisance parameters corresponding to a group of systematic uncertainties, and then evaluating $(\Delta \mu)^2$ by subtracting the resulting squared uncertainty of $\mu $ from its squared uncertainty found in the full fit. The same procedure is followed when quoting the effect of the ${t\bar {t}+{\geq }1b}$ normalisation. The total uncertainty is different from the sum in quadrature of the different components due to correlations between nuisance parameters existing in the fit.

Fraction (in percentages) of signal events, after SR and CR selections, originating from $b\bar {b}$, $WW$ and other remaining Higgs boson decay modes in the dilepton channel.

Fraction (in percentages) of signal events, after SR and CR selections, originating from $b\bar {b}$, $WW$ and other remaining Higgs boson decay modes in the single-lepton channels.

Predicted SM ${t\bar {t}H}$ cross-section in each of the five STXS ${\hat{p}_{{T}}^{H}}$ bins and signal acceptance times efficiency (including all event selection criteria) in each STXS bin as well as for the inclusive ${\hat{p}_{{T}}^{H}}$ range.

Number of expected signal events before the fit, after each selection requirement applied to enter the dilepton channel $SR^{\geq 4j}_{\geq 4b}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.

Number of expected signal events before the fit, after each selection requirement applied to enter the single-lepton channel resolved $SR^{\geq 6j}_{\geq 4b}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.

Number of expected signal events before the fit, after each selection requirement applied to enter the single-lepton channel boosted $SR_{boosted}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.

Fitted event-level b-tagging selection efficiency as a function of the dijet invariant mass. The uncertainties are evaluated by propagating the uncertainties on the fit parameters to the fitted efficiency.