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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

The $p_{T}$ dependencies of $v_{1}$ within $-0.4<y<-0.1$ for $t$ in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity and $p_{T}$ dependencies of $v_{1}$ for $^{3}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The $p_{T}$ dependencies of $v_{1}$ within $-0.1<y<0$ for $^4$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The $p_{T}$ dependencies of $v_{1}$ within $-0.4<y<-0.1$ for $^{4}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity and $p_{T}$ dependencies of $v_{2}$ for $p$ in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity and $p_{T}$ dependencies of $v_{2}$ for $d$ in 10-40% mid-central Au+Au collisions at 3 GeV.

The $p_{T}$ dependencies of $v_{2}$ within $-0.1<y<0$ for $t$ in 10-40% mid-central Au+Au collisions at 3 GeV.

The $p_{T}$ dependencies of $v_{2}$ within $-0.4<y<-0.1$ for $t$ in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity and $p_{T}$ dependencies of $v_{2}$ for $^{3}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The $p_{T}$ dependencies of $v_{2}$ within $-0.1<y<0$ for $^4$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The $p_{T}$ dependencies of $v_{2}$ within $-0.4<y<-0.1$ for $^{4}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity dependencies of $v_{1}$ for $p$, $d$, and $^{3}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity dependencies of $v_{1}$ for $t$ and $^{4}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity dependencies of $v_{1}$ from JAM plus coalescence calculations in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity dependencies of $v_{2}$ for $p$, $d$, and $^{3}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity dependencies of $v_{2}$ for $t$ and $^{4}$He in 10-40% mid-central Au+Au collisions at 3 GeV.

The rapidity dependencies of $v_{2}$ from JAM plus coalescence calculations in 10-40% mid-central Au+Au collisions at 3 GeV.

Light nucleus scaled $v_{1}$ slopes as a function os collision energy in 10-40 mid-cantral Au+Au collisions.

Light nucleus scaled $v_{1}$ slopes from JAM plus coalescence in 10-40 mid-cantral Au+Au collisions at 3 GeV.

Centrality dependence of the proton cumulant ratios for Au+Au collisions at $\sqrt{s_{NN}}$ = 3.0 GeV. Protons are from $-0.5 < y < 0$ and $0.4 < p_{T} < 2.0$ GeV/$c$. Systematic uncertainties are represented by gray bars. Statistical uncertainties are smaller than marker size. CBWC is applied to all cumulant ratios. While open squares represent the data without the VFC correction, blue triangles and red circles are the results with VFC using the $\langle N_{\rm{part}} \rangle$ distributions from the UrQMD and Glauber models, respectively. UrQMD model results are represented as gold dashed line.

Centrality dependence of the proton cumulant ratios for Au+Au collisions at $\sqrt{s_{NN}}$ = 3.0 GeV. Protons are from $-0.5 < y < 0$ and $0.4 < p_{T} < 2.0$ GeV/$c$. Systematic uncertainties are represented by gray bars. Statistical uncertainties are smaller than marker size. CBWC is applied to all cumulant ratios. While open squares represent the data without the VFC correction, blue triangles and red circles are the results with VFC using the $\langle N_{\rm{part}} \rangle$ distributions from the UrQMD and Glauber models, respectively. UrQMD model results are represented as gold dashed line.

Rapidity and transverse momentum dependence of the proton cumulant ratios for 0–5% central collisions. Black-squares, red-dots and blue-triangles stand for data without and with the VFC using Glauber and UrQMD, respectively.

Rapidity and transverse momentum dependence of the proton cumulant ratios for 0–5% central collisions. Black-squares, red-dots and blue-triangles stand for data without and with the VFC using Glauber and UrQMD, respectively.

Rapidity and transverse momentum dependence of the proton cumulant ratios for 0–5% central collisions. Black-squares, red-dots and blue-triangles stand for data without and with the VFC using Glauber and UrQMD, respectively.

Rapidity and transverse momentum dependence of the proton cumulant ratios for 0–5% central collisions. Black-squares, red-dots and blue-triangles stand for data without and with the VFC using Glauber and UrQMD, respectively.

