$v_4$ from three particle cumulant as a function of $p_T$ for the minimum bias Au+Au collisions. Background from secondary particles is expected to be less than 15%. Non-flow systematic uncertainties are approximately 20%. Fluctuations in initial geometry can lead to an effect of about a factor of 1.2 to 1.5.

$v_6$ with respect to the second harmonic event plane as a function of $p_T$ for the minimum bias Au+Au collisions. Background from secondary particles is expected to be less than 15%. Non-flow systematic uncertainties are approximately 20%. Fluctuations in initial geometry can lead to an effect of about a factor of 1.2 to 1.5.

$p_T$- and $\eta$-integrated $v_2$, $v_4$, $v_6$ as a function of centrality. Background from secondary particles is expected to be less than 15%. Non-flow systematic uncertainties are approximately 20%. Fluctuations in initial geometry can lead to an effect of about a factor of 1.2 to 1.5.

$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}^{+}$.

FIG. 4. Ridge yield Eq. (9) in $|\Delta \eta|<1.7$ and $2 \mathrm{GeV/c}<p^{assoc}_{t}$ < $p^{trig}_{t}$ as a function of $p_{t}^{t r i g}$. Solid lines are the systematic uncertainty due to $v_{2}$.

FIG. 5. Projection of $d N /\left.d \Delta \phi\right|_{a, b}$ for $0.7<|\Delta \eta|<1.4$ [Eq. (3)] in two trigger $p_{t}$ windows for $2<p_{t}^{\text {assoc }}<4 \mathrm{GeV} / \mathrm{c} .$ No background subtraction has been applied; note the suppressed zero on the vertical scale. The shaded band shows the fit of the function $k_{1}+k_{2}$. $\cos (2 \Delta \phi)$ to the recoil region $2<|\Delta \phi|<\pi$ for $4<p_{t}^{\text {trig }}<6 \mathrm{GeV} / c$. The width of the band indicates the fitting error. The solid curve represents the background estimate using the ZYAM normalization for $4<p_{t}^{\text {trig }}<6 \mathrm{GeV} / c .$ Systematic uncertainties are indicated by the light shaded band.

FIG. 6. (Color online) Differential $p_{t}$ spectrum for associated particles in central $\mathrm{Au}+$ Au collisions, with $4<p_{t}^{\text {trig }}<6$ and $6<p_{t}^{\text {trig }}<10 \mathrm{GeV} / c$ . The dash-dotted line is the inclusive hadron spectrum from central $\mathrm{Au}+$ Au collisions $[1] .$ Left panel: ridge spectrum; shaded bands show systematic uncertainty. Right panel: jet-like spectrum, also compared to $d+$ Au reference measurements. The lines in both panels show exponential fits to the data (see Table III). Data are offset horizontally for clarity.

FIG. 6. (Color online) Differential $p_{t}$ spectrum for associated particles in central $\mathrm{Au}+$ Au collisions, with $4<p_{t}^{\text {trig }}<6$ and $6<p_{t}^{\text {trig }}<10 \mathrm{GeV} / c$ . The dash-dotted line is the inclusive hadron spectrum from central $\mathrm{Au}+$ Au collisions $[1] .$ Left panel: ridge spectrum; shaded bands show systematic uncertainty. Right panel: jet-like spectrum, also compared to $d+$ Au reference measurements. The lines in both panels show exponential fits to the data (see Table III). Data are offset horizontally for clarity.

Non photonic electron yield in minimum bias D+AU collisions versus $p_{T}$.

Non photonic electron yield in AU+AU collisions versus PT, for a centrality range of 40-80% To obtain a differential cross-section in mb/(GeV2), multiply listed data by 30 Note that, in addition to the statistical and systematical errors, there is a normalization error on the value, given in the second column.

Non photonic electron yield in Au+Au collisions versus $p_{T}$, for a centrality range of 40-80%.

Non photonic electron yield in AU+AU collisions versus PT, for a centrality range of 10-40% To obtain a differential cross-section in mb/(GeV2), multiply listed data by 30 Note that, in addition to the statistical and systematical errors, there is a normalization error on the value, given in the second column.

Non photonic electron yield in Au+Au collisions versus $p_{T}$, for a centrality range of 10-40%.

