Measured natural width

The measured inclusive ratio of production cross sections

Summary of results for cross-section of non-prompt $\psi(2S)$ decaying to a muon pair for 13 TeV data in fb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $J/\psi$ decaying to a muon pair as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $\psi(2S)$ decaying to a muon pair as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Correction factors for an assumed angular dependence $\propto 1+\lambda_{\theta} \cos^2\theta $ in the helicity frame, for several values of the parameter $\lambda_{\theta}$. The correction factors were found to be the same (within 1% - 2%) for $J\psi$ and $\psi(2S)$ mesons, for prompt and non-prompt production mechanisms, and for the three rapidity intervals considered in this paper.

$p_{\rm T}$-distributions of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in pp collisions at \sqrt{s}$ = 5.02 TeV for 10-20\% multiplicity class.

$p_{\rm T}$-distributions of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in pp collisions at \sqrt{s}$ = 5.02 TeV for 20-30\% multiplicity class.

$p_{\rm T}$-distributions of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in pp collisions at \sqrt{s}$ = 5.02 TeV for 30-40\% multiplicity class.

$p_{\rm T}$-distributions of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in pp collisions at \sqrt{s}$ = 5.02 TeV for 40-50\% multiplicity class.

$p_{\rm T}$-distributions of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in pp collisions at \sqrt{s}$ = 5.02 TeV for 50-70\% multiplicity class.

$p_{\rm T}$-distributions of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in pp collisions at \sqrt{s}$ = 5.02 TeV for 70-100\% multiplicity class.

$p_{\rm T}$-distributions of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in pp collisions at \sqrt{s}$ = 5.02 TeV for 0-100\% multiplicity class.

$p_{\rm T}$-distributions of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in Xe Xe collisions at \sqrt{s_{NN}}$ = 5.44 TeV for 0-30\% centrality class.

$p_{\rm T}$-distributions of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in Xe Xe collisions at \sqrt{s_{NN}}$ = 5.44 TeV for 30-50\% centrality class.

$p_{\rm T}$-distributions of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in Xe Xe collisions at \sqrt{s_{NN}}$ = 5.44 TeV for 50-70\% centrality class.

$p_{\rm T}$-distributions of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in Xe Xe collisions at \sqrt{s_{NN}}$ = 5.44 TeV for 70-90\% centrality class.

$p_{\rm T}$-integrated yield d$N$/d$y$ of K$^{*}$ (average of particle and anti-particle) meson measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV.

Mean $p_{\rm T}$ of K$^{*}$ (average of particle and anti-particle) meson measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV.

$p_{\rm T}$-integrated yield d$N$/d$y$ of K$^{*}$ (average of particle and anti-particle) meson measured in pp collisions at $\sqrt{s}$ = 5.02 TeV.

Mean $p_{\rm T}$ of K$^{*}$ (average of particle and anti-particle) meson measured in pp collisions at $\sqrt{s}$ = 5.02 TeV.

$p_{\rm T}$-integrated yield ratio of $\rm{K}^{*}$ (average of particleand anti-particle) to K (average of particle and anti-particle) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV.

$p_{\rm T}$-integrated yield ratio of $\rm{K}^{*}$ (average of particleand anti-particle) to K (average of particle and anti-particle) measured in pp collisions at $\sqrt{s}$ = 5.02 TeV.

Lower limit of the hadronic phase lifetime measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV.

Lower limit of the hadronic phase lifetime measured in pp collisions at $\sqrt{s}$ = 5.02 TeV.

Kinetic freeze out temperature measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV.

$p_{\rm T}$ differential ratio $\rm{K}^{*}$ (average of particleand anti-particle) to K (average of particle and anti-particle) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV for $p_{\rm T}$ range 0.4-2.0 GeV/$\it{c}.

$p_{\rm T}$ differential ratio $\rm{K}^{*}$ (average of particleand anti-particle) to K (average of particle and anti-particle) measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV for $p_{\rm T}$ range 2.0-4.0 GeV/$\it{c}.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV for 0-30\% centrality.

RAA of $\rm{K}^{*}$ (average of particle and anti-particle) meson as a function of system size measured in Xe-Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV for $p_{\rm T}$ range 4.0-12.0 GeV/$\it{c}.

$p_{\rm T}$-distributions of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at \sqrt{s_{NN}}$ = 5.02 TeV for 40-60\% centrality.

$p_{\rm T}$-distributions of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at \sqrt{s_{NN}}$ = 5.02 TeV for 60-80\% centrality.

$p_{\rm T}$-integrated yield d$N$/d$y$ of K$^{*}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

Mean $p_{\rm T}$ of K$^{*}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{\rm T}$-integrated yield ratio of $\rm{K}^{*}$ (average of particleand anti-particle) to K (average of particle and anti-particle) measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{\rm T}$ differential ratio $(\rm{K}^{*\pm} + \rm{K}^{*\mp}) / (\rm{K}^{+}+\rm{K}^{-})$ measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for 0-10\% centrality interval.

$p_{\rm T}$ differential ratio $(\rm{K}^{*\pm} + \rm{K}^{*\mp}) / (\rm{K}^{+}+\rm{K}^{-})$ measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for 60-80\% centrality interval.

$p_{\rm T}$ differential ratio $(\rm{K}^{*\pm} + \rm{K}^{*\mp}) / (\rm{pi}^{+}+\rm{pi}^{-})$ measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for 0-10\% centrality interval.

$p_{\rm T}$ differential ratio $(\rm{K}^{*\pm} + \rm{K}^{*\mp}) / (\rm{pi}^{+}+\rm{pi}^{-})$ measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for 60-80\% centrality interval.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for 0-10\% centrality.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for 10-20\% centrality.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for 20-40\% centrality.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for 40-60\% centrality.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for 60-80\% centrality.

Measured ratio of $\mathcal{B}_{4\mu}/\mathcal{B}_{2\mu}$

Measured branching fraction $\mathcal{B}_{4\mu}$

Measured branching fraction $\mathcal{B}_{4\mu}$

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ and photon isolation cone radius $R=0.4$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ and photon isolation cone radius $R=0.4$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ and photon isolation cone radius $R=0.4$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ and photon isolation cone radius $R=0.2$.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ and isolation cone radius $0.2$ at NNLO QCD.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ at NNLO QCD.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $|\eta^{\gamma}|<0.6$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $0.6<|\eta^{\gamma}|<0.8$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $0.8<|\eta^{\gamma}|<1.37$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $1.56<|\eta^{\gamma}|<1.81$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $1.81<|\eta^{\gamma}|<2.01$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $2.01<|\eta^{\gamma}|<2.37$.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $|\eta^{\gamma}|<0.6$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $0.6<|\eta^{\gamma}|<0.8$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $0.8<|\eta^{\gamma}|<1.37$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $1.56<|\eta^{\gamma}|<1.81$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $1.81<|\eta^{\gamma}|<2.01$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $2.01<|\eta^{\gamma}|<2.37$ at NNLO QCD.

The profile likelihood scan as a function of Bs to mu+mu- and B0 to mu+mu- decay branching fractions in 2D.

Left pannel. The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 2.3%) is not graphically represented. The blue line represents the average for $p_T > 18$ GeV. For comparison, the LHCb measurement [10.1103/PhysRevLett.124.122002] is also shown. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (2.3%)

Left pannel. The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 2.3%) is not graphically represented. The blue line represents the average for $p_T > 18$ GeV. For comparison, the LHCb measurement [10.1103/PhysRevLett.124.122002] is also shown. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (2.3%)

Left pannel. The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 2.3%) is not graphically represented. The blue line represents the average for $p_T > 18$ GeV. For comparison, the LHCb measurement [10.1103/PhysRevLett.124.122002] is also shown. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (2.3%)

Right pannel. The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 2.3%) is not graphically represented. The blue line represents the average for $p_T > 18$ GeV. For comparison, the LHCb measurement [10.1103/PhysRevLett.124.122002] is also shown. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (2.3%)

Right pannel. The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 2.3%) is not graphically represented. The blue line represents the average for $p_T > 18$ GeV. For comparison, the LHCb measurement [10.1103/PhysRevLett.124.122002] is also shown. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (2.3%)

Right pannel. The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 2.3%) is not graphically represented. The blue line represents the average for $p_T > 18$ GeV. For comparison, the LHCb measurement [10.1103/PhysRevLett.124.122002] is also shown. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (2.3%)

Left panel.The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 4.9\%) is not graphically represented. The blue line represents the average of all the points. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (4.9%)

Left panel.The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 5.7%) is not graphically represented. The blue line represents the average of all the points. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (5.7%)

Left panel.The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 5.7%) is not graphically represented. The blue line represents the average of all the points. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (5.7%)

Right panel.The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 4.9\%) is not graphically represented. The blue line represents the average of all the points. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (4.9%)

Right panel.The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 5.7%) is not graphically represented. The blue line represents the average of all the points. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (5.7%)

Right panel.The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 5.7%) is not graphically represented. The blue line represents the average of all the points. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (5.7%)

Fraction of non-prompt J/$\psi$ in p--Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 8.16 TeV as a function of $p_\mathrm{T}$.

d$^2\sigma$/d$y$d$p_{\rm T}$ in bins of $p_{\mathrm{T}}^{J/\psi}$ for prompt J/$\psi$ in p--Pb collisions at $\sqrt{s_{NN}}$ = 8.16 TeV.

d$^2\sigma$/d$y$d$p_{\rm T}$ in bins of $p_{\mathrm{T}}^{J/\psi}$ for non-prompt J/$\psi$ in p--Pb collisions at $\sqrt{s_{NN}}$ = 8.16 TeV.

Nuclear modification factor ($R_{pPb}$) of prompt J/$\psi$ in p--Pb collisions at $\sqrt{s_{NN}}$ = 8.16 TeV at midrapidity.

Nuclear modification factor ($R_{pPb}$) of non-prompt J/$\psi$ in p--Pb collisions at $\sqrt{s_{NN}}$ = 8.16 TeV at midrapidity.

$p_{\rm T}$-differential density of $\Xi^{\pm}$ in jets and the underlying event in pp collisions at $\sqrt{s} = 13$ TeV.

$p_{\rm T}$-differential density of inclusive $\Omega^{\pm}$ in pp collisions at $\sqrt{s} = 13$ TeV.

$p_{\rm T}$-differential density of $\Omega^{\pm}$ in jets and the underlying event in pp collisions at $\sqrt{s} = 13$ TeV.

$p_{\rm T}$-differential $(\Lambda + \overline{\Lambda})/2{\rm K}_{\rm S}^{0}$ ratio of inclusive particles in pp collisions at $\sqrt{s} = 13$ TeV.

$p_{\rm T}$-differential $(\Lambda + \overline{\Lambda})/2{\rm K}_{\rm S}^{0}$ ratio in jets and the underlying event in pp collisions at $\sqrt{s} = 13$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / 2{\rm K}_{\rm S}^{0}$ and $(\Xi^{-} + \Xi^{+}) / (\Lambda + \overline{\Lambda})$ ratios of inclusive particles in pp collisions at $\sqrt{s} = 13$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / 2{\rm K}_{\rm S}^{0}$ and $(\Xi^{-} + \Xi^{+}) / (\Lambda + \overline{\Lambda})$ ratios in jets and the underlying event in pp collisions at $\sqrt{s} = 13$ TeV.

