The measured fiducial cross section vs $m(\ell\ell\gamma)$. The central values are provided along with the statistical and systematic uncertainties together with the sign information. The statistical and "Uncor" uncertainty should be treated as uncorrelated bin-to-bin, while the rest are correlated between bins, and they are written as signed NP variations. The parton-to-particle correction factor $C_{theory}$ is the ratio of the cross-section predicted by Sherpa LO samples at particle level within the fiducial phase-space region defined in Table 4 to the predicted cross-section at parton level within the same fiducial region but with the smooth-cone isolation prescription defined above replacing the particle-level photon isolation criterion, and with Born-level leptons in place of dressed leptons. This correction should be applied on fixed order parton-level calculations. The systematic uncertainty is evaluated from a comparison with the correction factor obtained using events generated with Sherpa 2.2.2 at NLO. The uncertainty is defined as Max(stat error, systematic difference between Sherpa LO and Sherpa 2.2.2 NLO), and cannot be considered correlated bin-to-bin. In the case that the calculations are valid for dressed leptons, a modified correction factor excluding the Born-to-dressed lepton correction should be applied instead. This correction only takes into account the particle-level isolation criteria, and is provided separately here. The Sherpa 2.2.8 NLO cross-sections given below include a small contribution from EW $Z\gamma jj$ production.

The measured fiducial cross section vs $p_{\mathrm{T}}^{\ell\ell\gamma}$. The central values are provided along with the statistical and systematic uncertainties together with the sign information. The statistical and "Uncor" uncertainty should be treated as uncorrelated bin-to-bin, while the rest are correlated between bins, and they are written as signed NP variations. The parton-to-particle correction factor $C_{theory}$ is the ratio of the cross-section predicted by Sherpa LO samples at particle level within the fiducial phase-space region defined in Table 4 to the predicted cross-section at parton level within the same fiducial region but with the smooth-cone isolation prescription defined above replacing the particle-level photon isolation criterion, and with Born-level leptons in place of dressed leptons. This correction should be applied on fixed order parton-level calculations. The systematic uncertainty is evaluated from a comparison with the correction factor obtained using events generated with Sherpa 2.2.2 at NLO. The uncertainty is defined as Max(stat error, systematic difference between Sherpa LO and Sherpa 2.2.2 NLO), and cannot be considered correlated bin-to-bin. In the case that the calculations are valid for dressed leptons, a modified correction factor excluding the Born-to-dressed lepton correction should be applied instead. This correction only takes into account the particle-level isolation criteria, and is provided separately here. The Sherpa 2.2.8 NLO cross-sections given below include a small contribution from EW $Z\gamma jj$ production.

The measured fiducial cross section vs $\Delta\phi(ll,\gamma)$. The central values are provided along with the statistical and systematic uncertainties together with the sign information. The statistical and "Uncor" uncertainty should be treated as uncorrelated bin-to-bin, while the rest are correlated between bins, and they are written as signed NP variations. The parton-to-particle correction factor $C_{theory}$ is the ratio of the cross-section predicted by Sherpa LO samples at particle level within the fiducial phase-space region defined in Table 4 to the predicted cross-section at parton level within the same fiducial region but with the smooth-cone isolation prescription defined above replacing the particle-level photon isolation criterion, and with Born-level leptons in place of dressed leptons. This correction should be applied on fixed order parton-level calculations. The systematic uncertainty is evaluated from a comparison with the correction factor obtained using events generated with Sherpa 2.2.2 at NLO. The uncertainty is defined as Max(stat error, systematic difference between Sherpa LO and Sherpa 2.2.2 NLO), and cannot be considered correlated bin-to-bin. In the case that the calculations are valid for dressed leptons, a modified correction factor excluding the Born-to-dressed lepton correction should be applied instead. This correction only takes into account the particle-level isolation criteria, and is provided separately here. The Sherpa 2.2.8 NLO cross-sections given below include a small contribution from EW $Z\gamma jj$ production.

