The minimum $\chi^2$ and the parameters $T_{fo}$ and $\beta^{max}_T$ for each of the five centrality selections. The best fit parameters are determined by averaging all parameter pairs within the 1$\sigma$ contour. The errors correspond to the standard deviation of the parameter pairs within the 1$\sigma$ $\chi^2$ contour.

Fit parameters for each particle species using Equations $\frac{dN}{dy}$ = $C \cdot$ ($N_{part}$)$^{\alpha_{part}}$ and $\frac{dN}{dy}$ = $C^{\prime} \cdot$ ($N_{coll}$)$^{\alpha_{coll}}$.

Values of the parameters $n_{pp}$ and $x$ from fitting Equation $R \equiv \frac{dN/dy}{N_{part}}$ = (1-$x$) $\cdot n_{pp} \frac{1}{2}$ + $x \cdot n_{pp} \frac{N_{coll}}{N_{part}}$ = $n_{pp}$ [$\frac{1}{2}$ + $x$($\frac{N_{coll}}{N_{part}}$ - $\frac{1}{2})$]to the observed $dN$/$dy$ per $N_{part}$.

Invariant yields for $\pi^{\pm}$ measured in minimum-bias events at midrapidity and normalized to one rapidity unit.

Invariant yields for $K^{\pm}$ measured in minimum-bias events at midrapidity and normalized to one rapidity unit.

Invariant yields for (anti)$p$ measured in minimum-bias events at midrapidity and normalized to one rapidity unit.

Pion invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

Pion invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

Kaon invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

Kaon invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

Minimum-bias invariant yields for $\pi^{\pm}$ in equal $p_T$ bins.

Minimum-bias invariant yields for $K^{\pm}$ in equal $p_T$ bins.

Minimum-bias invariant yields for (anti)$p$ in equal $p_T$ bins.

Pion invariant yields in each event centrality normalized to one rapidity unit at midrapidity.

Pion invariant yields in each event centrality normalized to one rapidity unit at midrapidity.

Kaon invariant yields in each event centrality normalized to one rapidity unit at midrapidity.

(Anti)proton invariant yields in each event centrality normalized to one rapidity unit at midrapidity.

(Anti)proton invariant yields in each event centrality normalized to one rapidity unit at midrapidity.

The integrated mean $p_T$ for pions, kaons, and (anti)protons produced in the five different classes of event centrality.

The expansion parameters $T_{fo}$ and $\beta^{max}_T$ as a function of the number of participants.

$v(Q)$, for central events, as a function of the azimuthal coverage of the detector for data. (Angles in degrees.)

The effect on $v(Q)$ of varying the acceptance in $\varphi_r$. The solid curve shows the expectation from global charge conservation. The 10% most central events are used for data and RQMD. (Angles in degrees.)

The $M_{e_T}$ distributions for four different centrality classes. The curves are the random baseline mixed event distributions. The source of differences in the data and mixed event distributions are addressed in the text.

The residual distribution between the data and mixed event $M_{e_T}$ spectra as a function of $M_{e_T}$ in units of standard deviation. The total ${\chi}^2$ and the number of degrees of freedom for the $0-5\%$, $0-10\%$, $10-20\%$, $20-30\%$ centrality classes are 310.0/32, 896.4/36, 678.7/47, and 553.9/53, respectively. A large fraction of the residual contributions are due to the effects of cluster merging.

Ratio of the electron data to the calculated background as a function of $p_T$ in minimum bias collisions (upper panel) and relative contributions to the background from various sources (lower panel). The curves for $\omega$ and $\phi$ show the sum of the Dalitz and the di-electron decay modes.

The background-subtracted electron spectra for central (0-10%) collisions compared with the expected contributions from open charm decays. Also shown, for central collisions only, are the expected contribution from bottom decays (dashed) and the conversion electron spectrum from a direct phonon prediction (dotted).

The background-subtracted electron spectra for minimum bias (0-92%) collisions compared with the expected contributions from open charm decays. Also shown, for central collisions only, are the expected contribution from bottom decays (dashed) and the conversion electron spectrum from a direct phonon prediction (dotted).

Single electron cross sections $\frac{d\sigma_e}{dy}|_{y=0}$ of this measurement and ISR experiments [9,11,31] are displayed (bottom of figure, right-hand scale) with charm decay contributions calculated with PYTHIA. Open and filled symbols for 1.0 < $p_T$ < 1.4 GeV/$c$ and $p_T$ > 1.4 GeV/$c$, respectively. The derived charm cross section of this measurement is compared with charm cross sections from SPS/FNAL experiments (top of figure, left-hand scale). The thick curve and the shaded band represent the charm cross section in the PYTHIA model and in a NLO pQCD calculation [30], respectively.

