Results of applying the BH algorithm to the mass spectra in the leading photon category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading electron category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading muon category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading small-R jet category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading bjet category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading large-R jet category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading photon category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading electron category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the leading muon category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the inclusive category, using the fitted background estimations from the initial step

Results of applying the BH algorithm to the mass spectra in the inclusive category, using the fitted background estimations from the initial step

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading small-R jet category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading bjet category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading large-R jet category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading photon category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading electron category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading muon category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading small-R jet category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading bjet category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading large-R jet category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading photon category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading electron category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the leading muon category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the inclusive category

Upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signal with a relative width value of 3% as a function of mass in the inclusive category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading small-R jet category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading bjet category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading large-R jet category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading photon category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading electron category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading muon category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading small-R jet category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading bjet category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading large-R jet category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading photon category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading electron category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the leading muon category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the inclusive category

Comparison of observed upper limits at 95% CL on the cross section times branching fraction times acceptance for a Gaussian-shaped signals with relative width values as a function of mass in the inclusive category

Acceptance times efficiency in an HVT model as a function of mass in the leading large-R jet category

Upper limits at 95% CL on the cross section times the branching fraction($W^{\prime} \to ZW$) for the HVT signal as functions of mass($m_{ZX}$) in the leading large-R jet category.The full(dashed) red curves correspond to the theoretical predictions of HVT model A(B)

Comparison of the identification efficiencies using standard and merged-ee reconstruction as a function of true PT(Z)

Background rejection factor as a function of signal efficiency. The red curve shows the BDT performance whereas the blue curve corresponds to that of a cut-based analysis relying only on the jet energyfraction deposited in the EM calorimeter

BH p-values of the 100 pseudo-experiments as a function of in the leading small-R jet category for an injected Gaussian-shaped signal with a relative width value of 3%. The fractions of the pseudo-experiments that have the correctly identified interval and BH p-values below the threshold of 0.01 are indicated. The background is derived by the background-only fit in the full fitting range.

BH p-values of the 100 pseudo-experiments as a function of in the leading small-R jet category for an injected Gaussian-shaped signal with a relative width value of 3%. The fractions of the pseudo-experiments that have the correctly identified interval and BH p-values below the threshold of 0.01 are indicated. The background is derived by the background-only fit in the range excluding the BH interval.

Fractions of pseudo-experiments in which the detected BH interval agrees with the injected mass point and the BH p-value is below 0.01 as a function of mass in the leading small-R jet category for Gaussian-shaped signal with a relative width of 3%.

Fractions of pseudo-experiments in which the detected BH interval agrees with the injected mass point and the BH p-value is below 0.01 as a function of mass in the leading small-R jet category for Gaussian-shaped signal with a relative width of 3%.

Distribution of exclusion upper limits on signal event yields at 95% CL from 1000 pseudo-experiments for Gaussian-shaped signals with relative width values of 3% at the low boundary of the limit-sensitive mass range for the ZX spectrum of the leading small-R jet category. The vertical line corresponds to the expected nominal exclusion limit at 95% CL.

Distribution of exclusion upper limits on signal event yields at 95% CL from 1000 pseudo-experiments for Gaussian-shaped signals with relative width values of 3% at the high boundary of the limit-sensitive mass range for the ZX spectrum of the leading small-R jet category. The vertical line corresponds to the expected nominal exclusion limit at 95% CL.

Distribution of exclusion upper limits on signal event yields at 95% CL from 1000 pseudo-experiments for Gaussian-shaped signals with relative width values of 5% at the low boundary of the limit-sensitive mass range for the ZX spectrum of the leading small-R jet category. The vertical line corresponds to the expected nominal exclusion limit at 95% CL.

Distribution of exclusion upper limits on signal event yields at 95% CL from 1000 pseudo-experiments for Gaussian-shaped signals with relative width values of 5% at the high boundary of the limit-sensitive mass range for the ZX spectrum of the leading small-R jet category. The vertical line corresponds to the expected nominal exclusion limit at 95% CL.

Distribution of exclusion upper limits on signal event yields at 95% CL from 1000 pseudo-experiments for Gaussian-shaped signals with relative width values of 10% at the low boundary of the limit-sensitive mass range for the ZX spectrum of the leading small-R jet category. The vertical line corresponds to the expected nominal exclusion limit at 95% CL.

