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

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

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

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

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

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

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

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

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

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

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

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

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

The rho-parameter, i.e. the ratio of the real to imaginary part of the elastic scattering amplitude extrapolated to t=0. The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter B from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter B from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter C from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter C from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter D from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter D from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The total elastic cross section measured inside the fiducial volume. The systematic uncertainty includes experimental uncertainties.

The total elastic cross section measured inside the fiducial volume. The systematic uncertainty includes experimental uncertainties.

The total elastic cross section obtained from the fitted parameters, extrapolated to full phase space using only the nuclear amplitude.

The total elastic cross section obtained from the fitted parameters, extrapolated to full phase space using only the nuclear amplitude.

The total inelastic cross section.

The total inelastic cross section.

The ratio of elastic to total cross section.

The ratio of elastic to total cross section.

The differential elastic cross section as function of t with statistical and systematic uncertainties. The systematic uncertainties are given as signed relative change for 20 sources of experimental uncertainty associated to nuisance parameters used in the fit for the extraction of physics parameters.

The differential elastic cross section as function of t with statistical and systematic uncertainties. The systematic uncertainties are given as signed relative change for 20 sources of experimental uncertainty associated to nuisance parameters used in the fit for the extraction of physics parameters.

Statistical covariance matrix for the measurement of the differential elastic cross section as function of t.

Statistical covariance matrix for the measurement of the differential elastic cross section as function of t.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Observed excluded cross-section at 95% CL for the WZ-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from Fig 8(aux).

Expected exclusion limits at 95% CL for the WZ-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Observed exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from from Fig 13(a) and from Fig 7 and Fig 10(aux).

Observed excluded cross-section at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from Fig 7(aux) and Fig 10(aux).

positive one $\sigma$ observed exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from from Fig 13(a) and from Fig 7 and Fig 10(aux).

negative one $\sigma$ observed exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from from Fig 13(a) and from Fig 7 and Fig 10(aux).

