Dependence of the parameter $A$ on $(t_1, t_2)$.

Dependence of the parameter $R$ on $(t_1, t_2)$.

Dependence of the parameter $c$ on $(t_1, t_2)$.

Dependence of the parameter $R$ on the two scattered proton transverse momenta.

Correlation coefficients among values of best parameters for the two-channel model, in the case of the exponential parametrisation of the proton-pomeron form factor.

Correlation coefficients among values of best parameters for the two-channel model, in the case of the Orear-type parametrisation of the proton-pomeron form factor.

Correlation coefficients among values of best parameters for the two-channel model, in the case of the power-law parametrisation of the proton-pomeron form factor.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and a $\tilde{\chi}_1^\pm$ that decays to $W^\pm\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction to $q_i\bar{q}_j W\pm\tilde{\chi}_1^0$.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to bottom quarks ($\tilde{g}\to b\bar{b}\tilde{\chi}_1^0$). Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to top quarks ($\tilde{g}\to t\bar{t}\tilde{\chi}_1^0$). Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for light-flavor squark pair production, where the squarks decay to $q\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay and 1-fold degeneracy in the light-flavor squarks (corresponding to the inner set of curves in the limit plot). To get the theory cross section for other N-fold degeneracy assumptions (e.g. 8-fold for the outer curves in the limit plot), just multiply by N.

Exclusion limits at 95% CL for bottom squark pair production, where the squarks decay to $b\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay to $t\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay to $b\tilde{\chi}_1^\pm$ and the $\tilde{\chi}_1^0$ decay to $W^\pm\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay either to $b\tilde{\chi}_1^\pm\to bW^\pm\tilde{\chi}_1^0$ or to $t\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for top squark pair production, where the squarks decay to $c\tilde{\chi}_1^0$. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, assuming unity branching fraction for the given decay.

Exclusion limits at 95% CL for the mono-$\phi$ model, in which a resonantly-produced colored scalar decays to a massive Dirac fermion and a quark. Signal cross sections are calculated at leading order in $\alpha_S$, assuming unity branching fraction for the given decay.

Cross section limits for $\mathrm{LQ}\to\mathrm{q}\nu$, where $q=u,\,d,\,s,\,\mathrm{or}\,c$. Limits are at the 95% confidence level. Theory cross sections are LO for vector LQ, and NLO for scalar LQ. Branching ratio is assumed to be 100% to $\mathrm{q}\nu$.

Cross section limits for $\mathrm{LQ}\to\mathrm{b}\nu$. Limits are at the 95% confidence level. Theory cross sections are LO for vector LQ, and NLO for scalar LQ. Branching ratio is assumed to be 100% to $\mathrm{b}\nu$.

Cross section limits for $\mathrm{LQ}\to\mathrm{t}\nu$. Limits are at the 95% confidence level. Theory cross sections are LO for vector LQ, and NLO for scalar LQ. Branching ratios are assumed to be $\mathcal{B}(\mathrm{LQ}\to\mathrm{t}\nu)=1-\beta$, and $\mathcal{B}(\mathrm{LQ}\to\mathrm{b}\tau)=\beta$.

Predictions and observations for monojet signal regions

Predictions and observations for signal regions with $250 \leq H_\mathrm{T} < 450$ GeV

Predictions and observations for signal regions with $450 \leq H_\mathrm{T} < 575$ GeV and $N_\mathrm{j}<7$

Predictions and observations for signal regions with $450 \leq H_\mathrm{T} < 575$ GeV and $N_\mathrm{j}\geq7$

Predictions and observations for signal regions with $575 \leq H_\mathrm{T} < 1200$ GeV and $N_\mathrm{j}^\mathrm{hi}<4$

Predictions and observations for signal regions with $575 \leq H_\mathrm{T} < 1200$ GeV and $4\leq N_\mathrm{j}^\mathrm{hi}<7$

Predictions and observations for signal regions with $575 \leq H_\mathrm{T} < 1200$ GeV and $N_\mathrm{j}\geq7$

Predictions and observations for signal regions with $1200 \leq H_\mathrm{T} < 1500$ GeV and $N_\mathrm{j}^\mathrm{hi}<4$

Predictions and observations for signal regions with $1200 \leq H_\mathrm{T} < 1500$ GeV and $4\leq N_\mathrm{j}^\mathrm{hi}<7$

Predictions and observations for signal regions with $1200 \leq H_\mathrm{T} < 1500$ GeV and $N_\mathrm{j}\geq7$

Predictions and observations for signal regions with $H_\mathrm{T} \geq 1500$ GeV and $N_\mathrm{j}<7$

Predictions and observations for signal regions with $H_\mathrm{T} \geq 1500$ GeV and $N_\mathrm{j}\geq7$

Covariance matrix for the 282 signal regions of the inclusive $M_\mathrm{T2}$ search

Correlation matrix for the 282 signal regions of the inclusive $M_\mathrm{T2}$ search

Bin number definitions for the $M_\mathrm{T2}$ covariance and correlation matrices

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 10$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 50$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 200$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct light squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 10$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark.

Exclusion limits at 95% CL for direct light squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 50$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark.

Exclusion limits at 95% CL for direct light squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino is long-lived with $c\tau_0 = 200$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino is long-lived with $c\tau_0 = 10$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino is long-lived with $c\tau_0 = 50$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino is long-lived with $c\tau_0 = 200$ cm and mass O(100) MeV greater than the neutralino's mass. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

The maximum chargino mass excluded at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the gluino mass is fixed to 1900 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

The maximum chargino mass excluded at 95% CL for direct squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the squark mass is fixed to 900 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, for a single light squark. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

The maximum chargino mass excluded at 95% CL for direct squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the squark mass is fixed to 1500 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, and the eight light squarks' masses are assumed to be degenerate. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 1$ to 2000 cm while the stop mass is fixed to 1000 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$. If all kinematically allowed chargino masses are excluded, the curves, including 68 and 95% expected, tend to overlap. At short decay lengths, horizontal exclusion lines are obtained from the inclusive analysis, as this is not affected by track reconstruction inefficiencies, which may arise when the chargino decays before the CMS tracker, and therefore shows better sensitivity to scenarios with very small lifetime compared to the disappearing track search, based on median expected limits.

Exclusion limits at 95% CL for direct gluino pair production, where the gluinos decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 5$ to 1000 cm while the gluino mass is fixed to 1600 GeV and the neutralino's mass is fixed to 1575 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Exclusion limits at 95% CL for direct squark pair production, where the squarks decay to light-flavor quarks and either the lightest neutralino, or the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 5$ to 1000 cm while the squark mass is fixed to 2000 GeV and the neutralino's mass is fixed to 1000 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$, and the eight light squarks' masses are assumed to be degenerate.

Exclusion limits at 95% CL for direct stop pair production, where the stops decay to either a top and the lightest neutralino, or a bottom and the lightest chargino, and the chargino mass is O(100) MeV greater than the neutralino's mass. The chargino's lifetime is varied from $c\tau_{0} = 5$ to 1000 cm while the stop mass is fixed to 1100 GeV and the neutralino's mass is fixed to 1000 GeV. Signal cross sections are calculated at approximately NNLO+NNLL order in $\alpha_S$.

Predictions and observations for 2016 disappearing track signal regions

Predictions and observations for 2017-2018 pixel track signal regions

Predictions and observations for 2017-2018 medium (M) and long (L) length track signal regions

Covariance matrix for the 68 signal regions of the disappearing tracks $M_\mathrm{T2}$ search

Correlation matrix for the 68 signal regions of the disappearing tracks $M_\mathrm{T2}$ search

Experimental Xi values and FJF predictions for the four NRQCD terms using BK LDME parameters

Experimental Xi values and FJF predictions for the four NRQCD terms using BCKL LDME parameters

Experimental Xi values and FJF predictions for the four NRQCD terms using BK LDME parameters

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.

Comparison of the $m_{miss}$ shapes for the simulated signal events within the fiducial region and those outside it, after including the effect of PU protons as describe in the text, for a generated $m_{X}$ mass of 1000 GeV. The distributions are shown for single(+z)-single(−z) proton reconstruction categories.

Product of the acceptance times the combined reconstruction and identification efficiency, as a function of $m_{X}$ , for events generated inside the fiducial volume defined in Table 1. The curves shown panel display the different final states. The definition of the fiducial region and signal model used to estimate the acceptance is provided in the text.

Product of the acceptance times the combined reconstruction and identification efficiency, as a function of $m_{X}$ , for different proton reconstruction categories in the $Z \rightarrow \mu\mu $ analysis, and for events generated inside the fiducial volume defined in Table 1. The curves shown panel display the different final states. The definition of the fiducial region and signal model used to estimate the acceptance is provided in the text.

Distributions of the reconstructed proton $\xi$ in the negative arm from the proton mixing procedure with simulated MC events, are compared to data. The ee final states are shown.The ee events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the positive arm from the proton mixing procedure with simulated MC events, are compared to data. The ee final states are shown.The ee events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed di-proton rapidity from the proton mixing procedure with simulated MC events, are compared to data. The ee final states are shown.The ee events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the negative arm from the proton mixing procedure with simulated MC events, are compared to data. The $\mu\mu$ final states are shown.The $\mu\mu$ events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the positive arm from the proton mixing procedure with simulated MC events, are compared to data. The $\mu\mu$ final states are shown.The $\mu\mu$ events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed di-proton rapidity from the proton mixing procedure with simulated MC events, are compared to data. The $\mu\mu$ final states are shown.The $\mu\mu$ events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the negative arm from the proton mixing procedure with simulated MC events, are compared to data. The $\gamma$ final states are shown.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the positive arm from the proton mixing procedure with simulated MC events, are compared to data. The $\gamma$ final states are shown.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed di-proton rapidity from the proton mixing procedure with simulated MC events, are compared to data. The $\gamma$ final states are shown.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Validation of the background modelling method, using the $e\mu$ control sample. Selected $e\mu$ events are mixed with protons from $Z\rightarrow \mu \mu$ events with $p_{T}(Z)<10$ GeV to simulate the combinatorial background shape, while the data points are unaltered $e\mu$ events. The proton $\xi$ distribution for the positive CT-PPS arm is shown.

