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.