Differential cross section as a function of $m_{\gamma\gamma}$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

Differential cross section as a function of $p_{T,\gamma\gamma}$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

Differential cross section as a function of $a_{T,\gamma\gamma}$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

Differential cross section as a function of $\phi_{\eta}*$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

Differential cross section as a function of $\pi-\Delta\phi_{\gamma\gamma}$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

Differential cross section as a function of $|cos\theta*|^{(CS)}$. The table contains the values measured in data and theory predictions from SHERPA, DIPHOX and NNLOJET.

The contribution from different systematic uncertainties to the measured $t\bar{t}t\bar{t}$ production cross section, grouped in categories.

The results of the fitted signal strength $\mu$ in the 1L/2LOS and 2LSS/3L combined channel.

The results of fitted inclusive ${t\bar{t}t\bar{t}}$ cross-section in the 1L/2LOS and 2LSS/3L combined channel.

Comparison between data and prediction for the distribution of the sum of the pseudo-continuous b-tagging score over the six jets with the highest score in the 1L,$\geq$9j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of the sum of the pseudo-continuous b-tagging score over the six jets with the highest score in the 1L,$\geq$9j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of the sum of the pseudo-continuous b-tagging score over the six jets with the highest score in the 2LOS,$\geq$7j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of the sum of the pseudo-continuous b-tagging score over the six jets with the highest score in the 2LOS,$\geq$7j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,$\geq$8j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,$\geq$8j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,$\geq$6j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,$\geq$6j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of number of jets in the 1L,$\geq$8j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of number of jets in the 1L,$\geq$8j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of number of jets in the 2LOS,$\geq$6j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of number of jets in the 2LOS,$\geq$6j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of b-jets multiplicity in the 1L,$\geq$8j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of b-jets multiplicity in the 1L,$\geq$8j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of b-jets multiplicity in the 2LOS,$\geq$6j,$\geq$3b region before the fit.

Comparison between data and prediction for the distribution of b-jets multiplicity in the 2LOS,$\geq$6j,$\geq$3b region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,4b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,4b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,$\geq$5b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,$\geq$5b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bL signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bL signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bH signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bH signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,4b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,4b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,$\geq$5b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,$\geq$5b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,7j,$\geq$4b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,7j,$\geq$4b signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bL signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bL signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bH signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bH signal region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,$\geq$4b signal region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,$\geq$4b signal region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,9j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 1L,$\geq$10j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,7j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,7j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bV validation region before the fit.

Comparison between data and prediction for the distribution of the BDT score in the 2LOS,$\geq$8j,3bV validation region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bL control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bL control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bH control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,3bH control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,9j,3bL control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,9j,3bL control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,9j,3bH control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,9j,3bH control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,4b control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,4b control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,$\geq$5b control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 1L,8j,$\geq$5b control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bL control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bL control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bH control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,3bH control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,7j,3bL control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,7j,3bL control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,7j,3bH control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,7j,3bH control region after the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,$\geq$4b control region before the fit.

Comparison between data and prediction for the distribution of the scalar sum of all jet and lepton pT in the event in the 2LOS,6j,$\geq$4b control region after the fit.

Summary of the observed and expected 95% CL upper limits on the cross section for $\mathrm{LQ}_{3}^{\mathrm{d}}$ pair production as a function of $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$ under the assumptions of B=1.

Summary of the observed and expected 95% CL upper limits on B as a function of $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$.

Summary of the observed and expected 95% CL upper limits on B as a function of $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$.

