The TPC mid-rapidity multiplicity distributions ($|\eta|$ < 0.5) for the corresponding E-FTPC selected centrality bins.

Uncorrected charged particle multiplicity distribution measured in the TPC within $|\eta|$ < 0.5 in pp collisions at 200 GeV.

The ratio of the transverse overlap area $S_{⊥}$ to $(N_{part}/2)^{2/3}$ versus $(N_{part}/2)^{2/3}$.

The ratio of the transverse overlap area $S_{⊥}$ to $(N_{part}/2)^{2/3}$ versus $(N_{part}/2)^{2/3}$.

The ratio of the transverse overlap area $S_{⊥}$ to $(N_{part}/2)^{2/3}$ versus $(N_{part}/2)^{2/3}$.

The ratio of the charged pion multiplicity to the transverse overlap area $dN_{\pi}/dy/S_{⊥}$ versus $N_{coll}/N_{part}$. Errors shown are total errors, dominated by systematic uncertainties. The systematic uncertainties are correlated between Npart, Ncoll, and S⊥, and are largely canceled in the plotted ratio quantities.

The ratio of the charged pion multiplicity to the transverse overlap area $dN_{\pi}/dy/S_{⊥}$ versus $N_{coll}/N_{part}$. Errors shown are total errors, dominated by systematic uncertainties. The systematic uncertainties are correlated between Npart, Ncoll, and S⊥, and are largely canceled in the plotted ratio quantities.

The ratio of the charged pion multiplicity to the transverse overlap area $dN_{\pi}/dy/S_{⊥}$ versus $N_{coll}/N_{part}$. Errors shown are total errors, dominated by systematic uncertainties. The systematic uncertainties are correlated between Npart, Ncoll, and S⊥, and are largely canceled in the plotted ratio quantities.

Energy loss effect for $\pi^{+-}$ (a), $K^{+-}$ (b), and p and pbar (c) at mid-rapidity (|y| < 0.1) as a function of particle momentum magnitude in 200 GeV pp and 62.4 GeV central 0-5% Au+Au collisions. Only negative particles are shown; energy loss for particles and antiparticles are the same. Errors shown are statistical only. The pion energy loss is already corrected by the track reconstruction algorithm.

Energy loss effect for $\pi^{+-}$ (a), $K^{+-}$ (b), and p and pbar (c) at mid-rapidity (|y| < 0.1) as a function of particle momentum magnitude in 200 GeV pp and 62.4 GeV central 0-5% Au+Au collisions. Only negative particles are shown; energy loss for particles and antiparticles are the same. Errors shown are statistical only. The pion energy loss is already corrected by the track reconstruction algorithm.

Energy loss effect for $\pi^{+-}$ (a), $K^{+-}$ (b), and p and pbar (c) at mid-rapidity (|y| < 0.1) as a function of particle momentum magnitude in 200 GeV pp and 62.4 GeV central 0-5% Au+Au collisions. Only negative particles are shown; energy loss for particles and antiparticles are the same. Errors shown are statistical only. The pion energy loss is already corrected by the track reconstruction algorithm.

Vertex-finding efficiency (ǫvtx) and fake vertex rate (δfake) as a function of the number of good global tracks (a) and the number of good primary tracks (b). Errors shown or smaller than the point size are statistical only.

Vertex-finding efficiency (ǫvtx) and fake vertex rate (δfake) as a function of the number of good global tracks (a) and the number of good primary tracks (b). Errors shown or smaller than the point size are statistical only.

The $p_{T}$ dependent correction to particle spectra due to fake vertex events, $\epsilon_{fake}(p_{T})$, in 200 GeV minimum bias pp and d+Au collisions. Errors shown are statistical only.

The $p_{T}$ dependent correction to particle spectra due to fake vertex events, $\epsilon_{fake}(p_{T})$, in 200 GeV minimum bias pp and d+Au collisions. Errors shown are statistical only.

Efficiency (product of tracking efficiency and detector acceptance) of π−, K−, and p in pp (a) and d+Au collisions (b) at 200 GeV as a function of input MC p⊥. Errors shown are binomial errors. The curves are parameterizations to the efficiency data and are used for corrections in the analysis.

Efficiency (product of tracking efficiency and detector acceptance) of π−, K−, and p in pp (a) and d+Au collisions (b) at 200 GeV as a function of input MC p⊥. Errors shown are binomial errors. The curves are parameterizations to the efficiency data and are used for corrections in the analysis.

