Measured differential yield of charged hadrons as a function of transverse momentum for pseudorapidities 0.1, 0.3, 0.5 and 0.7 for centre-of-mass energy 2360 GeV.

Measured differential yield of charged hadrons as a function of transverse momentum for pseudorapidities 0.9, 1.1, 1.3 and 1.5 for centre-of-mass energy 2360 GeV.

Measured differential yield of charged hadrons as a function of transverse momentum for pseudorapidities 1.7, 1.9, 2.1 and 2.3 for centre-of-mass energy 2360 GeV.

Measured differential yield of charged hadrons as a function of transverse momentum for the pseudorapidity range -2.4 to +2.4 for centre-of-mass energies 900 and 2360 GeV.

Reconstructed differential charged hadron yields as a function of pseudorapidity averaged over the cluster counting, tracklet and tracking methods at centre-of-mass energies of 900 and 2360 GeV.

Differential cross section of J/psi mesons from b hadrons as a function of dimuon pT in pp and PbPb collisions. The PbPb cross sections are normalised by TAA for direct comparison. Global uncertainties arise from the integrated luminosity uncertainty in pp collisions, and the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross section of prompt J/psi mesons as a function of dimuon rapidty in pp and PbPb collisions. The PbPb cross sections are normalised by TAA for direct comparison. Global uncertainties arise from the integrated luminosity uncertainty in pp collisions, and the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross section of J/psi mesons from b hadrons as a function of dimuon rapidty in pp and PbPb collisions. The PbPb cross sections are normalised by TAA for direct comparison. Global uncertainties arise from the integrated luminosity uncertainty in pp collisions, and the number of minimum bias events and TAA uncertainties for PbPb collisions.

Nuclear modification factor of prompt J/psi mesons as a function of dimuon rapidity. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of collision centrality. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of dimuon pT. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of dimuon pT in the mid-rapidity interval. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of dimuon pT in the most forward rapidity interval. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of collision centrality in the mid- and most forward rapidity intervals. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of pT in three centrality bins. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of collision centrality in the most forward rapidity range at moderate and high pT. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of collision centrality in the mid-rapidity range. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt psi(2S) mesons as a function of collision centrality in the mid-rapidity range. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

95% CL intervals on the nuclear modification factor of prompt psi(2S) mesons as a function of collision centrality in the mid-rapidity range. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt J/psi and psi(2S) mesons as a function of dimuon pT in the mid-rapidity range. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of collision centrality in the forward rapidity range. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt psi(2S) mesons as a function of collision centrality in the forward rapidity range. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

95% CL intervals on the nuclear modification factor of prompt psi(2S) mesons as a function of collision centrality in the forward rapidity range. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of dimuon pT in the forward rapidity range. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt psi(2S) mesons as a function of dimuon pT in the forward rapidity range. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

95% CL intervals on the nuclear modification factor of prompt psi(2S) mesons as a function of dimuon pT in the forward rapidity range. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of dimuon rapidity. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of collision centrality. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of dimuon pT. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of dimuon pT in the mid-rapidity interval. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of dimuon pT in the most forward rapidity interval. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of collision centrality in the mid- and most forward rapidity intervals. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of pT in three centrality bins. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of collision centrality in the most forward rapidity range at moderate and high pT. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Covariance matrix of the fitted parameters of interest.

Covariance matrix of the fitted parameters of interest.

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

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

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

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

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

The product of signal acceptance and efficiency in the 2l categories for the signal produced via vector boson fusion.

$m_{X}^{T}$ distribution in data in the 0l 2b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}^{T}$ distribution in data in the 0l $\leq$1b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}$ distribution in data in the 2e 2b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2e $\leq$1b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ 2b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ $\leq$1b non-VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}^{T}$ distribution in data in the 0l 2b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}^{T}$ distribution in data in the 0l $\leq$1b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}^{T}$ observed in the SR.

$m_{X}$ distribution in data in the 2e 2b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2e $\leq$1b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ 2b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

$m_{X}$ distribution in data in the 2$\mu$ $\leq$1b VBF category. The data is shown as Events/10 GeV. The distribution is shown up 4000 GeV, which corresponds to the event with the highest $m_{X}$ observed in the SR.