Cutflow of the preselection requirements (see Section 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the preselection requirements (see Section 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in 1$\ell+\geq 1\tau$ channel (see Table 3) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in 1$\ell+\geq 1\tau$ channel (see Table 3) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the 2$\ell$OS+$\geq 1\tau$ channel (see Table 4) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the 2$\ell$OS+$\geq 1\tau$ channel (see Table 4) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the $2\ell$SS/$3\ell+\geq 1\tau$ channel (see Table 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. In this channel, two signal regions (SR-L and SR-H) are defined based on $p_{\mathrm{T}, 1}^{\tau}$, with SR-L and SR-H requiring $125< p_{\mathrm{T}, 1}^{\tau} < 225$ GeV and $p_{\mathrm{T}, 1}^{\tau}>225$ GeV, respectively. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Cutflow of the signal region requirements in the $2\ell$SS/$3\ell+\geq 1\tau$ channel (see Table 5) for $\mathrm{LQ}_{3}^{\mathrm{d}}$ signals with $m_{\mathrm{LQ}_{3}^{\mathrm{d}}}$=0.9, 1.1, and 1.3 TeV, assuming B=1. Events that satisfy the preselection requirements are considered. In this channel, two signal regions (SR-L and SR-H) are defined based on $p_{\mathrm{T}, 1}^{\tau}$, with SR-L and SR-H requiring $125< p_{\mathrm{T}, 1}^{\tau} < 225$ GeV and $p_{\mathrm{T}, 1}^{\tau}>225$ GeV, respectively. The yields correspond to an integrated luminosity of 139 fb$^{-1}$.

Distributions of $s_{\text{T}}$ in the single-tau one-bin SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Distributions of $m_{\text{T}}(\tau)$ in the single-tau one-bin SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Distributions of $\Sigma m_{\text{T}}(b_{1,2})$ in the single-tau $p_{\text{T}}(\tau)$-binned SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Distributions of $p_{\text{T}}(\tau)$ in the single-tau $p_{\text{T}}(\tau)$-binned SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Observed event yields in data ('Observed') and expected event yields for SM background processes obtained from the background-only fit ('Total bkg.' and rows below) in the signal regions of the di-tau and single-tau channels. The quoted uncertainties include both the statistical and systematic uncertainties and are truncated at zero yield. By construction, no $t\bar{t}$ (2 real $\tau$) events can pass the selections in the single-tau channel. As the individual uncertainties are correlated, they do not add in quadrature to equal the total background uncertainty.

From left to right: upper limits at the 95% confidence level (CL) on the visible cross section ($\sigma_\text{vis}$) and on the number of signal events ($S_{\text{obs}}^{95}$). The third column ($S_{\text{exp}}^{95}$) shows the upper limit at the 95% CL on the number of signal events, given the expected number (and $\pm 1\,\sigma$ excursions on the expectation) of background events. The last two columns indicate the confidence level observed for the background-only hypothesis ($\text{CL}_{b}$), the discovery $p$-value ($p(s=0)$) and the significance ($Z$). In the di-tau SR, where fewer events are observed than predicted by the fitted background estimate, the $p$-value is capped at 0.5.

Expected and observed exclusion contours at the 95% confidence level for the vector third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{v}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ into a quark and a charged lepton. The plot shows the exclusion contour for the minimal-coupling scenario. The limits are derived from the binned single-tau signal region.

Expected and observed exclusion contours at the 95% confidence level for the vector third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{v}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ into a quark and a charged lepton. The plot shows the exclusion contour for the minimal-coupling scenario. The limits are derived from the binned single-tau signal region.

Expected and observed exclusion contours at the 95% confidence level for the vector third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{v}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ into a quark and a charged lepton. The plot shows the exclusion contour for vector leptoquarks with additional gauge couplings. The limits are derived from the binned single-tau signal region.

Expected and observed exclusion contours at the 95% confidence level for the vector third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{v}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ into a quark and a charged lepton. The plot shows the exclusion contour for vector leptoquarks with additional gauge couplings. The limits are derived from the binned single-tau signal region.

Exclusion contours at the 95% confidence level for the stop-stau signal model as a function of the masses of the top squark $m(\tilde{t}_{1})$ and of the tau slepton $m(\tilde{\tau}_{1})$. Expected and observed limits are shown for the present search in comparison to observed limits from previous ATLAS analyses based on data from Run-1 of the LHC at $\sqrt{s} = 8$ TeV [Eur. Phys. J. C 76 (2016)] and on a partial dataset from Run 2 at $\sqrt{s} = 13$ TeV [Phys. Rev. D 98 (2018) 032008]. The green band indicates the limit on the mass of the tau slepton (for a massless LSP) from the LEP experiments.

Exclusion contours at the 95% confidence level for the stop-stau signal model as a function of the masses of the top squark $m(\tilde{t}_{1})$ and of the tau slepton $m(\tilde{\tau}_{1})$. Expected and observed limits are shown for the present search in comparison to observed limits from previous ATLAS analyses based on data from Run-1 of the LHC at $\sqrt{s} = 8$ TeV [Eur. Phys. J. C 76 (2016)] and on a partial dataset from Run 2 at $\sqrt{s} = 13$ TeV [Phys. Rev. D 98 (2018) 032008]. The green band indicates the limit on the mass of the tau slepton (for a massless LSP) from the LEP experiments.

Expected and observed exclusion contours at the 95% confidence level for the scalar third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{u}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow q\ell)$ into a quark and a charged lepton. The plot shows the exclusion contour for up-type leptoquarks $\text{LQ}_{3}^{\text{u}})$ with charge $+2/3e$. The limits are derived from the binned single-tau signal region. Shown in gray for comparison are the observed exclusion-limit contours from the previous ATLAS publication that targets the same leptoquark models but is based on a subset of the Run-2 data [JHEP 06 (2019) 144]. In this previous publication five different analyses are considered that target not only the final state studied here but also the final states that correspond to a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow q\ell)$ of 0 or 1, leading to the concave shapes of the gray exclusion contours.