Non photonic electron yield in AU+AU collisions versus PT, for a centrality range of 0-5% To obtain a differential cross-section in mb/(GeV2), multiply listed data by 30 Note that, in addition to the statistical and systematical errors, there is a normalization error on the value, given in the second column.

Non photonic electron yield in Au+Au collisions versus $p_{T}$, for a centrality range of 0-5%.

Nuclear modification factor for non-photonic electrons in minimum bias D+AU reactions Note that, in addition to the statistical and systematical errors, there is a normalization error on the value, given in the second column.

Nuclear modification factor for non-photonic electrons in minimum bias d+Au reactions.

Nuclear modification factor for non-photonic electrons in AU+AU reactions in the centrality range of 0-5% Note that, in addition to the statistical and systematical errors, there is a normalization error on the value, given in the second column.

Nuclear modification factor for non-photonic electrons in Au+Au reactions in the centrality range of 0-5%.

$N_{part}$ dependence of the near-side associated-particle yield.

$N_{part}$ dependence of the near-side associated-particle yield.

$N_{part}$ dependence of the near-side associated-particle yield.

$N_{part}$ dependence of the near-side associated-particle yield.

$N_{part}$ dependence of the near-side associated-particle yield.

$N_{part}$ dependence of the near-side associated-particle yield.

$N_{part}$ dependence of the near-side associated-particle yield.

$N_{part}$ dependence of the near-side associated-particle yield.

$N_{part}$ dependence of the near-side associated-particle yield.

$N_{part}$ dependence of the near-side associated-particle yield.

$N_{part}$ dependence of the near-side associated-particle yield.

$N_{part}$ dependence of the near-side associated-particle yield.

$N_{part}$ dependence of the near-side associated-particle yield.

$N_{part}$ dependence of the near-side associated-particle yield.

$N_{part}$ dependence of the near-side associated-particle yield.

$N_{part}$ dependence of the near-side associated-particle yield.

$N_{part}$ dependence of the near-side associated-particle yield.

$N_{part}$ dependence of the near-side associated-particle yield.

$\Delta\phi$ distribution used to extract the away-side yield. The triangular pair acceptance correction in $\Delta\eta$ is not applied.

$N_{part}$ dependence of the away-side associated-particle yield.

$N_{part}$ dependence of the away-side associated-particle yield.

$N_{part}$ dependence of the away-side associated-particle yield.

$N_{part}$ dependence of the away-side associated-particle yield.

Away-side associated particle distribution.

$I_{AA}$.

$I_{AA}$.

$I_{AA}$.

$I_{AA}$.

Background-subtracted charge-independent (AAT ) correlated hadron pair density in 40-80% Au+Au collisions for 3<p(t)⊥<10 GeV/cand 1<p(a)⊥<3 GeV/c. The results are for near-side correlated hadrons within |∆φ1,2|<0.7, and corrected for the 3-particle ∆η-∆η acceptance. Statistical errors at (∆η1,∆η2)∼(0,0)are approximately 0.058 for Au+Au respectively..

Background-subtracted charge-independent (AAT ) correlated hadron pair density in 0-12% Au+Au collisions for 3<p(t)⊥<10 GeV/cand 1<p(a)⊥<3 GeV/c. The results are for near-side correlated hadrons within |∆φ1,2|<0.7, and corrected for the 3-particle ∆η-∆η acceptance. Statistical errors at (∆η1,∆η2)∼(0,0)are approximately 0.084 for Au+Au respectively..

The average correlated hadron pair density per trigger particle as a function of R for all charges in minimum 0-12% Au+Au collision. The average〈ˆP〉for AAT as a function of R = $\sqrt{∆η21+ ∆η22}$. The average density is peaked at R∼0 and decreases with R for all systems. In 0-12% Au+Au collisions, the average denstiy drops more slowly andbecomes approximately constant aboveR>1, indicatingthe presence of the ridge.

The average correlated hadron pair density per trigger particle as a function of R for all charges in minimum 40-80% Au+Au collision. The average〈ˆP〉for AAT as a function of R = $\sqrt{∆η21+ ∆η22}$. The average density is peaked at R∼0 and decreases with R for all systems. For 40-80% Au+Au collisions the average density at R>1 is consistent with zero, indicating no ridge contribution

