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

The $p_{T}$ spectra of $\pi^{-}$ measured at midrapidity (|y|<0.1) in Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV. Spectra are plotted for nine centrality classes, with some spectra multiplied by a scale factor to improve clarity, as indicatedin the legend

The $p_{T}$ spectra of $K^{+}$ measured at midrapidity (|y|<0.1) in Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV. Spectra are plotted for nine centrality classes, with some spectra multiplied by a scale factor to improve clarity, as indicatedin the legend

The $p_{T}$ spectra of $K^{-}$ measured at midrapidity (|y|<0.1) in Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV. Spectra are plotted for nine centrality classes, with some spectra multiplied by a scale factor to improve clarity, as indicatedin the legend

Average $p_{T}$ of $\pi^{+}$ as a function of number of participant for Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV.

Average $p_{T}$ of $\pi^{-}$ as a function of number of participant for Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV.

Average $p_{T}$ of $K^{+}$ as a function of number of participant for Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV.

Average $p_{T}$ of $K^{-}$ as a function of number of participant for Au+Au collisions at $\sqrt{s_{NN}}$= 14.5 GeV.

Average $p_{T}$ of p as a function of number of participant for Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV.

Average $p_{T}$ of p-bar as a function of number of participant for Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV.

dN/dy of $\pi^{+}$ scaled by 0.5*$N_{part}$ as a function of number of participant for Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV.

dN/dy of $\pi^{-}$ scaled by 0.5*$N_{part}$ as a function of number of participant for Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV.

dN/dy of $K^{+}$ scaled by 0.5*$N_{part}$ as a function of number of participant for Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV.

dN/dy of $K^{-}$ scaled by 0.5*$N_{part}$ as a function of number of participant for Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV.

dN/dy of proton scaled by 0.5*$N_{part}$ as a function of number of participant for Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV.

dN/dy of p-bar scaled by 0.5*$N_{part}$ as a function of number of participant for Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV.

Kinetic freeze-out temperature as a function of number of participant for Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV.

Velocity as a function of number of participant for Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV.

The event plane resolution calculated for Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV as a function of centrality.

Inclusive charged particle elliptic flow v2 at mid-pseudorapidity (|y| <1.0) as a function of $p_{T}$ for 10-20% centrality in Au + Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV.

Inclusive charged particle elliptic flow v2 at mid-pseudorapidity (|y| <1.0) as a function of $p_{T}$ for 20-30% centrality in Au + Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV.

Inclusive charged particle elliptic flow v2 at mid-pseudorapidity (|y| <1.0) as a function of $p_{T}$ for 30-40% centrality in Au + Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV.

Inclusive charged particle elliptic flow v2 at mid-pseudorapidity (|y| <1.0) as a function of transverse momentum $p_{T}$ for six centrality classes, obtained using the $\eta$-sub event plane method in Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV.

Inclusive charged particle elliptic flow v2 at mid-pseudorapidity (|y| <1.0) as a function of $p_{T}$-integrated v2($\eta$) for six centrality classes, obtained using the $\eta$-sub event plane method in Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV.

The ratio inclusive charged particle elliptic flow v2 over root-mean-square participant eccentricity $Epart_{2}$ at mid-pseudorapidity as a function of $p_{T}$ for 10–20%, 30–40%, and 50–60% collision centralities in Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV.

Summary of centrality bins, average number of participants $N_{part}$, number of binary collisions $N_{coll}$, reaction plane eccentricity eRP, participant eccentricity epart, root-mean-square of the participant eccentricity epart{2}, and transverse area $S_{part}$ from MC Glauber simulations at $\sqrt{s_{NN}}$ = 14.5 GeV.

The inclusive charged particle elliptic flow v2($\eta$-sub) versus pseudorapidity $\eta$ at mid-pseudorapidity for $\sqrt{s_{NN}}$ = 14.5 GeV.

Rapidity-odd charged particles directed flow v1 as a function of $p_{T}$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7 GeV for 0–10%, 10–40% and 40–80% centrality intervals.

Rapidity-odd charged particles directed flow v1 as a function of $p_{T}$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 11.5 GeV for 0–10%, 10–40% and 40–80% centrality intervals.

Rapidity-odd charged particles directed flow v1 as a function of $p_{T}$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV GeV for 0–10%, 10–40% and 40–80% centrality intervals.

Rapidity-odd charged particles directed flow v1 as a function of $p_{T}$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV for 0–10%, 10–40% and 40–80% centrality intervals.

Rapidity-odd charged particles directed flow v1 as a function of $p_{T}$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 27.0 GeV for 0–10%, 10–40% and 40–80% centrality intervals.

Rapidity-odd charged particles directed flow v1 as a function of $p_{T}$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 39.0 GeV for 0–10%, 10–40% and 40–80% centrality intervals.

Rapidity-odd charged particles directed flow v1 as a function of pseudorapidity $\eta$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7 GeV for 0–10%, 10–40% and 40–80% centrality intervals.

Rapidity-odd charged particles directed flow v1 as a function of pseudorapidity $\eta$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 11.5 GeV for 0–10%, 10–40% and 40–80% centrality intervals.

Rapidity-odd charged particles directed flow v1 as a function of pseudorapidity $\eta$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 14.5 GeV for 0–10%, 10–40% and 40–80% centrality intervals.

Rapidity-odd charged particles directed flow v1 as a function of pseudorapidity $\eta$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 19.6 GeV for 0–10%, 10–40% and 40–80% centrality intervals.

Rapidity-odd charged particles directed flow v1 as a function of pseudorapidity $\eta$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 27.0 GeV for 0–10%, 10–40% and 40–80% centrality intervals.

Rapidity-odd charged particles directed flow v1 as a function of pseudorapidity $\eta$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 39.0 GeV for 0–10%, 10–40% and 40–80% centrality intervals.

Rapidity-odd charged particles directed flow v1 as a function of pseudorapidity $\eta$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7 – 39 GeV for 30-60% centrality intervals.

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

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.