Expected and observed exclusion contours at the 95% confidence level for the scalar third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{u}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow q\ell)$ into a quark and a charged lepton. The plot shows the exclusion contour for up-type leptoquarks $\text{LQ}_{3}^{\text{u}})$ with charge $+2/3e$. The limits are derived from the binned single-tau signal region. Shown in gray for comparison are the observed exclusion-limit contours from the previous ATLAS publication that targets the same leptoquark models but is based on a subset of the Run-2 data [JHEP 06 (2019) 144]. In this previous publication five different analyses are considered that target not only the final state studied here but also the final states that correspond to a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow q\ell)$ of 0 or 1, leading to the concave shapes of the gray exclusion contours.

Expected and observed exclusion contours at the 95% confidence level for the scalar third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{d}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow q\ell)$ into a quark and a charged lepton. The plot shows the exclusion contour for down-type leptoquarks $\text{LQ}_{3}^{\text{d}})$ with charge $-1/3e$. The limits are derived from the binned single-tau signal region. Shown in gray for comparison are the observed exclusion-limit contours from the previous ATLAS publication that targets the same leptoquark models but is based on a subset of the Run-2 data [JHEP 06 (2019) 144]. In this previous publication five different analyses are considered that target not only the final state studied here but also the final states that correspond to a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow q\ell)$ of 0 or 1, leading to the concave shapes of the gray exclusion contours.

Expected and observed exclusion contours at the 95% confidence level for the scalar third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{d}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow q\ell)$ into a quark and a charged lepton. The plot shows the exclusion contour for down-type leptoquarks $\text{LQ}_{3}^{\text{d}})$ with charge $-1/3e$. The limits are derived from the binned single-tau signal region. Shown in gray for comparison are the observed exclusion-limit contours from the previous ATLAS publication that targets the same leptoquark models but is based on a subset of the Run-2 data [JHEP 06 (2019) 144]. In this previous publication five different analyses are considered that target not only the final state studied here but also the final states that correspond to a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow q\ell)$ of 0 or 1, leading to the concave shapes of the gray exclusion contours.

Upper limits on the signal cross section at the 95 % confidence level for the stop-stau signal model.

Upper limits on the signal cross section at the 95 % confidence level for the scalar third-generation leptoquark signal model with up-type leptoquarks.