The average correlated hadron pair density per trigger particle as a function of R for all charges in minimum d+Au collision. The average〈ˆP〉for AAT as a function of R = $\sqrt{∆η21+ ∆η22}$. The average density is peaked at R∼0 and decreases with R for all systems. For d+Au collisions the average density at R>1 is consistent with zero, indicating no ridge contribution

Radial dependence of average 3-particle correlation pair density for 0-12% Au+Au (like-sign triplets) data. The results are shown in Fig. 3(b) Indeed, no jet-like component is apparent in A^±A^±T^±.The A^±A^±T^∓ result contains both jet-like and ridge components.

Radial dependence of average 3-particle correlation pair density for 40-80% Au+Au (like-sign triplets) data. The contribution from other charge combinations, namely A^±A^∓ T^±, are simply the difference between AAT in Fig.3(a) and (A^±A^±T^± + A^±A^± T^∓) in Fig. 3(b). We found this to be equal to twice the A^±A^± T^∓contribution within errors.

Radial dependence of average 3-particle correlation pair density for d+Au (AAT) data. We expect the ridge contributions in the correlated pair density to be the same in all charge combinations. We verified this for large ∆η correlatedpair densities within our current statistics, as can be seen from Fig. 3(b). Therefore the total ridge particle pair density (ˆPrr) can be obtained as four times A^±A^±T^±. The remaining jet-like signal, the sum of jet-like correlated particle pairs (ˆPjj) and cross pairs of a jet-like and a ridge particle (ˆPjr), can then be obtained by subtracting the total ridge from AAT.

for same-sign associated particles (A^\pm A^\pm T^\pm and A^\pm A^\pm T^\mp) in 0-12% Au+Au collisions Systematic uncertainties are shown in the shaded boxes due to background normalization and in the solid curves due to flow.

The R dependence of the average〈ˆPrr〉and〈ˆPjj〉+〈ˆPjr〉in 0-12% Au+Au collisions.The ridge pair density is consistent with a constant 60.14±0.02 (χ2/ndf=5.8/7). Gaussian fits indicate a best fit value σ = 2.1 (χ2/ndf=4.8/6, solid curve) and σ >1.4 (dashed curve) with 84% confidence level. On the otherhand, the jet-like component is narrow with a Gaussianσ=0.34+0.13−0.09(χ2/ndf=0.8/6, dominated by statistical errors), comparing well to those from the correlated single hadron density.

The average correlated hadron pair density per trigger particle in 0-12% Au+Au collisions for the jet-like andridge components as a function of R. The solid curves are Gaussian fits. The dashed curve is a Gaussian fit with a fixed σ=1.4 to the ridge data. The systematic uncertainities on the ridge data are shown in shaded boxes due to background normalization and in open boxes due to flow.

The R dependence of the average in 0-12% Au+Au collisions.for the ridge as a function of ξ within R<1.4. The solid curves are Gaussian fits. To investigate possible structures in the ridge, we show in Fig. 4(b) the average ridge particle pair density as a function of ξ = arctan(∆η2/∆η1) within R<1.4. The data are consistent with a uniform distribution in ξ(χ2/ndf=1.7/7). This suggests that the ridge particles are uncorrelated in ∆η not only with the trigger particle but also between themselves. In other words,the ridge appears to be uniform in ∆η event-by-event.

FIG. 2: a) left side: The correlation data for the ratio of the histograms of same-event-pairs to mixed-event-pairs for unlike-sign charged pairs, shown in a two-dimensional (2-D) perspective plot $\Delta\phi$ - $\Delta\eta$. The plot was normalized to a mean of 1. b) right side: The similar correlation data for like-sign charge pairs.

FIG. 3: a) left side: Folded correlation data for unlike-sign charge pairs on $\Delta\phi$ and $\Delta\eta$ based on demonstrated symmetries in the data. This increases the statistics in a typical bin by a factor of four. Henceforth we will be dealing with folded data only. b) right side: Similar folded data for like-sign charge pairs.