$p_{\rm T}$-differential $(\Omega^{-} + \Omega^{+}) / 2{\rm K}_{\rm S}^{0}$, $(\Omega^{-} + \Omega^{+}) / (\Lambda + \overline{\Lambda})$ and $(\Omega^{-} + \Omega^{+}) / (\Xi^{-} + \Xi^{+})$ ratios of inclusive particles in pp collisions at $\sqrt{s} = 13$ TeV.

$p_{\rm T}$-differential $(\Omega^{-} + \Omega^{+}) / 2{\rm K}_{\rm S}^{0}$, $(\Omega^{-} + \Omega^{+}) / (\Lambda + \overline{\Lambda})$ and $(\Omega^{-} + \Omega^{+}) / (\Xi^{-} + \Xi^{+})$ ratios in jets and the underlying event in pp collisions at $\sqrt{s} = 13$ TeV.

$p_{\rm T}$-differential density of inclusive ${\rm K}_{\rm S}^{0}$ and $\Lambda$ ($\overline{\Lambda}$) in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential densities of ${\rm K}_{\rm S}^{0}$ and $\Lambda$ ($\overline{\Lambda}$) in jets and the underlying event in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of inclusive $\Xi^{\pm}$ in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of $\Xi^{\pm}$ in jets and the underlying event in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of inclusive $\Omega^{\pm}$ in p-Pb collisions at $\sqrt{s} = 5.02$ TeV.

$p_{\rm T}$-differential density of $\Omega^{\pm}$ in jets and the underlying event in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Lambda + \overline{\Lambda})/2{\rm K}_{\rm S}^{0}$ ratio in jets for various event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / 2{\rm K}_{\rm S}^{0}$ and $(\Xi^{-} + \Xi^{+}) / (\Lambda + \overline{\Lambda})$ ratios in jets for various event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Omega^{-} + \Omega^{+}) / 2{\rm K}_{\rm S}^{0}$, $(\Omega^{-} + \Omega^{+}) / (\Lambda + \overline{\Lambda})$ and $(\Omega^{-} + \Omega^{+}) / (\Xi^{-} + \Xi^{+})$ ratios in jets in MB p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of inclusive ${\rm K}_{\rm S}^{0}$ and $\Lambda$ ($\overline{\Lambda}$) for $0$-$10\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential densities of ${\rm K}_{\rm S}^{0}$ and $\Lambda$ ($\overline{\Lambda}$) in jets and the underlying event for $0$-$10\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of inclusive $\Xi^{\pm}$ for $0$-$10\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of $\Xi^{\pm}$ in jets and the underlying event for $0$-$10\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of inclusive ${\rm K}_{\rm S}^{0}$ and $\Lambda$ ($\overline{\Lambda}$) for $10$-$40\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential densities of ${\rm K}_{\rm S}^{0}$ and $\Lambda$ ($\overline{\Lambda}$) in jets and the underlying event for $10$-$40\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of inclusive $\Xi^{\pm}$ for $10$-$40\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of $\Xi^{\pm}$ in jets and the underlying event for $10$-$40\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of inclusive ${\rm K}_{\rm S}^{0}$ and $\Lambda$ ($\overline{\Lambda}$) for $40$-$100\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential densities of ${\rm K}_{\rm S}^{0}$ and $\Lambda$ ($\overline{\Lambda}$) in jets and the underlying event for $40$-$100\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of inclusive $\Xi^{\pm}$ for $40$-$100\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential density of $\Xi^{\pm}$ in jets and the underlying event for $40$-$100\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Lambda + \overline{\Lambda})/2{\rm K}_{\rm S}^{0}$ ratio of inclusive particles for various event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Lambda + \overline{\Lambda})/2{\rm K}_{\rm S}^{0}$ ratio in the underlying event for various event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Lambda + \overline{\Lambda})/2{\rm K}_{\rm S}^{0}$ ratio in jets for various event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / 2{\rm K}_{\rm S}^{0}$ ratio of inclusive particles for various event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / 2{\rm K}_{\rm S}^{0}$ ratio in the underlying event for $0$-$10\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / 2{\rm K}_{\rm S}^{0}$ ratio in the underlying event for $10$-$40\%$ and $40$-$100\%$ event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / 2{\rm K}_{\rm S}^{0}$ ratio in jets for various event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / (\Lambda + \overline{\Lambda})$ ratio of inclusive particles for various event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / (\Lambda + \overline{\Lambda})$ ratio in the underlying event for $0$-$10\%$ event activity class in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / (\Lambda + \overline{\Lambda})$ ratio in the underlying event for $10$-$40\%$ and $40$-$100\%$ event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

$p_{\rm T}$-differential $(\Xi^{-} + \Xi^{+}) / (\Lambda + \overline{\Lambda})$ ratio in jets for various event activity classes in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV.

Summary of the total uncertainty in the background prediction for each SR of the tt0L-low, tt0L-high, tt1L and tt2L analysis channels in the statistical combination. Their dominant contributions are indicated by individual lines. Individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

$E_{\text{T}}^{\text{miss}}$ distribution in SR0X for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

$E_{\text{T}}^{\text{miss}}$ distribution in SRWX for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

$E_{\text{T}}^{\text{miss}}$ distribution in SRTX for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Representative fit distribution in the different flavour leptons signal region for the tt2L analysis: each bin of such distribution, starting from the red arrow, corresponds to a single SR included in the fit. 'FNP' includes the contribution from fake/non-prompt lepton background arising from jets (mainly $\pi/K$, heavy-flavour hadron decays and photon conversion) misidentified as leptons, estimated in a purely data-driven way. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The total uncertainty in the SM expectation is represented with hatched bands and the expected distributions for selected signal models are shown as dashed lines.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$t\bar{t}$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$tW$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$tj$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$t\bar{t}$ model, defined as the number of selected reconstructed events divided by the acceptance.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$tW$ model, defined as the number of selected reconstructed events divided by the acceptance.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$tj$ model, defined as the number of selected reconstructed events divided by the acceptance.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Distribution of ${E}_{T}^{miss}$ after the fit of the multijet backgrounds, in the electron channel, in the $t\bar{t}$ 3-jets VR, without applying the cut on ${E}_{T}^{miss}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of ${E}_{T}^{miss}$ after the fit of the multijet backgrounds, in the electron channel, in the $t\bar{t}$ 4-jets VR, without applying the cut on ${E}_{T}^{miss}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the signal region, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $W$+jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $t\bar{t}$ 3-jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $t\bar{t}$ 4-jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Expected distributions of the MEM discriminant $P(S|X)$ in the SR, for the s-channel single-top signal, and for the $t\bar{t}$ and $W$+jets backgrounds, for MEM discriminant values larger than $2.0\times10^{-4}$. Each distribution is normalised to unity. The binning is the same as the optimised binning used in the signal extraction fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $W$+jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty band indicates the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $t\bar{t}$ 3-jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty band indicates the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $t\bar{t}$ 4-jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty bands indicate the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR before the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$. The lower panel of the figure shows the ratio of the data to the prediction, with different vertical axis ranges. The uncertainty band indicates the total uncertainties and their correlations in each bin. The uncertainties in the $t\bar{t}$ and $W$+jets normalisation factors, as well as in the s-channel signal cross-section, are not defined pre-fit and therefore not included. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR after the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$. The lower panel of the figure shows the ratio of the data to the prediction, with different vertical axis ranges. The uncertainty band indicates the total uncertainties and their correlations in each bin. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR after the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$, after subtraction of all backgrounds. The fitted distribution for the simulation of the signal is shown together with the post-fit uncertainty in the backgrounds. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Pre-fit and post-fit event yields in the SR, for MEM discriminant values larger than $2.0\times10^{-4}$. The central value of the event yield for each process is calculated by summing the values of the discriminant bin contents, using the nominal expected yield for the pre-fit value, and the best-fit estimate for the post-fit value. The error includes statistical and systematic uncertainties summed in quadrature. All sources of systematic uncertainties are included, taking into account correlations and anti-correlations in the post-fit case. The uncertainties in the $t\bar{t}$ and $W$+jets normalisation factors, as well as in the s-channel signal cross-section, are not defined pre-fit and therefore only included in the post-fit uncertainties.

Observed impact of the different sources of uncertainty on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. The 'Systematic uncertainties' category combines all sources of systematic uncertainties. The statistical uncertainty is obtained by repeating the fit after having fixed all nuisance parameters, including the $t\bar{t}$ and $W$+jets normalisation factors. 'Total' gives the total uncertainty on the measurement.

Observed impact of the different sources of $t\bar{t}$ modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, and 'ME/PS matching' to the matching of the ME to the parton shower.

Observed impact of the different sources of s-channel modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, as described in Section 7 in the paper.

Observed impact of the different sources of t-channel modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, as described in Section 7 in the paper.

Observed impact of the different sources of $tW$ modelling uncertainty on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, and '$t\bar{t}$ overlap' to the algorithm removing the overlap between $tW$ and $t\bar{t}$ production at NLO, as described in Section 7 in the paper.

Observed impact of the different sources of PDF uncertainties on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root.

Comparison between data and prediction after the fit to data in the signal region for the leading-jet $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the leading-jet $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the subleading-jet $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the subleading-jet $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the lepton $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the lepton $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the ${E}_{T}^{miss}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the $m_{T}^{W}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Nuisance parameters ranked according to their post-fit impacts on the best-fit value of the ratio $\mu$ of the measured cross-section to the predicted cross-section. In the figure, only the 20 nuisance parameters with the largest post-fit impacts are shown. The empty (solid) blue rectangles illustrate the pre-fit (post-fit) impact on $\mu$, corresponding to the upper axis. The pre-fit (post-fit) impact of each nuisance parameter, $\Delta\mu$, is calculated as the difference in the fitted value of $\mu$ between the nominal fit and the fit when fixing the corresponding nuisance parameter to $\hat{\theta}\pm\Delta\theta$ ($\hat{\theta}\pm\Delta\hat{\theta}$), where $\hat{\theta}$ is the best-fit value of the nuisance parameter and $\Delta\theta$ ($\Delta\hat{\theta}$) is its pre-fit (post-fit) uncertainty. Several systematic uncertainties are split into different nuisance parameters, which are indicated by NP. JES (JER) indicates jet energy scale (resolution), and $\gamma$ indicates a nuisance parameter associated to the MC statistics in one of the 18 bins numbered from 0 to 17. The black points show the best-fit values of the nuisance parameters, with the error bars representing the post-fit uncertainties. Each nuisance parameter is shown wrt. its nominal value, $\theta_0$, and in units of its pre-fit uncertainty, except the free-floating normalisation factors of the $t\bar{t}$ and $W$+jets backgrounds, and the parameters associated to the MC statistics in each bin, for which the post-fit values and uncertainties are shown.

Correlation matrix of the nuisance parameters and of the ratio $\mu$ of the measured cross-section to the predicted cross-section. The correlations are given after the fit to data. In the figure, only the parameters which have a correlation of at least 0.2 with any other parameter are shown.

Distribution of the MEM discriminant $P(S|X)$ in the SR for MEM discriminant values larger than $2.0\times10^{-4}$, for the collision data used for the measurement, and for 1000 pseudo-data replicas, generated using a bootstrapping technique, in order to assess the statistical correlations between this measurement and others, for the purpose of combinations. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the <a href="https://zenodo.org/record/5361038">BootstrapGenerator</a> software package , which implements a technique described in <a href="https://cds.cern.ch/record/2759945/">ATL-PHYS-PUB-2021-011</a>. The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. Each pseudo-data replica is assigned an index, ranging from 0 to 999, corresponding to the random number index used consistently for each observed data event.

