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.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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.

Distributions of $s_{\text{T}}$ in the single-tau one-bin SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Distributions of $m_{\text{T}}(\tau)$ in the single-tau one-bin SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Distributions of $\Sigma m_{\text{T}}(b_{1,2})$ in the single-tau $p_{\text{T}}(\tau)$-binned SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Distributions of $p_{\text{T}}(\tau)$ in the single-tau $p_{\text{T}}(\tau)$-binned SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Upper limits on the signal cross section at the 95 % confidence level for the stop-stau signal model.

Upper limits on the signal cross section at the 95 % confidence level for the scalar third-generation leptoquark signal model with up-type leptoquarks.

Upper limits on the signal cross section at the 95 % confidence level for the scalar third-generation leptoquark signal model with down-type leptoquarks.

Upper limits on the signal cross section at the 95 % confidence level for the vector third-generation leptoquark signal model with minimal coupling (MC).

Upper limits on the signal cross section at the 95 % confidence level for the vector third-generation leptoquark signal model with additional gauge couplings (YM).

Acceptance of the one-bin signal region of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the one-bin signal region of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the signal region of the di-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the one-bin signal region of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the signal region of the di-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the one-bin signal region of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the signal region of the di-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the signal region of the di-tau channel for pair production of top squarks with decays via tau sleptons.

Cutflow for the benchmark signal model $m(\tilde{t}_{1}) = 1350$ GeV, $m(\tilde{\tau}_{1}) = 1090$ GeV for the di-tau SR. The simulated sample contains 30,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\tilde{t}_{1}) = 1350$ GeV, $m(\tilde{\tau}_{1}) = 1090$ GeV for the single-tau one-bin SR. The simulated sample contains 30,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\tilde{t}_{1}) = 1350$ GeV, $m(\tilde{\tau}_{1}) = 1090$ GeV for the single-tau multi-bin SR. The simulated sample contains 30,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{u}}) = 1.2$ TeV, $\beta = 0.5$ for the di-tau SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{u}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau one-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{u}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau multi-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{d}}) = 1.2$ TeV, $\beta = 0.5$ for the di-tau SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{d}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau one-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{d}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau multi-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the minimal-coupling scenario for the di-tau SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the minimal-coupling scenario for the single-tau one-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the minimal-coupling scenario for the single-tau multi-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the Yang--Mills scenario for the di-tau SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the Yang--Mills scenario for the single-tau one-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the Yang--Mills scenario for the single-tau multi-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

CHI distribution for mass bin > 1200 GeV.

Centrality Ratio.

Signal region detector-level distribution for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$.

Signal region detector-level distribution for the observable $m_{e\mu}$.

Signal region detector-level distribution for the observable $p_{\mathrm{T}}^{e\mu}$.

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

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.

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.