Validation of the background modelling method, using the $e\mu$ control sample. Selected $e\mu$ events are mixed with protons from $Z\rightarrow \mu \mu$ events with $p_{T}(Z)<10$ GeV to simulate the combinatorial background shape, while the data points are unaltered $e\mu$ events. The proton $\xi$ distribution for the negative CT-PPS arm is shown.

Validation of the background modelling method, using the $e\mu$ control sample. Selected $e\mu$ events are mixed with protons from $Z\rightarrow \mu \mu$ events with $p_{T}(Z)<10$ GeV to simulate the combinatorial background shape, while the data points are unaltered $e\mu$ events. The di-proton invariant mass is shown.

Validation of the background modelling method, using the $e\mu$ control sample. Selected $e\mu$ events are mixed with protons from $Z\rightarrow \mu \mu$ events with $p_{T}(Z)<10$ GeV to simulate the combinatorial background shape, while the data points are unaltered $e\mu$ events. The $m_{miss}$ is shown.

$m_{miss}$ distributions in the $pp\rightarrow ppZeeX$, with the protons reconstructed with the multi(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZeeX$, with the protons reconstructed with the multi(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZeeX$, with the protons reconstructed with the single(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZeeX$, with the protons reconstructed with the single(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZ\mu\mu X$, with the protons reconstructed with the multi(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZ\mu\mu X$, with the protons reconstructed with the multi(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZ\mu\mu X$, with the protons reconstructed with the single(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZ\mu\mu X$, with the protons reconstructed with the single(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow pp\gamma X$, with the protons reconstructed with the multi(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 10 pb.

$m_{miss}$ distributions in the $pp\rightarrow pp\gamma X$, with the protons reconstructed with the multi(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 10 pb.

$m_{miss}$ distributions in the $pp\rightarrow pp\gamma X$, with the protons reconstructed with the single(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 10 pb.

$m_{miss}$ distributions in the $pp\rightarrow pp\gamma X$, with the protons reconstructed with the single(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 10 pb.

Upper limits on the $pp\rightarrow ppZ/\gamma X$ cross section at 95% CL, as a function of $m_X$. The 68 and 95% central intervals of the expected limits are represented by the green and yellow bands, respectively, while the observed limit is superimposed as a curve. The plot corresponds to the $Z\rightarrow ee$ analysis.

Upper limits on the $pp\rightarrow ppZ/\gamma X$ cross section at 95% CL, as a function of $m_X$. The 68 and 95% central intervals of the expected limits are represented by the green and yellow bands, respectively, while the observed limit is superimposed as a curve. The plot corresponds to the $Z\rightarrow \mu\mu$ analysis.

Upper limits on the $pp\rightarrow ppZ/\gamma X$ cross section at 95% CL, as a function of $m_X$. The 68 and 95% central intervals of the expected limits are represented by the green and yellow bands, respectively, while the observed limit is superimposed as a curve. The plot corresponds to the $Z$ analysis.

Upper limits on the $pp\rightarrow ppZ/\gamma X$ cross section at 95% CL, as a function of $m_X$. The 68 and 95% central intervals of the expected limits are represented by the green and yellow bands, respectively, while the observed limit is superimposed as a curve. The plot corresponds to the $\gamma$ analysis.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of Z pt on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of Z pt on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of |y(Z)| on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of |y(Z)| on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of $\phi^{\scriptscriptstyle *}_\eta$ on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the absolute cross section measurements in bins of $\phi^{\scriptscriptstyle *}_\eta$ on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of Z pt on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of Z pt on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of |y(Z)| on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of |y(Z)| on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of $\phi^{\scriptscriptstyle *}_\eta$ on dimuons.

The relative statistical and systematic uncertainties in % from various sources for the normalized cross section measurements in bins of $\phi^{\scriptscriptstyle *}_\eta$ on dielectrons.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dimuon final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dimuon final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dimuon final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dimuon final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dimuon final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dielectron final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dielectron final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dielectron final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dielectron final state.

The relative statistical and systematic uncertainties in % from various sources for the absolute double-differential cross section measurements in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dielectron final state.

The measured absolute cross sections in bins of |y(Z)|, using dressed level leptons. The cross sections are normalized by the bin width.

The measured absolute cross sections in bins of Z pt, using dressed level leptons. The cross sections are normalized by the bin width.

The measured absolute cross sections in bins of $\phi^{\scriptscriptstyle *}_\eta$, using dressed level leptons. The cross sections are normalized by the bin width.

The measured absolute cross sections in bins of Z pt different |y(Z)| bins, using dressed level leptons. The cross sections are normalized by the bin width.

The measured normalized cross sections in bins of Z pt, using dressed level leptons. The cross sections are divided by the bin width and then the distributions are normalized to 1.

The measured normalized cross sections in bins of $\phi^{\scriptscriptstyle *}_\eta$, using dressed level leptons. The cross sections are divided by the bin width and then the distributions are normalized to 1.

The measured normalized cross sections in bins of |y(Z)|, using dressed level leptons. The cross sections are divided by the bin width and then the distributions are normalized to 1.

The measured normalized cross sections in bins of Z pt in |y(Z)| bins, using dressed level leptons. The cross sections are divided by the bin width and then the distributions are normalized to 1.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of |y(Z)| in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dimuon final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of |y(Z)| in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dielectron final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of |y(Z)| in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dimuon final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of |y(Z)| in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dielectron final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of |y(Z)| in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dimuon final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of |y(Z)| in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dielectron final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of |y(Z)| in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dimuon final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of |y(Z)| in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dielectron final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dielectron final state.

The measured absolute cross sections in bins of |y(Z)|, using born level leptons. The cross sections are normalized by the bin width.

The measured absolute cross sections in bins of Z pt, using born level leptons. The cross sections are normalized by the bin width. The first bin (0-1 GeV) is not shown as large differences were observed in aMC@NLO and POWHEG predictions at Born level.

The measured absolute cross sections in bins of $\phi^{\scriptscriptstyle *}_\eta$, using born level leptons. The cross sections are normalized by the bin width.

The measured absolute cross sections in bins of Z pt different |y(Z)| bins, using born level leptons. The cross sections are normalized by the bin width. The first bin (0-1 GeV) is not shown as large differences were observed in aMC@NLO and POWHEG predictions at Born level.

The measured normalized cross sections in bins of Z pt, using born level leptons. The first bin (0-1 GeV) is not shown as large differences were observed in aMC@NLO and POWHEG predictions at Born level. The cross sections are divided by the bin width and then the distributions are normalized to 1.

The measured normalized cross sections in bins of $\phi^{\scriptscriptstyle *}_\eta$, using born level leptons. The cross sections are divided by the bin width and then the distributions are normalized to 1.

The measured normalized cross sections in bins of |y(Z)|, using born level leptons.

The measured normalized cross sections (left) in bins of Z pt in |y(Z)| bins, using born level leptons. The first bin (0-1 GeV) is not shown as large differences were observed in aMC@NLO and POWHEG predictions at Born level. The cross sections are divided by the bin width and then the distributions are normalized to 1.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of |y(Z)| in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dilepton final state.

Covariance matrix using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of |y(Z)| in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dilepton final state.

Covariance matrix using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of |y(Z)| in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dilepton final state.

Covariance matrix for normalized cross sections using dressed level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of |y(Z)| in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of $\phi^{\scriptscriptstyle *}_\eta$ in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0 < |y(Z)| < 0.4 bin in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.4 < |y(Z)| < 0.8 bin in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 0.8 < |y(Z)| < 1.2 bin in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.2 < |y(Z)| < 1.6 bin in the dilepton final state.

Covariance matrix for normalized cross sections using born level leptons for all bins used in bins of Z pt for the 1.6 < |y(Z)| < 2.4 bin in the dilepton final state.

Reconstruction level differential cross section spectla, obtained for the central acceptance, $|\eta| < 3$, for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology compared to the to the EPOS-LHC predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV

The number of high purity tracksin the first $\eta$ bin after a gap of $4.5 \leq \Delta\eta^F < 5.0$ for events with the Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV

The transeverse momentum of high purity tracksin the first $\eta$ bin after a gap of $4.5 \leq \Delta\eta^F < 5.0$ for events with the Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV

The total energy of all Particle Flow candidatesin the first $\eta$ bin after a gap of $4.5 \leq \Delta\eta^F < 5.0$ for events with the Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV

Unfolded diffraction enhanced differential cross section for events with Pomeron-Lead ($\mathrm{I\!P}\mathrm{Pb}$) topology, defined on the region $|\eta| < 5.2$, compared to hadron level predictions of the MC generators. The data are corrected for the contribution from events with undetectable energy in the HF calorimeter adjacent to the rapidity gap. The corrections are obtained using the EPOS-LHC MC sample. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffraction enhanced differential cross section for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology, defined on the region $|\eta| < 5.2$, compared to hadron level predictions of the MC generators. The data are corrected for the contribution from events with undetectable energy in the HF calorimeter adjacent to the rapidity gap. The corrections are obtained using the EPOS-LHC MC sample. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffractive enhanced differential cross section spectla for events with Pomeron-Lead ($\mathrm{I\!P}\mathrm{Pb}$) topology, compared to the to the hadron-level EPOS-LHC predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffractive enhanced differential cross section spectla for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology, compared to the to the hadron-level EPOS-LHC predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffractive enhanced differential cross section spectla for events with Pomeron-Lead ($\mathrm{I\!P}\mathrm{Pb}$) topology, compared to the to the hadron-level QGSJET II predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffractive enhanced differential cross section spectla for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology, compared to the to the hadron-level QGSJET II predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Reconstruction level diffraction enhanced differential cross section spectrum for Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology corrected for the contribution from events with undetectable energy in the HF calorimeter adjacent to the rapidity gap, compared with the spectrum obtained with events satisfying the ZDC veto requirement $E_{ZDC−} < 1 \mathrm{~TeV}$ which selects only the events without lead nuclear break up (without correction for HF undetectable energy). Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Elliptic-flow coefficients $v_2$ based on the six-particle correlations technique, as functions of transverse momentum and in bins of centrality. The results correspond to the range $|\eta| < 2.4$.

Elliptic-flow coefficients $v_2$ based on the eight-particle correlations technique, as functions of transverse momentum and in bins of centrality. The results correspond to the range $|\eta| < 2.4$.

Triangular-flow coefficients $v_3$ based on the two-particle correlations technique, as functions of transverse momentum and in bins of centrality. The results correspond to the range $|\eta| < 2.4$.

Triangular-flow coefficients $v_3$ based on the scalar-product technique, as functions of transverse momentum and in bins of centrality. The results correspond to the range $|\eta| < 0.8$.

Triangular-flow coefficients $v_3$ based on the four-particle correlations technique, as functions of transverse momentum and in bins of centrality. The results correspond to the range $|\eta| < 2.4$.

The $v_4$ coefficients based on the two-particle correlations technique, as functions of transverse momentum and in bins of centrality. The results correspond to the range $|\eta| < 2.4$.

The $v_4$ coefficients based on the scalar-product technique, as functions of transverse momentum and in bins of centrality. The results correspond to the range $|\eta| < 0.8$.

Centrality dependence of the spectrum-weighted $v_2$ flow harmonics with $0.3 < p_{\mathrm{T}} < 3.0~\mathrm{GeV}/c$. The $v_2$ results are shown for two-, four-, six-, and eight-particle correlations.

Centrality dependence of the spectrum-weighted $v_3$ flow harmonics with $0.3 < p_{\mathrm{T}} < 3.0~\mathrm{GeV}/c$. The results are shown for two- and four-particle correlations.

Centrality dependence of the spectrum-weighted $v_4$ flow harmonics with $0.3 < p_{\mathrm{T}} < 3.0~\mathrm{GeV}/c$. The results are shown for two-particle correlations.

Centrality dependence of $v_2\{4\}/v_2\{2\}$ ratios.

Centrality dependence of $v_2\{6\}/v_2\{4\}$ ratios.

Centrality dependence of $v_3\{4\}/v_3\{2\}$ ratios.

The $v_2$ results measured with two-particle correlations from PbPb collisions at $5.02~$TeV, shown as a function of $p_{\mathrm{T}}$ in eleven centrality bins.

The $v_3$ results measured with two-particle correlations from PbPb collisions at $5.02~$TeV, shown as a function of $p_{\mathrm{T}}$ in eleven centrality bins.

The $v_4$ results measured with two-particle correlations from PbPb collisions at $5.02~$TeV, shown as a function of $p_{\mathrm{T}}$ in eleven centrality bins.

Ratios of the $v_2$ harmonic coefficients from two-particle correlations in XeXe and PbPb collisions as functions of $p_{\mathrm{T}}$ in 11 centrality bins.

Ratios of the $v_3$ harmonic coefficients from two-particle correlations in XeXe and PbPb collisions as functions of $p_{\mathrm{T}}$ in 11 centrality bins.

Ratios of the $v_4$ harmonic coefficients from two-particle correlations in XeXe and PbPb collisions as functions of $p_{\mathrm{T}}$ in 11 centrality bins.

Centrality dependence of the spectrum-weighted $v_2$, $v_3$, and $v_4$ harmonic coefficients from two-particle correlations method for $0.3 < p_{\mathrm{T}} < 3.0 \mathrm{GeV}/c$ for PbPb collisions at $5.02$~TeV.

Ratios of the $v_2$, $v_3$, and $v_4$ harmonic coefficients from two-particle correlations in XeXe and PbPb collisions as functions or $0.3 < p_{\mathrm{T}} < 3.0~\mathrm{GeV}/c$ as a function of centrality.

Expected and observed upper limits on the AQGC operators $a^Z_C/\Lambda^2$, with no unitarization. The $y$ axis shows the limit on the ratio of the observed cross section to the cross section predicted for each anomalous coupling value ($\sigma_\mathrm{AQGC}$).

Limits on LEP-like dimension-6 anomalous quartic gauge coupling parameters, with and without unitarization via a clipping procedure.

Conversion of limits on $a^W_0$ to dimension-8 $f_{M,i}$ operators, using the assumption of vanishing $WWZ\gamma$ couplings to eliminate some parameters. When quoting limits on one of the operators, the other is fixed to zero. The results for $|f_{M,0}/\Lambda^{4}|$ and $|f_{M,4}/\Lambda^{4}|$ are shown with and without clipping of the signal model at 1.4 TeV, when the other parameter is fixed to the SM value of zero.

Conversion of limits on $a^W_0$ and $a^W_C$ to dimension-8 $f_{M,i}$ operators, using the assumption that all $f_{M,i}$ except one are equal to zero. The results are shown with and without clipping of the signal model at 1.4 TeV.

{Expected and observed limits in the two-dimensional plane of $a^W_0/\Lambda^2$ vs. $a^W_C/\Lambda^2$. The limits are described by analytical ellipses of equation $(x-x0)^2/a^2 + (y-y0)^2/b^2 = 1$ and rotated counter-clockwise by $\theta$ degrees, where $x$ and $y$ in the equation correspond to the $a_0^W$ and $a_C^W$ couplings, respectively.

{Expected and observed limits in the two-dimensional plane of $a^W_0/\Lambda^2$ vs. $a^W_C/\Lambda^2$ with unitarization imposed by clipping the signal model at 1.4 TeV. The limits are described by analytical ellipses of equation $(x-x0)^2/a^2 + (y-y0)^2/b^2 = 1$ and rotated counter-clockwise by $\theta$ degrees, where $x$ and $y$ in the equation correspond to the $a_0^W$ and $a_C^W$ couplings, respectively.

{Expected and observed limits in the two-dimensional plane of $a^Z_0/\Lambda^2$ vs. $a^Z_C/\Lambda^2$. The limits are described by analytical ellipses of equation $(x-x0)^2/a^2 + (y-y0)^2/b^2 = 1$ and rotated counter-clockwise by $\theta$ degrees, where $x$ and $y$ in the equation correspond to the $a_0^Z$ and $a_C^Z$ couplings, respectively.

Mean $p_T$ values of the $\Upsilon(1$S$)$, $\Upsilon(2$S$)$, and $\Upsilon(3$S$)$ states with $p_T\,>\,0\,GeV$ and $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $0<p_T<5\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $5<p_T<7\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $7<p_T<9\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $9<p_T<11\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $11<p_T<15\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $15<p_T<20\,$GeV range

Ratio $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $20<p_T<50\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $0<p_T<5\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $5<p_T<7\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $7<p_T<9\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $9<p_T<11\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $11<p_T<15\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $15<p_T<20\,$GeV range

Ratio $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ with $|y|\,<\,1.2$ as a function of track multiplicity $N_{track}$ in $20<p_T<50\,$GeV range

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of "forward" track multiplicity $N_{track}^{\Delta\phi}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$. Forward tracks are those with momentum direction in $\Delta\phi\,<\,\pi/3$ w.r.t. the $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of "transverse" track multiplicity $N_{track}^{\Delta\phi}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$. Transverse tracks are those with momentum direction in $\pi/3\,<\,\Delta\phi\,<\,2\pi/3$ w.r.t. the $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of "backward" track multiplicity $N_{track}^{\Delta\phi}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$. Backward tracks are those with momentum direction in $\Delta\phi\,>\,2\pi/3$ w.r.t. the $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$ and $N^{\Delta R}_{track}\,=\,0$ charged particles produced in $\Delta R\,<\,0.5$ cone around $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$ and $N^{\Delta R}_{track}\,=\,1$ charged particles produced in $\Delta R\,<\,0.5$ cone around $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$ and $N^{\Delta R}_{track}\,=\,2$ charged particles produced in $\Delta R\,<\,0.5$ cone around $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$ and $N^{\Delta R}_{track}\,>\,2$ charged particles produced in $\Delta R\,<\,0.5$ cone around $\Upsilon(n$S$)$ momentum direction.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ in transverse sphericity bin $0.00\,<\,S_T\,<\,0.55$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ in transverse sphericity bin $0.55\,<\,S_T\,<\,0.70$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ in transverse sphericity bin $0.70\,<\,S_T\,<\,0.85$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$.

Ratios $\Upsilon(2$S$)\,/\,\Upsilon(1$S$)$ and $\Upsilon(3$S$)\,/\,\Upsilon(1$S$)$ as functions of track multiplicity $N_{track}$ in transverse sphericity bin $0.85\,<\,S_T\,<\,1.00$ for $\Upsilon(n$S$)$ states with $p_T\,>\,7\,GeV$ and $|y|\,<\,1.2$.

Efficiency-corrected multiplicity bins used in the $\Upsilon(n$S$)$ ratio analysis and the corresponding mean number of charged particle tracks.

Charged particle yield distributions $Y(\Delta r)$ of b jets with $1 < p_{\text{T}}^{\text{trk}} < 12$ GeV are presented as functions of $\Delta r$.b jets with $p_T > 120$ GeV and charged particles with $1 < p^{\text{trk}}_{\text{T}} < 12$ GeV are used to construct the distributions as functions of $\Delta r$ differential $p_{\text{T}}^{\text{trk}}$ bins.

Charged particle yield differences of b to inclusive jets for $1 < p_{\text{T}}^{\text{trk}} < 12$ GeV are presented as functions of $\Delta r$.b jets with $p_T > 120$ GeV and charged particles with $1 < p^{\text{trk}}_{\text{T}} < 12$ GeV are used to construct the distributions as functions of $\Delta r$ differential $p_{\text{T}}^{\text{trk}}$ bins.

The jet shapes $\rho(\Delta r)$ of inclusive jets with $p_{\text{T}} > 120$ GeV and $p^{trk}_T > 1$ GeV are presented as functions of $\Delta r$ for differential $p_{\text{T}}^{\text{trk}}$ bin.

The jet shape distributions $\rho(\Delta r)$ of b jets with $p_{\text{T}} > 120$ GeV and $1< p^{\text{trk}}_{\text{T}} < 12$ GeV are presented as functions of $\Delta r$ for differential $p_{\text{T}}^{\text{trk}}$ bin.

The jet shapes $\rho (\Delta r)$ of inclusive jets with $p_{\text{T}} > 120$ GeV and $p^{\text{trk}}_{\text{T}} > 1$ GeV are presented as functions of $\Delta r.

The jet shapes $\rho (\Delta r)$ of b jets with $p_{\text{T}} > 120$ GeV and $p^{\text{trk}}_{\text{T}} > 1$ GeV are presented as functions of $\Delta r.

The b-to-inclusive jet shape ratios for $p_{\text{T}} > 120$ GeV and $p^{\text{trk}}_{\text{T}} > 1$ GeV are presented as functions of $\Delta r.

Invariant mass distribution of the diphoton pairs for the elastic selection region with events satisfying a < 0.005 (signal contribution is magnified by a factor 5000).

Predicted number of events having an elastic diphoton pair in association with a pair of protons observed within the range where the proton detectors have a radiation inefficiency less than 10%. The yields where the two-photon and two-proton systems mass and rapidity are matching at 2 and 3$\sigma$ are also quoted. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

Predicted number of events having an elastic diphoton pair in association with a pair of protons observed within the range where the proton detectors have a radiation inefficiency less than 10%. The yields where the two-photon and two-proton systems mass and rapidity are matching at 2 and 3$\sigma$ are also quoted. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

95% CL expected and observed frequentist upper limit for the LbL cross section within the $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$ search region.

95% CL expected and observed frequentist upper limit for the LbL cross section within the $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$ search region.

95% expected and observed one-dimensional limits on $\zeta_1$ and $\zeta_2$ anomalous LbyL production parameters, when the other parameter is set to zero. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

95% expected and observed one-dimensional limits on $\zeta_1$ and $\zeta_2$ anomalous LbyL production parameters, when the other parameter is set to zero. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

95% expected and observed limits on $\zeta_1$ and $\zeta_2$ anomalous LbyL production parameters. The parametric elliptic form is assumed: $a_0\zeta_1^2+a_1\zeta_1\zeta_2+a_2\zeta_2^2=1$. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

95% expected and observed limits on $\zeta_1$ and $\zeta_2$ anomalous LbyL production parameters. The parametric elliptic form is assumed: $a_0\zeta_1^2+a_1\zeta_1\zeta_2+a_2\zeta_2^2=1$. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

Measured normalized differential tt̄Z production cross section in the full phase space as a function of the transverse momentum of the Z boson, compared to the predictions obtained with the MadGraph5_aMC@NLO MC simulation, and to the theory prediction at NLO+NNLL accuracy (1905.07815). The distribution $1/σ\,Δσ$ is integrated over the bin, and $1/σ\,\mathrm{d}σ/\mathrm{d}p_{\mathrm{T}}(\mathrm{Z})$ is additionally divided by the bin width. The last bin includes the overflow contribution, but a finite bin width is used for the normalization.

Measured absolute differential tt̄Z production cross section in the full phase space as a function of $\cosθ^{*}_{\mathrm{Z}}$, compared to the predictions obtained with the MadGraph5_aMC@NLO MC simulation.

Measured absolute differential tt̄Z production cross section in the full phase space as a function of $\cosθ^{*}_{\mathrm{Z}}$, compared to the predictions obtained with the MadGraph5_aMC@NLO MC simulation.

Measured normalized differential tt̄Z production cross section in the full phase space as a function of $\cosθ^{*}_{\mathrm{Z}}$, compared to the predictions obtained with the MadGraph5_aMC@NLO MC simulation.

Measured normalized differential tt̄Z production cross section in the full phase space as a function of $\cosθ^{*}_{\mathrm{Z}}$, compared to the predictions obtained with the MadGraph5_aMC@NLO MC simulation.

Covariance matrix of the measured absolute differential tt̄Z production cross section in the full phase space as a function of the transverse momentum of the Z boson.

Covariance matrix of the measured absolute differential tt̄Z production cross section in the full phase space as a function of the transverse momentum of the Z boson.

Correlation matrix of the measured absolute differential tt̄Z production cross section in the full phase space as a function of the transverse momentum of the Z boson.

Correlation matrix of the measured absolute differential tt̄Z production cross section in the full phase space as a function of the transverse momentum of the Z boson.

Covariance matrix of the measured normalized differential tt̄Z production cross section in the full phase space as a function of the transverse momentum of the Z boson.

Covariance matrix of the measured normalized differential tt̄Z production cross section in the full phase space as a function of the transverse momentum of the Z boson.

Correlation matrix of the measured normalized differential tt̄Z production cross section in the full phase space as a function of the transverse momentum of the Z boson.

Correlation matrix of the measured normalized differential tt̄Z production cross section in the full phase space as a function of the transverse momentum of the Z boson.

Covariance matrix of the measured absolute differential tt̄Z production cross section in the full phase space as a function of $\cosθ^{*}_{\mathrm{Z}}$.

Covariance matrix of the measured absolute differential tt̄Z production cross section in the full phase space as a function of $\cosθ^{*}_{\mathrm{Z}}$.

Correlation matrix of the measured absolute differential tt̄Z production cross section in the full phase space as a function of $\cosθ^{*}_{\mathrm{Z}}$.

Correlation matrix of the measured absolute differential tt̄Z production cross section in the full phase space as a function of $\cosθ^{*}_{\mathrm{Z}}$.

Covariance matrix of the measured normalized differential tt̄Z production cross section in the full phase space as a function of $\cosθ^{*}_{\mathrm{Z}}$.

Covariance matrix of the measured normalized differential tt̄Z production cross section in the full phase space as a function of $\cosθ^{*}_{\mathrm{Z}}$.

Correlation matrix of the measured normalized differential tt̄Z production cross section in the full phase space as a function of $\cosθ^{*}_{\mathrm{Z}}$.

Correlation matrix of the measured normalized differential tt̄Z production cross section in the full phase space as a function of $\cosθ^{*}_{\mathrm{Z}}$.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the 1$^\text{st}$ jet $p_{\text{T}}$, $p_{\text{T}}(\text{j}_1)$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the 2$^\text{nd}$ jet $p_{\text{T}}$, $p_{\text{T}}(\text{j}_2)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the 2$^\text{nd}$ jet $p_{\text{T}}$, $p_{\text{T}}(\text{j}_2)$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the 3$^\text{rd}$ jet $p_{\text{T}}$, $p_{\text{T}}(\text{j}_3)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the 3$^\text{rd}$ jet $p_{\text{T}}$, $p_{\text{T}}(\text{j}_3)$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the 4$^\text{th}$ jet $p_{\text{T}}$, $p_{\text{T}}(\text{j}_4)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the 4$^\text{th}$ jet $p_{\text{T}}$, $p_{\text{T}}(\text{j}_4)$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the 5$^\text{th}$ jet $p_{\text{T}}$, $p_{\text{T}}(\text{j}_5)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the 5$^\text{th}$ jet $p_{\text{T}}$, $p_{\text{T}}(\text{j}_5)$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the 1$^\text{st}$ jet $|y|$, $|y(\text{j}_1)|$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the 1$^\text{st}$ jet $|y|$, $|y(\text{j}_1)|$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the 2$^\text{nd}$ jet $|y|$, $|y(\text{j}_2)|$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the 2$^\text{nd}$ jet $|y|$, $|y(\text{j}_2)|$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the 3$^\text{rd}$ jet $|y|$, $|y(\text{j}_3)|$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the 3$^\text{rd}$ jet $|y|$, $|y(\text{j}_3)|$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the 4$^\text{th}$ jet $|y|$, $|y(\text{j}_4)|$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the 4$^\text{th}$ jet $|y|$, $|y(\text{j}_4)|$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the 5$^\text{th}$ jet $|y|$, $|y(\text{j}_5)|$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the 5$^\text{th}$ jet $|y|$, $|y(\text{j}_5)|$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the Z boson $|y|$, $|y(\text{Z})|$, for events with at least one jet and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the Z boson $|y|$, $|y(\text{Z})|$, for events with at least one jet.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the Z boson $|y|$, $|y(\text{Z})|$ for events with at least one jet and $p_\text{T}(\text{Z}) > 150\,$GeV, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the Z boson $|y|$, $|y(\text{Z})|$ for events with at least one jet and $p_\text{T}(\text{Z}) > 150\,$GeV.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the Z boson $|y|$, $|y(\text{Z})|$ for events with at least $p_\text{T}(\text{Z}) > 300\,$GeV, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the Z boson $|y|$, $|y(\text{Z})|$ for events with at least $p_\text{T}(\text{Z}) > 300\,$GeV.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{diff}}$ of the Z boson and the leading jet, $y_{\text{diff}}(\text{Z}, \text{j}_1)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{diff}}$ of the Z boson and the leading jet, $y_{\text{diff}}(\text{Z}, \text{j}_1)$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{diff}}$ of the Z boson and the leading jet, $y_{\text{diff}}(\text{Z}, \text{j}_1)$ for events with $p_{\text{T}}(\text{Z}) > 150\,$GeV, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{diff}}$ of the Z boson and the leading jet, $y_{\text{diff}}(\text{Z}, \text{j}_1)$ for events with $p_{\text{T}}(\text{Z}) > 150\,$GeV.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{diff}}$ of the Z boson and the leading jet, $y_{\text{diff}}(\text{Z}, \text{j}_1)$ for events with $p_{\text{T}}(\text{Z}) > 300\,$GeV, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{diff}}$ of the Z boson and the leading jet, $y_{\text{diff}}(\text{Z}, \text{j}_1)$ for events with $p_{\text{T}}(\text{Z}) > 300\,$GeV.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{diff}}$ of the Z boson and the 2$^\text{nd}$leading jet, $y_{\text{diff}}(\text{Z}, \text{j}_2)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{diff}}$ of the Z boson and the 2$^\text{nd}$leading jet, $y_{\text{diff}}(\text{Z}, \text{j}_2)$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{diff}}$ of the Z boson and the 3$^\text{rd}$leading jet, $y_{\text{diff}}(\text{Z}, \text{j}_3)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{diff}}$ of the Z boson and the 3$^\text{rd}$leading jet, $y_{\text{diff}}(\text{Z}, \text{j}_3)$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{diff}}$ of the Z boson and the system constituted of the two leading jets, $y_{\text{diff}}(\text{Z}, \text{dijet})$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{diff}}$ of the Z boson and the system constituted of the two leading jets, $y_{\text{diff}}(\text{Z}, \text{dijet})$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{diff}}$ of the two leading jets, $y_{\text{diff}}(\text{j}_1, \text{j}_2)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{diff}}$ of the two leading jets, $y_{\text{diff}}(\text{j}_1, \text{j}_2)$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{sum}}$ of the two leading jets, $y_{\text{sum}}(\text{j}_1, \text{j}_2)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{sum}}$ of the two leading jets, $y_{\text{sum}}(\text{j}_1, \text{j}_2)$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{sum}}$ of the Z boson and the leading jet, $y_{\text{sum}}(\text{Z}, \text{j}_1)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{sum}}$ of the Z boson and the leading jet, $y_{\text{sum}}(\text{Z}, \text{j}_1)$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{sum}}$ of the Z boson and the leading jet, $y_{\text{sum}}(\text{Z}, \text{j}_1)$ for events with $p_{\text{T}}(\text{Z}) > 150\,$GeV, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{sum}}$ of the Z boson and the leading jet, $y_{\text{sum}}(\text{Z}, \text{j}_1)$ for events with $p_{\text{T}}(\text{Z}) > 150\,$GeV.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{sum}}$ of the Z boson and the leading jet, $y_{\text{sum}}(\text{Z}, \text{j}_1)$ for events with $p_{\text{T}}(\text{Z}) > 300\,$GeV, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{sum}}$ of the Z boson and the leading jet, $y_{\text{sum}}(\text{Z}, \text{j}_1)$ for events with $p_{\text{T}}(\text{Z}) > 300\,$GeV.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{sum}}$ of the Z boson and the 2$^\text{nd}$leading jet, $y_{\text{sum}}(\text{Z}, \text{j}_2)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{sum}}$ of the Z boson and the 2$^\text{nd}$leading jet, $y_{\text{sum}}(\text{Z}, \text{j}_2)$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{sum}}$ of the Z boson and the 3$^\text{rd}$leading jet, $y_{\text{sum}}(\text{Z}, \text{j}_3)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{sum}}$ of the Z boson and the 3$^\text{rd}$leading jet, $y_{\text{sum}}(\text{Z}, \text{j}_3)$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{sum}}$ of the Z boson and the system constituted of the two leading jets, $y_{\text{sum}}(\text{Z}, \text{dijet})$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the $y_{\text{sum}}$ of the Z boson and the system constituted of the two leading jets, $y_{\text{sum}}(\text{Z}, \text{dijet})$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of $H_\text{T}$ for events with at least one jet, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of $H_\text{T}$ for events with at least one jet.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of $H_\text{T}$ for events with at least two jets, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of $H_\text{T}$ for events with at least two jets.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of $H_\text{T}$ for events with at least three jets, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of $H_\text{T}$ for events with at least three jets.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of $H_\text{T}$ for events with at least four jets, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of $H_\text{T}$ for events with at least four jets.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of $H_\text{T}$ for events with at least five jets, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of $H_\text{T}$ for events with at least five jets.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the leading jet for events with at least ine jet, $\Delta\Phi(\text{Z},\text{j}_1)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the leading jet for events with at least ine jet, $\Delta\Phi(\text{Z},\text{j}_1)$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the leading jet for events with at least two jets, $\Delta\Phi(\text{Z},\text{j}_1)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the leading jet for events with at least two jets, $\Delta\Phi(\text{Z},\text{j}_1)$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the leading jet for events with at least three jets, $\Delta\Phi(\text{Z},\text{j}_1)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the leading jet for events with at least three jets, $\Delta\Phi(\text{Z},\text{j}_1)$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the 2$^\text{nd}$ leading jet for events with at least three jets, $\Delta\Phi(\text{Z},\text{j}_1)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the 2$^\text{nd}$ leading jet for events with at least three jets, $\Delta\Phi(\text{Z},\text{j}_1)$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the 3$^\text{rd}$ leading jet for events with at least three jets, $\Delta\Phi(\text{Z},\text{j}_1)$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the 3$^\text{rd}$ leading jet for events with at least three jets, $\Delta\Phi(\text{Z},\text{j}_1)$.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the leading jet$\Delta\Phi(\text{Z},\text{j}_1)$, for events with $p_{\text{T}}(\text{Z})>150\,$GeV and at least two jets, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the leading jet$\Delta\Phi(\text{Z},\text{j}_1)$, for events with $p_{\text{T}}(\text{Z})>150\,$GeV and at least two jets.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the leading jet$\Delta\Phi(\text{Z},\text{j}_1)$, for events with $p_{\text{T}}(\text{Z})>150\,$GeV and at least two jets, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the leading jet$\Delta\Phi(\text{Z},\text{j}_1)$, for events with $p_{\text{T}}(\text{Z})>150\,$GeV and at least two jets.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the leading jet$\Delta\Phi(\text{Z},\text{j}_1)$, for events with $p_{\text{T}}(\text{Z})>150\,$GeV and at least three jets, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the leading jet$\Delta\Phi(\text{Z},\text{j}_1)$, for events with $p_{\text{T}}(\text{Z})>150\,$GeV and at least three jets.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the 2$^{\text{nd}}$ leading jet$\Delta\Phi(\text{Z},\text{j}_2)$, for events with $p_{\text{T}}(\text{Z})>150\,$GeV and at least three jets, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the 2$^{\text{nd}}$ leading jet$\Delta\Phi(\text{Z},\text{j}_2)$, for events with $p_{\text{T}}(\text{Z})>150\,$GeV and at least three jets.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the 3$^{\text{rd}}$ leading jet$\Delta\Phi(\text{Z},\text{j}_3)$, for events with $p_{\text{T}}(\text{Z})>150\,$GeV and at least three jets, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the 3$^{\text{rd}}$ leading jet$\Delta\Phi(\text{Z},\text{j}_3)$, for events with $p_{\text{T}}(\text{Z})>150\,$GeV and at least three jets.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the leading jet, $\Delta\Phi(\text{Z},\text{j}_1)$, for events with $p_{\text{T}}(\text{Z})>300\,$GeV and at least one jet, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the leading jet, $\Delta\Phi(\text{Z},\text{j}_1)$, for events with $p_{\text{T}}(\text{Z})>300\,$GeV and at least one jet.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the leading jet, $\Delta\Phi(\text{Z},\text{j}_1)$, for events with $p_{\text{T}}(\text{Z})>300\,$GeV and at least two jets, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the leading jet, $\Delta\Phi(\text{Z},\text{j}_1)$, for events with $p_{\text{T}}(\text{Z})>300\,$GeV and at least two jets.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the leading jet, $\Delta\Phi(\text{Z},\text{j}_1)$, for events with $p_{\text{T}}(\text{Z})>300\,$GeV and at least three jets, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the leading jet, $\Delta\Phi(\text{Z},\text{j}_1)$, for events with $p_{\text{T}}(\text{Z})>300\,$GeV and at least three jets.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the 2$^{\text{nd}}$ leading jet, $\Delta\Phi(\text{Z},\text{j}_2)$, for events with $p_{\text{T}}(\text{Z})>300\,$GeV, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the 2$^{\text{nd}}$ leading jet, $\Delta\Phi(\text{Z},\text{j}_2)$, for events with $p_{\text{T}}(\text{Z})>300\,$GeV.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the 3$^{\text{rd}}$ leading jet, $\Delta\Phi(\text{Z},\text{j}_3)$, for events with $p_{\text{T}}(\text{Z})>300\,$GeV, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the 3$^{\text{rd}}$ leading jet, $\Delta\Phi(\text{Z},\text{j}_3)$, for events with $p_{\text{T}}(\text{Z})>300\,$GeV.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the leading jet, $\Delta\Phi(\text{Z},\text{j}_1)$, for events with at least three jets, $p_{\text{T}}(\text{Z})>150\,$GeV, and $H_{\text{T}}>300\,$GeV, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the leading jet, $\Delta\Phi(\text{Z},\text{j}_1)$, for events with at least three jets, $p_{\text{T}}(\text{Z})>150\,$GeV, and $H_{\text{T}}>300\,$GeV.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the 2$^{\text{nd}}$ leading jet, $\Delta\Phi(\text{Z},\text{j}_2)$, for events with at least three jets, $p_{\text{T}}(\text{Z})>150\,$GeV, and $H_{\text{T}}>300\,$GeV, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the 2$^{\text{nd}}$ leading jet, $\Delta\Phi(\text{Z},\text{j}_2)$, for events with at least three jets, $p_{\text{T}}(\text{Z})>150\,$GeV, and $H_{\text{T}}>300\,$GeV.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the 3$^{\text{rd}}$ leading jet, $\Delta\Phi(\text{Z},\text{j}_1)$, for events with at least three jets, $p_{\text{T}}(\text{Z})>150\,$GeV, and $H_{\text{T}}>300\,$GeV, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the Z boson and the 3$^{\text{rd}}$ leading jet, $\Delta\Phi(\text{Z},\text{j}_1)$, for events with at least three jets, $p_{\text{T}}(\text{Z})>150\,$GeV, and $H_{\text{T}}>300\,$GeV.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the two leading jets, $\Delta\Phi(\text{j}_1,\text{j}_2)$, for events with at least three jets, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the two leading jets, $\Delta\Phi(\text{j}_1,\text{j}_2)$, for events with at least three jets.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the two 1$^{\text{st}}$ and 3$^{\text{rd}}$ leading jets, $\Delta\Phi(\text{j}_1,\text{j}_3)$, for events with at least three jets, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the two 1$^{\text{st}}$ and 3$^{\text{rd}}$ leading jets, $\Delta\Phi(\text{j}_1,\text{j}_3)$, for events with at least three jets.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the two 2$^{\text{nd}}$ and 3$^{\text{rd}}$ leading jets, $\Delta\Phi(\text{j}_1,\text{j}_3)$, for events with at least three jets, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the two 2$^{\text{nd}}$ and 3$^{\text{rd}}$ leading jets, $\Delta\Phi(\text{j}_1,\text{j}_3)$, for events with at least three jets.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the two leading jets, $\Delta\Phi(\text{j}_1,\text{j}_2)$, for events with at least three jets and $p_{\text{T}}(\text{Z})>150\,$GeV, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the two leading jets, $\Delta\Phi(\text{j}_1,\text{j}_2)$, for events with at least three jets and $p_{\text{T}}(\text{Z})>150\,$GeV.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the 1$^{\text{st}}$ and 3$^{\text{rd}}$ leading jets, $\Delta\Phi(\text{j}_1,\text{j}_3)$, for events with at least three jets and $p_{\text{T}}(\text{Z})>150\,$GeV, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the 1$^{\text{st}}$ and 3$^{\text{rd}}$ leading jets, $\Delta\Phi(\text{j}_1,\text{j}_3)$, for events with at least three jets and $p_{\text{T}}(\text{Z})>150\,$GeV.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the 2$^{\text{nd}}$ and 3$^{\text{rd}}$ leading jets, $\Delta\Phi(\text{j}_2,\text{j}_3)$, for events with at least three jets and $p_{\text{T}}(\text{Z})>150\,$GeV, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the 2$^{\text{nd}}$ and 3$^{\text{rd}}$ leading jets, $\Delta\Phi(\text{j}_2,\text{j}_3)$, for events with at least three jets and $p_{\text{T}}(\text{Z})>150\,$GeV.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the 1$^{\text{st}}$ and 2$^{\text{nd}}$ leading jets, $\Delta\Phi(\text{j}_1,\text{j}_2)$, for events with at least three jets and $p_{\text{T}}(\text{Z})>300\,$GeV, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the 1$^{\text{st}}$ and 2$^{\text{nd}}$ leading jets, $\Delta\Phi(\text{j}_1,\text{j}_2)$, for events with at least three jets and $p_{\text{T}}(\text{Z})>300\,$GeV.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the 1$^{\text{st}}$ and 3$^{\text{rd}}$ leading jets, $\Delta\Phi(\text{j}_1,\text{j}_3)$, for events with at least three jets and $p_{\text{T}}(\text{Z})>300\,$GeV, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the 1$^{\text{st}}$ and 3$^{\text{rd}}$ leading jets, $\Delta\Phi(\text{j}_1,\text{j}_3)$, for events with at least three jets and $p_{\text{T}}(\text{Z})>300\,$GeV.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the 2$^{\text{nd}}$ and 3$^{\text{rd}}$ leading jets, $\Delta\Phi(\text{j}_1,\text{j}_3)$, for events with at least three jets and $p_{\text{T}}(\text{Z})>300\,$GeV, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the azimuthal angle between the 2$^{\text{nd}}$ and 3$^{\text{rd}}$ leading jets, $\Delta\Phi(\text{j}_1,\text{j}_3)$, for events with at least three jets and $p_{\text{T}}(\text{Z})>300\,$GeV.

The cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the mass of the system made of the two leading jets, $m_{\text{j}_1\text{j}_1}$, and breakdown of the relative uncertainty.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production measured as a function of the mass of the system made of the two leading jets, $m_{\text{j}_1\text{j}_1}$.

The cross section for Z($\rightarrow\ell\ell$) + jets production as a function of the leading jet transverse momentum and rapidity

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production as a function of the leading jet transverse momentum and rapidity

The cross section for Z($\rightarrow\ell\ell$) + jets production as a function of the Z boson and leading jet rapidities

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production as a function of the Z boson and leading jet rapidities

The cross section for Z($\rightarrow\ell\ell$) + jets production as a function of the rapidities of the Z boson and leading jet, and of the transverse momentum of the jet for the same configuration.

Bin-to-bin correlation in the the cross section for Z($\rightarrow\ell\ell$) + jets production as a function of the rapidities of the Z boson and leading jet, and of the transverse momentum of the jet for the same configuration.

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.

Differential cross section times dimuon branching fraction of Y(1S) as a function of $y^{Y}_{CM}$ in pPb collisions. The global uncertainty arises from the integrated luminosity uncertainty in pPb collisions.

Differential cross section times dimuon branching fraction of Y(2S) as a function of $y^{Y}_{CM}$ in pPb collisions. The global uncertainty arises from the integrated luminosity uncertainty in pPb collisions.