Observed and expected 95% CL upper limit on $\sigma \mathcal{B}$(Z'-> ZH) with all categories combined for the non-VBF signal, including all statistical and systematic uncertainties. The inner green band and the outer yellow band indicate the regions containing 68 and 95%, respectively, of the distribution of expected limits under the background-only hypothesis. The CMS search for a heavy resonance using 2016 data and the same final state [JHEP 11 (2018) 172] is shown as a comparison.

Observed and expected 95% CL upper limit on $\sigma \mathcal{B}$(Z'-> ZH) with all categories combined for the VBF signal, including all statistical and systematic uncertainties. The inner green band and the outer yellow band indicate the regions containing 68 and 95%, respectively, of the distribution of expected limits under the background-only hypothesis.

Observed exclusion limit in the space of the HVT model parameters [$g_{V}c_{H}$, $g^{2}c_{F}/g_{V}$] for mass hypotheses of 2 TeV for the non-VBF signal.

Observed exclusion limit in the space of the HVT model parameters [$g_{V}c_{H}$, $g^{2}c_{F}/g_{V}$] for mass hypotheses of 3 TeV for the non-VBF signal.

Observed exclusion limit in the space of the HVT model parameters [$g_{V}c_{H}$, $g^{2}c_{F}/g_{V}$] for mass hypotheses of 4 TeV for the non-VBF signal.

Three-jet events $p_{\mathrm{T}3}/p_{\mathrm{T}2}$ for large-angle radiation ($\Delta R_{23}$ > 1.0)

Three-jet events $\Delta R_{23}$ for soft radiation ($p_{\mathrm{T}3}/p_{\mathrm{T}2}$ < 0.3)

Three-jet events $\Delta R_{23}$ for soft radiation ($p_{\mathrm{T}3}/p_{\mathrm{T}2}$ < 0.3)

Three-jet events $\Delta R_{23}$ for hard radiation ($p_{\mathrm{T}3}/p_{\mathrm{T}2}$ > 0.6)

Three-jet events $\Delta R_{23}$ for hard radiation ($p_{\mathrm{T}3}/p_{\mathrm{T}2}$ > 0.6)

Three-jet events $p_{\mathrm{T}3}/p_{\mathrm{T}2}$ for small-angle radiation ($\Delta R_{23}$ < 1.0)

Three-jet events $p_{\mathrm{T}3}/p_{\mathrm{T}2}$ for small-angle radiation ($\Delta R_{23}$ < 1.0)

Three-jet events $p_{\mathrm{T}3}/p_{\mathrm{T}2}$ for large-angle radiation ($\Delta R_{23}$ > 1.0)

Three-jet events $p_{\mathrm{T}3}/p_{\mathrm{T}2}$ for large-angle radiation ($\Delta R_{23}$ > 1.0)

Three-jet events $\Delta R_{23}$ for soft radiation ($p_{\mathrm{T}3}/p_{\mathrm{T}2}$ < 0.3)

Three-jet events $\Delta R_{23}$ for soft radiation ($p_{\mathrm{T}3}/p_{\mathrm{T}2}$ < 0.3)

Three-jet events $\Delta R_{23}$ for hard radiation ($p_{\mathrm{T}3}/p_{\mathrm{T}2}$ > 0.6)

Three-jet events $\Delta R_{23}$ for hard radiation ($p_{\mathrm{T}3}/p_{\mathrm{T}2}$ > 0.6)

Z + two-jet events $p_{\mathrm{T}3}/p_{\mathrm{T}2}$ for small-angle radiation ($\Delta R_{23}$ < 1.0)

Z + two-jet events $p_{\mathrm{T}3}/p_{\mathrm{T}2}$ for small-angle radiation ($\Delta R_{23}$ < 1.0)

Z + two-jet events $p_{\mathrm{T}3}/p_{\mathrm{T}2}$ for large-angle radiation ($\Delta R_{23}$ > 1.0)

Z + two-jet events $p_{\mathrm{T}3}/p_{\mathrm{T}2}$ for large-angle radiation ($\Delta R_{23}$ > 1.0)

Z + two-jet events $\Delta R_{23}$ for soft radiation ($p_{\mathrm{T}3}/p_{\mathrm{T}2}$ < 0.3)

Z + two-jet events $\Delta R_{23}$ for soft radiation ($p_{\mathrm{T}3}/p_{\mathrm{T}2}$ < 0.3)

Z + two-jet events $\Delta R_{23}$ for hard radiation ($p_{\mathrm{T}3}/p_{\mathrm{T}2}$ > 0.6)

Z + two-jet events $\Delta R_{23}$ for hard radiation ($p_{\mathrm{T}3}/p_{\mathrm{T}2}$ > 0.6)

Product of acceptance and efficiency for dimuon pairs as a function of generated mass in simulated events. The DY samples are used to represent spin-1 particles, and RS graviton samples are used for spin-2 particles.

The invariant mass spectra of dielectron events. The points with error bars represent the observed yield. The histogram represents the expectations from the SM processes. The bins have equal width in logarithmic scale so that the width in GeV becomes larger with increasing mass.

The invariant mass spectra of dielectron events. The points with error bars represent the observed yield. The histogram represents the expectations from the SM processes. The bins have equal width in logarithmic scale so that the width in GeV becomes larger with increasing mass.

The invariant mass spectra of dimuon events. The points with error bars represent the observed yield. The histogram represents the expectations from the SM processes. The bins have equal width in logarithmic scale so that the width in GeV becomes larger with increasing mass.

The invariant mass spectra of dimuon events. The points with error bars represent the observed yield. The histogram represents the expectations from the SM processes. The bins have equal width in logarithmic scale so that the width in GeV becomes larger with increasing mass.

The invariant mass distribution of pairs of electrons observed in data and expected from simulated SM processes for $\cos\theta^* < 0$. The bin width gradually increases with mass. The quoted uncertainty represents the combined statistical and systematic uncertainties on the background. The signal contribution expected from the constructive LL CI model with $\Lambda = 16$ TeV is shown as a dashed green line.

The invariant mass distribution of pairs of electrons observed in data and expected from simulated SM processes for $\cos\theta^* < 0$. The bin width gradually increases with mass. The quoted uncertainty represents the combined statistical and systematic uncertainties on the background. The signal contribution expected from the constructive LL CI model with $\Lambda = 16$ TeV is shown as a dashed green line.

The invariant mass distribution of pairs of electrons observed in data and expected from simulated SM processes for $\cos\theta^* /geq 0$. The bin width gradually increases with mass. The quoted uncertainty represents the combined statistical and systematic uncertainties on the background. The signal contribution expected from the constructive LL CI model with $\Lambda = 16$ TeV is shown as a dashed green line.

The invariant mass distribution of pairs of electrons observed in data and expected from simulated SM processes for $\cos\theta^* /geq 0$. The bin width gradually increases with mass. The quoted uncertainty represents the combined statistical and systematic uncertainties on the background. The signal contribution expected from the constructive LL CI model with $\Lambda = 16$ TeV is shown as a dashed green line.

The invariant mass distribution of pairs of muons observed in data and expected from simulated SM processes for $\cos\theta^* < 0$. The bin width gradually increases with mass. The quoted uncertainty represents the combined statistical and systematic uncertainties on the background. The signal contribution expected from the constructive LL CI model with $\Lambda = 16$ TeV is shown as a dashed green line.

The invariant mass distribution of pairs of muons observed in data and expected from simulated SM processes for $\cos\theta^* < 0$. The bin width gradually increases with mass. The quoted uncertainty represents the combined statistical and systematic uncertainties on the background. The signal contribution expected from the constructive LL CI model with $\Lambda = 16$ TeV is shown as a dashed green line.

The invariant mass distribution of pairs of muons observed in data and expected from simulated SM processes for $\cos\theta^* /geq 0$. The bin width gradually increases with mass. The quoted uncertainty represents the combined statistical and systematic uncertainties on the background. The signal contribution expected from the constructive LL CI model with $\Lambda = 16$ TeV is shown as a dashed green line.

The invariant mass distribution of pairs of muons observed in data and expected from simulated SM processes for $\cos\theta^* /geq 0$. The bin width gradually increases with mass. The quoted uncertainty represents the combined statistical and systematic uncertainties on the background. The signal contribution expected from the constructive LL CI model with $\Lambda = 16$ TeV is shown as a dashed green line.

Observed and expected background yields for different mass ranges in the dielectron channel. The sum of all background contributions is shown as well as a breakdown into the three main categories. The quoted uncertainty includes both the statistical and the systematic components.

Observed and expected background yields for different mass ranges in the dielectron channel. The sum of all background contributions is shown as well as a breakdown into the three main categories. The quoted uncertainty includes both the statistical and the systematic components.

Observed and expected background yields for different mass ranges in the dimuon channel. The sum of all background contributions is shown as well as a breakdown into the three main categories. The quoted uncertainty includes both the statistical and the systematic components.

Observed and expected background yields for different mass ranges in the dimuon channel. The sum of all background contributions is shown as well as a breakdown into the three main categories. The quoted uncertainty includes both the statistical and the systematic components.

The observed upper limits at 95% CL on the product of production cross section and branching fraction for a spin-1 resonance with a width equal to 0.6% of the resonance mass, relative to the product of production cross section and branching fraction of a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The observed upper limits at 95% CL on the product of production cross section and branching fraction for a spin-1 resonance with a width equal to 0.6% of the resonance mass, relative to the product of production cross section and branching fraction of a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The expected upper limits together with the 68 and 95% quantiles at 95% CL on the product of production cross section and branching fraction for a spin-1 resonance with a width equal to 0.6% of the resonance mass, relative to the product of production cross section and branching fraction of a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The expected upper limits together with the 68 and 95% quantiles at 95% CL on the product of production cross section and branching fraction for a spin-1 resonance with a width equal to 0.6% of the resonance mass, relative to the product of production cross section and branching fraction of a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The observed upper limits at 95% CL on the product of production cross section and branching fraction for a spin-1 resonance, for widths equal to 0.6, 3, 5, and 10% of the resonance mass, relative to the product of production cross section and branching fraction for a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The observed upper limits at 95% CL on the product of production cross section and branching fraction for a spin-1 resonance, for widths equal to 0.6, 3, 5, and 10% of the resonance mass, relative to the product of production cross section and branching fraction for a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The expected upper limits at 95% CL on the product of production cross section and branching fraction for a spin-1 resonance, for widths equal to 0.6, 3, 5, and 10% of the resonance mass, relative to the product of production cross section and branching fraction for a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The expected upper limits at 95% CL on the product of production cross section and branching fraction for a spin-1 resonance, for widths equal to 0.6, 3, 5, and 10% of the resonance mass, relative to the product of production cross section and branching fraction for a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

Limits in the ( $c_{d}$ , $c_{u}$ ) plane obtained by recasting the combined limit at 95% CL on the $Z'$ boson cross section from dielectron and dimuon channels. The closed contour representing the GSM model class is composed of thick line segments. Each point on a segment corresponds to a particular model, and the location of the point gives the mass limit on the relevant $Z'$ boson. As indicated in the bottom left legend, the segment line styles correspond to ranges of the particular mixing angle for each considered model. The bottom right legend indicates the constituents of each model class.

Limits in the ( $c_{d}$ , $c_{u}$ ) plane obtained by recasting the combined limit at 95% CL on the $Z'$ boson cross section from dielectron and dimuon channels. The closed contour representing the GSM model class is composed of thick line segments. Each point on a segment corresponds to a particular model, and the location of the point gives the mass limit on the relevant $Z'$ boson. As indicated in the bottom left legend, the segment line styles correspond to ranges of the particular mixing angle for each considered model. The bottom right legend indicates the constituents of each model class.

Limits in the ( $c_{d}$ , $c_{u}$ ) plane obtained by recasting the combined limit at 95% CL on the $Z'$ boson cross section from dielectron and dimuon channels. The closed contour representing the LR model class is composed of thick line segments. Each point on a segment corresponds to a particular model, and the location of the point gives the mass limit on the relevant $Z'$ boson. As indicated in the bottom left legend, the segment line styles correspond to ranges of the particular mixing angle for each considered model. The bottom right legend indicates the constituents of each model class.

Limits in the ( $c_{d}$ , $c_{u}$ ) plane obtained by recasting the combined limit at 95% CL on the $Z'$ boson cross section from dielectron and dimuon channels. The closed contour representing the LR model class is composed of thick line segments. Each point on a segment corresponds to a particular model, and the location of the point gives the mass limit on the relevant $Z'$ boson. As indicated in the bottom left legend, the segment line styles correspond to ranges of the particular mixing angle for each considered model. The bottom right legend indicates the constituents of each model class.

Limits in the ( $c_{d}$ , $c_{u}$ ) plane obtained by recasting the combined limit at 95% CL on the $Z'$ boson cross section from dielectron and dimuon channels. The closed contour representing the E6 model class is composed of thick line segments. Each point on a segment corresponds to a particular model, and the location of the point gives the mass limit on the relevant $Z'$ boson. As indicated in the bottom left legend, the segment line styles correspond to ranges of the particular mixing angle for each considered model. The bottom right legend indicates the constituents of each model class.

Limits in the ( $c_{d}$ , $c_{u}$ ) plane obtained by recasting the combined limit at 95% CL on the $Z'$ boson cross section from dielectron and dimuon channels. The closed contour representing the E6 model class is composed of thick line segments. Each point on a segment corresponds to a particular model, and the location of the point gives the mass limit on the relevant $Z'$ boson. As indicated in the bottom left legend, the segment line styles correspond to ranges of the particular mixing angle for each considered model. The bottom right legend indicates the constituents of each model class.

The observed local p-value for the dielectron channel as a function of the dilepton invariant mass.

The observed local p-value for the dielectron channel as a function of the dilepton invariant mass.

The observed local p-value for the dimuon channel as a function of the dilepton invariant mass.

The observed local p-value for the dimuon channel as a function of the dilepton invariant mass.

The observed local p-value for the combined channel as a function of the dilepton invariant mass.

The observed local p-value for the combined channel as a function of the dilepton invariant mass.

The observed upper limits at 95% CL on the product of production cross section and branching fraction for a spin-2 resonance with a width equal to 0.6% of the resonance mass, relative to the product of production cross section and branching fraction of a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The observed upper limits at 95% CL on the product of production cross section and branching fraction for a spin-2 resonance with a width equal to 0.6% of the resonance mass, relative to the product of production cross section and branching fraction of a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The expected upper limits together with the 68 and 95% quantiles at 95% CL on the product of production cross section and branching fraction for a spin-2 resonance with a width equal to 0.6% of the resonance mass, relative to the product of production cross section and branching fraction of a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

The expected upper limits together with the 68 and 95% quantiles at 95% CL on the product of production cross section and branching fraction for a spin-2 resonance with a width equal to 0.6% of the resonance mass, relative to the product of production cross section and branching fraction of a Z boson, multiplied by the theoretical value of the cross section for $\sigma(pp \to \mathrm{Z}+X \to\ell\ell+X)$ of 1928 pb.

Limits at 95% confidence level for the masses of the DM particle, which is assumed to be Dirac fermion, and its associated mediator, in a simplified model of DM production via a vector mediator. The parameter exclusion is obtained by comparing the limits on product of the production cross section and the branching fraction for decay to a Z boson, with the values obtained from calculations in the simplified model. For each combination of the DM particle and mediator mass values, the width of the mediator is taken into account in the limit calculation. The lines with the hatching represents the excluded regions. The solid grey lines, marked as “Ωh 2 ≥ 0.12”, correspond to parameter regions that reproduce the observed DM relic density in the universe, with the hatched area indicating the region where the DM relic abundance exceeds the observed value.

Limits at 95% confidence level for the masses of the DM particle, which is assumed to be Dirac fermion, and its associated mediator, in a simplified model of DM production via a vector mediator. The parameter exclusion is obtained by comparing the limits on product of the production cross section and the branching fraction for decay to a Z boson, with the values obtained from calculations in the simplified model. For each combination of the DM particle and mediator mass values, the width of the mediator is taken into account in the limit calculation. The lines with the hatching represents the excluded regions. The solid grey lines, marked as “Ωh 2 ≥ 0.12”, correspond to parameter regions that reproduce the observed DM relic density in the universe, with the hatched area indicating the region where the DM relic abundance exceeds the observed value.

Limits at 95% confidence level for the masses of the DM particle, which is assumed to be Dirac fermion, and its associated mediator, in a simplified model of DM production via a vector mediator. The parameter exclusion is obtained by comparing the limits on product of the production cross section and the branching fraction for decay to a Z boson, with the values obtained from calculations in the simplified model. For each combination of the DM particle and mediator mass values, the width of the mediator is taken into account in the limit calculation. The lines with the hatching represents the excluded regions. The solid grey lines, marked as “Ωh 2 ≥ 0.12”, correspond to parameter regions that reproduce the observed DM relic density in the universe, with the hatched area indicating the region where the DM relic abundance exceeds the observed value.

Limits at 95% confidence level for the masses of the DM particle, which is assumed to be Dirac fermion, and its associated mediator, in a simplified model of DM production via a vector mediator. The parameter exclusion is obtained by comparing the limits on product of the production cross section and the branching fraction for decay to a Z boson, with the values obtained from calculations in the simplified model. For each combination of the DM particle and mediator mass values, the width of the mediator is taken into account in the limit calculation. The lines with the hatching represents the excluded regions. The solid grey lines, marked as “Ωh 2 ≥ 0.12”, correspond to parameter regions that reproduce the observed DM relic density in the universe, with the hatched area indicating the region where the DM relic abundance exceeds the observed value.

Limits at 95% confidence level for the masses of the DM particle, which is assumed to be Dirac fermion, and its associated mediator, in a simplified model of DM production via a axial vector. The parameter exclusion is obtained by comparing the limits on product of the production cross section and the branching fraction for decay to a Z boson, with the values obtained from calculations in the simplified model. For each combination of the DM particle and mediator mass values, the width of the mediator is taken into account in the limit calculation. The lines with the hatching represents the excluded regions. The solid grey lines, marked as “Ωh 2 ≥ 0.12”, correspond to parameter regions that reproduce the observed DM relic density in the universe, with the hatched area indicating the region where the DM relic abundance exceeds the observed value.

Limits at 95% confidence level for the masses of the DM particle, which is assumed to be Dirac fermion, and its associated mediator, in a simplified model of DM production via a axial vector. The parameter exclusion is obtained by comparing the limits on product of the production cross section and the branching fraction for decay to a Z boson, with the values obtained from calculations in the simplified model. For each combination of the DM particle and mediator mass values, the width of the mediator is taken into account in the limit calculation. The lines with the hatching represents the excluded regions. The solid grey lines, marked as “Ωh 2 ≥ 0.12”, correspond to parameter regions that reproduce the observed DM relic density in the universe, with the hatched area indicating the region where the DM relic abundance exceeds the observed value.

Limits at 95% confidence level for the masses of the DM particle, which is assumed to be Dirac fermion, and its associated mediator, in a simplified model of DM production via a axial vector. The parameter exclusion is obtained by comparing the limits on product of the production cross section and the branching fraction for decay to a Z boson, with the values obtained from calculations in the simplified model. For each combination of the DM particle and mediator mass values, the width of the mediator is taken into account in the limit calculation. The lines with the hatching represents the excluded regions. The solid grey lines, marked as “Ωh 2 ≥ 0.12”, correspond to parameter regions that reproduce the observed DM relic density in the universe, with the hatched area indicating the region where the DM relic abundance exceeds the observed value.

Limits at 95% confidence level for the masses of the DM particle, which is assumed to be Dirac fermion, and its associated mediator, in a simplified model of DM production via a axial vector. The parameter exclusion is obtained by comparing the limits on product of the production cross section and the branching fraction for decay to a Z boson, with the values obtained from calculations in the simplified model. For each combination of the DM particle and mediator mass values, the width of the mediator is taken into account in the limit calculation. The lines with the hatching represents the excluded regions. The solid grey lines, marked as “Ωh 2 ≥ 0.12”, correspond to parameter regions that reproduce the observed DM relic density in the universe, with the hatched area indicating the region where the DM relic abundance exceeds the observed value.

Exclusion limits at 95% CL on the ultraviolet cutoff for the dielectron channel, with $m_{\ell\ell}>1.8$ TeV ($m_{\ell\ell}>1.9$ TeV for 2016) in the GRW, Hewett, and HLZ conventions for the ADD model. Signal model cross sections are calculated up to LO, and an NNLO correction factor of 1.3 is applied.

Exclusion limits at 95% CL on the ultraviolet cutoff for the dielectron channel, with $m_{\ell\ell}>1.8$ TeV ($m_{\ell\ell}>1.9$ TeV for 2016) in the GRW, Hewett, and HLZ conventions for the ADD model. Signal model cross sections are calculated up to LO, and an NNLO correction factor of 1.3 is applied.

Exclusion limits at 95% CL on the ultraviolet cutoff for the dimuon channel, with $m_{\ell\ell}>1.8$ TeV ($m_{\ell\ell}>1.9$ TeV for 2016) in the GRW, Hewett, and HLZ conventions for the ADD model. Signal model cross sections are calculated up to LO, and an NNLO correction factor of 1.3 is applied.

Exclusion limits at 95% CL on the ultraviolet cutoff for the dimuon channel, with $m_{\ell\ell}>1.8$ TeV ($m_{\ell\ell}>1.9$ TeV for 2016) in the GRW, Hewett, and HLZ conventions for the ADD model. Signal model cross sections are calculated up to LO, and an NNLO correction factor of 1.3 is applied.

Exclusion limits at 95% CL on the ultraviolet cutoff for the combination of the dielectron and dimuon channel, with $m_{\ell\ell}>1.8$ TeV ($m_{\ell\ell}>1.9$ TeV for 2016) in the GRW, Hewett, and HLZ conventions for the ADD model. Signal model cross sections are calculated up to LO, and an NNLO correction factor of 1.3 is applied.

Exclusion limits at 95% CL on the ultraviolet cutoff for the combination of the dielectron and dimuon channel, with $m_{\ell\ell}>1.8$ TeV ($m_{\ell\ell}>1.9$ TeV for 2016) in the GRW, Hewett, and HLZ conventions for the ADD model. Signal model cross sections are calculated up to LO, and an NNLO correction factor of 1.3 is applied.

Dilepton lower exclusion limits at 95% CL on the CI scale ($\Lambda$) for the eight CI models considered, for the dielectron channel. The limits are obtained for $m_{\ell\ell}>400$ GeV.

Dilepton lower exclusion limits at 95% CL on the CI scale ($\Lambda$) for the eight CI models considered, for the dielectron channel. The limits are obtained for $m_{\ell\ell}>400$ GeV.

Dilepton lower exclusion limits at 95% CL on the CI scale ($\Lambda$) for the eight CI models considered, for the dimuon channel. The limits are obtained for $m_{\ell\ell}>400$ GeV.

Dilepton lower exclusion limits at 95% CL on the CI scale ($\Lambda$) for the eight CI models considered, for the dimuon channel. The limits are obtained for $m_{\ell\ell}>400$ GeV.

Dilepton lower exclusion limits at 95% CL on the CI scale ($\Lambda$) for the eight CI models considered, for the combination of the dielectron and dimuon channels. The limits are obtained for $m_{\ell\ell}>400$ GeV.

Dilepton lower exclusion limits at 95% CL on the CI scale ($\Lambda$) for the eight CI models considered, for the combination of the dielectron and dimuon channels. The limits are obtained for $m_{\ell\ell}>400$ GeV.

Ratio of the differential dilepton production cross section in the dimuon and dielectron channels $R_{\mu^{+}\mu^{-}/e^{+}e^{-}}$, as a function of $m_{\ell\ell}$ for events with central leptons. The ratio is obtained after correcting the reconstructed mass spectra to particle level.

Ratio of the differential dilepton production cross section in the dimuon and dielectron channels $R_{\mu^{+}\mu^{-}/e^{+}e^{-}}$, as a function of $m_{\ell\ell}$ for events with central leptons. The ratio is obtained after correcting the reconstructed mass spectra to particle level.

Ratio of the differential dilepton production cross section in the dimuon and dielectron channels $R_{\mu^{+}\mu^{-}/e^{+}e^{-}}$, as a function of $m_{\ell\ell}$ for events with forward leptons. The ratio is obtained after correcting the reconstructed mass spectra to particle level.

Ratio of the differential dilepton production cross section in the dimuon and dielectron channels $R_{\mu^{+}\mu^{-}/e^{+}e^{-}}$, as a function of $m_{\ell\ell}$ for events with forward leptons. The ratio is obtained after correcting the reconstructed mass spectra to particle level.

Invariant mass resolution for dielectron pairs as obtained from simulated DY samples as a function of generated mass.

Invariant mass resolution for dielectron pairs as obtained from simulated DY samples as a function of generated mass.

Invariant mass resolution for dimuon pairs as obtained from simulated DY samples as a function of generated mass.

Invariant mass resolution for dimuon pairs as obtained from simulated DY samples as a function of generated mass.

Observed data and total background estimate in all bins used in the limit calculation for the CI limits. The total event sample is split into the three years of data taking, the two lepton flavors, two categories in $\cos\theta^*$, and two categories in lepton pseudorapidity. Furthermore, each of these categories is split into 6 bins in dilepton invariant mass: [400,500,700,1100,1900,3500,10000] in GeV. The bin names indicate the $\cos\theta^*$ bin (cspos for positive values and csneg for negative values), the lepton flavor (dielectron or dimuon), the year of data taking (2016, 2017, 2018) with associated integrated luminosities of 35.9, 41.5, and 59.4 $\mathrm{fb}^{-1}$ (36.3, 42.1, 61.3 $\mathrm{fb}^{-1}$ in the dimuon channel), respectively. The lepton $|\eta|$ categories are labelled as BB (both leptons with $|\eta| < 1.442$ for electrons and $|\eta| < 1.2$ for muons) and BE (One electron with $|\eta| < 1.442$ and one with $|\eta| > 1.562$ or for muons at least one with $|\eta| > 1.2$). The mass bins are numbered in ascending order.

Observed data and total background estimate in all bins used in the limit calculation for the CI limits. The total event sample is split into the three years of data taking, the two lepton flavors, two categories in $\cos\theta^*$, and two categories in lepton pseudorapidity. Furthermore, each of these categories is split into 6 bins in dilepton invariant mass: [400,500,700,1100,1900,3500,10000] in GeV. The bin names indicate the $\cos\theta^*$ bin (cspos for positive values and csneg for negative values), the lepton flavor (dielectron or dimuon), the year of data taking (2016, 2017, 2018) with associated integrated luminosities of 35.9, 41.5, and 59.4 $\mathrm{fb}^{-1}$ (36.3, 42.1, 61.3 $\mathrm{fb}^{-1}$ in the dimuon channel), respectively. The lepton $|\eta|$ categories are labelled as BB (both leptons with $|\eta| < 1.442$ for electrons and $|\eta| < 1.2$ for muons) and BE (One electron with $|\eta| < 1.442$ and one with $|\eta| > 1.562$ or for muons at least one with $|\eta| > 1.2$). The mass bins are numbered in ascending order.

Observed data and total background estimate in all bins used in the limit calculation for the ADD limits. The total event sample is split into the three years of data taking, the two lepton flavors, two categories in $\cos\theta^*$, and two categories in lepton pseudorapidity. Furthermore, each of these categories is split into 5 bins in dilepton invariant mass: [1800,2200,2600,3000,3400,10000] ([1900,2200,2600,3000,3400,10000] for 2016) in GeV. The bin names indicate the $\cos\theta^*$ bin (cspos for positive values and csneg for negative values), the lepton flavor (dielectron or dimuon), the year of data taking (2016, 2017, 2018) with associated integrated luminosities of 35.9, 41.5, and 59.4 $\mathrm{fb}^{-1}$ (36.3, 42.1, 61.3 $\mathrm{fb}^{-1}$ in the dimuon channel), respectively. The lepton $|\eta|$ categories are labelled as BB (both leptons with $|\eta| < 1.442$ for electrons and $|\eta| < 1.2$ for muons) and BE (One electron with $|\eta| < 1.442$ and one with $|\eta| > 1.562$ or for muons at least one with $|\eta| > 1.2$). The mass bins are numbered in ascending order.

Observed data and total background estimate in all bins used in the limit calculation for the ADD limits. The total event sample is split into the three years of data taking, the two lepton flavors, two categories in $\cos\theta^*$, and two categories in lepton pseudorapidity. Furthermore, each of these categories is split into 5 bins in dilepton invariant mass: [1800,2200,2600,3000,3400,10000] ([1900,2200,2600,3000,3400,10000] for 2016) in GeV. The bin names indicate the $\cos\theta^*$ bin (cspos for positive values and csneg for negative values), the lepton flavor (dielectron or dimuon), the year of data taking (2016, 2017, 2018) with associated integrated luminosities of 35.9, 41.5, and 59.4 $\mathrm{fb}^{-1}$ (36.3, 42.1, 61.3 $\mathrm{fb}^{-1}$ in the dimuon channel), respectively. The lepton $|\eta|$ categories are labelled as BB (both leptons with $|\eta| < 1.442$ for electrons and $|\eta| < 1.2$ for muons) and BE (One electron with $|\eta| < 1.442$ and one with $|\eta| > 1.562$ or for muons at least one with $|\eta| > 1.2$). The mass bins are numbered in ascending order.

Bin-by-bin correlation for the limit calculation for the CI limits. The bin names indicate the $\cos\theta^*$ bin (cspos for positive values and csneg for negative values), the lepton flavor (dielectron or dimuon), the year of data taking (2016, 2017, 2018) with associated integrated luminosities of 35.9, 41.5, and 59.4 $\mathrm{fb}^{-1}$ (36.3, 42.1, 61.3 $\mathrm{fb}^{-1}$ in the dimuon channel), respectively. The lepton $|\eta|$ categories are labelled as BB (both leptons with $|\eta| < 1.442$ for electrons and $|\eta| < 1.2$ for muons) and BE (One electron with $|\eta| < 1.442$ and one with $|\eta| > 1.562$ or for muons at least one with $|\eta| > 1.2$). The mass bins are numbered in ascending order.

Bin-by-bin correlation for the limit calculation for the CI limits. The bin names indicate the $\cos\theta^*$ bin (cspos for positive values and csneg for negative values), the lepton flavor (dielectron or dimuon), the year of data taking (2016, 2017, 2018) with associated integrated luminosities of 35.9, 41.5, and 59.4 $\mathrm{fb}^{-1}$ (36.3, 42.1, 61.3 $\mathrm{fb}^{-1}$ in the dimuon channel), respectively. The lepton $|\eta|$ categories are labelled as BB (both leptons with $|\eta| < 1.442$ for electrons and $|\eta| < 1.2$ for muons) and BE (One electron with $|\eta| < 1.442$ and one with $|\eta| > 1.562$ or for muons at least one with $|\eta| > 1.2$). The mass bins are numbered in ascending order.

Bin-by-bin correlation for the limit calculation for the ADD limits. The bin names indicate the $\cos\theta^*$ bin (cspos for positive values and csneg for negative values), the lepton flavor (dielectron or dimuon), the year of data taking (2016, 2017, 2018) with associated integrated luminosities of 35.9, 41.5, and 59.4 $\mathrm{fb}^{-1}$ (36.3, 42.1, 61.3 $\mathrm{fb}^{-1}$ in the dimuon channel), respectively. The lepton $|\eta|$ categories are labelled as BB (both leptons with $|\eta| < 1.442$ for electrons and $|\eta| < 1.2$ for muons) and BE (One electron with $|\eta| < 1.442$ and one with $|\eta| > 1.562$ or for muons at least one with $|\eta| > 1.2$). The mass bins are numbered in ascending order.

Bin-by-bin correlation for the limit calculation for the ADD limits. The bin names indicate the $\cos\theta^*$ bin (cspos for positive values and csneg for negative values), the lepton flavor (dielectron or dimuon), the year of data taking (2016, 2017, 2018) with associated integrated luminosities of 35.9, 41.5, and 59.4 $\mathrm{fb}^{-1}$ (36.3, 42.1, 61.3 $\mathrm{fb}^{-1}$ in the dimuon channel), respectively. The lepton $|\eta|$ categories are labelled as BB (both leptons with $|\eta| < 1.442$ for electrons and $|\eta| < 1.2$ for muons) and BE (One electron with $|\eta| < 1.442$ and one with $|\eta| > 1.562$ or for muons at least one with $|\eta| > 1.2$). The mass bins are numbered in ascending order.

Integrated Fiducial Higgs cross section. The first uncertainty is the combined statistical uncertainty, the second is the combined systematic uncertainty. As described in the publication, the fiducial volume for 7 and 8 TeV is different than for 13 TeV.

Integrated Fiducial Higgs cross section individually with 2016, 2017 and 2018 dataset and with full Run 2 dataset. The first uncertainty is the combined statistical uncertainty, the second is the combined systematic uncertainty. Results are shown inclusively for all final states.

Integrated Fiducial Higgs cross section individually with 2016, 2017 and 2018 dataset and with full Run 2 dataset. The first uncertainty is the combined statistical uncertainty, the second is the combined systematic uncertainty. Results are shown inclusively for all final states.

Integrated Fiducial Higgs cross section individually with 2016, 2017 and 2018 dataset and with full Run 2 dataset. The first uncertainty is the combined statistical uncertainty, the second is the combined systematic uncertainty. Results are shown inclusively for all final states.

Integrated Fiducial Higgs cross section individually with 2016, 2017 and 2018 dataset and with full Run 2 dataset. The first uncertainty is the combined statistical uncertainty, the second is the combined systematic uncertainty. Results are shown inclusively for all final states.

Higgs fiducial cross section in bins of pT for the 4 leptons. The first uncertainty is statistical, the second is systematic uncertainties. The numbers in this HEP data entry are not divided by the bin width, and therefore the units are in fb.

Higgs fiducial cross section in bins of pT for the 4 leptons. The first uncertainty is statistical, the second is systematic uncertainties. The numbers in this HEP data entry are not divided by the bin width, and therefore the units are in fb.

Higgs fiducial cross section in bins of pT for the 4 leptons. The first uncertainty is statistical, the second is systematic uncertainties. The numbers in this HEP data entry are not divided by the bin width, and therefore the units are in fb.

Higgs fiducial cross section in bins of pT for the 4 leptons. The first uncertainty is statistical, the second is systematic uncertainties. The numbers in this HEP data entry are not divided by the bin width, and therefore the units are in fb.

Higgs fiducial cross section in bins of rapidity for the 4 leptons. The first uncertainty is statistical, the second is systematic uncertainties. The numbers in this HEP data entry are not divided by the bin width, and therefore the units are in fb.

Higgs fiducial cross section in bins of rapidity for the 4 leptons. The first uncertainty is statistical, the second is systematic uncertainties. The numbers in this HEP data entry are not divided by the bin width, and therefore the units are in fb.

Higgs fiducial cross section in bins of rapidity for the 4 leptons. The first uncertainty is statistical, the second is systematic uncertainties. The numbers in this HEP data entry are not divided by the bin width, and therefore the units are in fb.

Higgs fiducial cross section in bins of rapidity for the 4 leptons. The first uncertainty is statistical, the second is systematic uncertainties. The numbers in this HEP data entry are not divided by the bin width, and therefore the units are in fb.

Higgs fiducial cross section in bins of Jet Multiplicity The first uncertainty is statistical, the second is systematic uncertainty.

Higgs fiducial cross section in bins of Jet Multiplicity The first uncertainty is statistical, the second is systematic uncertainty.

Higgs fiducial cross section in bins of Jet Multiplicity The first uncertainty is statistical, the second is systematic uncertainty.

Higgs fiducial cross section in bins of Jet Multiplicity The first uncertainty is statistical, the second is systematic uncertainty.

Higgs fiducial cross section in bins of pT of leading jet. The first uncertainty is statistical, the second is the systematic uncertainty. The numbers in this HEP data entry are not divided by the bin width, and therefore the units are in fb. The 0-jet bin is also included here, which is not shown in the Figure.

Higgs fiducial cross section in bins of pT of leading jet. The first uncertainty is statistical, the second is the systematic uncertainty. The numbers in this HEP data entry are not divided by the bin width, and therefore the units are in fb. The 0-jet bin is also included here, which is not shown in the Figure.

Higgs fiducial cross section in bins of pT of leading jet. The first uncertainty is statistical, the second is the systematic uncertainty. The numbers in this HEP data entry are not divided by the bin width, and therefore the units are in fb. The 0-jet bin is also included here, which is not shown in the Figure.

Higgs fiducial cross section in bins of pT of leading jet. The first uncertainty is statistical, the second is the systematic uncertainty. The numbers in this HEP data entry are not divided by the bin width, and therefore the units are in fb. The 0-jet bin is also included here, which is not shown in the Figure.

Results of likelihood scans for the signal strength modifiers corresponding to the five main SM H boson production mechanisms. The uncertainties, corresponding to one standard deviation confidence intervals, include both statistical and systematic sources. The additional breakdown of the uncertainties into their separate statistical and systematic contributions is also shown.

Results of likelihood scans for the signal strength modifiers corresponding to the five main SM H boson production mechanisms. The uncertainties, corresponding to one standard deviation confidence intervals, include both statistical and systematic sources. The additional breakdown of the uncertainties into their separate statistical and systematic contributions is also shown.

Results of likelihood scans for the signal strength modifiers corresponding to the five main SM H boson production mechanisms. The uncertainties, corresponding to one standard deviation confidence intervals, include both statistical and systematic sources. The additional breakdown of the uncertainties into their separate statistical and systematic contributions is also shown.

Results of likelihood scans for the signal strength modifiers corresponding to the five main SM H boson production mechanisms. The uncertainties, corresponding to one standard deviation confidence intervals, include both statistical and systematic sources. The additional breakdown of the uncertainties into their separate statistical and systematic contributions is also shown.

Result of the 2D likelihood scan for the bosonic and fermionic signal strength modifiers. Contour for the 68% CL region.

Result of the 2D likelihood scan for the bosonic and fermionic signal strength modifiers. Contour for the 68% CL region.

Result of the 2D likelihood scan for the bosonic and fermionic signal strength modifiers. Contour for the 68% CL region.

Result of the 2D likelihood scan for the bosonic and fermionic signal strength modifiers. Contour for the 68% CL region.

Result of the 2D likelihood scan for the bosonic and fermionic signal strength modifiers. Contour for the 95% CL region.

Result of the 2D likelihood scan for the bosonic and fermionic signal strength modifiers. Contour for the 95% CL region.

Result of the 2D likelihood scan for the bosonic and fermionic signal strength modifiers. Contour for the 95% CL region.

Result of the 2D likelihood scan for the bosonic and fermionic signal strength modifiers. Contour for the 95% CL region.

The measured product of cross section times branching fraction for HZZ decay and the SM predictions for the Stage 0 STXS production bins and the inclusive measurement at mH= 125.38 GeV.

The measured product of cross section times branching fraction for HZZ decay and the SM predictions for the Stage 0 STXS production bins and the inclusive measurement at mH= 125.38 GeV.

The measured product of cross section times branching fraction for HZZ decay and the SM predictions for the Stage 0 STXS production bins and the inclusive measurement at mH= 125.38 GeV.

The measured product of cross section times branching fraction for HZZ decay and the SM predictions for the Stage 0 STXS production bins and the inclusive measurement at mH= 125.38 GeV.

The measured product of cross section times branching fraction for HZZ decay and the SM predictions for the merged Stage 1.2 STXS production bins at mH= 125.38 GeV.

The measured product of cross section times branching fraction for HZZ decay and the SM predictions for the merged Stage 1.2 STXS production bins at mH= 125.38 GeV.

The measured product of cross section times branching fraction for HZZ decay and the SM predictions for the merged Stage 1.2 STXS production bins at mH= 125.38 GeV.

The measured product of cross section times branching fraction for HZZ decay and the SM predictions for the merged Stage 1.2 STXS production bins at mH= 125.38 GeV.

Observed correlation matrice between the measured cross sections for the Stage 0 for HZZ decay.

Observed correlation matrice between the measured cross sections for the Stage 0 for HZZ decay.

Observed correlation matrice between the measured cross sections for the Stage 0 for HZZ decay.

Observed correlation matrix between the measured cross sections for the Stage 0 for HZZ decay.

Observed correlation matrice between the measured cross sections for the Stage 1.2 for HZZ decay.

Observed correlation matrice between the measured cross sections for the Stage 1.2 for HZZ decay.

Observed correlation matrice between the measured cross sections for the Stage 1.2 for HZZ decay.

Observed correlation matrix between the measured cross sections for the Stage 1.2 for HZZ decay.

Observed and expected correlations between the parameters in the production mode signal strength fit.

Results of the maximal merging scheme STXS fit. The best fit cross sections are shown together with the respective 68% C.L. intervals. The uncertainty is decomposed into the systematic and statistical components. The expected uncertainties on the fitted parameters are given in brackets. Also listed are the SM predictions for the cross sections and the theoretical uncertainty in those predictions.

The observed and expected impacts from the various sources of systematic uncertainty on the per-production mode signal strengths. The expected impacts are derived using an asimov dataset.

Observed and expected correlations between the parameters in the STXS stage 1.2 maximal merging fit.

Results of the maximal merging scheme STXS fit. The best fit cross sections are shown together with the respective 68% C.L. intervals. The uncertainty is decomposed into the systematic and statistical components. The expected uncertainties on the fitted parameters are given in brackets. Also listed are the SM predictions for the cross sections and the theoretical uncertainty in those predictions.

Results of the minimal merging scheme STXS fit. The best fit cross sections are shown together with the respective 68% C.L. intervals. The uncertainty is decomposed into the systematic and statistical components. The expected uncertainties on the fitted parameters are given in brackets. Also listed are the SM predictions for the cross sections and the theoretical uncertainty in those predictions.

Observed and expected correlations between the parameters in the STXS stage 1.2 maximal merging fit.

Observed and expected correlations between the parameters in the STXS stage 1.2 minimal merging fit.

Results of the minimal merging scheme STXS fit. The best fit cross sections are shown together with the respective 68% C.L. intervals. The uncertainty is decomposed into the systematic and statistical components. The expected uncertainties on the fitted parameters are given in brackets. Also listed are the SM predictions for the cross sections and the theoretical uncertainty in those predictions.

Observed likelihood surface.

Observed and expected correlations between the parameters in the STXS stage 1.2 minimal merging fit.

Expected likelihood surface.

Observed likelihood surface.

Observed likelihood surface.

Expected likelihood surface.

Expected likelihood surface.

Observed likelihood surface.

Results of the STXS stage 0 fit. The best fit cross sections are shown together with the respective 68% C.L. intervals. The uncertainty is decomposed into the systematic and statistical components. The expected uncertainties on the fitted parameters are given in brackets. Also listed are the SM predictions for the cross sections and the theoretical uncertainty in those predictions.

Expected likelihood surface.

Observed and expected correlations between the parameters in the STXS stage 0 fit.

Results of the STXS stage 0 fit. The best fit cross sections are shown together with the respective 68% C.L. intervals. The uncertainty is decomposed into the systematic and statistical components. The expected uncertainties on the fitted parameters are given in brackets. Also listed are the SM predictions for the cross sections and the theoretical uncertainty in those predictions.

Observed likelihood surface.

Observed and expected correlations between the parameters in the STXS stage 0 fit.

Expected likelihood surface.

Observed likelihood surface.

Observed likelihood surface.

Expected likelihood surface.

Expected likelihood surface.

Observed likelihood surface.

Expected likelihood surface.

The Z boson differential cross section as a function of |yZ|.

The Z boson differential cross section as a function of |yZ|.

The Z boson differential cross section as a function of |yZ|.

The Z boson differential cross section as a function of pTZ.

The Z boson differential cross section as a function of pTZ.

The Z boson differential cross section as a function of pTZ.

The TAA-normalized yields of Z bosons as a function of centrality.

The TAA-normalized yields of Z bosons as a function of centrality.

The TAA-normalized yields of Z bosons as a function of centrality.

The product of acceptance and selection efficiency within the fiducial region for the VBF $\mathrm{H}^{\pm\pm}\rightarrow\mathrm{W}^{\pm}\mathrm{W}^{\pm}\rightarrow 2\ell 2v$ and $\mathrm{H}^{\pm}\rightarrow\mathrm{W}^{\pm}\mathrm{Z}\rightarrow 3\ell v$ processes, as a function of $m_{\mathrm{H}_{5}}$. The combination of the statistical and systematic uncertainties is shown. The theoretical uncertainties in the acceptance are also included.

Expected and observed exclusion limits at 95\% confidence level for $\sigma_\mathrm{VBF}(\mathrm{H}^{\pm\pm}) \mathcal{B}(\mathrm{H}^{\pm\pm} \rightarrow \mathrm{W}^{\pm}\mathrm{W}^{\pm})$ as functions of $m_{\mathrm{H}^{\pm\pm}}$. The contribution of the $\mathrm{H}^{\pm}$ boson signal is set to zero for the derivation of the exclusion limits on the $\sigma_\mathrm{VBF}(\mathrm{H}^{\pm\pm}) \mathcal{B}(\mathrm{H}^{\pm\pm} \rightarrow \mathrm{W}^{\pm}\mathrm{W}^{\pm})$. Values above the curves are excluded.

Expected and observed exclusion limits at 95\% confidence level for $\sigma_\mathrm{VBF}(\mathrm{H}^{\pm}) \mathcal{B}(\mathrm{H}^{\pm} \rightarrow \mathrm{W}^{\pm}\mathrm{Z})$ as functions of $m_{\mathrm{H}^{\pm}}$. The contribution of the $\mathrm{H}^{\pm\pm}$ boson signal is set to zero for the derivation of the exclusion limits on the $\sigma_\mathrm{VBF}(\mathrm{H}^{\pm}) \mathcal{B}(\mathrm{H}^{\pm} \rightarrow \mathrm{W}^{\pm}\mathrm{Z})$. Values above the curves are excluded.

Expected and observed exclusion limits at 95\% confidence level for $s_{\mathrm{H}}$ as functions of $m_{\mathrm{H}_{5}}$ in the Georgi-Machacek model. Values above the curves are excluded. The exclusion limits for $s_{\mathrm{H}}$ are shown up to $m_{\mathrm{H}_{5}}=2000$ GeV, given the low sensitivity in the Georgi-Machacek model for values above that mass.

Covariance matrix for all the bins used in Figure 4. The covariance matrix from the background-only fit is shown.

The distribution of Top jet $ au_{3}/ au_{2}$ for the simulation of backgrounds and a 2400 GeV (left-handed) b* signal. The area of the total background contribution and area of the signal component are separately normalized to unity. All analysis selections are applied with the exception of the variable being plotted.

The distribution of W jet $ au_{2}/ au_{1}$ for the simulation of backgrounds and a 2400 GeV (left-handed) b* signal. The area of the total background contribution and area of the signal component are separately normalized to unity. All analysis selections are applied with the exception of the variable being plotted.

The distribution of m_{W} for the simulation of backgrounds and a 2400 GeV (left-handed) b* signal. The area of the total background contribution and area of the signal component are separately normalized to unity. All analysis selections are applied with the exception of the variable being plotted.

The distribution of Top jet DeepCSV score for the simulation of backgrounds and a 2400 GeV (left-handed) b* signal. The area of the total background contribution and area of the signal component are separately normalized to unity. All analysis selections are applied with the exception of the variable being plotted.

The distribution of m_{t} for the simulation of backgrounds and a 2400 GeV (left-handed) b* signal. The area of the total background contribution and area of the signal component are separately normalized to unity. All analysis selections are applied with the exception of the variable being plotted.

Distributions of $m_{tW}$ in the signal region for an interval of $m_{t}$ (65-105 GeV).

Distributions of $m_{tW}$ in the signal region for an interval of $m_{t}$ (105-225 GeV).

Distributions of $m_{tW}$ in the signal region for an interval of $m_{t}$ (225-285 GeV).

Distributions of $m_{tW}$ in themultijet control region of the signal region for an interval of $m_{t}$ (65-105 GeV).

Distributions of $m_{tW}$ in themultijet control region of the signal region for an interval of $m_{t}$ (105-225 GeV).

Distributions of $m_{tW}$ in themultijet control region of the signal region for an interval of $m_{t}$ (225-285 GeV).

Distributions of $m_{t}$ in the ttbar measurement region for an interval of $m_{tar{t}}$ (1200-1300 GeV).

Distributions of $m_{t}$ in the ttbar measurement region for an interval of $m_{tar{t}}$ (1300-1800 GeV).

Distributions of $m_{t}$ in the ttbar measurement region for an interval of $m_{tar{t}}$ (1800-3000 GeV).

Distributions of $m_{t}$ in themultijet control region of the ttbar measurement region for an interval of $m_{tar{t}}$ (1200-1300 GeV).

Distributions of $m_{t}$ in themultijet control region of the ttbar measurement region for an interval of $m_{tar{t}}$ (1300-1800 GeV).

Distributions of $m_{t}$ in themultijet control region of the ttbar measurement region for an interval of $m_{tar{t}}$ (1800-3000 GeV).

Dijet signal efficiencies as a function of the signal mass and lifetime for events satisfying all event and vertex requirements, with corrections based on systematic differences in the vertex reconstruction efficiency between data and simulation.

The distribution of $d_{\mathrm{BV}}$ for $\geq$5-track one-vertex events in data and three simulated multijet signal samples each with a mass of 1600 GeV. The production cross section for each signal model is assumed to be the lower limit excluded by CMS-EXO-17-018, corresponding to values of 0.8, 0.25, and 0.15 fb for the samples with $c\tau =$ 0.3, 1.0, and 10 mm, respectively. The last bin includes the overflow events. This bin includes one event in data with a vertex with large $d_{\mathrm{BV}}$ that appears to arise from tracks originating from separate pp interaction vertices, consistent with background.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data with two 3-track vertices. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data which have exactly one 4-track vertex and one 3-track vertex. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data with two 4-track vertices. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized using $\geq$5-track one-vertex event information. The two vertical red dashed lines separate the regions used in the fit.

Predicted yields for the background-only normalized template, predicted yields for three simulated multijet signals each with a mass of 1600 GeV, and the observed yield in each $d_{\mathrm{VV}}$ bin. The production cross section for each signal model is assumed to be the lower limit excluded by CMS-EXO-17-018, corresponding to values of 0.8, 0.25, and 0.15 fb for samples with $c\tau =$ 0.3, 1.0, and 10 mm, respectively. The uncertainties in the signal yields and the systematic uncertainties in the background prediction reflect the systematic uncertainties given in the text.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the dijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the top squark with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the dijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the top squark with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 300um in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 300um in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 1 mm in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 1 mm in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 10 mm in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 10 mm in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 800 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 800 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 1600 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 1600 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 2400 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals. For $m$ = 2400 GeV, the expected neutralino cross section is $\approx 8\times 10^{-5}$ fb and is not shown.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 2400 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Data-to-simulation efficiency correction factors, shown for multijet and dijet signal topologies in several ranges of $c\tau$. Note that these correction factors account for the two long-lived particles in the simulated events, and are therefore the total correction factors used to scale event yields rather than the correction factors one would apply to individual vertices.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 3-track one-vertex events in data, and is normalized to the number of 3-track two-vertex events in data.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 4-track and 3-track one-vertex events in data, and is normalized to the number of two-vertex events in data which have exactly one 4-track vertex and one 3-track vertex.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 4-track one-vertex events in data, and is normalized to the number of 4-track two-vertex events in data.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from $\geq$5-track one-vertex events in data, and is normalized using one-vertex event information. No $\geq$5-track two-vertex data events pass the selection.

The number of background events obtained from the 1- and 2-leg predictions derived from data, together with the direct observation in data, in bins in ChF, where either the leading or subleading jet has a ChF within the bin edges, and both have a ChF below the upper bin threshold. The bottom panel shows the ratios of the observation in data to the 1-leg and the 2-leg background predictions.

The expected and observed 95% CL upper limits on the production cross section for SIMPs with masses between 1 and 1000 GeV, with the assumption that the SIMP interaction in the detector can be approximated as neutron-like. The theoretical prediction of a simplified model incorporating this approximation and including a scalar mediator with couplings g_chi = -1 and g_q = 1 is also shown (red line). For masses above 100 GeV, where the modelling of the SIMP-nucleon interaction becomes more speculative, the obtained cross section upper limits are increasingly uncertain (shaded area).

Distribution of the transverse momentum of the diphoton system for the $\mathrm{Z}\gamma\gamma$ muon channel. The predicted yields are shown with their pre-fit normalisations. The observed data, the expected signal contribution and the background estimates are presented with error bars showing the corresponding statistical uncertainties.

Distribution of the transverse momentum of the diphoton system, obtained in the control region enriched in misidentified photons for the $\mathrm{W}\gamma\gamma$ electron channel. The predicted yields are shown with their pre-fit normalisations. The observed data, the expected signal contribution and the background estimates are presented with error bars showing the corresponding statistical uncertainties.

Distribution of the transverse momentum of the diphoton system, obtained in the control region enriched in misidentified photons for the $\mathrm{Z}\gamma\gamma$ electron channel. The predicted yields are shown with their pre-fit normalisations. The observed data, the expected signal contribution and the background estimates are presented with error bars showing the corresponding statistical uncertainties.

The measured $\mathrm{W}\gamma\gamma$ and $\mathrm{Z}\gamma\gamma$ fiducial cross sections and corresponding theoretical predictions from MadGraph5_aMC@NLO. The measured yields in the electron and muon channels are extrapolated to a common fiducial phase space determined from simulated signal events at the generated particle level. Generated particles are considered stable if their mean decay length is larger than $1\,\mathrm{cm}$. Generated leptons are required to have a $\mathrm{p}_\mathrm{t}>15\,\mathrm{GeV}$ and $|\eta|<2.5$. The momenta of photons in a cone of $\Delta R=0.1$, the same cone size as the one applied to reconstructed data, are added to the charged lepton momentum to correct for final-state radiation. Generated photons are required to have $\mathrm{p}_\mathrm{t}>15\,\mathrm{GeV}$ and $|\eta|<2.5$. Additionally, the candidate photons are required to have no selected leptons or photons in a cone of radius $\Delta R=0.4$ and no other stable particles, apart from photons and neutrinos, in a cone of radius $\Delta R=0.1$. Events are then selected in the $\mathrm{W}\gamma\gamma$ channel by requiring exactly one electron (muon) with $\mathrm{p}_\mathrm{t}>30\,\mathrm{GeV}$ and at least two photons with $\mathrm{p}_\mathrm{t}>20\,\mathrm{GeV}$. Events are selected in the $\mathrm{Z}\gamma\gamma$ channel by requiring two electrons (muons), at least one of them with $\mathrm{p}_\mathrm{t}>30\,\mathrm{GeV}$, and not less than two photons, each of them with $\mathrm{p}_\mathrm{t}>20\,\mathrm{GeV}$. Additionally, the invariant mass of the dilepton system is required to be $m_{\ell\ell}>55\,\mathrm{GeV}$.

Expected and observed 95% confidence level intervals for the different anomalous couplings parameters in both the $\ell\nu\gamma\gamma$ and $\ell\ell\gamma\gamma$ channels. All parameters are fixed to their SM values except the one that is fitted. No unitarity regularisation scheme is applied. Intervals with a NAN value are not measured in the corresponding channels.

Differential cross sections as function of transverse momentum imbalance between Z boson and dijet for Z+ ≥ 2 jets events.

Differential cross sections as function of transverse momentum imbalance between leading and subleading jet for Z+ ≥ 2 jets events.

Area normalized distribution as function of Delta Phi between Z boson and the leading jet for Z+ ≥ 1 jet events.

Area normalized distribution as function of transverse momentum imbalance between Z boson and the leading jet for Z+ ≥ 1 jet events.

Area normalized distribution as function of Delta Phi between Z boson and dijet for Z+ ≥ 2 jets events.

Area normalized distribution as function of transverse momentum imbalance between Z boson and dijet for Z+ ≥ 2 jets events.

Area normalized distribution as function of transverse momentum imbalance between leading and subleading jet for Z+ ≥ 2 jets events

Response matrix for Delta Phi between Z boson and the leading jet for Z+ ≥ 1 jet events.

Response matrix for transverse momentum imbalance between Z boson and the leading jet for Z+ ≥ 1 jet events.

Response matrix for Delta Phi between Z boson and dijet for Z+ ≥ 2 jets events.

Response matrix for transverse momentum imbalance between Z boson and dijet for Z+ ≥ 2 jets events.

Response matrix for transverse momentum imbalance between the leading jet and subleading jet for Z+ ≥ 2 jets events.

Correlation matrix for Delta Phi between Z boson and the leading jet for Z+ ≥ 1 jet events (for absolute differential cross section measurements).

Correlation matrix for transverse momentum imbalance between Z boson and the leading jet for Z+ ≥ 1 jet events (for absolute differential cross section measurements).

Correlation matrix for Delta Phi between Z boson and dijet for Z+ ≥ 2 jets events (for absolute differential cross section measurements).

Correlation matrix for transverse momentum imbalance between Z boson and dijet for Z+ ≥ 2 jets events (for absolute differential cross section measurements).

Correlation matrix for transverse momentum imbalance between the leading jet and subleading jet for Z+ ≥ 2 jets events (for absolute differential cross section measurements).

Correlation matrix for Delta Phi between Z boson and the leading jet for Z+ ≥ 1 jet events (for normalized differential cross section measurements).

Correlation matrix for transverse momentum imbalance between Z boson and the leading jet for Z+ ≥ 1 jet events (for normalized differential cross section measurements).

Correlation matrix for Delta Phi between Z boson and dijet for Z+ ≥ 2 jets events (for normalized differential cross section measurements).

Correlation matrix for transverse momentum imbalance between Z boson and dijet for Z+ ≥ 2 jets events (for normalized differential cross section measurements).

Correlation matrix for transverse momentum imbalance between the leading jet and subleading jet for Z+ ≥ 2 jets events (for normalized differential cross section measurements).

Expected and observed exclusion limits for heavy higgs boson mass 360 GeV

Expected and observed exclusion limits for heavy higgs boson mass 400 GeV

Expected and observed exclusion limits for heavy higgs boson mass 450 GeV

Expected and observed exclusion limits for heavy higgs boson mass 500 GeV

Expected and observed exclusion limits for heavy higgs boson mass 550 GeV

Expected and observed exclusion limits for heavy higgs boson mass 600 GeV

Expected and observed exclusion limits for heavy higgs boson mass 700 GeV

Expected and observed exclusion limits for heavy higgs boson mass 800 GeV

Expected and observed exclusion limits for heavy higgs boson mass 900 GeV

Expected and observed exclusion limits for heavy higgs boson mass 1000 GeV

Expected and observed exclusion limits for heavy higgs boson mass 1200 GeV

Expected and observed exclusion limits for heavy higgs boson mass 1400 GeV

Expected and observed exclusion limits for heavy higgs boson mass 1600 GeV

Expected and observed exclusion limits for heavy higgs boson mass 1800 GeV

Expected and observed exclusion limits for heavy higgs boson mass 2000 GeV

Expected and observed exclusion limits for heavy higgs boson mass 2500 GeV

Expected and observed exclusion limits for heavy higgs boson mass 3000 GeV

The measured single-differential cross sections in leading jet transverse momenta for the pure electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The measured single-differential cross sections in leading lepton transverse momenta for the pure electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The measured double-differential cross sections as functions of $m_{\mathrm{jj}}$ and $|\Delta\eta_{\mathrm{jj}}|$ for the pure electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The measured single-differential cross sections in photon transverse momenta for the combined QCD-induced and electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The measured single-differential cross sections in leading jet transverse momenta for the combined QCD-induced and electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The measured single-differential cross sections in leading lepton transverse momenta for the combined QCD-induced and electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The measured double-differential cross sections as functions of $m_{\mathrm{jj}}$ and $|\Delta\eta_{\mathrm{jj}}|$ for the combined QCD-induced and electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The expected and observed limits on the aQGC parameters at 95% confidence level. The last column presents the scattering energy values for which the amplitude would violate unitarity for the observed value of the aQGC parameter. All coupling parameter limits are set in TeV$^{-4}$, whereas the unitarity bounds are in TeV.

Distribution of the W boson transverse mass in the WZ 3l signal region. Events from conversions, and DY processes are grouped into the 'Others' category. The vertical error bars represent the statistical uncertainty in the data and the shaded band the uncertainty in the prediction. The signal contributions are scaled to the measured cross sections (postfit).

Results obtained in this analysis and other diboson production cross section measurements at different center-of-mass energies for the CMS, ATLAS, CDF, and D0 Collaborations are presented, and compared with the NNLO QCD x NLO EW and NLO predictions from MATRIX. The vertical error bars represent the uncertainty in the measured cross section.

Distribution of $m_{T}(W)$ in the $N_{jet}\geq 3$ signal region.

Distribution of $M_{3}$ in the $N_{jet}\geq 3$ signal region.

Distribution of $M_{3}$ in the $N_{jet}\geq 3$ signal region.

Distribution of $m(l,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Distribution of $m(l,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Distribution of $\Delta R(l,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Distribution of $\Delta R(l,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Distribution of $\Delta R(j,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Distribution of $\Delta R(j,\gamma)$ in the $N_{jet}\geq 3$ signal region.

Fit result of the multijet template obtained with loosely isolated leptons and the electroweak background to the measured $m_{T}(W)$ distribution with isolated leptons in the $N_{jet}=2$, $N_{b jet}=0$ selection for electrons.

Fit result of the multijet template obtained with loosely isolated leptons and the electroweak background to the measured $m_{T}(W)$ distribution with isolated leptons in the $N_{jet}=2$, $N_{b jet}=0$ selection for electrons.

Fit result of the multijet template obtained with loosely isolated leptons and the electroweak background to the measured $m_{T}(W)$ distribution with isolated leptons in the $N_{jet}=2$, $N_{b jet}=0$ selection for muons.

Fit result of the multijet template obtained with loosely isolated leptons and the electroweak background to the measured $m_{T}(W)$ distribution with isolated leptons in the $N_{jet}=2$, $N_{b jet}=0$ selection for muons.

Distribution of the invariant mass of the lepton and the photon ($m(l,\gamma)$) in the $N_{jet}\geq 3$, $N_{b jet}=0$ selection for the e channel.

Distribution of the invariant mass of the lepton and the photon ($m(l,\gamma)$) in the $N_{jet}\geq 3$, $N_{b jet}=0$ selection for the e channel.

Distribution of the invariant mass of the lepton and the photon ($m(l,\gamma)$) in the $N_{jet}\geq 3$, $N_{b jet}=0$ selection for the $\mu$ channel.

Distribution of the invariant mass of the lepton and the photon ($m(l,\gamma)$) in the $N_{jet}\geq 3$, $N_{b jet}=0$ selection for the $\mu$ channel.

Extracted scale factors for the contribution from misidentified electrons for the three data-taking periods, and the Z$\gamma$, W$\gamma$ simulations.

Extracted scale factors for the contribution from misidentified electrons for the three data-taking periods, and the Z$\gamma$, W$\gamma$ simulations.

Predicted and observed yields in the control regions in the $N_{jet}= 3$ and $\geq 4$ seletions using the post-fit values of the nuisance parameters.

Predicted and observed yields in the control regions in the $N_{jet}= 3$ and $\geq 4$ seletions using the post-fit values of the nuisance parameters.

Predicted and observed yields in the signal regions in the $N_{jet}= 3$ and $\geq 4$ seletions using the post-fit values of the nuisance parameters.

Predicted and observed yields in the signal regions in the $N_{jet}= 3$ and $\geq 4$ seletions using the post-fit values of the nuisance parameters.

The measured inclusive ttgamma cross section in the fiducial phase space compared to the prediction from simulation using Madgraph_aMC@NLO at a center-of-mass energy of 13 TeV.

The measured inclusive ttgamma cross section in the fiducial phase space compared to the prediction from simulation using Madgraph_aMC@NLO at a center-of-mass energy of 13 TeV.

Summary of the measured cross section ratios with respect to the NLO cross section prediction for signal regions binned in the electron channel, muon channel and the combined single lepton measurement.

Summary of the measured cross section ratios with respect to the NLO cross section prediction for signal regions binned in the electron channel, muon channel and the combined single lepton measurement.

The unfolded differential cross sections for $p_{T}(\gamma)$ and the comparison to simulations.

The unfolded differential cross sections for $p_{T}(\gamma)$ and the comparison to simulations.

The unfolded differential cross sections for $|\eta(\gamma)|$ and the comparison to simulations.

The unfolded differential cross sections for $|\eta(\gamma)|$ and the comparison to simulations.

The unfolded differential cross sections for $\Delta R(l,\gamma)$ and the comparison to simulations.

The unfolded differential cross sections for $\Delta R(l,\gamma)$ and the comparison to simulations.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The covariance matrix of systematic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The covariance matrix of statistic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The correlation matrix of statistical uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $p_{T}(\gamma)$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $|\eta(\gamma)|$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

The correlation matrix of systematic uncertainties for the unfolded differential measurement for $\Delta R(l,\gamma)$.

Summary of the one-dimensional intervals at 68 and 95% CL.

Summary of the one-dimensional intervals at 68 and 95% CL.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR3 signal region for the electron channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR3 signal region for the electron channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR3 signal region for the muon channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR3 signal region for the muon channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR4p signal region for the electron channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR4p signal region for the electron channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR4p signal region for the muon channel.

The observed and predicted post-fit yields for the combined Run 2 data set in the SR4p signal region for the muon channel.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional profiled scan for the Wilson coefficient $c_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional profiled scan for the Wilson coefficient $c_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional profiled scan for the Wilson coefficient $c^{I}_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional profiled scan for the Wilson coefficient $c^{I}_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional scan for the Wilson coefficient $c_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional scan for the Wilson coefficient $c_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional scan for the Wilson coefficient $c^{I}_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the one-dimensional scan for the Wilson coefficient $c^{I}_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the two-dimensional scan for the Wilson coefficients $c_{tZ}$ and $c^{I}_{tZ}$.

Negative log-likelihood ratio values with respect to the best fit value of the two-dimensional scan for the Wilson coefficients $c_{tZ}$ and $c^{I}_{tZ}$.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 7 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 7 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The cluster efficiency in bins of hadronic and EM energy in region A. Region A is defined as 391 cm $< r <$ 695.5 cm and 400 cm $< |z| <$ 671 cm. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.

The cluster efficiency in bins of hadronic and EM energy in region A. Region A is defined as 391 cm $< r <$ 695.5 cm and 400 cm $< |z| <$ 671 cm. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.

The cluster efficiency in bins of hadronic and EM energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.

The cluster efficiency in bins of hadronic and EM energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.

The efficiency of $N_{station} > 1$ requirement in bins of hadronic energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bin corresponds to LLPs that decayed leptonically with 0 hadronic energy. The efficiency is calculated with respect to clusters that pass all cluster-level cuts described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 10%.

The efficiency of $N_{station} > 1$ requirement in bins of hadronic energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bin corresponds to LLPs that decayed leptonically with 0 hadronic energy. The efficiency is calculated with respect to clusters that pass all cluster-level cuts described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 10%.

The geometric signal acceptance as a function of $c\tau$. The acceptance is shown for four different LLP mass hypotheses: 7, 15, 40, and 55 GeV. The acceptance is defined by requiring at least 1 LLP to decay in the region defined as 400 cm $< |z| <$ 1100 cm, $r < $ 695.5 cm, and $|\eta| <$ 2.4.

${p}_{T}^{miss}$ distribution of data and MC events in the signal region with the signal stacked on top of the background prediction for a mass hypothesis of ${m}_{stop} = 225GeV$ and ${m}_{LSP} = 50 GeV$. Events from ttW, ttZ, DY, non-prompt leptons, and diboson processes are grouped into the 'Other' category. The lower panel contains the data-to-prediction ratio. The uncertainty band includes statistical, background normalization and all systematic uncertainties.

${p}_{T}^{miss}$ normalized distribution for signal points with different top squark and neutralino masses and SM ttbar events are compared.

${m}_{T2} (e\mu)$ normalized distribution for signal points with different top squark and neutralino masses and SM ttbar events are compared.

$\Delta\eta (e\mu)$ normalized distribution for signal points with different top squark and neutralino masses and SM ttbar events are compared.

$\Delta\phi (e\mu)$ normalized distribution for signal points with different top squark and neutralino masses and SM ttbar events are compared.

Normalized DNN score distribution in the signal region for two mass hypotheses to compare signal and ttbar background: $m_{LSP}=$ 50 (100) GeV and $m_{stop}=$ 225 (275) GeV.

Post-fit DNN score distributions in the signal region with the signal superimposed and scaled by a factor of 20 for different mass hypotheses of $({m}_{stop}, {m}_{LSP})=$ (225, 50) GeV. The uncertainty band includes statistical, background normalization, and all systematic uncertainties. Events from ttW, ttZ, DY, non-prompt leptons, and diboson processes are grouped into the 'Other' category. The lower panel contains the data-to-prediction ratio before the fit (line) and after (dots), each of them with their corresponding band of uncertainties, and the ratio between the sum of the signal and background predictions and the background prediction. The masses of the signal model correspond to the values of the DNN mass parameters in each distribution.

Post-fit DNN score distributions in the signal region with the signal superimposed and scaled by a factor of 20 for different mass hypotheses of $({m}_{stop}, {m}_{LSP})=$ (275, 100) GeV. The uncertainty band includes statistical, background normalization, and all systematic uncertainties. Events from ttW, ttZ, DY, non-prompt leptons, and diboson processes are grouped into the 'Other' category. The lower panel contains the data-to-prediction ratio before the fit (line) and after (dots), each of them with their corresponding band of uncertainties, and the ratio between the sum of the signal and background predictions and the background prediction. The masses of the signal model correspond to the values of the DNN mass parameters in each distribution.

Post-fit DNN score distributions in the signal region with the signal superimposed and scaled by a factor of 20 for different mass hypotheses of $({m}_{stop}, {m}_{LSP})=$ (275, 70) GeV. The uncertainty band includes statistical, background normalization, and all systematic uncertainties. Events from ttW, ttZ, DY, non-prompt leptons, and diboson processes are grouped into the 'Other' category. The lower panel contains the data-to-prediction ratio before the fit (line) and after (dots), each of them with their corresponding band of uncertainties, and the ratio between the sum of the signal and background predictions and the background prediction. The masses of the signal model correspond to the values of the DNN mass parameters in each distribution.

Post-fit DNN score distributions in the signal region with the signal superimposed and scaled by a factor of 20 for different mass hypotheses of $({m}_{stop}, {m}_{LSP})=$ (245, 100) GeV. The uncertainty band includes statistical, background normalization, and all systematic uncertainties. Events from ttW, ttZ, DY, non-prompt leptons, and diboson processes are grouped into the 'Other' category. The lower panel contains the data-to-prediction ratio before the fit (line) and after (dots), each of them with their corresponding band of uncertainties, and the ratio between the sum of the signal and background predictions and the background prediction. The masses of the signal model correspond to the values of the DNN mass parameters in each distribution.

Upper limit at 95% CL on the signal cross section as a function of the top squark and neutralino masses. The model is excluded for all of the colored region inside the black boundary.

Upper limit at 95% CL on the signal cross section as a function of the top squark mass for mass difference of 175 GeV. The green and yellow bands represent the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The purple dotted line indicates the approximate NNLO + NNLL productioncross section.

Upper limit at 95% CL on the signal cross section as a function of the top squark mass for mass difference of 185 GeV. The green and yellow bands represent the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The purple dotted line indicates the approximate NNLO + NNLL productioncross section.

Upper limit at 95% CL on the signal cross section as a function of the top squark mass for mass difference of 165 GeV. The green and yellow bands represent the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The purple dotted line indicates the approximate NNLO + NNLL productioncross section.

Summary of the contributions of the experimental uncertainties in the DNN score distribution for signal and the ttbar background. The values represent the relative variation in the number of expected events across different signal models in the signal region.

Summary of the contribution of each modeling uncertainty source to the DNN score distribution for the ttbar background.

95% CL expected upper limit on top squark cross section for the direct decay model.

95% CL observed upper limit on top squark cross section for the direct decay model.

Contour of the region excluded with 95% CL (expected).

Contour of the region excluded with 95% CL (observed).

95% CL expected upper limit on top squark cross section for the cascade decay model.

95% CL observed upper limit on top squark cross section for the cascade decay model.

Contour of the region excluded with 95% CL (expected).

Contour of the region excluded with 95% CL (observed).

95% CL expected upper limit on top squark cross section for the mixed decay model.

95% CL observed upper limit on top squark cross section for the mixed decay model.

Contour of the region excluded with 95% CL (expected).

Contour of the region excluded with 95% CL (observed).

95% CL expected limit on $\sigma/\sigma_{theory}$ for a fermionic DM particle with m = 1 GeV, as a function of the mediator mass for a scalar.

95% CL observed limit on $\sigma/\sigma_{theory}$ for a fermionic DM particle with m = 1 GeV, as a function of the mediator mass for a scalar.

95% CL expected limit on $\sigma/\sigma_{theory}$ for a fermionic DM particle with m = 1 GeV, as a function of the mediator mass for a pseudoscalar.

95% CL observed limit on $\sigma/\sigma_{theory}$ for a fermionic DM particle with m = 1 GeV, as a function of the mediator mass for a pseudoscalar.

The correlation matrix for the Higgs pT measurements, both for the unregularized and regularized fits. The last bin is inclusive.

The correlation matrix for the jet multiplicity measurements, both for the unregularized and regularized fits. The last bin is inclusive.

The correlation matrix for the leading jet pT measurements, both for the unregularized and regularized fits. The last bin is inclusive.

The fiducial integrated signal strength and cross section.