Efficiency (product of tracking efficiency and detector acceptance) of π−, K−, and p in 70-80% peripheral Au+Au (a) and 0-5% central Au+Au collisions (b) at 62.4 GeV as a function of input MC p⊥. Errors shown are binomial errors. The curves are parameterizations to the efficiency data and are used for corrections in the analysis.

Efficiency (product of tracking efficiency and detector acceptance) of π−, K−, and p in 70-80% peripheral Au+Au (a) and 0-5% central Au+Au collisions (b) at 62.4 GeV as a function of input MC p⊥. Errors shown are binomial errors. The curves are parameterizations to the efficiency data and are used for corrections in the analysis.

Efficiency (product of tracking efficiency and detector acceptance) of π−, K−, and p in 70-80% peripheral Au+Au (a) and 0-5% central Au+Au collisions (b) at 62.4 GeV as a function of input MC p⊥. Errors shown are binomial errors. The curves are parameterizations to the efficiency data and are used for corrections in the analysis.

Efficiency (product of tracking efficiency and detector acceptance) of π−, K−, and p in 70-80% peripheral Au+Au (a) and 0-5% central Au+Au collisions (b) at 62.4 GeV as a function of input MC p⊥. Errors shown are binomial errors. The curves are parameterizations to the efficiency data and are used for corrections in the analysis.

The dca distributions of protons and antiprotons for 0.40 < p⊥ < 0.45 GeV/c in 200 GeV minimum bias d+Au. Errors shown are statistical only.

The dca distributions of protons and antiprotons for 0.70 < p⊥ < 0.75 GeV/c in 200 GeV minimum bias d+Au. Errors shown are statistical only.

The dca distributions of protons and antiprotons for 0.40 < p⊥ < 0.45 GeV/c in 62.4 GeV 0-5% central Au+Au collisions. Errors shown are statistical only.

The dca distributions of protons and antiprotons for 0.70 < p⊥ < 0.75 GeV/c in 62.4 GeV 0-5% central Au+Au collisions. Errors shown are statistical only.

Pion background fraction from weak decays ($\Lambda, K_S^0$) and $\mu^{+-}$ contamination as a function of p⊥ in minimum bias d+Au collisions at 200 GeV. Errors shown are statistical only.

Pion background fraction from weak decays ($\Lambda, K_S^0$) and $\mu^{+-}$ contamination as a function of p⊥ in minimum bias d+Au collisions at 200 GeV. Errors shown are statistical only.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity identified antiproton spectra in 200 GeV minimum bias d+Au (a) and 62.4 GeV central Au+Au collisions (b) measured by dE/dx together with those by TOF [46, 47]. The dE/dx data are from |y| < 0.1 and the TOF data are from |y| < 0.5. The curves are various fits to the dE/dx data for extrapolation. The quadratic sum of statistical errors and point-to-point systematic errors are plotted, but are smaller than the point size.

Mid-rapidity (|y| < 0.1) identified particle spectra in d+Au collisions at 200 GeV. The p and pbar spectra are inclusive, including weak decay products. Spectra are plotted for three centrality bins and for minimum bias events. Spectra from top to bottom are for 0-20% scaled by 4, 20-40% scaled by 2, minimum bias not scaled, and 40-100% scaled by 1/2. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the minimum bias data; the normalizations of the curves are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified particle spectra in d+Au collisions at 200 GeV. The p and pbar spectra are inclusive, including weak decay products. Spectra are plotted for three centrality bins and for minimum bias events. Spectra from top to bottom are for 0-20% scaled by 4, 20-40% scaled by 2, minimum bias not scaled, and 40-100% scaled by 1/2. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the minimum bias data; the normalizations of the curves are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified particle spectra in d+Au collisions at 200 GeV. The p and pbar spectra are inclusive, including weak decay products. Spectra are plotted for three centrality bins and for minimum bias events. Spectra from top to bottom are for 0-20% scaled by 4, 20-40% scaled by 2, minimum bias not scaled, and 40-100% scaled by 1/2. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the minimum bias data; the normalizations of the curves are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified particle spectra in Au+Au collisions at 62.4 GeV. The p and p spectra are inclusive, including weak decay products. Spectra are plotted for nine centrality bins, from top to bottom, 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 60-70%, and 70-80%. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the spectra; the normalizations of the curves in (a,b) are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified particle spectra in Au+Au collisions at 62.4 GeV. The p and p spectra are inclusive, including weak decay products. Spectra are plotted for nine centrality bins, from top to bottom, 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 60-70%, and 70-80%. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the spectra; the normalizations of the curves in (a,b) are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified particle spectra in Au+Au collisions at 62.4 GeV. The p and p spectra are inclusive, including weak decay products. Spectra are plotted for nine centrality bins, from top to bottom, 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 60-70%, and 70-80%. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the spectra; the normalizations of the curves in (a,b) are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified particle spectra in Au+Au collisions at 62.4 GeV. The p and p spectra are inclusive, including weak decay products. Spectra are plotted for nine centrality bins, from top to bottom, 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 60-70%, and 70-80%. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but are smaller than the point size. The curves are the blast-wave model fits to the spectra; the normalizations of the curves in (a,b) are fixed by the corresponding negative particle spectra.