Differential cross section times dimuon branching fraction of Y(3S) as a function of $y^{Y}_{CM}$ in pPb collisions. The global uncertainty arises from the integrated luminosity uncertainty in pPb collisions.

Differential cross section times dimuon branching fraction of Y(1S) as a function of pT in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(2S) as a function of pT in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(3S) as a function of pT in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(1S) as a function of $|y^{Y}_{CM}|$ in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(2S) as a function of $|y^{Y}_{CM}|$ in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(3S) as a function of $|y^{Y}_{CM}|$ in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Nuclear modification factor of Y(1S) as a function of pT. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) as a function of pT. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) as a function of pT. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(1S) as a function of $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) as a function of $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) as a function of $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(1S) at forward and backward $y^{Y}_{CM}$ for pT < 6 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) at forward and backward $y^{Y}_{CM}$ for pT < 6 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) at forward and backward $y^{Y}_{CM}$ for pT < 6 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(1S) at forward and backward $y^{Y}_{CM}$ for 6 < pT < 30 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) at forward and backward $y^{Y}_{CM}$ for 6 < pT < 30 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) at forward and backward $y^{Y}_{CM}$ for 6 < pT < 30 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

RFB of Y(1S) versus $N^{|\eta_{lab}|<2.4}_{tracks}$.

