Fraction of events where both bosons are longitudinally polarized in the region with $p_T^Z>200$ GeV using two unconstrained parameters.

Numbers of observed and expected events in the 00-enhanced signal regions, before the fit. The total uncertainties are quoted.

Summary of the relative uncertainties in the measured longitudinal-longitudinal fractions $f_{00}$. The uncertainties are reported as percentages.

Numbers of observed and expected events in the 00-enhanced signal regions, after the fit. The total uncertainties are quoted.

$|\Delta Y(l_{W}Z)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 20$ GeV of the TT state.

$|\Delta Y(l_{W}Z)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 40$ GeV of the TT state.

$|\Delta Y(l_{W}Z)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 70$ GeV of the TT state.

$|\Delta Y(l_{W}Z)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 20$ GeV of of the sum of all polarization.

$|\Delta Y(l_{W}Z)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 40$ GeV of of the sum of all polarization.

$|\Delta Y(l_{W}Z)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 70$ GeV of of the sum of all polarization.

$|\Delta Y(WZ)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 20$ GeV of the TT state.

$|\Delta Y(WZ)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 40$ GeV of the TT state.

$|\Delta Y(WZ)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 70$ GeV of the TT state.

$|\Delta Y(WZ)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 20$ GeV of of the sum of all polarization.

$|\Delta Y(WZ)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 40$ GeV of of the sum of all polarization.

$|\Delta Y(WZ)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 70$ GeV of of the sum of all polarization.

The $\mathcal{D}$ value for the unfolded $|\Delta Y(l_{W}Z)|$ distributions of the TT polarization state as a function of the $p_T^{WZ}$ cut value.

The $\mathcal{D}$ value for the unfolded $|\Delta Y(WZ)|$ distributions of the TT polarization state as a function of the $p_T^{WZ}$ cut value.

Measured normalized differential $|\Delta Y(l_{W}Z)|$ cross-section of the TT state in the RAZ Region $p_{T}^{WZ}< 20$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(l_{W}Z)|$ cross-section of the TT state in the RAZ Region $p_{T}^{WZ}< 40$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(l_{W}Z)|$ cross-section of the TT state in the RAZ Region $p_{T}^{WZ}< 70$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(WZ)|$ cross-section of the TT state in the RAZ Region $p_{T}^{WZ}< 20$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(WZ)|$ cross-section of the TT state in the RAZ Region $p_{T}^{WZ}< 40$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(WZ)|$ cross-section of the TT state in the RAZ Region $p_{T}^{WZ}< 70$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(l_{W}Z)|$ cross-section of the sum of polarization states in the RAZ Region $p_{T}^{WZ}< 20$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(l_{W}Z)|$ cross-section of the sum of polarization states in the RAZ Region $p_{T}^{WZ}< 40$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(l_{W}Z)|$ cross-section of the sum of polarization states in the RAZ Region $p_{T}^{WZ}< 70$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(WZ)|$ cross-section of the sum of polarization states in the RAZ Region $p_{T}^{WZ}< 20$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(WZ)|$ cross-section of the sum of polarization states in the RAZ Region $p_{T}^{WZ}< 40$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(WZ)|$ cross-section of the sum of polarization states in the RAZ Region $p_{T}^{WZ}< 70$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Summary of the relative uncertainties in the measured depth of the TT state in the RAZ Region $p_{T}^{WZ}<20$ GeV using $|\Delta Y(l_{W}Z)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the TT state in the RAZ Region $p_{T}^{WZ}<40$ GeV using $|\Delta Y(l_{W}Z)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the TT state in the RAZ Region $p_{T}^{WZ}<70$ GeV using $|\Delta Y(l_{W}Z)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the TT state in the RAZ Region $p_{T}^{WZ}<20$ GeV using $|\Delta Y(WZ)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the TT state in the RAZ Region $p_{T}^{WZ}<40$ GeV using $|\Delta Y(WZ)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the TT state in the RAZ Region $p_{T}^{WZ}<70$ GeV using $|\Delta Y(WZ)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the sum of polarization states in the RAZ Region $p_{T}^{WZ}<20$ GeV using $|\Delta Y(l_{W}Z)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the sum of polarization states in the RAZ Region $p_{T}^{WZ}<40$ GeV using $|\Delta Y(l_{W}Z)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the sum of polarization states in the RAZ Region $p_{T}^{WZ}<70$ GeV using $|\Delta Y(l_{W}Z)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the sum of polarization states in the RAZ Region $p_{T}^{WZ}<20$ GeV using $|\Delta Y(WZ)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the sum of polarization states in the RAZ Region $p_{T}^{WZ}<40$ GeV using $|\Delta Y(WZ)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the sum of polarization states in the RAZ Region $p_{T}^{WZ}<70$ GeV using $|\Delta Y(WZ)|$. The uncertainties are reported as percentages.