Mid-rapidity (|y| < 0.1) identified pion spectra in Au+Au collisions at 130 GeV. Spectra are plotted for eight centrality bins, from top to bottom, 0-6%, 6-11%, 11-18%, 18- 26%, 26-34%, 34-45%, 45-58%, and 58-85%. Errors plotted are statistical and point-to-point systematic errors added in quadrature, but they are smaller than the data point size. The curves are the Bose-Einstein fits to the spectra; the normalizations of the curves are fixed by the corresponding negative particle spectra.

Estimate of the product of the Bjorken energy density and the formation time ($\epsilon_{B_j}\tau$) as a function of centrality Npart. Errors shown are the quadratic sum of statistical and systematic uncertainties.

Estimate of the product of the Bjorken energy density and the formation time ($\epsilon_{B_j}\tau$) as a function of centrality Npart. Errors shown are the quadratic sum of statistical and systematic uncertainties.

Estimate of the product of the Bjorken energy density and the formation time ($\epsilon_{B_j}\tau$) as a function of centrality Npart. Errors shown are the quadratic sum of statistical and systematic uncertainties.

The K+/π+ and K−/π− ratios as a function of the collision energy in pp and central heavy-ion collisions.

The K+/K− ratio as a function of the collision energy in central heavy-ion collisions.

The K−/π− ratio as a function of the number of participants Npart in heavy-ion collisions at RHIC. Errors shown are the quadratic sum of statistical and systematic uncertainties.

The K−/π− ratio as a function of the number of participants Npart in heavy-ion collisions at RHIC. Errors shown are the quadratic sum of statistical and systematic uncertainties.

The K−/π− ratio as a function of the number of participants Npart in heavy-ion collisions at RHIC. Errors shown are the quadratic sum of statistical and systematic uncertainties.

The K−/π− ratio as a function of dNπ/dy/S⊥ in heavy-ion collisions at RHIC. Errors shown are the quadratic sum of statistical and systematic uncertainties.

The K−/π− ratio as a function of dNπ/dy/S⊥ in heavy-ion collisions at RHIC. Errors shown are the quadratic sum of statistical and systematic uncertainties.

The K−/π− ratio as a function of dNπ/dy/S⊥ in heavy-ion collisions at RHIC. Errors shown are the quadratic sum of statistical and systematic uncertainties.

Baryon chemical potential extracted for central heavy-ion collisions as a function of the collision energy. Errors shown are the total statistical and systematic errors.

The extracted chemical freeze-out temperatures for central heavy-ion collisions as a function of the collision energy. Errors shown are the total statistical and systematic errors.

The extracted chemical freeze-out temperatures for central heavy-ion collisions as a function of the collision energy. Errors shown are the total statistical and systematic errors.

The extracted kinetic freeze-out temperatures for central heavy-ion collisions as a function of the collision energy. Errors shown are the total statistical and systematic errors.

Average transverse radial flow velocity extracted from the blast-wave model for central heavy-ion collisions as a function of the collision energy. Errors shown are the total statistical and systematic errors.

Phase diagram plot of chemical freezeout temperature versus baryon chemical potential extracted from chemical equilibrium models. Errors shown are the total statistical and systematic errors.

Phase diagram plot of chemical freezeout temperature versus baryon chemical potential extracted from chemical equilibrium models. Errors shown are the total statistical and systematic errors.

Differential cross-sections obtained from the optical and MC Glauber calculations for Au+Au collisions at 200 GeV. Statistical errors are smaller than the point size.