RFB of Y(2S) versus $N^{|\eta_{lab}|<2.4}_{tracks}$.

RFB of Y(3S) versus $N^{|\eta_{lab}|<2.4}_{tracks}$.

RFB of Y(1S) versus $E^{|\eta_{lab}|>4}_{T}$.

RFB of Y(2S) versus $E^{|\eta_{lab}|>4}_{T}$.

RFB of Y(3S) versus $E^{|\eta_{lab}|>4}_{T}$.

Nuclear modification factor of Y(1S), Y(2S), and Y(3S) integrated over pT and $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Exclusion limit for BrHZX_BrXee

Exclusion limit for BrHZX_BrXmumu

Exclusion limit for BrHZX_BrXll

Exclusion limit for cah

Exclusion limit for cZh

Exclusion limit for kappa

Expected and observed yields from the standard model processeses in the WW signal region. The combination of the statistical and systematic uncertainties are shown. The predicted yields are shown with their best-fit normalizations from the simultaneous fit.

Expected and observed yields from the standard model processeses in the WW signal region. The combination of the statistical and systematic uncertainties are shown. The predicted yields are shown with their best-fit normalizations from the simultaneous fit.

Expected and observed yields from the standard model processeses in the WW signal region. The combination of the statistical and systematic uncertainties are shown. The predicted yields are shown with their best-fit normalizations from the simultaneous fit.