The measured fiducial cross section vs $p_{\mathrm{T}}^{\ell\ell\gamma}/m(\ell\ell\gamma)$. The central values are provided along with the statistical and systematic uncertainties together with the sign information. The statistical and "Uncor" uncertainty should be treated as uncorrelated bin-to-bin, while the rest are correlated between bins, and they are written as signed NP variations. The parton-to-particle correction factor $C_{theory}$ is the ratio of the cross-section predicted by Sherpa LO samples at particle level within the fiducial phase-space region defined in Table 4 to the predicted cross-section at parton level within the same fiducial region but with the smooth-cone isolation prescription defined above replacing the particle-level photon isolation criterion, and with Born-level leptons in place of dressed leptons. This correction should be applied on fixed order parton-level calculations. The systematic uncertainty is evaluated from a comparison with the correction factor obtained using events generated with Sherpa 2.2.2 at NLO. The uncertainty is defined as Max(stat error, systematic difference between Sherpa LO and Sherpa 2.2.2 NLO), and cannot be considered correlated bin-to-bin. In the case that the calculations are valid for dressed leptons, a modified correction factor excluding the Born-to-dressed lepton correction should be applied instead. This correction only takes into account the particle-level isolation criteria, and is provided separately here. The Sherpa 2.2.8 NLO cross-sections given below include a small contribution from EW $Z\gamma jj$ production.

Expected and observed 95% CL upper limits on $\sigma(pp\to H)B(H\to Z(\eta_c~\text{or}~J/\psi))/$pb. The smaller (larger) quoted ranges around the expected limits represent $\pm 1\sigma$ ($\pm 2\sigma$) fluctuations.

Lower limits at 95$\%$ CL on $\Lambda$ for the dimuon channel for different signal chiralities in the (constructive/destructive interference) SRs of the analysis.

Lower limits at 95$\%$ CL on $\Lambda$ for the dilepton combination for different signal chiralities in the (constructive/destructive interference) SRs of the analysis.

The dielectron data events passing the full selection.

The dimuon data events passing the full selection.

The absolute differential cross-section measured in the fiducial phase-space as a function of the minimum $\Delta R$ between the photon and the leptons in the electron-muon channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute differential cross-section measured in the fiducial phase-space as a function of the $\Delta\phi$ between the two leptons in the electron-muon channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The absolute differential cross-section measured in the fiducial phase-space as a function of the $|\Delta\eta|$ between the two leptons in the electron-muon channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised differential cross-section measured in the fiducial phase-space as a function of the photon pT in the electron-muon channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised differential cross-section measured in the fiducial phase-space as a function of the photon $|\eta|$ in the electron-muon channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised differential cross-section measured in the fiducial phase-space as a function of the minimum $\Delta R$ between the photon and the leptons in the electron-muon channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised differential cross-section measured in the fiducial phase-space as a function of the $\Delta\phi$ between the two leptons in the electron-muon channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The normalised differential cross-section measured in the fiducial phase-space as a function of the $|\Delta\eta|$ between the two leptons in the electron-muon channel. The uncertainty is decomposed into four components which are the signal modelling uncertainty, the background modelling uncertainty, the experimental uncertainty, and the data statistical uncertainty.

The total correlation matrix of the absolute differential cross-section measured in the fiducial phase-space as a function of the photon pT in the electron-muon channel. The individual systematic uncertainties are symmetrized before deriving the correlation matrix.

The total correlation matrix of the absolute differential cross-section measured in the fiducial phase-space as a function of the photon $|\eta|$ in the electron-muon channel. The individual systematic uncertainties are symmetrized before deriving the correlation matrix.

The total correlation matrix of the absolute differential cross-section measured in the fiducial phase-space as a function of the minimum $\Delta R$ between the photon and the leptons in the electron-muon channel. The individual systematic uncertainties are symmetrized before deriving the correlation matrix.

The total correlation matrix of the absolute differential cross-section measured in the fiducial phase-space as a function of the $\Delta\phi$ between the two leptons in the electron-muon channel. The individual systematic uncertainties are symmetrized before deriving the correlation matrix.

The total correlation matrix of the absolute differential cross-section measured in the fiducial phase-space as a function of the $|\Delta\eta|$ between the two leptons in the electron-muon channel. The individual systematic uncertainties are symmetrized before deriving the correlation matrix.

The total correlation matrix of the normalised differential cross-section measured in the fiducial phase-space as a function of the photon pT in the electron-muon channel. The individual systematic uncertainties are symmetrized before deriving the correlation matrix.