Distribution of exclusion upper limits on signal event yields at 95% CL from 1000 pseudo-experiments for Gaussian-shaped signals with relative width values of 10% at the high boundary of the limit-sensitive mass range for the ZX spectrum of the leading small-R jet category. The vertical line corresponds to the expected nominal exclusion limit at 95% CL.

Data yields of the six event categories in the $Z\to e^+e^-$ and $\mu^+\mu^-$ decay channels. The merged-$e^+e^-$ events are included in the $e^+e^-$ channel, increasing the event yield, mainly in the leading large-$R$-jet category, by 0.6 %.

A list of mass spectra, event categories and their corresponding fit ranges, functional forms, numbers of free parameters and global $\chi^2$ $p$-values from background-only fits. The fit range values are rounded to the nearest 5 GeV. Here $f_1(x)=p_0\left(\mathrm{e}^{-p_1x}+p_2\mathrm{e}^{-(p_1+p_3)x}+p_4\mathrm{e}^{-(p_1+p_3+p_5)x}+\cdots\right)$ and $f_2(x)=p_0(1-x)^{p_1}x^{p_2+p_3\ln\!x+p_4\ln^2\!x+\cdots}$

A list of mass spectra, event categories and their corresponding signal-sensitive mass ranges. The initial BH $p$-value is obtained by using the background derived from the background-only fit in the full fit range, whereas the new BH $p$-value uses the background derived from the background-only fit in the range excluding the initial BH interval.

A list of mass spectra, event categories, relative width values of Gaussian-shaped signals and limit-sensitive mass ranges and fractions, corresponding to the mass values in the previous column, of pseudo-experiments having exclusion upper limits higher than the nominal exclusion limit at 95% CL

Cutflow of HVT model $A$ signals ($W^\prime \to ZW \to \ell\ell qq$) with $m_{W^\prime} = 1$ TeV and $m_{W^\prime} = 4$ TeV based on MC simulations.

Acceptance times efficiency ($\mathcal{A} \times \epsilon$) values in % in the dominant event category for $p^Z_{\mathrm{T}} > 100$ GeV in the $Z \to \ell^{+}\ell^{-}$ decay 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.

Particle-level TEEC results for the third HT2 bin

Particle-level TEEC results for the fourth HT2 bin

Particle-level TEEC results for the fifth HT2 bin

Particle-level TEEC results for the sixth HT2 bin

Particle-level TEEC results for the seventh HT2 bin

Particle-level TEEC results for the eighth HT2 bin

Particle-level TEEC results for the ninth HT2 bin

Particle-level TEEC results for the tenth HT2 bin

Particle-level ATEEC results

Particle-level ATEEC results for the first HT2 bin

Particle-level ATEEC results for the second HT2 bin

Particle-level ATEEC results for the third HT2 bin

Particle-level ATEEC results for the fourth HT2 bin

Particle-level ATEEC results for the fifth HT2 bin

Particle-level ATEEC results for the sixth HT2 bin

Particle-level ATEEC results for the seventh HT2 bin

Particle-level ATEEC results for the eighth HT2 bin

Particle-level ATEEC results for the ninth HT2 bin

Particle-level ATEEC results for the tenth HT2 bin

Particle-level TEEC predictions

Particle-level TEEC predictions for the first HT2 bin

Particle-level TEEC predictions for the second HT2 bin

Particle-level TEEC predictions for the third HT2 bin

Particle-level TEEC predictions for the fourth HT2 bin

Particle-level TEEC predictions for the fifth HT2 bin

Particle-level TEEC predictions for the sixth HT2 bin

Particle-level TEEC predictions for the seventh HT2 bin

Particle-level TEEC predictions for the eighth HT2 bin

Particle-level TEEC predictions for the ninth HT2 bin

Particle-level TEEC predictions for the tenth HT2 bin

Particle-level ATEEC predictions

Particle-level ATEEC predictions for the first HT2 bin

Particle-level ATEEC predictions for the second HT2 bin

Particle-level ATEEC predictions for the third HT2 bin

Particle-level ATEEC predictions for the fourth HT2 bin

Particle-level ATEEC predictions for the fifth HT2 bin

Particle-level ATEEC predictions for the sixth HT2 bin

Particle-level ATEEC predictions for the seventh HT2 bin

Particle-level ATEEC predictions for the eighth HT2 bin