Expected and observed yields from the standard model processeses in the WZ signal region. The combination of the statistical and systematic uncertainties are shown. The predicted yields are shown with their best-fit normalizations from the simultaneous fit.

Expected and observed yields from the standard model processeses in the WZ signal region. The combination of the statistical and systematic uncertainties are shown. The predicted yields are shown with their best-fit normalizations from the simultaneous fit.

Expected and observed yields from the standard model processeses in the WZ signal region. The combination of the statistical and systematic uncertainties are shown. The predicted yields are shown with their best-fit normalizations from the simultaneous fit.

The measured inclusive fiducial cross section measurements. The WW fiducial region is defined by requiring two same-sign leptons with $p_{T}>20$, $|\eta|<2.5$, and $m_{ll}>20$, and two jets with $m_{jj}>500$ and $|\Delta \eta_{jj}|>2.5$. The jets at generator level are clustered from stable particles, excluding neutrinos, using the anti-kt clustering algorithm with R = 0.4, and are required to have $p_{T}>50$ and $|\eta|<4.7$. The jets within $\Delta R<0.4$ of the selected charged leptons are not included. The WZ fiducial region is defined by requiring three leptons with $p_{T}>20$, $|\eta|<2.5$, a pair of opposite charge same-flavor lepton pair with $|m_{ll}-m_{Z}|<15$, and two jets with $m_{jj}>500$ and $|\Delta \eta_{jj}|>2.5$.

The measured inclusive fiducial cross section measurements. The WW fiducial region is defined by requiring two same-sign leptons with $p_{T}>20$, $|\eta|<2.5$, and $m_{ll}>20$, and two jets with $m_{jj}>500$ and $|\Delta \eta_{jj}|>2.5$. The jets at generator level are clustered from stable particles, excluding neutrinos, using the anti-kt clustering algorithm with R = 0.4, and are required to have $p_{T}>50$ and $|\eta|<4.7$. The jets within $\Delta R<0.4$ of the selected charged leptons are not included. The WZ fiducial region is defined by requiring three leptons with $p_{T}>20$, $|\eta|<2.5$, a pair of opposite charge same-flavor lepton pair with $|m_{ll}-m_{Z}|<15$, and two jets with $m_{jj}>500$ and $|\Delta \eta_{jj}|>2.5$.

The measured inclusive fiducial cross section measurements. The WW fiducial region is defined by requiring two same-sign leptons with $p_{T}>20$, $|\eta|<2.5$, and $m_{ll}>20$, and two jets with $m_{jj}>500$ and $|\Delta \eta_{jj}|>2.5$. The jets at generator level are clustered from stable particles, excluding neutrinos, using the anti-kt clustering algorithm with R = 0.4, and are required to have $p_{T}>50$ and $|\eta|<4.7$. The jets within $\Delta R<0.4$ of the selected charged leptons are not included. The WZ fiducial region is defined by requiring three leptons with $p_{T}>20$, $|\eta|<2.5$, a pair of opposite charge same-flavor lepton pair with $|m_{ll}-m_{Z}|<15$, and two jets with $m_{jj}>500$ and $|\Delta \eta_{jj}|>2.5$.

Distributions of $m_{jj}$ in the WW signal region.

Distributions of $m_{jj}$ in the WW signal region.

Distributions of $m_{jj}$ in the WW signal region.

Distributions of $m_{ll}$ in the WW signal region.

Distributions of $m_{ll}$ in the WW signal region.

Distributions of $m_{ll}$ in the WW signal region.

Distributions of $m_{jj}$ in the WZ signal region.

Distributions of $m_{jj}$ in the WZ signal region.

Distributions of $m_{jj}$ in the WZ signal region.

Distributions of BDT score in the WZ signal region.

Distributions of BDT score in the WZ signal region.

Distributions of BDT score in the WZ signal region.

Absolute WW cross section in $m_{jj}$ bins.

Absolute WW cross section in $m_{jj}$ bins.

Absolute WW cross section in $m_{jj}$ bins.

Normalized WW cross section in $m_{jj}$ bins.

Normalized WW cross section in $m_{jj}$ bins.

Normalized WW cross section in $m_{jj}$ bins.

Absolute WW cross section in $m_{ll}$ bins.

Absolute WW cross section in $m_{ll}$ bins.

Absolute WW cross section in $m_{ll}$ bins.

Normalized WW cross section in $m_{ll}$ bins.

Normalized WW cross section in $m_{ll}$ bins.

Normalized WW cross section in $m_{ll}$ bins.

Absolute WW cross section in $p_{T}^{l max}$ bins.

Absolute WW cross section in $p_{T}^{l max}$ bins.

Absolute WW cross section in $p_{T}^{l max}$ bins.

Normalized WW cross section in $p_{T}^{l max}$ bins.

