Observed values of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the single-Higgs and double-Higgs analyses combination, with all other coupling modifiers fixed to unity.

Observed values of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the double-Higgs analyses, with all other coupling modifiers fixed to unity.

Observed values of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the single-Higgs analyses, with all other coupling modifiers fixed to unity.

Observed values of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the single-Higgs and double-Higgs combination for the generic model (free floating $\kappa_t$, $\kappa_b$, $\kappa_V$ and $\kappa_\tau$). The observed best-fit value of $\kappa_\lambda$ for the generic model is shifted slightly relative to the other models because of its correlation with the best-fit values of the $\kappa_b$, $\kappa_t$ and $\kappa_\tau$ parameters, which are slightly below, but compatible with unity.

Expected values of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the single-Higgs and double-Higgs analyses combination derived from the combined single-Higgs and double-Higgs analyses, with all other coupling modifiers fixed to unity.

Expected values of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the double-Higgs analyses.

Expected values of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the single-Higgs analyses, with all other coupling modifiers fixed to unity.

Expected values of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the single-Higgs and double-Higgs analyses for the generic model (free floating $\kappa_t$, $\kappa_b$, $\kappa_V$ and $\kappa_\tau$).

Observed constraints in the $\kappa_\lambda$–$\kappa_t$ plane from single-Higgs and double-Higgs combination. The solid lines show the 68% CL contours. The observed constraint for the single- and double-Higgs combination for $\kappa_t$ values below unity is slightly less stringent than that for the single-Higgs fit alone due to the slightly higher best-fit value for this coupling modifier.

Observed constraints in the $\kappa_\lambda$–$\kappa_t$ plane from single-Higgs and double-Higgs combination. The dashed lines show the 95% CL contours. The observed constraint for the single- and double-Higgs combination for $\kappa_t$ values below unity is slightly less stringent than that for the single-Higgs fit alone due to the slightly higher best-fit value for this coupling modifier.

Observed constraints in the $\kappa_\lambda$–$\kappa_t$ plane from double-Higgs analysis. The solid lines show the 68% CL contours. The double-Higgs contours are shown for values of $\kappa_t$ smaller than 1.2.

Observed constraints in the $\kappa_\lambda$–$\kappa_t$ plane from double-Higgs analysis. The dashed lines show the 95% CL contours. The double-Higgs contours are shown for values of $\kappa_t$ smaller than 1.2.

Observed constraints in the $\kappa_\lambda$–$\kappa_t$ plane from single-Higgs analysis. The solid lines show the 68% CL contours.

Observed constraints in the $\kappa_\lambda$–$\kappa_t$ plane from single-Higgs analysis. The dashed lines show the 95% CL contours.

Expected constraints in the $\kappa_\lambda$–$\kappa_t$ plane from single-Higgs and double-Higgs combination. The solid lines show the 68% CL contours. The double-Higgs contours are shown for values of $\kappa_t$ smaller than 1.2.

Expected constraints in the $\kappa_\lambda$–$\kappa_t$ plane from single-Higgs and double-Higgs combination. The dashed lines show the 95% CL contours. The double-Higgs contours are shown for values of $\kappa_t$ smaller than 1.2.

Expected constraints in the $\kappa_\lambda$–$\kappa_t$ plane from double-Higgs analyses. The solid lines show the 68% CL contours. The double-Higgs contours are shown for values of $\kappa_t$ smaller than 1.2.

Expected constraints in the $\kappa_\lambda$–$\kappa_t$ plane from double-Higgs analyses. The dashed lines show the 95% CL contours. The double-Higgs contours are shown for values of $\kappa_t$ smaller than 1.2.

Expected constraints in the $\kappa_\lambda$–$\kappa_t$ plane from single-Higgs analyses. The solid lines show the 68% CL contours.

Expected constraints in the $\kappa_\lambda$–$\kappa_t$ plane from single-Higgs analyses. The dashed lines show the 95% CL contours.

