Nuclear modification factor of $\Upsilon(1\mathrm{S})$ as a function of the average number of participants $\langle N_{\mathrm{part}} \rangle$ or as a function of the collision centrality.

Nuclear modification factor of $\Upsilon(2\mathrm{S})$ as a function of the average number of participants $\langle N_{\mathrm{part}} \rangle$ or as a function of the collision centrality.

Ratio of $\Upsilon(2\mathrm{S})$ and $\Upsilon(1\mathrm{S})$ yields (cf. equation 2 in the article for the definition) as a function of the average number of participants $\langle N_{\mathrm{part}} \rangle$ or as a function of the collision centrality. The global uncertainty is the quadratic sum of the branching ratio uncertainties.

Relative nuclear modification factor or double yield ratio between $\Upsilon(2\mathrm{S})$ and $\Upsilon(1\mathrm{S})$ as a function of the average number of participants $\langle N_{\mathrm{part}} \rangle$ or as a function of the collision centrality. The global uncertainty corresponds to the systematic uncertainty on the cross-section ratio in proton–proton collisions.

Nuclear modification factor of $\Upsilon(1\mathrm{S})$ as a function of transverse momentum for the 0–90% centrality interval.

Nuclear modification factor of $\Upsilon(1\mathrm{S})$ as a function of rapidity for the 0–90% centrality interval.

Nuclear modification factor of $\Upsilon(2\mathrm{S})$ as a function of rapidity for the 0–90% centrality interval.