Measured values of the signal cross-section and of the $t\bar{t}$ and $W$+jets normalisation factors, obtained by statistical-only fits to the collision data used for the measurement, and to 1000 pseudo-data replicas, generated using a bootstrapping technique, in order to assess the statistical correlations between this measurement and others, for the purpose of combinations. The central values and their statistical uncertainties are obtained by repeating the fit after having fixed all nuisance parameters, except the $t\bar{t}$ and $W$+jets normalisation factors, which are let free-floating (unlike for the statistical uncertainty on the cross-section quoted in the paper). The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the <a href="https://zenodo.org/record/5361038">BootstrapGenerator</a> software package , which implements a technique described in <a href="https://cds.cern.ch/record/2759945/">ATL-PHYS-PUB-2021-011</a>. The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. Each pseudo-data replica is assigned an index, ranging from 0 to 999, corresponding to the random number index used consistently for each observed data event.

$p_{\rm{T}}$-differential yield of $\Sigma^{*-}$ + cc in Pb-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (0-10% multiplicity class).

$p_{\rm{T}}$-differential yield of $\Sigma^{*-}$ + cc in Pb-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (30-50% multiplicity class).

$p_{\rm{T}}$-differential yield of $\Sigma^{*-}$ + cc in Pb-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (50-90% multiplicity class).

$p_{\rm T}$-integrated yield d$N$/d$y$ of ($\Sigma^{*+}+\overline{\Sigma}^{*-}$) in Pb-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV for each V0A multiplicity event class.

$p_{\rm T}$-integrated yield d$N$/d$y$ of ($\Sigma^{*-}+\overline{\Sigma}^{*+}$) in Pb-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV for each V0A multiplicity event class.

Mean $p_{\rm T}$ of ($\Sigma^{*+}+\overline{\Sigma}^{*-}$) in Pb-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV for each V0A multiplicity event class.

Mean $p_{\rm T}$ of ($\Sigma^{*-}+\overline{\Sigma}^{*+}$) in Pb-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV for each V0A multiplicity event class.

Integrated yield ratio.

The measured branching fraction ratios

Transverse momentum distribution for B-candidate

Transverse momentum distribution for B-candidate

Transverse momentum distribution for B-candidate

Transverse momentum distribution for B-candidate

The measured branching fractions

The measured branching fractions

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Comparison of predicted and observed event yields in each of the control and signal regions in the dilepton channel after the fit to the data. The uncertainty band includes all uncertainties and their correlations.

Comparison of predicted and observed event yields in each of the control and signal regions in the single-lepton channels after the fit to the data. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $120\le p_T^H<200$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $200\le p_T^H<300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $120\le p_T^H<200$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $200\le p_T^H<300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 450$ GeV (yield only). The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton boosted SRs after the inclusive fit to the data for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton boosted SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for ${\Delta R^{{avg}}_{bb}}$ after the inclusive fit to the data in the single-lepton $CR^{5j}_{{\geq}4b\ lo}$ control region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for ${\Delta R^{{avg}}_{bb}}$ after the inclusive fit to the data in the single-lepton $CR^{5j}_{{\geq}4b\ hi}$ control region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit yields of signal ($S$) and total background ($B$) as a function of $\log (S/B)$, compared with data. Final-discriminant bins in all dilepton and single-lepton analysis regions are combined into bins of $\log (S/B)$, with the signal normalised to the SM prediction used for the computation of $\log (S/B)$. The signal is then shown normalised to the best-fit value and the SM prediction. The lower frame reports the ratio of data to background, and this is compared with the expected ${t\bar {t}H}$-signal-plus-background yield divided by the background-only yield for the best-fit signal strength (solid red line) and the SM prediction (dashed orange line).

Comparison between data and prediction for the reconstruction BDT score for the Higgs boson candidate identified using Higgs boson information, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the average $\Delta \eta $ between $b$-tagged jets, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the likelihood discriminant, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the average $\Delta R$ for all possible combinations of $b$-tagged jet pairs, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the DNN $P(H)$ output for the Higgs boson candidate after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the DNN $P(H)$ output for the Higgs boson candidate after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the reconstructed Higgs boson candidate mass for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit distribution of the reconstructed Higgs boson candidate mass for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit distribution of the reconstructed Higgs boson candidate mass for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Fitted values of the ${t\bar {t}H}$ signal strength parameter in the individual channels and in the inclusive signal-strength measurement.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the fit. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Pre-fit distribution of the number of jets in the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Pre-fit distribution of the number of jets in the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Pre-fit distribution of the number of jets in the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Post-fit distribution of the number of jets in the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the number of jets in the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the number of jets in the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Signal-strength measurements in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive signal strength.

95% CL simplified template cross-section upper limits in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive limit. The observed limits are shown (solid black lines), together with the expected limits both in the background-only hypothesis (dotted black lines) and in the SM hypothesis (dotted red lines). In the case of the expected limits in the background-only hypothesis, one- and two-standard-deviation uncertainty bands are also shown. The hatched uncertainty bands correspond to the theory uncertainty in the fiducial cross-section prediction in each bin.

The ratios $S/B$ (black solid line, referring to the vertical axis on the left) and $S/\sqrt{B}$ (red dashed line, referring to the vertical axis on the right) for each category in the inclusive analysis in the dilepton channel (left) and in the single-lepton channels (right), where $S$ ($B$) is the number of selected signal (background) events predicted by the simulation and normalised to a luminosity of 139 fb$^{-1}$ .

Comparison between data and prediction for the $\Delta R$ between the Higgs candidate and the ${t\bar {t}}$ candidate system, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the number of $b$-tagged jet pairs with an invariant mass within 30 GeV of 125 GeV, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the reconstruction BDT score for the Higgs boson candidate identified using Higgs boson information, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the $\Delta R$ between the two highest ${p_{{T}}}$ $b$-tagged jets, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the sum of $b$-tagging discriminants of jets from Higgs, hadronic top and leptonic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the sum of $b$-tagging discriminants of jets from Higgs, hadronic top and leptonic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the hadronic top candidate invariant mass, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the hadronic top candidate invariant mass, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the fraction of the sum of $b$-tagging discriminants of all jets not associated to the Higgs or hadronic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the fraction of the sum of $b$-tagging discriminants of all jets not associated to the Higgs or hadronic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $0\le {\hat{p}_{{T}}^{H}}<120$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $120\le {\hat{p}_{{T}}^{H}}<200$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $200\le {\hat{p}_{{T}}^{H}}<300$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $300\le {\hat{p}_{{T}}^{H}}<450$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for ${\hat{p}_{{T}}^{H}}\ge 450$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

95% confidence level upper limits on signal-strength measurements in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive signal-strength limit, after the fit used to extract multiple signal-strength parameters. The observed limits are shown (solid black lines), together with the expected limits both in the background-only hypothesis (dotted black lines) and in the SM hypothesis (dotted red lines). In the case of the expected limits in the background-only hypothesis, one- and two-standard-deviation uncertainty bands are also shown.

Post-fit correlation matrix (in percentages) between the $\mu $ values obtained in the STXS bins.

Performance of the Higgs boson reconstruction algorithms. For each row of `truth' ${\hat{p}_{{T}}^{H}}$, the matrix shows (in percentages) the fraction of Higgs boson candidates which are truth-matched to ${b\bar {b}}$ decays, with reconstructed $p_T^H$ in the various bins of the dilepton (left), single lepton resolved (middle) and boosted (right) channels.

Pre-fit event yields in the dilepton signal regions and control regions. All uncertainties are included except the $k({t\bar {t}+{\geq }1b})$ uncertainty that is not defined pre-fit. For the ${t\bar {t}H}$ signal, the pre-fit yield values correspond to the theoretical prediction and corresponding uncertainties. `Other sources' refers to s-channel, t-channel, $tW$, $tWZ$, $tZq$, $Z+$ jets and diboson events.

Post-fit event yields in the dilepton signal regions and control regions, after the inclusive fit in all channels. All uncertainties are included, taking into account correlations. For the ${t\bar {t}H}$ signal, the post-fit yield and uncertainties correspond to those in the inclusive signal-strength measurement. `Other sources' refers to s-channel, t-channel, $tW$, $tWZ$, $tZq$, $Z+$ jets and diboson events.

Pre-fit event yields in the single-lepton resolved and boosted signal regions and control regions. All uncertainties are included except the $k({t\bar {t}+{\geq }1b})$ uncertainty that is not defined pre-fit. For the ${t\bar {t}H}$ signal, the pre-fit yield values correspond to the theoretical prediction and corresponding uncertainties. `Other top sources' refers to s-channel, t-channel, $tWZ$ and $tZq$ events.