Corrected $m(HH)$ distribution in the resolved $4b$ signal region (dots), after the fit under the background-only hypothesis. The error bars on the $4b$ points represent the Poisson uncertainties corresponding to their event yields. The background model (teal histogram) is shown with its total post-fit uncertainty (gray band). The final bin includes overflow. Representative spin-0 signal hypotheses (dashed, dotted, and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. The bottom panel shows the difference between the $4b$ distribution and the background model, relative to the background model. No significant excess of data relative to the SM background is observed.

Corrected $m(HH)$ distribution in the resolved $4b$ signal region (dots), after the fit under the background-only hypothesis. The error bars on the $4b$ points represent the Poisson uncertainties corresponding to their event yields. The background model (teal histogram) is shown with its total post-fit uncertainty (gray band). The final bin includes overflow. Representative spin-2 signal hypotheses (dashed, dotted, and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. The bottom panel shows the difference between the $4b$ distribution and the background model, relative to the background model. No significant excess of data relative to the SM background is observed.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The steps up to the $b$-tag categorization are shown for the spin-0.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The steps up to the $b$-tag categorization are shown for the spin-2.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The efficiencies of the three b-tag categories are shown for the spin-0 scenario; this efficiency is obtained after the other selection steps including the SR definition. The signal efficiency in the 4b region has a maximum around 1.5 TeV. Above that value the track-jets start to merge together, and for the highest resonance masses the 2b category becomes the most efficient.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The efficiencies of the three b-tag categories are shown for the spin-2 scenario; this efficiency is obtained after the other selection steps including the SR definition. The signal efficiency in the 4b region has a maximum around 1.5 TeV. Above that value the track-jets start to merge together, and for the highest resonance masses the 2b category becomes the most efficient.

Comparison of the background model (stacked histograms) with data (dots) in the $2b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

Comparison of the background model (stacked histograms) with data (dots) in the $3b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

Comparison of the background model (stacked histograms) with data (dots) in the $4b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

The $m(HH)$ distributions in the boosted $2b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $2b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $3b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $3b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $4b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $4b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

Expected (dashed black lines) and observed (solid black lines) 95% CL upper limits on the cross-section of resonant $HH$ production in the spin-0 signal models. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits (colored bands) are shown. Expected limits using each of the resolved and boosted channels individually (dashed colored lines) are shown. The nominal $H\rightarrow b\bar{b}$ branching ratio is taken as 0.582.

Expected (dashed black lines) and observed (solid black lines) 95% CL upper limits on the cross-section of resonant $HH$ production in the spin-2 signal models. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits (colored bands) are shown. Expected limits using each of the resolved and boosted channels individually (dashed colored lines) are shown. The theoretical prediction for the bulk RS model with $k/\bar{M}_{\text{Pl}} = 1$ (solid red line) is shown; the decrease below 350 GeV is due to a sharp reduction in the $G^{*}_{\text{KK}} \rightarrow HH$ branching ratio. The nominal $H\rightarrow b\bar{b}$ branching ratio is taken as 0.582.

Measured differential cross section as a function of observable $ \Delta R (l,l)$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ N_{jets}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{Jet1}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{Jet2}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{Jet2}/p_{T}^{Jet1}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ m_{jj}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ m_{ll\gamma j}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ H_{T}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{\gamma} / \sqrt{H_{T}}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ \Delta \phi (Jet,\gamma)$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma j}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ \phi_{CS}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ \cos \theta_{CS}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma} / m_{ll\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma} / m_{ll\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma} / m_{ll\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll} - p_{T}^{\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll} - p_{T}^{\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll} - p_{T}^{\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma j}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma j}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma j}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{\gamma} / \sqrt{H_{T}}$ (Fig. 8 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ H_{T}$ (Fig. 8 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ \Delta \phi (Jet,\gamma)$ (Fig. 8 (c))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ \Delta R (l,l)$ (Fig. 5 (d))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll} - p_{T}^{\gamma}$ (Fig. 5 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll} - p_{T}^{\gamma} \textrm{ in bin } p_{T}^{ll} + p_{T}^{\gamma} < 200 GeV$ (Fig. 11 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll} - p_{T}^{\gamma} \textrm{ in bin } 200 GeV < p_{T}^{ll} + p_{T}^{\gamma} < 300 GeV$ (Fig. 11 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll} - p_{T}^{\gamma} \textrm{ in bin } p_{T}^{ll} + p_{T}^{\gamma} > 300 GeV$ (Fig. 11 (c))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ m_{jj}$ (Fig. 7 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ m_{ll\gamma j}$ (Fig. 7 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ N_{jets}$ (Fig. 6 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{Jet1}$ (Fig. 6 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{Jet2}$ (Fig. 6 (c))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{Jet2}/p_{T}^{Jet1}$ (Fig. 6 (d))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll}$ (Fig. 5 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma j}$ (Fig. 8 (d))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma j} \textrm{ in bin } p_{T}^{ll\gamma} < 50 GeV$ (Fig. 12 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma j} \textrm{ in bin } 50 GeV < p_{T}^{ll\gamma} < 75 GeV$ (Fig. 12 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma j} \textrm{ in bin } p_{T}^{ll\gamma} > 75 GeV$ (Fig. 12 (c))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma} / m_{ll\gamma} \textrm{ in bin } 125 GeV < m_{ll\gamma} < 200 GeV$ (Fig. 10 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma} / m_{ll\gamma} \textrm{ in bin } 200 GeV < m_{ll\gamma} < 300 GeV$ (Fig. 10 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma} / m_{ll\gamma} \textrm{ in bin } m_{ll\gamma} > 300 GeV$ (Fig. 10 (c))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll} + p_{T}^{\gamma}$ (Fig. 5 (c))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ \cos \theta_{CS}$ (Fig. 9 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ \phi_{CS}$ (Fig. 9 (a))

Normalisation factors for the major background processes and ttW ($\Delta \eta < 0$) from the fit using only statistical uncertainties (fixing Nuisance Parameters (NP) to their fitted values) to data in all analysis regions.

Normalisation factors for the major background processes and ttW ($\Delta \eta < 0$) from the fit using using both statistical and systematic uncertainties to data in all analysis regions.

Correlation matrix between the Normalisation Factors in the fit using only statistical uncertainties (fixing Nuisance Parameters (NP) to their fitted values) to data in all analysis regions.

The most relevant systematic uncertainties ranked by their impact on the leptonic charge asymmetry ($A_c^{\ell}$) parameter at reconstructed level. The impact of the uncertainties is shown before and after the combined profile-likelihood fit to data in the signal and control regions. Pulls introduced by the fitting procedure are also shown. The entries shown in bold are the uncertainties of the freely floating background normalisations. ME stands for "matrix element", PS for "parton shower" and JER for "jet energy resolution". The gamma-uncertainties refer to the MC statistical uncertainties in a specific region and bin.

The most relevant systematic uncertainties ranked by their impact on the ttW Normalisation Factor ($\Delta \eta < 0$) parameter at reconstructed level. The impact of the uncertainties is shown before and after the combined profile-likelihood fit to data in the signal and control regions. Pulls introduced by the fitting procedure are also shown. ME stands for "matrix element", PS for "parton shower" and JER for "jet energy resolution". The gamma-uncertainties refer to the MC statistical uncertainties in a specific region and bin.

The predicted and observed numbers of events in the control and signal regions. The predictions are shown before the fit to data. The indicated uncertainties consider statistical as well as all experimental and theoretical systematic uncertainties. Background categories with event yields with a value of 1.0e-6 have a negligible contribution to the region and are manually set to that value.

The predicted and observed numbers of events in the control and signal regions. The predictions are shown after the fit to data. The indicated uncertainties consider statistical as well as all experimental and theoretical systematic uncertainties. Background categories with event yields with a value of 1.0e-6 have a negligible contribution to the region and are manually set to that value.

A summary of the uncertainties in the background estimates for SR-Gtt-0L-C. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-1L-B. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-1L-M1. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-1L-M2. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-1L-C. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-B. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-M. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-C. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtb-B. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtb-M. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtb-C. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-2100-1. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-1800-1. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-2300-1200. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gtt-1900-1400. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-2800-1400. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-2300-1000. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-2100-1600. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

A summary of the uncertainties in the background estimates for SR-Gbb-2000-1800. The individual experimental and theoretical uncertainties are assumed to be uncorrelated and are combined by adding in quadrature.

Results of the background-only fit extrapolated to SR_Gtt_0L_B in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_0L_M1 in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_0L_M2 in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_0L_C in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_1L_B in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_1L_M1 in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_1L_M2 in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_1L_C in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_B in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_M in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_C in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtb_B in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtb_M in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtb_C in the CC analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_2100_1 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_1800_1 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_2300_1200 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gtt_1900_1400 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_2800_1400 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_2300_1000 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_2100_1600 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Results of the background-only fit extrapolated to SR_Gbb_2000_1800 in the NN analysis, for both the total expected background yields and the main contributing background processes. The quoted uncertainties include both experimental and theoretical systematics. The data in the SRs are not included in the fit. The background category $t\bar{t}+X$ includes $t\bar{t} W/Z$, $t\bar{t} H$ and $t\bar{t} t\bar{t}$ events. The row ``Pre-fit background'' provides the total background prediction when the $t\bar{t}$ and $Z+$jets normalisations are obtained from theoretical calculation, taking into account the kinematic weights described in Section 5.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the NN analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed (left) 95\%~CL exclusion limits on the gluino mass as a function of BR$(\tilde{g} \to b\bar{b}\tilde\chi^{0}_{1}$) (vertical) and BR$(\tilde{g} \to t\bar{t}\tilde\chi^{0}_{1}$) (horizontal) for Gtb models with $m(\tilde\chi^{0}_{1}) = 1$~GeV, obtained from the CC analysis.

Expected (right) 95\%~CL exclusion limits on the gluino mass as a function of BR$(\tilde{g} \to b\bar{b}\tilde\chi^{0}_{1}$) (vertical) and BR$(\tilde{g} \to t\bar{t}\tilde\chi^{0}_{1}$) (horizontal) for Gtb models with $m(\tilde\chi^{0}_{1}) = 1$~GeV, obtained from the CC analysis.

Observed (left) 95\%~CL exclusion limits on the gluino mass as a function of BR$(\tilde{g} \to b\bar{b}\tilde\chi^{0}_{1}$) (vertical) and BR$(\tilde{g} \to t\bar{t}\tilde\chi^{0}_{1}$) (horizontal) for Gtb models with $m(\tilde\chi^{0}_{1}) = 600$~GeV, obtained from the CC analysis.

Expected (right) 95\%~CL exclusion limits on the gluino mass as a function of BR$(\tilde{g} \to b\bar{b}\tilde\chi^{0}_{1}$) (vertical) and BR$(\tilde{g} \to t\bar{t}\tilde\chi^{0}_{1}$) (horizontal) for Gtb models with $m(\tilde\chi^{0}_{1}) = 600$~GeV, obtained from the CC analysis.

Observed (left) 95\%~CL exclusion limits on the gluino mass as a function of BR$(\tilde{g} \to b\bar{b}\tilde\chi^{0}_{1}$) (vertical) and BR$(\tilde{g} \to t\bar{t}\tilde\chi^{0}_{1}$) (horizontal) for Gtb models with $m(\tilde\chi^{0}_{1}) = 1$~TeV, obtained from the CC analysis.

Expected (right) 95\%~CL exclusion limits on the gluino mass as a function of BR$(\tilde{g} \to b\bar{b}\tilde\chi^{0}_{1}$) (vertical) and BR$(\tilde{g} \to t\bar{t}\tilde\chi^{0}_{1}$) (horizontal) for Gtb models with $m(\tilde\chi^{0}_{1}) = 1$~TeV, obtained from the CC analysis.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Observed exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Expected exclusion limit in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb models obtained from the CC analysis. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm 1 \sigma$ of its theoretical uncertainty.

Upper limit at 95\% CL on the cross-section times branching ratio (fb) in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb (right) models obtained from the CC analysis. The numbers give the observed 95\% CL upper limit on the cross section in fb, with the label colour matching the associated best-expected region. Only a lower limit on the excluded cross section (>0.7 fb) is given at some points due to the very small number events expected and observed in the chosen SR. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background theoretical uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm1\sigma$ of its theoretical uncertainty.

Upper limit at 95\% CL on the cross-section times branching ratio (fb) in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gbb (right) models obtained from the NN analysis. The numbers give the observed 95\% CL upper limit on the cross section in fb, with the label colour matching the associated best-expected region. Only a lower limit on the excluded cross section (>0.7 fb) is given at some points due to the very small number events expected and observed in the chosen SR. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background theoretical uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm1\sigma$ of its theoretical uncertainty.

Upper limit at 95\% CL on the cross-section times branching ratio (fb) in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt (left) models obtained from the CC analysis. The numbers give the observed 95\% CL upper limit on the cross section in fb, with the label colour matching the associated best-expected region. Only a lower limit on the excluded cross section (>0.7 fb) is given at some points due to the very small number events expected and observed in the chosen SR. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background theoretical uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm1\sigma$ of its theoretical uncertainty.

Upper limit at 95\% CL on the cross-section times branching ratio (fb) in the $\tilde{g}$--$\tilde\chi^0_1$ mass plane for the Gtt (left) models obtained from the NN analysis. The numbers give the observed 95\% CL upper limit on the cross section in fb, with the label colour matching the associated best-expected region. Only a lower limit on the excluded cross section (>0.7 fb) is given at some points due to the very small number events expected and observed in the chosen SR. The dashed and solid bold lines show the 95\% CL expected and observed limits, respectively. The shaded bands around the expected limits show the impact of the experimental and background theoretical uncertainties. The dotted lines show the impact on the observed limit of the variation of the nominal signal cross-section by $\pm1\sigma$ of its theoretical uncertainty.

Acceptance for SR-Gtt-0L-B and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-0L-B and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-0L-M1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-0L-M1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-0L-M2 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-0L-M2 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-0L-C and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-0L-C and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-1L-B and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-1L-B and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-1L-M1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-1L-M1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-1L-M2 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-1L-M2 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-1L-C and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-1L-C and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-B and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-B and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-M and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-M and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-C and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-C and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-2100-1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-2100-1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-1800-1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-1800-1 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-2300-1200 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-2300-1200 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gtt-1900-1400 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gtt-1900-1400 and the $\tilde{g}\rightarrow t\bar{t}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-2800-1400 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-2800-1400 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-2300-1000 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-2300-1000 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-2100-1600 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-2100-1600 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Acceptance for SR-Gbb-2000-1800 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Efficiency for SR-Gbb-2000-1800 and the $\tilde{g}\rightarrow b\bar{b}\tilde\chi^0_1$ signal process.

Cutflow for the SR-Gtt-0L-B for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-0L-M1 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-0L-M2 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-0L-C for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-1L-B for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-1L-M1 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-1L-M2 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-1L-C for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-B for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-M for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-C for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtb-B for a representative Gtb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtb-M for a representative Gtb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtb-C for a representative Gtb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-2100-1 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-1800-1 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-2300-1200 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gtt-1900-1400 for a representative Gtt signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-2800-1400 for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-2300-1000 for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-2100-1600 for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Cutflow for the SR-Gbb-2000-1800 for a representative Gbb signal. Signal was generated with 30000 events. Expected yields are normalised to a luminosity of 139~fb$^{-1}$.

Expected 95% CL exclusion contours in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for boosted.

Observed 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for boosted.

Expected 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for boosted.

Observed 95% CL exclusion contours for the assumpition of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for boosted.

Expected 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for boosted.

Observed 95% CL exclusion contours for the assumption of Majorana type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for resolved.

Expected 95% CL exclusion contours for the assumption of Majorana type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for resolved.

Observed 95% CL exclusion contours for the assumption of Majorana type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for resolved.

Expected 95% CL exclusion contours for the assumption of Majorana type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for resolved.

Observed 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for resolved.

Expected 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for resolved.

Observed 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for resolved.

Expected 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for resolved.

The $m_{eJ}$ distribution in the boosted one electron channel.

The $m_{eeJ}$ distribution in the boosted two electron channel.

The $m_{\mu\mu J}$ distribution in the boosted two muon channel.

The $m_{jj}$ distribution in the resolved electron channel.

The $m_{jj}$ distribution in the resolved muon channel.

The $m_{eejj}$ distribution in the resolved electron channel.

The $m_{\mu\mu jj}$ distribution in the resolved muon channel.

The $h_{T}$ distribution in the resolved electron channel.

The $h_{T}$ distribution in the resolved muon channel.

Acceptance $\times$ efficiency for bSR1e.

Acceptance $\times$ efficiency for bSR2e.

Acceptance $\times$ efficiency for bSR2mu.

Acceptance $\times$ efficiency for rSROS2e.

Acceptance $\times$ efficiency for rSROS2mu.

Acceptance $\times$ efficiency for rSRSS2e.

Acceptance $\times$ efficiency for rSRSS2mu.

Expected and observed 95% CL upper limits for $\sigma(pp\to W_{R})\times B(W_{R}\to eeqq)$ as a function of $m(W_{R})$ for the special case of $m(W_{R}) = 2\times m(N_{R})$. The theoretical cross section corresponds to either the Dirac or Majorana neutrino interpretation by ignoring the charge of the leptons.

Expected and observed 95% CL upper limits for $\sigma(pp\to W_{R})\times B(W_{R} \to \mu\mu qq)$ as a function of $m(W_{R})$ for the special case of $m(W_{R}) = 2\times m(N_{R})$. The theoretical cross section corresponds to either the Dirac or Majorana neutrino interpretation by ignoring the charge of the leptons.

Expected and observed 95% CL upper limits for $\sigma(pp\to W_{R})\times B(W_{R} \to eeqq)$ as a function of $m(W_{R})$ for the special case of $m(N_{R})=200$ GeV. The theoretical cross section corresponds to either the Dirac or Majorana neutrino interpretation by ignoring the charge of the leptons.

Expected and observed 95% CL upper limits for $\sigma(pp\to W_{R})\times B(W_{R} \to \mu \mu qq)$ as a function of $m(W_{R})$ for the special case of $m(N_{R})=200$ GeV. The theoretical cross section corresponds to either the Dirac or Majorana neutrino interpretation by ignoring the charge of the leptons.

Expected and observed 95% CL upper limits for $\sigma(pp \to W_{R}) \times B(W_{R} \to eeqq)$ as a function of $m(W_{R})$ for the Majorana neutrino interpretations for the special case of $m(W_{R}) = 2\times m(N_{R})$.

Expected and observed 95% CL upper limits for $\sigma(pp \to W_{R}) \times B(W_{R} \to \mu\mu qq)$ as a function of $m(W_{R})$ for the Majorana neutrino interpretations for the special case of $m(W_{R}) = 2 \times m(N_{R})$.

Expected and observed 95% CL upper limits for $\sigma(pp \to W_{R}) \times B(W_{R} \to eeqq)$ as a function of $m(W_{R})$ for the Dirac neutrino interpretations for the special case of $m(W_{R}) = 2 \times m(N_{R})$.

Expected and observed 95% CL upper limits for $\sigma(pp \to W_{R}) \times B(W_{R} \to \mu \mu qq)$ as a function of $m(W_{R})$ for the Dirac neutrino interpretations for the special case of $m(W_{R}) = 2 \times m(N_{R})$.

Observed 95% CL upper limits on the cross section for all masses and charges of photon-fusion pair-produced spin-1/2 monopoles as a function of mass for magnetic charges $|g|=1g_D$ and $|g|=2g_D$.

Observed 95% CL upper limits on the cross section for all masses and charges of Drell-Yan spin-0 HECOs production as a function of mass for various values of electric charge in the range $20\le|z|\le100$.

Observed 95% CL upper limits on the cross section for all masses and charges of Drell-Yan spin-1/2 HECOs production as a function of mass for various values of electric charge in the range $20\le|z|\le100$.

Observed 95% CL upper limits on the cross section for all masses and charges of photon-fusion pair-produced spin-0 HECOs as a function of mass for various values of electric charge in the range $20\le|z|\le100$.

Observed 95% CL upper limits on the cross section for all masses and charges of photon-fusion pair-produced spin-1/2 HECOs as a function of mass for various values of electric charge in the range $20\le|z|\le100$.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 4000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 4000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 4000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 4000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 4000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 4000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 4000 GeV.

Predicted and observed yields as a function of $m_{4\ell}$ in the VBS-Suppressed region. Overflow events are included in the last bin of the distribution.

Relative systematic uncertainties in the unfolded differential cross section for 4ljj production in the VBS-Enhanced region as a function of m_{jj}

Relative systematic uncertainties in the unfolded differential cross section for 4ljj production in the VBS-Enhanced region as a function of m_{4\ell}

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $m_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $m_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $m_{4\ell}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $m_{4\ell}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $ p_{T,4l}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $ p_{T,4l}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $\Delta\phi_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $\Delta\phi_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $|\Delta y_{jj}|$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $|\Delta y_{jj}|$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $cos(\theta^{*}_{12})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $cos(\theta^{*}_{12})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $cos(\theta^{*}_{34})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $cos(\theta^{*}_{34})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $p_{T,jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $p_{T,jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $p_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $p_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $S_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $S_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Statistical correlation between different bins of the unfolded cross-section for all the measured observables in the VBS-Enhanced region.

Statistical correlation between different bins of the unfolded cross-section for all the measured observables in the VBS-Suppressed region.

Breakdown of systematic uncertainties in the measured inclusive cross-section. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Data bootstrap post unfolding for the inclusive cross-section. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb].

Differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The last bin of the distribution contains the overflow.

Differential absolute cross-section for $|\eta^{l}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $|\eta^{l}|$.

Data bootstrap post unfolding for the differential absolute cross-section for $|\eta^{l}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb].

Differential absolute cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential absolute cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The last bin of the distribution contains the overflow.

Differential absolute cross-section for $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $m^{e\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential absolute cross-section for $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The last bin of the distribution contains the overflow.

Differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The last bin of the distribution contains the overflow.

Differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The last bin of the distribution contains the overflow.

Differential absolute cross-section for $|\Delta\phi^{e\mu}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $|\Delta\phi^{e\mu}|$.

Data bootstrap post unfolding for the differential absolute cross-section for $|\Delta\phi^{e\mu}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/rad].

Differential absolute cross-section for $|y^{e\mu}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $|y^{e\mu}|$.

Data bootstrap post unfolding for the differential absolute cross-section for $|y^{e\mu}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb].

Double differential absolute cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest invariant mass bin of the two distributions contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential absolute cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The highest invariant mass bin of the two distributions contains the overflows.

Data bootstrap post unfolding for the double differential absolute cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The highest invariant mass bin of the two distributions contains the overflows.

Double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest invariant mass bin of the two distributions contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The highest invariant mass bin of the two distributions contains the overflows.

Data bootstrap post unfolding for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV/rad]. The highest invariant mass bin of the two distributions contains the overflows.

