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

The total coherent $\mathrm{J}/\psi$ photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$, measured in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$ values used correspond to the center of each rapidity range. The theoretical uncertainties is due to the uncertainties in the photon flux.

The nuclear gluon suppression factor $R_{\mathrm{g}}^{\mathrm{Pb}}$ as a function of Bjorken $x$ extracted from the CMS measurement of the coherent $\mathrm{J}/\psi$ photoproduction in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $x$ values are evaluated at the centers of their corresponding rapidity ranges. The theoretical uncertainties are due to the uncertainties in the photon flux and the impulse approximation model.

The nuclear gluon suppression factor $R_{\mathrm{g}}^{\mathrm{Pb}}$ as a function of Bjorken $x$ extracted from the CMS measurement of the coherent $\mathrm{J}/\psi$ photoproduction in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $x$ values are evaluated at the centers of their corresponding rapidity ranges. The theoretical uncertainties are due to the uncertainties in the photon flux and the impulse approximation model.

The total covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The covariance matrix includes both the experimental and theoretical (photon flux) uncertainties. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The photon flux values from STARLight as a function of rapidity, in different neutron multiplicity classes: 0n0n, 0nXn, and XnXn.

The experimental covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The covariance matrix for the flux in the 0n0n neutron multiplicity class.

The theoretical (photon flux) covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The covariance matrix for the flux in the 0nXn neutron multiplicity class.

The covariance matrix for the flux in the XnXn neutron multiplicity class.

The total covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The covariance matrix includes both the experimental and theoretical (photon flux) uncertainties. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The experimental covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The theoretical (photon flux) covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($-0.9 < y < -0.6$).

$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($-1.2 < y < -0.9$).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($0.0 < y < 0.3$).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($-0.3 < y < 0.0$).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($-0.6 < y < -0.3$).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($-0.9 < y < -0.6$).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($-1.2 < y < -0.9$).

d$N$/d$y$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (V0A multiplicity classes).

d$N$/d$y$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (V0A multiplicity classes).

$\langle p_\mathrm{T}\rangle$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (V0A multiplicity classes).

$\langle p_\mathrm{T}\rangle$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (V0A multiplicity classes).

(d$N$/$\mathrm{d}y$)/(d$N$/$\mathrm{d}y)_{y=0}$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (0-100 %).

(d$N$/$\mathrm{d}y$)/(d$N$/$\mathrm{d}y)_{y=0}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (0-100 %).

$\langle p_\mathrm{T}\rangle_\mathrm{y}$/$\langle p_\mathrm{T}\rangle_\mathrm{y=0}$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (0-100 %).

$\langle p_\mathrm{T}\rangle_\mathrm{y}$/$\langle p_\mathrm{T}\rangle_\mathrm{y=0}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (0-100%).

$y_{asym}$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (V0A multiplicity classes).

$y_{asym}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (V0A multiplicity classes).

$Q_{CP}$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (0-10 %/40-100 %).

$Q_{CP}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (0-10 %/40-100%).

$Q_{CP}$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (10-40 %/40-100 %).

$Q_{CP}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (10-40 %/40-100%).

Transverse mass $m_T$ spectrum at midrapidity of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 158 GeV/c. $m_0$ is the PDG mass of the $\phi$ meson.

Transverse mass $m_T$ spectrum at midrapidity of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 80 GeV/c. $m_0$ is the PDG mass of the $\phi$ meson.

Dependence of the slope parameter $T$ of the transverse momentum spectrum on rapidity $y$, of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 158 GeV/c.

Dependence of the slope parameter $T$ of the transverse momentum spectrum on rapidity $y$, of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 80 GeV/c.

Dependence of the width $\sigma_y$ of the rapidity distributions on $p_T$, of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 158 GeV/c.

Dependence of the width $\sigma_y$ of the rapidity distributions on $p_T$, of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 80 GeV/c.

Rapidity spectrum of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 158 GeV/c.

Rapidity spectrum of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 80 GeV/c.

Rapidity spectrum of $\phi$ mesons produced in minimum bias p+p collisions at beam momentum of 40 GeV/c.

Energy dependence of ratios of total yields of $\phi$ mesons to mean total yields of pions produced in in minimum bias p+p collisions at CERN SPS energies.

Energy dependence of double ratios of total yields of $\phi$ mesons to mean total yields of pions in central Pb+Pb collisions over minimum bias p+p collisions at CERN SPS energies.

Energy dependence of total yields of $\phi$ mesons produced in minimum bias p+p collisions at CERN SPS energies.

Energy dependence of midrapidity yields of $\phi$ mesons produced in minimum bias p+p collisions at CERN SPS energies.