Observed and expected 95% CL upper limits on the sum of the ggF HH and VBF HH production cross-section from the bbbb, bb$\tau\tau$ and bb$\gamma\gamma$ decay channels, and their statistical combination. The value $m_H$=125.09 GeV is assumed when deriving the predicted SM cross section. The expected limit and the corresponding error bands are derived assuming the absence of the HH process with all nuisance parameters profiled to the observed data. The SM prediction together with its theoretical uncertainty is also shown (red vertical band).

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the HH to bbbb analysis. All other coupling modifiers are fixed to their SM value.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the HH to bb$\tau\tau$ analysis. All other coupling modifiers are fixed to their SM value.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the HH to bb$\gamma\gamma$ analysis. All other coupling modifiers are fixed to their SM value.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for the double-Higgs combination. All other coupling modifiers are fixed to their SM value.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for HH to bbbb analysis. All other coupling modifiers are fixed to their SM value.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for HH to bb$\tau\tau$ analysis. All other coupling modifiers are fixed to their SM value.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for HH to bb$\gamma\gamma$ analysis. All other coupling modifiers are fixed to their SM value.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_\lambda$ parameter for double-Higgs combination. All other coupling modifiers are fixed to their SM value.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for the HH to bbbb analysis. All other coupling modifiers are fixed to their SM value.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for the HH to bb$\tau\tau$ analysis. All other coupling modifiers are fixed to their SM value.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for the HH to bb$\gamma\gamma$ analysis. All other coupling modifiers are fixed to their SM value.

Observed value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for the double-Higgs combination. All other coupling modifiers are fixed to their SM value.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for the HH to bbbb analysis. All other coupling modifiers are fixed to their SM value.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for the HH to bb$\tau\tau$ analysis. All other coupling modifiers are fixed to their SM value.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for the HH to bb$\gamma\gamma$ analysis. All other coupling modifiers are fixed to their SM value.

Expected value of the test statistic (-2ln$\Lambda$), as a function of the $\kappa_{2V}$ parameter for the double-Higgs combination. All other coupling modifiers are fixed to their SM value.

Observed constraints in the $\kappa_{2V}$–$\kappa_{V}$ plane from double-Higgs combination. The solid lines show the 68% (95%) CL contours.

Observed constraints in the $\kappa_{2V}$–$\kappa_{V}$ plane from double-Higgs combination. The dashed lines show the 68% (95%) CL contours.

Expected constraints in the $\kappa_{2V}$-$\kappa_{V}$ plane from double-Higgs combination. The solid lines show the 68% CL contours.

Expected constraints in the $\kappa_{2V}$-$\kappa_{V}$ plane from double-Higgs combination. The dashed lines show the 95% CL contours.

Acceptance times efficiency as a function of resonant mass for each event selection step in the search for a neutral heavy scalar resonance produced in association with a W boson decaying to a charged lepton and a neutrino.

Acceptance times efficiency as a function of resonant mass for each event selection step in the search for a neutral heavy scalar resonance produced in association with a Z boson decaying to charged leptons.