Collision energy dependence of the ratios of cumulants, $C_4/C_2$, for proton (squares) and net-proton (red circles) from top 0-5% Au+Au collisions at RHIC. The points for protons are shifted horizontally for clarity. The new result for proton from $\sqrt{s_{NN}}$ = 3.0 GeV collisions is shown as a filled square. HADES data of $\sqrt{s_{NN}}$ = 2.4 GeV 0-10% collisions is also shown. The vertical black and gray bars are the statistical and systematic uncertainties, respectively. In addition, results from the HRG model, based on both Canonical Ensemble (CE) and Grand-Canonical Ensemble (GCE), and transport model UrQMD are presented.

The correlation functions (corrected for nonuniform detector efficiency in $\phi$; not corrected for the absolute detection efficiency) vs. azimuthal angle difference between forward ($2.6<\eta<4.0$) $\pi^{0}$s in $p$+$p$ collisions at $\sqrt{s_{\mathrm{_{NN}}}}=200$ GeV at high $p_{T}$ ($p^{trig}_{T}$=2.5-3 GeV/c, $p^{asso}_{T}$=2-2.5 GeV/c)

The correlation functions (corrected for nonuniform detector efficiency in $\phi$; not corrected for the absolute detection efficiency) vs. azimuthal angle difference between forward ($2.6<\eta<4.0$) $\pi^{0}$s in $p$+Al collisions at $\sqrt{s_{\mathrm{_{NN}}}}=200$ GeV at high $p_{T}$ ($p^{trig}_{T}$=2.5-3 GeV/c, $p^{asso}_{T}$=2-2.5 GeV/c)

The correlation functions (corrected for nonuniform detector efficiency in $\phi$; not corrected for the absolute detection efficiency) vs. azimuthal angle difference between forward ($2.6<\eta<4.0$) $\pi^{0}$s in $p$+$Au collisions at $\sqrt{s_{\mathrm{_{NN}}}}=200$ GeV at high $p_{T}$ ($p^{trig}_{T}$=2.5-3 GeV/c, $p^{asso}_{T}$=2-2.5 GeV/c)

Relative area of MinBias $pAu/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

Relative area of MinBias $pAl/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

Relative area of $pAu/pp$ at b=0 of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from prediction

Relative width of MinBias $pAu/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

Relative width of MinBias $pAl/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

Relative pedestal of MinBias $pAu/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

Relative pedestal of MinBias $pAl/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

Relative area of MinBias $pp/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$), the ratio is unity, error is zero

Relative area of MinBias $pA/pp$ of back-to-back di-$\pi^0$ correlations at forward pseudorapidities ($2.6<\eta<4.0$) from data

The efficiency for simulated VBF events to pass each analysis cut.

The efficiency for simulated VH events to pass each analysis cut.

The efficiency for simulated ttH events to pass each analysis cut.

The fractional contribution of each production mode to a given analysis region around the Higgs boson peak, defined as 105 < mJ < 140 GeV. The fraction is given with respect to the total signal yield in the analysis region in question.

Signal acceptance times efficiency within the fiducial volume used in the fiducial region.

Signal acceptance times efficiency for the STXS volumes in the differential measurement. Along with the pTH requirements shown, |yH | < 2 is required. For events with pTH < 300 GeV, the acceptance times efficiency is less than 0.1 × 10−2.

Event yields and associated uncertainties after the global likelihood fit in the pT(H) > 450 GeV fiducial region. The total background yield can differ from the sum of the individual components due to rounding.

Event yields and associated uncertainties after the global likelihood fit of the differential regions. The total background yield can differ from the sum of the individual components due to rounding. The five STXS volumes are labeled pT0–pT4 corresponding to pTH < 300 GeV, 300 – 450 GeV, 450 – 650 GeV, 650 – 1000 GeV, and > 1000 GeV, respectively. The pT0 event yield is constrained to its SM value within the theoretical and experimental uncertainties and free parameters act independently on the remaining four volumes.

The JER for R = 0.2 jets in Pb+Pb collisions as a function of $2|\Psi_2-\phi^{reco}|$ for centrality selections of 0-5%, 5-10%, 10-20%, 20-40% and 40-60%.

The systematic uncertainties in v2 for 20-40% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v2 for 5-10% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v3 for 20-40% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v3 for 5-10% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v4 for 20-40% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v4 for 5-10% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v2 for pT = 71--398 GeV jets as a function of centrality.

The systematic uncertainties in v3 for pT = 71--398 GeV jets as a function of centrality.

The systematic uncertainties in v4 for pT = 71--398 GeV jets as a function of centrality.