Particle-level ATEEC predictions for the ninth HT2 bin

Particle-level ATEEC predictions for the tenth HT2 bin

Fitted values for the strong coupling constant extracted from TEEC with MMHT 2014 PDF

Fitted values for the strong coupling constant extracted from TEEC with NNPDF 3.0

Fitted values for the strong coupling constant extracted from TEEC with CT14 PDF

Fitted values for the strong coupling constant extracted from ATEEC with MMHT 2014 PDF

Fitted values for the strong coupling constant extracted from ATEEC with NNPDF 3.0

Fitted values for the strong coupling constant extracted from ATEEC with CT14 PDF

The unfolded differential charge asymmetry as a function of the transverse momentum of the top pair system. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NNLO in QCD and NLO in EW theory are listed. The scale uncertainty is obtained by varying renormalisation and factorisation scales independently by a factor of 2 or 0.5 around $\mu_0$ to calculate the maximum and minimum value of the asymmetry, respectively. The nominal value $\mu_0$ is chosen as $H_T/4$. The variations in which one scale is multiplied by 2 while the other scale is divided by 2 are excluded. Finally, the scale and MC integration uncertainties are added in quadrature.

The unfolded differential charge asymmetry as a function of the longitudinal boost of the top pair system. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NNLO in QCD and NLO in EW theory are listed. The scale uncertainty is obtained by varying renormalisation and factorisation scales independently by a factor of 2 or 0.5 around $\mu_0$ to calculate the maximum and minimum value of the asymmetry, respectively. The nominal value $\mu_0$ is chosen as $H_T/4$. The variations in which one scale is multiplied by 2 while the other scale is divided by 2 are excluded. Finally, the scale and MC integration uncertainties are added in quadrature.

The unfolded inclusive leptonic asymmetry. The unfolded $A_C^{\ell\bar{\ell}}$ is obtained in the reduced phase-space defined by the requirement $|\Delta |\eta_{\ell\bar{\ell}}||<2.5$. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NLO in QCD and NLO in EW theory are listed. The theory uncertainty is obtained by varying both scales by a factor of 0.5 or 2.0 to calculate the minimum and maximum value of the asymmetry, respectively.

The unfolded differential leptonic asymmetry as a function of the invariant mass of the di-lepton pair. The unfolded $A_C^{\ell\bar{\ell}}$ is obtained in the reduced phase-space defined by the requirement $|\Delta |\eta_{\ell\bar{\ell}}||<2.5$. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NLO in QCD and NLO in EW theory are listed. The theory uncertainty is obtained by varying both scales by a factor of 0.5 or 2.0 to calculate the minimum and maximum value of the asymmetry, respectively.

The unfolded differential leptonic asymmetry as a function of the transverse momentum of the di-lepton pair. The unfolded $A_C^{\ell\bar{\ell}}$ is obtained in the reduced phase-space defined by the requirement $|\Delta |\eta_{\ell\bar{\ell}}||<2.5$. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NLO in QCD and NLO in EW theory are listed. The theory uncertainty is obtained by varying both scales by a factor of 0.5 or 2.0 to calculate the minimum and maximum value of the asymmetry, respectively.

The unfolded differential leptonic asymmetry as a function of the longitudinal boost of the di-lepton pair. The unfolded $A_C^{\ell\bar{\ell}}$ is obtained in the reduced phase-space defined by the requirement $|\Delta |\eta_{\ell\bar{\ell}}||<2.5$. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NLO in QCD and NLO in EW theory are listed. The theory uncertainty is obtained by varying both scales by a factor of 0.5 or 2.0 to calculate the minimum and maximum value of the asymmetry, respectively.

Individual 68% and 95% CL bounds on the relevant Wilson coefficients of the SM Effective Field Theory in units of $\text{TeV}^{-2}$. The bounds are derived from the $A_C^{t\bar{t}}$ inclusive measurement. The experimental uncertainties are accounted for, in the form of the complete covariance matrix that keeps track of correlations between bins for the differential measurement. The theory uncertainty from the NNLO QCD + NLO EW calculation is included by explicitly varying the renormalization and factorization scales, or the parton density functions, in the calculation and registering the variations in the intervals.