Acceptance times efficiency as a function of resonant mass for each event selection step in the search for a neutral heavy scalar resonance produced in association with a Z boson decaying to charged leptons.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy narrow-width pseudoscalar resonance that decays to a 260 GeV scalar resonance and a Z boson, which decays to neutrinos.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy narrow-width pseudoscalar resonance that decays to a 260 GeV scalar resonance and a Z boson, which decays to neutrinos.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy narrow-width pseudoscalar resonance that decays to a 300 GeV scalar resonance and a Z boson, which decays to neutrinos.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy narrow-width pseudoscalar resonance that decays to a 300 GeV scalar resonance and a Z boson, which decays to neutrinos.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy narrow-width pseudoscalar resonance that decays to a 400 GeV scalar resonance and a Z boson, which decays to neutrinos.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy narrow-width pseudoscalar resonance that decays to a 400 GeV scalar resonance and a Z boson, which decays to neutrinos.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy large-width pseudoscalar resonance that decays to a 260 GeV scalar resonance and a Z boson, which decays to neutrinos.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy large-width pseudoscalar resonance that decays to a 260 GeV scalar resonance and a Z boson, which decays to neutrinos.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy large-width pseudoscalar resonance that decays to a 300 GeV scalar resonance and a Z boson, which decays to neutrinos.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy large-width pseudoscalar resonance that decays to a 300 GeV scalar resonance and a Z boson, which decays to neutrinos.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy large-width pseudoscalar resonance that decays to a 400 GeV scalar resonance and a Z boson, which decays to neutrinos.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy large-width pseudoscalar resonance that decays to a 400 GeV scalar resonance and a Z boson, which decays to neutrinos.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy narrow-width pseudoscalar resonance that decays to a 260 GeV scalar resonance and a Z boson, which decays to charged leptons.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy narrow-width pseudoscalar resonance that decays to a 260 GeV scalar resonance and a Z boson, which decays to charged leptons.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy narrow-width pseudoscalar resonance that decays to a 300 GeV scalar resonance and a Z boson, which decays to charged leptons.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy narrow-width pseudoscalar resonance that decays to a 300 GeV scalar resonance and a Z boson, which decays to charged leptons.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy narrow-width pseudoscalar resonance that decays to a 400 GeV scalar resonance and a Z boson, which decays to charged leptons.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy narrow-width pseudoscalar resonance that decays to a 400 GeV scalar resonance and a Z boson, which decays to charged leptons.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy large-width pseudoscalar resonance that decays to a 260 GeV scalar resonance and a Z boson, which decays to charged leptons.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy large-width pseudoscalar resonance that decays to a 260 GeV scalar resonance and a Z boson, which decays to charged leptons.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy large-width pseudoscalar resonance that decays to a 300 GeV scalar resonance and a Z boson, which decays to charged leptons.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy large-width pseudoscalar resonance that decays to a 300 GeV scalar resonance and a Z boson, which decays to charged leptons.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy large-width pseudoscalar resonance that decays to a 400 GeV scalar resonance and a Z boson, which decays to charged leptons.

Acceptance times efficiency as a function of pseudoscalar resonant mass for each event selection step in the search for a neutral heavy large-width pseudoscalar resonance that decays to a 400 GeV scalar resonance and a Z boson, which decays to charged leptons.

Expected and observed 95% CL upper limits on the cross-section of resonant $H\to 4b$ production in association with a W boson. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits are shown.

Expected and observed 95% CL upper limits on the cross-section of resonant $H\to 4b$ production in association with a W boson. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits are shown.

Expected and observed 95% CL upper limits on the cross-section of resonant $H\to 4b$ production in association with a Z boson. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits are shown.

Expected and observed 95% CL upper limits on the cross-section of resonant $H\to 4b$ production in association with a Z boson. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits are shown.

Expected 95% CL upper limits on the cross-section of a heavy narrow-width pseudoscalar resonance decaying to a Z boson and a heavy scalar resonance decaying to $H\to 4b$. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits are shown.

Expected 95% CL upper limits on the cross-section of a heavy narrow-width pseudoscalar resonance decaying to a Z boson and a heavy scalar resonance decaying to $H\to 4b$. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits are shown.

Observed 95% CL upper limits on the cross-section of a heavy narrow-width pseudoscalar resonance decaying to a Z boson and a heavy scalar resonance decaying to $H\to 4b$. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits are shown.

Observed 95% CL upper limits on the cross-section of a heavy narrow-width pseudoscalar resonance decaying to a Z boson and a heavy scalar resonance decaying to $H\to 4b$. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits are shown.

Expected 95% CL upper limits on the cross-section of a heavy large-width pseudoscalar resonance decaying to a Z boson and a heavy scalar resonance decaying to $H\to 4b$. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits are shown.

Expected 95% CL upper limits on the cross-section of a heavy large-width pseudoscalar resonance decaying to a Z boson and a heavy scalar resonance decaying to $H\to 4b$. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits are shown.

Expected 95% CL upper limits on the cross-section of a heavy large-width pseudoscalar resonance decaying to a Z boson and a heavy scalar resonance decaying to $H\to 4b$. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits are shown.

Observed 95% CL upper limits on the cross-section of a heavy large-width pseudoscalar resonance decaying to a Z boson and a heavy scalar resonance decaying to $H\to 4b$. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits are shown.

