$\Upsilon(1S)$ production cross-section in $p$Pb and Pb$p$, as a function of $y^*$. The uncertainty is the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ production cross-section in $p$Pb and Pb$p$, as a function of $p_{T}$. The uncertainty is the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ production cross-section in $p$Pb and Pb$p$, as a function of $y^*$. The uncertainty is the sum in quadrature of the statistical and systematic components.

Scaled $pp$ differential cross-section in $p_{T}$ at $\sqrt{s_{NN}}=8.16$ TeV. The first uncertainty is statistical, the second is systematic, which includes the systematic uncertainty from the $pp$ measurement and that estimated by changing the interpolation function.

Scaled $pp$ differential cross-section in $y$ at $\sqrt{s_{NN}}=8.16$ TeV. The first uncertainty is statistical, the second is systematic, which includes the systematic uncertainty from the $pp$ measurement and that estimated by changing the interpolation function.

$\Upsilon(1S)$ nuclear modification factor, $R_{p{\rm Pb}}^{\Upsilon(1S)}$, in $p$Pb and Pb$p$ as a function of $p_{T}$ integrated over $y^*$ in the range $1.5 < y^* < 4.0$ for $p$Pb and $-5.0 < y^* < -2.5$ for Pb$p$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(1S)$ nuclear modification factor, $R_{p{\rm Pb}}^{\Upsilon(1S)}$, in $p$Pb and Pb$p$ as a function of $y^*$ integrated over $p_{T}$ in the range $0 < p_{T} < 25$ GeV$/c$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ nuclear modification factor, $R_{p{\rm Pb}}^{\Upsilon(2S)}$, in $p$Pb and Pb$p$ as a function of $p_{T}$ integrated over $y^*$ in the range $1.5 < y^* < 4.0$ for $p$Pb and $-5.0 < y^* < -2.5$ for Pb$p$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ nuclear modification factor, $R_{p{\rm Pb}}^{\Upsilon(2S)}$, in $p$Pb and Pb$p$ as a function of $y^*$ integrated over $p_{T}$ in the range $0 < p_{T} < 25$ GeV$/c$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(1S)$ forward-to-backward ratio, $R_{{\rm FB}}^{\Upsilon(1S)}$, as a function of $p_{T}$ integrated over $\vert y^* \vert$ in the range $2.5 < \vert y^* \vert < 4.0$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(1S)$ forward-to-backward ratio, $R_{{\rm FB}}^{\Upsilon(1S)}$, as a function of $\vert y^* \vert$ integrated over $p_{T}$ in the range $0 < p_{T} < 25$ GeV/c. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ to $\Upsilon(1S)$ ratio, in $p$Pb and Pb$p$ as a function of $p_{T}$ integrated over $y^*$ in the range $1.5 < y^* < 4.0$ for $p$Pb and $-5.0 < y^* < -2.5$ for Pb$p$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ to $\Upsilon(1S)$ ratio, in $p$Pb and Pb$p$ as a function of $y^*$ integrated over $p_{T}$ in the range $0 < p_{T} < 25$ GeV/c. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(1S)$ to non-prompt $J/\psi$, in $p$Pb and Pb$p$ as a function of $y^*$ integrated over $p_{T}$ in the range $0 < p_{T} < 25$ GeV/c. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

Spectrum in transverse momentum of electrons created in open heavy flavor decays, for 0-10% centrality.

Spectrum in transverse momentum of electrons created in open heavy flavor decays, for 10-20% centrality.

Spectrum in transverse momentum of electrons created in open heavy flavor decays, for 20-40% centrality.

Spectrum in transverse momentum of electrons created in open heavy flavor decays, for 40-60% centrality.

Spectra in transverse momentum of electrons created in open heavy flavor decays, for 60-92% centrality.

Non photonic electrons differential yield scaled by the number of collisions, as a function of centrality. PT belongs to 0.8-4.0 GeV/c.

Non photonic electrons differential yield scaled by the number of collisions, for minimum bias collisions. PT belongs to 0.8-4.0 GeV/c.

Number of collisions as a fonction of centrality There are only systematic errors due to the calculation of Ncoll.

Number of collisions for minimum bias events There are only systematic errors due to the calculation of Ncoll.

Nuclear overlap function TAA as a function of centrality The axis is labelled with 1/SIG, because the unit of the nuclear overlap function is barn-1 There are only systematic errors due to the calculation of TAA.

Nuclear overlap function TAA for minimum bias events The axis is labelled with 1/SIG, as the unit of the nuclear overlap function is barn-1 There are only systematic errors due to the calculation of TAA.

Differential charm cross section scaled by the number of collisions, at mid rapidity, as a function of centrality The differential cross section per number of collisions D(SIG)/Ncoll is calculated as following : D(SIG)/Ncoll = D(N(cc_bar))/TAA.

Differential charm cross section scaled by the number of collisions, at mid rapidity, for minimum bias events The differential cross section per number of collisions D(SIG)/Ncoll is calculated as following : D(SIG)/Ncoll = D(N(cc_bar))/TAA.