Differential cross-sections obtained from the optical and MC Glauber calculations for Au+Au collisions at 200 GeV. Statistical errors are smaller than the point size.

(Color online) Dynamical net charge fluctuations, $\nu_{+−,dyn}$, of particles produced with pseudorapidity $|\eta|$ < 0.5 scaled by (a) the multiplicity, $dN_{ch}/d\eta$. The dashed line corresponds to charge conservation effect and the solid line to the prediction for a resonance gas, (b) the number of participants, and (c) the number of binary collisions.

(Color online) Dynamical net charge fluctuations, $\nu_{+−,dyn}$, of particles produced with pseudorapidity $|\eta|$ < 0.5 scaled by (a) the multiplicity, $dN_{ch}/d\eta$. The dashed line corresponds to charge conservation effect and the solid line to the prediction for a resonance gas, (b) the number of participants, and (c) the number of binary collisions.

(Color online) Dynamical fluctuations $\nu_{+−,dyn}$, normalized to their value for |$\eta$| < 1, as function of the integrated pseudorapidity range. (a) data for $Au + Au$ collisions at $\sqrt{s_{NN}}$ = 62.4, 200 GeV (0-5%) along with data for $Cu + Cu$ collisions at $\sqrt{s_{NN}}$ = 62.4, 200 GeV (0-10%), are compared to inclusive $p + p$ data at $\sqrt{s}$ = 200 GeV, and (b) data for $Au + Au$ collisions at $\sqrt{s_{NN}}$ = 62.4, 200 GeV (30-40%) along with data for $Cu + Cu$ collisions at $\sqrt{s_{NN}}$ = 62.4, 200 GeV (0-10%), are compared to inclusive $p + p$ collision data at $\sqrt{}s$ = 200 GeV.

(Color online) Dynamical fluctuations $\nu_{+−,dyn}$, normalized to their value for |$\eta$| < 1, as function of the integrated pseudorapidity range. (a) data for $Au + Au$ collisions at $\sqrt{s_{NN}}$ = 62.4, 200 GeV (0-5%) along with data for $Cu + Cu$ collisions at $\sqrt{s_{NN}}$ = 62.4, 200 GeV (0-10%), are compared to inclusive $p + p$ data at $\sqrt{s}$ = 200 GeV, and (b) data for $Au + Au$ collisions at $\sqrt{s_{NN}}$ = 62.4, 200 GeV (30-40%) along with data for $Cu + Cu$ collisions at $\sqrt{s_{NN}}$ = 62.4, 200 GeV (0-10%), are compared to inclusive $p + p$ collision data at $\sqrt{}s$ = 200 GeV.

(Color online) Dynamical fluctuations $\nu_{+−,dyn}$, as a function of the integrated azimuthal range $\phi$ for selected collision centralities for (a) $Au + Au$ collisions at $\sqrt{s_{NN}}$ = 200 GeV, and (b) $Cu + Cu$ collisions at $\sqrt{s_{NN}}$ = 200 GeV.

(Color online) Dynamical fluctuations $\nu_{+−,dyn}$, as a function of the integrated azimuthal range $\phi$ for selected collision centralities for (a) $Au + Au$ collisions at $\sqrt{s_{NN}}$ = 200 GeV, and (b) $Cu + Cu$ collisions at $\sqrt{s_{NN}}$ = 200 GeV.

(Color online) Dynamical fluctuations $\nu_{+−,dyn}$, as a function of the integrated azimuthal range $\phi$ for similar number of participating nucleons for $Au + Au$ and $Cu + Cu$ $\sqrt{s_{NN}}$ = 200 GeV.

Charged particle $v_1(\eta)$ for 40-80 % centrality in Au+Au collisions at 200 GeV.

Charged particle $v_1(p_t)$ in 200 GeV Au+Au for three centralities and two pseudorapidity ranges.

Charged particle $v_1(\eta)$ for mid-central (30-60%) Au+Au and Cu+Cu collisions at 200 GeV and 62.4 GeV.

Charged particle $v_1$ versus centrality for Au+Au and Cu+Cu collisions at 200 GeV and 62.4 GeV in two pseudorapidity ranges.

Charged particle $v_1$ versus $\eta-y_{beam}$, for 30 − 60% Au+Au collisions at 200 GeV.

Charged particle $v_1$ versus $\eta-y_{beam}$, for 30 − 60% Cu+Cu collisions at 200 GeV.

Charged particle $v_1$ versus $\eta-y_{beam}$, for 30 − 60% Au+Au collisions at 62.4 GeV.