Upper limits on the signal cross section at the 95 % confidence level for the scalar third-generation leptoquark signal model with down-type leptoquarks.

Upper limits on the signal cross section at the 95 % confidence level for the vector third-generation leptoquark signal model with minimal coupling (MC).

Upper limits on the signal cross section at the 95 % confidence level for the vector third-generation leptoquark signal model with additional gauge couplings (YM).

Acceptance of the one-bin signal region of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the one-bin signal region of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the signal region of the di-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the one-bin signal region of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the signal region of the di-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the one-bin signal region of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the signal region of the di-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the signal region of the di-tau channel for pair production of top squarks with decays via tau sleptons.

Cutflow for the benchmark signal model $m(\tilde{t}_{1}) = 1350$ GeV, $m(\tilde{\tau}_{1}) = 1090$ GeV for the di-tau SR. The simulated sample contains 30,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\tilde{t}_{1}) = 1350$ GeV, $m(\tilde{\tau}_{1}) = 1090$ GeV for the single-tau one-bin SR. The simulated sample contains 30,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\tilde{t}_{1}) = 1350$ GeV, $m(\tilde{\tau}_{1}) = 1090$ GeV for the single-tau multi-bin SR. The simulated sample contains 30,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{u}}) = 1.2$ TeV, $\beta = 0.5$ for the di-tau SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{u}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau one-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{u}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau multi-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{d}}) = 1.2$ TeV, $\beta = 0.5$ for the di-tau SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{d}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau one-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{d}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau multi-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the minimal-coupling scenario for the di-tau SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the minimal-coupling scenario for the single-tau one-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the minimal-coupling scenario for the single-tau multi-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the Yang--Mills scenario for the di-tau SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the Yang--Mills scenario for the single-tau one-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the Yang--Mills scenario for the single-tau multi-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Expected +- 1 sigma 95% CL exclusion limit for the Zprime-2HDM model.

Expected +- 2 sigma 95% CL exclusion limit for the Zprime-2HDM model.

Observed 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected +- 1 sigma 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected +- 2 sigma 95% CL exclusion limit for the 2HDM+a model ggF production.

Observed 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected +- 1 sigma 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected +- 2 sigma 95% CL exclusion limit for the 2HDM+a model bbA production.

Observed 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected +- 1 sigma 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected +- 2 sigma 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected and observed upper limits at 95% CL on cross-section for Zprime-2HDM model.

Expected and observed upper limits at 95% CL on cross-section for ggF producton in the 2HDM+a model.

Expected and observed upper limits at 95% CL on cross-section for bbA producton in the 2HDM+a model.

Model-independent upper limits on the visible cross-section $σ_{vis, $h(\bar{b})+DM} ≡ σ_{h+DM} \times B(h \to b\bar{b}) \times \mathcal{A} \times \epsilon$ in the different signal regions.

Theory cross-section for Zprime-2HDM model.

Theory cross-section for bbA production in the 2HDM+a model.

Theory cross-section for ggF production in the 2HDM+a model.

Distribution of Higgs boson candidate mass in 2b region with MET=150-200 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=200-350 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=350-500 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=500-750 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET > 750 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=150-200 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=200-350 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=350-500 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET > 500 GeV.

Yields in 1-lepton control region.

Yields in 2-lepton control region.

MET distribution in 2b region of the 0-lepton channel.

MET distribution in 3b region of the 0-lepton channel.

Expected signal yields after certain selection cuts in 2b region with MET=150-200 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=200-350 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=350-500 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=500-750 GeV.

Expected signal yields after certain selection cuts in 2b region with MET > 750 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=150-200 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=200-350 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=350-500 GeV.

Expected signal yields after certain selection cuts in 3b region with MET > 500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET > 750 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET>500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET > 750 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET > 500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET > 750 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET > 500 GeV.