Angular distribution of jets as a function of the observed psi2 plane for jets with 71 < pT < 79 GeV in the 10-20% centrality bin.

Angular distribution of jets as a function of the observed psi3 plane for jets with 71 < pT < 79 GeV in the 10-20% centrality bin.

Angular distribution of jets as a function of the observed psi4 plane for jets with 71 < pT < 79 GeV in the 10-20% centrality bin.

The v2 values for R = 0.2 jets as a function of centrality for jets in several pT ranges.

The v2 values for R = 0.2 jets as a function of pT for 0-5%, 5-10%, and 20-40% centrality collisions.

The v2, v3, and v4 as a function of centrality for jets with pT = 71-398 GeV.

The v3 values for R = 0.2 jets as a function of centrality for jets in several pT ranges.

The v4 values for R = 0.2 jets as a function of centrality for jets in several pT ranges.

The v2 as a function of pT for jets in 10-20% centrality collisions.

The v3 as a function of pT for jets in 10-20% centrality collisions.

The v2 as a function of pT for jets in 20-40% centrality collisions.

The v3 as a function of pT for jets in 20-40% centrality collisions.

The v2 for jets in 10-20% centrality collisions.

The v3 for jets in 10-20% centrality collisions.

R2max as a function of pT (filled circles). Also shown is 1 - 4v2/(1+2v2).

R3max as a function of pT (filled circles). Also shown is 1 - 4v3/(1+2v3).

The systematic uncertainties in v2 for 40-60% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v3 for 40-60% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v4 for 40-60% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v2 for 10-20% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v3 for 10-20% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v4 for 10-20% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v2 for 0-5% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v3 for 0-5% centrality Pb+Pb collisions as a function of pT.

The systematic uncertainties in v4 for 0-5% centrality Pb+Pb collisions as a function of pT.

The v2 as a function of pT for jets in 40-60% centrality collisions.

The v2 as a function of pT for jets in 5-10% centrality collisions.

The v2 as a function of pT for jets in 0-5% centrality collisions.

The v3 as a function of pT for jets in 40-60% centrality collisions.

The v3 as a function of pT for jets in 5-10% centrality collisions.

The v3 as a function of pT for jets in 0-5% centrality collisions.

Data from Figure 11a

Model-independent per-channel efficiencies calculated in the fiducial volume of the High Mass channel for the $2e2\mu$ final state.

Model-independent per-channel efficiencies calculated in the fiducial volume of the High Mass channel for the $4e$ final state.

Model-dependent per-channel fiducial region acceptances for the $H\to aa\to 4\mu$ process. The regions of 2 GeV to 4.4 GeV, and 8 GeV to 12 GeV are each excluded by the quarkonia veto. Both, the Low Mass and the High Mass channel, include $m_X = $ 15 GeV at their signal region's edges, hence the dual entry for $m_X = $ 15 GeV.

Model-dependent per-channel fiducial region acceptances for the $H\to Z_dZ_d\to 4\ell$ process. The regions of 2 GeV to 4.4 GeV, and 8 GeV to 12 GeV are each excluded by the quarkonia veto. Both, the Low Mass and the High Mass channel, include $m_X = $ 15 GeV at their signal region's edges, hence the dual entry for $m_X = $ 15 GeV.

Model-dependent per-channel fiducial region acceptances for the $H\to Z_dZ_d\to 2e2\mu$ process.

Model-dependent per-channel fiducial region acceptances for the $H\to Z_dZ_d\to 4e$ process.

Upper limits at 95% CL on fiducial cross-sections for the $H\rightarrow XX \rightarrow 4\mu$ process. The step change at $m_X$ = 15 GeV is due to the change in efficiency caused by the change in fiducial phase-space definition. The regions of 2 GeV to 4.4 GeV, and 8 GeV to 12 GeV are each excluded by the quarkonia veto.

Upper limits at 95% CL on fiducial cross-sections for the $H\rightarrow XX \rightarrow 4e$ process.

Upper limits at 95% CL on fiducial cross-sections for the $H\rightarrow XX \rightarrow 2e2\mu$ process.

Observed and expected upper limits at 95% CL for the cross section of the $H\to Z_dZ_d\to 4\ell$ process, assuming SM Higgs boson production via the gluon-gluon fusion process. All final states are combined.