Charged particle $v_1$ versus $\eta-y_{beam}$, for 30 − 60% Cu+Cu collisions at 62.4 GeV.

FIG. 1: (a) Raw two-particle correlation signal $Y_2$ (red), background $aB_{inc}F_2$ (solid histogram), and background systematic uncertainty from a (dashed histograms). (b) Background-subtracted two-particle correlation $\hat{Y}_2$ (red), and systematic uncertainties due to a (dashed histograms) and flow (blue histograms). (c) Raw three-particle correlation $Y_3$. (d) $ba^2Y_{inc}^2$ . (e) Sum of trig-corr-bkgd and trigger flow. Data are from 12% central Au+Au collisions. Statistical errors in (a,b) are smaller than the point size.

FIG. 1: (a) Raw two-particle correlation signal $Y_2$ (red), background $aB_{inc}F_2$ (solid histogram), and background systematic uncertainty from a (dashed histograms). (b) Background-subtracted two-particle correlation $\hat{Y}_2$ (red), and systematic uncertainties due to a (dashed histograms) and flow (blue histograms). (c) Raw three-particle correlation $Y_3$. (d) $ba^2Y_{inc}^2$ . (e) Sum of trig-corr-bkgd and trigger flow. Data are from 12% central Au+Au collisions. Statistical errors in (a,b) are smaller than the point size.

FIG. 2: Background subtracted three-particle correlations, $\hat{Y}_3$, for (a) pp, (b) d+Au, (c) 80-50% Au+Au, (d) 50-30% Au+Au, (e) 30-10% Au+Au, and (f) 12% central Au+Au. Statistical errors per bin are approximately $\pm$0.012 in (a) and $\pm$0.006 in (b), both at ($\pi$, $\pi$), and are $\pm$0.022, $\pm$0.049, $\pm$0.099 and $\pm$0.077 from (c) to (f), similar for all bins.

FIG. 2: Background subtracted three-particle correlations, $\hat{Y}_3$, for (a) pp, (b) d+Au, (c) 80-50% Au+Au, (d) 50-30% Au+Au, (e) 30-10% Au+Au, and (f) 12% central Au+Au. Statistical errors per bin are approximately $\pm$0.012 in (a) and $\pm$0.006 in (b), both at ($\pi$, $\pi$), and are $\pm$0.022, $\pm$0.049, $\pm$0.099 and $\pm$0.077 from (c) to (f), similar for all bins.

FIG. 2: Background subtracted three-particle correlations, $\hat{Y}_3$, for (a) pp, (b) d+Au, (c) 80-50% Au+Au, (d) 50-30% Au+Au, (e) 30-10% Au+Au, and (f) 12% central Au+Au. Statistical errors per bin are approximately $\pm$0.012 in (a) and $\pm$0.006 in (b), both at ($\pi$, $\pi$), and are $\pm$0.022, $\pm$0.049, $\pm$0.099 and $\pm$0.077 from (c) to (f), similar for all bins.

FIG. 2: Background subtracted three-particle correlations, $\hat{Y}_3$, for (a) pp, (b) d+Au, (c) 80-50% Au+Au, (d) 50-30% Au+Au, (e) 30-10% Au+Au, and (f) 12% central Au+Au. Statistical errors per bin are approximately $\pm$0.012 in (a) and $\pm$0.006 in (b), both at ($\pi$, $\pi$), and are $\pm$0.022, $\pm$0.049, $\pm$0.099 and $\pm$0.077 from (c) to (f), similar for all bins.

FIG. 2: Background subtracted three-particle correlations, $\hat{Y}_3$, for (a) pp, (b) d+Au, (c) 80-50% Au+Au, (d) 50-30% Au+Au, (e) 30-10% Au+Au, and (f) 12% central Au+Au. Statistical errors per bin are approximately $\pm$0.012 in (a) and $\pm$0.006 in (b), both at ($\pi$, $\pi$), and are $\pm$0.022, $\pm$0.049, $\pm$0.099 and $\pm$0.077 from (c) to (f), similar for all bins.

FIG. 2: Background subtracted three-particle correlations, $\hat{Y}_3$, for (a) pp, (b) d+Au, (c) 80-50% Au+Au, (d) 50-30% Au+Au, (e) 30-10% Au+Au, and (f) 12% central Au+Au. Statistical errors per bin are approximately $\pm$0.012 in (a) and $\pm$0.006 in (b), both at ($\pi$, $\pi$), and are $\pm$0.022, $\pm$0.049, $\pm$0.099 and $\pm$0.077 from (c) to (f), similar for all bins.

