Acceptance times efficiency for signal grid in merged two-prong signal region.

Acceptance times efficiency for signal grid in resolved two-prong signal region.

Acceptance times efficiency for signal grid in resolved two-prong signal region.

The obtained p-value from comparing data to background estimation across all bins, defined by windows in the X and Y particle masses, in the anomaly signal region.

The obtained p-value from comparing data to background estimation across all bins, defined by windows in the X and Y particle masses, in the anomaly signal region.

The expected 95% CL limits on the cross-section $\sigma(pp \rightarrow Y \rightarrow XH \rightarrow q\bar{q}b\bar{b}$) in pb in the two-dimensional space of $m_Y$ versus $m_X$, obtained from a simultaneous fit of both merged and resolved two-prong signal regions with all statistical and systematic uncertainties.

The expected 95% CL limits on the cross-section $\sigma(pp \rightarrow Y \rightarrow XH \rightarrow q\bar{q}b\bar{b}$) in pb in the two-dimensional space of $m_Y$ versus $m_X$, obtained from a simultaneous fit of both merged and resolved two-prong signal regions with all statistical and systematic uncertainties.

The observed 95% CL limits on the cross-section $\sigma(pp \rightarrow Y \rightarrow XH \rightarrow q\bar{q}b\bar{b}$) in pb in the two-dimensional space of $m_Y$ versus $m_X$, obtained from a simultaneous fit of both merged and resolved two-prong signal regions with all statistical and systematic uncertainties.

The observed 95% CL limits on the cross-section $\sigma(pp \rightarrow Y \rightarrow XH \rightarrow q\bar{q}b\bar{b}$) in pb in the two-dimensional space of $m_Y$ versus $m_X$, obtained from a simultaneous fit of both merged and resolved two-prong signal regions with all statistical and systematic uncertainties.

Average timing distributions for SR data and the estimated background as determined by the background-only fit, in each of the five exclusive $\rho$ categories. For comparison, the expected timing shapes for a few different signal models are superimposed, with each model labeled by the values of the $\tilde\chi^0_1$ mass and lifetime, as well as decay mode. To provide some indication of the variations in signal yield and shape, three signal models are shown for each of the $\tilde\chi^0_1$ decay modes, namely $\tilde\chi^0_1$ $\rightarrow$ $H \tilde G$ and $\tilde\chi^0_1$ $\rightarrow$ $Z \tilde G$. The models shown include a rather low $\tilde\chi^0_1$ mass value of 135 GeV for lifetimes of either 2 ns or 10 ns, and a higher $\tilde\chi^0_1$ mass value which is near the 95% CL exclusion limit for each decay mode for a lifetime of 2 ns. Each signal model is shown with the signal normalization corresponding to a BR value of unity for the decay mode in question.

Average timing distributions for SR data and the estimated background as determined by the background-only fit, in each of the five exclusive $\rho$ categories. For comparison, the expected timing shapes for a few different signal models are superimposed, with each model labeled by the values of the $\tilde\chi^0_1$ mass and lifetime, as well as decay mode. To provide some indication of the variations in signal yield and shape, three signal models are shown for each of the $\tilde\chi^0_1$ decay modes, namely $\tilde\chi^0_1$ $\rightarrow$ $H \tilde G$ and $\tilde\chi^0_1$ $\rightarrow$ $Z \tilde G$. The models shown include a rather low $\tilde\chi^0_1$ mass value of 135 GeV for lifetimes of either 2 ns or 10 ns, and a higher $\tilde\chi^0_1$ mass value which is near the 95% CL exclusion limit for each decay mode for a lifetime of 2 ns. Each signal model is shown with the signal normalization corresponding to a BR value of unity for the decay mode in question.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ branching ratio to the SM Higgs boson, where the assumed cross-section is for higgsino production, and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 - $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$). Several signal hypotheses are overlaid that are labelled by the $\tilde\chi^0_1$ mass, all with a fixed $\tilde\chi^0_1$ lifetime of 2 ns.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ branching ratio to the SM Higgs boson, where the assumed cross-section is for higgsino production, and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 - $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$). Several signal hypotheses are overlaid that are labelled by the $\tilde\chi^0_1$ mass, all with a fixed $\tilde\chi^0_1$ lifetime of 2 ns.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ branching ratio to the SM Higgs boson, where the assumed cross-section is for higgsino production, and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 - $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$). Several signal hypotheses are overlaid that are labelled by the $\tilde\chi^0_1$ mass, all with a fixed $\tilde\chi^0_1$ lifetime of 2 ns.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ branching ratio to the SM Higgs boson, where the assumed cross-section is for higgsino production, and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 - $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$). Several signal hypotheses are overlaid that are labelled by the $\tilde\chi^0_1$ mass, all with a fixed $\tilde\chi^0_1$ lifetime of 2 ns.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

Cutflow for an example higgsino signal with mass 225 GeV and lifetime 10 ns, in the H decay mode. Acceptance is defined at truth level, and efficiency compares the events passing at reconstruction level with respect to truth.

Cutflow for an example higgsino signal with mass 225 GeV and lifetime 10 ns, in the Z decay mode. Acceptance is defined at truth level, and efficiency compares the events passing at reconstruction level with respect to truth.

Cutflow for an example higgsino signal with mass 225 GeV and lifetime 2 ns, in the H decay mode. Acceptance is defined at truth level, and efficiency compares the events passing at reconstruction level with respect to truth.

Cutflow for an example higgsino signal with mass 225 GeV and lifetime 2 ns, in the Z decay mode. Acceptance is defined at truth level, and efficiency compares the events passing at reconstruction level with respect to truth.