List of relative uncertainties of the measured inclusive $t\bar{t}\text{Z}$ cross section from the combined fit. The uncertainties are symmetrised for presentation and grouped into the categories described in the text. The quadratic sum of the individual uncertainties is not equal to the total uncertainty due to correlations introduced by the fit.

The definitions of the trilepton signal regions: for the inclusive measurement, a combination of the regions with pseudo-continuous $b$-tagging 3$\ell$-Z-1$b$4$j$-PCBT and 3$\ell$-Z-2$b$3$j$-PCBT is used, whereas for the differential measurement, only the region 3$\ell$-Z-2$b$3$j$, with a fixed $b$-tagging WP is employed.

The definitions of the trilepton signal regions: for the inclusive measurement, a combination of the regions with pseudo-continuous $b$-tagging 3$\ell$-Z-1$b$4$j$-PCBT and 3$\ell$-Z-2$b$3$j$-PCBT is used, whereas for the differential measurement, only the region 3$\ell$-Z-2$b$3$j$, with a fixed $b$-tagging WP is employed.

The definitions of the four tetralepton signal regions. The regions are defined to target different $b$-jet multiplicities and flavour combinations of the non-Z leptons.

The definitions of the four tetralepton signal regions. The regions are defined to target different $b$-jet multiplicities and flavour combinations of the non-Z leptons.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|$ in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the absolute parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the normalised parton-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel.

The total correlation matrix of the absolute particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}$ of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the absolute value of rapidity of the $Z$ boson in the 3$\ell$+4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{l \textrm{non-}Z}$ in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta y (Z, t_{\textrm{lep}})|/\pi$ in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 3$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (l_{t}^{+}, l_{\bar{t}}^{-})|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $|\Delta \phi (t\bar{t}, Z)|/\pi$ in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the $p_{\textrm{T}}^{t\bar{t}}$ in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel.

The total correlation matrix of the normalised particle-level differential cross-section measured in the fiducial phase-space as a function of the number of jets in the 4$\ell$ channel.

Summary of the total uncertainty in the background prediction for each SR of the tt0L-low, tt0L-high, tt1L and tt2L analysis channels in the statistical combination. Their dominant contributions are indicated by individual lines. Individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

$E_{\text{T}}^{\text{miss}}$ distribution in SR0X for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

$E_{\text{T}}^{\text{miss}}$ distribution in SRWX for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

$E_{\text{T}}^{\text{miss}}$ distribution in SRTX for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Representative fit distribution in the different flavour leptons signal region for the tt2L analysis: each bin of such distribution, starting from the red arrow, corresponds to a single SR included in the fit. 'FNP' includes the contribution from fake/non-prompt lepton background arising from jets (mainly $\pi/K$, heavy-flavour hadron decays and photon conversion) misidentified as leptons, estimated in a purely data-driven way. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The total uncertainty in the SM expectation is represented with hatched bands and the expected distributions for selected signal models are shown as dashed lines.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$t\bar{t}$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$tW$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$tj$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$t\bar{t}$ model, defined as the number of selected reconstructed events divided by the acceptance.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$tW$ model, defined as the number of selected reconstructed events divided by the acceptance.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$tj$ model, defined as the number of selected reconstructed events divided by the acceptance.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Measured fiducial differential cross sections for the jet multiplicity in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $\Delta R_{Z,j}^{min}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $r_{Z,j}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $r_{Z,j}$ in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $r_{Z,j}$ in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the H$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $\Delta R_{Z,j}^{min}$ in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The statistical, systematic, and luminosity uncertainties are given.