Single electron cross sections $\frac{d\sigma_e}{dy}|_{y=0}$ of this measurement and ISR experiments [9,11,31] are displayed (bottom of figure, right-hand scale) with charm decay contributions calculated with PYTHIA. Open and filled symbols for 1.0 < $p_T$ < 1.4 GeV/$c$ and $p_T$ > 1.4 GeV/$c$, respectively. The derived charm cross section of this measurement is compared with charm cross sections from SPS/FNAL experiments (top of figure, left-hand scale). The thick curve and the shaded band represent the charm cross section in the PYTHIA model and in a NLO pQCD calculation [30], respectively.

$v_2$ vs Fixed $p_T$ for several centrality selections. [F] and [A] follow the notation Fig. 2. Systematic errors are estimated to be $\sim 5$%.

$v_2$ vs Assorted $p_T$ for several centrality selections. [F] and [A] follow the notation Fig. 2. Systematic errors are estimated to be $\sim 5$%.

$v_2$ vs Assorted $p_T$ for several centrality selections. [F] and [A] follow the notation Fig. 2. Systematic errors are estimated to be $\sim 5$%.

The centrality dependence of $A_2$ for two different $p_T$ selections. [F] and [A] follow the notation in Fig. 2. Systematic errors are estimated to be $\sim 10$%, dominated by the normalization of the correction function and the model determination of$ \varepsilon$.

The centrality dependence of $A_2$ for two different $p_T$ selections. [F] and [A] follow the notation in Fig. 2. Systematic errors are estimated to be $\sim 10$%, dominated by the normalization of the correction function and the model determination of$ \varepsilon$.

$\Delta\phi$ distributions for meson triggers with 2.5 < $p_T$ < 4.0 GeV/$c$ and associated charged hadrons with 1.7 < $p_T$ < 2.5 GeV/$c$ for 60-70% centrality in Au+Au collisions.

$\Delta\phi$ distributions for baryon triggers with 2.5 < $p_T$ < 4.0 GeV/$c$ and associated charged hadrons with 1.7 < $p_T$ < 2.5 GeV/$c$ for 60-70% centrality in Au+Au collisions.

Combinatorial Background Levels from Figure 1 Top Panels.

Combinatorial Background Levels from Figure 1 Bottom Panels.

Trigger Particle $v_2$ values.

Partner Particle $v_2$ values.

Yield per trigger for associated charged hadrons between 1.7 < $p_T$ < 2.5 GeV/$c$ for the near- (top) and away- (bottom) side jets.

Yield per trigger for associated charged hadrons between 1.7 < $p_T$ < 2.5 GeV/$c$ for the near- (top) and away- (bottom) side jets.

Uncertainties in Figure 2 as a function of $N_{part}$ and trigger type for AuAu and $d$+Au.

Uncertainties in Figure 2 as a function of $N_{part}$ and trigger type for AuAu and $d$+Au.

$p_T$ spectra of the near side associated charged hadrons corrected to the full jet yield for meson triggers at 2.5 < $p_T$ < 4.0 GeV/$c$ and $|\eta|$ < 0.35 for six centralities in Au+Au and $d$+Au collisions.

$p_T$ spectra of the near side associated charged hadrons corrected to the full jet yield for meson triggers at 2.5 < $p_T$ < 4.0 GeV/$c$ and $|\eta|$ < 0.35 for six centralities in Au+Au and $d$+Au collisions.

$p_T$ spectra of the near side associated charged hadrons corrected to the full jet yield for baryon triggers at 2.5 < $p_T$ < 4.0 GeV/$c$ and $|\eta|$ < 0.35 for six centralities in Au+Au, $d$+Au and $p$+$p$ collisions.

Inverse slopes from the fits in Figure 3.

Inverse slopes from the fits in Figure 3.

Inverse slopes from the fits in Figure 3.

Differential anisotropy $v_2$($p_T$) for charged hadrons in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV obtained via reaction plane methods

Differential anisotropy as a function of centrality for $\sqrt{s_{NN}}$ = 62.4 GeV/$c$

Differential anisotropy as a function of centrality for $\sqrt{s_{NN}}$ = 130 GeV/$c$