Acceptance across the H decay mode signal grid, calculated using truth information. The selection applied corresponds to the model-independent signal region (i.e. the standard SR with $t_{\text{avg}$ > 0.9 ns).

Acceptance across the Z decay mode signal grid, calculated using truth information. The selection applied corresponds to the model-independent signal region (i.e. the standard SR with $t_{\text{avg}$ > 0.9 ns).

Efficiency across the H decay mode signal grid, calculated using reco information. The selection applied corresponds to the model-independent signal region (i.e. the standard SR with $t_{\text{avg}$ > 0.9 ns). Here, the numerator is the signal yield passing the reco selection and the denominator is the signal yield passing the truth selection.

Efficiency across the Z decay mode signal grid, calculated using reco information. The selection applied corresponds to the model-independent signal region (i.e. the standard SR with $t_{\text{avg}$ > 0.9 ns). Here, the numerator is the signal yield passing the reco selection and the denominator is the signal yield passing the truth selection.

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the High-pT jet SR for observed data events

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the High-pT jet SR for expected signal events in the strong gluino pair pair production model with m(gluino)=1.8 TeV, m(chi0)=0.2 TeV, tau(chi0)=0.1 ns

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the Trackless jet SR for observed data events

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the Trackless jet SR for expected signal events in the electroweak pair production model

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Exclusion limits at 95% CL on the production cross section in the electroweak pair production model.

Exclusion limits at 95% CL on the production cross section in the strong gluino pair production models and m(gluino)=2.4 TeV

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Expected exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Observed exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Expected exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Observed exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Exclusion limits at 95% CL on the production cross section in the strong gluino pair production models and m($ ilde{\chi}^0_1$)=1.25 TeV

Acceptance cutflow for the High-pT SR for representative points in the strong gluino pair production model. See additional resources for more information.

Acceptance cutflow for the Trackless SR for representative points in the electroweak pair production model. See additional resources for more information.

Acceptance cutflow for the Trackless SR for representative points in the electroweak pair production model with heavy-flavor quarks final state. See additional resources for more information.

Acceptance cutflow for the High-pT SR for representative points in the electroweak pair production model with heavy-flavor quarks final state. See additional resources for more information.

Reinterpretation Material: Event-level Efficiency for HighPt SR selections, R < 1150 mm

Reinterpretation Material: Event-level Efficiency for HighPt SR selections, R [1150, 3870] mm

Reinterpretation Material: Event-level Efficiency for HighPt SR selections, R > 3870 mm

Reinterpretation Material: Event-level Efficiency for Trackless SR selections, R < 1150 mm

Reinterpretation Material: Event-level Efficiency for Trackless SR selections, R [1150, 3870] mm

Reinterpretation Material: Event-level Efficiency for Trackless SR selections, R > 3870 mm

Reinterpretation Material: Vertex-level Efficiency for R < 22 mm

Reinterpretation Material: Vertex-level Efficiency for R [22, 25] mm

Reinterpretation Material: Vertex-level Efficiency for R [25, 29] mm

Reinterpretation Material: Vertex-level Efficiency for R [29, 38] mm

Reinterpretation Material: Vertex-level Efficiency for R [38, 46] mm

Reinterpretation Material: Vertex-level Efficiency for R [46, 73] mm

Reinterpretation Material: Vertex-level Efficiency for R [73, 84] mm

Reinterpretation Material: Vertex-level Efficiency for R [84, 111] mm

Reinterpretation Material: Vertex-level Efficiency for R [111, 120] mm

Reinterpretation Material: Vertex-level Efficiency for R [120, 145] mm

Reinterpretation Material: Vertex-level Efficiency for R [145, 180] mm

Reinterpretation Material: Vertex-level Efficiency for R [180, 300] mm

Cutflow (acceptance x efficiency) for the High-pT SR for representative points in the strong gluino pair production model. See additional resources for more information.

Cutflow (acceptance x efficiency) for the Trackless SR for representative points in the electroweak pair production model. See additional resources for more information.

Cutflow (acceptance x efficiency) for the Trackless SR for representative points in the electroweak pair production model with heavy-flavor quarks. See additional resources for more information.

Cutflow (acceptance x efficiency) for the High-pT SR for representative points in the electroweak pair production model with heavy-flavor quarks. See additional resources for more information.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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.

Observed (points with error bars) and expected background (solid histograms) distributions for $E_{T}^{miss}$ in the signal region (a) SRL, (b) SRM and (c) SRH after the background-only fit applied to the CRs. The predicted signal distributions for the two models with a gluino mass of 2000 GeV and neutralino mass of 250 GeV (SRL), 1050 GeV (SRM) or 1950 GeV (SRH) are also shown for comparison. The uncertainties in the SM background are only statistical.

Observed and expected exclusion limit in the gluino-neutralino mass plane at 95% CL combined using the signal region with the best expected sensitivity at each point, for the full Run-2 dataset corresponding to an integrated luminosity of $139~\mathrm{fb}^{-1}$, for $\gamma/Z$ (a) and $\gamma/h$ (b) signal models. The black solid line corresponds to the expected limits at 95% CL, with the light (yellow) bands indicating the 1$\sigma$ exclusions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves, the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties. For each point in the higgsino-bino parameter space, the labels indicate the best-expected signal region, where L, M and H mean SRL, SRM and SRH, respectively.

Observed and expected exclusion limit in the gluino-neutralino mass plane at 95% CL combined using the signal region with the best expected sensitivity at each point, for the full Run-2 dataset corresponding to an integrated luminosity of $139~\mathrm{fb}^{-1}$, for $\gamma/Z$ (a) and $\gamma/h$ (b) signal models. The black solid line corresponds to the expected limits at 95% CL, with the light (yellow) bands indicating the 1$\sigma$ exclusions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves, the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties. For each point in the higgsino-bino parameter space, the labels indicate the best-expected signal region, where L, M and H mean SRL, SRM and SRH, respectively.

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/Z$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Acceptance (left) and efficiency (right) for the $\gamma/h$ model signal grid for SRL (top), SRM (middle) and SRH (bottom).

Cutflow for the SRL selection, for two relevant signal points for both $\gamma/Z$ and $\gamma/h$ models, where the gluinos have mass of 2000 GeV and the neutralinos have a mass of 250 GeV (10000 generated events). The numbers are normalized to a luminosity of 139 $fb^{-1}$.

Cutflow for the SRM selection, for two relevant signal points for both $\gamma/Z$ and $\gamma/h$ models, where the gluinos have mass of 2000 GeV and the neutralinos have a mass of 1050 GeV (10000 generated events). The numbers are normalized to a luminosity of 139 $fb^{-1}$.

Cutflow for the SRH selection, for two relevant signal points for both $\gamma/Z$ and $\gamma/h$ models, where the gluinos have mass of 2000 GeV and the neutralinos have a mass of 1950 GeV (10000 generated events). The numbers are normalized to a luminosity of 139 $fb^{-1}$.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected exclusion limits in the gluino–neutralino mass plane at 95% CL for the full Run-2 dataset corresponding to an integrated luminosity of 139 fb−1 , for the (a) $\gamma/Z$ and (b) $\gamma/h$ signal models. They are obtained by combining limits from the signal region with the best expected sensitivity at each point. The dashed (black) line corresponds to the expected limits at 95% CL, with the light (yellow) band indicating the $\pm 1\sigma$ excursions due to experimental and background-theory uncertainties. The observed limits are indicated by medium (red) curves: the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross section by the theoretical scale and PDF uncertainties.

Observed and expected 95% CL lower limits on the $T$-quark mass as a function of the $T$-quark width-to-mass ratio and the branching fraction of the $T \rightarrow Ht$ decay ($\Gamma_{T}$ is the $T$-quark width).

Cutflow table listing the number of events passing each criterion for a $T$-quark hypothesis with a mass of 1.6 TeV and $\kappa_{T} = 0.5$. The initial signal event yield is the predicted number of $T$-quark events inclusive in the Higgs-boson and top-quark decays for 139 fb$^{-1}$.

Observed 95% CL upper limits on the single $T$-quark production cross-section as a function of the $T$-quark coupling $\kappa_{T}$ and $m_{T}$.

Expected 95% CL upper limits on the single $T$-quark production cross-section as a function of the $T$-quark coupling $\kappa_{T}$ and $m_{T}$.

Observed and expected 95% CL lower limits on the $T$-quark mass as a function of the $T$-quark width-to-mass ratio and the branching fraction of the $T \rightarrow Wb$ decay ($\Gamma_{T}$ is the $T$-quark width).

The BDT distribution of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the low mass region ($m^{*}_{b\bar{b}\gamma\gamma} < 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.881 in the low mass region are discarded.

The BDT distribution (with x-axis zoomed in) of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the low mass region ($m^{*}_{b\bar{b}\gamma\gamma} < 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.881 in the low mass region are discarded. The range of BDT scores is from 0.8 to 1.

The BDT distribution (with x-axis zoomed in) of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the low mass region ($m^{*}_{b\bar{b}\gamma\gamma} < 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.881 in the low mass region are discarded. The range of BDT scores is from 0.8 to 1.

The BDT distribution (with x-axis zoomed in) of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the low mass region ($m^{*}_{b\bar{b}\gamma\gamma} < 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.881 in the low mass region are discarded. The range of BDT scores is from 0.8 to 1.

The BDT distribution (with x-axis zoomed in) of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the low mass region ($m^{*}_{b\bar{b}\gamma\gamma} < 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.881 in the low mass region are discarded. The range of BDT scores is from 0.8 to 1.

The BDT distribution of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the high mass region ($m^{*}_{b\bar{b}\gamma\gamma} > 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.857 in the high mass region are discarded.

The BDT distribution of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the high mass region ($m^{*}_{b\bar{b}\gamma\gamma} > 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.857 in the high mass region are discarded.

The BDT distribution of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the high mass region ($m^{*}_{b\bar{b}\gamma\gamma} > 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.857 in the high mass region are discarded.

The BDT distribution of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the high mass region ($m^{*}_{b\bar{b}\gamma\gamma} > 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.857 in the high mass region are discarded.

The BDT distribution (with x-axis zoomed in) of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the high mass region ($m^{*}_{b\bar{b}\gamma\gamma} > 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.857 in the high mass region are discarded. The range of BDT scores is from 0.8 to 1.

The BDT distribution (with x-axis zoomed in) of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the high mass region ($m^{*}_{b\bar{b}\gamma\gamma} > 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.857 in the high mass region are discarded. The range of BDT scores is from 0.8 to 1.

The BDT distribution (with x-axis zoomed in) of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the high mass region ($m^{*}_{b\bar{b}\gamma\gamma} > 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.857 in the high mass region are discarded. The range of BDT scores is from 0.8 to 1.

The BDT distribution (with x-axis zoomed in) of the di-Higgs ggF signal for two different values of $\kappa_{\lambda}$ and the main backgrounds in the high mass region ($m^{*}_{b\bar{b}\gamma\gamma} > 350$ GeV). Distributions are normalized to unit area. The dotted lines denote the category boundaries. Events with a BDT score below 0.857 in the high mass region are discarded. The range of BDT scores is from 0.8 to 1.

The $BDT_{tot}$ score for the benchmark signal $m_{X}$ = 300 GeV and for the main backgrounds. Distributions are normalized to unit area. The dotted line denotes the event selection threshold. Events with a $BDT_{tot}$ score below 0.85 for $m_{X}$ = 300 GeV are discarded.

The $BDT_{tot}$ score for the benchmark signal $m_{X}$ = 300 GeV and for the main backgrounds. Distributions are normalized to unit area. The dotted line denotes the event selection threshold. Events with a $BDT_{tot}$ score below 0.85 for $m_{X}$ = 300 GeV are discarded.

The $BDT_{tot}$ score for the benchmark signal $m_{X}$ = 300 GeV and for the main backgrounds. Distributions are normalized to unit area. The dotted line denotes the event selection threshold. Events with a $BDT_{tot}$ score below 0.85 for $m_{X}$ = 300 GeV are discarded.

The $BDT_{tot}$ score for the benchmark signal $m_{X}$ = 300 GeV and for the main backgrounds. Distributions are normalized to unit area. The dotted line denotes the event selection threshold. Events with a $BDT_{tot}$ score below 0.85 for $m_{X}$ = 300 GeV are discarded.

The $BDT_{tot}$ score for the benchmark signal $m_{X}$ = 500 GeV and for the main backgrounds. Distributions are normalized to unit area. The dotted line denotes the event selection threshold. Events with a $BDT_{tot}$ score below 0.75 for $m_{X}$ = 500 GeV are discarded.

The $BDT_{tot}$ score for the benchmark signal $m_{X}$ = 500 GeV and for the main backgrounds. Distributions are normalized to unit area. The dotted line denotes the event selection threshold. Events with a $BDT_{tot}$ score below 0.75 for $m_{X}$ = 500 GeV are discarded.

The $BDT_{tot}$ score for the benchmark signal $m_{X}$ = 500 GeV and for the main backgrounds. Distributions are normalized to unit area. The dotted line denotes the event selection threshold. Events with a $BDT_{tot}$ score below 0.75 for $m_{X}$ = 500 GeV are discarded.

The $BDT_{tot}$ score for the benchmark signal $m_{X}$ = 500 GeV and for the main backgrounds. Distributions are normalized to unit area. The dotted line denotes the event selection threshold. Events with a $BDT_{tot}$ score below 0.75 for $m_{X}$ = 500 GeV are discarded.

Distributions of $m_{\gamma\gamma}$ in high mass BDT tight category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in high mass BDT tight category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in high mass BDT tight category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in high mass BDT tight category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in high mass BDT loose category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in high mass BDT loose category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in high mass BDT loose category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in high mass BDT loose category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in low mass BDT tight category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in low mass BDT tight category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in low mass BDT tight category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in low mass BDT tight category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in low mass BDT loose category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in low mass BDT loose category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in low mass BDT loose category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ in low mass BDT loose category for the nonresonant $HH$ search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband.

Distributions of $m_{\gamma\gamma}$ for the selections used for the resonance mass point $m_{X}$ = 300 GeV for the resonant search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband. The scalar resonance signal is scaled to a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ = 370 fb.

Distributions of $m_{\gamma\gamma}$ for the selections used for the resonance mass point $m_{X}$ = 300 GeV for the resonant search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband. The scalar resonance signal is scaled to a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ = 370 fb.

Distributions of $m_{\gamma\gamma}$ for the selections used for the resonance mass point $m_{X}$ = 300 GeV for the resonant search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband. The scalar resonance signal is scaled to a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ = 370 fb.

Distributions of $m_{\gamma\gamma}$ for the selections used for the resonance mass point $m_{X}$ = 300 GeV for the resonant search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband. The scalar resonance signal is scaled to a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ = 370 fb.

Distributions of $m_{\gamma\gamma}$ for the selections used for the resonance mass point $m_{X}$ = 500 GeV for the resonant search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband. The scalar resonance signal is scaled to a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ = 67 fb.

Distributions of $m_{\gamma\gamma}$ for the selections used for the resonance mass point $m_{X}$ = 500 GeV for the resonant search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband. The scalar resonance signal is scaled to a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ = 67 fb.

Distributions of $m_{\gamma\gamma}$ for the selections used for the resonance mass point $m_{X}$ = 500 GeV for the resonant search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband. The scalar resonance signal is scaled to a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ = 67 fb.

Distributions of $m_{\gamma\gamma}$ for the selections used for the resonance mass point $m_{X}$ = 500 GeV for the resonant search. The data-derived fractions of nonresonant $\gamma\gamma$, $\gamma$-jet or jet-$\gamma$, and dijet background are applied and the total background is normalized to the data sideband. The scalar resonance signal is scaled to a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ = 67 fb.

The number of data events observed in the 120 GeV < $m_{\gamma\gamma}$ < 130 GeV window, the number of $HH$ signal events expected for $\kappa_{\lambda}$ = 1 and for $\kappa_{\lambda}$ = 10, and events expected for single Higgs boson production (estimated using MC simulation), as well as for continuum background. For the single Higgs boson, Rest includes VBF, $WH$, $tHqb$, and $tHW$. The values are obtained from a fit of the Asimov data set generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$ = 1. The continuum background component of the Asimov data set is obtained from the fit of the data sideband. The uncertainties in $HH$ signals and single Higgs boson background include the systematic uncertainties discussed in Section 6. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious-signal uncertainty.

The number of data events observed in the 120 GeV < $m_{\gamma\gamma}$ < 130 GeV window, the number of $HH$ signal events expected for $\kappa_{\lambda}$ = 1 and for $\kappa_{\lambda}$ = 10, and events expected for single Higgs boson production (estimated using MC simulation), as well as for continuum background. For the single Higgs boson, Rest includes VBF, $WH$, $tHqb$, and $tHW$. The values are obtained from a fit of the Asimov data set generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$ = 1. The continuum background component of the Asimov data set is obtained from the fit of the data sideband. The uncertainties in $HH$ signals and single Higgs boson background include the systematic uncertainties discussed in Section 6. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious-signal uncertainty.

The number of data events observed in the 120 GeV < $m_{\gamma\gamma}$ < 130 GeV window, the number of $HH$ signal events expected for $\kappa_{\lambda}$ = 1 and for $\kappa_{\lambda}$ = 10, and events expected for single Higgs boson production (estimated using MC simulation), as well as for continuum background. For the single Higgs boson, Rest includes VBF, $WH$, $tHqb$, and $tHW$. The values are obtained from a fit of the Asimov data set generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$ = 1. The continuum background component of the Asimov data set is obtained from the fit of the data sideband. The uncertainties in $HH$ signals and single Higgs boson background include the systematic uncertainties discussed in Section 6. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious-signal uncertainty.

The number of data events observed in the 120 GeV < $m_{\gamma\gamma}$ < 130 GeV window, the number of $HH$ signal events expected for $\kappa_{\lambda}$ = 1 and for $\kappa_{\lambda}$ = 10, and events expected for single Higgs boson production (estimated using MC simulation), as well as for continuum background. For the single Higgs boson, Rest includes VBF, $WH$, $tHqb$, and $tHW$. The values are obtained from a fit of the Asimov data set generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$ = 1. The continuum background component of the Asimov data set is obtained from the fit of the data sideband. The uncertainties in $HH$ signals and single Higgs boson background include the systematic uncertainties discussed in Section 6. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious-signal uncertainty.

Observed and expected limits at 95% CL on the cross section of nonresonant Higgs boson pair production as a function of the Higgs boson self-coupling modifier $\kappa_{\lambda}= \lambda_{HHH}/\lambda^{\textrm{SM}}_{HHH}$. The expected constraints on $\kappa_{\lambda}$ are obtained with a background hypothesis excluding $pp \rightarrow HH$ production. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit due to statistical and systematic uncertainties are also shown. The theory prediction curve represents the scenario where all parameters and couplings are set to their SM values except for $\kappa_{\lambda}$. The uncertainty band of the theory prediction curve shows the cross-section uncertainty.

Observed and expected limits at 95% CL on the cross section of nonresonant Higgs boson pair production as a function of the Higgs boson self-coupling modifier $\kappa_{\lambda}= \lambda_{HHH}/\lambda^{\textrm{SM}}_{HHH}$. The expected constraints on $\kappa_{\lambda}$ are obtained with a background hypothesis excluding $pp \rightarrow HH$ production. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit due to statistical and systematic uncertainties are also shown. The theory prediction curve represents the scenario where all parameters and couplings are set to their SM values except for $\kappa_{\lambda}$. The uncertainty band of the theory prediction curve shows the cross-section uncertainty.

Observed and expected limits at 95% CL on the cross section of nonresonant Higgs boson pair production as a function of the Higgs boson self-coupling modifier $\kappa_{\lambda}= \lambda_{HHH}/\lambda^{\textrm{SM}}_{HHH}$. The expected constraints on $\kappa_{\lambda}$ are obtained with a background hypothesis excluding $pp \rightarrow HH$ production. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit due to statistical and systematic uncertainties are also shown. The theory prediction curve represents the scenario where all parameters and couplings are set to their SM values except for $\kappa_{\lambda}$. The uncertainty band of the theory prediction curve shows the cross-section uncertainty.

Observed and expected limits at 95% CL on the cross section of nonresonant Higgs boson pair production as a function of the Higgs boson self-coupling modifier $\kappa_{\lambda}= \lambda_{HHH}/\lambda^{\textrm{SM}}_{HHH}$. The expected constraints on $\kappa_{\lambda}$ are obtained with a background hypothesis excluding $pp \rightarrow HH$ production. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit due to statistical and systematic uncertainties are also shown. The theory prediction curve represents the scenario where all parameters and couplings are set to their SM values except for $\kappa_{\lambda}$. The uncertainty band of the theory prediction curve shows the cross-section uncertainty.

Values of the negative log-profile-likelihood ratio ($-2ln\Lambda$) as a function of $\kappa_{\lambda}$ evaluated for the combination of all the categories of the nonresonant search. The coupling of the Higgs boson to fermions and gauge bosons is set to SM values in the profile likelihood calculation. The Asimov data set is generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$= 1. All systematic uncertainties, including the theoretical uncertainties in the di-Higgs boson production cross section, are included. The intersections of the solid curves and the horizontal dashed lines indicate the 1$\sigma$ and 2$\sigma$ confidence-level intervals.

Values of the negative log-profile-likelihood ratio ($-2ln\Lambda$) as a function of $\kappa_{\lambda}$ evaluated for the combination of all the categories of the nonresonant search. The coupling of the Higgs boson to fermions and gauge bosons is set to SM values in the profile likelihood calculation. The Asimov data set is generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$= 1. All systematic uncertainties, including the theoretical uncertainties in the di-Higgs boson production cross section, are included. The intersections of the solid curves and the horizontal dashed lines indicate the 1$\sigma$ and 2$\sigma$ confidence-level intervals.

Values of the negative log-profile-likelihood ratio ($-2ln\Lambda$) as a function of $\kappa_{\lambda}$ evaluated for the combination of all the categories of the nonresonant search. The coupling of the Higgs boson to fermions and gauge bosons is set to SM values in the profile likelihood calculation. The Asimov data set is generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$= 1. All systematic uncertainties, including the theoretical uncertainties in the di-Higgs boson production cross section, are included. The intersections of the solid curves and the horizontal dashed lines indicate the 1$\sigma$ and 2$\sigma$ confidence-level intervals.

Values of the negative log-profile-likelihood ratio ($-2ln\Lambda$) as a function of $\kappa_{\lambda}$ evaluated for the combination of all the categories of the nonresonant search. The coupling of the Higgs boson to fermions and gauge bosons is set to SM values in the profile likelihood calculation. The Asimov data set is generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$= 1. All systematic uncertainties, including the theoretical uncertainties in the di-Higgs boson production cross section, are included. The intersections of the solid curves and the horizontal dashed lines indicate the 1$\sigma$ and 2$\sigma$ confidence-level intervals.

The number of events observed in the 120 < $m_{\gamma\gamma}$ < 130 GeV window in data, the number of events expected for scalar resonance signals of masses $m_{X}$ = 300 GeV and $m_{X}$ = 500 GeV assuming a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ equal to the observed exclusion limits of Figure 15, and events expected for SM $HH$ and single Higgs boson production (estimated using MC simulation), as well as for continuum background. The values are obtained from a fit of the Asimov data set generated under the signal-plus-background hypothesis. The continuum background component of the Asimov data set is obtained from the fit of the data sideband. The uncertainties in the resonant signals and the SM $HH$ and single-Higgs-boson backgrounds include the systematic uncertainties discussed in Section 6. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious-signal uncertainty.

The number of events observed in the 120 < $m_{\gamma\gamma}$ < 130 GeV window in data, the number of events expected for scalar resonance signals of masses $m_{X}$ = 300 GeV and $m_{X}$ = 500 GeV assuming a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ equal to the observed exclusion limits of Figure 15, and events expected for SM $HH$ and single Higgs boson production (estimated using MC simulation), as well as for continuum background. The values are obtained from a fit of the Asimov data set generated under the signal-plus-background hypothesis. The continuum background component of the Asimov data set is obtained from the fit of the data sideband. The uncertainties in the resonant signals and the SM $HH$ and single-Higgs-boson backgrounds include the systematic uncertainties discussed in Section 6. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious-signal uncertainty.

The number of events observed in the 120 < $m_{\gamma\gamma}$ < 130 GeV window in data, the number of events expected for scalar resonance signals of masses $m_{X}$ = 300 GeV and $m_{X}$ = 500 GeV assuming a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ equal to the observed exclusion limits of Figure 15, and events expected for SM $HH$ and single Higgs boson production (estimated using MC simulation), as well as for continuum background. The values are obtained from a fit of the Asimov data set generated under the signal-plus-background hypothesis. The continuum background component of the Asimov data set is obtained from the fit of the data sideband. The uncertainties in the resonant signals and the SM $HH$ and single-Higgs-boson backgrounds include the systematic uncertainties discussed in Section 6. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious-signal uncertainty.

The number of events observed in the 120 < $m_{\gamma\gamma}$ < 130 GeV window in data, the number of events expected for scalar resonance signals of masses $m_{X}$ = 300 GeV and $m_{X}$ = 500 GeV assuming a total production cross section $\sigma(pp \rightarrow X \rightarrow HH)$ equal to the observed exclusion limits of Figure 15, and events expected for SM $HH$ and single Higgs boson production (estimated using MC simulation), as well as for continuum background. The values are obtained from a fit of the Asimov data set generated under the signal-plus-background hypothesis. The continuum background component of the Asimov data set is obtained from the fit of the data sideband. The uncertainties in the resonant signals and the SM $HH$ and single-Higgs-boson backgrounds include the systematic uncertainties discussed in Section 6. The uncertainty in the continuum background is given by the sum in quadrature of the statistical uncertainty from the fit to the data and the spurious-signal uncertainty.

Observed and expected limits at 95% CL on the production cross section of a narrow-width scalar resonance $X$ as a function of the mass $m_{X}$ of the hypothetical scalar particle. The black solid line represents the observed upper limits. The dashed line represents the expected upper limits. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit due to statistical and systematic uncertainties are also shown.

Observed and expected limits at 95% CL on the production cross section of a narrow-width scalar resonance $X$ as a function of the mass $m_{X}$ of the hypothetical scalar particle. The black solid line represents the observed upper limits. The dashed line represents the expected upper limits. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit due to statistical and systematic uncertainties are also shown.

Observed and expected limits at 95% CL on the production cross section of a narrow-width scalar resonance $X$ as a function of the mass $m_{X}$ of the hypothetical scalar particle. The black solid line represents the observed upper limits. The dashed line represents the expected upper limits. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit due to statistical and systematic uncertainties are also shown.

Observed and expected limits at 95% CL on the production cross section of a narrow-width scalar resonance $X$ as a function of the mass $m_{X}$ of the hypothetical scalar particle. The black solid line represents the observed upper limits. The dashed line represents the expected upper limits. The $\pm 1\sigma$ and $\pm 2\sigma$ variations about the expected limit due to statistical and systematic uncertainties are also shown.

Breakdown of the dominant systematic uncertainties. The impact of the uncertainties is defined according to the statistical analysis described in Section 7. It corresponds to the relative variation of the expected upper limit on the cross section when re-evaluating the profile likelihood ratio after fixing the nuisance parameter in question to its best-fit value, while all remaining nuisance parameters remain free to float. The impact is shown in %. Only systematic uncertainties with an impact of at least 0.2% are shown. Uncertainties of the "Norm. + Shape" type affect both the normalization and the parameters of the functional form. The rest of the uncertainties affect only the yields.

Breakdown of the dominant systematic uncertainties. The impact of the uncertainties is defined according to the statistical analysis described in Section 7. It corresponds to the relative variation of the expected upper limit on the cross section when re-evaluating the profile likelihood ratio after fixing the nuisance parameter in question to its best-fit value, while all remaining nuisance parameters remain free to float. The impact is shown in %. Only systematic uncertainties with an impact of at least 0.2% are shown. Uncertainties of the "Norm. + Shape" type affect both the normalization and the parameters of the functional form. The rest of the uncertainties affect only the yields.

Breakdown of the dominant systematic uncertainties. The impact of the uncertainties is defined according to the statistical analysis described in Section 7. It corresponds to the relative variation of the expected upper limit on the cross section when re-evaluating the profile likelihood ratio after fixing the nuisance parameter in question to its best-fit value, while all remaining nuisance parameters remain free to float. The impact is shown in %. Only systematic uncertainties with an impact of at least 0.2% are shown. Uncertainties of the "Norm. + Shape" type affect both the normalization and the parameters of the functional form. The rest of the uncertainties affect only the yields.

Breakdown of the dominant systematic uncertainties. The impact of the uncertainties is defined according to the statistical analysis described in Section 7. It corresponds to the relative variation of the expected upper limit on the cross section when re-evaluating the profile likelihood ratio after fixing the nuisance parameter in question to its best-fit value, while all remaining nuisance parameters remain free to float. The impact is shown in %. Only systematic uncertainties with an impact of at least 0.2% are shown. Uncertainties of the "Norm. + Shape" type affect both the normalization and the parameters of the functional form. The rest of the uncertainties affect only the yields.

Cutflow for nonresonant di-Higgs ggF signal sample, yields are normalized to 139 $fb^{-1}$.

Cutflow for nonresonant di-Higgs ggF signal sample, yields are normalized to 139 $fb^{-1}$.

Cutflow for nonresonant di-Higgs ggF signal sample, yields are normalized to 139 $fb^{-1}$.

Cutflow for nonresonant di-Higgs ggF signal sample, yields are normalized to 139 $fb^{-1}$.

Cutflow for resonant signal sample, with $m_{X}$ = 300 GeV, yields are normalized to 139 $fb^{-1}$.

Cutflow for resonant signal sample, with $m_{X}$ = 300 GeV, yields are normalized to 139 $fb^{-1}$.

Cutflow for resonant signal sample, with $m_{X}$ = 300 GeV, yields are normalized to 139 $fb^{-1}$.

Cutflow for resonant signal sample, with $m_{X}$ = 300 GeV, yields are normalized to 139 $fb^{-1}$.

Cutflow for resonant signal sample, with $m_{X}$ = 500 GeV, yields are normalized to 139 $fb^{-1}$.

Cutflow for resonant signal sample, with $m_{X}$ = 500 GeV, yields are normalized to 139 $fb^{-1}$.

Cutflow for resonant signal sample, with $m_{X}$ = 500 GeV, yields are normalized to 139 $fb^{-1}$.

Cutflow for resonant signal sample, with $m_{X}$ = 500 GeV, yields are normalized to 139 $fb^{-1}$.

Comparison of $m_{b\bar{b}}$ distributions when applying the specific b-jet energy calibration and the nominal jet energy calibration. The distributions are fitted using a Bukin function, and the values of the peak position, resolution and the relative improvement are reported in the legend.

Comparison of $m_{b\bar{b}}$ distributions when applying the specific b-jet energy calibration and the nominal jet energy calibration. The distributions are fitted using a Bukin function, and the values of the peak position, resolution and the relative improvement are reported in the legend.

Comparison of $m_{b\bar{b}}$ distributions when applying the specific b-jet energy calibration and the nominal jet energy calibration. The distributions are fitted using a Bukin function, and the values of the peak position, resolution and the relative improvement are reported in the legend.

Comparison of $m_{b\bar{b}}$ distributions when applying the specific b-jet energy calibration and the nominal jet energy calibration. The distributions are fitted using a Bukin function, and the values of the peak position, resolution and the relative improvement are reported in the legend.

Fit results of $m_{b\bar{b}}$ distributions when applying the specific b-jet energy calibration and the nominal jet energy calibration. The distributions are fitted using a Bukin function, and the values of the peak position, resolution and the relative improvement are reported in the legend.

Fit results of $m_{b\bar{b}}$ distributions when applying the specific b-jet energy calibration and the nominal jet energy calibration. The distributions are fitted using a Bukin function, and the values of the peak position, resolution and the relative improvement are reported in the legend.

Fit results of $m_{b\bar{b}}$ distributions when applying the specific b-jet energy calibration and the nominal jet energy calibration. The distributions are fitted using a Bukin function, and the values of the peak position, resolution and the relative improvement are reported in the legend.

Fit results of $m_{b\bar{b}}$ distributions when applying the specific b-jet energy calibration and the nominal jet energy calibration. The distributions are fitted using a Bukin function, and the values of the peak position, resolution and the relative improvement are reported in the legend.

The relative amount (purity) of expected events from SM $HH$ and single Higgs boson production processes for each of the four categories of the nonresonant search. The Higgs boson pair production with $k_{\lambda} = 1$ is considered as signal in (a), while the case with $k_{\lambda} = 10$ is considered as signal in (b).

The relative amount (purity) of expected events from SM $HH$ and single Higgs boson production processes for each of the four categories of the nonresonant search. The Higgs boson pair production with $k_{\lambda} = 1$ is considered as signal in (a), while the case with $k_{\lambda} = 10$ is considered as signal in (b).

The relative amount (purity) of expected events from SM $HH$ and single Higgs boson production processes for each of the four categories of the nonresonant search. The Higgs boson pair production with $k_{\lambda} = 1$ is considered as signal in (a), while the case with $k_{\lambda} = 10$ is considered as signal in (b).

The relative amount (purity) of expected events from SM $HH$ and single Higgs boson production processes for each of the four categories of the nonresonant search. The Higgs boson pair production with $k_{\lambda} = 1$ is considered as signal in (a), while the case with $k_{\lambda} = 10$ is considered as signal in (b).

The expected significance in each of the four categories of the nonresonant search. The Higgs boson pair production with $k_{\lambda} = 1$ is considered as signal in (a), while the case with $k_{\lambda} = 10$ is considered as signal in (b). The single Higgs boson processes and the di-photon continuum spectrum are considered as background.

The expected significance in each of the four categories of the nonresonant search. The Higgs boson pair production with $k_{\lambda} = 1$ is considered as signal in (a), while the case with $k_{\lambda} = 10$ is considered as signal in (b). The single Higgs boson processes and the di-photon continuum spectrum are considered as background.

The expected significance in each of the four categories of the nonresonant search. The Higgs boson pair production with $k_{\lambda} = 1$ is considered as signal in (a), while the case with $k_{\lambda} = 10$ is considered as signal in (b). The single Higgs boson processes and the di-photon continuum spectrum are considered as background.

The expected significance in each of the four categories of the nonresonant search. The Higgs boson pair production with $k_{\lambda} = 1$ is considered as signal in (a), while the case with $k_{\lambda} = 10$ is considered as signal in (b). The single Higgs boson processes and the di-photon continuum spectrum are considered as background.

Spurious signal result for the exponential pdf for the various ggF nonresonant di-Higgs categories. In each category, the spurious signal value ($n_{sp}$) and its ratio to the expected statistical error ($Z_{spur}$) from data are shown.

Spurious signal result for the exponential pdf for the various ggF nonresonant di-Higgs categories. In each category, the spurious signal value ($n_{sp}$) and its ratio to the expected statistical error ($Z_{spur}$) from data are shown.

Spurious signal result for the exponential pdf for the various ggF nonresonant di-Higgs categories. In each category, the spurious signal value ($n_{sp}$) and its ratio to the expected statistical error ($Z_{spur}$) from data are shown.

Spurious signal result for the exponential pdf for the various ggF nonresonant di-Higgs categories. In each category, the spurious signal value ($n_{sp}$) and its ratio to the expected statistical error ($Z_{spur}$) from data are shown.

Spurious signal result for the exponential pdf as function of the resonant di-Higgs signal mass. The spurious signal value and its ratio to the expected statistical error from data are shown.

Spurious signal result for the exponential pdf as function of the resonant di-Higgs signal mass. The spurious signal value and its ratio to the expected statistical error from data are shown.

Spurious signal result for the exponential pdf as function of the resonant di-Higgs signal mass. The spurious signal value and its ratio to the expected statistical error from data are shown.

Spurious signal result for the exponential pdf as function of the resonant di-Higgs signal mass. The spurious signal value and its ratio to the expected statistical error from data are shown.

Distributions of the signal efficiency as a function of $\kappa_{\lambda}$, for the di-Higgs boson ggF nonresonant production mode. The range of $\kappa_{\lambda}$ in the table is from -10 to 10.

Distributions of the signal efficiency as a function of $\kappa_{\lambda}$, for the di-Higgs boson ggF nonresonant production mode. The range of $\kappa_{\lambda}$ in the table is from -10 to 10.

Distributions of the signal efficiency as a function of $\kappa_{\lambda}$, for the di-Higgs boson ggF nonresonant production mode. The range of $\kappa_{\lambda}$ in the table is from -10 to 10.

Distributions of the signal efficiency as a function of $\kappa_{\lambda}$, for the di-Higgs boson ggF nonresonant production mode. The range of $\kappa_{\lambda}$ in the table is from -10 to 10.

Distributions of the signal efficiency as a function of $\kappa_{\lambda}$, for the di-Higgs boson VBF nonresonant production mode. The range of $\kappa_{\lambda}$ in the table is from -10 to 10.

Distributions of the signal efficiency as a function of $\kappa_{\lambda}$, for the di-Higgs boson VBF nonresonant production mode. The range of $\kappa_{\lambda}$ in the table is from -10 to 10.

Distributions of the signal efficiency as a function of $\kappa_{\lambda}$, for the di-Higgs boson VBF nonresonant production mode. The range of $\kappa_{\lambda}$ in the table is from -10 to 10.

Distributions of the signal efficiency as a function of $\kappa_{\lambda}$, for the di-Higgs boson VBF nonresonant production mode. The range of $\kappa_{\lambda}$ in the table is from -10 to 10.

Values of the negative log-profile-likelihood ($-2ln\Lambda$) as a function of $\kappa_{\lambda}$ evaluated for the combination of all the categories of the nonresonant search. The coupling of the Higgs boson to fermions and gauge bosons is set to SM values in the profile likelihood calculation. The Asimov data set is generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$= 1. All systematic uncertainties, including the theoretical uncertainties on the di-Higgs boson production cross section, are included. The intersections of the solid curves and the horizontal dashed lines indicate the 1$\sigma$ and 2$\sigma$ confidence level intervals. The best fit value corresponds to $\kappa_{\lambda}$ = $2.8^{+2.0}_{-2.2}(^{+3.8}_{-4.3})$ for the 1$\sigma$(2$\sigma$) confidence interval. The expected value corresponds to $\kappa_{\lambda}$ = $1.0^{+5.5}_{-2.4}(^{+7.3}_{-4.2})$ for the 1$\sigma$(2$\sigma$) confidence interval. The dashed curves represent values of the negative log-profile-likelihood where the Higgs boson branching fractions and the cross section of the production modes are varied as a function of $\kappa_{\lambda}$. In this case,the best fit value corresponds to $\kappa_{\lambda}$ = $2.7^{+2.0}_{-2.2}(^{+3.8}_{-4.3})$ and the expected value corresponds to $\kappa_{\lambda}$ = $1.0^{+5.4}_{-2.5}(^{+7.3}_{-4.3})$ for the 1$\sigma$(2$\sigma$) confidence interval.

Values of the negative log-profile-likelihood ($-2ln\Lambda$) as a function of $\kappa_{\lambda}$ evaluated for the combination of all the categories of the nonresonant search. The coupling of the Higgs boson to fermions and gauge bosons is set to SM values in the profile likelihood calculation. The Asimov data set is generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$= 1. All systematic uncertainties, including the theoretical uncertainties on the di-Higgs boson production cross section, are included. The intersections of the solid curves and the horizontal dashed lines indicate the 1$\sigma$ and 2$\sigma$ confidence level intervals. The best fit value corresponds to $\kappa_{\lambda}$ = $2.8^{+2.0}_{-2.2}(^{+3.8}_{-4.3})$ for the 1$\sigma$(2$\sigma$) confidence interval. The expected value corresponds to $\kappa_{\lambda}$ = $1.0^{+5.5}_{-2.4}(^{+7.3}_{-4.2})$ for the 1$\sigma$(2$\sigma$) confidence interval. The dashed curves represent values of the negative log-profile-likelihood where the Higgs boson branching fractions and the cross section of the production modes are varied as a function of $\kappa_{\lambda}$. In this case,the best fit value corresponds to $\kappa_{\lambda}$ = $2.7^{+2.0}_{-2.2}(^{+3.8}_{-4.3})$ and the expected value corresponds to $\kappa_{\lambda}$ = $1.0^{+5.4}_{-2.5}(^{+7.3}_{-4.3})$ for the 1$\sigma$(2$\sigma$) confidence interval.

Values of the negative log-profile-likelihood ($-2ln\Lambda$) as a function of $\kappa_{\lambda}$ evaluated for the combination of all the categories of the nonresonant search. The coupling of the Higgs boson to fermions and gauge bosons is set to SM values in the profile likelihood calculation. The Asimov data set is generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$= 1. All systematic uncertainties, including the theoretical uncertainties on the di-Higgs boson production cross section, are included. The intersections of the solid curves and the horizontal dashed lines indicate the 1$\sigma$ and 2$\sigma$ confidence level intervals. The best fit value corresponds to $\kappa_{\lambda}$ = $2.8^{+2.0}_{-2.2}(^{+3.8}_{-4.3})$ for the 1$\sigma$(2$\sigma$) confidence interval. The expected value corresponds to $\kappa_{\lambda}$ = $1.0^{+5.5}_{-2.4}(^{+7.3}_{-4.2})$ for the 1$\sigma$(2$\sigma$) confidence interval. The dashed curves represent values of the negative log-profile-likelihood where the Higgs boson branching fractions and the cross section of the production modes are varied as a function of $\kappa_{\lambda}$. In this case,the best fit value corresponds to $\kappa_{\lambda}$ = $2.7^{+2.0}_{-2.2}(^{+3.8}_{-4.3})$ and the expected value corresponds to $\kappa_{\lambda}$ = $1.0^{+5.4}_{-2.5}(^{+7.3}_{-4.3})$ for the 1$\sigma$(2$\sigma$) confidence interval.

Values of the negative log-profile-likelihood ($-2ln\Lambda$) as a function of $\kappa_{\lambda}$ evaluated for the combination of all the categories of the nonresonant search. The coupling of the Higgs boson to fermions and gauge bosons is set to SM values in the profile likelihood calculation. The Asimov data set is generated under the SM signal-plus-background hypothesis, $\kappa_{\lambda}$= 1. All systematic uncertainties, including the theoretical uncertainties on the di-Higgs boson production cross section, are included. The intersections of the solid curves and the horizontal dashed lines indicate the 1$\sigma$ and 2$\sigma$ confidence level intervals. The best fit value corresponds to $\kappa_{\lambda}$ = $2.8^{+2.0}_{-2.2}(^{+3.8}_{-4.3})$ for the 1$\sigma$(2$\sigma$) confidence interval. The expected value corresponds to $\kappa_{\lambda}$ = $1.0^{+5.5}_{-2.4}(^{+7.3}_{-4.2})$ for the 1$\sigma$(2$\sigma$) confidence interval. The dashed curves represent values of the negative log-profile-likelihood where the Higgs boson branching fractions and the cross section of the production modes are varied as a function of $\kappa_{\lambda}$. In this case,the best fit value corresponds to $\kappa_{\lambda}$ = $2.7^{+2.0}_{-2.2}(^{+3.8}_{-4.3})$ and the expected value corresponds to $\kappa_{\lambda}$ = $1.0^{+5.4}_{-2.5}(^{+7.3}_{-4.3})$ for the 1$\sigma$(2$\sigma$) confidence interval.

Minimum BDT value of the events passing the selection criteria of the resonant search. The combined BDT score is formed using as coefficients $C_{1}$ = 0.65 and $C_{2}$ = 1 − $C_{1}$. The selection efficiency for the resonant $X \rightarrow HH$ signal is also shown.

Minimum BDT value of the events passing the selection criteria of the resonant search. The combined BDT score is formed using as coefficients $C_{1}$ = 0.65 and $C_{2}$ = 1 − $C_{1}$. The selection efficiency for the resonant $X \rightarrow HH$ signal is also shown.

Minimum BDT value of the events passing the selection criteria of the resonant search. The combined BDT score is formed using as coefficients $C_{1}$ = 0.65 and $C_{2}$ = 1 − $C_{1}$. The selection efficiency for the resonant $X \rightarrow HH$ signal is also shown.

Minimum BDT value of the events passing the selection criteria of the resonant search. The combined BDT score is formed using as coefficients $C_{1}$ = 0.65 and $C_{2}$ = 1 − $C_{1}$. The selection efficiency for the resonant $X \rightarrow HH$ signal is also shown.