Systematic uncertainties for the Z boson p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the leading jet p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $\Delta R_{Z,j}^{min}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $r_{Z,j}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $r_{Z,j}$ in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $r_{Z,j}$ in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the H$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $\Delta R_{Z,j}^{min}$ in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the Z boson p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the leading jet p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the jet multiplicity in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the jet multiplicity in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the $\Delta R_{Z,j}^{min}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the $r_{Z,j}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the $r_{Z,j}$ in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the $r_{Z,j}$ in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the jet multiplicity in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the jet multiplicity in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the H$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the $\Delta R_{Z,j}^{min}$ in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with leptons at the Born-level to the cross section calculated with dressed leptons as a function of the jet multiplicity in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the Z boson p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the leading jet p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the jet multiplicity in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the jet multiplicity in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the $\Delta R_{Z,j}^{min}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the $r_{Z,j}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the $r_{Z,j}$ in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the $r_{Z,j}$ in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the jet multiplicity in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the jet multiplicity in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the H$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the $\Delta R_{Z,j}^{min}$ in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Correction scale factor from the cross section calculated with an overlap removal with jets of pT greater than 100 GeV to the cross section calculated with an overlap removal with jets of pT greater than 30 GeV as a function of the jet multiplicity in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, averaging the electron and muon channels, derived with Sherpa2.2.11. The systematic uncertainty is obtained with an enveloppe around scale factors computed from Sherpa2.2.1 and MG5_aMC+Py8 CKKWL.

Measured fiducial differential cross sections for the Z boson p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the leading jet p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $\Delta R_{Z,j}^{min}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $r_{Z,j}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $r_{Z,j}$ in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $r_{Z,j}$ in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the H$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the $\Delta R_{Z,j}^{min}$ in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Measured fiducial differential cross sections for the jet multiplicity in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The statistical, systematic, and luminosity uncertainties are given.

Systematic uncertainties for the Z boson p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the leading jet p$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $\Delta R_{Z,j}^{min}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $r_{Z,j}$ in the high-p$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $r_{Z,j}$ in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $r_{Z,j}$ in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the collinear region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the back-to-back region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the H$_{\mathrm{T}}$ in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the $\Delta R_{Z,j}^{min}$ in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Systematic uncertainties for the jet multiplicity in the high-S$_{\mathrm{T}}$ region in Z($\to \ell^{+} \ell^{-}$) + high p$_{\mathrm{T}}$ jets events, where the EW Zjj contribution is treated as signal and not subtracted as background. The uncertainties are presented as a percentage of the measured cross-section for the upward variation of each source of uncertainty in each bin.

Average timing distributions for SR data and the estimated background as determined by the background-only fit, in each of the five exclusive $\rho$ categories. For comparison, the expected timing shapes for a few different signal models are superimposed, with each model labeled by the values of the $\tilde\chi^0_1$ mass and lifetime, as well as decay mode. To provide some indication of the variations in signal yield and shape, three signal models are shown for each of the $\tilde\chi^0_1$ decay modes, namely $\tilde\chi^0_1$ $\rightarrow$ $H \tilde G$ and $\tilde\chi^0_1$ $\rightarrow$ $Z \tilde G$. The models shown include a rather low $\tilde\chi^0_1$ mass value of 135 GeV for lifetimes of either 2 ns or 10 ns, and a higher $\tilde\chi^0_1$ mass value which is near the 95% CL exclusion limit for each decay mode for a lifetime of 2 ns. Each signal model is shown with the signal normalization corresponding to a BR value of unity for the decay mode in question.