Total charm cross section scaled by the number of collisions, as a function of centrality The total cross section per number of collisions SIG/Ncoll is calculated as following : SIG/Ncoll = N(cc_bar)/TAA.

Total charm cross section scaled by the number of collisions, for minimum bias events The total cross section per number of collisions SIG/Ncoll is calculated as following : SIG/Ncoll = N(cc_bar)/TAA.

Nuclear modification factor as a function of PT, for 10-20% central reactions Note that the systematic error given is related to the the uncertainties in the p+p measurements.An additional systematic error, symmetrical on the + and - side, related to the uncertainties in the Au+Au measurement, is given in the second column. Another, PT-independant, 13%systematic error due to the uncertainty on the overlap function and the Pi0 yield normalization is to add.

Nuclear modification factor as a function of PT, for 20-40% central reactions Note that the systematic error given is related to the the uncertainties in the p+p measurements.An additional systematic error, symmetrical on the + and - side, related to the uncertainties in the Au+Au measurement, is given in the second column. Another, PT-independant, 13%systematic error due to the uncertainty on the overlap function and the Pi0 yield normalization is to add.

Nuclear modification factor as a function of PT, for 40-60% central reactions Note that the systematic error given is related to the the uncertainties in the p+p measurements.An additional systematic error, symmetrical on the + and - side, related to the uncertainties in the Au+Au measurement, is given in the second column. Another, PT-independant, 13%systematic error due to the uncertainty on the overlap function and the Pi0 yield normalization is to add.

Centrality dependence of the $p_{T}$ distribution for $K^{-}$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Errors are statistical only.

Centrality dependence of the $p_{T}$ distribution for protons in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Errors are statistical only.

Centrality dependence of the $p_{T}$ distribution for antiprotons in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Errors are statistical only.

Centrality dependence of the $p_{T}$ distribution for antiprotons in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Errors are statistical only.

Inverse slope parameters for $\pi^{+}$, $\pi^{-}$, $K^{+}$, $K^{-}$, proton, and antiproton for the 0-5%, 40-50%, and 60-92% centrality bins. Errors are statistical only.

Fit paramaters $T_{0}^{+}$ and $T_{0}^{-}$. Errors are statistical only.

Fit paramaters $u_{t}^{+}$ and $u_{t}^{-}$. Errors are statistical only.

Relationship between $N_{part}$ and centrality.

Mean transverse momentum as a function of centrality for pions, kaons, protons,and antiprotons in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Errors are systematic only.

dN/dy as a function of centrality for pions, kaons, protons,and antiprotons in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Errors are systematic only.

Particle ratios of $\pi^{-}$/$\pi^{+}$ for centrality 0-5% and 60-92%. Errors are statistical and systematic.

Particle ratios of $K^{-}$/$K^{+}$ for centrality 0-5% and 60-92%. Errors are statistical and systematic.

Particle ratios of antiprotons/protons for centrality 0-5%, 60-92%, and minimum bias. Errors are statistical and systematic.

Particle ratios of kaons/pions for centrality 0-5% and 60-92%. Errors are statistical and systematic.

Particle ratios of proton/$\pi^{+}$ and antiproton/$\pi^{-}$ for centrality 0-10%, 20-30%, and 60-92%. Errors are statistical and systematic.

Particle ratios of proton/$\pi^{0}$ and antiproton/$\pi^{0}$ for centrality 0-10%, 20-30%, and 60-92%. Errors are statistical and systematic.

Centrality dependence of particle ratios for $\pi^{-}$/$\pi^{+}$, $K^{-}$/$K^{+}$, and antiprotons/protons. Errors are statistical and systematic.

Rcp as a function of centrality 0-10% and 60-92% for ($\pi^{+}$ + $\pi^{-}$)/2. Errors are statistical only.

Rcp as a function of centrality 0-10% and 60-92% for ($K^{+}$ + $K^{-}$)/2. Errors are statistical only.

Rcp as a function of centrality 0-10% and 60-92% for charged (p+pbar)/2. Errors are statistical only.

Rcp as a function of centrality 0-10% and 60-92% for neutral pions. Errors are statistical only.

Integrated Rcp above $p_{T}$ = 1.5 GeV/c as a function of centrality 60-92% for pions, kaons, and (p+pbar)/2. Errors are statistical only.

Bertsch-Pratt radius parameters versus the cube root of the number of participants ($N^{1/3}_{part}$) for $\pi^+\pi^+$ (blue squares) and $\pi^-\pi^-$ (red triangle) for the fits to a core-halo structure by Eq. (2) with statistical error bars and systematic error bands, in the 0.2 GeV/$c$ < $k_T$ < 2.0 GeV/$c$ range, with $<k_T>$ ~ 0.45 GeV/$c$. The solid line shows to $p_0$+$p_1$*$N^{1/3}_{part}$ for both signs.

Bertsch-Pratt radius parameters versus the cube root of the number of participants ($N^{1/3}_{part}$) for $\pi^+\pi^+$ (blue squares) and $\pi^-\pi^-$ (red triangle) for the fits to a core-halo structure by Eq. (2) with statistical error bars and systematic error bands, in the 0.2 GeV/$c$ < $k_T$ < 2.0 GeV/$c$ range, with $<k_T>$ ~ 0.45 GeV/$c$. The solid line shows to $p_0$+$p_1$*$N^{1/3}_{part}$ for both signs.

