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.

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

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 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 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 three exclusive signal regions.

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.

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 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 ttW and ttZ.

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

Charged particle mean multiplicities measured in 4x4 kinematic bins of $Q^2$ and $y$.

Charged particle mean multiplicities measured with the additional restriction to select only particles from the current region of the Breit frame $0<\eta^{*}<4$, in 4x4 kinematic bins of $Q^2$ and $y$.

Variance of the charged particle multiplicity measured in 4x4 kinematic bins of $Q^2$ and $y$.

Variance of the charged particle multiplicity measured with the additional restriction to select only particles from the current region of the Breit frame $0<\eta^{*}<4$, in 4x4 kinematic bins of $Q^2$ and $y$.

KNO function for charged particles neasured with the additional restriction to select only particles from the current region of the Breit frame $0<\eta^{*}<4$, in 4x4 kinematic bins of $Q^2$ and $y$.

Hadron entropy for charged particles selected in a sliding pseudorapidity window of 1.4 units size with transverse momenta $P_{T lab}>150$ MeV, in the accessible kinematic bins of momentum transfer $Q^2$ and $y$. The data are presented as a function of $Q^2$ and the corresponding average $x$-Bjorken.

Hadron entropy for charged particles neasured with the additional restriction to select only particles from the current region of the Breit frame $0<\eta^{*}<4$, in 4x4 kinematic bins of $Q^2$ and $y$. The data are presented as a function of $Q^2$ and the corresponding average $x$-Bjorken.

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

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

The ratios (&mu;) of the 95&#37; C.L. upper limits on the combined s&rarr; W^&plusmn;W^&#8723; and s&rarr; ZZ cross section to simplified model expectations for the m_Z'=1 TeV scenario, for various m_s hypotheses. The observed limits (solid line) are consistent with the expectation under the SM-only hypothesis (dashed line) within uncertainties (filled band), except for a small excess for m_s=160 GeV, discussed in the text.

The ratios (&mu;) of the 95&#37; C.L. upper limits on the combined s&rarr; W^&plusmn;W^&#8723; and s&rarr; ZZ cross section to simplified model expectations for the m_Z'=1.7 TeV scenario, for various m_s hypotheses. The observed limits (solid line) are consistent with the expectation under the SM-only hypothesis (dashed line) within uncertainties (filled band), except for a small excess for m_s=160 GeV, discussed in the text.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) × B(s&rarr; VV) for m_Z'=0.5 TeV signal points. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) × B(s&rarr; VV) process for m_Z'=0.5 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) × B(s&rarr; VV) for m_Z'=1 TeV signal points. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) × B(s&rarr; VV) process for m_Z'=1 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) × B(s&rarr; VV) for m_Z'=1.7 TeV signal points. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) × B(s&rarr; VV) process for m_Z'=1.7 TeV, shown in dashed blue.

SM background post-fit yields stacked in each SR and CR category and E_T^miss bin and data overlaid with the maximum likelihood estimators set to the conditional values of the combined signal and control region fit. The hatched uncertainty band shown includes simulation statistics uncertainties, experimental systematic uncertainties, and V+jets theory modelling systematic uncertainties. Pre-fit uncertainties cover differences between the data and pre-fit background prediction.

Cumulative efficiencies for the merged category for signal samples with m_s=160 GeV (a), m_s=235 GeV (b) and m_s=310 GeV (c), each with m_Z'=1 TeV. The dark Higgs candidate selection includes stringent jet substructure requirements and typically at most one candidate is present in signal events. Here, &Delta; &phi;_{jets_1,2,3 E_T^miss} is the smallest azimuthal angle between the E_T^miss and any of the three highest-p_T (leading) small-R jets.

Cumulative efficiencies for the intermediate category for signal samples with m_s=160 GeV (a), m_s=235 GeV (b) and m_s=310 GeV (c), each with m_Z'=1 TeV. The TAR+Comb algorithm reconstructs the dark Higgs candidate from a TAR jet with m^TAR&gt;60 GeV that is supplemented by up to two additional small-R jets within &Delta;R_cone=2.5 of the TAR jet. Here, &Delta; &phi;_{jets_1,2,3 E_T^miss} is the smallest azimuthal angle between the E_T^miss and any of the three highest-p_T (leading) small-R jets. For details see text.

The product of acceptance and efficiency (A &times; &varepsilon;), defined as the number of signal events satisfying the full set of selection criteria in the merged or intermediate signal regions, divided by the total number of generated signal events, for the s(W^&plusmn;W^&#8723;) dark Higgs signal points with dark Higgs boson mass m_s and Z' boson mass m_Z'.

The product of acceptance and efficiency (A &times; &varepsilon;), defined as the number of signal events satisfying the full set of selection criteria in the merged or intermediate signal regions, divided by the total number of generated signal events, for the s(ZZ) dark Higgs signal points with dark Higgs boson mass m_s and Z' boson mass m_Z'.