The total correlation matrix of the normalised differential cross-section measured in the fiducial phase-space as a function of the photon $|\eta|$ in the electron-muon channel. The individual systematic uncertainties are symmetrized before deriving the correlation matrix.

The total correlation matrix of the normalised differential cross-section measured in the fiducial phase-space as a function of the minimum $\Delta R$ between the photon and the leptons in the electron-muon channel. The individual systematic uncertainties are symmetrized before deriving the correlation matrix.

The total correlation matrix of the normalised differential cross-section measured in the fiducial phase-space as a function of the $\Delta\phi$ between the two leptons in the electron-muon channel. The individual systematic uncertainties are symmetrized before deriving the correlation matrix.

The total correlation matrix of the normalised differential cross-section measured in the fiducial phase-space as a function of the $|\Delta\eta|$ between the two leptons in the electron-muon channel. The individual systematic uncertainties are symmetrized before deriving the correlation matrix.

The statistical correlation matrix of all the absolute differential cross-sections measured in the fiducial phase-space in the electron-muon channel.

The statistical correlation matrix of all the normalised differential cross-sections measured in the fiducial phase-space in the electron-muon channel.

Fiducial region definition.

Impact of different categories of systematic uncertainties on the fiducial and inclusive measurements. The quoted values are obtained by repeating the fit, fixing a set of nuisance parameters of the sources corresponding to the considered category, and subtracting in quadrature the resulting uncertainty from the total uncertainty of the nominal fit. The total uncertainty is different from the sum in quadrature of the different components due to correlations between nuisance parameters built by the fit.

Fiducial region definition

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

Expected and observed events in the 16 discovery regions along with the according control regions.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the MS-agnostic R-hadron search.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector R-hadron search.

p0-values and model-independent upper limits on cross-section x acceptance x efficiency for the 16 discovery regions.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the MS-agnostic search for metastable gluino R-hadrons.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector direct-stau search.

Expected signal yield and acceptance x efficiency, estimated background and observed number of events in data for the full range of simulated masses in the full-detector chargino search.

Upper cross-section limit in gluino R-hadron search.

Upper cross-section limit in sbottom R-hadron search.

Upper cross-section limit in stop R-hadron search.

Upper cross-section limit in stau search.

Upper cross-section limit in chargino search.

Lower mass limit as function of gluino lifetime.

Acceptance x efficiency, acceptance and efficiency for the full range of simulated masses in the MS-agnostic R-hadron search.

Upper cross-section limit in meta-stable gluino R-hadron search.

Flavor composition of 800 GeV stop R-hadrons simulated using the generic model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV anti-stop R-hadrons simulated using the generic model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV stop R-hadrons simulated using the Regge model as a function of radial distance from the interaction point.

Flavor composition of 800 GeV anti-stop R-hadrons simulated using the Regge model as a function of radial distance from the interaction point.

ETmiss trigger efficiency as function of true ETmiss (EtmissTurnOn).

Single-muon trigger efficiency as function of $|\eta|$ and $\beta$ (SingleMuTurnOn).

Candidate reconstruction efficiency for ID+Calo selection (IDCaloEff).

Candidate reconstruction efficiency for loose selection (LooseEff).

Efficiency for a loose candidate to be promoted to a tight candidate (TightPromotionEff).

Resolution and average of reconstructed dE/dx mass for a given simulated mass for ID+calo candidates.

Resolution and average of reconstructed ToF mass for a given simulated mass for ID+calo candidates.

Resolution and average of reconstructed ToF mass for a given simulated mass for FullDet candidates.

Track multiplicity $n_{Tracks}$ for preselected DVs in muon-triggered events with at least one muon passing the full selection. Along with the data shown with black markers, the stacked filled histograms represent the background estimates, and predictions for signal scenarios are overlaid with dashed lines. The errors include statistical and systematic uncertainties and are indicated by hatched bands. The DV full selection requirements, $n_{Tracks} \geq 3$ and $m_{DV} > 20$ GeV are visualized with a black arrow.