Data and post-fit signal and background from S+B fit for 315 GeV resonant $H\to 4b$ production in association with a W boson.

Data and post-fit signal and background from S+B fit for 315 GeV resonant $H\to 4b$ production in association with a W boson.

Data and post-fit signal and background from S+B fit for 400 GeV resonant $H\to 4b$ production in association with a W boson.

Data and post-fit signal and background from S+B fit for 400 GeV resonant $H\to 4b$ production in association with a W boson.

Data and post-fit signal and background from S+B fit for 550 GeV resonant $H\to 4b$ production in association with a Z boson.

Data and post-fit signal and background from S+B fit for 550 GeV resonant $H\to 4b$ production in association with a Z boson.

Data and post-fit signal and background from S+B fit for 400 GeV resonant $H\to 4b$ production in association with a Z boson.

Data and post-fit signal and background from S+B fit for 400 GeV resonant $H\to 4b$ production in association with a Z boson.

Data and post-fit signal and background from S+B fit for a 790 GeV narrow-width pseudoscalar resonance decaying to a Z boson and a 300 GeV scalar resonance decaying to $H\to 4b$.

Data and post-fit signal and background from S+B fit for a 790 GeV narrow-width pseudoscalar resonance decaying to a Z boson and a 300 GeV scalar resonance decaying to $H\to 4b$.

Data and post-fit signal and background from S+B fit for a 420 GeV large-width pseudoscalar resonance decaying to a Z boson and a 320 GeV scalar resonance decaying to $H\to 4b$.

Data and post-fit signal and background from S+B fit for a 420 GeV large-width pseudoscalar resonance decaying to a Z boson and a 320 GeV scalar resonance decaying to $H\to 4b$.

Data and post-fit signal and background from S+B fit for a 700 GeV large-width pseudoscalar resonance decaying to a Z boson and a 380 GeV scalar resonance decaying to $H\to 4b$.

Data and post-fit signal and background from S+B fit for a 700 GeV large-width pseudoscalar resonance decaying to a Z boson and a 380 GeV scalar resonance decaying to $H\to 4b$.

Data and post-fit signal and background from S+B fit for SM VHH production, with each Higgs boson decaying to $2b$.

Data and post-fit signal and background from S+B fit for SM VHH production, with each Higgs boson decaying to $2b$.

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$.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Naming convention for the observables at different levels of the analysis. At the background-subtracted level the contributions of tracks from pile-up collisions and tracks from secondary vertices are subtracted. At the corrected level the tracking-efficiency correction (TEC) is applied. The observables at particle level are the analysis results.

The total pile-up scale-factor relative uncertainty parameterised in $\mu$ and $n_\text{trk,out}$ and expressed in percent.

The $\chi^2$ and NDF for measured normalised differential cross-sections obtained by comparing the different predictions with the unfolded data. Global($n_\text{ch},\Sigma_{n_{\text{ch}}} p_{\text{T}}$) denotes the scenario in which the covariance matrix is built including the correlations of systematic uncertainties between the two observables $n_{\text{ch}}$ and $\Sigma_{n_{\text{ch}}} p_{\text{T}}$

Normalised differential cross-section as a function of $n_\text{ch}$.

Normalised differential cross-section as a function of $\sum_{n_{\text{ch}}} p_{\text{T}}$.

Normalised double-differential cross-section as a function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$.

Normalised double-differential cross-section as a function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$.

Normalised double-differential cross-section as a function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$.

Normalised double-differential cross-section as a function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$.

Normalised double-differential cross-section as a function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n\text{ch}$ in $ n_\text{ch} \geq 80$.

The $\chi^2$ and NDF for the measured normalised differential cross-sections obtained by comparing the different predictions with the unfolded data. The values corresponding to the `Total' scenario are given. In this scenario, only uncertainties described in Section 8 of the paper are included, while the extra theory uncertainties defined in Section 9 are omitted.

The $\chi^2$ and NDF for the measured normalised differential cross-sections obtained by comparing the different predictions with the unfolded data. The values corresponding to the `De-correlate modelling' scenario are given. This scenario uses the detector covariance matrix, and adding modelling uncertainties, scale variations in the matrix element and parton shower as well as the $h_{\text{damp}}$ variation only to the diagonal elements of the covariance matrix.