FIG. 3: a) left side: Folded correlation data for unlike-sign charge pairs on $\Delta\phi$ and $\Delta\eta$ based on demonstrated symmetries in the data. This increases the statistics in a typical bin by a factor of four. Henceforth we will be dealing with folded data only. b) right side: Similar folded data for like-sign charge pairs.

FIG. 4: “(Color online)”a) left side: The 2 dimensional (2-D) approximately gaussian signal shape equation (4) from the fit to the normalized and folded unlike-sign charge pairs signal data. b) right side: The unlike-sign charge pairs signal data corresponding to the adjoining signal function on the left side.

FIG. 4: “(Color online)”a) left side: The 2 dimensional (2-D) approximately gaussian signal shape equation (4) from the fit to the normalized and folded unlike-sign charge pairs signal data. b) right side: The unlike-sign charge pairs signal data corresponding to the adjoining signal function on the left side.

FIG. 5 :“(Color online)”a) left side: A perspective plot of the fit to the 2-D like-sign charge pairs signal shape of equation (5). b) right side: Like-sign charge pairs signal data corresponding to the adjacent signal fit on the left side.

FIG. 5: “(Color online)”a) left side: A perspective plot of the fit to the 2-D like-sign charge pairs signal shape of equation (5). b) right side: Like-sign charge pairs signal data corresponding to the adjacent signal fit on the left side.

FIG. 6: “(Color online)”a) left side: The charge dependent (CD) signal shape is a 2 dimensional (2-D) approximate gaussian which is symmetrical. The average gaussian width is 27.5$^{\circ}$ $\pm$ 3.2$^{\circ}$ (0.48 $\pm$ 0.056 rad) both in $\Delta\phi$ and $\Delta\eta$. Thus we are observing a greater probability for unlike-sign charge pairs of particles than for like-sign charge pairs emitted randomly on 2-D $\eta\phi$ directions. b) right side: The CD signal data corresponding to the adjacent signal fit on the left.

FIG. 6: “(Color online)”a) left side: The charge dependent (CD) signal shape is a 2 dimensional (2-D) approximate gaussian which is symmetrical. The average gaussian width is 27.5$^{\circ}$ $\pm$ 3.2$^{\circ}$ (0.48 $\pm$ 0.056 rad) both in $\Delta\phi$ and $\Delta\eta$. Thus we are observing a greater probability for unlike-sign charge pairs of particles than for like-sign charge pairs emitted randomly on 2-D $\eta\phi$ directions. b) right side: The CD signal data corresponding to the adjacent signal fit on the left.

FIG. 7: “(Color online)”a) left side: The fit to the charge independent (CI) signal shape which is the sum of the fits of like-sign plus unlike-sign charge pairs signals. The 2-D gaussian equivalent RMS signal has a $\Delta\phi$ width of about 32$^{\circ}$ (0.52 rad), and $\Delta\eta$ width of about 66$^{\circ}$ (1.15 rad) which is about double the $\Delta\phi$ width. b) right side: The CI signal data corresponding to the adjacent fit on the left.

FIG. 7: “(Color online)”a) left side: The fit to the charge independent (CI) signal shape which is the sum of the fits of like-sign plus unlike-sign charge pairs signals. The 2-D gaussian equivalent RMS signal has a $\Delta\phi$ width of about 32$^{\circ}$ (0.52 rad), and $\Delta\eta$ width of about 66$^{\circ}$ (1.15 rad) which is about double the $\Delta\phi$ width. b) right side: The CI signal data corresponding to the adjacent fit on the left.

FIG. 8: a) left side: The charge independent (CI) signal data with the lower range of elliptic flow amplitude and the best $\chi^2$ (the same as Figure 7b) plotted as a 2-D perspective plot on $\Delta\phi$ vs. $\Delta\eta$. b) right side: The charge independent (CI) signal data with our maximum elliptic flow amplitude used in our analysis ($\chi^2$ was worse by 1$\sigma$).

FIG. 8: a) left side: The charge independent (CI) signal data with the lower range of elliptic flow amplitude and the best $\chi^2$ (the same as Figure 7b) plotted as a 2-D perspective plot on $\Delta\phi$ vs. $\Delta\eta$. b) right side: The charge independent (CI) signal data with our maximum elliptic flow amplitude used in our analysis ($\chi^2$ was worse by 1$\sigma$).

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.