Observed and expected upper limits at 95% CL for the cross sections of the $H\to Z_dZ_d\to 4\mu$ processes, assuming SM Higgs boson production via the gluon-gluon fusion process. The regions of 2 GeV to 4.4 GeV, and 8 GeV to 12 GeV are each excluded by the quarkonia veto. The step changes at $m_{Z_d}= 15 \,\text{GeV}$ are due to the change in selection from the Low Mass to the High Mass analysis.

Observed and expected upper limits at 95% CL for the cross sections of the $H\to aa\to 4\mu$ processes, assuming SM Higgs boson production via the gluon-gluon fusion process. The regions of 2 GeV to 4.4 GeV, and 8 GeV to 12 GeV are each excluded by the quarkonia veto. The step changes at $m_{Z_d}= 15 \,\text{GeV}$ are due to the change in selection from the Low Mass to the High Mass analysis.

Model-independent efficiencies for the $H\to ZX$ process for the $4\mu$ and $2e2\mu$ final states combined calculated in the $ZX$ fiducial volume.

Model-independent efficiencies for the $H\to ZX$ process for the $4\mu$, $2e2\mu$, $2\mu 2e$, and $4e$ final states combined calculated in the $ZX$ fiducial volume.

Model-independent efficiencies for the $H\to ZX$ process for the $2\mu 2e$, and $4e$ final states combined calculated in the $ZX$ fiducial volume.

Model-dependent per-channel fiducial region acceptances for the $H\to ZZ_d\to 4\ell$ process for the $4\mu$, and $2e2\mu$ final states.

Model-dependent per-channel fiducial region acceptances for the $H\to ZZ_d\to 4\ell$ process for the $4\mu$, $2e2\mu$, $2\mu 2e$, and $4e$ final states.

Model-dependent per-channel fiducial region acceptances for the $H\to ZZ_d\to 4\ell$ process for the $2\mu 2e$, and $4e$ final states.

Per-channel upper limit at 95% CL on the fiducial cross-section for the $H\rightarrow ZX \to 2\ell2e$ process.

Per-channel upper limit at 95% CL on the fiducial cross-section for the $H\rightarrow ZX \to 2\ell2\mu$ process.

Observed and expected upper limits at 95% CL for the cross sections of the $H\to ZZ_d\to 4\ell$ process, assuming SM Higgs~boson production via the gluon-gluon fusion process.

Observed and expected upper limits at 95% CL for the cross sections of the $H\to Za\to 2\ell2\mu$ process, assuming SM Higgs~boson production via the gluon-gluon fusion process.

95% CL upper limits on the cross section times the model-dependent branching ratio divided by the SM Higgs boson production cross section for the $H\rightarrow Z_dZ_d$ process for the benchmark HAHM. The regions of 2 GeV to 4.4 GeV, and 8 GeV to 12 GeV are each excluded by the quarkonia veto. The step changes at $m_{Z_d}= 15 \,\text{GeV}$ are due to the change in selection from the Low Mass to the High Mass analysis.

95% CL upper limits on the cross section times the model-dependent branching ratio divided by the SM Higgs boson production cross section for the $H\rightarrow aa$ process for the benchmark 2HDM+S model. The regions of 2 GeV to 4.4 GeV, and 8 GeV to 12 GeV are each excluded by the quarkonia veto. The step changes at $m_{Z_d}= 15 \,\text{GeV}$ are due to the change in selection from the Low Mass to the High Mass analysis.

Upper limit at 95% CL on the effective Higgs mixing parameter $\kappa^\prime=\kappa \, m_H^2 / |m_H^2-m_S^2|$. The regions of 2 GeV to 4.4 GeV, and 8 GeV to 12 GeV are each excluded by the quarkonia veto. The step changes at $m_{Z_d}= 15 \,\text{GeV}$ are due to the change in selection from the Low Mass to the High Mass analysis.

Upper limit at 95% CL on the $Z_d$ kinetic mixing parameter $\epsilon$, assuming the SM Higgs boson production cross section.

Upper limit at 95% CL on the $Z$-$Z_d$ mass mixing parameter $\delta^2 \times \textrm{BR}(Z_d\to \ell\ell)$, assuming the SM Higgs boson production cross section.

Parton dijet $M_{inv}$ vs $A_{LL}$ values with associated uncertainties, for topology C.

Parton dijet $M_{inv}$ vs $A_{LL}$ values with associated uncertainties, for topology D.