Double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Data bootstrap post unfolding for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV/rad]. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Data bootstrap post unfolding for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV/rad]. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV]. The last bin of the distribution contains the overflow.

Differential normalised cross-section for $|\eta^{l}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $|\eta^{l}|$.

Data bootstrap post unfolding for the differential normalised cross-section for $|\eta^{l}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry.

Differential normalised cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential normalised cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry. The last bin of the distribution contains the overflow.

Differential normalised cross-section for $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $m^{e\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential normalised cross-section for $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry. The last bin of the distribution contains the overflow.

Differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV]. The last bin of the distribution contains the overflow.

Differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV]. The last bin of the distribution contains the overflow.

Differential normalised cross-section for $|\Delta\phi^{e\mu}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $|\Delta\phi^{e\mu}|$.

Data bootstrap post unfolding for the differential normalised cross-section for $|\Delta\phi^{e\mu}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/rad].

Differential normalised cross-section for $|y^{e\mu}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $|y^{e\mu}|$.

Data bootstrap post unfolding for the differential normalised cross-section for $|y^{e\mu}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry.

Double differential normalised cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest invariant mass bin of the two distributions contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential normalised cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The highest invariant mass bin of the two distributions contains the overflows.

Data bootstrap post unfolding for the double differential normalised cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV]. The highest invariant mass bin of the two distributions contains the overflows.

Double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest invariant mass bin of the two distributions contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The highest invariant mass bin of the two distributions contains the overflows.

Data bootstrap post unfolding for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV/rad]. The highest invariant mass bin of the two distributions contains the overflows.

Double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Data bootstrap post unfolding for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV/rad]. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Data bootstrap post unfolding for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV/rad]. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.