Differential anisotropy as a function of centrality for $\sqrt{s_{NN}}$ = 200 GeV/$c$

Differential anisotropy $v_2$($p_T$) for centrality selection 13 - 26%.

Differential anisotropy $v_2$($p_T$) for centrality selection 13 - 26%.

Differential $v_2$ vs. $\sqrt{S_{NN}}$ for charged hadrons in nucleus-nucleus collisions. Results are shown for the centrality cut of 13 - 26% and $p_T$ selections of 1.75 GeV/$c$ and 0.65 GeV/$c$.

Centrality dependence of $⟨p^{trunc}_T⟩$, the average $p_T$ of charged particles with $p_T$ above a threshold $p^{min}_T$ minus the threshold $p^{min}_T$. Shown are values for two $p^{min}_T$ cuts, one at $p_T >$ 0.5 GeV/$c$ representing all data presented in Fig. 2 and the other one at $p_T >$ 2 GeV/$c$.

Nuclear modification factor ($R_{AA}$) for the 60-80%, 30-60%, 15-30%, and 5-15% centrality selections compared to the one for the most central sample (0-5%). Due to insufficient statistics $R_{AA}$ is not shown for the 80-92% sample.

The top panel gives the nuclear modification factor $R_{AA}$ for three exclusive $p_T$ regions as a function of the centrality of the collision. The lower panel shows essentially the same quantity but normalized to the number of participant pairs rather than to the number of binary collisions. The dotted line indicates the expectation for scaling with the number of binary collisions (top) or with the number of participants (bottom). Only statistical errors are shown. The systematic error on the scale and the centrality dependence are identical to the errors shown in Fig. 5. These errors are correlated, i.e. take their maximum or minimum value simultaneously for all centrality and $p_T$ selections. In addition, there are also $p_T$ dependent systematic errors, which are given in Table 1. The systematic errors do not alter the trends in the data.

The top panel shows the measured $R_{side}$ from identical pions for PHENIX. The bottom panel shows the ratio $R_{out}/R_{side}$ as a function of $k_T$. Longitudinal Co-Moving System (LCMS) frame for $\pi^-$

The top panel shows the measured $R_{side}$ from identical pions for PHENIX. The bottom panel shows the ratio $R_{out}/R_{side}$ as a function of $k_T$. Pair Center-of-Mass System (PCMS) frame for $\pi^+$

The top panel shows the measured $R_{side}$ from identical pions for PHENIX. The bottom panel shows the ratio $R_{out}/R_{side}$ as a function of $k_T$. Pair Center-of-Mass System (PCMS) frame for $\pi^-$

Pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in pp collisions at $\sqrt{s} = 5020~\mathrm{GeV}$.

Pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in Pb-p collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in pp collisions at various LHC energies.

$\langle N_{\rm \gamma} \rangle$ measured within $2.3<\eta_{\rm lab}<3.9$ as a function of collision energy in pp collisions.

Centrality-dependent pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$ (Centrality estimator - CL1).

Centrality-dependent pseudorapidity distribution of inclusive photons measured within $2.3<\eta_{\rm lab}<3.9$ in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$ (Centrality estimator - ZNA).

$\langle N_{\rm \gamma} \rangle$ measured within $2.3<\eta_{\rm lab}<3.9$ as a function of $\langle N_{\rm part} \rangle$ for CL1 centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

$\langle N_{\rm \gamma} \rangle$ measured within $2.3<\eta_{\rm lab}<3.9$ as a function of $\langle N_{\rm part} \rangle$ for ZNA (high-$p_{\rm T}$) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

$\langle N_{\rm \gamma} \rangle$ measured within $2.3<\eta_{\rm lab}<3.9$ as a function of $\langle N_{\rm part} \rangle$ for ZNA (Pb-side) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm part} \rangle$ for CL1 centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm part} \rangle$ for ZNA (high-$p_{\rm T}$) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm part} \rangle$ for ZNA (Pb-side) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm \gamma} \rangle$ for CL1 centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm \gamma} \rangle$ for ZNA (high-$p_{\rm T}$) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.

Values of $\langle N_{\rm \gamma} \rangle$/$\langle N_{\rm part} \rangle$ as a function of $\langle N_{\rm \gamma} \rangle$ for ZNA (Pb-side) centrality estimator in p-Pb collisions at $\sqrt{s_\mathrm{NN}} = 5020~\mathrm{GeV}$.