Individual 68% and 95% CL bounds on the relevant Wilson coefficients of the SM Effective Field Theory in units of $\text{TeV}^{-2}$. The bounds are derived from the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement. The experimental uncertainties are accounted for, in the form of the complete covariance matrix that keeps track of correlations between bins for the differential measurement. The theory uncertainty from the NNLO QCD + NLO EW calculation is included by explicitly varying the renormalization and factorization scales, or the parton density functions, in the calculation and registering the variations in the intervals.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ inclusive measurement. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement for $\beta_{z,t\bar{t}}$ $\in$ [0,0.3]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement for $\beta_{z,t\bar{t}}$ $\in$ [0.3,0.6]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement for $\beta_{z,t\bar{t}}$ $\in$ [0.6,0.8]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement for $\beta_{z,t\bar{t}}$ $\in$ [0.8,1]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ < 500 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ $\in$ [500,750] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ $\in$ [750,1000] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ $\in$ [1000,1500] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ > 1500 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement for $p_{T,t\bar{t}}$ < 30 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement for $p_{T,t\bar{t}}$ $\in$ [30,120] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement for $p_{T,t\bar{t}}$ > 120 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ inclusive measurement. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement for $\beta_{z,\ell\bar{\ell}}$ $\in$[0,0.3]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement for $\beta_{z,\ell\bar{\ell}}$ $\in$[0.3,0.6]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement for $\beta_{z,\ell\bar{\ell}}$ $\in$[0.6,0.8]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement for $\beta_{z,\ell\bar{\ell}}$ $\in$[0.8,1]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement for $m_{\ell\bar{\ell}}$ < 200 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement for $m_{\ell\bar{\ell}}$ $\in$ [200,300] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement for $m_{\ell\bar{\ell}}$ $\in$ [300,400] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement for $m_{\ell\bar{\ell}}$ > 400 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement for $p_{T,\ell\bar{\ell}}$ < 20 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement for $p_{T,\ell\bar{\ell}}$ $\in$ [20, 70] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement for $p_{T,\ell\bar{\ell}}$ > 70 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{t\bar{t}}$ inclusive measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{\ell\ell}$ inclusive measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Covariance matrix for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement.

Covariance matrix for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement.

Covariance matrix for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement.

Covariance matrix for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement.

Covariance matrix for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement.

Covariance matrix for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement.

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.

Measured unfolded differential cross-section as a function of the dilepton transverse momentum $p^{ll}_{\mathrm{T}}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the the four-body transverse momentum $p^{ll\gamma\gamma}_{\mathrm{T}}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the diphoton invariant mass $m_{\gamma\gamma}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the four-body invariant mass $m_{ll\gamma\gamma}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Expected and observed $95\%$ confidence intervals for the coupling parameters $f_{T,j}/\Lambda^{4}$ of transverse dimension-8 operators. All parameter values outside of the stated range are excluded at the chosen confidence level. No unitarity constraints are applied.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,8}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,0}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,1}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,2}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,5}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,6}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,7}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,9}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

The $m_{jj}$ distributions for the CRs and the BDT classifier response distribution for the SR after the fit in all regions. The dashed line shows the total background distribution before the fit. The vertical error bars on the data points correspond to the data's statistical uncertainty. Overflows are included in the last bin. The uncertainty band corresponds to the combination of the MC statistical uncertainty and systematic uncertainties obtained in the fit.

The $m_{jj}$ distributions for the CRs and the BDT classifier response distribution for the SR after the fit in all regions. The dashed line shows the total background distribution before the fit. The vertical error bars on the data points correspond to the data's statistical uncertainty. Overflows are included in the last bin. The uncertainty band corresponds to the combination of the MC statistical uncertainty and systematic uncertainties obtained in the fit.

The $m_{jj}$ distributions for the CRs and the BDT classifier response distribution for the SR after the fit in all regions. The dashed line shows the total background distribution before the fit. The vertical error bars on the data points correspond to the data's statistical uncertainty. Overflows are included in the last bin. The uncertainty band corresponds to the combination of the MC statistical uncertainty and systematic uncertainties obtained in the fit.

Observed and expected event yields for the signal and all of the background processes considered in this analysis after the fit to the data in all of the regions. The uncertainty on the expected yield is the combination of statistical and systematic uncertainties obtained in the fit. The individual uncertainties can be correlated and do not necessarily add in quadrature to equal the total background uncertainty.

