Total efficiencies for the selection of signal events after categorization and application of the tighter photon $E_{\mathrm{T}}^{\gamma}$ selection used to optimize the signal significance spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$. In addition to the total efficiency, contributions to the signal selection from each of the separate event categories are shown. The efficiencies calculated from MC samples with $W/Z$ hadronic decays are shown as the points on each curve. The line presents interpolated results.

Total efficiencies for the selection of signal events after categorization and application of the tighter photon $E_{\mathrm{T}}^{\gamma}$ selection used to optimize the signal significance spin-2 $q\bar{q}\to X^0 \to Z\gamma$. In addition to the total efficiency, contributions to the signal selection from each of the separate event categories are shown. The efficiencies calculated from MC samples with $W/Z$ hadronic decays are shown as the points on each curve. The line presents interpolated results.

Total efficiencies for the selection of signal events after categorization and application of the tighter photon $E_{\mathrm{T}}^{\gamma}$ selection used to optimize the signal significance spin-2 $gg\to X^0 \to Z\gamma$. In addition to the total efficiency, contributions to the signal selection from each of the separate event categories are shown. The efficiencies calculated from MC samples with $W/Z$ hadronic decays are shown as the points on each curve. The line presents interpolated results.

The $m_{J\gamma}$ distributions of data events selected for the spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$ search in the WMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-0 $gg \to X^{0} \to Z\gamma$ search in the BTAG category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-0 $gg \to X^{0} \to Z\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-0 $gg \to X^{0} \to Z\gamma$ search in the ZMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $gg \to X^{0} \to Z\gamma$ search in the BTAG category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $gg \to X^{0} \to Z\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $gg \to X^{0} \to Z\gamma$ search in the ZMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $q\bar{q} \to X^{0} \to Z\gamma$ search in the BTAG category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $q\bar{q} \to X^{0} \to Z\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $q\bar{q} \to X^{0} \to Z\gamma$ search in the ZMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-0 $gg \to X^{0} \to Z\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-2 $gg \to X^{0} \to Z\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-2 $q\bar{q} \to X^{0} \to Z\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The jet mass distribution of large-$R$ jets originating from the hadronic decay of $W$ and $Z$ bosons produced from the decay of BSM bosons with mass $m_X = 2000$ GeV. The decays simulated are for the production models $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$ with a spin-1 resonance $X^{\pm}$ and $gg\to X^0 \to Z\gamma$ with a spin-0 resonance $X^{0}$.

The photon $E_{T}$ selection criteria as a function of resonance mass for D2 categories. Selections for four different signal channels are shown together.

The photon $E_{T}$ selection criteria as a function of resonance mass for ZMASS or WMASS categories. Selections for four different signal channels are shown together.

Signal selection efficiency ($\epsilon$) times acceptance ($A$) as a function of $H^{\pm}$. The estimated $\epsilon\times A$ arises from the lepton selection and triggering (∼30%) as well as jet selection and flavour tagging (∼10% or lower). The decrease of $\epsilon\times A$ for $m_{H^{\pm}}$ = 120 GeV and higher is expected from the kinematic constraint on the $H^{\pm}$ decay products due to the top-quark mass.

The performance of jet energy resolution (JER) for jets with |y| < 2.1 evaluated as a function of pT_truth in different centrality bins. Simulated hard scatter events were overlaid onto events from a dedicated sample of minimum-bias Xe+Xe data. The fit parameters are listed in a sperate table (Extras 1)

The relative magnitude of systematic uncertainties for per-pair normalized xJ distributions in 0-10% Xe+Xe centrality

The relative magnitude of systematic uncertainties for absolutely normalized xJ distributions in 0-10% Xe+Xe centrality

The relative magnitude of systematic uncertainties for rho distributions for leading jets in 0-10% Xe+Xe centrality

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

Parameter a,b, and c from JER fits in Figure 1b.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Distribution of the NN output in the $\geq$1fj$\,$SR in data and the expected contribution of the signal and background processes after the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions considering the correlations of the uncertainties as obtained by the fit.

Distribution of the NN output in the $t\bar{t}\gamma\,$CR in data and the expected contribution of the signal and background processes after the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions considering the correlations of the uncertainties as obtained by the fit.

Total event yield in the $W\gamma\,$CR in data and the expected contribution of the signal and background processes after the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions considering the correlations of the uncertainties as obtained by the fit.

Distribution of the scalar sum of the jet transverse momenta in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions. The first and last bins include the underflow and overflow, respectively.

Distribution of the $\eta$ of the $b$-tagged jet in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions.

Distribution of the reconstructed top-quark mass in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions. The first and last bins include the underflow and overflow, respectively.

Distribution of the $p_{\mathrm{T}}$ of the top-quark+photon system in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the photon $p_{\mathrm{T}}$ in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the photon $\eta$ in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions.

Distribution of the scalar sum of the jet transverse momenta in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the invariant mass of the $b$-tagged jet and the highest-$p_{\mathrm{T}}$ forward jet in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of $p_{\mathrm{T}}$ of the highest-$p_{\mathrm{T}}$ forward jet in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the difference in $\eta$ between the highest-$p_{\mathrm{T}}$ forward jet and the photon in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the energy of the system formed by the highest-$p_{\mathrm{T}}$ forward jet and the photon in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions. The first and last bins include the underflow and overflow, respectively.

Distribution of the $\eta$ of the $b$-tagged jet in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions.

Distribution of the reconstructed top-quark mass in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions. The first and last bins include the underflow and overflow, respectively.

Distribution of the $p_{\mathrm{T}}$ of the top-quark and photon system in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the photon $p_{\mathrm{T}}$ in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the photon $\eta$ in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions.

Ordered list of the 30 systematic uncertainties with the largest impact on the measured signal normalisation in the fit to data in the parton-level measurement considered as nuisance parameters (NPs) in the profile-likelihood fit. The column "NP value, error" corresponds to the nominal best-fit values and the corresponding uncertainties. The impact of each NP, $\Delta\sigma$/$\sigma_{\mathrm{pred}}$, is computed by comparing the nominal best-fit value of the POI ($\sigma$/$\sigma_{\mathrm{pred}}$) with the result of the fits when fixing the considered NP to its best-fit value shifted by its pre-fit and post-fit uncertainties. The corresponding impacts are listed in the "POI impact prefit high/low" and "POI impact high/low" columns, respectively. The "MC stat." NPs represent the MC statistical uncertainty and they enter the likelihood with a Poisson term, while all the other NPs enter the likelihood via a Gaussian term.

Ordered list of the 30 systematic uncertainties with the largest impact on the measured signal normalisation in the fit to data in the particle-level measurement considered as nuisance parameters (NPs) in the profile-likelihood fit. The column "NP value, error" corresponds to the nominal best-fit values and the corresponding uncertainties. The impact of each NP, $\Delta\sigma$/$\sigma_{\mathrm{pred}}$, is computed by comparing the nominal best-fit value of the POI ($\sigma$/$\sigma_{\mathrm{pred}}$) with the result of the fits when fixing the considered NP to its best-fit value shifted by its pre-fit and post-fit uncertainties. The corresponding impacts are listed in the "POI impact prefit high/low" and "POI impact high/low" columns, respectively. The "MC stat." NPs represent the MC statistical uncertainty and they enter the likelihood with a Poisson term, while all the other NPs enter the likelihood via a Gaussian term.