Projections of away-side three-particle correlations along the diagonal $\Sigma$ within 0 < $\Delta$ < 0.35 (squares) and along the off-diagonal $\Delta$ within |$\Sigma$| < 0.35 (points) in (a) d+Au and (b) 12% central Au+Au collisions. The shaded areas indicate systematic uncertainties on the off-diagonal projections. The histogram in (a) is the near-side off-diagonal projection. The histogram in (b) is the away-side off-diagonal projection of our result with a = b = 1. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Projections of away-side three-particle correlations along the diagonal $\Sigma$ within 0 < $\Delta$ < 0.35 (squares) and along the off-diagonal $\Delta$ within |$\Sigma$| < 0.35 (points) in (a) d+Au and (b) 12% central Au+Au collisions. The shaded areas indicate systematic uncertainties on the off-diagonal projections. The histogram in (a) is the near-side off-diagonal projection. The histogram in (b) is the away-side off-diagonal projection of our result with a = b = 1. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Three-particle correlation strength per radian$^2$ versus ($N_{part}$/2)$^{1/3}$ where $N_{part}$ is the number of participants. Some of the data points have been displaced in $N_{part}$ for clarity.

Analyzing powers A_N(-x_F) in x_F bins at < eta > =3.3 and x_F < 0.

Analyzing powers versus pi^0 transverse momentum (p_T) for events with scaled pi^0 longitudinal momentum x_F > 0.4.

Analyzing powers versus pi^0 transverse momentum (p_T) for events with scaled pi^0 longitudinal momentum x_F < -0.4.

Analyzing powers versus pi^0 transverse momentum (p_T) in 0.25 < x_F < 0.30 bin.

Analyzing powers versus pi^0 transverse momentum (p_T) in 0.30 < x_F < 0.35 bin.

Analyzing powers versus pi^0 transverse momentum (p_T) in 0.35 < x_F < 0.40 bin.

Analyzing powers versus pi^0 transverse momentum (p_T) in 0.40 < x_F < 0.47 bin.

Analyzing powers versus pi^0 transverse momentum (p_T) in 0.47 < x_F < 0.56 bin.

Analyzing powers versus pi^0 transverse momentum (p_T) in 0.56 < x_F < 0.80 bin.

$v_{2}$ as a function of $p_{T}$ for charged hadrons with $|\eta| < 1.0$ in 10–40$%$ $Au+Au$ collisions, at $\sqrt{s_{NN}} = 200$ GeV, from the Lee-Yang Zero method (solid circles) with Sum Generating Function.

$v_{2}$ as a function of $p_{T}$ for charged hadrons with $|\eta| < 1.0$ in 0–80$%$ $Au+Au$ collisions, at $\sqrt{s_{NN}} = 200$ GeV, from the Event Plane method.

$p_{T}$ integrated charged hadron $v_{2}$ in the TPC as a function of geometrical cross section. Shown are the Event Plane method ($v_{2}${EP}) (open circles). All data are from $Au+Au$ collisions, at $\sqrt{s_{NN}} = 200$ GeV. For the TPC, $|\eta|< 1.0$ was used.

$p_{T}$ integrated charged hadron $v_{2}$ in the TPC as a function of geometrical cross section. Shown are the 4 Cumulant method ($v_{2}${4}) (solid squares). All data are from $Au+Au$ collisions, at $\sqrt{s_{NN}} = 200$ GeV. For the TPC, $|\eta|< 1.0$ was used.

$p_{T}$ integrated charged hadron $v_{2}$ in the TPC as a function of geometrical cross section. Shown are the Lee-Yang Zero method with Sum Generating Function ($v_{2}$) (solid circles). All data are from $Au+Au$ collisions, at $\sqrt{s_{NN}} = 200$ GeV. For the TPC, $|\eta|< 1.3$ was used.

$p_{T}$ integrated charged hadron $v_{2}$ in the TPC as a function of geometrical cross section. Shown are the Lee-Yang Zero method with Product Generating Function ($v_{2}$) (open stars). All data are from $Au+Au$ collisions, at $\sqrt{s_{NN}} = 200$ GeV. For the TPC, $|\eta|< 1.3$ was used.