The correlation matrix for the point-to-point uncertainties (statistical and systematics) for the inclusive jet measurements.The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (statistical and systematics) for the inclusive jet measurements coupling with the forward-forward dijet measurements (topology A). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (statistical and systematics) for the inclusive jet measurements coupling with the forward-forward dijet measurements (topology B). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (statistical and systematics) for the inclusive jet measurements coupling with the forward-forward dijet measurements (topology C). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (statistical and systematics) for the inclusive jet measurements coupling with the forward-forward dijet measurements (topology D). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) for forward-forward dijet measurements (topology A). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology A) with forward-central dijet measurements (topology B). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology A) with forward-central dijet measurements (topology C). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology A) with forward-central dijet measurements (topology D). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) for forward-forward dijet measurements (topology B). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology B) with forward-central dijet measurements (topology C). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology B) with forward-central dijet measurements (topology D). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) for forward-forward dijet measurements (topology C). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) coupling forward-forward dijet measurements (topology C) with forward-central dijet measurements (topology D). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

The correlation matrix for the point-to-point uncertainties (systematics only) for forward-forward dijet measurements (topology D). The relative luminosity and beam polarization uncertainties are not included because they are the same for all points.

Correlation coefficients $\rho_{i,j}$ for the statistical and systematic uncertainties between the $i$-th and $j$-th bin of the differential $A_E$ measurement as a function of the jet scattering angle $\theta_j$

The effects on the energy asymmetry of $1\sigma$ variations in its influencing nuisance parameters for the three $\theta_j$ bins. These are extracted from the samples of the posterior distribution with $\sigma_i^{(j)} = c_{ij}/\sqrt{c_{jj}}$ being the estimated shift of bin $i$ in conjunction with a shift $\Delta\theta_j$ of nuisance parameter $j$. The data statistical (Data stat.) uncertainty is obtained from running the unfolding with all nuisance parameters being fixed to their post-marginalised values, the MC statistical uncertainty on the response matrix ($t\bar{t}$ response MC stat.) is evaluated using a bootstrapping method from the covariance matrix of the ensemble of repeated unfolding results with varied response matrices. The $\gamma$ variations denote the statistical uncertainties of the background predictions in the corresponding bin of the $\Delta E$ vs $\theta_{j}$ distribution. The numbers appended to the $W$+jets PDF variations denote the corresponding NNPDF3.0 PDF sets.

The effects on the energy asymmetry of $1\sigma$ variations in its influencing nuisance parameters for the three $\theta_j$ bins. These are extracted from the samples of the posterior distribution with $\sigma_i^{(j)} = c_{ij}/\sqrt{c_{jj}}$ being the estimated shift of bin $i$ in conjunction with a shift $\Delta\theta_j$ of nuisance parameter $j$. The data statistical (Data stat.) uncertainty is obtained from running the unfolding with all nuisance parameters being fixed to their post-marginalised values, the MC statistical uncertainty on the response matrix ($t\bar{t}$ response MC stat.) is evaluated using a bootstrapping method from the covariance matrix of the ensemble of repeated unfolding results with varied response matrices. The $\gamma$ variations denote the statistical uncertainties of the background predictions in the corresponding bin of the $\Delta E$ vs $\theta_{j}$ distribution. The numbers appended to the $W$+jets PDF variations denote the corresponding NNPDF3.0 PDF sets.

Covariances $c_{ij}$ for the highest ranked systematic uncertainties in the energy asymmetry measurement differential in $\theta_j$.

Covariances $c_{ij}$ for the highest ranked systematic uncertainties in the energy asymmetry measurement differential in $\theta_j$.

Definition of the fiducial phase space with lepton candidate ($\ell$), hadronic top candidate ($j_h$), leptonic top candidate ($j_l$) and associated jet candidate ($j_a$).

Definition of the fiducial phase space with lepton candidate ($\ell$), hadronic top candidate ($j_h$), leptonic top candidate ($j_l$) and associated jet candidate ($j_a$).

Bounds on individual Wilson coefficients $C($TeV$/\Lambda)^2$ from the energy asymmetry, obtained from a combined fit to its values in all three $\theta_j$ bins.

Bounds on individual Wilson coefficients $C($TeV$/\Lambda)^2$ from the energy asymmetry, obtained from a combined fit to its values in all three $\theta_j$ bins.