Normalized WW cross section in $p_{T}^{l max}$ bins.

Normalized WW cross section in $p_{T}^{l max}$ bins.

Absolute WZ cross section in $m_{jj}$ bins.

Absolute WZ cross section in $m_{jj}$ bins.

Absolute WZ cross section in $m_{jj}$ bins.

Normalized WZ cross section in $m_{jj}$ bins.

Normalized WZ cross section in $m_{jj}$ bins.

Normalized WZ cross section in $m_{jj}$ bins.

Distributions of $m_{T}^{WW}$ in the WW signal region.

Distributions of $m_{T}^{WW}$ in the WW signal region.

Distributions of $m_{T}^{WW}$ in the WW signal region.

Distributions of $m_{T}^{WZ}$ in the WZ signal region.

Distributions of $m_{T}^{WZ}$ in the WZ signal region.

Distributions of $m_{T}^{WZ}$ in the WZ signal region.

Observed and expected lower and upper 95\% confidence level limits in TeV$^{-4}$ on the parameters of the quartic, obtained without using any unitarization procedure.

Observed and expected lower and upper 95\% confidence level limits in TeV$^{-4}$ on the parameters of the quartic, obtained without using any unitarization procedure.

Observed and expected lower and upper 95\% confidence level limits in TeV$^{-4}$ on the parameters of the quartic, obtained without using any unitarization procedure.

Observed and expected lower and upper 95\% confidence level limits in TeV$^{-4}$ on the parameters of the quartic by cutting the EFT expansion at the unitarity limit.

Observed and expected lower and upper 95\% confidence level limits in TeV$^{-4}$ on the parameters of the quartic by cutting the EFT expansion at the unitarity limit.

Observed and expected lower and upper 95\% confidence level limits in TeV$^{-4}$ on the parameters of the quartic by cutting the EFT expansion at the unitarity limit.

Measured diferential cross sections as a function of the transverse momentum of the lepton from the W decay.

Measured diferential cross sections ratio $R=\sigma(W^+ + \overline{c}) / \sigma(W^- + c)$ as a function of the absolute value of the pseudorapidity of the lepton from the W decay.

Measured diferential cross sections ratio $R=\sigma(W^+ + \overline{c}) / \sigma(W^- + c)$ as a function of the transverse momentum of the lepton from the W decay.

The acceptance- and efficiency-corrected yields for the $\mathrm{B}^{+}$ mesons, scaled by $T_{\mathrm{AA}}$ and $N_{\text{MB}}$ in PbPb collisions at $\sqrt{s_{NN}}=5.02$ TeV. The results are shown as a function of the event $\langle N_{part} \rangle$.

The acceptance- and efficiency-corrected $\mathrm{B}^{0}_{\mathrm{s}}/\mathrm{B}^{+}$ yields in PbPb collisions at $\sqrt{s_{NN}}=5.02$. TeV. The results are shown as a function of the meson $p_T$.

The acceptance- and efficiency-corrected $\mathrm{B}^{0}_{\mathrm{s}}/\mathrm{B}^{+}$ yields in PbPb collisions at $\sqrt{s_{NN}}=5.02$. The results are shown as a function of the event $\langle N_{part} \rangle$.

The nuclear modification factor R$_{\rm AA}$, as a function of $z$ in the 0-20 and 20-30 collisions centrality ranges.

Distributions of $\Delta\phi_{\mathrm{trk,Z}}$ in 30-50% centrality PbPb collisions at 5.02 TeV.

Distributions of $\Delta\phi_{\mathrm{trk,Z}}$ in 0-30% centrality PbPb collisions at 5.02 TeV.

PbPb - pp difference of $\Delta\phi_{\mathrm{trk,Z}}$ distributions for 70-90% centrality PbPb collisions at 5.02 TeV.

PbPb - pp difference of $\Delta\phi_{\mathrm{trk,Z}}$ distributions for 50-70% centrality PbPb collisions at 5.02 TeV.

PbPb - pp difference of $\Delta\phi_{\mathrm{trk,Z}}$ distributions for 30-50% centrality PbPb collisions at 5.02 TeV.

PbPb - pp difference of $\Delta\phi_{\mathrm{trk,Z}}$ distributions for 0-30% centrality PbPb collisions at 5.02 TeV.

Distributions of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ in pp collisions at 5.02 TeV.

Distributions of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ in 70-90% centrality PbPb collisions at 5.02 TeV.

Distributions of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ in 50-70% centrality PbPb collisions at 5.02 TeV.

Distributions of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ in 30-50% centrality PbPb collisions at 5.02 TeV.

Distributions of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ in 0-30% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ distributions for 70-90% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ distributions for 50-70% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ distributions for 30-50% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ distributions for 0-30% centrality PbPb collisions at 5.02 TeV.

Distributions of p$^{\mathrm{trk}}_{\mathrm{T}}$ in pp collisions at 5.02 TeV.

Distributions of p$^{\mathrm{trk}}_{\mathrm{T}}$ in pp collisions at 5.02 TeV.

Distributions of p$^{\mathrm{trk}}_{\mathrm{T}}$ in 70-90% centrality PbPb collisions at 5.02 TeV.

Distributions of p$^{\mathrm{trk}}_{\mathrm{T}}$ in 50-70% centrality PbPb collisions at 5.02 TeV.

Distributions of p$^{\mathrm{trk}}_{\mathrm{T}}$ in 30-50% centrality PbPb collisions at 5.02 TeV.

Distributions of p$^{\mathrm{trk}}_{\mathrm{T}}$ in 0-30% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of p$^{\mathrm{trk}}_{\mathrm{T}}$ distributions for 70-90% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of p$^{\mathrm{trk}}_{\mathrm{T}}$ distributions for 50-70% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of p$^{\mathrm{trk}}_{\mathrm{T}}$ distributions for 30-50% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of p$^{\mathrm{trk}}_{\mathrm{T}}$ distributions for 0-30% centrality PbPb collisions at 5.02 TeV.

Differential cross section in bins of pT(Z). Values are expressed as a fraction of the total cross section. The inclusive final state is shown.

Differential cross section in bins of pT(Z). Values are expressed as a fraction of the total cross section. The emm and mmm final states are shown.

Differential cross section in bins of mass of the WZ system. Values are expressed as a fraction of the total cross section.

Measured WZ production cross sections computed separately in each of the flavour categories.

Measured fiducial cross sections and their corresponding uncertainties for each of the individual flavour categories, as well as for the combination of the four. The combined value is the result of a symultaneous fit to the four categories, therefore both the central value and its total uncertaintydiffer from the sum of the central values and the quadratic sum of the uncertainties respectively, because of correlations among sources of uncertainty in the categorized values.

Differential cross section in bins of pT(leading jet). Values are expressed as a fraction of the total cross section. The inclusive final state is shown.

Expected and observed yields for each of the relevant processes and flavour categories. Combined statistical and systematic uncertainties are shown for each case except for the observed data yields for which only statistical uncertainties are presented. All expected yields correspond to quantities estimated after the maximum likelihood fit. Uncertainties are computed taking into account the full correlation matrix between sources of uncertainty, processes, and flavour categories.

Differential cross section in bins of pT(leading jet). Values are expressed as a fraction of the total cross section. The eee, eem, emm, and mmm final states are shown.

Summary of the total postfit impact of each uncertainty source on the uncertainty in the signal strength measurement, for the four flavour categories and their combination. Theoretical uncertainties are only included in the signal acceptance during the extrapolation to the total phase space, so they are not included in the likelihood fit. The values are percentages and correspond to half the difference between the up and down variation of each systematic uncertainty component.

Efficiencies estimated as transfer factors from the fiducial region to the signal region using generator-level information, for an integrated luminosity L of 35.9 inverse femtobarn. Statistical, scale, and PDF uncertainties are later propagated to the final cross section measurement.

Measured WZ production cross section computed in the inclusive final state and split by W boson charge, and asymmetry ratio.

Average hadronic energy reconstructed in the region −6.6 < eta < −5.2 as a function of the number of reconstructed tracks for abs(eta)<2.

Ratio of the average electromagnetic and hadronic energies reconstructed in the region −6.6 < eta < −5.2 as a function of the number of reconstructed tracks for abs(eta)<2.

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

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

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

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

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

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

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

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

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

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

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

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

Invariant $p_{T}$-differential spectra of $\Lambda + \bar{\Lambda}$ in p+p and p+Pb at $\sqrt{s}$=5.02 TeV in various $y_{CM}$ ranges

Invariant $p_{T}$-differential spectra of $\Xi- + \bar{\Xi+}$ in p+p and p+Pb at $\sqrt{s}$=5.02 TeV in various |$y_{CM}$| ranges

Invariant $p_{T}$-differential spectra of $\Xi- + \bar{\Xi+}$ in p+p and p+Pb at $\sqrt{s}$=5.02 TeV in various $y_{CM}$ ranges

Invariant $p_{T}$-differential spectra of $\Omega- + \bar{\Omega+}$ in p+p and p+Pb at $\sqrt{s}$=5.02 TeV in |$y_{CM}$| < 1.8

Invariant $p_{T}$-differential spectra of $\Omega- + \bar{\Omega+}$ in p+p and p+Pb at $\sqrt{s}$=5.02 TeV in |$y_{CM}$| < 1.8

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in |$y_{CM}$| < 1.8

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in |$y_{CM}$| < 1.8

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in -1.8 < $y_{CM}$ < 0

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in |$y_{CM}$| < 1.8

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0 < $y_{CM}$ < 1.8

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in |$y_{CM}$| < 1.8

Invariant $p_{T}$-differential spectra of ${K_{0}}^{S}$ in p+Pb at $\sqrt{s}$=5.02 TeV in various $y_{CM}$ ranges

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in |$y_{CM}$| < 1.8

Invariant $p_{T}$-differential spectra of $\Lambda + \bar{\Lambda}$ in p+Pb at $\sqrt{s}$=5.02 TeV in various $y_{CM}$ ranges

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in -1.8 < $y_{CM}$ < 0

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0.3 < $y_{CM}$ < 0.8

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in -1.8 < $y_{CM}$ < 0

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0.8 < $y_{CM}$ < 1.3

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in -1.8 < $y_{CM}$ < 0

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 1.3 < $y_{CM}$ < 1.8

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0 < $y_{CM}$ < 1.8

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0 < $y_{CM}$ < 1.8

$R_{pPb}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0 < $y_{CM}$ < 1.8