$v_{2}$ as a function of $p_{T}$ for $10–40\%$ centrality using $\eta$-subevent method (open crosses), are shown in (a) for $\Lambda$. All data are from $\sqrt{s_{NN}} = 200$ GeV $Au+Au$ collisions.

$v_{2}$ as a function of $p_{T}$ for $10–40\%$ centrality using Event Plane method (open circles), are shown in (a) for $\Lambda$. All data are from $\sqrt{s_{NN}} = 200$ GeV $Au+Au$ collisions.

$v_{2}$ as a function of $p_{T}$ for $10–40\%$ centrality using $\eta$-subevent method (open crosses), are shown in (b) for $K_{S}^{0}$. All data are from $\sqrt{s_{NN}} = 200$ GeV $Au+Au$ collisions.

$v_{2}$ as a function of $p_{T}$ for $10–40\%$ centrality using Event Plane method (open circles), are shown in (b) for $K_{S}^{0}$. All data are from $\sqrt{s_{NN}} = 200$ GeV $Au+Au$ collisions.

$v_{2}$ as a function of $p_{T}$ for $10–40\%$ centrality using Lee-Yang Zero method with Sum Generating Function (solid circles), are shown in (a) and (b) for $\Lambda$ and $K_{S}^{0}$ respectively. All data are from $\sqrt{s_{NN}} = 200$ GeV $Au+Au$ collisions.

$v_{2}$ of $K_{S}^{0}$ (open circles) as a function of $p_{T}$ for (a) $0–80\%$, (b) $40–80\%$, (c) $10–40\%$ and (d) $0–10\%$ centralities in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Lambda+\bar{\Lambda}$ (open squares) as a function of $p_{T}$ for (a) $0–80\%$, (b) $40–80\%$, (c) $10–40\%$ and (d) $0–10\%$ centralities in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Xi^{-}+\bar{\Xi}^{+}$ (filled triangles) as a function of $p_{T}$ for (a) $0–80\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Xi^{-}+\bar{\Xi}^{+}$ (filled triangles) as a function of $p_{T}$ for (b) $40–80\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Xi^{-}+\bar{\Xi}^{+}$ (filled triangles) as a function of $p_{T}$ for (c) $10–40\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Xi^{-}+\bar{\Xi}^{+}$ (filled triangles) as a function of $p_{T}$ for (d) $0–10\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Omega^{-}+\bar{\Omega}^{+}$ (filled circles) as a function of $p_{T}$ for (a) $0–80\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $\Omega^{-}+\bar{\Omega}^{+}$ (filled circles) as a function of $p_{T}$ for (d) $0–10\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

$v_{2}$ of $p+\bar{p}$ (filled squares) and $\pi^{+}+\pi^{-}$ (filled stars) as a function of $p_{T}$ for (a) $0–80\%$ centrality in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV.

Fig12d inset $\Lambda$ (open squares) expansion at low $p_{T}$

Fig12d inset for $K_{S}^{0}$ (filled circles) expansion at low $p_{T}$

Fig13 Centrality dependence of $v_{2}$ versus number of participants $N_{part}$ for charged hadrons (crosses) in $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV. For identified particle's value, $v_{2}$ in Tab. II scaled by participant eccentricity in Tab. I versus number of participant in Tab. I.

Fig14 $p_{T}$ dependence of $v_{2}$ of $K_{S}^{0}$ at 200 GeV.

Fig14 $p_{T}$ dependence of $v_{2}$ of $\Lambda$ at 200 GeV.

The invariant mass distribution for the coherently produced $\rho^{0}$ candidates obtained from the topology sample with the cut on the $\rho^{0}$ transverse momentum $p_{T}$ < 150 MeV/c. The hatched area is the contribution from the combinatorial background. The solid line corresponds to Eq. 3 which encompasses the Breit-Wigner (dashed), the mass independent contribution from the direct $\pi^{+}\pi^{-}$ production (dash-dotted), and the interference term(dotted).

The ratio |B/A| as a function of $y_{\rho^{0}}$ for the minimum bias data, obtained by fitting Eq.3 to the invariant mass distributions in bins of $y_{\rho^{0}}$.

Coherent $\rho^{0}$ production cross-section for the minimum bias data set as a function of $y_{\rho^{0}}$ (black dots) overlaid by the $d\sigma/dy$ distribution predicted by the KN model [22] (solid line).

$\rho^{0}$ production cross-section as a function of the momentum transfer squared $t$, together with the fit of Eq. 5. The fit parameters are shown in Table I.