Interpolated production cross sections and cross-section ratio for $\Upsilon(1\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ and $\Upsilon(2\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ in proton–proton collisions, integrated over the rapidity interval 2.5–4.0 and $p_{\mathrm{T}} <$ 15 GeV/$c$. Results used for the determination of the integrated nuclear modification factors and of the nuclear modification factors as a function of centrality (Figures 4 and 5 (right)).

Interpolated $p_{\mathrm{T}}$-differential cross section for $\Upsilon(1\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ in proton–proton collisions. Results used for the determination of the nuclear modification factor of $\Upsilon(1\mathrm{S})$ as a function of its transverse momentum (Figure 6).

Interpolated rapidity-differential cross section for $\Upsilon(1\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ in proton–proton collisions. Results used for the determination of the nuclear modification factor of $\Upsilon(1\mathrm{S})$ as a function of rapidity (Figure 7).

Interpolated rapidity-differential cross section for $\Upsilon(2\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ in proton–proton collisions. Results used for the determination of the nuclear modification factor of $\Upsilon(2\mathrm{S})$ as a function of rapidity (Figure 7).

Distribution of $E^{miss}_T$ in data and for the expected SM background in the Single-Muon CR after performing the 'simplified shape fit'. Overflows are included in the fourth bin of each distribution. The $E^{miss}_T$ calculation in this CR does not include the muon contribution. The error bars are statistical, and the dashed band includes statistical and systematic uncertainties determined by the fit. The lower panel shows the ratio of data to expected background event yields.

Distribution of $E^{miss}_T$ in data and for the expected SM background in the Two-Muon CR after performing the 'simplified shape fit'. The $E^{miss}_T$ calculation in this CR does not include the muon contribution. The error bars are statistical, and the dashed band includes statistical and systematic uncertainties determined by the fit. The lower panel shows the ratio of data to expected background event yields.

Distribution of $E^{miss}_T$ in data and for the expected SM background in the Two-Muon CR after performing the 'simplified shape fit'. Overflows are included in the fourth bin of each distribution. The $E^{miss}_T$ calculation in this CR does not include the muon contribution. The error bars are statistical, and the dashed band includes statistical and systematic uncertainties determined by the fit. The lower panel shows the ratio of data to expected background event yields.

Distribution of $E^{miss}_T$ in data and for the expected SM background in the Two-Electron CR after performing the 'simplified shape fit'. The $E^{miss}_$ calculation in this CR does not include electron contribution. The error bars are statistical, and the dashed band includes statistical and systematic uncertainties determined by the fit. The lower panel shows the ratio of data to expected background event yields.

Distribution of $E^{miss}_T$ in data and for the expected SM background in the Two-Electron CR after performing the 'simplified shape fit'. Overflows are included in the fourth bin of each distribution. The $E^{miss}_T$ calculation in this CR does not include electron contribution. The error bars are statistical, and the dashed band includes statistical and systematic uncertainties determined by the fit. The lower panel shows the ratio of data to expected background event yields.

The observed exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The observed exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The expected exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The expected exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The expected + 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The expected + 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The expected - 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The expected - 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The observed exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.1.

The observed exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.1.

The expected exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.1.

The expected exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.1.

The expected + 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.1.

The expected + 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.1.

The expected - 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.1.

The expected - 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a simplified DM model involving an axial-vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.1.

The observed exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The observed exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The expected exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The expected exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The expected + 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The expected + 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The expected - 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The expected - 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0.

The observed exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings and $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.01

The observed exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings and $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.01

The expected exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.01.

The expected exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.01.

The expected + 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.01.

The expected + 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.01.

The expected - 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.01.

The expected - 1 sigma exclusion contours on the $m_{\chi} - m_{med}$ plane for a vector mediator with couplings $g_{\chi}$ = 1, $g_q$ = 0.1 and $g_l$ = 0.01.

The 90% CL exclusion limit on $\chi$-proton spin-dependent scattering cross section in an axial-vector model with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 as a function of the dark-matter mass $m_{\chi}$.

The 90% CL exclusion limit on $\chi$-proton spin-dependent scattering cross section in an axial-vector model with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 as a function of the dark-matter mass $m_{\chi}$. The 90% CL exclusion limit is displayed only in the regime where the theory is perturbative.

The 90% CL exclusion limit on $\chi$-neutron spin-dependent scattering cross section in an axial-vector model with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 as a function of the dark-matter mass $m_{\chi}$.

The 90% CL exclusion limit on $\chi$-neutron spin-dependent scattering cross section in an axial-vector model with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 as a function of the dark-matter mass $m_{\chi}$. The 90% CL exclusion limit is displayed only in the regime where the theory is perturbative.

The 90% CL exclusion limit on $\chi$-nucleon spin-independent scattering cross section in a vector model with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 as a function of the dark-matter mass $m_{\chi}$.

The 90% CL exclusion limit on $\chi$-nucleon spin-independent scattering cross section in a vector model with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 as a function of the dark-matter mass $m_{\chi}$.

Observed and expected exclusion at 95% CL on the coupling $c_W$ as a function of the effective scale $f_a$ for an ALP mass of 1 MeV.

Observed and expected exclusion at 95% CL on the coupling $c_W$ as a function of the effective scale $f_a$ for an ALP mass of 1 MeV.

The observed 95% CL upper limits on signal strength $\mu$, the ratio of the experimental limit to the predicted signal cross section, for a simplified model of dark matter production involving an axial-vector operator and couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 as a function of the mass of the dark matter, $m_{\chi}$ and the mass of the axial-vector mediator, $m_{med}$.

The observed 95% CL upper limits on signal strength $\mu$, the ratio of the experimental limit to the predicted signal cross section, for a simplified model of dark matter production involving an axial-vector operator and couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 as a function of the mass of the dark matter, $m_{\chi}$ and the mass of the axial-vector mediator, $m_{med}$.

The observed 95% CL upper limits on cross section for a simplified model of dark matter production involving an axial-vector operator and couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 as a function of the mass of the dark matter, $m_{\chi}$ and the mass of the axial-vector mediator, $m_{med}$.

The observed 95% CL upper limits on cross section for a simplified model of dark matter production involving an axial-vector operator and couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 as a function of the mass of the dark matter, $m_{\chi}$ and the mass of the axial-vector mediator, $m_{med}$.

Acceptance in % for a simplified model of dark matter production involving an axial-vector operator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 for different values of $m_{\chi}$ and $m_{med}$ in the four inclusive signal regions.

Acceptance in % for a simplified model of dark matter production involving an axial-vector operator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 for different values of $m_{\chi}$ and $m_{med}$ in the four inclusive signal regions.

Acceptance in % for a simplified model of dark matter production involving an axial-vector operator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 for different values of $m_{\chi}$ and $m_{med}$ in the three exclusive signal regions.

Acceptance in % for a simplified model of dark matter production involving an axial-vector operator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 for different values of $m_{\chi}$ and $m_{med}$ in the three exclusive signal regions.

Efficiencies for a simplified model of dark matter production involving an axial-vector operator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 for different values of $m_{\chi}$ and $m_{med}$ in the four inclusive signal regions.

Efficiencies for a simplified model of dark matter production involving an axial-vector operator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 for different values of $m_{\chi}$ and $m_{med}$ in the four inclusive signal regions.

Efficiencies for a simplified model of dark matter production involving an axial-vector operator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 for different values of $m_{\chi}$ and $m_{med}$ in the three exclusive signal regions.

Efficiencies for a simplified model of dark matter production involving an axial-vector operator with couplings $g_{\chi}$ = 1, $g_q$ = 0.25 and $g_l$ = 0 for different values of $m_{\chi}$ and $m_{med}$ in the three exclusive signal regions.

Neutron multiplicity dependence of acoplanarity ($\alpha$) from process $\gamma\gamma$ to $\mu^+\mu^-$ in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV.

Neutron multiplicity dependence of acoplanarity ($\alpha$) from process $\gamma\gamma$ to $\mu^+\mu^-$ in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV.

Neutron multiplicity dependence of acoplanarity ($\alpha$) from process $\gamma\gamma$ to $\mu^+\mu^-$ in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV.

Neutron multiplicity dependence of average core acoplanarity ($<\alpha^{core}>$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV.

Neutron multiplicity dependence of average invariant mass ($<m_{\mu\mu}>$) from process $\gamma\gamma$ to $\mu^+\mu^-$ in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 0n1n class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing larger neutron multiplicity.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 0nXn class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing larger neutron multiplicity.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 1nXn class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing larger neutron multiplicity.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 0n1n class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing smaller neutron multiplicity.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 0nXn class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing smaller neutron multiplicity.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 1nXn class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing smaller neutron multiplicity.

The measured total fiducial cross sections. The values labeled 'stat.' and 'eff.' represent the statistical uncertainty and the systematic uncertainty estimated from the efficiencies, respectively. The value 'sys.' includes all remaining systematic uncertainties, with the exception of the luminosity. The 9\% uncertainty associated with the luminosity measurement is labeled as 'lumi'.

The measured total cross sections. The values labeled 'stat.' and 'eff.' represent the statistical uncertainty and the systematic uncertainty estimated from the efficiencies, respectively. The value 'sys.' includes all remaining systematic uncertainties, with the exception of the luminosity. The 9\% uncertainty associated with the luminosity measurement is labeled as 'lumi'.

Measurements for the ratio of the fiducial cross sections for production of $W^+$ and $W^-$ bosons in bins of the decay charged lepton pseudorapidity.

The measured $(W^+ + W^-)/Z$

Summary of the sources of systematic and statistical uncertainties and their impact on the measurement of the ttH and tH signal rates, and the measured value of the unconstrained nuisance parameters. The quantity $\Delta r_{x}/r_{x}$ corresponds to the change in uncertainty when fixing the nuisance parameters associated with that uncertainty in the fit. Under the label "MC and sideband statistical uncertainty" are the uncertainties associated with the limited number of simulated MC events and the amount of data events in the application region of the MP method.

Number of events selected in the 2$\ell$ss+0$\tau_{\mathrm{h}}$, 3$\ell$+0$\tau_{\mathrm{h}}$, and 2$\ell$ss+1$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Number of events selected in the 2$\ell$ss+0$\tau_{\mathrm{h}}$, 3$\ell$+0$\tau_{\mathrm{h}}$, and 2$\ell$ss+1$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Number of events selected in the 1$\ell$+1$\tau_{\mathrm{h}}$, 0$\ell$+2$\tau_{\mathrm{h}}$, 2$\ell$os+1$\tau_{\mathrm{h}}$, and 1$\ell$+2$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Number of events selected in the 1$\ell$+1$\tau_{\mathrm{h}}$, 0$\ell$+2$\tau_{\mathrm{h}}$, 2$\ell$os+1$\tau_{\mathrm{h}}$, and 1$\ell$+2$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Number of events selected in the 4$\ell$+0$\tau_{\mathrm{h}}$, 3$\ell$+1$\tau_{\mathrm{h}}$, and 2$\ell$+2$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Number of events selected in the 4$\ell$+0$\tau_{\mathrm{h}}$, 3$\ell$+1$\tau_{\mathrm{h}}$, and 2$\ell$+2$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Dependence of the likelihood function L as a function of $\kappa_{\mathrm{t}}$ and $\kappa_{\mathrm{V}}$.

Dependence of the likelihood function L as a function of $\kappa_{\mathrm{t}}$ and $\kappa_{\mathrm{V}}$.

$-2\Delta log( L)$ as a function of the production rate of ttH and tH.

$-2\Delta log( L)$ as a function of the production rate of ttH and tH.

$-2\Delta log( L)$ as a function of the production rate of ttH and ttW.

$-2\Delta log( L)$ as a function of the production rate of ttH and ttW.

$-2\Delta log( L)$ as a function of the production rate of ttH and ttZ.

$-2\Delta log( L)$ as a function of the production rate of ttH and ttZ.

$-2\Delta log( L)$ as a function of the production rate of ttW and ttZ.

$-2\Delta log( L)$ as a function of the production rate of ttW and ttZ.

Covariance matrix of the fitted parameters of interest.

Covariance matrix of the fitted parameters of interest.

Production rate of the ttH signal measured in each of the ten channels individually and for the combination of all channels

Production rate of the ttH signal measured in each of the ten channels individually and for the combination of all channels

Production rate of the tH signal measured in each of the ten channels individually and for the combination of all channels

Production rate of the tH signal measured in each of the ten channels individually and for the combination of all channels

Probability for ttH signal events to pass the event selection in each of the analysis categories employed.

Observed metSig distributions in signal regions MB-SSd. The histograms show the MC background predictions normalised by the background-only fit. The hatched (red) error bands indicate experimental and MC statistical uncertainties. Expected distributions for benchmark signal model points, normalised using the approximate NNLO+NNLL cross-section times integrated luminosity, are also shown for comparison (masses in GeV).

Observed $m_{\mathrm{eff}}$ distributions in signal regions MB-GGd. The histograms show the MC background predictions normalised by the background-only fit. The hatched (red) error bands indicate experimental and MC statistical uncertainties. Expected distributions for benchmark signal model points, normalised using the approximate NNLO+NNLL cross-section times integrated luminosity, are also shown for comparison (masses in GeV).

Observed $m_{\mathrm{eff}}$ distributions in signal regions MB-GGd. The histograms show the MC background predictions normalised by the background-only fit. The hatched (red) error bands indicate experimental and MC statistical uncertainties. Expected distributions for benchmark signal model points, normalised using the approximate NNLO+NNLL cross-section times integrated luminosity, are also shown for comparison (masses in GeV).

Observed metSig distributions in signal regions MB-GGd. The histograms show the MC background predictions normalised by the background-only fit. The hatched (red) error bands indicate experimental and MC statistical uncertainties. Expected distributions for benchmark signal model points, normalised using the approximate NNLO+NNLL cross-section times integrated luminosity, are also shown for comparison (masses in GeV).

Observed metSig distributions in signal regions MB-GGd. The histograms show the MC background predictions normalised by the background-only fit. The hatched (red) error bands indicate experimental and MC statistical uncertainties. Expected distributions for benchmark signal model points, normalised using the approximate NNLO+NNLL cross-section times integrated luminosity, are also shown for comparison (masses in GeV).

Observed $m_{\mathrm{eff}}$ distributions in signal regions MB-C. The histograms show the MC background predictions normalised by the background-only fit. The hatched (red) error bands indicate experimental and MC statistical uncertainties. Expected distributions for benchmark signal model points, normalised using the approximate NNLO+NNLL cross-section times integrated luminosity, are also shown for comparison (masses in GeV).

Observed $m_{\mathrm{eff}}$ distributions in signal regions MB-C. The histograms show the MC background predictions normalised by the background-only fit. The hatched (red) error bands indicate experimental and MC statistical uncertainties. Expected distributions for benchmark signal model points, normalised using the approximate NNLO+NNLL cross-section times integrated luminosity, are also shown for comparison (masses in GeV).

Observed metSig distributions in signal regions MB-C. The histograms show the MC background predictions normalised by the background-only fit. The hatched (red) error bands indicate experimental and MC statistical uncertainties. Expected distributions for benchmark signal model points, normalised using the approximate NNLO+NNLL cross-section times integrated luminosity, are also shown for comparison (masses in GeV).

Observed metSig distributions in signal regions MB-C. The histograms show the MC background predictions normalised by the background-only fit. The hatched (red) error bands indicate experimental and MC statistical uncertainties. Expected distributions for benchmark signal model points, normalised using the approximate NNLO+NNLL cross-section times integrated luminosity, are also shown for comparison (masses in GeV).

Observed BDT-GGd1 score distributions in signal regions GGd1. The histograms show the MC background predictions normalised by the background-only fit. The hatched (red) error bands indicate experimental and MC statistical uncertainties. Expected distributions for benchmark signal model points, normalised using the approximate NNLO+NNLL cross-section times integrated luminosity, are also shown for comparison (masses in GeV).

Observed BDT-GGd1 score distributions in signal regions GGd1. The histograms show the MC background predictions normalised by the background-only fit. The hatched (red) error bands indicate experimental and MC statistical uncertainties. Expected distributions for benchmark signal model points, normalised using the approximate NNLO+NNLL cross-section times integrated luminosity, are also shown for comparison (masses in GeV).

Observed BDT-GGo1 score distributions in signal regions GGo1. The histograms show the MC background predictions normalised by the background-only fit. The hatched (red) error bands indicate experimental and MC statistical uncertainties. Expected distributions for benchmark signal model points, normalised using the approximate NNLO+NNLL cross-section times integrated luminosity, are also shown for comparison (masses in GeV).

Observed BDT-GGo1 score distributions in signal regions GGo1. The histograms show the MC background predictions normalised by the background-only fit. The hatched (red) error bands indicate experimental and MC statistical uncertainties. Expected distributions for benchmark signal model points, normalised using the approximate NNLO+NNLL cross-section times integrated luminosity, are also shown for comparison (masses in GeV).

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and neutralino in SR BDT-GGd1

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and neutralino in SR BDT-GGd1

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and neutralino in SR BDT-GGd2

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and neutralino in SR BDT-GGd2

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and neutralino in SR BDT-GGd3

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and neutralino in SR BDT-GGd3

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and neutralino in SR BDT-GGd4

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and neutralino in SR BDT-GGd4

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and chargino in SR BDT-GGo1

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and chargino in SR BDT-GGo1

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and chargino in SR BDT-GGo2

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and chargino in SR BDT-GGo2

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and chargino in SR BDT-GGo3

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and chargino in SR BDT-GGo3

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and chargino in SR BDT-GGo4

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and chargino in SR BDT-GGo4

Signal region acceptance for simplified model with squark pair production and decays to a quark and neutralino in SR 2j-1600

Signal region acceptance for simplified model with squark pair production and decays to a quark and neutralino in SR 2j-1600

Signal region acceptance for simplified model with squark pair production and decays to a quark and neutralino in SR 2j-2200

Signal region acceptance for simplified model with squark pair production and decays to a quark and neutralino in SR 2j-2200

Signal region acceptance for simplified model with squark pair production and decays to a quark and neutralino in SR 2j-2800

Signal region acceptance for simplified model with squark pair production and decays to a quark and neutralino in SR 2j-2800

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and neutralino in SR 4j-1000

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and neutralino in SR 4j-1000

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and neutralino in SR 4j-2200

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and neutralino in SR 4j-2200

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and neutralino in SR 4j-3400

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and neutralino in SR 4j-3400

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and neutralino in SR 5j-1600

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and neutralino in SR 5j-1600

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and chargino in SR 6j-1000

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and chargino in SR 6j-1000

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and chargino in SR 6j-2200

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and chargino in SR 6j-2200

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and chargino in SR 6j-3400

Signal region acceptance for simplified model with gluino pair production and decays to two quarks and chargino in SR 6j-3400

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and neutralino in SR BDT-GGd1. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and neutralino in SR BDT-GGd1. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and neutralino in SR BDT-GGd2. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and neutralino in SR BDT-GGd2. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and neutralino in SR BDT-GGd3. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and neutralino in SR BDT-GGd3. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and neutralino in SR BDT-GGd4. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and neutralino in SR BDT-GGd4. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and chargino in SR BDT-GGo1. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and chargino in SR BDT-GGo1. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and chargino in SR BDT-GGo2. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and chargino in SR BDT-GGo2. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and chargino in SR BDT-GGo3. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and chargino in SR BDT-GGo3. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and chargino in SR BDT-GGo4. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and chargino in SR BDT-GGo4. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with squark pair production and decays to a quark and neutralino in SR 2j-1600. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with squark pair production and decays to a quark and neutralino in SR 2j-1600. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with squark pair production and decays to a quark and neutralino in SR 2j-2200. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with squark pair production and decays to a quark and neutralino in SR 2j-2200. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with squark pair production and decays to a quark and neutralino in SR 2j-2800. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with squark pair production and decays to a quark and neutralino in SR 2j-2800. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and neutralino in SR 4j-1000. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and neutralino in SR 4j-1000. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and neutralino in SR 4j-2200. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and neutralino in SR 4j-2200. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and neutralino in SR 4j-3400. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and neutralino in SR 4j-3400. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and neutralino in SR 5j-1600. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and neutralino in SR 5j-1600. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and chargino in SR 6j-1000. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and chargino in SR 6j-1000. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and chargino in SR 6j-2200. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and chargino in SR 6j-2200. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and chargino in SR 6j-3400. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Signal region efficiency for simplified model with gluino pair production and decays to two quarks and chargino in SR 6j-3400. Efficiencies on signal points with low statistics are not reported. The efficiency value -1.0 in the table corresponds to the case where efficiency cannot be calculated due to the null acceptance of the model point.

Exclusion limits in the mass plane of the lightest neutralino and first- and second-generation squarks assuming squark pair production and direct decays obtained by using the signal region with the best expected sensitivity at each point. The expected limits are indicated with a dark dashed curve, with the light (yellow) band indicating the $1\sigma$ excursions due to experimental and background-only theoretical uncertainties.

Exclusion limits in the mass plane of the lightest neutralino and first- and second-generation squarks assuming squark pair production and direct decays obtained by using the signal region with the best expected sensitivity at each point. The expected limits are indicated with a dark dashed curve, with the light (yellow) band indicating the $1\sigma$ excursions due to experimental and background-only theoretical uncertainties.

Exclusion limits in the mass plane of the lightest neutralino and first- and second-generation squarks assuming squark pair production and direct decays obtained by using the signal region with the best expected sensitivity at each point. The observed limits are indicated by the medium dark (maroon) curve where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalisation and factorisation scale and PDF uncertainties.

Exclusion limits in the mass plane of the lightest neutralino and first- and second-generation squarks assuming squark pair production and direct decays obtained by using the signal region with the best expected sensitivity at each point. The observed limits are indicated by the medium dark (maroon) curve where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalisation and factorisation scale and PDF uncertainties.

Exclusion limits in the mass plane of the lightest neutralino and non degenerated squark pair production and direct decays obtained by using the signal region with the best expected sensitivity at each point. The expected limits are indicated with a dark dashed curve, with the light (yellow) band indicating the $1\sigma$ excursions due to experimental and background-only theoretical uncertainties.

Exclusion limits in the mass plane of the lightest neutralino and non degenerated squark pair production and direct decays obtained by using the signal region with the best expected sensitivity at each point. The expected limits are indicated with a dark dashed curve, with the light (yellow) band indicating the $1\sigma$ excursions due to experimental and background-only theoretical uncertainties.

Exclusion limits in the mass plane of the lightest neutralino and non degenerated squark pair production and direct decays obtained by using the signal region with the best expected sensitivity at each point. The observed limits are indicated by the medium dark (maroon) curve where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalisation and factorisation scale and PDF uncertainties.

Exclusion limits in the mass plane of the lightest neutralino and non degenerated squark pair production and direct decays obtained by using the signal region with the best expected sensitivity at each point. The observed limits are indicated by the medium dark (maroon) curve where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalisation and factorisation scale and PDF uncertainties.

Exclusion limits in the mass plane of the lightest neutralino and the gluino for gluino pair production with direct decays obtained by using the signal region with the best expected sensitivity at each point. The expected limits are indicated with a dark dashed curve, with the light (yellow) band indicating the $1\sigma$ excursions due to experimental and background-only theoretical uncertainties.

Exclusion limits in the mass plane of the lightest neutralino and the gluino for gluino pair production with direct decays obtained by using the signal region with the best expected sensitivity at each point. The expected limits are indicated with a dark dashed curve, with the light (yellow) band indicating the $1\sigma$ excursions due to experimental and background-only theoretical uncertainties.

Exclusion limits in the mass plane of the lightest neutralino and the gluino for gluino pair production with direct decays obtained by using the signal region with the best expected sensitivity at each point. The observed limits are indicated by the medium dark (maroon) curve where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalisation and factorisation scale and PDF uncertainties.

Exclusion limits in the mass plane of the lightest neutralino and the gluino for gluino pair production with direct decays obtained by using the signal region with the best expected sensitivity at each point. The observed limits are indicated by the medium dark (maroon) curve where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalisation and factorisation scale and PDF uncertainties.

Exclusion limits in the mass plane of the lightest neutralino and squarks. The expected limits are indicated with dark dashed curves, with the light (yellow) bands indicating the $1\sigma$ excursions due to experimental and background-only theoretical uncertainties.

Exclusion limits in the mass plane of the lightest neutralino and squarks. The expected limits are indicated with dark dashed curves, with the light (yellow) bands indicating the $1\sigma$ excursions due to experimental and background-only theoretical uncertainties.

Exclusion limits in the mass plane of the lightest neutralino and squarks. The observed limits are indicated by the medium dark (maroon) curves where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalisation and factorisation scale and PDF uncertainties.

Exclusion limits in the mass plane of the lightest neutralino and squarks. The observed limits are indicated by the medium dark (maroon) curves where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalisation and factorisation scale and PDF uncertainties.

Exclusion limits for squark pair production with a one-step decay via an intermediate chargino. The neutralino mass is fixed at 60 GeV and exclusion limits are given for mass difference ratio, $X$, as a function of the squark mass. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The expected limits are indicated with dark dashed curves, with the light (yellow) bands indicating the $1\sigma$ excursions due to experimental and background-only theoretical uncertainties.

Exclusion limits for squark pair production with a one-step decay via an intermediate chargino. The neutralino mass is fixed at 60 GeV and exclusion limits are given for mass difference ratio, $X$, as a function of the squark mass. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The expected limits are indicated with dark dashed curves, with the light (yellow) bands indicating the $1\sigma$ excursions due to experimental and background-only theoretical uncertainties.

Exclusion limits for squark pair production with a one-step decay via an intermediate chargino. The neutralino mass is fixed at 60 GeV and exclusion limits are given for mass difference ratio, $X$, as a function of the squark mass. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The observed limits are indicated by the medium dark (maroon) curves where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalisation and factorisation scale and PDF uncertainties.

Exclusion limits for squark pair production with a one-step decay via an intermediate chargino. The neutralino mass is fixed at 60 GeV and exclusion limits are given for mass difference ratio, $X$, as a function of the squark mass. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The observed limits are indicated by the medium dark (maroon) curves where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalisation and factorisation scale and PDF uncertainties.

Exclusion limits in the mass plane of the lightest neutralino and gluinos. The expected limits are indicated with dark dashed curves, with the light (yellow) bands indicating the $1\sigma$ excursions due to experimental and background-only theoretical uncertainties.

Exclusion limits in the mass plane of the lightest neutralino and gluinos. The expected limits are indicated with dark dashed curves, with the light (yellow) bands indicating the $1\sigma$ excursions due to experimental and background-only theoretical uncertainties.

Exclusion limits in the mass plane of the lightest neutralino and gluinos. The observed limits are indicated by the medium dark (maroon) curve where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalisation and factorisation scale and PDF uncertainties.

Exclusion limits in the mass plane of the lightest neutralino and gluinos. The observed limits are indicated by the medium dark (maroon) curve where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalisation and factorisation scale and PDF uncertainties.

Exclusion limits for gluino pair production with a one-step decay via an intermediate chargino. The neutralino mass is fixed at 60 GeV and exclusion limits are given for mass difference ratio, $X$, as a function of the gluino mass. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The expected limits are indicated with dark dashed curves, with the light (yellow) bands indicating the $1\sigma$ excursions due to experimental and background-only theoretical uncertainties.

Exclusion limits for gluino pair production with a one-step decay via an intermediate chargino. The neutralino mass is fixed at 60 GeV and exclusion limits are given for mass difference ratio, $X$, as a function of the gluino mass. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The expected limits are indicated with dark dashed curves, with the light (yellow) bands indicating the $1\sigma$ excursions due to experimental and background-only theoretical uncertainties.

Exclusion limits for gluino pair production with a one-step decay via an intermediate chargino. The neutralino mass is fixed at 60 GeV and exclusion limits are given for mass difference ratio, $X$, as a function of the gluino mass. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The observed limits are indicated by the medium dark (maroon) curves where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalisation and factorisation scale and PDF uncertainties.

Exclusion limits for gluino pair production with a one-step decay via an intermediate chargino. The neutralino mass is fixed at 60 GeV and exclusion limits are given for mass difference ratio, $X$, as a function of the gluino mass. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The observed limits are indicated by the medium dark (maroon) curves where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalisation and factorisation scale and PDF uncertainties.

Exclusion limits for the model with combined production of squark pairs, gluino pairs, and of squark--gluino pairs. The neutralino mass is fixed at 0 GeV. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The expected limits are indicated with dark dashed curves, with the light (yellow) bands indicating the $1\sigma$ excursions due to experimental and background-only theoretical uncertainties.

Exclusion limits for the model with combined production of squark pairs, gluino pairs, and of squark--gluino pairs. The neutralino mass is fixed at 0 GeV. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The expected limits are indicated with dark dashed curves, with the light (yellow) bands indicating the $1\sigma$ excursions due to experimental and background-only theoretical uncertainties.

Exclusion limits for the model with combined production of squark pairs, gluino pairs, and of squark--gluino pairs. The neutralino mass is fixed at 0 GeV. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The observed limits are indicated by the medium dark (maroon) curves where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalisation and factorisation scale and PDF uncertainties.

Exclusion limits for the model with combined production of squark pairs, gluino pairs, and of squark--gluino pairs. The neutralino mass is fixed at 0 GeV. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The observed limits are indicated by the medium dark (maroon) curves where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalisation and factorisation scale and PDF uncertainties.

Exclusion limits for the model with combined production of squark pairs, gluino pairs, and of squark--gluino pairs. The neutralino mass is fixed at 995 GeV. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The expected limits are indicated with dark dashed curves, with the light (yellow) bands indicating the $1\sigma$ excursions due to experimental and background-only theoretical uncertainties.

Exclusion limits for the model with combined production of squark pairs, gluino pairs, and of squark--gluino pairs. The neutralino mass is fixed at 995 GeV. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The expected limits are indicated with dark dashed curves, with the light (yellow) bands indicating the $1\sigma$ excursions due to experimental and background-only theoretical uncertainties.

Exclusion limits for the model with combined production of squark pairs, gluino pairs, and of squark--gluino pairs. The neutralino mass is fixed at 995 GeV. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The observed limits are indicated by the medium dark (maroon) curves where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalisation and factorisation scale and PDF uncertainties.

Exclusion limits for the model with combined production of squark pairs, gluino pairs, and of squark--gluino pairs. The neutralino mass is fixed at 995 GeV. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The observed limits are indicated by the medium dark (maroon) curves where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalisation and factorisation scale and PDF uncertainties.

Exclusion limits for the model with combined production of squark pairs, gluino pairs, and of squark--gluino pairs. The neutralino mass is fixed at 1495 GeV. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The expected limits are indicated with dark dashed curves, with the light (yellow) bands indicating the $1\sigma$ excursions due to experimental and background-only theoretical uncertainties.

Exclusion limits for the model with combined production of squark pairs, gluino pairs, and of squark--gluino pairs. The neutralino mass is fixed at 1495 GeV. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The expected limits are indicated with dark dashed curves, with the light (yellow) bands indicating the $1\sigma$ excursions due to experimental and background-only theoretical uncertainties.

Exclusion limits for the model with combined production of squark pairs, gluino pairs, and of squark--gluino pairs. The neutralino mass is fixed at 1495 GeV. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The observed limits are indicated by the medium dark (maroon) curves where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalisation and factorisation scale and PDF uncertainties.

Exclusion limits for the model with combined production of squark pairs, gluino pairs, and of squark--gluino pairs. The neutralino mass is fixed at 1495 GeV. Exclusion limits are obtained by using the signal region with the best expected sensitivity at each point. The observed limits are indicated by the medium dark (maroon) curves where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the signal cross-section by the renormalisation and factorisation scale and PDF uncertainties.

The observed upper limits on signal cross section corresponding to the best expected signal region in each mass point for squark-pair production with direct decays.

The observed upper limits on signal cross section corresponding to the best expected signal region in each mass point for squark-pair production with direct decays.

The observed upper limits on signal cross section corresponding to the best expected signal region in each mass point for gluino-pair production with direct decays

The observed upper limits on signal cross section corresponding to the best expected signal region in each mass point for gluino-pair production with direct decays

The observed upper limits on signal cross section corresponding to the best expected signal region in each mass point for squark-pair production with a one-step decay via an intermediate chargino.

The observed upper limits on signal cross section corresponding to the best expected signal region in each mass point for squark-pair production with a one-step decay via an intermediate chargino.

The observed upper limits on signal cross section corresponding to the best expected signal region in each mass point for squark-pair production with a one-step decay via an intermediate chargino. The neutralino mass is fixed at 60 GeV and exclusion limits are given for mass difference ratio, $X$, as a function of the squark mass.

The observed upper limits on signal cross section corresponding to the best expected signal region in each mass point for squark-pair production with a one-step decay via an intermediate chargino. The neutralino mass is fixed at 60 GeV and exclusion limits are given for mass difference ratio, $X$, as a function of the squark mass.

The observed upper limits on signal cross section corresponding to the best expected signal region in each mass point for gluino-pair production with a one-step decay via an intermediate chargino.

The observed upper limits on signal cross section corresponding to the best expected signal region in each mass point for gluino-pair production with a one-step decay via an intermediate chargino.

The observed upper limits on signal cross section corresponding to the best expected signal region in each mass point for gluino-pair production with a one-step decay via an intermediate chargino. The neutralino mass is fixed at 60~GeV and exclusio limits are given for mass difference ratio, $X$, as a function of the gluino mass.

The observed upper limits on signal cross section corresponding to the best expected signal region in each mass point for gluino-pair production with a one-step decay via an intermediate chargino. The neutralino mass is fixed at 60~GeV and exclusio limits are given for mass difference ratio, $X$, as a function of the gluino mass.

Cut-flow for model-independent search regions targeting squarks for SS direct model points. Expected yields are normalized to a luminosity of 139 fb$^{-1}$.

Cut-flow for model-independent search regions targeting squarks for SS direct model points. Expected yields are normalized to a luminosity of 139 fb$^{-1}$.

Cut-flow for model-independent search regions targeting gluinos for GG direct model points. Expected yields are normalized to a luminosity of 139 fb$^{-1}$.

Cut-flow for model-independent search regions targeting gluinos for GG direct model points. Expected yields are normalized to a luminosity of 139 fb$^{-1}$.

Cut-flow for model-independent search regions targeting squarks and gluinos in models with one-step decay. Expected yields are normalized to a luminosity of 139 $fb^{-1}$.

Cut-flow for model-independent search regions targeting squarks and gluinos in models with one-step decay. Expected yields are normalized to a luminosity of 139 $fb^{-1}$.

Cut-flow for BDT search regions targeting gluinos in models with one-step decays. Expected yields are normalized to a luminosity of 139 $fb^{-1}$.

Cut-flow for BDT search regions targeting gluinos in models with one-step decays. Expected yields are normalized to a luminosity of 139 $fb^{-1}$.

Cut-flow for BDT search regions targeting gluinos in models with direct decays. Expected yields are normalized to a luminosity of 139 $fb^{-1}$.

Cut-flow for BDT search regions targeting gluinos in models with direct decays. Expected yields are normalized to a luminosity of 139 $fb^{-1}$.

The \(m_{bb}\) distributions in the MediumBDT categories, overlaid with contributions from the \(H\gamma jj\) signal as well as the resonant \(Z\gamma jj\) and non-resonant \(b\bar{b} \gamma jj\) background fits. The combined \(\chi^2\) per degree of freedom is \(45.2/45\). The bottom panel in each plot presents the significance of the Higgs boson signal relative to the non-resonant \(b\bar{b} \gamma jj\) background in each bin.

The \(m_{bb}\) distributions in the MediumBDT categories, overlaid with contributions from the \(H\gamma jj\) signal as well as the resonant \(Z\gamma jj\) and non-resonant \(b\bar{b} \gamma jj\) background fits. The combined \(\chi^2\) per degree of freedom is \(45.2/45\). The bottom panel in each plot presents the significance of the Higgs boson signal relative to the non-resonant \(b\bar{b} \gamma jj\) background in each bin.

Comparisons of data and simulated event distributions of the BDT input variables \(\Delta R (b1, \gamma)\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(m_{jj}\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(\Delta R (b2, \gamma)\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(\text{centrality}(\gamma, j1, j2)\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(\Delta \phi(b\bar{b}, jj)\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(p_{\mathrm{T}}^{jj}\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(\cos \theta\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Comparisons of data and simulated event distributions of the BDT input variables \(\Delta R(b1,j1)\) in the two \(m_{bb}\) sidebands after kinematic reweighting of the non-resonant \(b\bar{b}\gamma jj\) background. The data are shown as black points, and the background contributions are stacked in coloured histograms. The Higgs boson signal contribution is scaled up and represented by the dashed red line. The bottom panel in each plot shows the ratio between the data and the SM prediction, where the uncertainty band corresponds to the statistical uncertainty only.

Numerical results of $b \bar{b}$- and $c \bar{c}$-dijet cross-sections, $c \bar{c}$/$b \bar{b}$ dijet cross-section ratios and their total uncertainties as a function of the leading jet $p_T$.

Numerical results of $b \bar{b}$- and $c \bar{c}$-dijet cross-sections, $c\bar{c}$/$b \bar{b}$ dijet cross-section ratios and their total uncertainties as a function of $m_{jj}$ (dijet invariant mass).

Covariance matrix, corresponding to the total uncertainties, obtained between the leading jet eta intervals of the $b \bar{b}$-dijet differential cross sections. The unit of all the elements of the matrix is nb$^2$.

Covariance matrix, corresponding to the total uncertainties, obtained between the leading jet eta intervals of the $c \bar{c}$-dijet differential cross sections. The unit of all the elements of the matrix is nb$^2$.

Covariance matrix, corresponding to the total uncertainties, obtained between the leading jet $\eta$ intervals of the $b \bar{b}$ (horizontal) and $c \bar{c}$ (vertical) differential cross sections. The unit of all the elements of the matrix is nb$^2$.

Covariance matrix, corresponding to the total uncertainties, obtained between the $\Delta y^*$ intervals of the $b \bar{b}$-dijet differential cross sections. The unit of all the elements of the matrix is nb$^2$.

Covariance matrix, corresponding to the total uncertainties, obtained between the $\Delta y^*$ intervals of the $c \bar{c}$-dijet differential cross sections. The unit of all the elements of the matrix is nb$^2$.

Covariance matrix, corresponding to the total uncertainties, obtained between the $\Delta y^*$ intervals of the $b \bar{b}$ (horizontal) and $c \bar{c}$ (vertical) differential cross sections. The unit of all the elements of the matrix is nb$^2$.

Covariance matrix, corresponding to the total uncertainties, obtained between the leading jet $p_T$ intervals of the $b \bar{b}$-dijet differential cross sections. The unit of all the elements of the matrix is (nb GeV$/c)^2$ and the $p_T$ intervals are given in GeV$/c$.

Covariance matrix, corresponding to the total uncertainties, obtained between the leading jet $p_T$ intervals of the $c \bar{c}$-dijet differential cross sections. The unit of all the elements of the matrix is (nb GeV$/c)^2$ and the $p_T$ intervals are given in GeV$/$c .

Covariance matrix, corresponding to the total uncertainties, obtained between the leading jet $p_T$ intervals of the $b \bar{b}$ (horizontal) and $c \bar{c}$ (vertical) differential cross sections. The unit of all the elements of the matrix is (nb GeV$/c)^2$ and the $p_T$ intervals are given in GeV$/c$.

Covariance matrix, corresponding to the total uncertainties, obtained between the $m_{jj}$ intervals of the $b \bar{b}$-dijet differential cross sections. The unit of all the elements of the matrix is (nb GeV$/c^2)^2$ and the mass intervals are given in GeV$/c^2$.

Covariance matrix, corresponding to the total uncertainties, obtained between the $m_{jj}$ intervals of the $c \bar{c}$-dijet differential cross sections. The unit of all the elements of the matrix is (nb GeV$/c^2)^2$ and the mass intervals are given in GeV$/c^2$.

Covariance matrix, corresponding to the total uncertainties, obtained between the $m_{jj}$ intervals of the $b \bar{b}$ (horizontal) and $c \bar{c}$ (vertical) differential cross sections. The unit of all the elements of the matrix is $($nb GeV$/c^2)^2$ and the mass intervals are given in GeV$/c^2$.

ALP width dominantly invisible, see caption of Fig 7.

scalar production search, see caption of Fig 8.