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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Integrated away-side $\gamma_{dir}$-h per-trigger $I_{AA}$ and $I_{dA}$ of Au+Au, d+Au, and p+p, as a function of $\xi$.

$I_{AA}$ vs $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7 GeV/$c$, 7-9 GeV/$c$, and 9-12 GeV/$c$.

$I_{AA}$ vs $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7 GeV/$c$, 7-9 GeV/$c$, and 9-12 GeV/$c$.

$I_{AA}$ vs $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7 GeV/$c$, 7-9 GeV/$c$, and 9-12 GeV/$c$.

$I_{AA}$ as a function of $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7, 7-9, and 9-12 GeV/$c$. Three away-side integration ranges are chosen to calculate the per-trigger yield and the corresponding $I_{AA}$, $|\Delta\phi-\pi|<\pi/2$, $|\Delta\phi-\pi|<\pi/3$, and $|\Delta\phi-\pi|<\pi/6$.

$I_{AA}$ as a function of $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7, 7-9, and 9-12 GeV/$c$. Three away-side integration ranges are chosen to calculate the per-trigger yield and the corresponding $I_{AA}$, $|\Delta\phi-\pi|<\pi/2$, $|\Delta\phi-\pi|<\pi/3$, and $|\Delta\phi-\pi|<\pi/6$.

$I_{AA}$ as a function of $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7, 7-9, and 9-12 GeV/$c$. Three away-side integration ranges are chosen to calculate the per-trigger yield and the corresponding $I_{AA}$, $|\Delta\phi-\pi|<\pi/2$, $|\Delta\phi-\pi|<\pi/3$, and $|\Delta\phi-\pi|<\pi/6$.

Ratios of $I_{AA}$ as a function of direct photon $p_T$ for three different away-side integration ranges.

Measured $I_{AA}$ for direct photon $p_T$ of 5-7, 7-9, and 9-12 GeV/$c$, as a function of $\xi$, are compared with theoretical model calculations.

Measured $I_{AA}$ for direct photon $p_T$ of 5-7, 7-9, and 9-12 GeV/$c$, as a function of $\xi$, are compared with theoretical model calculations.

Measured $I_{AA}$ for direct photon $p_T$ of 5-7, 7-9, and 9-12 GeV/$c$, as a function of $\xi$, are compared with theoretical model calculations.

$p_{\rm T}$-differential yields of $\Lambda$(1520) (sum of particle and anti-particle states) in p--Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ $\mathrm{=}$ 5.02 TeV in multiplicity interval 20--40\%. The uncertainty 'sys,$p_{\rm T}$-correlated' indicates the systematic uncertainty after removing the contributions of $p_{\rm T}$-uncorrelated uncertainty.

$p_{\rm T}$-differential yields of $\Lambda$(1520) (sum of particle and anti-particle states) in p--Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ $\mathrm{=}$ 5.02 TeV in multiplicity interval 40--60\%. The uncertainty 'sys,$p_{\rm T}$-correlated' indicates the systematic uncertainty after removing the contributions of $p_{\rm T}$-uncorrelated uncertainty.

$p_{\rm T}$-differential yields of $\Lambda$(1520) (sum of particle and anti-particle states) in p--Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ $\mathrm{=}$ 5.02 TeV in multiplicity interval 60--100\%. The uncertainty 'sys,$p_{\rm T}$-correlated' indicates the systematic uncertainty after removing the contributions of $p_{\rm T}$-uncorrelated uncertainty.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to charged $\pi$ ($\pi^{+} + \pi^{-}$) at midrapidity as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in inelastic pp collisions at $\sqrt{s}$ $\mathrm{=}$ 7 TeV.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to charged $\pi$ ($\pi^{+} + \pi^{-}$) as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in p--Pb collisions at $\sqrt{s}$ $\mathrm{=}$ 5.02 TeV. The uncertainty 'sys,uncorrelated' indicates the multiplicity uncorrelated systematic uncertainty.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to K$^{-}$ at midrapidity as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in inelastic pp collisions at $\sqrt{s}$ $\mathrm{=}$ 7 TeV.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to K$^{-}$ as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in p--Pb collisions at $\sqrt{s}$ $\mathrm{=}$ 5.02 TeV. The uncertainty 'sys,uncorrelated' indicates the multiplicity uncorrelated systematic uncertainty.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to its ground state particle $\Lambda$ ($\Lambda + \bar{\Lambda}$) at midrapidity as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in inelastic pp collisions at $\sqrt{s}$ $\mathrm{=}$ 7 TeV.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to its ground state particle $\Lambda$ ($\Lambda + \bar{\Lambda}$) as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in p--Pb collisions at $\sqrt{s}$ $\mathrm{=}$ 5.02 TeV. The uncertainty 'sys,uncorrelated' indicates the multiplicity uncorrelated systematic uncertainty.