Summary of the event yield for processes in all regions, after the fit over all regions. The dashed line shows the total background distribution before the fit. The vertical error bars on the data points correspond to the data's statistical uncertainty. The uncertainty band corresponds to the combination of the MC statistical uncertainty and systematic uncertainties obtained in the fit.

Impact of different components of systematic uncertainty on the measured cross section, without taking into account the correlations. The impact is calculated by fixing the value of the corresponding nuisance parameters to the values obtained in the fit used to measure the cross section, performing the fit, estimating the signal strength uncertainty and subtracting it from the nominal uncertainty in quadrature.

Fitted POI values for this analysis, the previous ATLAS analysis, and their combination. The first and second columns present the values obtained in the individual analyses. The third column presents the values obtained in the combination.

Observed and expected one-dimensional limits on dimension 8 aQGC parameters. Limits are obtained by setting all aQGCs parameters except one to zero. Unitarity is not preserved.

Observed and expected one-dimensional limits on dimension 8 aQGC parameters in the region, where unitarity is preserved. Cutoff scales in MC, $E_{c}$, are given for each parameter. Limits are obtained by setting all aQGCs parameters except one to zero.

The $E_{T}^{\gamma}$ distribution in the SR after the fit in the control regions. The red (green) line shows the expected number of events in the case of non-zero EFT coefficient $f_{T0}/\Lambda^4$ ($f_{M0}/\Lambda^4$) with the value shown in the legend. The vertical error bars on the data points correspond to the data statistical uncertainty. Overflows are included in the last bin. The uncertainty band corresponds to the combination of the MC statistical uncertainty and systematic uncertainties obtained in the fit.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{T0}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 600$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{T5}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 600$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{T8}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 600$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{T9}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 600$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{M0}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 400$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{M1}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 400$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{M2}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 400$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Correlation matrix of the parameters for the fit in all regions. Only parameters with absolute value of the correlation coefficient larger than 10 are shown.

Distribution of ${E}_{T}^{miss}$ after the fit of the multijet backgrounds, in the electron channel, in the $t\bar{t}$ 3-jets VR, without applying the cut on ${E}_{T}^{miss}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of ${E}_{T}^{miss}$ after the fit of the multijet backgrounds, in the electron channel, in the $t\bar{t}$ 4-jets VR, without applying the cut on ${E}_{T}^{miss}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the signal region, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $W$+jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $t\bar{t}$ 3-jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $t\bar{t}$ 4-jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Expected distributions of the MEM discriminant $P(S|X)$ in the SR, for the s-channel single-top signal, and for the $t\bar{t}$ and $W$+jets backgrounds, for MEM discriminant values larger than $2.0\times10^{-4}$. Each distribution is normalised to unity. The binning is the same as the optimised binning used in the signal extraction fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $W$+jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty band indicates the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $t\bar{t}$ 3-jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty band indicates the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $t\bar{t}$ 4-jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty bands indicate the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR before the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$. The lower panel of the figure shows the ratio of the data to the prediction, with different vertical axis ranges. The uncertainty band indicates the total uncertainties and their correlations in each bin. The uncertainties in the $t\bar{t}$ and $W$+jets normalisation factors, as well as in the s-channel signal cross-section, are not defined pre-fit and therefore not included. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR after the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$. The lower panel of the figure shows the ratio of the data to the prediction, with different vertical axis ranges. The uncertainty band indicates the total uncertainties and their correlations in each bin. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR after the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$, after subtraction of all backgrounds. The fitted distribution for the simulation of the signal is shown together with the post-fit uncertainty in the backgrounds. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Pre-fit and post-fit event yields in the SR, for MEM discriminant values larger than $2.0\times10^{-4}$. The central value of the event yield for each process is calculated by summing the values of the discriminant bin contents, using the nominal expected yield for the pre-fit value, and the best-fit estimate for the post-fit value. The error includes statistical and systematic uncertainties summed in quadrature. All sources of systematic uncertainties are included, taking into account correlations and anti-correlations in the post-fit case. The uncertainties in the $t\bar{t}$ and $W$+jets normalisation factors, as well as in the s-channel signal cross-section, are not defined pre-fit and therefore only included in the post-fit uncertainties.

Observed impact of the different sources of uncertainty on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. The 'Systematic uncertainties' category combines all sources of systematic uncertainties. The statistical uncertainty is obtained by repeating the fit after having fixed all nuisance parameters, including the $t\bar{t}$ and $W$+jets normalisation factors. 'Total' gives the total uncertainty on the measurement.