The $\chi^2$ and NDF for measured absolute differential cross-sections obtained by comparing the different predictions with the unfolded data. Global($n_\text{ch},\Sigma_{n_{\text{ch}}} p_{\text{T}}$) denotes the scenario in which the covariance matrix is built including the correlations of systematic uncertainties between the two observables $n_{\text{ch}}$ and $\Sigma_{n_{\text{ch}}} p_{\text{T}}$

The $\chi^2$ and NDF for the measured absolute differential cross-sections obtained by comparing the different predictions with the unfolded data. The values corresponding to the `Total' scenario are given. In this scenario, only uncertainties described in Section 8 of the paper are included, while the extra theory uncertainties defined in Section 9 are omitted.

The $\chi^2$ and NDF for the measured normalised differential cross-sections obtained by comparing the different predictions with the unfolded data. The values corresponding to the `De-correlate modelling' scenario are given. This scenario uses the detector covariance matrix, and adding modelling uncertainties, scale variations in the matrix element and parton shower as well as the $h_{\text{damp}}$ variation only to the diagonal elements of the covariance matrix.

Absolute differential cross-section as a function of $n_\text{ch}$.

Absolute differential cross-section as a function of $\sum_{n_{\text{ch}}} p_{\text{T}}$.

Absolute double-differential cross-section as a function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$.

Absolute double-differential cross-section as a function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$.

Absolute double-differential cross-section as a function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$.

Absolute double-differential cross-section as a function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$.

Absolute double-differential cross-section as a function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n\text{ch}$ in $ n_\text{ch} \geq 80$.

Covariance matrix of the normalised differential cross-section as function of $n_\text{ch}$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix of the normalised differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ n_\text{ch} < 20$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $40 \leq n_\text{ch} < 60$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $60 \leq n_\text{ch} < 80$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $40 \leq n_\text{ch} < 60$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $60 \leq n_\text{ch} < 80$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $40 \leq n_\text{ch} < 60$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $60 \leq n_\text{ch} < 80$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $40 \leq n_\text{ch} < 60$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $60 \leq n_\text{ch} < 80$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $40 \leq n_\text{ch} < 60$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $60 \leq n_\text{ch} < 80$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Global covariance matrix of the normalised differential cross-section as function of $n_\text{ch}$ and $\sum_{n_{\text{ch}}} p_{\text{T}}$ in $n_\text{ch}$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Global covariance matrix of the normalised differential cross-section as function of $n_\text{ch}$ and $\sum_{n_{\text{ch}}} p_{\text{T}}$ in $n_\text{ch}$ vs. $\sum_{n_{\text{ch}}} p_{\text{T}}$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Global covariance matrix of the normalised differential cross-section as function of $n_\text{ch}$ and $\sum_{n_{\text{ch}}} p_{\text{T}}$ in $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Global covariance matrix of the normalised differential cross-section as function of $n_\text{ch}$ and $\sum_{n_{\text{ch}}} p_{\text{T}}$ in $\sum_{n_{\text{ch}}} p_{\text{T}}$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix of the normalised differential cross-section as function of $n_\text{ch}$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix of the normalised differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ n_\text{ch} < 20$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $40 \leq n_\text{ch} < 60$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $60 \leq n_\text{ch} < 80$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $40 \leq n_\text{ch} < 60$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $60 \leq n_\text{ch} < 80$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $40 \leq n_\text{ch} < 60$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $60 \leq n_\text{ch} < 80$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $40 \leq n_\text{ch} < 60$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $60 \leq n_\text{ch} < 80$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $40 \leq n_\text{ch} < 60$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $60 \leq n_\text{ch} < 80$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ and the normalised double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Global covariance matrix of the normalised differential cross-section as function of $n_\text{ch}$ and $\sum_{n_{\text{ch}}} p_{\text{T}}$ in $n_\text{ch}$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Global covariance matrix of the normalised differential cross-section as function of $n_\text{ch}$ and $\sum_{n_{\text{ch}}} p_{\text{T}}$ in $n_\text{ch}$ vs. $\sum_{n_{\text{ch}}} p_{\text{T}}$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Global covariance matrix of the normalised differential cross-section as function of $n_\text{ch}$ and $\sum_{n_{\text{ch}}} p_{\text{T}}$ in $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Global covariance matrix of the normalised differential cross-section as function of $n_\text{ch}$ and $\sum_{n_{\text{ch}}} p_{\text{T}}$ in $\sum_{n_{\text{ch}}} p_{\text{T}}$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix of the absolute differential cross-section as function of $n_\text{ch}$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix of the absolute differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ n_\text{ch} < 20$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $40 \leq n_\text{ch} < 60$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $60 \leq n_\text{ch} < 80$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $40 \leq n_\text{ch} < 60$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $60 \leq n_\text{ch} < 80$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $40 \leq n_\text{ch} < 60$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $60 \leq n_\text{ch} < 80$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $40 \leq n_\text{ch} < 60$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $60 \leq n_\text{ch} < 80$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $40 \leq n_\text{ch} < 60$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $60 \leq n_\text{ch} < 80$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Global covariance matrix of the absolute differential cross-section as function of $n_\text{ch}$ and $\sum_{n_{\text{ch}}} p_{\text{T}}$ in $n_\text{ch}$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Global covariance matrix of the absolute differential cross-section as function of $n_\text{ch}$ and $\sum_{n_{\text{ch}}} p_{\text{T}}$ in $n_\text{ch}$ vs. $\sum_{n_{\text{ch}}} p_{\text{T}}$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Global covariance matrix of the absolute differential cross-section as function of $n_\text{ch}$ and $\sum_{n_{\text{ch}}} p_{\text{T}}$ in $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Global covariance matrix of the absolute differential cross-section as function of $n_\text{ch}$ and $\sum_{n_{\text{ch}}} p_{\text{T}}$ in $\sum_{n_{\text{ch}}} p_{\text{T}}$ at particle level, accounting for the statistical, systematic uncertainties, and uncertainties in the theoretical predictions.