Invariant $p_{T}$-differential spectra of ${K_{0}}^{S}$ in p+Pb at $\sqrt{s}$=5.02 TeV in various $y_{CM}$ ranges

Invariant $p_{T}$-differential spectra of $\Lambda + \bar{\Lambda}$ in p+Pb at $\sqrt{s}$=5.02 TeV in various $y_{CM}$ ranges

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0.3 < |$y_{CM}$| < 0.8

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0.3 < |$y_{CM}$| < 0.8

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0.3 < |$y_{CM}$| < 0.8

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0.8 < |$y_{CM}$| < 1.3

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0.8 < |$y_{CM}$| < 1.3

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 0.8 < |$y_{CM}$| < 1.3

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 1.3 < |$y_{CM}$| < 1.8

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 1.3 < |$y_{CM}$| < 1.8

$Y_{asym}$ in p+Pb at $\sqrt{s}$=5.02 TeV in 1.3 < |$y_{CM}$| < 1.8

Results of elliptic flow, corrected for short range correlations, for prompt neutral D mesons, as a function of multiplicity for $|y_{lab}|< 1$, with 6$ < p_{T} < $8 GeV in pp collisions at 13 TeV.

Results of elliptic flow, corrected for short range correlations, for prompt neutral D mesons, as a function of multiplicity for $|y_{lab}|< 1$, with 2$ < p_{T} < $4 GeV in pPb collisions at 8.16 TeV.

Results of elliptic flow, corrected for short range correlations, for prompt neutral D mesons, as a function of multiplicity for $|y_{lab}|< 1$, with 4$ < p_{T} < $6 GeV in pPb collisions at 8.16 TeV.

Results of elliptic flow, corrected for short range correlations, for prompt neutral D mesons, as a function of multiplicity for $|y_{lab}|< 1$, with 6$ < p_{T} < $8 GeV in pPb collisions at 8.16 TeV.

Results of elliptic flow, corrected for short range correlations, for nonprompt neutral D mesons, as a function of transverse momenta for $|y_{lab}|< 1$, with $185 \leq N^{offline}_{trk} < 250$ in pPb collisions at 8.16 TeV.

Prompt $D^0$ meson $v_3$ as a function of $p_T$ in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_3$ as a function of $p_T$ in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_3$ as a function of $p_T$ in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_2$ as a function of $p_T$ in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_2$ as a function of $p_T$ in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_2$ as a function of $p_T$ in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_3$ as a function of $p_T$ in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_3$ as a function of $p_T$ in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_3$ as a function of $p_T$ in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_2$ as a function of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_2$ as a function of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_3$ as a function of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_3$ as a function of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_2$ as a function of rapidity in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $v_3$ as a function of rapidity in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

Prompt $D^0$ meson $\Delta v_2$ as a function of rapidity in PbPb collisions at $\sqrt{s_{NN}}=5.02~TeV$.

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)

The intercept parameter as a function of kT using the HCS method in pp collisions at 13 TeV.

1/Rinv^2 as a function of mT using the HCS method in pp collisions at 13 TeV.

The depth of the anticorrelation as function of multiplicity for kT integrated values.

Correlation function from the double ratio technique, integrated in the range 0 < kT < 1 GeV, for 0=<Ntrkoffline<=4. The error bars represent statistical uncertainties.

Correlation function from the double ratio technique, integrated in the range 0 < kT < 1 GeV, for 10=<Ntrkoffline<=12. The error bars represent statistical uncertainties.

Correlation function from the double ratio technique, integrated in the range 0 < kT < 1 GeV, for 31=<Ntrkoffline<=33. The error bars represent statistical uncertainties.

Correlation function from the double ratio technique, integrated in the range 0 < kT < 1 GeV, for 80=<Ntrkoffline<=84. The error bars represent statistical uncertainties.

Correlation function from the double ratio technique, integrated in the range 0 < kT < 1 GeV, for 105=<Ntrkoffline<=109. The error bars represent statistical uncertainties.

Correlation function from the double ratio technique, integrated in the range 0 < kT < 1 GeV, for 130=<Ntrkoffline<=250. The error bars represent statistical uncertainties.

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

Single vector leptoquark (LQv) k = 1 total selection efficiency, accounting for both the decay branching fraction and the event selection, for events that pass the signal region requirements and any of the top quark or b jet categories defined in the search.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 50-90 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in pp collisions. The results are shown in different dijet momentum balance bins.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 0-10 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $0.7 < p_{\mathrm{T}}^{\mathrm{ch}} < 1$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 10-30 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $0.7 < p_{\mathrm{T}}^{\mathrm{ch}} < 1$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 30-50 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $0.7 < p_{\mathrm{T}}^{\mathrm{ch}} < 1$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 50-90 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $0.7 < p_{\mathrm{T}}^{\mathrm{ch}} < 1$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in pp collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $0.7 < p_{\mathrm{T}}^{\mathrm{ch}} < 1$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 0-10 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $1 < p_{\mathrm{T}}^{\mathrm{ch}} < 2$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 10-30 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $1 < p_{\mathrm{T}}^{\mathrm{ch}} < 2$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 30-50 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $1 < p_{\mathrm{T}}^{\mathrm{ch}} < 2$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 50-90 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $1 < p_{\mathrm{T}}^{\mathrm{ch}} < 2$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in pp collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $1 < p_{\mathrm{T}}^{\mathrm{ch}} < 2$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 0-10 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $2 < p_{\mathrm{T}}^{\mathrm{ch}} < 3$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 10-30 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $2 < p_{\mathrm{T}}^{\mathrm{ch}} < 3$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 30-50 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $2 < p_{\mathrm{T}}^{\mathrm{ch}} < 3$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 50-90 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $2 < p_{\mathrm{T}}^{\mathrm{ch}} < 3$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in pp collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $2 < p_{\mathrm{T}}^{\mathrm{ch}} < 3$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 0-10 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $3 < p_{\mathrm{T}}^{\mathrm{ch}} < 4$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 10-30 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $3 < p_{\mathrm{T}}^{\mathrm{ch}} < 4$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 30-50 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $3 < p_{\mathrm{T}}^{\mathrm{ch}} < 4$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 50-90 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $3 < p_{\mathrm{T}}^{\mathrm{ch}} < 4$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in pp collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $3 < p_{\mathrm{T}}^{\mathrm{ch}} < 4$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 0-10 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $4 < p_{\mathrm{T}}^{\mathrm{ch}} < 8$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 10-30 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $4 < p_{\mathrm{T}}^{\mathrm{ch}} < 8$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 30-50 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $4 < p_{\mathrm{T}}^{\mathrm{ch}} < 8$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 50-90 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $4 < p_{\mathrm{T}}^{\mathrm{ch}} < 8$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in pp collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $4 < p_{\mathrm{T}}^{\mathrm{ch}} < 8$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 0-10 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $8 < p_{\mathrm{T}}^{\mathrm{ch}} < 12$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 10-30 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $8 < p_{\mathrm{T}}^{\mathrm{ch}} < 12$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 30-50 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $8 < p_{\mathrm{T}}^{\mathrm{ch}} < 12$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in the 50-90 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $8 < p_{\mathrm{T}}^{\mathrm{ch}} < 12$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the leading jets as a function of $\Delta\eta$ in pp collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $8 < p_{\mathrm{T}}^{\mathrm{ch}} < 12$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 0-10 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 10-30 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 30-50 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 50-90 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in pp collisions. The results are shown in different dijet momentum balance bins.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 0-10 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $0.7 < p_{\mathrm{T}}^{\mathrm{ch}} < 1$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 10-30 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $0.7 < p_{\mathrm{T}}^{\mathrm{ch}} < 1$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 30-50 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $0.7 < p_{\mathrm{T}}^{\mathrm{ch}} < 1$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 50-90 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $0.7 < p_{\mathrm{T}}^{\mathrm{ch}} < 1$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in pp collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $0.7 < p_{\mathrm{T}}^{\mathrm{ch}} < 1$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 0-10 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $1 < p_{\mathrm{T}}^{\mathrm{ch}} < 2$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 10-30 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $1 < p_{\mathrm{T}}^{\mathrm{ch}} < 2$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 30-50 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $1 < p_{\mathrm{T}}^{\mathrm{ch}} < 2$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 50-90 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $1 < p_{\mathrm{T}}^{\mathrm{ch}} < 2$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in pp collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $1 < p_{\mathrm{T}}^{\mathrm{ch}} < 2$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 0-10 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $2 < p_{\mathrm{T}}^{\mathrm{ch}} < 3$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 10-30 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $2 < p_{\mathrm{T}}^{\mathrm{ch}} < 3$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 30-50 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $2 < p_{\mathrm{T}}^{\mathrm{ch}} < 3$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 50-90 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $2 < p_{\mathrm{T}}^{\mathrm{ch}} < 3$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in pp collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $2 < p_{\mathrm{T}}^{\mathrm{ch}} < 3$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 0-10 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $3 < p_{\mathrm{T}}^{\mathrm{ch}} < 4$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 10-30 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $3 < p_{\mathrm{T}}^{\mathrm{ch}} < 4$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 30-50 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $3 < p_{\mathrm{T}}^{\mathrm{ch}} < 4$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 50-90 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $3 < p_{\mathrm{T}}^{\mathrm{ch}} < 4$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in pp collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $3 < p_{\mathrm{T}}^{\mathrm{ch}} < 4$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 0-10 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $4 < p_{\mathrm{T}}^{\mathrm{ch}} < 8$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 10-30 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $4 < p_{\mathrm{T}}^{\mathrm{ch}} < 8$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 30-50 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $4 < p_{\mathrm{T}}^{\mathrm{ch}} < 8$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 50-90 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $4 < p_{\mathrm{T}}^{\mathrm{ch}} < 8$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in pp collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $4 < p_{\mathrm{T}}^{\mathrm{ch}} < 8$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 0-10 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $8 < p_{\mathrm{T}}^{\mathrm{ch}} < 12$ GeV.

The distribution of charged particle yields within $|\Delta\varphi| < 1.0$ correlated with the subleading jets as a function of $\Delta\eta$ in the 10-30 % centrality bin in PbPb collisions. The results are shown in different dijet momentum balance bins for the charged particle $p_{\mathrm{T}}$ bin $8 < p_{\mathrm{T}}^{\mathrm{ch}} < 12$ GeV.