Observed impact of the different sources of $t\bar{t}$ modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, and 'ME/PS matching' to the matching of the ME to the parton shower.

Observed impact of the different sources of s-channel modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, as described in Section 7 in the paper.

Observed impact of the different sources of t-channel modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, as described in Section 7 in the paper.

Observed impact of the different sources of $tW$ modelling uncertainty on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, and '$t\bar{t}$ overlap' to the algorithm removing the overlap between $tW$ and $t\bar{t}$ production at NLO, as described in Section 7 in the paper.

Observed impact of the different sources of PDF uncertainties on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root.

Comparison between data and prediction after the fit to data in the signal region for the leading-jet $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the leading-jet $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the subleading-jet $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the subleading-jet $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the lepton $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the lepton $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the ${E}_{T}^{miss}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the $m_{T}^{W}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Nuisance parameters ranked according to their post-fit impacts on the best-fit value of the ratio $\mu$ of the measured cross-section to the predicted cross-section. In the figure, only the 20 nuisance parameters with the largest post-fit impacts are shown. The empty (solid) blue rectangles illustrate the pre-fit (post-fit) impact on $\mu$, corresponding to the upper axis. The pre-fit (post-fit) impact of each nuisance parameter, $\Delta\mu$, is calculated as the difference in the fitted value of $\mu$ between the nominal fit and the fit when fixing the corresponding nuisance parameter to $\hat{\theta}\pm\Delta\theta$ ($\hat{\theta}\pm\Delta\hat{\theta}$), where $\hat{\theta}$ is the best-fit value of the nuisance parameter and $\Delta\theta$ ($\Delta\hat{\theta}$) is its pre-fit (post-fit) uncertainty. Several systematic uncertainties are split into different nuisance parameters, which are indicated by NP. JES (JER) indicates jet energy scale (resolution), and $\gamma$ indicates a nuisance parameter associated to the MC statistics in one of the 18 bins numbered from 0 to 17. The black points show the best-fit values of the nuisance parameters, with the error bars representing the post-fit uncertainties. Each nuisance parameter is shown wrt. its nominal value, $\theta_0$, and in units of its pre-fit uncertainty, except the free-floating normalisation factors of the $t\bar{t}$ and $W$+jets backgrounds, and the parameters associated to the MC statistics in each bin, for which the post-fit values and uncertainties are shown.

Correlation matrix of the nuisance parameters and of the ratio $\mu$ of the measured cross-section to the predicted cross-section. The correlations are given after the fit to data. In the figure, only the parameters which have a correlation of at least 0.2 with any other parameter are shown.

Distribution of the MEM discriminant $P(S|X)$ in the SR for MEM discriminant values larger than $2.0\times10^{-4}$, for the collision data used for the measurement, and for 1000 pseudo-data replicas, generated using a bootstrapping technique, in order to assess the statistical correlations between this measurement and others, for the purpose of combinations. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the <a href="https://zenodo.org/record/5361038">BootstrapGenerator</a> software package , which implements a technique described in <a href="https://cds.cern.ch/record/2759945/">ATL-PHYS-PUB-2021-011</a>. The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. Each pseudo-data replica is assigned an index, ranging from 0 to 999, corresponding to the random number index used consistently for each observed data event.

Measured values of the signal cross-section and of the $t\bar{t}$ and $W$+jets normalisation factors, obtained by statistical-only fits to the collision data used for the measurement, and to 1000 pseudo-data replicas, generated using a bootstrapping technique, in order to assess the statistical correlations between this measurement and others, for the purpose of combinations. The central values and their statistical uncertainties are obtained by repeating the fit after having fixed all nuisance parameters, except the $t\bar{t}$ and $W$+jets normalisation factors, which are let free-floating (unlike for the statistical uncertainty on the cross-section quoted in the paper). The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the <a href="https://zenodo.org/record/5361038">BootstrapGenerator</a> software package , which implements a technique described in <a href="https://cds.cern.ch/record/2759945/">ATL-PHYS-PUB-2021-011</a>. The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. Each pseudo-data replica is assigned an index, ranging from 0 to 999, corresponding to the random number index used consistently for each observed data event.

Expected 2-D likelihood scan of the signal strength $\mu_{\tau\tau}$ versus the $CP$-mixing angle $\phi_{\tau}$.

Free-floating parameters in the measurement.