Covariance matrix of the absolute differential cross-section as function of $n_\text{ch}$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix of the absolute differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ n_\text{ch} < 20$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $40 \leq n_\text{ch} < 60$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $60 \leq n_\text{ch} < 80$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $40 \leq n_\text{ch} < 60$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $60 \leq n_\text{ch} < 80$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $40 \leq n_\text{ch} < 60$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $60 \leq n_\text{ch} < 80$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 40 \leq n_\text{ch} < 60$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $40 \leq n_\text{ch} < 60$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $60 \leq n_\text{ch} < 80$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 60 \leq n_\text{ch} < 80$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} < 20$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $ 20 \leq n_\text{ch} < 40$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $40 \leq n_\text{ch} < 60$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $60 \leq n_\text{ch} < 80$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Covariance matrix between the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ and the absolute double-differential cross-section as function of $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ in $n_\text{ch} \geq 80$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Global covariance matrix of the absolute differential cross-section as function of $n_\text{ch}$ and $\sum_{n_{\text{ch}}} p_{\text{T}}$ in $n_\text{ch}$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Global covariance matrix of the absolute differential cross-section as function of $n_\text{ch}$ and $\sum_{n_{\text{ch}}} p_{\text{T}}$ in $n_\text{ch}$ vs. $\sum_{n_{\text{ch}}} p_{\text{T}}$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Global covariance matrix of the absolute differential cross-section as function of $n_\text{ch}$ and $\sum_{n_{\text{ch}}} p_{\text{T}}$ in $\sum_{n_{\text{ch}}} p_{\text{T}}$ vs. $n_\text{ch}$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

Global covariance matrix of the absolute differential cross-section as function of $n_\text{ch}$ and $\sum_{n_{\text{ch}}} p_{\text{T}}$ in $\sum_{n_{\text{ch}}} p_{\text{T}}$ at particle level, accounting for the statistical and systematic uncertainties, but excluding uncertainties in the theoretical predictions.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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