Direct $\gamma$ invariant yields as a function of transverse momentum for 9 centrality selections and minimum bias Au+AU collisions at $\sqrt{s_{NN}}$ = 200 GeV. Data with no errors represents 90% confidence level upper limit. The bin range is not an uncertainty in the x-axis because the actual uncertainty by having the finite bin width is corrected for by the bin-shift correction. These bins were constructed using the corrected finite values as centers.

Ratio of Au+Au yield to $p$+$p$ yield normalized by the number of binary nucleon collisions as a function of centrality given by $N_{part}$ for direct $\gamma$ and $\pi^0$ yields integrated above 6 GeV/$c$. Data with no errors represents 90% confidence level upper limit.

Correlation function, C(q) for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

Correlation function, C(q) for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

Correlation function, C(q) for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

Correlation functions, C(q), for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

Correlation functions, C(q), for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

Correlation functions, C(q), for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

Correlation functions, C(q), for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

Correlation functions, C(q), for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

Correlation functions, C(q), for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

Correlation functions, C(q), for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

Correlation functions, C(q), for $\pi^+\pi^+$ and $\pi^-\pi^-$ pairs.

$k_T$ dependence of extracted spheroidal source parameters for pion source functions (a) $R_T$ , (b) $a$, and (c) $\lambda$ for most central collisions & perhipheral collisions.

$k_T$ dependence of extracted spheroidal source parameters for pion source functions (a) $R_T$ , (b) $a$, and (c) $\lambda$ for most central collisions & perhipheral collisions.

$k_T$ dependence of extracted spheroidal source parameters for pion source functions (a) $R_T$ , (b) $a$, and (c) $\lambda$ for most central collisions & perhipheral collisions.

$\bar{d}/d$ ratio vs. $p_T$ for minimum bias data.

Transverse momentum dependepnce of $v_2$ for identified particle $\pi^-$.

Transverse momentum dependepnce of $v_2$ for identified particle $K^-$.

Transverse momentum dependepnce of $v_2$ for identified particle $\bar{p}$.

Transverse momentum dependepnce of $v_2$ for identified particle $\pi^+$.

Transverse momentum dependepnce of $v_2$ for identified particle $K$.

Transverse momentum dependepnce of $v_2$ for identified particle $p$.

Transverse momentum dependepnce of $v_2$ for combined particles $\pi^+$ and $\pi^-$.

Transverse momentum dependepnce of $v_2$ for combined particles $K^+$ and $K^-$.

Transverse momentum dependepnce of $v_2$ for combined particles $p$ and $\bar{p}$.

Transverse momentum dependepnce of $v_2$ for combined particles $\pi^+$ and $\pi^-$ scaled by quarks for each particle as motivated by a quark coalescence model.

Transverse momentum dependepnce of $v_2$ for combined particles $K^+$ and $K^-$ scaled by quarks for each particle as motivated by a quark coalescence model.

Transverse momentum dependepnce of $v_2$ for combined particles $p$ and $\bar{p}$ scaled by quarks for each particle as motivated by a quark coalescence model.

Transverse momentum dependence of $v_2$ for negative charged particle distributions at 0-20% centrality.

Transverse momentum dependence of $v_2$ for combined $\pi^-$ and $K^-$ at 0-20% centrality.

Transverse momentum dependence of $v_2$ for $\bar{p}$ at 0-20% centrality.

Transverse momentum dependence of $v_2$ for positive charged particle distributions at 0-20% centrality.

Transverse momentum dependence of $v_2$ for combined $\pi^+$ and $K^+$ at 0-20% centrality.

Transverse momentum dependence of $v_2$ for $p$ at 0-20% centrality.

Transverse momentum dependence of $v_2$ for negative charged particle distributions at 20-40% centrality.

Transverse momentum dependence of $v_2$ for combined $\pi^-$ and $K^-$ at 20-40% centrality.

Transverse momentum dependence of $v_2$ for $\bar{p}$ at 20-40% centrality.

Transverse momentum dependence of $v_2$ for positive charged particle distributions at 20-40% centrality.

Transverse momentum dependence of $v_2$ for combined $\pi^+$ and $K^+$ at 20-40% centrality.

Transverse momentum dependence of $v_2$ for $p$ at 20-40% centrality.

Transverse momentum dependence of $v_2$ for negative charged particle distributions at 40-60% centrality.

Transverse momentum dependence of $v_2$ for combined $\pi^-$ and $K^-$ at 40-60% centrality.

Transverse momentum dependence of $v_2$ for $\bar{p}$ at 40-60% centrality.

Transverse momentum dependence of $v_2$ for positive charged particle distributions at 40-60% centrality.

Transverse momentum dependence of $v_2$ for combined $\pi^+$ and $K^+$ at 40-60% centrality.

Transverse momentum dependence of $v_2$ for $p$ at 40-60% centrality.