Invariant mass $m_{DV}$ for the highest-mass preselected DV with at least three associated tracks in MET-triggered events with at least one muon passing the full selection. Along with the data shown with black markers, the stacked filled histograms represent the background estimates, and predictions for signal scenarios are overlaid with dashed lines. The errors include statistical and systematic uncertainties and are indicated by hatched bands. The DV full selection requirements, $n_{Tracks} \geq 3$ and $m_{DV} > 20$ GeV are visualized with a black arrow.

Invariant mass $m_{DV}$ for the highest-mass preselected DV with at least three associated tracks in muon-triggered events with at least one muon passing the full selection. Along with the data shown with black markers, the stacked filled histograms represent the background estimates, and predictions for signal scenarios are overlaid with dashed lines. The errors include statistical and systematic uncertainties and are indicated by hatched bands. The DV full selection requirements, $n_{Tracks} \geq 3$ and $m_{DV} > 20$ GeV are visualized with a black arrow.

The observed event yields in the control, validation and signal regions are shown for the MET Trigger selections, along with the predicted background yields. The bottom panel shows the ratio of observed events to the total background yields. The errors represent the total uncertainty of the backgrounds prediction, including the statistical and systematic uncertainties added in quadrature.

The observed event yields in the control, validation and signal regions are shown for the Muon Trigger selections, along with the predicted background yields. The bottom panel shows the ratio of observed events to the total background yields. The errors represent the total uncertainty of the backgrounds prediction, including the statistical and systematic uncertainties added in quadrature.

Expected exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Expected (1 sigma band) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Expected (2 sigma band) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Observed exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Observed (+1 sigma) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Observed (-1 sigma) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Exclusion limits on the production cross section as a function of m($\tilde{t}$) are shown for several values of $\tau(\tilde{t})$ along with the nominal signal production cross section and its theoretical uncertainty.

Parameterized event selection efficiencies for the $E_{T}^{miss}$ Trigger SR. The event-level efficiencies for each SR are extracted for all events passing the acceptance of the corresponding SR.

Parameterized event selection efficiencies for the Muon Trigger SR. The event-level efficiencies for each SR are extracted for all events passing the acceptance of the corresponding SR.

Parameterized muon-level reconstruction efficiencies as a function of the muon $p_{T}$ and $d_{0}$. The muon-level efficiencies are extracted using muons passing the muon acceptance criteria.

Parameterized vertex-level reconstruction efficiencies as a function of the radial position of the truth vertex. The efficiency is calculated independent of the muons originating from this truth vertex.

Parameterized vertex-level reconstruction efficiencies as a function of the radial position of the truth vertex. The efficiency is calculated only for truth vertices which have a muon originating from them which is matched to a reconstructed muon.

The $p_{T}$ distribution of all muons originating from LLP decays in the samples used to calculate and validate the efficiencies.

The invariant mass and multiplicity of selected decay products of all truth vertices used in the calculation and validation of the reconstructed efficiencies.

Observed $m_{t\bar{t}}^{reco}$ distributions in data for SR$2b$, shown together with the result of the fit with the three-shape-parameter function. The error bars indicate the effect of the fit parameter uncertainty on the background prediction. The bin width of the distributions is chosen to be the same as that used in the background parameterization.

Observed and expected upper limits on the cross-section times branching fraction of the $Z'_{TC2} \rightarrow t\bar{t}$ as a function of the $Z'_{TC2}$ mass ($\Gamma/m = 1.2$%).

Fit parameter values obtained from fits to data in the two signal regions. The fit function is $f(x) = p_{0} (1-x)^{p_{1}} x^{p_{2} + p_{3} \log(x)}$, where $x = m^{reco}_{t\bar{t}}/\sqrt{s}$.

Distribution of $v_{3}$ from MBT events plotted as a function of the A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of $v_{3}$ from $p_{T}^{jet}>75$ GeV events plotted as a function of the A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of $v_{3}$ from $p_{T}^{jet}>100$ GeV events plotted as a function of the A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of $v_{2}$ from MBT events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>75$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>100$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of $v_{2}$ from MBT events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>75$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>100$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of $v_{2}$ from MBT events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>75$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>100$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative UE-UE pair fractions from MBT events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-HS pair fractions from MBT events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-UE pair fractions from MBT events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-HS pair fractions from MBT events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-HS pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-HS pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-HS pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-HS pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative HS-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative UE-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative HS-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative UE-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative HS-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).