Post-fit event yields in the single-lepton resolved and boosted signal regions and control regions, after the inclusive fit in all channels. All uncertainties are included, taking into account correlations. For the ${t\bar {t}H}$ signal, the post-fit yield and uncertainties correspond to those in the inclusive signal-strength measurement. `Other top sources' refers to s-channel, t-channel, $tWZ$ and $tZq$ events.

Breakdown of the contributions to the uncertainties in $\mu$. The contributions from the different sources of uncertainty are evaluated after the fit. The $\Delta \mu $ values are obtained by repeating the fit after having fixed a certain set of nuisance parameters corresponding to a group of systematic uncertainties, and then evaluating $(\Delta \mu)^2$ by subtracting the resulting squared uncertainty of $\mu $ from its squared uncertainty found in the full fit. The same procedure is followed when quoting the effect of the ${t\bar {t}+{\geq }1b}$ normalisation. The total uncertainty is different from the sum in quadrature of the different components due to correlations between nuisance parameters existing in the fit.

Fraction (in percentages) of signal events, after SR and CR selections, originating from $b\bar {b}$, $WW$ and other remaining Higgs boson decay modes in the dilepton channel.

Fraction (in percentages) of signal events, after SR and CR selections, originating from $b\bar {b}$, $WW$ and other remaining Higgs boson decay modes in the single-lepton channels.

Predicted SM ${t\bar {t}H}$ cross-section in each of the five STXS ${\hat{p}_{{T}}^{H}}$ bins and signal acceptance times efficiency (including all event selection criteria) in each STXS bin as well as for the inclusive ${\hat{p}_{{T}}^{H}}$ range.

Number of expected signal events before the fit, after each selection requirement applied to enter the dilepton channel $SR^{\geq 4j}_{\geq 4b}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.

Number of expected signal events before the fit, after each selection requirement applied to enter the single-lepton channel resolved $SR^{\geq 6j}_{\geq 4b}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.

Number of expected signal events before the fit, after each selection requirement applied to enter the single-lepton channel boosted $SR_{boosted}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.

Comparison of the pT spectrum for the fourth-leading from data to different PYTHIA8 (P8),HERWIG++ (H++),and HERWIG7 (H7) tunes.

Comparison of the eta spectrum for the leading jet from data to different PYTHIA8 (P8),HERWIG++ (H++),and HERWIG7 (H7) tunes.

Comparison of the eta spectrum for the sub-leading from data to different PYTHIA8 (P8),HERWIG++ (H++),and HERWIG7 (H7) tunes.

Comparison of the eta spectrum for the third-leading from data to different PYTHIA8 (P8),HERWIG++ (H++),and HERWIG7 (H7) tunes.

Comparison of the eta spectrum for the fourth-leading from data to different PYTHIA8 (P8),HERWIG++ (H++),and HERWIG7 (H7) tunes.

Comparison of the DeltaPhiSoft distribution from data to different PYTHIA8 (P8),HERWIG++ (H++), and HERWIG7 (H7) tunes. All distributions have been normalized to regions where a reduced DPS contribution is expected.

Comparison of the DeltaPhiMin distribution from data to different PYTHIA8 (P8),HERWIG++ (H++), and HERWIG7 (H7) tunes. All distributions have been normalized to regions where a reduced DPS contribution is expected.

Comparison of the DeltaY distribution from data to different PYTHIA8 (P8),HERWIG++ (H++), and HERWIG7 (H7) tunes. All distributions have been normalized to regions where a reduced DPS contribution is expected.

Comparison of the phi_ij distribution from data to different PYTHIA8 (P8),HERWIG++ (H++), and HERWIG7 (H7) tunes. All distributions have been normalized to regions where a reduced DPS contribution is expected.

Comparison of the Deltap_{T,Soft} distribution from data to different PYTHIA 8 (P8), HERWIG ++ (H++), and HERWIG 7 (H7) tunes. All distributions have been normalized to regions where a reduced DPS contribution is expected.

Comparison of the DeltaS distribution from data to different PYTHIA 8 (P8), HERWIG ++ (H++), and HERWIG 7 (H7) tunes. All distributions have been normalized to regions where a reduced DPS contribution is expected.

Comparison of the pT spectrum for the leading jet from data with different KATIE(KT), MADGRAPH5aMC@NLO(MG5), and POWHEG(PW) models

Comparison of the pT spectrum for the sub-leading from data with different KATIE(KT), MADGRAPH5aMC@NLO(MG5), and POWHEG(PW) models

Comparison of the pT spectrum for the third-leading from data to different PYTHIA8 (P8),HERWIG++ (H++),and HERWIG7 (H7) tunes.

Comparison of the pT spectrum for the fourth-leading from data with different KATIE(KT), MADGRAPH5aMC@NLO(MG5), and POWHEG(PW) models

Comparison of the eta spectrum for the leading jet from data with different KATIE(KT), MADGRAPH5aMC@NLO(MG5), and POWHEG(PW) models

Comparison of the eta spectrum for the sub-leading from data with different KATIE(KT), MADGRAPH5aMC@NLO(MG5), and POWHEG(PW) models

Comparison of the eta spectrum for the third-leading from data with different KATIE(KT), MADGRAPH5aMC@NLO(MG5), and POWHEG(PW) models

Comparison of the eta spectrum for the fourth-leading from data with different KATIE(KT), MADGRAPH5aMC@NLO(MG5), and POWHEG(PW) models

Comparison of the DeltaPhiSoft distribution from data to differentKATIE(KT), MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations. All distributions have been normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the DeltaPhiMin, distribution from data to differentKATIE(KT), MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations. All distributions have been normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the DeltaY and distribution from data to differentKATIE(KT), MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations. All distributions have been normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the phi_ij distribution from data to differentKATIE(KT), MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations. All distributions have been normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the Deltap_{T,Soft} distribution from data to different KATIE(KT),MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations. All distributions havebeen normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the DeltaS distribution from data to different KATIE(KT),MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations. All distributions havebeen normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the pT spectrum for the leading jet from data with different SPS+DPS KATIE( KT) and PYTHIA8 (P8) models.

Comparison of the pT spectrum for the sub-leading from data with different SPS+DPS KATIE( KT) and PYTHIA8 (P8) models.

Comparison of the pT spectrum for the third-leading from data with different SPS+DPS KATIE( KT) and PYTHIA8 (P8) models.

Comparison of the pT spectrum for the fourth-leading from data with different SPS+DPS KATIE( KT) and PYTHIA8 (P8) models.

Comparison of the eta spectrum for the leading jet from data with different SPS+DPS KATIE( KT) and PYTHIA8 (P8) models.

Comparison of the eta spectrum for the sub-leading from data with different SPS+DPS KATIE( KT) and PYTHIA8 (P8) models.

Comparison of the eta spectrum for the third-leading from data with different SPS+DPS KATIE( KT) and PYTHIA8 (P8) models.

Comparison of the eta spectrum for the fourth-leading from data with different SPS+DPS KATIE( KT) and PYTHIA8 (P8) models.

Comparison of the DeltaPhiSoft distribution from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models. All distributions have been normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the DeltaPhiMin distribution from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models. All distributions have been normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the DeltaY distribution from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models. All distributions have been normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the phi_ij distribution from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models. All distributions have been normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the Deltap_{T,Soft} distributions from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models. All distributions havebeen normalized to regions where a reduced DPS sensitivity is expected.

Comparison of the DeltaS distributions from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models. All distributions have been normalized to regions where a reduced DPS contribution is expected.

The DeltaS distribution obtained from the mixed data sample compared to predictions from the pure DPS sample in PYTHIA 8 (P8) and KATIE (KT). The distributions are normalized to unity.

The results of the template fit for the POWHEG (PW) NLO 2 -> 2 model without the hard MPI removed. As the DeltaS distribution obtained from the mixed data sample carries a statistical and systematic uncertainty, so does the total fitted sample.

Comparison of the DeltaPhiSoft distribution from data to different PYTHIA8 (P8),HERWIG++ (H++), and HERWIG7 (H7) tunes.

Comparison of the DeltaPhiMin distribution from data to different PYTHIA8 (P8),HERWIG++ (H++), and HERWIG7 (H7) tunes.

Comparison of the DeltaY distribution from data to different PYTHIA8 (P8),HERWIG++ (H++), and HERWIG7 (H7) tunes.

Comparison of the phi_ij distribution from data to different PYTHIA8 (P8),HERWIG++ (H++), and HERWIG7 (H7) tunes.

Comparison of the Deltap_{T,Soft} distribution from data to different PYTHIA 8 (P8), HERWIG ++ (H++), and HERWIG 7 (H7) tunes.

Comparison of the DeltaS distribution from data to different PYTHIA 8 (P8), HERWIG ++ (H++), and HERWIG 7 (H7) tunes.

Comparison of the DeltaPhiSoft distribution from data to different KATIE(KT),MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations.

Comparison of the DeltaPhiMin distribution from data to different KATIE(KT),MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations.

Comparison of the DeltaY distribution from data to different KATIE(KT),MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations.

Comparison of the phi_ij distribution from data to different KATIE(KT),MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations.

Comparison of the Deltap_{T,Soft} distribution from data to different KATIE(KT),MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations.

Comparison of the DeltaS distribution from data to different KATIE(KT),MADGRAPH5aMC@NLO(MG5), andPOWHEG(PW) implementations.

Comparison of the DeltaPhiSoft distribution from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models.

Comparison of the DeltaPhiMin distribution from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models.

Comparison of the DeltaY distribution from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models.

Comparison of the phi_ij distribution from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models.

Comparison of the Deltap_{T,Soft} distribution from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models.

Comparison of the DeltaS distributions from data to different SPS+DPS KATIE (KT) and PYTHIA 8 (P8) models.

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.

Lower mass requirement for signal regions.

Lower mass requirement for signal regions.

Lower mass requirement for signal regions.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Upper cross-section limit in stau search.

Upper cross-section limit in stau search.

Upper cross-section limit in chargino search.

Upper cross-section limit in chargino search.

Lower mass limit as function of gluino lifetime.

Lower mass limit as function of gluino lifetime.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Candidate reconstruction efficiency for loose selection (LooseEff).

Candidate reconstruction efficiency for loose selection (LooseEff).

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

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

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

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

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

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

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

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

Predicted ratio of cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<1.37$.

Measured ratio of cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$.

Predicted ratio of cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$.

Measured ratio of cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.37$.

Predicted ratio of cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.37$.

Measured double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$.

Predicted double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$.

Measured double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<1.37$.

Predicted double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<1.37$.

Measured double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$.

Predicted double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$.

Measured double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.37$.

Predicted double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.37$.

v_{2}{EP} vs. p_{T} of average prompt D0, D+, D*+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the rapidity interval |y|<0.8 The first (sys) error is the systematic uncertainty from the other sources The second (sys) error is the systematic uncertainty from the B feed-down contribution.

v_{2}{EP} vs. p_{T} of prompt D0 mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the 60% small-q_{2}^{TPC} sample in the rapidity interval |y|<0.8 The first (sys) error is the systematic uncertainty from the other sources The second (sys) error is the systematic uncertainty from the B feed-down contribution.

v_{2}{EP} vs. p_{T} of prompt D0 mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the 20% large-q_{2}^{TPC} sample in the rapidity interval |y|<0.8 The first (sys) error is the systematic uncertainty from the other sources The second (sys) error is the systematic uncertainty from the B feed-down contribution.

v_{2}{EP} vs. p_{T} of prompt D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the 60% small-q_{2}^{TPC} sample in the rapidity interval |y|<0.8 The first (sys) error is the systematic uncertainty from the other sources The second (sys) error is the systematic uncertainty from the B feed-down contribution.

v_{2}{EP} vs. p_{T} of prompt D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the 20% large-q_{2}^{TPC} sample in the rapidity interval |y|<0.8 The first (sys) error is the systematic uncertainty from the other sources The second (sys) error is the systematic uncertainty from the B feed-down contribution.

v_{2}{EP} vs. p_{T} of prompt D0 mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the 60% small-q_{2}^{TPC} sample in the rapidity interval |y|<0.8 The first (sys) error is the systematic uncertainty from the other sources The second (sys) error is the systematic uncertainty from the B feed-down contribution.

v_{2}{EP} vs. p_{T} of prompt D0 mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the 20% large-q_{2}^{TPC} sample in the rapidity interval |y|<0.8 The first (sys) error is the systematic uncertainty from the other sources The second (sys) error is the systematic uncertainty from the B feed-down contribution.

v_{2}{EP} vs. p_{T} of prompt D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the 60% small-q_{2}^{TPC} sample in the rapidity interval |y|<0.8 The first (sys) error is the systematic uncertainty from the other sources The second (sys) error is the systematic uncertainty from the B feed-down contribution.

v_{2}{EP} vs. p_{T} of prompt D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the 20% large-q_{2}^{TPC} sample in the rapidity interval |y|<0.8 The first (sys) error is the systematic uncertainty from the other sources The second (sys) error is the systematic uncertainty from the B feed-down contribution.

v_{2}{EP} vs. p_{T} of average prompt D0, D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the 60% small-q_{2}^{TPC} sample in the rapidity interval |y|<0.8 The first (sys) error is the systematic uncertainty from the other sources The second (sys) error is the systematic uncertainty from the B feed-down contribution.

v_{2}{EP} vs. p_{T} of average prompt D0, D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the 20% large-q_{2}^{TPC} sample in the rapidity interval |y|<0.8 The first (sys) error is the systematic uncertainty from the other sources The second (sys) error is the systematic uncertainty from the B feed-down contribution.

v_{2}{EP} vs. p_{T} of average prompt D0, D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the 60% small-q_{2}^{TPC} sample in the rapidity interval |y|<0.8 The first (sys) error is the systematic uncertainty from the other sources The second (sys) error is the systematic uncertainty from the B feed-down contribution.

v_{2}{EP} vs. p_{T} of average prompt D0, D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the 20% large-q_{2}^{TPC} sample in the rapidity interval |y|<0.8 The first (sys) error is the systematic uncertainty from the other sources The second (sys) error is the systematic uncertainty from the B feed-down contribution.

Ratio of v_{2}{EP} in the 60% small-q_{2}^{TPC} sample to the unbiased sample vs. p_{T} of average prompt D0, D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the rapidity interval |y|<0.8.

Ratio of v_{2}{EP} in the 20% large-q_{2}^{TPC} sample to the unbiased sample vs. p_{T} of average prompt D0, D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the rapidity interval |y|<0.8.

Ratio of v_{2}{EP} in the 60% small-q_{2}^{TPC} sample to the unbiased sample vs. p_{T} of average prompt D0, D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the rapidity interval |y|<0.8.

Ratio of v_{2}{EP} in the 20% large-q_{2}^{TPC} sample to the unbiased sample vs. p_{T} of average prompt D0, D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the rapidity interval |y|<0.8.

v_{2}{EP} vs. p_{T} of average prompt D0, D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the 60% small-q_{2}^{V0A} sample in the rapidity interval |y|<0.8 The first (sys) error is the systematic uncertainty from the other sources The second (sys) error is the systematic uncertainty from the B feed-down contribution.

v_{2}{EP} vs. p_{T} of average prompt D0, D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the 20% large-q_{2}^{V0A} sample in the rapidity interval |y|<0.8 The first (sys) error is the systematic uncertainty from the other sources The second (sys) error is the systematic uncertainty from the B feed-down contribution.

v_{2}{EP} vs. p_{T} of average prompt D0, D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the 60% small-q_{2}^{V0A} sample in the rapidity interval |y|<0.8 The first (sys) error is the systematic uncertainty from the other sources The second (sys) error is the systematic uncertainty from the B feed-down contribution.

v_{2}{EP} vs. p_{T} of average prompt D0, D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the 20% large-q_{2}^{V0A} sample in the rapidity interval |y|<0.8 The first (sys) error is the systematic uncertainty from the other sources The second (sys) error is the systematic uncertainty from the B feed-down contribution.

Ratio of v_{2}{EP} in the 60% small-q_{2}^{V0A} sample to the unbiased sample vs. p_{T} of average prompt D0, D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the rapidity interval |y|<0.8.

Ratio of v_{2}{EP} in the 20% large-q_{2}^{V0A} sample to the unbiased sample vs. p_{T} of average prompt D0, D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the rapidity interval |y|<0.8.

Ratio of v_{2}{EP} in the 60% small-q_{2}^{V0A} sample to the unbiased sample vs. p_{T} of average prompt D0, D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the rapidity interval |y|<0.8.

Ratio of v_{2}{EP} in the 20% large-q_{2}^{V0A} sample to the unbiased sample vs. p_{T} of average prompt D0, D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the rapidity interval |y|<0.8.

v_{2}{SP} vs. p_{T} of charged particles in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the rapidity interval |eta|<0.8.

v_{2}{SP} vs. p_{T} of charged particles in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the 60% small-q_{2}^{V0A} sample in the rapidity interval |eta|<0.8.

v_{2}{SP} vs. p_{T} of charged particles in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the 20% large-q_{2}^{V0A} sample in the rapidity interval |eta|<0.8.

v_{2}{SP} vs. p_{T} of charged particles in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the rapidity interval |eta|<0.8.

v_{2}{SP} vs. p_{T} of charged particles in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the 60% small-q_{2}^{V0A} sample in the rapidity interval |eta|<0.8.

v_{2}{SP} vs. p_{T} of charged particles in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the 20% large-q_{2}^{V0A} sample in the rapidity interval |eta|<0.8.

Ratio of v_{2}{SP} in the 60% small-q_{2}^{V0A} sample to the unbiased sample vs. p_{T} of charged particles in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the rapidity interval |eta|<0.8.

Ratio of v_{2}{SP} in the 20% large-q_{2}^{V0A} sample to the unbiased sample vs. p_{T} of charged particles in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the rapidity interval |eta|<0.8.

Ratio of v_{2}{SP} in the 60% small-q_{2}^{V0A} sample to the unbiased sample vs. p_{T} of charged particles in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the rapidity interval |eta|<0.8.

Ratio of v_{2}{SP} in the 20% large-q_{2}^{V0A} sample to the unbiased sample vs. p_{T} of charged particles in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the rapidity interval |eta|<0.8.

Ratio of the p_{T} distribution in the 60% small-q_{2}^{TPC} sample to the unbiased sample of prompt D0 mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the rapidity interval |y|<0.8.

Ratio of the p_{T} distribution in the 60% small-q_{2}^{TPC} sample to the unbiased sample of prompt D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the rapidity interval |y|<0.8.

Ratio of the p_{T} distribution in the 20% large-q_{2}^{TPC} sample to the unbiased sample of prompt D0 mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the rapidity interval |y|<0.8.

Ratio of the p_{T} distribution in the 20% large-q_{2}^{TPC} sample to the unbiased sample of prompt D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the rapidity interval |y|<0.8.

Ratio of the p_{T} distribution in the 60% small-q_{2}^{TPC} sample to the unbiased sample of prompt D0 mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the rapidity interval |y|<0.8.

Ratio of the p_{T} distribution in the 60% small-q_{2}^{TPC} sample to the unbiased sample of prompt D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the rapidity interval |y|<0.8.

Ratio of the p_{T} distribution in the 20% large-q_{2}^{TPC} sample to the unbiased sample of prompt D0 mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the rapidity interval |y|<0.8.

Ratio of the p_{T} distribution in the 20% large-q_{2}^{TPC} sample to the unbiased sample of prompt D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the rapidity interval |y|<0.8.

Ratio of the p_{T} distribution in the 60% small-q_{2}^{TPC} sample to the unbiased sample of average prompt D0, D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the rapidity interval |y|<0.8.

Ratio of the p_{T} distribution in the 20% large-q_{2}^{TPC} sample to the unbiased sample of average prompt D0, D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 10-30% in the rapidity interval |y|<0.8.

Ratio of the p_{T} distribution in the 60% small-q_{2}^{TPC} sample to the unbiased sample of average prompt D0, D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the rapidity interval |y|<0.8.

Ratio of the p_{T} distribution in the 20% large-q_{2}^{TPC} sample to the unbiased sample of average prompt D0, D+ mesons in Pb-Pb collisions at sqrt{s_{NN}}=5.02 TeV in the centrality class 30-50% in the rapidity interval |y|<0.8.

Invariant mass distribution of the selected $\mathrm{B}^+\mathrm{K}^-$ candidates in data

Invariant mass distribution of the selected $\mathrm{B}^0 \mathrm{K}^{0}_{\mathrm{S}}$ candidates in data

The observed signal yields ($N$), natural widths ($\Gamma$), and mass differences from the fits to the $m_{\mathrm{BK}}$ distributions in data, for the two channels: $\mathrm{B}^+\mathrm{K}^-$ and $\mathrm{B}^0 \mathrm{K}^{0}_{\mathrm{S}}$. The uncertainties are statistical only.

The measured ratios of branching fractions

The measured ratios of cross sections times branching fractions

The measured mass differences $M(\mathrm{B}^{(*)0}_{\mathrm{s1,2}})-M^{\mathrm{PDG}}_{\mathrm{B}}-M^{\mathrm{PDG}}_{\mathrm{K}}$

The measured masses of $\mathrm{B}^{*}_{\mathrm{s2}}(5840)^0$ and $\mathrm{B}_{\mathrm{s1}}(5830)^0$ mesons

The measured mass differences between neutral and charged B mesons

The measured $\mathrm{B}^{*}_{\mathrm{s2}}(5840)^0$ natural width

Acceptance x efficiency versus resonance mass for both lephad and hadhad channels in the RS bulk model with k/MPl = 1

Acceptance x efficiency versus resonance mass for both lephad and hadhad channels in the RS bulk model with k/MPl = 2

Acceptance x efficiency versus resonance mass for both lephad and hadhad channels in the scalar model

Upper limits on the production cross-section times the HH to bbtautau branching ratio for non-resonant HH at 95% CLS and their interpretation as multiples of the SM prediction

Upper limits on the production cross-section times the HH to bbtautau branching ratio divided by the SM prediction for non-resonant HH at 95% CL

Post-fit expected number of signal and background events and observed number of data events after applying the selection criteria and requiring exactly 2 b-tagged jets and assuming a background-only hypothesis

Post-fit expected number of signal and background events and observed number of data events in the last two bins of the non-resonant BDT score distribution of the SM signal after applying the selection criteria and requiring exactly 2 b-tagged jets and assuming a background-only hypothesis

Transverse momentum spectra of charged particles in V0M IV multiplicity class

Transverse momentum spectra of charged particles in V0M V multiplicity class

Transverse momentum spectra of charged particles in V0M VI multiplicity class

Transverse momentum spectra of charged particles in V0M VII multiplicity class

Transverse momentum spectra of charged particles in V0M VIII multiplicity class

Transverse momentum spectra of charged particles in V0M IX multiplicity class

Transverse momentum spectra of charged particles in V0M X multiplicity class

Transverse momentum spectra of $\pi^{+} + \pi^{-}$ in V0M multiplicity class I

Transverse momentum spectra of $\pi^{+} + \pi^{-}$ in V0M multiplicity class II

Transverse momentum spectra of $\pi^{+} + \pi^{-}$ in V0M multiplicity class III

Transverse momentum spectra of $\pi^{+} + \pi^{-}$ in V0M multiplicity class IV

Transverse momentum spectra of $\pi^{+} + \pi^{-}$ in V0M multiplicity class V

Transverse momentum spectra of $\pi^{+} + \pi^{-}$ in V0M multiplicity class VI

Transverse momentum spectra of $\pi^{+} + \pi^{-}$ in V0M multiplicity class VII

Transverse momentum spectra of $\pi^{+} + \pi^{-}$ in V0M multiplicity class VIII

Transverse momentum spectra of $\pi^{+} + \pi^{-}$ in V0M multiplicity class IX

Transverse momentum spectra of $\pi^{+} + \pi^{-}$ in V0M multiplicity class X

Transverse momentum spectra of $K^{+} + K^{-}$ in V0M multiplicity class I

Transverse momentum spectra of $K^{+} + K^{-}$ in V0M multiplicity class II

Transverse momentum spectra of $K^{+} + K^{-}$ in V0M multiplicity class III

Transverse momentum spectra of $K^{+} + K^{-}$ in V0M multiplicity class IV

Transverse momentum spectra of $K^{+} + K^{-}$ in V0M multiplicity class V

Transverse momentum spectra of $K^{+} + K^{-}$ in V0M multiplicity class VI

Transverse momentum spectra of $K^{+} + K^{-}$ in V0M multiplicity class VII

Transverse momentum spectra of $K^{+} + K^{-}$ in V0M multiplicity class VIII

Transverse momentum spectra of $K^{+} + K^{-}$ in V0M multiplicity class IX

Transverse momentum spectra of $K^{+} + K^{-}$ in V0M multiplicity class X

Transverse momentum spectra of $p + \bar{p}$ in V0M multiplicity class I

Transverse momentum spectra of $p + \bar{p}$ in V0M multiplicity class II

Transverse momentum spectra of $p + \bar{p}$ in V0M multiplicity class III

Transverse momentum spectra of $p + \bar{p}$ in V0M multiplicity class IV

Transverse momentum spectra of $p + \bar{p}$ in V0M multiplicity class V

Transverse momentum spectra of $p + \bar{p}$ in V0M multiplicity class VI

Transverse momentum spectra of $p + \bar{p}$ in V0M multiplicity class VII

Transverse momentum spectra of $p + \bar{p}$ in V0M multiplicity class VIII

Transverse momentum spectra of $p + \bar{p}$ in V0M multiplicity class IX

Transverse momentum spectra of $p + \bar{p}$ in V0M multiplicity class X

Transverse momentum spectra of $K* + \bar{K}*$ in V0M multiplicity class I

Transverse momentum spectra of $K* + \bar{K}*$ in V0M multiplicity class II

Transverse momentum spectra of $K* + \bar{K}*$ in V0M multiplicity class III

Transverse momentum spectra of $K* + \bar{K}*$ in V0M multiplicity class IV + V

Transverse momentum spectra of $K* + \bar{K}*$ in V0M multiplicity class VI

Transverse momentum spectra of $K* + \bar{K}*$ in V0M multiplicity class VII

Transverse momentum spectra of $K* + \bar{K}*$ in V0M multiplicity class VIII

Transverse momentum spectra of $K* + \bar{K}*$ in V0M multiplicity class IX

Transverse momentum spectra of $K* + \bar{K}*$ in V0M multiplicity class X

Transverse momentum spectra of $\phi$ in V0M multiplicity class I

Transverse momentum spectra of $\phi$ in V0M multiplicity class II

Transverse momentum spectra of $\phi$ in V0M multiplicity class III

Transverse momentum spectra of $\phi$ in V0M multiplicity class IV + V

Transverse momentum spectra of $\phi$ in V0M multiplicity class VI

Transverse momentum spectra of $\phi$ in V0M multiplicity class VII

Transverse momentum spectra of $\phi$ in V0M multiplicity class VIII

Transverse momentum spectra of $\phi$ in V0M multiplicity class IX

Transverse momentum spectra of $\phi$ in V0M multiplicity class X

Unidentified particle spectral modification in V0M Class IX pp collisions.

$\pi^{+}+\pi^{-}$ spectral modification in V0M Class IX pp collisions.

$K^{+}+K^{-}$ spectral modification in V0M Class IX pp collisions.

$p+\bar{p}$ spectral modification in V0M Class IX pp collisions.

$\Lambda + \bar{\Lambda}$ spectral modification in V0M Class IX pp collisions.

$\Xi^{-} + \bar{\Xi}^{+}$ spectral modification in V0M Class IX pp collisions.

$K* + \bar{K}*$ spectral modification in V0M Class IX pp collisions.

$\phi$ spectral modification in V0M Class IX pp collisions.

Unidentified particle spectral modification in V0M Class I pp collisions.

$\pi^{+}+\pi^{-}$ spectral modification in V0M Class I pp collisions.

$K^{+}+K^{-}$ spectral modification in V0M Class I pp collisions.

$p+\bar{p}$ spectral modification in V0M Class I pp collisions.

$\Lambda + \bar{\Lambda}$ spectral modification in V0M Class I pp collisions.

$\Xi^{-} + \bar{\Xi}^{+}$ spectral modification in V0M Class I pp collisions.

$K* + \bar{K}*$ spectral modification in V0M Class I pp collisions.

$\phi$ spectral modification in V0M Class I pp collisions.

Transverse momentum dependence of $(p + \bar{p})/\phi$ yield ratios in low multiplicity pp collisions at $\sqrt{s} = 7$ (TeV).

Transverse momentum dependence of $(K^{+} + K^{-})/(\pi^{+} + \pi^{-})$ yield ratios in low multiplicity pp collisions at $\sqrt{s} = 7$ (TeV).

Transverse momentum dependence of $(p + \bar{p})/(\pi^{+} + \pi^{-})$ yield ratios in low multiplicity pp collisions at $\sqrt{s} = 7$ (TeV).

Transverse momentum dependence of $\Lambda/K^{0}_{s}$ yield ratios in low multiplicity pp collisions at $\sqrt{s} = 7$ (TeV).

Transverse momentum dependence of $(p + \bar{p})/\phi$ yield ratios in high multiplicity pp collisions at $\sqrt{s} = 7$ (TeV).

Transverse momentum dependence of $(K^{+} + K^{-})/(\pi^{+} + \pi^{-})$ yield ratios in high multiplicity pp collisions at $\sqrt{s} = 7$ (TeV).

Transverse momentum dependence of $(p + \bar{p})/(\pi^{+} + \pi^{-})$ yield ratios in high multiplicity pp collisions at $\sqrt{s} = 7$ (TeV).

Transverse momentum dependence of $\Lambda/K^{0}_{s}$ yield ratios in high multiplicity pp collisions at $\sqrt{s} = 7$ (TeV).

Transverse momentum dependence of $(p + \bar{p})/\phi$ ratio in low-$p_{T}$ interval in pp collisions as a function of $\langle dN_{ch}/d\eta \rangle_{|\eta|<0.5}$.

Transverse momentum dependence of $(K^{+} + K^{-})/(\pi^{+} + \pi^{-})$ ratio in low-$p_{T}$ interval in pp collisions as a function of $\langle dN_{ch}/d\eta \rangle_{|\eta|<0.5}$.

Transverse momentum dependence of $(p + \bar{p})/(\pi^{+} + \pi^{-})$ ratio in low-$p_{T}$ interval in pp collisions as a function of $\langle dN_{ch}/d\eta \rangle_{|\eta|<0.5}$.

Transverse momentum dependence of $\Lambda/K^{0}_{s}$ ratio in low-$p_{T}$ interval in pp collisions as a function of $\langle dN_{ch}/d\eta \rangle_{|\eta|<0.5}$.

Transverse momentum dependence of $(p + \bar{p})/\phi$ ratio in mid-$p_{T}$ interval in pp collisions as a function of $\langle dN_{ch}/d\eta \rangle_{|\eta|<0.5}$.

Transverse momentum dependence of $(K^{+} + K^{-})/(\pi^{+} + \pi^{-})$ ratio in mid-$p_{T}$ interval in pp collisions as a function of $\langle dN_{ch}/d\eta \rangle_{|\eta|<0.5}$.

Transverse momentum dependence of $(p + \bar{p})/(\pi^{+} + \pi^{-})$ ratio in mid-$p_{T}$ interval in pp collisions as a function of $\langle dN_{ch}/d\eta \rangle_{|\eta|<0.5}$.

Transverse momentum dependence of $\Lambda/K^{0}_{s}$ ratio in mid-$p_{T}$ interval in pp collisions as a function of $\langle dN_{ch}/d\eta \rangle_{|\eta|<0.5}$.

Transverse momentum dependence of $(p + \bar{p})/\phi$ ratio in high-$p_{T}$ interval in pp collisions as a function of $\langle dN_{ch}/d\eta \rangle_{|\eta|<0.5}$.

Transverse momentum dependence of $(K^{+} + K^{-})/(\pi^{+} + \pi^{-})$ ratio in high-$p_{T}$ interval in pp collisions as a function of $\langle dN_{ch}/d\eta \rangle_{|\eta|<0.5}$.

Transverse momentum dependence of $(p + \bar{p})/(\pi^{+} + \pi^{-})$ ratio in high-$p_{T}$ interval in pp collisions as a function of $\langle dN_{ch}/d\eta \rangle_{|\eta|<0.5}$.

Transverse momentum dependence of $\Lambda/K^{0}_{s}$ ratio in high-$p_{T}$ interval in pp collisions as a function of $\langle dN_{ch}/d\eta \rangle_{|\eta|<0.5}$.

Integrated $(K^{+} + K^{-})$/$(\pi^{+} + \pi^{-})$ ratios as a function of $\langle dN_{ch}/d\eta\rangle$ in pp collisions at $\sqrt{s} = 7$ TeV.

Integrated $(K^{*} + \bar{K}^{*})$/$(\pi^{+} + \pi^{-})$ ratios as a function of $\langle dN_{ch}/d\eta\rangle$ in pp collisions at $\sqrt{s} = 7$ TeV.

Integrated $\phi$/$(\pi^{+} + \pi^{-})$ ratios as a function of $\langle dN_{ch}/d\eta\rangle$ in pp collisions at $\sqrt{s} = 7$ TeV.

$\langle\beta_{T}\rangle$ from blast-wave fit to PbPb collisions at $\sqrt{s}$ = 2.76 TeV

$T_{kin}$ from blast-wave fit to PbPb collisions at $\sqrt{s}$ = 2.76 TeV

$n$ from blast-wave fit to PbPb collisions at $\sqrt{s}$ = 2.76 TeV

$\langle\beta_{T}\rangle$ from blast-wave fit to pPb collisions at $\sqrt{s}$ = 5.02 TeV

$T_{kin}$ from blast-wave fit to pPb collisions at $\sqrt{s}$ = 5.02 TeV

$n$ from blast-wave fit to pPb collisions at $\sqrt{s}$ = 5.02 TeV

$\langle\beta_{T}\rangle$ from blast-wave fit to pp collisions at $\sqrt{s}$ = 7 TeV

$T_{kin}$ from blast-wave fit to pp collisions at $\sqrt{s}$ = 7 TeV

$n$ from blast-wave fit to pp collisions at $\sqrt{s}$ = 7 TeV

Dijet mass spectra after the background only fit with the background prediction in the high-mass region with two b-tags.

Dijet mass spectra after the background only fit with the background prediction in the low-mass region with two b-tags.

The online b-tagging efficiency with respect to the offline b-tagging efficiency as a function of pT. The b-tagging online and offline working points correspond to an efficiency of 60% and 70%, respectively.

Observed and expected 95% credibility-level upper limits on the cross-section for the b* model in the high-mass region with inclusive b-jet selection.

Observed and expected 95% credibility-level upper limits on the cross-section times branching ratio for the SSM and leptophobic Z' models in the low- and high-mass region with two b-tags selection.

Observed and expected 95% credibility-level upper limits on the cross-section for DM Z' models in the low-mass region with two b-tags selection. The Z' is expected to decay to all five quark flavors other than the top quark and the mediator to SM quark coupling (gSM) equal to 0.1 is assumed.

Observed and expected 95% credibility-level upper limits on the cross-section times branching ratio for DM Z'->bb models in the high-mass region with two b-tags selection. The Z' is expected to decay to bb only and the mediator to SM quark coupling (gSM) equal to 0.25 is assumed.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances with intrinsic width smaller than the detector resolution. The width of the Gaussian reconstructed shape is dominated by the dijet mass resolution. The table shows the limits obtained from the high-mass inclusive one b-tag selection.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances with intrinsic width smaller than the detector resolution. The width of the Gaussian reconstructed shape is dominated by the dijet mass resolution. The table shows the limits obtained from the combined low- and high-mass two b-tags selection.

The mass distributions for the inclusive one b-tagged selection and two b-tagged selection using an integrated luminosity of 36.1 $fb^{-1}$. The inclusive one b-tagged Pythia8 MC distribution is normalized to the inclusive one b-tagged data. The two b-tagged Pythia8 MC distribution is normalized to the two b-tagged data. The systematic uncertainty band is dominated by the b-tagging scale factor and the b-jet energy scale.

Signal acceptance times efficiency for inclusive 1 b-tag and 2 b-tag categories as a function of the simulated b* and the Z' masses.

Signal acceptance times efficiency for 2 b-tag categories as a function of the simulated Z' masses.

The flavor composition of the simulated dijet background as a function of dijet mass before tagging. The fraction of the six combinations of the b-jet , c-jet and light-flavor jet are shown. All offline selections are applied.

The flavor composition of the simulated dijet background as a function of dijet mass with inclusive one b-tag. The fraction of the six combinations of the b-jet , c-jet and light-flavor jet are shown. All offline selections are applied.

The flavor composition of the simulated dijet background as a function of dijet mass with two b-tags. The fraction of the six combinations of the b-jet , c-jet and light-flavor jet are shown. All offline selections are applied.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the inclusive b-jet selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 3%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the inclusive b-jet selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 7%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the inclusive b-jet selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 10%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the inclusive b-jet selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 15%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the combined low- and high-mass two b-tags selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 3%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the combined low- and high-mass two b-tags selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 7%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the combined low- and high-mass two b-tags selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 10%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the combined low- and high-mass two b-tags selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 15%.

d$\sigma$/dy for non-prompt J/$\psi$ in p-Pb collisions at $\sqrt{s_NN}=5.02 TeV$

d$^2\sigma$/d$y$d$p_{\rm T}$ in bins of $p_T^{J/\psi}$ for non-prompt J/$\psi$ in p-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. The data point in the first $p_{\rm T}$ bin is un upper limit at 95% Confidence Level.

Nuclear modification factor ($R_{pPb}$) in bins of $p_T^{J/\psi}$ for non-prompt J/$\psi$ in p-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. The data point in the first $p_{\rm T}$ bin is un upper limit at 95% Confidence Level.

Nuclear modification factor ($R_{pPb}$) for non-prompt J/$\psi$ integrated over $p_{\rm T}$ in p-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV at mid-rapidity.

d$^2\sigma$/d $p_{\rm{T}}$dy versus $p_T^{J/\psi}$ for prompt J/$\psi$ in p-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

Nuclear modification factor ($R_{pPb}$) in bins of $p_T^{J/\psi}$ for prompt J/$\psi$ in p-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

Nuclear modification factor ($R_{pPb}$) for prompt J/$\psi$ integrated over $p_{\rm T}$ in p-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV at mid-rapidity.

Beauty production cross section at mid-rapidity (d$sigma$/d$y$) in p-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

d$\sigma$/dy for prompt J/$\psi$ in p-Pb collisions at $\sqrt{s_NN}=5.02 TeV$ integrated over $p_{\rm T}$>0.

Ratio of the measured ETA to PI0 ratio to the one obtained via m_T scaling in NSD p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV.

The nuclear modification factor RPPB of PI0 as a function of pT in NSD p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV, the overall normalization uncertainty of 5.96 % is not included.

The nuclear modification factor RPPB of ETA as a function of pT in NSD p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV, the overall normalization uncertainty of 5.96 % is not included.

Measured B$^\pm$ differential cross-sections (in units of $\mu$b) at 7 TeV and 13 TeV as functions of $y$ in the $p_T$ range $0<p_T<40$GeV$/c$. The cross-section ratio between 13 TeV and 7 TeV is also presented.

Invariant differential yields of tritons and antitritons in inelastic pp collisions at $\sqrt{s}$ = 7 TeV. The uncertainties of $_{-2.0}^{+5.0}$% due to the extrapolation to inelastic pp collisions are not included in the systematic uncertainties.

Invariant differential yields of $^{3}$He and $^{3}\overline{\text{He}}$ nuclei in inelastic pp collisions at $\sqrt{s}$ = 7 TeV. The uncertainties of $_{-2.0}^{+5.0}$% due to the extrapolation to inelastic pp collisions are not included in the systematic uncertainties.

Coalescence parameter ($B_{2}$) of deuterons and antideuterons in inelastic pp collisions at $\sqrt{s}$ = 0.9 TeV.

Coalescence parameter ($B_{2}$) of deuterons and antideuterons in inelastic pp collisions at $\sqrt{s}$ = 2.76 TeV.

Coalescence parameter ($B_{2}$) of deuterons and antideuterons in inelastic pp collisions at $\sqrt{s}$ = 7 TeV.

Coalescence parameter ($B_{3}$) of tritons and antitritons in inelastic pp collisions at $\sqrt{s}$ = 7 TeV.

Coalescence parameter ($B_{3}$) of $^{3}$He and $^{3}\overline{\text{He}}$ nuclei in inelastic pp collisions at $\sqrt{s}$ = 7 TeV.

Integrated deuteron-to-proton (d/p) and antideuteron-to-antiproton ($\overline{\text{d}}/\overline{\text{p}}$) ratios in inelastic pp collisions as a function of the average charged particle multiplicity for different center-of-mass energies.

RFB of inclusive JPSI(1S) in p-Pb collisions as a function of relative charged-particle pseudorapidity density.

Relative mean transverse momentum of JPSI(1S) mesons for backward rapidity as a function of the relative charged-particle pseudorapidity density.

Relative mean transverse momentum of JPSI(1S) mesons for forward rapidity as a function of the relative charged-particle pseudorapidity density.

Definition of the fiducial phase space.

The seven input variables to the NN ordered by their discriminating power. The jet that is not $b$-tagged is referred to as $\textit{untagged}~$jet.

The seven input variables to the NN ordered by their discriminating power. The jet that is not $b$-tagged is referred to as $\textit{untagged}~$jet.

Event yields for the different processes estimated with the fit to the $O_\mathrm{NN}$ distribution compared to the numbers of observed events. Only the statistical uncertainties are quoted. The $Z,VV+\mathrm{jets}$ contributions and the multijet background are fixed in the fit; therefore no uncertainty is quoted for these processes.

Event yields for the different processes estimated with the fit to the $O_\mathrm{NN}$ distribution compared to the numbers of observed events. Only the statistical uncertainties are quoted. The $Z,VV+\mathrm{jets}$ contributions and the multijet background are fixed in the fit; therefore no uncertainty is quoted for these processes.

Detailed list of the contribution from each source of uncertainty to the total uncertainty in the measured values of $\sigma_{\mathrm{fid}}(tq)$ and $\sigma_{\mathrm{fid}}(\bar tq)$. The estimation of the systematic uncertainties has a statistical uncertainty of $0.3\%$. Uncertainties contributing less than $0.5\%$ are marked with ‘<0.5’.

Detailed list of the contribution from each source of uncertainty to the total uncertainty in the measured values of $\sigma_{\mathrm{fid}}(tq)$ and $\sigma_{\mathrm{fid}}(\bar tq)$. The estimation of the systematic uncertainties has a statistical uncertainty of $0.3\%$. Uncertainties contributing less than $0.5\%$ are marked with ‘<0.5’.

Significant contributions to the total relative uncertainty in the measured value of $R_{t}$. The estimation of the systematic uncertainties has a statistical uncertainty of $0.3~\%$. Uncertainties contributing less than $0.5~\%$ are not shown.

Significant contributions to the total relative uncertainty in the measured value of $R_{t}$. The estimation of the systematic uncertainties has a statistical uncertainty of $0.3~\%$. Uncertainties contributing less than $0.5~\%$ are not shown.

Slopes $a$ of the mass dependence of the measured cross$-$sections.

Slopes $a$ of the mass dependence of the measured cross$-$sections.

Predicted (post-fit) and observed event yields for the signal region (SR), after the requirement on the neural network discriminant, $O_{\mathrm{NN}}~>~0.8$. The multijet background prediction is obtained from the fit to the $E_{\mathrm{T}}^{\mathrm{miss}}$ distribution described in Section 6, while all the other predictions and uncertainties are derived from the total cross$-$section measurement. In some cases there is no uncertainty quoted. In these cases the uncertainty is < 0.5.

Predicted (post-fit) and observed event yields for the signal region (SR), after the requirement on the neural network discriminant, $O_{\mathrm{NN}}~>~0.8$. The multijet background prediction is obtained from the fit to the $E_{\mathrm{T}}^{\mathrm{miss}}$ distribution described in Section 6, while all the other predictions and uncertainties are derived from the total cross$-$section measurement. In some cases there is no uncertainty quoted. In these cases the uncertainty is < 0.5.

Predicted (post-fit) and observed event yields for the signal region (SR), after the requirement on the second neural network discriminant, $O_{\mathrm{NN2}}~>~0.8$. The multijet background prediction is obtained from the fit to the $E_{\mathrm{T}}^{\mathrm{miss}}$ distribution described in Section 6, while all the other predictions and uncertainties are derived from the total cross$-$section measurement. In some cases there is no uncertainty quoted. In these cases the uncertainty is < 0.5.

Predicted (post-fit) and observed event yields for the signal region (SR), after the requirement on the second neural network discriminant, $O_{\mathrm{NN2}}~>~0.8$. The multijet background prediction is obtained from the fit to the $E_{\mathrm{T}}^{\mathrm{miss}}$ distribution described in Section 6, while all the other predictions and uncertainties are derived from the total cross$-$section measurement. In some cases there is no uncertainty quoted. In these cases the uncertainty is < 0.5.

Migration matrix for $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at the particle level. The pseudo top quark is shown on the $y$-axis and the reconstructed variable is shown on the $x$-axis.

Migration matrix for $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at the particle level. The pseudo top quark is shown on the $y$-axis and the reconstructed variable is shown on the $x$-axis.

Migration matrix for $p_{\mathrm{T}}(t)$ at the parton level. The parton-level quark is shown on the $y$-axis and the reconstructed variable is shown on the $x$-axis.

Migration matrix for $p_{\mathrm{T}}(t)$ at the parton level. The parton-level quark is shown on the $y$-axis and the reconstructed variable is shown on the $x$-axis.

Migration matrix for $|y(\hat{t\hspace{-0.2mm}})|$ at the particle level. The pseudo top quark is shown on the $y$-axis and the reconstructed variable is shown on the $x$-axis.

Migration matrix for $|y(\hat{t\hspace{-0.2mm}})|$ at the particle level. The pseudo top quark is shown on the $y$-axis and the reconstructed variable is shown on the $x$-axis.

Migration matrix for $|y(t)|$ at the parton level. The parton-level quark is shown on the $y$-axis and the reconstructed variable is shown on the $x$-axis.

Migration matrix for $|y(t)|$ at the parton level. The parton-level quark is shown on the $y$-axis and the reconstructed variable is shown on the $x$-axis.

Uncertainties in the normalisations of the different backgrounds for all processes, as derived from the total cross-section measurement.

Uncertainties in the normalisations of the different backgrounds for all processes, as derived from the total cross-section measurement.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at particle level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at particle level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at particle level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at particle level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ at particle level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ at particle level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ at particle level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ at particle level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ at particle level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ at particle level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ at particle level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ at particle level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ at particle level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ at particle level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ at particle level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ at particle level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $p_{\mathrm{T}}(t)$ at parton level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $p_{\mathrm{T}}(t)$ at parton level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $p_{\mathrm{T}}(t)$ at parton level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $p_{\mathrm{T}}(t)$ at parton level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $|y(t)|$ at parton level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $|y(t)|$ at parton level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $|y(t)|$ at parton level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $|y(t)|$ at parton level.

Statistical correlation matrix for the absolute differential cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ for $tq$ events(at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ for $ \bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ for $\bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ for $ \bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ for $\bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ for $\bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ for $\bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ for $\bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ for $\bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $p_{\mathrm{T}}(t)$ for $tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $p_{\mathrm{T}}(t)$ for $ \bar tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $p_{\mathrm{T}}(t)$ for $tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $p_{\mathrm{T}}(t)$ for $ \bar tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $|y(t)|$ for $tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $|y(t)|$ for $\bar tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $|y(t)|$ for $tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $|y(t)|$ for $\bar tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Fiducial acceptance $A_{\mathrm{fid}}$ for different $t$-channel single top-quark MC samples. $^{\mathrm{(a)}}$ Calculation taken from AcerMC $+$ $\mathrm{P{\scriptsize YTHIA}6}$. $^{\mathrm{(b)}}$ Calculation taken from $\mathrm{P{\scriptsize OWHEG}}$-$\mathrm{B{\scriptsize OX}}$ $+$ $\mathrm{P{\scriptsize YTHIA}6}$.

Fiducial acceptance $A_{\mathrm{fid}}$ for different $t$-channel single top-antiquark MC samples. $^{\mathrm{(a)}}$ Calculation taken from AcerMC $+$ $\mathrm{P{\scriptsize YTHIA}6}$. $^{\mathrm{(b)}}$ Calculation taken from $\mathrm{P{\scriptsize OWHEG}}$-$\mathrm{B{\scriptsize OX}}$ $+$ $\mathrm{P{\scriptsize YTHIA}6}$.

Uncertainties for the absolute differential $tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at particle level per bin ([0,35,50,75,100,150,200,300] GeV) in percent of $\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at particle level per bin ([0,35,50,75,100,150,200,300] GeV) in percent of $\left( \dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $\bar tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at particle level per bin ([0,35,50,75,100,150,300] GeV) in percent of $\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})}$.

Uncertainties for the normalised differential $\bar tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at particle level per bin ([0,35,50,75,100,150,300] GeV) in percent of $\left( \dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})}$.

Uncertainties for the absolute differential $tq$ cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ at particle level per bin ([0,0.15,0.3,0.45,0.7,1.0,1.3,2.2]) in percent of $\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(\hat{t\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $tq$ cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ at particle level per bin ([0,0.15,0.3,0.45,0.7,1.0,1.3,2.2]) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(\hat{t\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $\bar tq$ cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ at particle level per bin ([0,0.15,0.3,0.45,0.7,1.0,1.3,2.2]) in percent of $\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(\hat{t\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $\bar tq$ cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ at particle level per bin ([0,0.15,0.3,0.45,0.7,1.0,1.3,2.2]) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(\hat{t\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ at particle level per bin ([30,45,60,75,100,150,300] GeV) in percent of $\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ at particle level per bin ([30,45,60,75,100,150,300] GeV) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $\bar tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ at particle level per bin ([30,45,60,75,100,150,300] GeV) in percent of $\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $\bar tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ at particle level per bin ([30,45,60,75,100,150,300] GeV) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $tq$ cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ at particle level per bin ([0.0, 1.2, 1.7, 2.2, 2.7, 3.3, 4.5]) in percent of $\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(\hat{j\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $tq$ cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ at particle level per bin ([0.0, 1.2, 1.7, 2.2, 2.7, 3.3, 4.5]) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(\hat{j\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $\bar tq$ cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ at particle level per bin ([0.0, 1.2, 1.7, 2.2, 2.7, 3.3, 4.5]) in percent of $\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(\hat{j\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $\bar tq$ cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ at particle level per bin ([0.0, 1.2, 1.7, 2.2, 2.7, 3.3, 4.5]) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(\hat{j\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $tq$ cross-section as a function of $p_{\mathrm{T}}(t)$ at parton level per bin ([0,50,100,150,200,300] GeV) in percent of $\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}p_{\mathrm{T}}(t)}$.

Uncertainties for the normalised differential $tq$ cross-section as a function of $p_{\mathrm{T}}(t)$ at parton level per bin ([0,50,100,150,200,300] GeV) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}p_{\mathrm{T}}(t)}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $\bar tq $ cross-section as a function of $p_{\mathrm{T}}(t)$ at parton level per bin ([0,50,100,150,300] GeV) in percent of $\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}p_{\mathrm{T}}(t)}$.

Uncertainties for the normalised differential $\bar tq $ cross-section as a function of $p_{\mathrm{T}}(t)$ at parton level per bin ([0,50,100,150,300] GeV) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}p_{\mathrm{T}}(t)}$.

Uncertainties for the absolute differential $ tq $ cross-section as a function of $|y(t)|$ at parton level per bin ([0,0.3,0.7,1.3,2.2]) in percent of $\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(t)|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $ tq $ cross-section as a function of $|y(t)|$ at parton level per bin ([0,0.3,0.7,1.3,2.2]) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(t)|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $ \bar tq $ cross-section as a function of $|y(t)|$ at parton level per bin ([0,0.3,0.7,1.3,2.2]) in percent of $\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(t)|}$.

Uncertainties for the normalised differential $ \bar tq $ cross-section as a function of $|y(t)|$ at parton level per bin ([0,0.3,0.7,1.3,2.2]) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(t)|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

$p_{\rm{T}}$-differential yield of $\Sigma^{*+}$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (NSD).

$p_{\rm{T}}$-differential yield of $\Sigma^{*+}$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (0-20% multiplicity class).

$p_{\rm{T}}$-differential yield of $\Sigma^{*+}$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (0-20% multiplicity class).

$p_{\rm{T}}$-differential yield of $\Sigma^{*+}$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (20-60% multiplicity class).

$p_{\rm{T}}$-differential yield of $\Sigma^{*+}$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (20-60% multiplicity class).

$p_{\rm{T}}$-differential yield of $\Sigma^{*+}$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (60-100% multiplicity class).

$p_{\rm{T}}$-differential yield of $\Sigma^{*+}$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (60-100% multiplicity class).

$p_{\rm{T}}$-differential yield of ($\Xi^{*0}+\overline{\Xi}^{*0}$)/2 in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (NSD).

$p_{\rm{T}}$-differential yield of ($\Xi^{*0}+\overline{\Xi}^{*0}$)/2 in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (NSD).

$p_{\rm{T}}$-differential yield of ($\Xi^{*0}+\overline{\Xi}^{*0}$)/2 in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (0-20\% multiplicity class).

$p_{\rm{T}}$-differential yield of ($\Xi^{*0}+\overline{\Xi}^{*0}$)/2 in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (0-20\% multiplicity class).

$p_{\rm{T}}$-differential yield of ($\Xi^{*0}+\overline{\Xi}^{*0}$)/2 in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (20-40\% multiplicity class).

$p_{\rm{T}}$-differential yield of ($\Xi^{*0}+\overline{\Xi}^{*0}$)/2 in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (20-40\% multiplicity class).

$p_{\rm{T}}$-differential yield of ($\Xi^{*0}+\overline{\Xi}^{*0}$)/2 in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (40-60\% multiplicity class).

$p_{\rm{T}}$-differential yield of ($\Xi^{*0}+\overline{\Xi}^{*0}$)/2 in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (40-60\% multiplicity class).

$p_{\rm{T}}$-differential yield of ($\Xi^{*0}+\overline{\Xi}^{*0}$)/2 in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (60-100\% multiplicity class).

$p_{\rm{T}}$-differential yield of ($\Xi^{*0}+\overline{\Xi}^{*0}$)/2 in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (60-100\% multiplicity class).

$p_{\rm T}$-integrated yield d$N$/d$y$ of ($\Sigma^{*\pm}+\overline{\Sigma}^{*\mp}$)/4 in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV. Yields are normalised to the visible cross section for each V0A multiplicity event class and to NSD events for the minimum bias case.

$p_{\rm T}$-integrated yield d$N$/d$y$ of ($\Sigma^{*\pm}+\overline{\Sigma}^{*\mp}$)/4 in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV. Yields are normalised to the visible cross section for each V0A multiplicity event class and to NSD events for the minimum bias case.

Mean $p_{\rm T}$ of $\Sigma^{*\pm}$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV for each V0A multiplicity event class.

Mean $p_{\rm T}$ of $\Sigma^{*\pm}$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV for each V0A multiplicity event class.

$p_{\rm T}$-integrated yield d$N$/d$y$ of ($\Xi^{*0}+\overline{\Xi}^{0})$)/2 in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV. Yields are normalised to the visible cross section for each V0A multiplicity event class and to NSD events for the minimum bias case.

$p_{\rm T}$-integrated yield d$N$/d$y$ of ($\Xi^{*0}+\overline{\Xi}^{0})$)/2 in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV. Yields are normalised to the visible cross section for each V0A multiplicity event class and to NSD events for the minimum bias case.

Mean $p_{\rm T}$ of $\Xi^{*0}$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV for each V0A multiplicity event class.

Mean $p_{\rm T}$ of $\Xi^{*0}$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV for each V0A multiplicity event class.

Integrated yield ratio.

Integrated yield ratio.

Integrated yield ratio.

Integrated yield ratio.

Integrated yield ratio.

Integrated yield ratio.

Integrated yield ratio.

Integrated yield ratio.

Breakdown of systematic uncertainties in percent for the measured fiducial cross section times leptonic branching ratio for Z/gamma* production in the Z/gamma* -> mu+mu- final state at 13TeV.

Measured fiducial cross section times leptonic branching ratio for Z/gamma* production in the combined Z/gamma* -> l+l- (l = e, mu) final state.

Measured total cross section times leptonic branching ratio for Z/gamma* production in the combined Z/gamma* -> l+ l- (l = e, mu) final state.

Measured total cross-section for ttbar events using e-mu events with b-tagged jets.

Correlation coefficients amongst the combined (Z/gamma^* -> l+ l-) where (l = e, mu) fiducial cross section and ttbar total cross section at 13, 8, and 7TeV.

Correlation coefficients amongst the combined (Z/gamma^* -> l+ l-) where (l = e, mu) fiducial cross section and ttbar total cross section at 13, 8, and 7TeV, omitting the common normalisation uncertainty due to the luminosity calibration and beam energy.

Z-boson and ttbar production cross-section single ratios at 13, 8, 7TeV. The beam-energy uncertainty, which largely cancels in the ratios, is included in the systematic uncertainty. In the first column, the total (TOT) cross sections are used for both numerator and denominator. In the second column, the total cross section is used for the numerator while the fiducial (FID) cross section is used for the denominator. In the third column, the fiducial cross sections are used for both numerator and denominator. Apart for the MZ requirement, the kinematic fiducial requirements listed below apply only to the fiducial cross sections. The luminosity uncertainty almost entirely cancels in some of the ratio measurements.

Z-boson and ttbar production cross-section double ratios at 13, 8, 7TeV. The beam-energy uncertainty, which largely cancels in the ratios, is included in the systematic uncertainty. In the first column, the total (TOT) cross sections are used for both ttbar and Z processes. In the second column, the total cross section is used for ttbar while the fiducial (FID) cross section is used for Z. In the third column, the fiducial cross sections are used for both ttbar and Z. Apart for the MZ requirement, the kinematic fiducial requirements listed below apply only to the fiducial cross sections.