This table lists the kinematic requirements on parton-level objects used to define of the fiducial phase space for the parton-level measurement. Frixione isolation ($\href{https://arxiv.org/abs/hep-ph/9801442}{\text{hep-ph/9801442}}$) with a chosen radius of $\Delta R = 0.2$ is applied to photons ($\gamma$). The measured fiducial parton-level cross section is $\sigma_{tq\gamma}\times\mathcal{B}(t\rightarrow l\nu b) = 688\pm 23(\text{stat.})^{+75}_{-71}(\text{syst.})\,$fb.

This table lists the kinematic requirements on particle-level objects used to define of the fiducial phase space for the particle-level measurement. The particle level objects are photons ($\gamma$) not from a hadron decay, neutrinos not from a hadron decay ($\nu$), prompt electrons and muons ($\ell$) "dressed" by adding close-by ($\Delta R < 0.1$) photons, and anti-$k_t$ $R = 0.4$ jets built from stable particles ($\tau > 30\,$ps) and tau leptons excluding neutrinos and prompt dressed muons. Jets are $b$-tagged ($b$-jet) using ghost-matched $b$-hadrons with $p_{\text{T}} > 5\,$GeV. Apart from the kinematic requirements, isolation and overlap removal criteria are applied. Jets within $\Delta R = 0.4$ of a photon are removed if the $p_{\text{T}}$ of charged particles within $\Delta R = 0.3$ of the photon is smaller than $10\,\%$ of its $p_{\text{T}}$. Jets within $\Delta R = 0.4$ of a lepton are removed. Events are removed where a photon is close ($\Delta R < 0.4$) to a lepton or a surviving jet. The measured fiducial particle-level cross section is $\sigma_{tq\gamma}\times\mathcal{B}(t\rightarrow l\nu b)+\sigma_{t(\rightarrow l\nu b\gamma)q} = 303\pm 9(\text{stat.})^{+33}_{-32}(\text{syst.})\,$fb.

Comparison between data and background prediction before the fit (Pre-Fit) for the mass of the FCNC top-quark candidate in SR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The four FCNC LH signals are also shown separately, normalized to five times the cross-section corresponding to the most stringent observed branching ratio limits. The first (last) bin in all distributions includes the underflow (overflow). The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction before the fit (Pre-Fit) for the mass of the SM top-quark candidate in SR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The four FCNC LH signals are also shown separately, normalized to five times the cross-section corresponding to the most stringent observed branching ratio limits. The first (last) bin in all distributions includes the underflow (overflow). The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction before the fit (Pre-Fit) for the transverse momentum of the Z boson candidate in SR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The four FCNC LH signals are also shown separately, normalized to five times the cross-section corresponding to the most stringent observed branching ratio limits. The first (last) bin in all distributions includes the underflow (overflow). The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{1}$ discriminant in the mass sideband CR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also separately shown, normalized to 500 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{2}^{u}$ discriminant in the mass sideband CR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also separately shown, normalized to 500 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{1}$ discriminant in SR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also separately shown, normalized to 50 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{2}^{u}$ discriminant in SR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also separately shown, normalized to 50 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZc LH coupling extraction. The distribution is for the $D_{1}$ discriminant in the mass sideband CR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZc LH signals are also separately shown, normalized to 500 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZc LH coupling extraction. The distribution is for the $D_{2}^{c}$ discriminant in the mass sideband CR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZc LH signals are also separately shown, normalized to 500 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZc LH coupling extraction. The distribution is for the $D_{1}$ discriminant in SR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZc LH signals are also separately shown, normalized to 50 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZc LH coupling extraction. The distribution is for the $D_{2}^{c}$ discriminant in SR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZc LH signals are also separately shown, normalized to 50 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the leading lepton $p_{T}$ in the $t\bar{t}$ CR. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also shown separately, normalized to $10^{3}$ times the best fit of the signal yield. The last bin includes the overflow. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the third lepton $p_{T}$ in the $t\bar{t}$ CR. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also shown separately, normalized to $10^{3}$ times the best fit of the signal yield. The last bin includes the overflow. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the leading lepton $p_{T}$ in the $t\bar{t}Z$ CR. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also shown separately, normalized to $10^{3}$ times the best fit of the signal yield. The last bin includes the overflow. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{1}$ discriminant in the $t\bar{t}Z$ CR. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also shown separately, normalized to $10^{3}$ times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Bayesian upper exclusion limits on the ratio $g_{W'}/g_{W}$ for an SSM-like W$\prime$ boson are shown. The unity coupling ratio (blue dotted curve) corresponds to the SSM common benchmark. The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian lower exclusion limits on the NUGIM G(221) mixing angle $\cot(\theta_{E})$ are shown as a function of the W$\prime$ boson mass. The theoretically excluded region is shaded in grey. The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian upper exclusion limits at 95% CL on the product of the production cross section and branching fraction of a QBH in an associated $ au$ lepton and neutrino final state. Minimum threshold masses $m_{th}$ of up to 6.6 TeV are excluded at 95% CL. The observed limit (solid line) is obtained from the intersection with the LO QBH cross section (blue dotted curve). The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian upper limits at 95% CL on the cross section of the process $pp\rightarrow\tau\nu$ mediated via LQ exchange in the t-channel. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The predicted LQ cross section at LO in the three coupling benchmark scenarios is depicted in different colors for $g_{U}=1$. The uncertainty bandy correspond to the sum in quadrature of PDF and scale variations. The first benchmark scenario considers only couplings to left-handed SM fermions (i.e. $\beta_{\text{R}}^{ij}=0$) and is referred to as "best fit LH". The second benchmark, referred to as "best fit LH+RH", considers $|\beta_{\text{R}}^{\text{b}\tau}|=1$ and all other $\beta_{\text{R}}^{ij}=0$. In the third "democratic" benchmark, equal couplings only to LH fermions are assumed, i.e. $\beta_{\text{L}}^{ij}=1$ and $\beta_{\text{R}}^{ij}=0$ for all $i$ and $j$.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the LH scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the LH+RH scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the democratic scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Bayesian upper exclusion limits at 95% CL on each of the Wilson coefficients described by the EFT model based on 2016-2018 data. The three different coupling types represent a left-handed vector coupling ($\epsilon^{cb}_{L}$), tensor-like coupling ($\epsilon^{cb}_{T}$), and scalar-tensor-like coupling ($\epsilon^{cb}_{S_{L}}$). The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Summary of exclusion limits (expected and observed) calculated at 95% CL for full Run-2 CMS data.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "year_binnumber", following the binning from Figure 4.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "year_binnumber", following the binning used in Figure 4.

Predicted signal yields for the 2017 data-taking period, corresponding to 41.3 fb$^{-1}$, after the application of each search requirement (cumulative) for various signal hypothesis. The requirements listed are presented as total efficiencies w.r.t. the previous selection step.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

The distributions of $m_{\mathrm{b1},\mathrm{W-tagged}}$ in the 0L inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $m_{\mathrm{T}}^{\mathrm{b,E_{\mathrm{T}^{\mathrm{miss}}}}}$ in the 0L inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $N_{\mathrm{W-tagged}}$ in the 0L inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $m_{\mathrm{b1},\mathrm{\cancel{b1}}}$ in the hadronic top inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $m_{\mathrm{b1},\mathrm{\cancel{b1}}}$ in the leptonic top inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $m_{\mathrm{b1},\mathrm{\cancel{b1}}}$ in the leptonic top inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Cutflow for the reference point $(\it{m}_{\mathrm{H^{\pm}}}, \it{m}_{a}, tan\beta, sin\theta )=$ (500,100,1,0.7) , (800,150,20,0.7), (600,250,30,0.7), (1000,400,1,0.7) in 0L regions. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(\it{m}_{\mathrm{H^{\pm}}}, \it{m}_{a}, tan\beta, sin\theta )=$ (500,100,1,0.7) , (800,150,20,0.7), (600,250,30,0.7), (1000,400,1,0.7) in 1L leptonic top regions. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(\it{m}_{\mathrm{H^{\pm}}}, \it{m}_{a}, tan\beta, sin\theta )=$ (500,100,1,0.7) , (800,150,20,0.7), (600,250,30,0.7), (1000,400,1,0.7) in 1L hadronic top regions. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Distribution of $m(e^{\pm},\mu^{\pm})_{\mathrm{lead}}$ in the electron-muon signal region after the background-only fit.

Distribution of $m(\mu^{\pm},\mu^{\pm})_{\mathrm{lead}}$ in the muon-muon signal region after the background-only fit.

Distribution of $m(\ell^{\pm},\ell^{\prime \pm})_{\mathrm{lead}}$ in the three-lepton signal region after the background-only fit.

Event yield in the four-lepton signal region after the background-only fit.

Two-lepton signal region weighted cutflows for electron-electron channel. Data and a few representative signal samples are presented alongside the main background contributions. Yield stands for the observed yield in data and the expected yields based on $\sigma \cdot \mathcal{L}$ for the signal and background samples. Top quark and multiboson process backgrounds are merged under the 'other' category.

Two-lepton signal region weighted cutflows for electron-muon channel. Data and a few representative signal samples are presented alongside the main background contributions. Yield stands for the observed yield in data and the expected yields based on $\sigma \cdot \mathcal{L}$ for the signal and background samples. Top quark and multiboson process backgrounds are merged under the 'other' category.

Two-lepton signal region weighted cutflows for muon-muon channel. Data and a few representative signal samples are presented alongside the main background contributions. Yield stands for the observed yield in data and the expected yields based on $\sigma \cdot \mathcal{L}$ for the signal and background samples. Top quark and multiboson process backgrounds are merged under the 'other' category.

Three-lepton signal region weighted cutflow for $\ell \ell \ell$ channel. Data and a few representative signal samples are presented alongside the main background contributions. Yield stands for the observed yield in data and the expected yields based on $\sigma \cdot \mathcal{L}$ for the signal and background samples. Top quark and multiboson process backgrounds are merged under the 'other' category.

Four-lepton signal region weighted cutflow for $\ell \ell \ell \ell$ channel. Data and a few representative signal samples are presented alongside the main background contributions. Yield stands for the observed yield in data and the expected yields based on $\sigma \cdot \mathcal{L}$ for the signal and background samples. Top quark, Drell-Yan, and multiboson process backgrounds are merged under the 'other' category.

Probability of finding at least one TAR jet, where the p_T-leading TAR jet passes the m_Wcand and D₂^&beta;=1 requirements, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=2100 GeV, with the preselections applied.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=500 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=1000 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=1700 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=2100 GeV, with the preselections applied that do not pass the requirements of the merged category.

Observed exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected+1&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected-1&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected+2&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected-2&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=0.5 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=0.5 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=1 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=1 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=1.7 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=1.7 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=2.1 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=2.1 TeV, shown in dashed blue.

Data overlaid on SM background yields stacked in each SR and CR category after the fit to data ('Post-fit'). The yields in the SR are broken down into their contributions to the individual bins. The maximum-likelihood estimators are set to the conditional values of the CR-only fit, and propagated to SR and CRs.

Dominant sources of uncertainty for three dark Higgs scenarios after the fit to data. The uncertainties are quantified in terms of their contribution to the fitted signal uncertainty that is expressed relative to the theory prediction. Three representative dark Higgs signal scenarios with g_q=0.25, g_&chi;=1.0, sin&theta;=0.01 and m_&chi;=200 GeV are considered, which are indicated using the (m_Z', m_s) format in units of GeV in the table columns.

Cumulative efficiencies in the merged category for three representative dark Higgs signal scenarios with g_q=0.25, g_&chi;=1.0, sin&theta;=0.01, m_Z' = 1 TeV, and m_&chi;=200 GeV considering s&rarr;W(&ell;&nu;)W(qq) decays only.

Cumulative efficiencies in the resolved category for three representative dark Higgs signal scenarios with g_q=0.25, g_&chi;=1.0, sin&theta;=0.01, m_Z' = 1 TeV, and m_&chi;=200 GeV considering s&rarr;W(&ell;&nu;)W(qq) decays only.

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={300, 350, 400, 500, 750} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={1000, 1700} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={2100, 2500, 2900, 3300} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Parametrization of the $A_{X}$ factor, defined as the fraction of diphoton resonances satisfying the fiducial acceptance at the particle level, as function of $m_{X}$ obtained from the narrow width simulated signal samples produced in gluon fusion.

The correction factor, $C_{X}$, defined as the ratio between the number of reconstructed signal events passing the analysis cuts and the number of signal events at the particle level generated within the fiducial volume, and acceptance correction factor, $A_{X}$, defined as the fraction of diphoton resonances satisfying the fiducial acceptance at the particle level. Both are computed for NWA spin-0 models as a function of $m_{X}$.

Effect of event selections on a scalar MC signal sample generated for $m_{X}$ = 15 GeV and on the data. For the MC sample, the efficiencies are shown after applying event weights and a truth level filter that requires two photons with $p^{\gamma\gamma}_{T}>40$ GeV; for the data, the absolute yields are shown. The initial yields for data include a trigger preselection that is the OR of a list of single photon and diphoton triggers. The "2 $loose$ photons" step includes the kinematic acceptance cuts.

Parameterization of the Double Sided Crystal Ball function parameters describing the scalar mass resolution model as a function of $m_{X}$ [GeV].

Pre-fit comparison between data and background in the baseline SR for two of the variables used as input for the BSM pBDT: HT.

Data and post-fit background comparison for the distributions of the discriminant variables fitted in the control regions obtained with the background-only fit: CR HF e.

Data and post-fit background comparison for the distributions of the discriminant variables fitted in the control regions obtained with the background-only fit: CR HF mu.

Data and post-fit background comparison for the distributions of the discriminant variables fitted in the control regions obtained with the background-only fit: CR Conv.

Data and post-fit background comparison for the distributions of the discriminant variables fitted in the control regions obtained with the background-only fit: CR ttW.

Data and post-fit background comparison for the distributions of the discriminant variables fitted in the control regions obtained with the background-only fit: CR lowBDT.

Data and post-fit background comparison obtained with the background-only fit in the validation region defined for ttW+jets events: number of jets.

Data and post-fit background comparison obtained with the background-only fit in the validation region defined for ttW+jets events: SM BDT.

Data and post-fit background comparison obtained with the background-only fit to the BSM SR for the BSM pBDT distribution used for mH = 400 GeV.

Data and post-fit background comparison obtained with the background-only fit to the BSM SR for the BSM pBDT distribution used for mH = 1000 GeV.

Observed (black solid line) and expected (black dashed line) 95% CL upper limits on the ttH/A cross-section times branching fraction of H/A->tt as a function of mH/A.

Observed (red line) and expected (black dashed line) exclusion regions at 95% CL in the tan(beta) versus mass plane assuming that both, a heavy scalar H and pseudo-scalar A, contribute to the tttt final state and have the same mass mH=mA.

Observed (red line) and expected (black dashed line) exclusion regions at 95% CL in the tan(beta) versus mass plane assuming that only the scalar H contributes to the tttt final state.

Data and post-fit background comparison obtained with the background-only fit to the BSM SR for the BSM pBDT distribution used for mH = 500 GeV.

Data and post-fit background comparison obtained with the background-only fit to the BSM SR for the BSM pBDT distribution used for mH = 600 GeV.

Data and post-fit background comparison obtained with the background-only fit to the BSM SR for the BSM pBDT distribution used for mH = 700 GeV.

Data and post-fit background comparison obtained with the background-only fit to the BSM SR for the BSM pBDT distribution used for mH = 800 GeV.

Data and post-fit background comparison obtained with the background-only fit to the BSM SR for the BSM pBDT distribution used for mH = 900 GeV.

Observed (red line) and expected (black dashed line) exclusion regions at 95% CL in the tan(beta) versus mass plane assuming that only the pseudo-scalar A contributes to the tttt final state.

The product of acceptance and efficiency (A×ε), defined as the number of signal events satisfying the full set of selection criteria of the BSM SR, divided by the total number of generated signal events, for the tt̄H(→ tt̄ ) process in a type-II 2HDM for tanβ = 1 for the different heavy Higgs mass hypotheses.

Observed 95% CL limits on the production cross-section for scalar up-type LQs decaying into a top quark and a neutrino or a bottom quark and an electron.

Expected 95% CL limits on the production cross-section for scalar up-type LQs decaying into a top quark and a neutrino or a bottom quark and an electron.

Observed 95% CL limits on the production cross-section for scalar down-type LQs decaying into a bottom quark and a neutrino or a top quark and a muon.

Expected 95% CL limits on the production cross-section for scalar down-type LQs decaying into a bottom quark and a neutrino or a top quark and a muon.

Observed 95% CL limits on the production cross-section for scalar down-type LQs decaying into a bottom quark and a neutrino or a top quark and an electron.

Expected 95% CL limits on the production cross-section for scalar down-type LQs decaying into a bottom quark and a neutrino or a top quark and an electron.

Observed 95% CL limits on the production cross-section for vector up-type LQs in the Yang-Mills coupling scenario decaying into a top quark and a neutrino or a bottom quark and a muon.

Expected 95% CL limits on the production cross-section for vector up-type LQs in the Yang-Mills coupling scenario decaying into a top quark and a neutrino or a bottom quark and a muon.

Observed 95% CL limits on the production cross-section for vector up-type LQs in the Yang-Mills coupling scenario decaying into a top quark and a neutrino or a bottom quark and an electron.

Expected 95% CL limits on the production cross-section for vector up-type LQs in the Yang-Mills coupling scenario decaying into a top quark and a neutrino or a bottom quark and an electron.

Observed 95% CL limits on the production cross-section for vector up-type LQs in the minimal coupling scenario decaying into a top quark and a neutrino or a bottom quark and a muon.

Expected 95% CL limits on the production cross-section for vector up-type LQs in the minimal coupling scenario decaying into a top quark and a neutrino or a bottom quark and a muon.

Observed 95% CL limits on the production cross-section for vector up-type LQs in the minimal coupling scenario decaying into a top quark and a neutrino or a bottom quark and an electron.

Expected 95% CL limits on the production cross-section for vector up-type LQs in the minimal coupling scenario decaying into a top quark and a neutrino or a bottom quark and an electron.

Product of signal acceptance and efficiency in the training region for scalar up-type LQs decaying into top quarks and neutrinos or bottom quarks and muons depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(LQ_{mix}^u\rightarrow b\mu)$.

Product of signal acceptance and efficiency in the training region for scalar up-type LQs decaying into top quarks and neutrinos or bottom quarks and electrons depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(LQ_{mix}^u\rightarrow be)$.

Product of signal acceptance and efficiency in the training region for scalar down-type LQs decaying into top quarks and muons or bottom quarks and neutrinos depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(LQ_{mix}^d\rightarrow t\mu)$.

Product of signal acceptance and efficiency in the training region for scalar down-type LQs decaying into top quarks and electrons or bottom quarks and neutrinos depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(LQ_{mix}^d\rightarrow te)$.

Product of signal acceptance and efficiency in the training region for vector LQs in the Yang-Mills coupling scenario decaying into top quarks and neutrinos or bottom quarks and muons depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(vLQ_{mix}^{YM}\rightarrow b\mu)$.

Product of signal acceptance and efficiency in the training region for vector LQs in the Yang-Mills coupling scenario decaying into top quarks and neutrinos or bottom quarks and electrons depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(vLQ_{mix}^{YM}\rightarrow be)$.

Product of signal acceptance and efficiency in the training region for vector LQs in the minimal coupling scenario decaying into top quarks and neutrinos or bottom quarks and muons depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(vLQ_{mix}^{min}\rightarrow b\mu)$.

Product of signal acceptance and efficiency in the training region for vector LQs in the minimal coupling scenario decaying into top quarks and neutrinos or bottom quarks and electrons depending on the signal mass, $m_{LQ}$, and $\mathcal{B}(vLQ_{mix}^{min}\rightarrow be)$.

Cut-flow for the four models containing scalar LQs at a mass of $m_{LQ}=1300$ GeV and $\mathcal{B}(LQ_{mix}\rightarrow q_3\ell)=0.5$.

Cut-flow for the four models containing vector LQs at a mass of $m_{LQ}=1300$ GeV and $\mathcal{B}(LQ_{mix}\rightarrow q_3\ell)=0.5$.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[110,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[110,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[120,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[120,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[130,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[130,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[100,105)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[100,105)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[105,110)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[105,110)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[110,115)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[110,115)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[115,120)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[115,120)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[120,125)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[120,125)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[125,130)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[125,130)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[130,140)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[130,140)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[140,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-0J $m_{\mathrm{T2}}^{100} \in[140,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[100,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[100,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[110,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[110,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[120,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[120,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[130,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[130,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[100,105)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[100,105)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[105,110)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[105,110)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[110,115)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[110,115)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[115,120)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[115,120)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[120,125)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[120,125)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[125,130)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[125,130)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[130,140)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[130,140)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[140,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the slepton pair production model, in the SR-1J $m_{\mathrm{T2}}^{100} \in[140,\infty)$ region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

Cutflow table for the slepton signal sample with $m(\tilde{\ell},\tilde{\chi}_1^0) = (100,70)$ GeV, in the SR-0J $m_{\mathrm{T2}}^{100} \in [100,\infty)$ region. The yields include the process cross section and are weighted to the 139 fb$^{-1}$ luminosity. 246000 events were generated for the sample.

Cutflow table for the slepton signal sample with $m(\tilde{\ell},\tilde{\chi}_1^0) = (100,70)$ GeV, in the SR-1J $m_{\mathrm{T2}}^{100} \in [100,\infty)$ region. The yields include the process cross section and are weighted to the 139 fb$^{-1}$ luminosity. 246000 events were generated for the sample.

Observed and expected exclusion limits on SUSY simplified models, with observed upper limits on signal cross-section (fb) overlaid, for slepton-pair production in the $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ plane. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the (a) $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\ell})-\Delta m(\tilde{\ell},\tilde{\chi}_1^0)$ planes. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP for $\tilde{\mu}_{\textup{R}}$ and by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for direct selectron production in the (a) $m(\tilde{e})-m(\tilde{\chi}_1^0)$ and (c) $m(\tilde{e})-\Delta m(\tilde{e},\tilde{\chi}_1^0)$ planes, and for direct smuon production in the (b) $m(\tilde{\mu})-m(\tilde{\chi}_1^0)$ and (d) $m(\tilde{\mu})-\Delta m(\tilde{\mu},\tilde{\chi}_1^0)$ planes. In Figure (a) and (c) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{e}_{\textup{L,R}}$ and for $\tilde{e}_{\textup{L}}$ and $\tilde{e}_{\textup{R}}$. In Figure (b) and (d) the observed (solid thick lines) and expected (dashed lines) exclusion contours are indicated for combined $\tilde{\mu}_{\textup{L,R}}$ and for $\tilde{\mu}_{\textup{L}}$. No unique sensitivity to $\tilde{\mu}_{\textup{R}}$ is observed. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown in the shaded areas.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ plane. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The red contour shows the exclusion limits obtained using both the SR-0J and SR-1J region, as presented in Figure 6. The blue and green contours correspond to the result obtained considering only SR-0J and SR-1J region respectively. All limits are computed at 95% CL. The observed limits obtained by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ plane. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The red contour shows the exclusion limits obtained using both the SR-0J and SR-1J region, as presented in Figure 6. The blue and green contours correspond to the result obtained considering only SR-0J and SR-1J region respectively. All limits are computed at 95% CL. The observed limits obtained by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ plane. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The red contour shows the exclusion limits obtained using both the SR-0J and SR-1J region, as presented in Figure 6. The blue and green contours correspond to the result obtained considering only SR-0J and SR-1J region respectively. All limits are computed at 95% CL. The observed limits obtained by the ATLAS experiment in previous searches are also shown.

Observed and expected exclusion limits on SUSY simplified models for slepton-pair production in the $m(\tilde{\ell})-m(\tilde{\chi}_1^0)$ plane. Only $\tilde{e}$ and $\tilde{\mu}$ are considered. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The red contour shows the exclusion limits obtained using both the SR-0J and SR-1J region, as presented in Figure 6. The blue and green contours correspond to the result obtained considering only SR-0J and SR-1J region respectively. All limits are computed at 95% CL. The observed limits obtained by the ATLAS experiment in previous searches are also shown.

The upper panel shows the observed number of events in each of the binned SRs defined in Table 3, together with the expected SM backgrounds obtained after applying the efficiency correction method to compute the number of expected FSB events. `Others' include the non-dominant background sources, e.g. $t \bar{t}$+$V$, Higgs boson and Drell--Yan events. The uncertainty band includes systematic and statistical errors from all sources. The distributions of two signal points with mass splittings $\Delta m(\tilde{\ell},\tilde{\chi}_1^0) = m(\tilde{\ell})-m(\tilde{\chi}_1^0) = 30$ GeV and $\Delta m(\tilde{\ell},\tilde{\chi}_1^0) = m(\tilde{\ell})-m(\tilde{\chi}_1^0) = 50$ GeV are overlaid. The lower panel shows the significance as defined in Ref. [115].

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR$^{\text{-DF BDT-signal}\in(0.81,1]}_{\text{-SF BDT-signal}\in(0.77,1]}$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR$^{\text{-DF BDT-signal}\in(0.81,1]}_{\text{-SF BDT-signal}\in(0.77,1]}$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.81,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.81,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.82,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.82,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.83,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.83,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.84,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.84,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.85,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.85,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.81,0.8125]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.81,0.8125]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8125,0.815]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8125,0.815]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.815,0.8175]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.815,0.8175]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8175,0.82]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8175,0.82]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.82,0.8225]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.82,0.8225]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8225,0.825]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8225,0.825]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.825,0.8275]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.825,0.8275]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8275,0.83]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8275,0.83]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.83,0.8325]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.83,0.8325]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8325,0.835]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8325,0.835]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.835,0.8375]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.835,0.8375]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8375,0.84]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.8375,0.84]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.84,0.845]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.84,0.845]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.845,0.85]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.845,0.85]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.85,0.86]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.85,0.86]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.86,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-DF BDT-signal$\in(0.86,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.77,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.77,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.78,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.78,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.79,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.79,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.80,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.80,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.77,0.775]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.77,0.775]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.775,0.78]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.775,0.78]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.78,0.785]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.78,0.785]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.785,0.79]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.785,0.79]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.79,0.795]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.79,0.795]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.795,0.80]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.795,0.80]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.80,0.81]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.80,0.81]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.81,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

The figure shows the signal acceptance (a) and efficiency (b) plots for the $\tilde{\chi}_1^+\tilde{\chi}_1^-$ production with $W$-boson-mediated decay model, in the SR-SF BDT-signal$\in(0.81,1]$ inclusive region. Acceptance is calculated by applying the signal region requirements to particle-level objects, which do not suffer from identification inefficiencies or mismeasurements. The efficiency is calculated with fully reconstructed objects with the acceptance divided out. Large acceptance and efficiency differences in neighbouring points are due to statistical fluctuations.

Cutflow table for the chargino signal sample with $m\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0=(125,25)$ GeV, in the SR-SF BDT-signal$\in (0.77,1]$ and SR-DF BDT-signal$\in (0.81,1]$ regions. The yields include the process cross-section and are weighted to the 139 fb$^{-1}$ luminosity. 170000 events were generated for the sample.

Observed and expected exclusion limits on SUSY simplified models, with observed upper limits on signal cross-section (fb) overlaid, for chargino-pair production with $W$-boson-mediated decays in the $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ plane. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

Observed and expected exclusion limits on SUSY simplified models for chargino-pair production with $W$-boson-mediated decays in the (a) $m(\tilde{\chi}_1^{\pm})-m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{\chi}_1^{\pm})-\Delta m(\tilde{\chi}_1^{\pm},\tilde{\chi}_1^0)$ planes. The observed (solid thick line) and expected (thin dashed line) exclusion contours are indicated. The shaded band around the dashed line corresponds to the $\pm 1 \sigma$ variations in the expected limit, including all uncertainties except theoretical uncertainties in the signal cross-section. The dotted lines around the observed limit illustrate the change in the observed limit as the nominal signal cross-section is scaled up and down by the theoretical uncertainty. All limits are computed at 95% CL. The observed limits obtained at LEP and by the ATLAS experiment in previous searches are also shown. In case of the search performed on ATLAS Run 1 data at $\sqrt{s} = 8$ TeV no sensitivity was expected for the exclusion in the mass plane.

The upper panel shows the observed number of events in the SRs defined in Table 3, together with the expected SM backgrounds obtained after the background fit in the CRs. `Others' include the non-dominant background sources, e.g.$t \bar{t}$+$V$, Higgs boson and Drell--Yan events. The uncertainty band includes systematic and statistical errors from all sources. Distributions for three benchmark signal points are overlaid for comparison. The lower panel shows the significance as defined in Ref. [115].

Cumulative efficiency for simulated signal events to pass each stage of the event selection in the hadhad category. The efficiencies are calculated with respect to HH to bbtautau decays in which both tau-leptons decay hadronically. The ‘Pre-selection’ consists of basic requirements, including that at least two tau-had-vis pass loose kinematic requirements, at least one of the tau-had-vis candidate passes identification requirement, and that events do not contain an electron or muon. The ‘Object selections’ require exactly two tau-had-vis candidates, and at least two jets with pT > 25 GeV and abs(eta) < 2.5. The ‘Loose tau ID’ requires that both tau-had-vis candidates pass identification requirements. The ‘DTT offline jet cuts’ are cuts placed on the pT of the reconstructed jet or tau-had-vis that are geometrically matched to the HLT objects, to ensure the efficiencies of the HLT objects reach the plateau region.

Cumulative efficiency for simulated signal events to pass each stage of the event selection in the lephad SLT category. The efficiencies are calculated with respect to HH to bbtautau decays in which one tau-lepton decays hadronically and one decays leptonically. The ‘Pre-selection’ consists of basic requirements, including that at least one tau-had-vis candidate and one lepton pass loose kinematic requirements. The ‘Object selections’ require exactly one tau-had-vis candidate, and at least two jets with pT > 25 GeV and abs(eta) < 2.5. The ‘Trigger specific offline pT cuts’ are cuts placed on the pT of the reconstructed jet or tau-had-vis that are geometrically matched to the HLT objects, to ensure the efficiencies of the HLT objects reach the plateau region.

Cumulative efficiency for simulated signal events to pass each stage of the event selection in the lephad LTT category. The efficiencies are calculated with respect to HH to bbtautau decays in which one tau-lepton decays hadronically and one decays leptonically. The ‘Pre-selection’ consists of basic requirements, including that at least one tau-had-vis candidate and one lepton pass loose kinematic requirements. The ‘Object selections’ require exactly one tau-had-vis candidate, and at least two jets with pT > 25 GeV and abs(eta) < 2.5. The ‘Trigger specific offline pT cuts’ are cuts placed on the pT of the reconstructed jet or tau-had-vis that are geometrically matched to the HLT objects, to ensure the efficiencies of the HLT objects reach the plateau region.

Post-fit expected number of signal and background events and observed number of data events after applying the selection criteria and requiring exactly 2 b-tagged jets and assuming a background-only hypothesis.

Observed and expected limits at 95% CL on the cross-section of HH production, for the non-resonant ggF+VBF HH search, and the resonant HH search for four values of the resonance mass mX.

Acceptance times efficiency for the full analysis selections as a function of the resonance mass mX in the hadhad, lephad SLT and lephad LTT trigger categories, and the combined lephad. The acceptance times efficiency is evaluated for HH to bbtautau decays, with respect to the targeted tau-lepton decay modes (lephad or hadhad).

Post-fit distribution of mHH in the hadhad channel.

Post-fit distribution of mHH in the lephad SLT channel.

Post-fit distribution of mHH in the lephad LTT channel.

Post-fit distribution of mtautau MMC in the hadhad channel.

Post-fit distribution of mtautau MMC in the hadhad channel.

Post-fit distribution of mtautau MMC in the hadhad channel.

Post-fit distribution of of di-b-jet mass in the hadhad channel.

Post-fit distribution of di-b-jet mass in the LTT lephad channel.

Post-fit distribution of di-b-jet mass in the LTT lephad channel.

BDT for SM HH in the hadhad channel

NN for SM HH in the lephad SLT channel

NN for SM HH in the lephad LTT channel

PNN for mX = 500 GeV resonant HH in the hadhad channel

PNN for mX = 500 GeV resonant HH in the lephad SLT channel

PNN for mX = 500 GeV resonant HH in the lephad LTT channel

PNN for mX = 1000 GeV resonant HH in the hadhad channel

PNN for mX = 1000 GeV resonant HH in the lephad SLT channel

PNN for mX = 1000 GeV resonant HH in the lephad LTT channel

Event yields as a function of log10(S/B) for data, background and non-resonant HH signal. Final discriminant bins from the hadhad, lephad SLT and lephad LTT categories are combined into bins of log10(S/B). The B is the fitted background yield assuming background-only hypothesis, and the signal S is scaled to the SM expected cross-section.

Observed and expected limits at 95% CL on the cross-section of the resonant HH production as a function of the scalar resonance mass mX.

Post-fit distribution of delta-R between the taus in the hadhad channel.

Post-fit distribution of delta-R between the b-tagged jets in the hadhad channel.

Post-fit distribution of delta-R between the taus in the lephad SLT channel.

Post-fit distribution of delta-R between the b-tagged jets in the lephad SLT channel.

Post-fit distribution of delta-pT between the tau and lepton in the lephad SLT channel.

Post-fit distribution of pT of the subleading b-tagged jet in the lephad SLT channel.

Post-fit distribution of MTW in the lephad SLT channel.

Post-fit distribution of missing transverse momentum in the lephad SLT channel.

Post-fit distribution of missing transverse momentum centrality in the lephad SLT channel.

Post-fit distribution of delta-phi between the Higgs boson candidates in the lephad SLT channel.

Post-fit distribution of delta-pT between the tau and lepton in the lephad LTT channel.

Post-fit distribution of delta-R between the taus in the lephad LTT channel.

Post-fit distribution of delta-phi between the lepton and the missing transverse momentum in the lephad LTT channel.

Post-fit distribution of delta-phi between the Higgs boson candidates in the lephad LTT channel.

Post-fit distribution of the total transverse momentum s in the lephad LTT channel.

PNN for mX = 300 GeV resonant HH in the hadhad channel

PNN for mX = 300 GeV resonant HH in the lephad SLT channel

PNN for mX = 300 GeV resonant HH in the lephad LTT channel

PNN for mX = 1600 GeV resonant HH in the hadhad channel

PNN for mX = 1600 GeV resonant HH in the lephad SLT channel

PNN for mX = 1600 GeV resonant HH in the lephad LTT channel

Local p-value of the background-only hypothesis as a function of the resonance mass.

Correlation matrix for the individual sources of systematic uncertainty considered in the analysis, for the combined opposite sign (OS) and same sign (SS) binned-template profile likelihood fit to data. The OS or SS refers to the charge signs of the primary lepton and the soft muon. The gamma parameters are NPs used to describe the effect of the limited statistics of the samples.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Leading tau Pt distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{h}\tau_{lep}\tau_{had}$-2j region. Other MC includes single top, V+jets, and other small backgrounds. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Leading tau Pt distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{h}\tau_{lep}\tau_{had}$-3j region. Other MC includes single top, V+jets, and other small backgrounds. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Leading tau Pt distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{\ell}\tau_{had}\tau_{had}$SS region. Other MC includes single top, V+jets, and other small backgrounds. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Leading tau Pt distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{h}\tau_{had}\tau_{had}$-2j region. Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Leading tau Pt distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{h}\tau_{had}\tau_{had}$-3j region. Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Leading tau Pt distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{h}\tau_{had}\tau_{had}$-3jSS region. Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Di-tau mass distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{\ell}\tau_{had}\tau_{had}$ region. Other MC includes single top, V+jets, and other small backgrounds. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Di-tau mass distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{\ell}\tau_{had}$-1j region. Other MC includes single top, V+jets, and other small backgrounds. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Di-tau mass distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{\ell}\tau_{had}$-2j region. Other MC includes single top, V+jets, and other small backgrounds. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Di-tau mass distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{h}\tau_{lep}\tau_{had}$-2j region. Other MC includes single top, V+jets, and other small backgrounds. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Di-tau mass distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{h}\tau_{lep}\tau_{had}$-3j region. Other MC includes single top, V+jets, and other small backgrounds. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Di-tau mass distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{\ell}\tau_{had}\tau_{had}$SS region. Other MC includes single top, V+jets, and other small backgrounds. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Di-tau mass distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{h}\tau_{had}\tau_{had}$-2j region. Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Di-tau mass distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{h}\tau_{had}\tau_{had}$-3j region. Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

Di-tau mass distributions obtained before the fit to data (Pre-Fit) showing the expected background and tuH signals after applying fake factors in the $t_{h}\tau_{had}\tau_{had}$-3jSS region. Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. The tuH signals with nominal branching ratio of 0.1% are scaled using normalization factors of 2 to 50. Statistical and systematic uncertainties are included in the "Total background".

BDT output distributions obtained from a signal+background fit to the data for the tuH search in the $t_{\ell}\tau_{had}\tau_{had}$ region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(uH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tuH search in the $t_{\ell}\tau_{had}$-1j region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(uH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tuH search in the $t_{\ell}\tau_{had}$-2j region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(uH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tuH search in the $t_{h}\tau_{lep}\tau_{had}$-2j region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(uH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tuH search in the $t_{h}\tau_{lep}\tau_{had}$-3j region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(uH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tuH search in the $t_{\ell}\tau_{had}\tau_{had}$SS region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(uH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tuH search in the $t_{h}\tau_{had}\tau_{had}$-2j region, Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(uH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tuH search in the $t_{h}\tau_{had}\tau_{had}$-3j region, Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(uH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tuH search in the $t_{h}\tau_{had}\tau_{had}$-3jSS region, Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(uH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tcH search in the $t_{\ell}\tau_{had}\tau_{had}$ region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(cH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tcH search in the $t_{\ell}\tau_{had}$-1j region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(cH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tcH search in the $t_{\ell}\tau_{had}$-2j region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(cH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tcH search in the $t_{h}\tau_{lep}\tau_{had}$-2j region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(cH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tcH search in the $t_{h}\tau_{lep}\tau_{had}$-3j region,Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(cH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tcH search in the $t_{\ell}\tau_{had}\tau_{had}$SS region, Other MC includes single top, V+jets, and other small backgrounds. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(cH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tcH search in the $t_{h}\tau_{had}\tau_{had}$-2j region, Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(cH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tcH search in the $t_{h}\tau_{had}\tau_{had}$-3j region, Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(cH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$ of 0.1%.

BDT output distributions obtained from a signal+background fit to the data for the tcH search in the $t_{h}\tau_{had}\tau_{had}$-3jSS region, Rare includes single top, V+jets, and other small backgrounds. $\tau_{sub}$ real includes the contribution of fakes for which the sub-leading tau is real. Statistical and systematic uncertainties are included in the "Total background". The signal shapes of tt(cH), tH, and their sum are also shown using a normalisation of 2 x $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$ of 0.1%.

95% CL upper limits on $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$ for the individual searches as well as their combination, assuming $\mathcal{B}(\mathrm{t}\to\mathrm{uH}) = 0$. The observed limits are compared with the expected (median) limits under the background-only hypothesis. The surrounding shaded bands correspond to the 68% and 95% CL intervals around the expected limits, denoted by $\pm 1\sigma$ and $\pm 2\sigma$, respectively.

95% CL upper limits on $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ for the individual searches as well as their combination, assuming $\mathcal{B}(\mathrm{t}\to\mathrm{cH}) = 0$. The observed limits are compared with the expected (median) limits under the background-only hypothesis. The surrounding shaded bands correspond to the 68% and 95% CL intervals around the expected limits, denoted by $\pm 1\sigma$ and $\pm 2\sigma$, respectively.

Observed upper limits at 95% CL on the branching fractions in the plane of $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ and $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$.

Expected upper limits at 95% CL on the branching fractions in the plane of $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ and $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$.

Expected $+2\sigma$ upper limits at 95% CL on the branching fractions in the plane of $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ and $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$.

Expected $+1\sigma$ upper limits at 95% CL on the branching fractions in the plane of $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ and $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$.

Expected $-1\sigma$ upper limits at 95% CL on the branching fractions in the plane of $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ and $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$.

Expected $-2\sigma$ upper limits at 95% CL on the branching fractions in the plane of $\mathcal{B}(\mathrm{t}\to\mathrm{uH})$ and $\mathcal{B}(\mathrm{t}\to\mathrm{cH})$.

Observed upper limits at 95% CL on the anomalous couplings in the plane of $C_{\mathrm{c\phi}}$ and $C_{\mathrm{u\phi}}$.

Expected upper limits at 95% CL on the anomalous couplings in the plane of $C_{\mathrm{c\phi}}$ and $C_{\mathrm{u\phi}}$.

Expected $+2\sigma$ upper limits at 95% CL on the anomalous couplings in the plane of $C_{\mathrm{c\phi}}$ and $C_{\mathrm{u\phi}}$.

Expected $+1\sigma$ upper limits at 95% CL on the anomalous couplings in the plane of $C_{\mathrm{c\phi}}$ and $C_{\mathrm{u\phi}}$.

Expected $-1\sigma$ upper limits at 95% CL on the anomalous couplings in the plane of $C_{\mathrm{c\phi}}$ and $C_{\mathrm{u\phi}}$.

Expected $-2\sigma$ upper limits at 95% CL on the anomalous couplings in the plane of $C_{\mathrm{c\phi}}$ and $C_{\mathrm{u\phi}}$.

Predicted and observed yields in each of the analysis regions considered in leptonic channel.

Predicted and observed yields in each of the analysis regions considered in hadronic channel.

Absolute uncertainties on $\mathcal{B}(\mathrm{t}\to\mathrm{qH})$ obtained from the combined fit to the data.

Summary of 95% CL upper limits on $\mathcal{B}(\mathrm{t}\to\mathrm{qH})$, significance and best-fit branching ratio in the signal regions. The values in the tables are in the form of observed(expected).

Comparisons of the data and the expected background distributions of the WZ invariant mass in the ANN-based VBF signal region. The background predictions are obtained through a background-only simultaneous fit to the VBF signal region and the WZ-QCD and ZZ VBF control regions. The yields are normalized to the bin width

Comparisons of the observed data and the expected background distributions of the WZ invariant mass using the cut-based VBF selection. The background predictions are obtained through a background-only simultaneous fit to the VBF cut-based signal region and the WZ-QCD and ZZ VBF control regions. The yields are normalized to the bin width.

Comparisons of the observed data and the expected background distributions of the WZ invariant mass using the cut-based VBF selection. The background predictions are obtained through a background-only simultaneous fit to the VBF cut-based signal region and the WZ-QCD and ZZ VBF control regions. The yields are normalized to the bin width.

Drell-Yan signal region selection cutflow for a simulated W' in the HVT model A with m_W' = 1 TeV. The unweighted number of events is shown.

Drell-Yan signal region selection cutflow for a simulated W' in the HVT model A with m_W' = 1 TeV. The unweighted number of events is shown.

VBF signal region selection cutflow for a simulated W' in the HVT model C with m_W' = 500 GeV. The unweighted number of events is shown.

VBF signal region selection cutflow for a simulated W' in the HVT model C with m_W' = 500 GeV. The unweighted number of events is shown.

VBF signal region selection cutflow for a simulated H5+ in the GM model with m_H5+ = 450 GeV. The unweighted number of events is shown.

VBF signal region selection cutflow for a simulated H5+ in the GM model with m_H5+ = 450 GeV. The unweighted number of events is shown.

The acceptancetimes efficiencyof the HVT W' in the Drell-Yan signal region for different mass points and for the individual channels and the sum of all channels. The uncertainty includes both statistical and experimental systematic components.

The acceptancetimes efficiencyof the HVT W' in the Drell-Yan signal region for different mass points and for the individual channels and the sum of all channels. The uncertainty includes both statistical and experimental systematic components.

The acceptancetimes efficiencyof VBF H5+ selection after the ANN-based VBF selection at different mass points for the individual channels and the sum of all channels. The uncertainty includes both statistical and experimental systematic components.

The acceptancetimes efficiencyof VBF H5+ selection after the ANN-based VBF selection at different mass points for the individual channels and the sum of all channels. The uncertainty includes both statistical and experimental systematic components.

The acceptancetimes efficiencyof VBF HVT W' selection after the ANN-based VBF selection at different mass points for the individual channels and the sum of all channels. The uncertainty includes both statistical and experimental systematic components.

The acceptancetimes efficiencyof VBF HVT W' selection after the ANN-based VBF selection at different mass points for the individual channels and the sum of all channels. The uncertainty includes both statistical and experimental systematic components.

The acceptancetimes efficiencyof VBF H5+ selection after the cut-based VBF selection at different mass points for the individual channels and the sum of all channels. The uncertainty includes both statistical and experimental systematic components.

The acceptancetimes efficiencyof VBF H5+ selection after the cut-based VBF selection at different mass points for the individual channels and the sum of all channels. The uncertainty includes both statistical and experimental systematic components.

The acceptancetimes efficiencyof VBF HVT W' selection after the cut-based VBF selection at different mass points for the individual channels and the sum of all channels. The uncertainty includes both statistical and experimental systematic components.

The acceptancetimes efficiencyof VBF HVT W' selection after the cut-based VBF selection at different mass points for the individual channels and the sum of all channels. The uncertainty includes both statistical and experimental systematic components.

Observed and expected 95% CL exclusion upper limits on sigma * B(W' -> WZ) for the Drell-Yan production of a W' boson in the HVT model as a function of its mass. The LO theory predictions for HVT Model A with g_V=1 and Model B with g_V=$ are also shown.

Observed and expected 95% CL exclusion upper limits on sigma * B(W' -> WZ) for the Drell-Yan production of a W' boson in the HVT model as a function of its mass. The LO theory predictions for HVT Model A with g_V=1 and Model B with g_V=3 are also shown.

Using the ANN VBF selection, observed and expected 95% CL upper limits on sigma * B(W' -> WZ) for the VBF production of a W' boson in the HVT with parameter c_F=0, as a function of its mass. The LO theory predictions for HVT VBF model with different values of the coupling parameters g_V and c_H are also shown.

Using the ANN VBF selection, observed and expected 95% CL upper limits on sigma * B(W' -> WZ) for the VBF production of a W' boson in the HVT with parameter c_F=0, as a function of its mass. The LO theory predictions for HVT VBF model with different values of the coupling parameters g_V and c_H are also shown.

Using the ANN VBF selection, observed and expected 95% CL upper limits on sigma * B(W' -> WZ) of the GM model as a function of m_H_5.

Using the ANN VBF selection, observed and expected 95% CL upper limits on sigma * B(W' -> WZ) of the GM model as a function of m_H_5.

Using the ANN VBF selection, observed and expected 95% CL upper limits on sin(thetaH) of the GM model as a function of m_H_5.

Using the ANN VBF selection, observed and expected 95% CL upper limits on sin(thetaH) of the GM model as a function of m_H_5.

Using the cut-based VBF selection, observed and expected 95% CL upper limits on sigma * B(W' -> WZ) for the VBF production of a W' boson in the HVT with parameter c_F=0, as a function of its mass. The LO theory predictions for HVT VBF model with different values of the coupling parameters g_V and c_H are also shown.

Using the cut-based VBF selection, observed and expected 95% CL upper limits on sigma * B(W' -> WZ) for the VBF production of a W' boson in the HVT with parameter c_F=0, as a function of its mass. The LO theory predictions for HVT VBF model with different values of the coupling parameters g_V and c_H are also shown.

Using the cut-based VBF selection, observed and expected 95% CL upper limits on sigma * B(W' -> WZ) of the GM model as a function of m_H_5.

Using the cut-based VBF selection, observed and expected 95% CL upper limits on sigma * B(W' -> WZ) of the GM model as a function of m_H_5.

Using the cut-based VBF selection, observed and expected 95% CL upper limits on sin(thetaH) of the GM model as a function of m_H_5.

Using the cut-based VBF selection, observed and expected 95% CL upper limits on sin(thetaH) of the GM model as a function of m_H_5.

Correlation matrix for the measurement of STXS parameters in each of the 28 regions considered.

Fitted values for kappa_g and kappa_y.

Correlation matrix for the measurement of kappa_g and kappa_y.

Summary of the 68% CL confidence intervals for individual measurements of SMEFT parameters observed in data. In each case, SMEFT parameters other than the one measured are fixed to 0.

Results of the EV parameter measurement in data, in the linear and linear+quadratic parameterizations of the SMEFT. All the EVs parameters are free to vary in the fits. The ranges correspond to 68% CL confidence intervals.

Observed linear correlation coefficients of the EVs parameters in the linear parameterization.

Observed linear correlation coefficients of the EVs parameters in the linear+quadratic parameterization.

Cross-sections times H->yy branching ratio for ggF +bbH, VBF, VH, ttH, and tH production. The values are obtained from a simultaneous fit to all categories. The theory uncertainties in the predictions include uncertainties due to missing higher-order terms in the perturbative QCD calculations and choices of parton distribution functions and value of alpha_s, as well as the H->yy branching ratio uncertainty.

Best-fit values and uncertainties for STXS parameters in each of the 28 regions considered. The values for the gg->H process also include the contributions from bbH production.

Correlation matrix for the measurement of STXS parameters in each of the 33 regions considered.

Acceptance * reconstruction efficiency for the P P --> Zprime --> Zh --> llbb/cc signals in the 2-lepton channel.

Acceptance * reconstruction efficiency for the P P --> bbA --> Zh --> vvbb signals in the 0-lepton channel.

Acceptance * reconstruction efficiency for the P P --> bbA --> Zh --> vvbb signals in the 0-lepton channel.

Acceptance * reconstruction efficiency for the P P --> bbA --> Zh --> llbb signals in the 2-lepton channel.

Acceptance * reconstruction efficiency for the P P --> bbA --> Zh --> llbb signals in the 2-lepton channel.

Acceptance * reconstruction efficiency for the P P --> Wprime --> Zh --> lvbb/cc signals in the 0-lepton channel.

Acceptance * reconstruction efficiency for the P P --> Wprime --> Zh --> lvbb/cc signals in the 0-lepton channel.

Acceptance * reconstruction efficiency for the P P --> Wprime --> Zh --> lvbb/cc signals in the 1-lepton channel.

Acceptance * reconstruction efficiency for the P P --> Wprime --> Zh --> lvbb/cc signals in the 1-lepton channel.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the resolved 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the resolved 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the resolved 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the resolved 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the merged 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the merged 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the merged 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the merged 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the resolved 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the resolved 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the resolved 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the resolved 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the merged 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the merged 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the merged 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the merged 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 3+ b-tag signal region. The background prediction is shown after a background-only maximum-likelihood bbA fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 3+ b-tag signal region. The background prediction is shown after a background-only maximum-likelihood bbA fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the resolved 3+ b-tag signal region. The background prediction is shown after a background-only maximum-likelihood bbA fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the resolved 3+ b-tag signal region. The background prediction is shown after a background-only maximum-likelihood bbA fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 2 b-tag signal region with additional b-tagged track jets not associated with the large-R jet. The background prediction is shown after a background-only maximum-likelihood bbA fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 2 b-tag signal region with additional b-tagged track jets not associated with the large-R jet. The background prediction is shown after a background-only maximum-likelihood bbA fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the merged 1+2 b-tag signal region with additional b-tagged track jets not associated with the large-R jet. The background prediction is shown after a background-only maximum-likelihood bbA fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the merged 1+2 b-tag signal region with additional b-tagged track jets not associated with the large-R jet. The background prediction is shown after a background-only maximum-likelihood bbA fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the resolved 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the resolved 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the resolved 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the resolved 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the merged 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the merged 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the merged 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the merged 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the resolved top control region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the resolved top control region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Upper limits on Zprime to Z h production cross section times branching fraction in pb.

Upper limits on Zprime to Z h production cross section times branching fraction in pb.

Upper limits on Wprime to W h production cross section times branching fraction in pb.

Upper limits on Wprime to W h production cross section times branching fraction in pb.

Upper limits on ggA to Z h production cross section times branching fraction in pb.

Upper limits on ggA to Z h production cross section times branching fraction in pb.

Upper limits on bbA to Z h production cross section times branching fraction in pb.

Upper limits on bbA to Z h production cross section times branching fraction in pb.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 220 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 220 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 260 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 260 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 300 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 300 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 340 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 340 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 380 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 380 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 400 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 400 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 420 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 420 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 440 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 440 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 460 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 460 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 500 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 500 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 600 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 600 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 700 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 700 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 800 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 800 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 900 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 900 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 1000 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 1000 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 1200 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 1200 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 1400 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 1400 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 1600 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 1600 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 2000 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 2000 GeV.

Acceptance * reconstruction efficiency for the P P --> A --> Zh --> vvbb signal in the 0-lepton channel.

Acceptance * reconstruction efficiency for the P P --> A --> Zh --> vvbb signal in the 0-lepton channel.

Acceptance * reconstruction efficiency for the P P --> A --> Zh --> llbb signal in the 2-lepton channel.

Acceptance * reconstruction efficiency for the P P --> A --> Zh --> llbb signal in the 2-lepton channel.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Distributions of expected upper limits at 95% confidence level on the cross section of P P --> A --> Zh as a function of bbA fraction an signal mass.

Distributions of expected upper limits at 95% confidence level on the cross section of P P --> A --> Zh as a function of bbA fraction an signal mass.

Distributions of observed upper limits at 95% confidence level on the cross section of P P --> A --> Zh as a function of bbA fraction an signal mass.

Distributions of observed upper limits at 95% confidence level on the cross section of P P --> A --> Zh as a function of bbA fraction an signal mass.

Signal region detector-level distribution for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$.

Signal region detector-level distribution for the observable $m_{e\mu}$.

Signal region detector-level distribution for the observable $p_{\mathrm{T}}^{e\mu}$.

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

The reconstruction efficiency for caloDPJs produced by the decay of &gamma;_d into e⁺e^- or qq&#772;. The reconstruction efficiency is shown for &gamma;_d with 0&lt;|&eta;|&lt;1.1 as a function of their transverse momentum in events where the &gamma;_d L_xy is between 2&nbsp;m and 4&nbsp;m.