Comparison of theoretical predictions to the measured differential cross-section for coherent $\rho^{0}$ production. The statistical errors are shown by the solid vertical line at each data point. The sum of the statistical and systematic error bars is shown by the grey band.

Projections of the two dimensional efficiency corrected $\Phi_{h}$ vs $cos(\Theta_{h})$ distributions obtained with the minimum bias data set. The solid line shows the result of the two-dimensional fit to the data with Eq. 6 and the coefficients given in Tab. III.

Projections of the two dimensional efficiency corrected $\Phi_{h}$ vs $cos(\Theta_{h})$ distributions obtained with the minimum bias data set. The solid line shows the result of the two-dimensional fit to the data with Eq. 6 and the coefficients given in Tab. III.

$R_{AA}$ as a function of $p_{T}$ from $0-5\%$ and $60-80\%$ central Au+Au events for $\Lambda$. Stat. and syst. errors added in quadrature.The band at unity shows the systematic uncertainty on $⟨N_{bin}⟩$. The dashed line below unity shows the expected value of RAA should the yields scale with $⟨N_{part}⟩$, and the band around it shows the systematic uncertainty on ⟨Npart ⟩. $N_{bin}$ Scaling 1 $\pm$ 0.16 and $N_{part}$ Scaling 0.168 $\pm$ 0.02 .

$R_{AA}$ as a function of $p_{T}$ from $0-5\%$ and $60-80\%$ central Au+Au events for $\Xi^{-} + \bar{\Xi}^{+}$. Stat. and syst. errors added in quadrature. The band at unity shows the systematic uncertainty on $⟨N_{bin}⟩$. The dashed line below unity shows the expected value of RAA should the yields scale with $⟨N_{part}⟩$, and the band around it shows the systematic uncertainty on ⟨Npart ⟩. $N_{bin}$ Scaling 1 $\pm$ 0.16 and $N_{part}$ Scaling 0.168 $\pm$ 0.02 .

$R_{AA}$ as a function of $p_{T}$ from $0-5\%$ and $60-80\%$ central Au+Au events for inclusive $p+\bar{p}$. Stat. and syst. errors added in quadrature. The band at unity shows the systematic uncertainty on $⟨N_{bin}⟩$. The dashed line below unity shows the expected value of RAA should the yields scale with $⟨N_{part}⟩$, and the band around it shows the systematic uncertainty on ⟨Npart ⟩. $N_{bin}$ Scaling 1 $\pm$ 0.16 and $N_{part}$ Scaling 0.168 $\pm$ 0.02 .

$R_{AA}$ from HIJING as a function of $p_{T}$ from $0-5\%$ central Au+Au events for inclusive $p+\bar{p}$.

$R_{AA}$ from HIJING as a function of $p_{T}$ from $0-5\%$ central Au+Au events for $\Lambda$.

$R_{AA}$ from HIJING as a function of $p_{T}$ from $0-5\%$ central Au+Au events for $\Xi^{-}+\bar{\Xi}^{+}$.

$R_{AA}$ from EPOS as a function of $p_{T}$ from $0-5\%$ central Au+Au events for inclusive $p+\bar{p}$.

$R_{AA}$ from EPOS as a function of $p_{T}$ from $0-5\%$ central Au+Au events for $\Lambda$.

$R_{AA}$ from EPOS as a function of $p_{T}$ from $0-5\%$ central Au+Au events for $\Xi^{-}+\bar{\Xi}^{+}$.

$R_{AA}$ from HIJING as a function of $p_{T}$ from $60-80\%$ central Au+Au events for inclusive $p+\bar{p}$.

$R_{AA}$ from HIJING as a function of $p_{T}$ from $60-80\%$ central Au+Au events for $\Lambda$.

$R_{AA}$ from HIJING as a function of $p_{T}$ from $60-80\%$ central Au+Au events for $\Xi^{-}+\bar{\Xi}^{+}$.

$R_{AA}$ from EPOS as a function of $p_{T}$ from $60-80\%$ central Au+Au events for inclusive $p+\bar{p}$.

$R_{AA}$ from EPOS as a function of $p_{T}$ from $60-80\%$ central Au+Au events for $\Lambda$.

$R_{AA}$ from EPOS as a function of $p_{T}$ from $60-80\%$ central Au+Au events for $\Xi^{-}+\bar{\Xi}^{+}$.