Expected exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Expected exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Expected exclusion limits at 95% CL for the Wh-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{bRPV}_{2l-SS}$. in a susy scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV, tan$\beta$=5. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{bRPV}_{3l}$. in a susy scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV, tan$\beta$=5. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{WZ}_{high-m_{T2}}$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 150 GeV, $m(\tilde{\chi}^{0} _{1})$ = 50 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region $SR^{WZ}_{low-m_{T2}}$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 150 GeV, $m(\tilde{\chi}^{0} _{1})$ = 50 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the low mass $SR^{RPV}_{2l1b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the medium mass $SR^{RPV}_{2l1b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the low mass $SR^{RPV}_{2l2b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the medium mass $SR^{RPV}_{2l2b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the high mass $SR^{RPV}_{2l2b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the low mass $SR^{RPV}_{2l3b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the medium mass $SR^{RPV}_{2l3b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the high mass $SR^{RPV}_{2l3b}$, where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{0} _{1}/\tilde{\chi}^{0} _{2})$ = 200 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the $SR^{Wh}_{low-m_{T2} }$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Higgs bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 300 GeV, $m(\tilde{\chi}^{0} _{1})$ = 100 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the $SR^{Wh}_{high-m_{T2} }$. The wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Higgs bosons. The masses of the superpartners involved in the process are set to $m(\tilde{\chi}^{\pm} _{1}/\tilde{\chi}^{0} _{2})$ = 300 GeV, $m(\tilde{\chi}^{0} _{1})$ = 100 GeV. Only statistical uncertainties are shown.

Signal Hepdataeptance for $SR^{bRPV}_{2l-SS}$ signal region from Fig 13(a)(aux) in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal Hepdataeptance for $SR^{bRPV}_{3l}$ signal region from Fig 13(b)(aux) in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal acceptance for $SR^{WZ}_{high-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal acceptance for $SR^{WZ}_{low-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal acceptance for $SR^{RPV}_{2l1b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l1b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l2b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l2b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l2b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l3b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l3b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{RPV}_{2l3b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{bRPV}_{2l-SS}$ signal region in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal efficiency for $SR^{bRPV}_{3l}$ signal region in a SUSY scenario where $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are produced in pairs and decay to all possible allowed bRPV decays.

Signal efficiency for $SR^{WZ}_{high-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal efficiency for $SR^{WZ}_{low-m_{T2}}$ in a SUSY scenario where the wino-like doublet pair ($\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$) were produced and then decays into $bino-like \tilde{\chi}^{0} _{1}$ which is the lightest SUSY particle (LSP) accompanied by mass on-shell or mass off-shell W and Z bosons.

Signal efficiency for $SR^{RPV}_{2l1b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l1b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l2b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l2b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l2b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l3b}-L$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l3b}-M$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal efficiency for $SR^{RPV}_{2l3b}-H$ signal region in a SUSY scenario where the $\tilde{\chi}^{0} _{1} and \tilde{\chi}^{0} _{2}$ are directly produced and undergoes prompt RPV decays.

Signal acceptance for $SR^{Wh}_{high-m_{T2} }$ signal region from Fig 11(a)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Signal acceptance for $SR^{Wh}_{low-m_{T2} }$ signal region from Fig 11(b)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Signal efficiency for $SR^{Wh}_{high-m_{T2} }$ signal region from Fig 15(a)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Signal efficiency for $SR^{Wh}_{low-m_{T2} }$ signal region from Fig 15(b)(aux) in a SUSY scenario where direct production of a lightest $\tilde{\chi}^{\pm} _{1} and \tilde{\chi}^{0} _{2}$ , decay with 100% branching ratio to a final state with a same sign light lepton (e or $\mu$) pair and two lightest neutralino1, via the on-shell emission of SM W and Higgs bosons,

Observed 95% X-section upper limits as a function of higgsino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{1}/\tilde{\chi}^{0}_{2}$ mass in the bilinear RPV model from Fig 14.

Observed 95% X-section upper limits as a function of higgsino $\tilde{\chi}^{0}_{1}/\tilde{\chi}^{0}_{2}$ mass in the UDD RPV model from Fig 18.

Observed 95% X-section upper limits as a function of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ mass in the WZ-mediated simplified model of wino $\tilde{\chi}^{\pm}_{1}/\tilde{\chi}^{0}_{2}$ production from Fig 9(aux).

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{WZ}_{high-m_{T2}}$ from publication's Figure 11(a) . The last bin is inclusive.

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{WZ}_{low-m_{T2}}$ from publication's Figure 11(b) . The last bin is inclusive.

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{bRPV}_{2l-SS}$ from publication's Figure 11(c) . The last bin is inclusive.

N-1 distributions for $m_{T2}$ of observed data and expected background towards $SR^{bRPV}_{3l}$ from publication's Figure 11(d) . The last bin is inclusive.

N-1 distributions for $\sum p^{b-jet}_{T}/\sum p^{jet}_{T}$ of observed data and expected background towards $SR^{RPV}_{2l1b}-L$ from publication's Figure 16(a) . The last bin is inclusive.

N-1 distributions for $\sum p^{b-jet}_{T}/\sum p^{jet}_{T}$ of observed data and expected background towards $SR^{RPV}_{2l2b}-M$ from publication's Figure 16(b) . The last bin is inclusive.

N-1 distributions for $\sum p^{b-jet}_{T}/\sum p^{jet}_{T}$ of observed data and expected background towards $SR^{RPV}_{2l3b}-H$ from publication's Figure 16(c) . The last bin is inclusive.

N-1 distribution for $E_{T}^{miss}$ in $SR^{Wh}_{high-m_{T2} }$ in ee channel

N-1 distribution for $E_{T}^{miss}$ in $SR^{Wh}_{high-m_{T2} }$ in e$\mu$ channel

N-1 distribution for $E_{T}^{miss}$ in $SR^{Wh}_{high-m_{T2} }$ in $\mu\mu$ channel

N-1 distribution for $\mathcal{S}(E_{T}^{miss})$ in $SR^{Wh}_{low-m_{T2} }$ in ee channel

N-1 distribution for $\mathcal{S}(E_{T}^{miss})$ in $SR^{Wh}_{low-m_{T2} }$ in e$\mu$ channel

N-1 distribution for $\mathcal{S}(E_{T}^{miss})$ in $SR^{Wh}_{low-m_{T2} }$ in $\mu\mu$ channel

Summary of results for cross-section of non-prompt $\psi(2S)$ decaying to a muon pair for 13 TeV data in fb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $J/\psi$ decaying to a muon pair as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $\psi(2S)$ decaying to a muon pair as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Correction factors for an assumed angular dependence $\propto 1+\lambda_{\theta} \cos^2\theta $ in the helicity frame, for several values of the parameter $\lambda_{\theta}$. The correction factors were found to be the same (within 1% - 2%) for $J\psi$ and $\psi(2S)$ mesons, for prompt and non-prompt production mechanisms, and for the three rapidity intervals considered in this paper.

Expected exclusion limits at 95% CL from Fig 7(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Expected exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Expected exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRGGWZ-H.

N-1 distribution for $E_{\mathrm{T}}^{\mathrm{miss}}$of observed data and expected background in SRGGSlep-M.

N-1 distribution for $\sum{p_{\mathrm{T}}^\mathrm{jet}}$of observed data and expected background in SRUDD-ge2b.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRLQD.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRSSWZ-H.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRSSSlep-H(loose).

Signal acceptance for SRGGWZ-H signal region from Fig 10(c) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-H signal region from Fig 15(c) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGWZ-M signal region from Fig 10(b) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-M signal region from Fig 15(b) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGWZ-L signal region from Fig 10(a) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-L signal region from Fig 15(a) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-L signal region from Fig 12(a) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-L signal region from Fig 17(a) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-M signal region from Fig 12(b) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-M signal region from Fig 17(b) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-H signal region from Fig 12(c) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-H signal region from Fig 17(c) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRUDD-1b signal region from Fig 14(b) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-1b signal region from Fig 19(b) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-2b signal region from Fig 14(c) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-2b signal region from Fig 19(c) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-ge2b signal region from Fig 14(d) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-ge2b signal region from Fig 19(d) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-ge3b signal region from Fig 14(e) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-ge3b signal region from Fig 19(e) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRLQD signal region from Fig 14(a) in a SUSY scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Signal efficiency for SRLQD signal region from Fig 19(a) in a SUSY scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Signal acceptance for SRSSWZ-L signal region from Fig 11(a) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-L signal region from Fig 16(a) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-ML signal region from Fig 11(b) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-ML signal region from Fig 16(b) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-MH signal region from Fig 11(c) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-MH signal region from Fig 16(c) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-H signal region from Fig 11(d) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-H signal region from Fig 16(d) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-H signal region from Fig 13(d) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-H signal region from Fig 18(d) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-MH signal region from Fig 13(c) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-MH signal region from Fig 18(c) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-L signal region from Fig 13(a) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-L signal region from Fig 18(a) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-ML signal region from Fig 13(b) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-ML signal region from Fig 18(b) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-H(loose) signal region from Fig 13(e) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-H(loose) signal region from Fig 18(e) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGWZ-H in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGWZ-M in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGWZ-L in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGSlep-L in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGSlep-M in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGSlep-H in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-1b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-2b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-ge2b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-ge3b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRLQD in a susy scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2200 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1870 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-L in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-ML in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-MH in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-H in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-H in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-MH in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-L in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-ML in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-H(loose) in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Cross-section upper limits at 95% CL from Fig1(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Cross-section upper limits at 95% CL from Fig1(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Cross-section upper limits at 95% CL from Fig1(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Expected 95% CL exclusion contour for the Gbb signal.

Observed 95% CL exclusion contour for the Gbb signal.

Expected 95% CL exclusion contour for the Gtt combination.

Observed 95% CL exclusion contour for the Gtt combination.

Acceptances for the Gbb model in SR-Gbb-A. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gbb model in SR-Gbb-B. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gbb model in SR-Gbb-C. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gtt model in SR-Gtt-0L-A. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gtt model in SR-Gtt-0L-B. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gtt model in SR-Gtt-0L-C. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gtt model in SR-Gtt-1L-A. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptances for the Gtt model in SR-Gtt-1L-B. Acceptance is evaluated at truth level, with only leptons from heavy bosons and taus considered, and no further quality or isolation criteria applied in their selection.

Acceptance times efficiency for the Gbb model in SR-Gbb-A.

Acceptance times efficiency for the Gbb model in SR-Gbb-B.

Acceptance times efficiency for the Gbb model in SR-Gbb-C.

Acceptance times efficiency for the Gtt model in SR-Gtt-0L-A.

Acceptance times efficiency for the Gtt model in SR-Gtt-0L-B.

Acceptance times efficiency for the Gtt model in SR-Gtt-0L-C.

Acceptance times efficiency for the Gtt model in SR-Gtt-1L-A.

Acceptance times efficiency for the Gtt model in SR-Gtt-1L-B.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gbb model in SR-Gbb-A.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gbb model in SR-Gbb-B.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gbb model in SR-Gbb-C.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gtt model in SR-Gtt-0L-A.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gtt model in SR-Gtt-0L-B.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gtt model in SR-Gtt-0L-C.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gtt model in SR-Gtt-1L-A.

95% CL upper limit on the cross-section times branching ratio (in fb) for the Gtt model in SR-Gtt-1L-B.

Signal region yielding the best expected sensitivity for each point of the parameter space in the Gbb model.

Signal region yielding the best expected sensitivity for each point of the parameter space in the Gtt model for the 0-lepton channel.

Signal region yielding the best expected sensitivity for each point of the parameter space in the Gtt model for the 1-lepton channel.

Combination of two 0-lepton and 1-lepton signal regions yielding the best expected sensitivity for each point of the parameter space in the Gtt model.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 10j50 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

Observed 95% CL limit for the pMSSM grid.

Observed 95% CL limit for the pMSSM grid when the signal cross section is increased by one standard deviation.

Observed 95% CL limit for the pMSSM grid when the signal cross section is decreased by one standard deviation.

Expected 95% CL limit for the pMSSM grid.

+1 sigma excursion of the expected 95% CL limit for the pMSSM grid.

-1 sigma excursion of the expected 95% CL limit for the pMSSM grid.

Observed 95% CL limit for the 2Step grid.

Observed 95% CL limit for the 2Step grid when the signal cross section is increased by one standard deviation.

Observed 95% CL limit for the 2Step grid when the signal cross section is decreased by one standard deviation.

Expected 95% CL limit for the 2Step grid.

+1 sigma excursion of the expected 95% CL limit for the 2Step grid.

-1 sigma excursion of the expected 95% CL limit for the 2Step grid.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j50 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j50 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j50 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 9j50 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 9j50 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 9j50 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 10j50 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 10j50 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 10j50 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 7j80 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 7j80 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 7j80 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region 8j80 2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

Degree of multijet closure for signal and vaidation regions with at no b-jet requirement. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The bins labelled in bold are signal regions, while the others are validation regions. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Degree of multijet closure for signal and vaidation regions with at least 1 b-jet. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The bins labelled in bold are signal regions, while the others are validation regions. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Degree of multijet closure for signal and vaidation regions with at least 2 b-jets. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The bins labelled in bold are signal regions, while the others are validation regions. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Summary of all 15 signal regions (post-fit).

Signal region yielding the best-expected CLs value, the best expected CLs value, and the corresponding observed CLs value for the 2Step grid.

Signal region yielding the best-expected CLs value, the best expected CLs value, and the corresponding observed CLs value for the pMSSM grid.

95% CLs observed upper limit on model cross-section for 2-step signal points for the best-expected signal region.

Performance of the 8j50-0b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j50-1b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j50-2b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-0b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-1b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-2b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-0b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-1b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-2b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-0b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-1b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-2b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-0b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-1b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-2b selection for the pMSSM grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j50-0b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j50-1b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j50-2b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-0b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-1b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 9j50-2b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-0b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-1b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 10j50-2b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-0b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-1b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 7j80-2b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-0b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-1b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Performance of the 8j80-2b selection for the 2Step grid: number of generated signal events; total signal cross-section; acceptance; efficiency (fractional); observed CL using this region alone; expected CL using this region alone.

Charged-particle multiplicities in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of transverse momentum for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicity distribution in proton-proton collisions at a centre-of-mass energy of 13000 GeV for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Average transverse momentum in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of the number of charged particles in the event for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Extrapolated charged-particle multiplicities in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of pseudorapidity for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Extrapolated charged-particle multiplicities in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of transverse momentum for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Extrapolated charged-particle multiplicity distributions in proton-proton collisions at a centre-of-mass energy of 13000 GeV for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Extrapolated average transverse momentum in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of the number of charged particles in the event for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of pseudorapidity for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <0.8.

Charged-particle multiplicities in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of transverse momentum for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <0.8.

Charged-particle multiplicity distribution in proton-proton collisions at a centre-of-mass energy of 13000 GeV for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <0.8.

Average transverse momentum in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of the number of charged particles in the event for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <0.8.

Extrapolated charged-particle multiplicities in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of pseudorapidity for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <0.8.

Extrapolated charged-particle multiplicities in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of transverse momentum for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <0.8.

Extrapolated charged-particle multiplicity distributions in proton-proton collisions at a centre-of-mass energy of 13000 GeV for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <0.8.

Extrapolated average transverse momentum in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of the number of charged particles in the event for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <0.8.

Comparison of data with estimated backgrounds in the $\tilde{E}_\text{T}^\text{miss}$ distribution with the TZCR1 event selection except for the requirement on $\tilde{E}_\text{T}^\text{miss}$. The variables $\tilde{E}_\text{T}^\text{miss}$ and $\tilde{m}_\text{T}$ are constructed in the same way as $E_\text{T}^\text{miss}$ and $m_\text{T}$ but treating the leading photon transverse momentum as invisible. The predicted backgrounds are scaled with normalization factors. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin includes overflow.

Comparison of data with estimated backgrounds in the $\tilde{m}_\text{T}$ distribution with the TZCR1 event selection except for the requirement on $\tilde{m}_\text{T}$. The variables $\tilde{E}_\text{T}^\text{miss}$ and $\tilde{m}_\text{T}$ are constructed in the same way as $E_\text{T}^\text{miss}$ and $m_\text{T}$ but treating the leading photon transverse momentum as invisible. The predicted backgrounds are scaled with normalization factors. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin includes overflow.

Comparison of the observed data ($n_\text{obs}$) with the predicted background ($n_\text{exp}$) in the validation and signal regions. The background predictions are obtained using the background-only fit configuration. The bottom panel shows the significance of the difference between data and predicted background, where the significance is based on the total uncertainty ($\sigma_\text{tot}$).

Jet multiplicity distributions for events where exactly two signal leptons are selected. No correction factors are included in the background normalizations. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin includes overflow.

Jet multiplicity distributions for events where exactly one lepton plus one $\tau$ candidate are selected. No correction factors are included in the background normalizations. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin includes overflow.

The $E_\text{T}^\text{miss}$ distribution in SR1. In the plot, the full event selection in the corresponding signal region is applied, except for the requirement on $E_\text{T}^\text{miss}$. The predicted backgrounds are scaled with normalization factors. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin contains the overflow. Benchmark signal models are overlaid for comparison. The benchmark models are specified by the gluino and stop masses, given in TeV in the table.

The $m_\text{T}$ distribution in SR1. In the plot, the full event selection in the corresponding signal region is applied, except for the requirement on $m_\text{T}$. The predicted backgrounds are scaled with normalization factors. The uncertainty band includes statistical and all experimental systematic uncertainties. The last bin contains the overflow. Benchmark signal models are overlaid for comparison. The benchmark models are specified by the gluino and stop masses, given in TeV in the table.

Expected (black dashed) 95% excluded regions in the plane of $m_{\tilde{g}}$ versus $m_{\tilde{t}_1}$ for gluino-mediated stop production.

Observed (red solid) 95% excluded regions in the plane of $m_{\tilde{g}}$ versus $m_{\tilde{t}_1}$ for gluino-mediated stop production.

Expected (black dashed) 95% excluded regions in the plane of $m_{\tilde{t}_1}$ versus $m_{\tilde{\chi}_1^0}$ for direct stop production.

Observed (red solid) 95% excluded regions in the plane of $m_{\tilde{t}_1}$ versus $m_{\tilde{\chi}_1^0}$ for direct stop production.

The expected upper limits on $T$ quark pair production times the squared branching ratio for $T \rightarrow tZ$ as a function of the $T$ quark mass.

The observed upper limits on $T$ quark pair production times the squared branching ratio for $T \rightarrow tZ$ as a function of the $T$ quark mass.

The expected limits on $T$ quarks as a function of the branching ratios $B\left(T \rightarrow bW\right)$ and $B\left(T \rightarrow tH\right)$ for a $T$ quark with a mass of 800 GeV. The $T$ is assumed to decay in three possible ways: $T \to tZ$, $T \to tH$, and $T \to bW$.

The observed limits on $T$ quarks as a function of the branching ratios $B\left(T \rightarrow bW\right)$ and $B\left(T \rightarrow tH\right)$ for a $T$ quark with a mass of 800 GeV. The $T$ is assumed to decay in three possible ways: $T \to tZ$, $T \to tH$, and $T \to bW$.

The $m_\text{T}$ distribution in the WVR2-tail validation region which has the same preselection and jet $p_\text{T}$ requirements as SR2.

The $am_\text{T2}$ distribution in the WVR2-tail validation region which has the same preselection and jet $p_\text{T}$ requirements as SR2.

Large-radius jet mass ($R=1.2$), decomposed into the number of small-radius jet constituents. The lower panel shows the ratio of the total data to the total prediction (summed over all jet multiplicities). Events are required to have one lepton, four jets with $p_\text{T}>80,50,40,40$ GeV, at least one $b$-tagged jet, $E_\text{T}^\text{miss}>200$ GeV, and $m_\text{T}>30$ GeV.

Distribution of $m_\text{T2}^\tau$ in data for a selection enriched in $t\bar{t}$ events with one hadronically decaying $\tau$. Events that have no hadronic $\tau$ candidate (that passes the Loose identification criteria, as well as other requirements) are not shown in the plot.

Upper limits on the model cross-section in units of pb for the gluino-mediated stop models.

Upper limits on the model cross-section in units of pb for the models with direct stop pair production.

Illustration of the best expected signal region per signal grid point for the gluino-mediated stop models. This mapping is used for the final combined exclusion limits.

Illustration of the best expected signal region per signal grid point for models with direct stop pair production. This mapping is used for the final combined exclusion limits.

Expected $CL_s$ values for the gluino-mediated stop models.

Observed $CL_s$ values for the gluino-mediated stop models.

Expected $CL_s$ values for the direct stop pair production models.

Observed $CL_s$ values for the direct stop pair production models.

Expected limit using SR1 for models with direct stop pair production and an unpolarized stop (and bino LSP).

Expected limit using SR1 for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Expected limit using SR1 for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Observed limit using SR1 for models with direct stop pair production and an unpolarized stop (and bino LSP).

Observed limit using SR1 for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Observed limit using SR1 for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Expected limit using SR2 for models with direct stop pair production and an unpolarized stop (and bino LSP).

Expected limit using SR2 for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Expected limit using SR2 for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Observed limit using SR2 for models with direct stop pair production and an unpolarized stop (and bino LSP).

Observed limit using SR2 for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Observed limit using SR2 for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Expected limit using SR1+SR2 (best expected) for models with direct stop pair production and an unpolarized stop (and bino LSP).

Expected limit using SR1+SR2 (best expected) for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Expected limit using SR1+SR2 (best expected) for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Observed limit using SR1+SR2 (best expected) for models with direct stop pair production and an unpolarized stop (and bino LSP).

Observed limit using SR1+SR2 (best expected) for models with direct stop pair production with $\tilde{t}_1=\tilde{t}_L$ (and bino LSP).

Observed limit using SR1+SR2 (best expected) for models with direct stop pair production with $\tilde{t}_1\sim\tilde{t}_R$ (and bino LSP).

Acceptance for SR1 in the gluino-mediated stop models. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance for SR1 in the direct stop pair production. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance for SR2 in the gluino-mediated stop models. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance for SR2 in the direct stop pair production. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance for SR3 in the gluino-mediated stop models. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Acceptance for SR3 in the direct stop pair production. The acceptance is defined as the fraction of signal events that pass the analysis selection performed on generator-level objects, therefore emulating an ideal detector with perfect particle identification and no measurement resolution effects.

Efficiency for SR1 in the gluino-mediated stop models. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency for SR1 in the direct stop pair production. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency for SR2 in the gluino-mediated stop models. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency for SR2 in the direct stop pair production. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency for SR3 in the gluino-mediated stop models. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).

Efficiency for SR3 in the direct stop pair production. The efficiency is the ratio between the expected signal rate calculated with simulated data passing all the reconstruction level cuts applied to reconstructed objects, and the signal rate for an ideal detector (with perfect particle identification and no measurement resolution effects).