Direct photon spectra(Nucl. Part. Phys. 23, A1 (1997)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in p+p at $\sqrt{s_{NN}}$= 63 GeV.

Direct photon spectra(Sov. J. Nucl. Phys. 51, 836 (1990)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in p+p at $\sqrt{s_{NN}}$= 63 GeV.

Direct photon spectra normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Au+Au at $\sqrt{s_{NN}}$= 62.4 GeV.

Direct photon spectra normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Au+Au at $\sqrt{s_{NN}}$= 39.0 GeV.

Direct photon spectra(Physical Review C91, 064904 (2015)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Au+Au at $\sqrt{s_{NN}}$= 200 GeV.

Direct photon spectra(Physical Review Letters 104, 132301 (2010)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Au+Au at $\sqrt{s_{NN}}$= 200 GeV.

Direct photon spectra(Physical Review Letters 109, 152302 (2012)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Au+Au at $\sqrt{s_{NN}}$= 200 GeV.

Direct photon spectra(Physical Review C98, 054902 (2018)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Cu+Cu at $\sqrt{s_{NN}}$= 200 GeV.

Direct photon spectra(Physics Letters B754 235 (2016)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Pb+Pb at $\sqrt{s_{NN}}$= 2760 GeV.

Number of binary collisions N_{coll} vs. Charged particle multiplicity in Au+Au at various cm energies.

Charged particle multiplicity vs. centrality in Au+Au at various cm energies.

Number of binary collisions N_{coll} vs. centrality in Au+Au at various cm energies.

Integrated direct photon yield vs. Charged particle multiplicity in various systems at various cm energies.

Charged particle multiplicity vs. centrality in various systems at various cm energies.

Integrated direct photon yield vs. centrality in various systems at various cm energies.

Integrated direct photon yield vs. Charged particle multiplicity in various systems at various cm energies.

Charged particle multiplicity vs. centrality in various systems at various cm energies.

Integrated direct photon yield vs. centrality in various systems at various cm energies.

Integrated per-trigger yields as a function of associated particle angle relative to $\psi_2$ event plane $\phi^a - \psi_2$ integrated over the near side $|\Delta\phi|<\pi/4$.

Integrated per-trigger yields as a function of associated particle angle relative to $\psi_2$ event plane $\phi^a - \psi_2$ integrated over the away side $|\Delta\phi|<\pi/4$.

Integrated per-trigger yields as a function of associated particle angle relative to $\psi_2$ event plane for near-side $|\Delta\phi|<\pi/4$ and for away-side $|\Delta\phi|<\pi/4$.

Integrated per-trigger yields as a function of associated particle angle relative to $\psi_2$ event plane for near-side $|\Delta\phi|<\pi/4$ and for away-side $|\Delta\phi|<\pi/4$.

Per-trigger yields $Y(\Delta\phi)$ of dihadrons pairs measured in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200~{\rm GeV}$ after subtracting the underlying event model with several $p_T$ selections. Systematic uncertainties due to track matching are given.

Per-trigger yields $Y(\Delta\phi)$ of dihadrons pairs measured in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200~{\rm GeV}$ after subtracting the underlying event model with several $p_T$ selections. Systematic uncertainties due to track matching are given.

$\Psi_{2}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $1<p_{T}^{assoc}<2$ GeV/c.

$\Psi_{2}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{2}$-dependent $C(\Delta\phi)$ for $4<p_{T}^{trig}<10$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $1<p_{T}^{assoc}<2$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $4<p_{T}^{trig}<10$ GeV/c and $1<p_{T}^{assoc}<2$ GeV/c.

$\Psi_{2}$-dependent $Y(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $1<p_{T}^{assoc}<2$ GeV/c.

$\Psi_{2}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{2}$-dependent $C(\Delta\phi)$ for $4<p_{T}^{trig}<10$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $1<p_{T}^{assoc}<2$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $4<p_{T}^{trig}<10$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_2$-dependent C for 2 $< p_T^t <$ 4 and 1 $< p_T^a <$ 2 GeV/c.

$\Psi_2$-dependent C for 2 $< p_T^t <$ 4 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_2$-dependent C for 4 $< p_T^t <$ 10 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_2$-dependent Y for 2 $< p_T^t <$ 4 and 1 $< p_T^a <$ 2 GeV/c.

$\Psi_2$-dependent Y for 2 $< p_T^t <$ 4 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_2$-dependent Y for 4 $< p_T^t <$ 10 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_3$-dependent C for 2 $< p_T^t <$ 4 and 1 $< p_T^a <$ 2 GeV/c.

$\Psi_3$-dependent C for 2 $< p_T^t <$ 4 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_3$-dependent C for 4 $< p_T^t <$ 10 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_3$-dependent Y for 2 $< p_T^t <$ 4 and 1 $< p_T^a <$ 2 GeV/c.

$\Psi_3$-dependent Y for 2 $< p_T^t <$ 4 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_3$-dependent Y for 4 $< p_T^t <$ 10 and 2 $< p_T^a <$ 4 GeV/c.

Multiplicity distribution in the pseudorapidity region -3.4 to 3.4 for NSD collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 5.0 for NSD collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -2.0 to 2.0 for NSD collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -2.4 to 2.4 for NSD collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -3.0 to 3.0 for NSD collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 3.4 for NSD collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 5.0 for NSD collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -2.0 to 2.0 for NSD collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -2.4 to 2.4 for NSD collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -3.0 to 3.0 for NSD collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 3.4 for NSD collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 5.0 for NSD collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -2.0 to 2.0 for INEL collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -2.4 to 2.4 for INEL collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -3.0 to 3.0 for INEL collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 3.4 for INEL collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 5.0 for INEL collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -2.0 to 2.0 for INEL collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -2.4 to 2.4 for INEL collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -3.0 to 3.0 for INEL collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 3.4 for INEL collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 5.0 for INEL collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -2.0 to 2.0 for INEL collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -2.4 to 2.4 for INEL collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -3.0 to 3.0 for INEL collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 3.4 for INEL collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 5.0 for INEL collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -2.0 to 2.0 for INEL>0 collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -2.4 to 2.4 for INEL>0 collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -3.0 to 3.0 for INEL>0 collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 3.4 for INEL>0 collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 5.0 for INEL>0 collisions at a centre-of-mass energy of 900 GeV.

Multiplicity distribution in the pseudorapidity region -2.0 to 2.0 for INEL>0 collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -2.4 to 2.4 for INEL>0 collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -3.0 to 3.0 for INEL>0 collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 3.4 for INEL>0 collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 5.0 for INEL>0 collisions at a centre-of-mass energy of 7000 GeV.

Multiplicity distribution in the pseudorapidity region -2.0 to 2.0 for INEL>0 collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -2.4 to 2.4 for INEL>0 collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -3.0 to 3.0 for INEL>0 collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 3.4 for INEL>0 collisions at a centre-of-mass energy of 8000 GeV.

Multiplicity distribution in the pseudorapidity region -3.4 to 5.0 for INEL>0 collisions at a centre-of-mass energy of 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.0 to 2.0 in NSD collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.4 to 2.4 in NSD collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.0 to 3.0 in NSD collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 3.4 in NSD collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 5.0 in NSD collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.0 to 2.0 in NSD collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.4 to 2.4 in NSD collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.0 to 3.0 in NSD collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 3.4 in NSD collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 5.0 in NSD collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.0 to 2.0 in NSD collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.4 to 2.4 in NSD collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.0 to 3.0 in NSD collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 3.4 in NSD collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 5.0 in NSD collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.0 to 2.0 in INEL collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.4 to 2.4 in INEL collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.0 to 3.0 in INEL collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 3.4 in INEL collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 5.0 in INEL collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.0 to 2.0 in INEL collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.4 to 2.4 in INEL collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.0 to 3.0 in INEL collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 3.4 in INEL collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 5.0 in INEL collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.0 to 2.0 in INEL collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.4 to 2.4 in INEL collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.0 to 3.0 in INEL collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 3.4 in INEL collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 5.0 in INEL collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.0 to 2.0 in INEL>0 collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.4 to 2.4 in INEL>0 collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.0 to 3.0 in INEL>0 collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 3.4 in INEL>0 collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 5.0 in INEL>0 collisions at center-of-mass energy 900 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.0 to 2.0 in INEL>0 collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.4 to 2.4 in INEL>0 collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.0 to 3.0 in INEL>0 collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 3.4 in INEL>0 collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 5.0 in INEL>0 collisions at center-of-mass energy 7000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.0 to 2.0 in INEL>0 collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -2.4 to 2.4 in INEL>0 collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.0 to 3.0 in INEL>0 collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 3.4 in INEL>0 collisions at center-of-mass energy 8000 GeV.

Set of fit parameters for the multiplicity distribution in pseudorapidity range -3.4 to 5.0 in INEL>0 collisions at center-of-mass energy 8000 GeV.

The double-differential cross sections for $J/\psi$ production from $b$-decays as a function of transverse momentum for the rapidity range $1.5 < y^* < 4.0$ in the nucleon-nucleon centre-of-mass frame. The first quoted uncertainty indicates the bin-by-bin correlated systematic uncertainty and the second is the bin-by-bin uncorrelated systematic uncertainty.

The double-differential cross sections for prompt $J/\psi$ production, assuming no polarisation, as a function of transverse momentum for the rapidity range $-5.0 < y^* < -2.5$ in the nucleon-nucleon centre-of-mass frame. The first quoted uncertainty indicates the bin-by-bin correlated systematic uncertainty and the second is the bin-by-bin uncorrelated systematic uncertainty.

The double-differential cross sections for $J/\psi$ production from $b$-hadron decays as a function of transverse momentum for the rapidity range $-5.0 < y^* < -2.5$ in the nucleon-nucleon centre-of-mass frame. The first quoted uncertainty indicates the bin-by-bin correlated systematic uncertainty and the second is the bin-by-bin uncorrelated systematic uncertainty.

Fraction of $J/\psi$ from $b$-hadron decay in bins of transverse momentum and rapidity in the nucleon-nucleon rest frame, assuming no polarisation, in the proton-lead beam configuration corresponding to rapidity range $1.5< y^* <4.0$ in the nucleon-nucleon centre-of-mass frame. The uncertainty is uncorrelated between bins.

Fraction of $J/\psi$ from $b$-hadron decay in bins of transverse momentum and rapidity, assuming no polarisation, in the proton-lead beam configuration corresponding to the rapidity range $-5.0 < y^* < -2.5$ in the nucleon-nucleon centre-of-mass frame. The uncertainty is uncorrelated between bins.

The nuclear modification factor $R_{p \mathrm{Pb}}$ of prompt $J/\psi$ for transverse momentum smaller than 14 GeV/c as a function of transverse momentum integrated over the rapidity ranges $1.5 < y^* < 4.0$ for $p$Pb and $-5.0 < y^* < -2.5$ for Pb$p$. The quoted uncertainties are the quadratic sums of statistical and systematic uncertainties.

The nuclear modification factor $R_{p\mathrm{Pb}}$ of prompt $J/\psi$ integrated over transverse momentum smaller than 14 GeV/$c$ as a function of rapidity in the nucleon-nucleon centre-of-mass frame ($y^*$). The quoted uncertainties are the quadratic sums of statistical and systematic uncertainties.

The nuclear modification factor $R_{p\mathrm{Pb}}$ of $J/psi$ from the decay of $b$-hadrons for transverse momentum smaller than 14 GeV/$c$ as a function of transverse momentum integrated over the rapidity range $1.5 < y^* < 4.0$ for $p$Pb and $-5.0 < y^* < -2.5$ for Pb$p$. The quoted uncertainties are the quadratic sums of statistical and systematic uncertainties.

The nuclear modification factor $R_{p\mathrm{Pb}}$ of $J/\psi$ from decay of $b$-hadrons integrated over transverse momentum smaller than 14 GeV/$c$ as a function of rapidity in the nucleon-nucleon centre-of-mass frame. The quoted uncertainties are the quadratic sums of statistical and systematic uncertainties.

The FORWARD/BACKWARD production ratio $R_{\mathrm{FB}}$ of prompt $J/\psi$ as a function of transverse momentum integrated over $|y^*|$ in the range $2.5 < |y^*| < 4.0$. The quoted uncertainties are the quadratic sums of statistical and systematic uncertainties.

The FORWARD/BACKWARD production ratio $R_{\mathrm{FB}}$ of prompt $J/\psi$ as a function of rapidity $y^*$ in the nucleon-nucleon centre-of-mass frame integrated over transverse momentum smaller than 14 GeV/$c$. The quoted uncertainties are the quadratic sums of statistical and systematic uncertainties.

The FORWARD/BACKWARD production ratio $R_{\mathrm{FB}}$ of $J/\psi$ from $b$-hadron decays as a function of transverse momentum integrated over $|y^*|$ in the range $2.5 < |y^*| < 4.0$. The quoted uncertainties are the quadratic sums of statistical and systematic uncertainties.

The FORWARD/BACKWARD production ratio $R_{\mathrm{FB}}$ of $J/\psi$ from $b$-hadron decays as a function of rapidity $y^*$ in the nucleon-nucleon centre-of-mass frame integrated over transverse momentum smaller than 14 GeV/$c$. The quoted uncertainties are the quadratic sums of statistical and systematic uncertainties.