Average timing distributions for SR data and the estimated background as determined by the background-only fit, in each of the five exclusive $\rho$ categories. For comparison, the expected timing shapes for a few different signal models are superimposed, with each model labeled by the values of the $\tilde\chi^0_1$ mass and lifetime, as well as decay mode. To provide some indication of the variations in signal yield and shape, three signal models are shown for each of the $\tilde\chi^0_1$ decay modes, namely $\tilde\chi^0_1$ $\rightarrow$ $H \tilde G$ and $\tilde\chi^0_1$ $\rightarrow$ $Z \tilde G$. The models shown include a rather low $\tilde\chi^0_1$ mass value of 135 GeV for lifetimes of either 2 ns or 10 ns, and a higher $\tilde\chi^0_1$ mass value which is near the 95% CL exclusion limit for each decay mode for a lifetime of 2 ns. Each signal model is shown with the signal normalization corresponding to a BR value of unity for the decay mode in question.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ branching ratio to the SM Higgs boson, where the assumed cross-section is for higgsino production, and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 - $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$). Several signal hypotheses are overlaid that are labelled by the $\tilde\chi^0_1$ mass, all with a fixed $\tilde\chi^0_1$ lifetime of 2 ns.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ branching ratio to the SM Higgs boson, where the assumed cross-section is for higgsino production, and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 - $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$). Several signal hypotheses are overlaid that are labelled by the $\tilde\chi^0_1$ mass, all with a fixed $\tilde\chi^0_1$ lifetime of 2 ns.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ branching ratio to the SM Higgs boson, where the assumed cross-section is for higgsino production, and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 - $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$). Several signal hypotheses are overlaid that are labelled by the $\tilde\chi^0_1$ mass, all with a fixed $\tilde\chi^0_1$ lifetime of 2 ns.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ branching ratio to the SM Higgs boson, where the assumed cross-section is for higgsino production, and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 - $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$). Several signal hypotheses are overlaid that are labelled by the $\tilde\chi^0_1$ mass, all with a fixed $\tilde\chi^0_1$ lifetime of 2 ns.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

Cutflow for an example higgsino signal with mass 225 GeV and lifetime 10 ns, in the H decay mode. Acceptance is defined at truth level, and efficiency compares the events passing at reconstruction level with respect to truth.

Cutflow for an example higgsino signal with mass 225 GeV and lifetime 10 ns, in the Z decay mode. Acceptance is defined at truth level, and efficiency compares the events passing at reconstruction level with respect to truth.

Cutflow for an example higgsino signal with mass 225 GeV and lifetime 2 ns, in the H decay mode. Acceptance is defined at truth level, and efficiency compares the events passing at reconstruction level with respect to truth.

Cutflow for an example higgsino signal with mass 225 GeV and lifetime 2 ns, in the Z decay mode. Acceptance is defined at truth level, and efficiency compares the events passing at reconstruction level with respect to truth.

Acceptance across the H decay mode signal grid, calculated using truth information. The selection applied corresponds to the model-independent signal region (i.e. the standard SR with $t_{\text{avg}$ > 0.9 ns).

Acceptance across the Z decay mode signal grid, calculated using truth information. The selection applied corresponds to the model-independent signal region (i.e. the standard SR with $t_{\text{avg}$ > 0.9 ns).

Efficiency across the H decay mode signal grid, calculated using reco information. The selection applied corresponds to the model-independent signal region (i.e. the standard SR with $t_{\text{avg}$ > 0.9 ns). Here, the numerator is the signal yield passing the reco selection and the denominator is the signal yield passing the truth selection.

Efficiency across the Z decay mode signal grid, calculated using reco information. The selection applied corresponds to the model-independent signal region (i.e. the standard SR with $t_{\text{avg}$ > 0.9 ns). Here, the numerator is the signal yield passing the reco selection and the denominator is the signal yield passing the truth selection.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ and photon isolation cone radius $R=0.4$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ and photon isolation cone radius $R=0.4$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ and photon isolation cone radius $R=0.4$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ and photon isolation cone radius $R=0.2$.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ and isolation cone radius $0.2$ at NNLO QCD.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ at NNLO QCD.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $|\eta^{\gamma}|<0.6$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $0.6<|\eta^{\gamma}|<0.8$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $0.8<|\eta^{\gamma}|<1.37$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $1.56<|\eta^{\gamma}|<1.81$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $1.81<|\eta^{\gamma}|<2.01$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $2.01<|\eta^{\gamma}|<2.37$.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $|\eta^{\gamma}|<0.6$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $0.6<|\eta^{\gamma}|<0.8$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $0.8<|\eta^{\gamma}|<1.37$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $1.56<|\eta^{\gamma}|<1.81$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $1.81<|\eta^{\gamma}|<2.01$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $2.01<|\eta^{\gamma}|<2.37$ at NNLO QCD.