Widths of rapidity distributions of $\phi$ mesons produced in minimum bias p+p collisions at CERN SPS energies, as a function of beam rapidity.

$R_{pPb}$ as a function of $p_{T}$ integrated over rapidity range $-1.8<y*<1.3$ for the 0--90% centrality interval for the three geometric models: (a) Glauber, (b) Glauber-Gribov with $\omega_{\sigma}=0.11$ and (c) Glauber-Gribov with $\omega_{\sigma}=0.2$.

$R_{pPb}$ values as a function of $p_{T}$ in the 0--90% centrality interval averaged over $|y*|<0.5$.The $\langle T_{Pb} \rangle$ value for the ATLAS centrality correction is calculated with the Glauber model. The total systematic uncertainties, which include the uncertainty in $\langle T_{Pb} \rangle$, are indicated by lines of the same colour. Strict quantitative agreement is not expected as each measurement uses different rapidity intervals for the centrality determination and apply different event selection criteria to reject diffractive collisions.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $-1.8<y*<1.3$ for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $-1.8<y*<1.3$ for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $-1.8<y*<1.3$ for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{CP}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $|\eta|<2.3$ for seven centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{CP}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $|\eta|<2.3$ for seven centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{CP}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $|\eta|<2.3$ for seven centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

Double-differential yields, $\frac{d^{2}n}{x_{_F}p_{_T}}$ and $f_n(x_{_F},p_{T})$, for $x_{_F}>0$.

$p_{_T}$ integrated yield $\frac{dn}{dy}$, the inverse slope parameter $T$ and the mean transverse mass $\langle m_{_T}\rangle - m_{\Lambda}$.

$p_{_T}$ integrated yield $\frac{dn}{x_{_F}}$ and the invariant cross section $F(x_{_F})$.

Transverse momentum spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 20 GeV in various rapidity ranges.

Transverse momentum spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 31 GeV in various rapidity ranges.

Transverse momentum spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 31 GeV in various rapidity ranges.

Transverse momentum spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 31 GeV in various rapidity ranges.

Transverse momentum spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 31 GeV in various rapidity ranges.

Transverse momentum spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 40 GeV in various rapidity ranges.

Transverse momentum spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 40 GeV in various rapidity ranges.

Transverse momentum spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 40 GeV in various rapidity ranges.

Transverse momentum spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 40 GeV in various rapidity ranges.

Transverse momentum spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 80 GeV in various rapidity ranges.

Transverse momentum spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 80 GeV in various rapidity ranges.

Transverse momentum spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 80 GeV in various rapidity ranges.

Transverse momentum spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 80 GeV in various rapidity ranges.

Transverse momentum spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 80 GeV in various rapidity ranges.

Transverse momentum spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 158 GeV in various rapidity ranges.

Transverse momentum spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 158 GeV in various rapidity ranges.

Transverse momentum spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 158 GeV in various rapidity ranges.

Transverse momentum spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 158 GeV in various rapidity ranges.

Transverse momentum spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 158 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 20 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 20 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 20 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 20 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 31 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 31 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 31 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 31 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 40 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 40 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 40 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 40 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 80 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 80 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 80 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 80 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 158 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 158 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 158 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 158 GeV in various rapidity ranges.

Transverse mass spectra of $\pi^-$ mesons produced in inelastic $p p$ interactions at 158 GeV in various rapidity ranges.

Inverse slope parameter T as a function of rapidity for the 5 beam momenta.

Rapidity spectra for the 5 beam momenta.

Mean transverse mass as a function of rapidity for the 5 beam momenta.

Numerical values of the parameters fitted to the rapidiy distributions.

The corrected per-event transverse momentum distribution of Z bosons in the centrality region 10-20%.

The corrected per-event transverse momentum distribution of Z bosons in the centrality region 20-40%.

The corrected per-event transverse momentum distribution of Z bosons in the centrality region 40-80%.

Combined results for the centrality (Npart) dependence of Z boson yields divided by Ncoll for the PT range > 0 GeV/c. The systematic error includes the uncertainty in Ncoll.

Combined results for the centrality (Npart) dependence of Z boson yields divided by Ncoll for the PT range 0 to 10 GeV/c. The systematic error includes the uncertainty in Ncoll.

Combined results for the centrality (Npart) dependence of Z boson yields divided by Ncoll for the PT range 10 to 30 GeV/c. The systematic error includes the uncertainty in Ncoll.

Combined results for the centrality (Npart) dependence of Z boson yields divided by Ncoll for the PT range > 30 GeV/c. The systematic error includes the uncertainty in Ncoll.

CHI distribution for mass bin > 1200 GeV.

Centrality Ratio.