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.

The spin transfer $D_{LL}$ to (a) $\Lambda$ and (b) $\bar{\Lambda}$ hyperons produced at positive pseudorapidity with respect to the polarized proton beam from $MB$, $JP$, and $HT$ data versus hyperon transverse momenta $p_{T}$. The sizes of the statistical and systematic uncertainties are indicated by the vertical bars and bands, respectively. For clarity, the HT data points have been shifted slightly in $p_{T}$. The dotted vertical lines indicate the $p_{T}$ intervals in the analysis of HT and JP data.

Comparison of $\Lambda$ and $\bar{\Lambda}$ spin transfer $D_{LL}$ in polarized proton-proton collisions at $\sqrt{s} = 200 GeV$ for (a) positive and (b) negative $\eta$ versus $p_{T}$. The vertical bars and bands indicate the sizes of the statistical and systematic uncertainties, respectively. The $\bar{\Lambda}$ data points have been shifted slightly in $p_{T}$ for clarity. The dotted vertical lines indicate the $p_{T}$ intervals in the analysis of HT and JP data. The horizontal lines show model predictions evaluated at $\eta$ and largest $p_{T}$ of the data.

Comparison of $\Lambda$ and $\bar{\Lambda}$ spin transfer $D_{LL}$ in polarized proton-proton collisions at $\sqrt{s} = 200$ GeV for (a) positive and (b) negative $\eta$ versus $p_{T}$. The vertical bars and bands indicate the sizes of the statistical and systematic uncertainties, respectively. The $\bar{\Lambda}$ data points have been shifted slightly in $p_{T}$ for clarity. The dotted vertical lines indicate the $p_{T}$ intervals in the analysis of HT and JP data. The horizontal lines show model predictions evaluated at $\eta$ and largest $p_{T}$ of the data.

The spin transfer $D_{LL}$ to (a) $\Lambda$ and (b) $\bar{\Lambda}$ hyperons, calculated using the new $\alpha_\Lambda=0.732$ updated on PDG2020, produced at positive pseudorapidity with respect to the polarized proton beam from $MB$, $JP$, and $HT$ data versus hyperon transverse momenta $p_{T}$. The sizes of the statistical and systematic uncertainties are indicated by the vertical bars and bands, respectively. For clarity, the HT data points have been shifted slightly in $p_{T}$. The dotted vertical lines indicate the $p_{T}$ intervals in the analysis of HT and JP data.

The spin transfer $D_{LL}$ to (a) $\Lambda$ and (b) $\bar{\Lambda}$ hyperons, calculated using the new $\alpha_\Lambda=0.732$ updated on PDG2020, produced at positive pseudorapidity with respect to the polarized proton beam from $MB$, $JP$, and $HT$ data versus hyperon transverse momenta $p_{T}$. The sizes of the statistical and systematic uncertainties are indicated by the vertical bars and bands, respectively. For clarity, the HT data points have been shifted slightly in $p_{T}$. The dotted vertical lines indicate the $p_{T}$ intervals in the analysis of HT and JP data.

Comparison of $\Lambda$ and $\bar{\Lambda}$ spin transfer $D_{LL}$, calculated using the new $\alpha_\Lambda=0.732$ updated on PDG2020, in polarized proton-proton collisions at $\sqrt{s} = 200$ GeV for (a) positive and (b) negative $\eta$ versus $p_{T}$. The vertical bars and bands indicate the sizes of the statistical and systematic uncertainties, respectively. The $\bar{\Lambda}$ data points have been shifted slightly in $p_{T}$ for clarity. The dotted vertical lines indicate the $p_{T}$ intervals in the analysis of HT and JP data. The horizontal lines show model predictions evaluated at $\eta$ and largest $p_{T}$ of the data.

Transverse momentum spectra for charged pions at midrapidity.

Transverse momentum spectra for protons at midrapidity.

Transverse momentum spectra for antiprotons at midrapidity.

Transverse momentum spectra for kaons at midrapidity.

dN/dy of pi+ and proton, normalised by N_part/2.

dN/dy of pi+ and proton, normalised by N_part/2.

dN/dy of pi+ and proton, normalised by N_part/2.

dN/dy of K+ and K-, normalised by N_part/2.

dN/dy of K+ and K-, normalised by N_part/2.

dN/dy of K+ and K-, normalised by N_part/2.

Transverse momentum of pi+, K+ and proton.

Transverse momentum of pi+, K+ and proton.

Transverse momentum of pi+, K+ and proton.

K-/K+, K-/pi-, p/pi+ and K+/pi+ vs. N_part.

K-/K+, K-/pi-, p/pi+ and K+/pi+ vs. N_part.

K-/K+, K-/pi-, p/pi+ and K+/pi+ vs. N_part.

dN/d_eta, normalised by N_part/2, as a function of sqrt(S_NN).

dN/dy, normalised by N_part/2; and m_t - m of pi+ and pi- as functions of sqrt(S_NN).

dN/dy, normalised by N_part/2; and m_t - m of K+ and K- as functions of sqrt(S_NN).

pi-/pi+, pbar/p, K-/K+, K+/pi+ and K-/pi- as functions of sqrt(S_NN).

v1 as a function of eta and eta/y-beam from different experiments and methods.

v2 as a function of p_T for different particle species and experiments.

v2 as a function of s_NN from different experiments.

T-mu_B plot for heavy-ion collisions at different s_NN.

$\langle cos(\phi_{\alpha}+\phi_{\beta}−2\phi_{c})\rangle$ as a function of reference multiplicity for different charge combinations, after corrections for acceptance effects. In the legend the signs indicate the charge of particles $\alpha$, $\beta$, and c. The results shown are for Au+Au collisions at 200 GeV obtained in the Full Field.

$\langle cos(\phi_{\alpha}-\phi_{\beta})\rangle$ as a function of centrality for different charge combinations and FF and RFF configurations. The data points corresponding to different charge and field configurations are slightly shifted in the horizontal direction with respect to each other for clarity. The error bars are statistical. Also shown are model predictions described in Section VII.

A comparison of the correlations obtained by selecting the third particle from the main TPC or from the Forward TPCs. The TPC and FTPC points are shifted horizontally relative to one another for clarity purposes. The error bars are statistical.

A comparison of the correlations obtained by selecting the third particle from the main TPC or from the Forward TPCs after scaling by the flow of the third particle. The shaded areas represent the uncertainty from $v_{2,c}$ scaling (see text for details). The TPC and FTPC points are shifted horizontally relative to one another for clarity purposes. The error bars are statistical.

$\langle cos(\phi_{\alpha}+\phi_{\beta}−2\Psi_{RP})\rangle$ in Au+Au and Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV calculated using Eq. 7. The error-bars show the statistical errors. The shaded area reflects the uncertainty in the elliptic flow values used in calculations, with lower (in magnitude) limit obtained with elliptic flow from two-particle correlations and upper limit from four-particle cumulants. For details, see Section IV. Thick solid (Au+Au) and dashed (Cu+Cu) lines represent possible non-reaction-plane dependent contribution from many-particle clusters as estimated by HIJING, see Section VII A.

$\langle cos(\phi_{\alpha}+\phi_{\beta}−2\Psi_{RP})\rangle$ in Au+Au and Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62 GeV calculated using Eq. 7. The error-bars show the statistical errors. The shaded area reflects the uncertainty in the elliptic flow values used in calculations. For details, see Section IV. Thick solid (Au+Au) and dashed (Cu+Cu) lines represent possible non-reaction-plane dependent contribution from many-particle clusters as estimated by HIJING, see Section VII A.

Au+Au and Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV. The correlations are scaled with the number of participants and are plotted as function of centrality. The error-bars indicate the statistical errors. The shaded area reflects the uncertainty in the elliptic flow values used in calculations. For details, see Section IV. Thick solid (Au+Au) and dashed (Cu+Cu) lines represent possible non-reaction-plane dependent contribution from many-particle clusters as estimated by HIJING, see Section VII A.

Au+Au and Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV. The correlations are scaled with the number of participants and are plotted as function of number of participants. The error-bars indicate the statistical errors. The shaded area reflects the uncertainty in the elliptic flow values used in calculations. For details, see Section IV. Thick solid (Au+Au) and dashed (Cu+Cu) lines represent possible non-reaction-plane dependent contribution from many-particle clusters as estimated by HIJING, see Section VII A.

Au+Au at 200 GeV. The correlations dependence on pseudorapidity separation $\Delta\eta$ = |$\eta_{\alpha}$ − $\eta_{\beta}$| for centrality 30-50%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Au+Au at 200 GeV. The correlations dependence on pseudorapidity separation $\Delta\eta$ = |$\eta_{\alpha}$ − $\eta_{\beta}$| for centrality 10-30%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Au+Au at 200 GeV. The correlations dependence on (p$_{t,\alpha}$ + p$_{t,\beta}$)/2 for centrality 30-50%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Au+Au at 200 GeV. The correlations dependence on (p$_{t,\alpha}$ + p$_{t,\beta}$)/2 for centrality 10-30%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Au+Au at 200 GeV. The correlations dependence on |p$_{t,\alpha}$ - p$_{t,\beta}$| for centrality 30-50%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Au+Au at 200 GeV. The correlations dependence on |p$_{t,\alpha}$ - p$_{t,\beta}$| for centrality 10-30%. The shaded band indicates uncertainty associated with $v_{2}$ measurements and has been calculated using two- and four-particle cumulant results as the limits.

Three-particle correlator in Au+Au and Cu+Cu collisions compared to HIJING calculations shown as thick lines. All three particles are taken in the main TPC region, |$\eta$| < 1.0. The correlator has been scaled with number of participants and are plotted versus number of participants. In this representation HIJING results for Au+Au and Cu+Cu collisions coincide in the region of overlap.

$\langle cos(\phi_{\alpha} + \phi_{\beta} − 2\Psi_{RP})\rangle$ calculated for 200 GeV Au+Au events with event generators HIJING (with and without an “elliptic flow afterburner”), UrQMD, and MEVSIM. Blue symbols mark opposite-charge correlations, and red are same-charge. Solid stars represent the values from the data to facilitate comparison. Acceptance cuts of 0.15 < p$_{t}$ < 2 GeV/c and |$\eta$| < 1.0 were used in all cases. For MEVSIM, HIJING, and UrQMD points the true reaction plane from the generated event was used for $\Psi_{RP}$ . Thick solid lighter colored lines represent possible non-reactionplane dependent contribution from many-particle clusters as estimated by HIJING and discussed in Section VII A. Corresponding estimates from UrQMD are about factor of two smaller.

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.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Photon conversion point radius distribution from real data and MC simulation in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Conversion probability correction factor for $\pi^0$ as a function of $p_{T}$

The overall detection efficiency of $\pi^0$ from embedding studt in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The overall detection efficiency of $\pi^0$ from embedding studt in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The overall detection efficiency of $\pi^0$ from embedding studt in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Invariant yield of STAR $\pi^0$ as a function of $p_{T}$ at midrapidity for different collision centrality bins in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Invariant yield of STAR $\pi^0$ as a function of $p_{T}$ at midrapidity for different collision centrality bins in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Invariant yield of STAR $\pi^0$ as a function of $p_{T}$ at midrapidity for different collision centrality bins in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The ratios of STAR $\pi^0$ spectra over $\pi^{\pm}$ from STAR and $\pi^0$ from PHENIX for different collision centrality bins in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The nuclear modification factor $R_{CP}$ as a functon of $p_{T}$ of STAR $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The nuclear modification factor $R_{AA}$ as a functon of $p_{T}$ of STAR $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Fig. 2. (Color online.) Event-by-event photon multiplicity distributions (solid circles) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV} .$ The distributions for top $0-5 \%$ central $\mathrm{Au}+$ Au collisions and top $0-10 \%$ central $\mathrm{Cu}+\mathrm{Cu}$ collisions are also shown (open circles). The photon multiplicity distributions for central collisions are observed to be Gaussian (solid line). Only statistical errors are shown. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 2. (Color online.) Event-by-event photon multiplicity distributions (solid circles) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV} .$ The distributions for top $0-5 \%$ central $\mathrm{Au}+$ Au collisions and top $0-10 \%$ central $\mathrm{Cu}+\mathrm{Cu}$ collisions are also shown (open circles). The photon multiplicity distributions for central collisions are observed to be Gaussian (solid line). Only statistical errors are shown. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 2. (Color online.) Event-by-event photon multiplicity distributions (solid circles) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV} .$ The distributions for top $0-5 \%$ central $\mathrm{Au}+$ Au collisions and top $0-10 \%$ central $\mathrm{Cu}+\mathrm{Cu}$ collisions are also shown (open circles). The photon multiplicity distributions for central collisions are observed to be Gaussian (solid line). Only statistical errors are shown. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several centrality classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several centrality classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several centrality classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 4. (Color online.) Top panel: The number of photons divided by $\left\langle N_{\text{part}}\right\rangle / 2$ as a function of average number of participating nucleons for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s} _{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$ for $-3.7<\eta<-2.3 .$ Errors shown are systematic only and include those for $\left\langle N_{\text {part }}\right\rangle .$ Results from HIJING are shown as lines (solid for $\mathrm{Au}+\mathrm{Au}$ and dashed for $\mathrm{Cu}+\mathrm{Cu}$ ). Bottom panel: Same as above, for both photons and charged particles from PHOBOS [10]. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 4. (Color online.) Top panel: The number of photons divided by $\left\langle N_{\text{part}}\right\rangle / 2$ as a function of average number of participating nucleons for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s} _{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$ for $-3.7<\eta<-2.3 .$ Errors shown are systematic only and include those for $\left\langle N_{\text {part }}\right\rangle .$ Results from HIJING are shown as lines (solid for $\mathrm{Au}+\mathrm{Au}$ and dashed for $\mathrm{Cu}+\mathrm{Cu}$ ). Bottom panel: Same as above, for both photons and charged particles from PHOBOS [10]. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 5. (Color online.) Photon pseudorapidity distributions normalized by the average number of participating nucleon pairs for different collision centralities are plotted as a function of pseudorapidity shifted by the beam rapidity (-5.36 for $200 \mathrm{GeV}$ and -4.19 for $62.4 \mathrm{GeV}$ ) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ collisions at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$. Errors are systematic only, statistical errors are negligible in comparison. For clarity of presentation, results for only four centralities are shown. The $\mathrm{Cu}+$ Cu data are shifted by 0.1 unit in $\eta-y_{\text {beam }} .$ The solid line is a second order polynomial fit to the data (see text for details).

Fig. 5. (Color online.) Photon pseudorapidity distributions normalized by the average number of participating nucleon pairs for different collision centralities are plotted as a function of pseudorapidity shifted by the beam rapidity (-5.36 for $200 \mathrm{GeV}$ and -4.19 for $62.4 \mathrm{GeV}$ ) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ collisions at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$. Errors are systematic only, statistical errors are negligible in comparison. For clarity of presentation, results for only four centralities are shown. The $\mathrm{Cu}+$ Cu data are shifted by 0.1 unit in $\eta-y_{\text {beam }} .$ The solid line is a second order polynomial fit to the data (see text for details).

xT scaling (with xT.

J/psi nuclear modification factor RAA from sqrt(s).

J/psi nuclear modification factor RAA from sqrt(s).

J/psi-hadron azimuthal correlations in sqrt(s).

Contribution from B-hadron feed-down to inclusive J/psi production in sqrt(s).

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

HBT parameters for all centralities in 62.4 GeV Cu+Cu.

HBT parameters for all centralities in 200 GeV Cu+Cu.

HBT parameters for 200 GeV Au+Au, 0-5% and Cu+Cu, 0-10%

HBT parameters for 200 GeV and 62.4 GeV Cu+Cu, 0-10% and 62.4 GeV Au+Au, 0-5%

HBT radii for radii ratios in 200 GeV and 62.4 GeV Cu+Cu, 0-10% and 200 GeV & 62.4 GeV Au+Au, 0-5%. Note that on Figure 8 in the paper the radii ratios are shown while here we list the radii values.

Energy dependence of the pi- HBT Radii for Volume estimate (Rside*Rside*Rlong) in central Au+Au, Pb+Pb, and Pb+Au collisions (AGS,SPS and RHIC) at midrapidity and k_T ~ 0.2-0.3 GeV/c. Note that on Figure 9 in the paper Rside*Rside*Rlong is shown while here we list the Rside and Rlong values.

HBT Radii for Volume estimate(Rside*Rside*Rlong & Rout*Rout*Rlong) comparisons for 200 GeV and 62.4 GeV Au+Au with Npart and dNch/deta. Note that on Figure 10 in the paper specific combinations of the HBT radii are shown while here we list the radii values.

HBT Radii for Volume estimate(Rside*Rside*Rlong & Rout*Rout*Rlong) comparisons for 200 GeV and 62.4 GeV Au+Au & Cu+Cu with dNch/deta. Note that on Figure 11 in the paper specific combinations of the HBT radii are shown while here we list the radii values.

HBT radii comparisons for 200 GeV and 62.4 GeV Au+Au & Cu+Cu with dNch/deta1/3

(Color online) Measured dynamical $K/\pi$ fluctuations in terms of $\nu_{dyn,K\pi}$ for 62.4 and 200 GeV Au+Au compared with $\sigma^{2}_{dyn}$ from central Pb+Pb collisions at 6.3, 7.6, 8.8, 12.3, and 17.3 GeV from NA49 [7]. Statistical errors are shown for the STAR data. Statistical and systematic errors are shown for the NA49 results. The solid line corresponds to a fit to the STAR data of the form $c + d/(dN/d\eta)$.

(Color online) The dN/dη scaled dynamical $K/\pi$ fluctuations for summed charges (stars), same signs (circles), and opposite signs (squares) as a function of $dN/d\eta$. The errors shown are statistical. The open and filled symbols refer to Au+Au collisions at 62.4 GeV and 200 GeV respectively. The dash-dot, dotted, and dashed lines represents HIJING calculations for summed charges, same signs, and opposite signs respectively.

Efficiency corrected $t_{\perp}$ spectrum for $\rho^{0}$ from (top) minium bias and (bottom) topology data, for mid-rapidity (left) and larger rapidity (right) samples. The points are the data, while the solid lines are the results of fits to Eq. (3).

Efficiency corrected $t_{\perp}$ spectrum for $\rho^{0}$ from (top) minium bias and (bottom) topology data, for mid-rapidity (left) and larger rapidity (right) samples. The points are the data, while the solid lines are the results of fits to Eq. (3).

Efficiency corrected $t_{\perp}$ spectrum for $\rho^{0}$ from (top) minium bias and (bottom) topology data, for mid-rapidity (left) and larger rapidity (right) samples. The points are the data, while the solid lines are the results of fits to Eq. (3).

Efficiency corrected $t_{\perp}$ spectrum for $\rho^{0}$ from (top) minium bias and (bottom) topology data, for mid-rapidity (left) and larger rapidity (right) samples. The points are the data, while the solid lines are the results of fits to Eq. (3).

Upper panels. $N_{\scriptsize{\mbox{part}}}$ scaled ($R^{N_{\scriptsize{\mbox{part}}}}_{AA}$) nuclear modification factors as a function of $p_{T}$ of $\phi$ mesons for $0-10\%$ and $20-30\%$ $Cu+Cu$ and $Au+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV. Lower panel. Same as above for $N_{\scriptsize{\mbox{bin}}}$ scaled ($R^{N_{\scriptsize{\mbox{bin}}}}_{AA}$) nuclear modification factor. The error bars represent the statistical and systematic errors added in quadrature. The shaded band in upper panel around 1 at $p_{T}=4.5-5.5$ GeV/$c$ in the right side reflects the uncertainty in $N_{\scriptsize{\mbox{part}}}$ and that on the lower panel for $N_{\scriptsize{\mbox{bin}}}$ calculation for central $Au+Au$ collisions. The respective uncertainties for central $Cu+Cu$ collisions are of similar order.

Upper panels. $N_{\scriptsize{\mbox{part}}}$ scaled ($R^{N_{\scriptsize{\mbox{part}}}}_{AA}$) nuclear modification factors as a function of $p_{T}$ of $\phi$ mesons for $0-10\%$ and $20-30\%$ $Cu+Cu$ and $Au+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV. Lower panel. Same as above for $N_{\scriptsize{\mbox{bin}}}$ scaled ($R^{N_{\scriptsize{\mbox{bin}}}}_{AA}$) nuclear modification factor. The error bars represent the statistical and systematic errors added in quadrature. The shaded band in upper panel around 1 at $p_{T}=4.5-5.5$ GeV/$c$ in the right side reflects the uncertainty in $N_{\scriptsize{\mbox{part}}}$ and that on the lower panel for $N_{\scriptsize{\mbox{bin}}}$ calculation for central $Au+Au$ collisions. The respective uncertainties for central $Cu+Cu$ collisions are of similar order.

Upper panels. $N_{\scriptsize{\mbox{part}}}$ scaled ($R^{N_{\scriptsize{\mbox{part}}}}_{AA}$) nuclear modification factors as a function of $p_{T}$ of $\phi$ mesons for $0-10\%$ and $20-30\%$ $Cu+Cu$ and $Au+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV. Lower panel. Same as above for $N_{\scriptsize{\mbox{bin}}}$ scaled ($R^{N_{\scriptsize{\mbox{bin}}}}_{AA}$) nuclear modification factor. The error bars represent the statistical and systematic errors added in quadrature. The shaded band in upper panel around 1 at $p_{T}=4.5-5.5$ GeV/$c$ in the right side reflects the uncertainty in $N_{\scriptsize{\mbox{part}}}$ and that on the lower panel for $N_{\scriptsize{\mbox{bin}}}$ calculation for central $Au+Au$ collisions. The respective uncertainties for central $Cu+Cu$ collisions are of similar order.

Upper panel. The ratio of the yields of $K^{−}$, $\phi$, $\bar{\Lambda}$ and $\Xi+\bar{\Xi}$ normalized to $N_{\scriptsize{\mbox{part}}}$ in nucleus-nucleus collisions to corresponding yields in inelastic $p+p$ collisions as a function of $N_{\scriptsize{\mbox{part}}}$ at 200 GeV. Lower panel. Same as above for $\phi$ mesons in $Cu+Cu$ collisions at 200 and 62.4 GeV. The $p+p$ collision data at 200 GeV are from STAR [5] and at 62.4 GeV from ISR [29]. The error bars shown represent the statistical and systematic errors added in quadrature.

Upper panel. The ratio of the yields of $K^{−}$, $\phi$, $\bar{\Lambda}$ and $\Xi+\bar{\Xi}$ normalized to $N_{\scriptsize{\mbox{part}}}$ in nucleus-nucleus collisions to corresponding yields in inelastic $p+p$ collisions as a function of $N_{\scriptsize{\mbox{part}}}$ at 200 GeV. Lower panel. Same as above for $\phi$ mesons in $Cu+Cu$ collisions at 200 and 62.4 GeV. The $p+p$ collision data at 200 GeV are from STAR [5] and at 62.4 GeV from ISR [29]. The error bars shown represent the statistical and systematic errors added in quadrature.

Upper panels: same-event (full points) and mixed-event (solid line) $K^{+}K^{-}$ invariant mass distributions at 0.6 < $p_{T}$ < 1.4 GeV/c in p + p 200 GeV collisions (a), 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 62.4 GeV collisions (60–80%) (c), and 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 200 GeV collisions (0–10%) (e). Lower panels: the corresponding $\phi$ meson mass peaks after subtracting the background. Dashed curves show a Breit-Wigner + linear background function fit in (b), (d). In (f), both linear and quadratic backgrounds are shown as dashed and dot-dashed lines, respectively.

Upper panels: same-event (full points) and mixed-event (solid line) $K^{+}K^{-}$ invariant mass distributions at 0.6 < $p_{T}$ < 1.4 GeV/c in p + p 200 GeV collisions (a), 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 62.4 GeV collisions (60–80%) (c), and 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 200 GeV collisions (0–10%) (e). Lower panels: the corresponding $\phi$ meson mass peaks after subtracting the background. Dashed curves show a Breit-Wigner + linear background function fit in (b), (d). In (f), both linear and quadratic backgrounds are shown as dashed and dot-dashed lines, respectively.

Upper panels: same-event (full points) and mixed-event (solid line) $K^{+}K^{-}$ invariant mass distributions at 0.6 < $p_{T}$ < 1.4 GeV/c in p + p 200 GeV collisions (a), 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 62.4 GeV collisions (60–80%) (c), and 0.8 < $p_{T}$ < 1.2 GeV/c in Au + Au 200 GeV collisions (0–10%) (e). Lower panels: the corresponding $\phi$ meson mass peaks after subtracting the background. Dashed curves show a Breit-Wigner + linear background function fit in (b), (d). In (f), both linear and quadratic backgrounds are shown as dashed and dot-dashed lines, respectively.

Reconstruction efficiency including acceptance of $\phi$ meson as a function of $p_{T}$ in several centrality bins of Au + Au, d + Au, and p + p 200 GeV collisions.

Reconstruction efficiency including acceptance of $\phi$ meson as a function of $p_{T}$ in several centrality bins of Au + Au, d + Au, and p + p 200 GeV collisions.

Reconstruction efficiency including acceptance of $\phi$ meson as a function of $p_{T}$ in several centrality bins of Au + Au, d + Au, and p + p 200 GeV collisions.

Event plane $\Phi_{2}$ resolution as a function of centrality in Au + Au 200 GeV collisions, where the vertical axis starts from 0.3 for clarity.

$\phi−\Phi_{2}$ distribution for $\phi$ meson at 1.5 < $p_{T}$ < 2.0 GeV/c in Au + Au collisions (0–80$\%$) at 200 GeV. The line is the fitting result. Error bars are statistical only.

<cos2($\phi_{K^{+}K^{-}}$-$\Phi_{2}$)> (full red points) and <sin2($\phi_{K^{+}K^{-}}$-$\Phi_{2}$)> (open blue points) as a function of $m_{inv}$ of $K^{+}K^{-}$ pairs at 0.5 < $p_{T}$ < 1.0 GeV/c in Au + Au 200 GeV collisions (0–80\%), where the solid curve is the result of fitting by Eq. (10). The arrow shows the position of the $\phi$ invariant mass peak. The dashed line shows the zero horizontal line.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Masses and widths (FWHMs) of $\phi$ as a function of $p_{T}$ in p + p 200 GeV (NSD), d + Au 200 GeV (0–20$\%$), Au + Au 62.4 GeV (0–20$\%$), and Au + Au 200 GeV (0–5$\%$) collisions, with the corresponding $PDG$ values.

Invariant mass distributions of $\phi$ meson at 0.6 < $p_{T}$ < 1.0 GeV/c in p + p 200 GeV (NSD) and Au + Au 200 GeV (0–5$\%$) collisions. Solid symbols: experimental data. Open symbols: MC simulation. Curves are the results of a Breit-Wigner function fit. Note: Two sets of MC data are shown for Au + Au 200 GeV, and see text for details.

Invariant mass distributions of $\phi$ meson at 0.6 < $p_{T}$ < 1.0 GeV/c in p + p 200 GeV (NSD) and Au + Au 200 GeV (0–5$\%$) collisions. Solid symbols: experimental data. Open symbols: MC simulation. Curves are the results of a Breit-Wigner function fit. Note: Two sets of MC data are shown for Au + Au 200 GeV, and see text for details.

Invariant mass distributions of $\phi$ meson at 0.6 < $p_{T}$ < 1.0 GeV/c in p + p 200 GeV (NSD) and Au + Au 200 GeV (0–5$\%$) collisions. Solid symbols: experimental data. Open symbols: MC simulation. Curves are the results of a Breit-Wigner function fit. Note: Two sets of MC data are shown for Au + Au 200 GeV, and see text for details.

Invariant mass distributions of $\phi$ meson at 0.6 < $p_{T}$ < 1.0 GeV/c in p + p 200 GeV (NSD) and Au + Au 200 GeV (0–5$\%$) collisions. Solid symbols: experimental data. Open symbols: MC simulation. Curves are the results of a Breit-Wigner function fit. Note: Two sets of MC data are shown for Au + Au 200 GeV, and see text for details.

Invariant mass distributions of $\phi$ meson at 0.6 < $p_{T}$ < 1.0 GeV/c in p + p 200 GeV (NSD) and Au + Au 200 GeV (0–5$\%$) collisions. Solid symbols: experimental data. Open symbols: MC simulation. Curves are the results of a Breit-Wigner function fit. Note: Two sets of MC data are shown for Au + Au 200 GeV, and see text for details.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

$\phi$ meson transverse mass distributions for different collision systems and different energies. For clarity, distributions for some centrality bins have been scaled by factors indicated in the figure. Curves represent the exponential (solid) and Levy (dashed) function fits to the distributions. Error bars are statistical only. Note that a scale factor of 1.09 is applied to the $\phi$ meson spectra for Au + Au collisions at 200 GeV in Ref. [50] to correct for the kaon identification efficiency effect missed previously.

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

Comparison of transverse momentum spectra shape among different 200 GeV collision systems: Au + Au (0–5$\%$), d + Au (0–20$\%$), and p + p (inelastic). The spectra are normalized by $N_{bin}$ (top panel) and $N_{part}/2$ (bottom panel).

$N_{part}$ dependence of $(dN/dy)/0.5N_{part}$ in five different collision systems: Au + Au 62.4, 130, 200 GeV; p + p 200 GeV (inelastic); and d + Au 200 GeV. Statistical and systematic errors are included.

$N_{part}$ dependence of $(dN/dy)/0.5N_{part}$ in five different collision systems: Au + Au 62.4, 130, 200 GeV; p + p 200 GeV (inelastic); and d + Au 200 GeV. Statistical and systematic errors are included.

$N_{part}$ dependence of $(dN/dy)/0.5N_{part}$ in five different collision systems: Au + Au 62.4, 130, 200 GeV; p + p 200 GeV (inelastic); and d + Au 200 GeV. Statistical and systematic errors are included.

$N_{part}$ dependence of $(dN/dy)/0.5N_{part}$ in five different collision systems: Au + Au 62.4, 130, 200 GeV; p + p 200 GeV (inelastic); and d + Au 200 GeV. Statistical and systematic errors are included.

$N_{part}$ dependence of $(dN/dy)/0.5N_{part}$ in five different collision systems: Au + Au 62.4, 130, 200 GeV; p + p 200 GeV (inelastic); and d + Au 200 GeV. Statistical and systematic errors are included.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: $N_{part}$ dependence of $<p_{T}>_{\phi}$ in different collision systems. Bottom panel: Hadron mass dependence of $<p_{T}>$ in central Au + Au collisions at 62.4 and 200 GeV. The band and curve show two hydrodynamic model calculations for central Au + Au collisions at 200 GeV. Note: Hadron masses for the Au + Au 62.4 GeV data are shifted slightly in the $x$-axis direction for clarity, and systematic errors are included for $\phi$.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of the ratio $\phi/\pi^{-}$ in A + A (full points) and p + p (open points) collisions. Stars are data from the STARexperiment at RHIC.Bottom panel: $N_{part}$ dependence of ratio $\phi/\pi^{-}$ in different collision systems. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

Top panel: Energy dependence of ratio $\phi/K^{-}$ in $A + A$ and elementary ($e + e$ and $p + p$) collisions. Stars are data from STAR experiments at RHIC. Bottom panel: $N_{part}$ dependence of ratio $\phi/K^{-}$ in different collision systems. The dashed line shows results from $UrQMD$ model calculations. Systematic errors are included for the STAR data points.

The $\Omega/\phi$ ratio vs $p_{T}$ for three centrality bins in $\sqrt{s_{NN}} = 200 GeV$ Au + Au collisions, where the data points for 40–60$\%$ are shifted slightly for clarity. As shown in the legend, the lines represent results from Hwa and Yang [96], Ko et al. [97], and, for SCF, Refs. [98,99].

The $\Omega/\phi$ ratio vs $p_{T}$ for three centrality bins in $\sqrt{s_{NN}} = 200 GeV$ Au + Au collisions, where the data points for 40–60$\%$ are shifted slightly for clarity. As shown in the legend, the lines represent results from Hwa and Yang [96], Ko et al. [97], and, for SCF, Refs. [98,99].

$p_{T}$ dependence of the nuclear modification factor $R_{cp}$ in Au + Au 200 GeV collisions. The top and bottom panels present $R_{cp}$ from midperipheral and most-peripheral collisions, respectively. See legend for symbol and line designations. The rectangular bands showthe uncertainties of binary and participant scalings. Statistical and systematic errors are included.

$p_{T}$ dependence of the nuclear modification factor $R_{cp}$ in Au + Au 200 GeV collisions. The top and bottom panels present $R_{cp}$ from midperipheral and most-peripheral collisions, respectively. See legend for symbol and line designations. The rectangular bands showthe uncertainties of binary and participant scalings. Statistical and systematic errors are included.

$p_{T}$ dependence of the nuclear modification factor $R_{cp}$ in Au + Au 62.4 and 200 GeV and d + Au 200 GeV collisions, where rectangular bands represent the uncertainties of binary and participant scalings (see legend). Statistical and systematic errors are included.

$p_{T}$ dependence of the nuclear modification factor $R_{cp}$ in Au + Au 62.4 and 200 GeV and d + Au 200 GeV collisions, where rectangular bands represent the uncertainties of binary and participant scalings (see legend). Statistical and systematic errors are included.

$p_{T}$ dependence of the nuclear modification factor $R_{cp}$ in Au + Au 62.4 and 200 GeV and d + Au 200 GeV collisions, where rectangular bands represent the uncertainties of binary and participant scalings (see legend). Statistical and systematic errors are included.

$p_{T}$ dependence of the nuclear modification factor $R_{AB}$ for $\phi$ in Au + Au 200 GeV and d + Au 200 GeV collisions. For comparison, data points for $\pi^{+}+\pi^{-}$ in d + Au 200 GeV and $p+\bar{p}$ in d + Au 200 GeV are also shown (see legend). Rectangular bands show the uncertainties of binary (solid line) and participant (dot-dash line) scalings. Systematic errors are included for $\phi$, $\pi^{+}+\pi^{-}$, and $p+\bar{p}$.

$p_{T}$ dependence of the nuclear modification factor $R_{AB}$ for $\phi$ in Au + Au 200 GeV and d + Au 200 GeV collisions. For comparison, data points for $\pi^{+}+\pi^{-}$ in d + Au 200 GeV and $p+\bar{p}$ in d + Au 200 GeV are also shown (see legend). Rectangular bands show the uncertainties of binary (solid line) and participant (dot-dash line) scalings. Systematic errors are included for $\phi$, $\pi^{+}+\pi^{-}$, and $p+\bar{p}$.

$p_{T}$ dependence of the elliptic flow $v_{2}$ of $\phi$, $\Lambda$, and $K_{s}^{0}$ in Au + Au collisions (0–80$\%$) at 200 GeV. Data points for $\phi$ are from the reaction plane method (full up-triangles) and invariant mass method (full circles), where data points from the reaction plane method are shifted slightly along the $x$ axis for clarity. Vertical error bars represent statistical errors, while the square bands represent systematic uncertainties. The magenta curved band represents the $v_{2}$ of $\phi$ meson from the AMPT model with a string melting mechanism [116]. The dash and dot curves represent parametrizations inspired by number-of-quark scaling ideas from Ref. [117] for NQ = 2 and NQ = 3, respectively.

Elliptic flow $v_{2}$ as a function of $p_{T}$ , $v_{2}$($p_{T}$), for the $\phi$ meson from different centralities. The vertical error bars represent the statistical errors, while the square bands represent the systematic uncertainties. For clarity, data points of 10–40$\%$ are shifted in the $p_{T}$ direction slightly.

Elliptic flow $v_{2}$ as a function of $p_{T}$ , $v_{2}$($p_{T}$), for the $\phi$ meson from different centralities. The vertical error bars represent the statistical errors, while the square bands represent the systematic uncertainties. For clarity, data points of 10–40$\%$ are shifted in the $p_{T}$ direction slightly.

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.

Delta++ reconstruction efficiency times detector acceptance as a function of the invariant mass for minimum bias d+Au collisions. The error shown is due to the available statistics in the simulation.

Delta++ reconstruction efficiency times detector acceptance as a function of the invariant mass for p+p collisions. The error shown is due to the available statistics in the simulation.

Lambda* reconstruction efficiency times detector acceptance as a function of pT for minimum bias d+Au collisions. The error shown is due to the available statistics in the simulation.

Sigma*+- reconstruction efficiency times detector acceptance as a function of pT for minimum bias d+Au collisions. The error shown is due to the available statistics in the simulation.

Sigma*+- reconstruction efficiency times detector acceptance as a function of pT for d+Au at 40-100% centrality. The error shown is due to the available statistics in the simulation.

Sigma*+- reconstruction efficiency times detector acceptance as a function of pT for d+Au at 20-40% centrality . The error shown is due to the available statistics in the simulation.

Sigma*+- reconstruction efficiency times detector acceptance as a function of pT for d+Au at 0-20% centrality . The error shown is due to the available statistics in the simulation.

Raw pi+pi- invariant mass distribution at midrapidity after subtraction of the like-sign reference distribution for minimum bias d+Au collisions with the hadronic cocktail fit.

rho0 mass as a function of pT at |y| < 0.5 for d+A at 0-20 % centrality.

rho0 mass as a function of pT at |y| < 0.5 for d+A at 20-40 % centrality.

rho0 mass as a function of pT at |y| < 0.5 for d+A at 40-100 % centrality.

rho0 mass as a function of pT at |y| < 0.5 for minimum bias d+A collisions.

rho0 mass as a function of pT at |y| < 0.5 for p+p collisions.

rho0 mass as a function of pT at |y| < 0.5 for high multiplicity p+p collisions.

rho0 mass as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Invariant pi+pi- mass distribution after background subtraction for minimum bias p+p collisions measured by NA27.

Mixed-event background subtracted K pi raw invariant mass distribution for a particular K*0 pT bin at |y| < 0.5 for minimum bias d+Au interactions.

Mixed-event background subtracted KS*0 pi raw invariant mass distribution integrated over nvariant mass distribution for a particular K*+- pT bin at |y| < 0.5 for minimum bias d+Au interactions.

K*0 mass as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

K*+- mass as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

K*0 width as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

K*+- width as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Mixed-event background subtracted p pi raw invariant mass distribution for a particular pT bin at |y| < 0.5 for minimum bias d+Au interactions.

Mixed-event background subtracted p pi raw invariant mass distribution for a particular pT bin at |y| < 0.5 for minimum bias p+p collisions.

Delta++ mass as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Delta++ width as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Delta++ mass as a function of pT at |y| < 0.5 for minimum bias p+p collisions.

Delta++ width as a function of pT at |y| < 0.5 for minimum bias p+p collisions.

Mixed-event background subtracted Lambda pi raw invariant mass distribution integrated over the Sigma* pT at |y| < 0.5 for minimum bias d+Au interactions.

Mixed-event background subtracted p K raw invariant mass distribution integrated over the Lambda* pT at |y| < 0.5 for minimum bias d+Au interactions.

Sigma* mass as a function of pT at |y| < 0.75 for minimum bias d+Au collisions.

Lambda* mass as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Lambda* width as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

rho0 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 40-100% centrality.

rho0 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 20-40% centrality.

rho0 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 0-20% centrality.

rho0 invariant yields as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

(K*0+K*0bar)/2 invariant yields as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

(K*0+K*0bar)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 0-20% centrality.

(K*0 + K*0bar)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 20-40% centrality.

(K*0 + K*0bar)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 40-100% centrality.

(K*+ + K*-)/2 invariant yields as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

(K*+ +K*-)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 0-20% centrality.

(K*+ +K*-)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 20-40% centrality.

(K*+ +K*-)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 40-100% centrality.

(Delta++ + Delta--bar)/2 invariant yields as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

(Delta++ + Delta--bar)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 0-20% centrality.

(Delta++ + Delta--bar)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 20-40% centrality.

(Delta++ + Delta--bar)/2 invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 40-100% centrality.

(Delta++ + Delta--bar)/2 invariant yields as a function of pT at |y| < 0.5 for minimum bias p+p collisions. The errors shown also include the systematic uncertainties.

Sigma* invariant yields as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Sigma* invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 0-20% centrality.

Sigma* invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 20-40% centrality.

Sigma* invariant yields as a function of pT at |y| < 0.5 for d+Au collisions at 40-100% centrality.

Lambda* invariant yields as a function of pT at |y| < 0.5 for minimum bias d+Au collisions.

Rescaled spectra of pion minimum bias d+Au.

Rescaled spectra of p minimum bias d+Au.

Rescaled spectra of rho0 minimum bias d+Au.

Rescaled spectra of K* minimum bias d+Au.

Rescaled spectra of Delta minimum bias d+Au.

Rescaled spectra of Sigma* minimum bias d+Au.

Rescaled spectra of Lambda* minimum bias d+Au.

pion as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

proton as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

K as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

rho0 as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

K* as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

pion as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

proton as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

K as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

Delta++ as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

Sigma* as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

Lambda* as a function of < dNch > /deta for minimum bias p+p, minimum bias d+Au, and 0-20%, 20-40%, 40-100% of the total d+Au cross section.

rho0/pi ratios in p+p, various centralities in d+Au, and in peripheral Au+Au collisions as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

K*/k- ratios in p+p, various centralities in d+Au and Au+Au collisions as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

Delta++/p ratios in p+p, various centralities in d+Au as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

Sigma*/Lambda ratios in p+p, various centralities in d+Au, and in central Au+Au collisions as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

Lambda*/Lambda ratios in p+p, minimum bias d+Au, and two different centralities in Au+Au collisions as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

Delta++/p ratios normalized by their ratio measured in p+p collisions at the same beam energy as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

K*/K- ratios normalized by their ratio measured in p+p collisions at the same beam energy as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

rho0/pi- ratios normalized by their ratio measured in p+p collisions at the same beam energy as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

Lambda*/Lambda ratios normalized by their ratio measured in p+p collisions at the same beam energy as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

Sigma*/Lambda ratios normalized by their ratio measured in p+p collisions at the same beam energy as a function of dNch/deta. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

RdAu for pion.

RdAu for charged hadrons.

RdAu for charged rho0.

RdAu for pion.

RdAu for charged hadrons.

RdAu for charged K*0.

RdAu for charged K*+-.

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.

(a) The solid curve represents the fraction of $\Delta G$ for GRSV-std that has $x > x_{min}$ for scale $Q^{2}$ = 100 $GeV^{2}/c^{2}$. The histograms show the $x_{gluon}$ sampled in the lowest and highest jet $p_{T}$ bins. (b) Confidence levels for several gluon polarization distributions, characterized by their $\Delta G$ at an input scale of 0.4 $GeV^{2}/c^{2}$ [4, 14, 23].

(a) The solid curve represents the fraction of $\Delta G$ for GRSV-std that has $x > x_{min}$ for scale $Q^{2}$ = 100 $GeV^{2}/c^{2}$. The histograms show the $x_{gluon}$ sampled in the lowest and highest jet $p_{T}$ bins. (b) Confidence levels for several gluon polarization distributions, characterized by their $\Delta G$ at an input scale of 0.4 $GeV^{2}/c^{2}$ [4, 14, 23].

(Color online) $\Lambda$ and $\overline{\Lambda}$ spectra on the deuteron and on the gold side in $d$ + Au minimum bias collisions. The data points on the gold side are multiplied by 2 for better visibility. The statistical errors are smaller than the points marking the measurements. The curves show a fit with a Boltzmann function in transverse mass to the data points.

(Color online) (a) Comparison of the measured $\overline{\Lambda}$ yield with model calculations. Statistical errors are shown as vertical error bars, the vertical caps show the quadratic sum of statistical and systematic er- rors including the overall normalization uncertainty. In both panels the target and projectile beam rapidities are indicated by arrows.

(Color online) (b) Comparison of the net $\Lambda$ yield with model calculations. Statistical errors are shown as vertical error bars, the vertical caps show the quadratic sum of statistical and systematic er- rors including the overall normalization uncertainty. In both panels the target and projectile beam rapidities are indicated by arrows.

(Color online) Comparison of $\overline{\Lambda}$ and net $\Lambda$ yields to model calculations for all three centrality classes. Statistical errors are shown as vertical error bars, the vertical caps show the quadratic sum of statistical and systematic errors. Beam rapidity is indicated by arrows.

(Color online) Comparison of $\overline{\Lambda}$ and net $\Lambda$ yields to model calculations for all three centrality classes. Statistical errors are shown as vertical error bars, the vertical caps show the quadratic sum of statistical and systematic errors. Beam rapidity is indicated by arrows.

(Color online) Comparison of $\overline{\Lambda}$ and net $\Lambda$ yields to model calculations for all three centrality classes. Statistical errors are shown as vertical error bars, the vertical caps show the quadratic sum of statistical and systematic errors. Beam rapidity is indicated by arrows.

(Color online) Comparison of $\overline{\Lambda}$ and net $\Lambda$ yields to model calculations for all three centrality classes. Statistical errors are shown as vertical error bars, the vertical caps show the quadratic sum of statistical and systematic errors. Beam rapidity is indicated by arrows.

(Color online) Comparison of $\overline{\Lambda}$ and net $\Lambda$ yields to model calculations for all three centrality classes. Statistical errors are shown as vertical error bars, the vertical caps show the quadratic sum of statistical and systematic errors. Beam rapidity is indicated by arrows.

(Color online) Comparison of $\overline{\Lambda}$ and net $\Lambda$ yields to model calculations for all three centrality classes. Statistical errors are shown as vertical error bars, the vertical caps show the quadratic sum of statistical and systematic errors. Beam rapidity is indicated by arrows.

(Color online) Minimum bias $\overline{\Lambda}/\Lambda$ ratio compared to model calculations. On the deuteron side HIJING/BB shows the best agreement with the results, while on the Au side only AMPT and EPOS give a satisfactory description of the data.

$\overline{\Lambda}/\Lambda$ ratio and net $\Lambda$ and $\overline{\Lambda}$ yields as a function of collision centrality on both the deuteron (left) and the gold side (right). On the deuteron side, centrality is expressed by the number of collisions per deuteron participant, while on the gold side the number of Au participants is chosen. Only statistical errors are shown. The increase in baryon number transport with centrality, shown by the net $\Lambda$ yield, is matched by the increase of $\overline{\Lambda}$-$\Lambda$ pair production, thus keeping the $\overline{\Lambda}/\Lambda$ ratio constant over a wide centrality range.

$\overline{\Lambda}/\Lambda$ ratio and net $\Lambda$ and $\overline{\Lambda}$ yields as a function of collision centrality on both the deuteron (left) and the gold side (right). On the deuteron side, centrality is expressed by the number of collisions per deuteron participant, while on the gold side the number of Au participants is chosen. Only statistical errors are shown. The increase in baryon number transport with centrality, shown by the net $\Lambda$ yield, is matched by the increase of $\overline{\Lambda}$-$\Lambda$ pair production, thus keeping the $\overline{\Lambda}/\Lambda$ ratio constant over a wide centrality range.

$\overline{\Lambda}/\Lambda$ ratio and net $\Lambda$ and $\overline{\Lambda}$ yields as a function of collision centrality on both the deuteron (left) and the gold side (right). On the deuteron side, centrality is expressed by the number of collisions per deuteron participant, while on the gold side the number of Au participants is chosen. Only statistical errors are shown. The increase in baryon number transport with centrality, shown by the net $\Lambda$ yield, is matched by the increase of $\overline{\Lambda}$-$\Lambda$ pair production, thus keeping the $\overline{\Lambda}/\Lambda$ ratio constant over a wide centrality range.

$\overline{\Lambda}/\Lambda$ ratio and net $\Lambda$ and $\overline{\Lambda}$ yields as a function of collision centrality on both the deuteron (left) and the gold side (right). On the deuteron side, centrality is expressed by the number of collisions per deuteron participant, while on the gold side the number of Au participants is chosen. Only statistical errors are shown. The increase in baryon number transport with centrality, shown by the net $\Lambda$ yield, is matched by the increase of $\overline{\Lambda}$-$\Lambda$ pair production, thus keeping the $\overline{\Lambda}/\Lambda$ ratio constant over a wide centrality range.

(Color online) Net $\Lambda dN/dy$ for central, mid-central and peripheral events on both the deuteron and the Au side of the collision. The data are compared to calculations of the distribution of net baryons obtained with the Multichain model [16] with $\alpha$ = 2.9, scaled by 0.1 to account for the conversion from net baryons to net $\Lambda$. An overall scale uncertainty of 10% on the model curves from this conversion is not shown. See text for details.

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

(Color online) Global polarization of $\Lambda$–hyperons as a function of centrality given as fraction of the total inelastic hadronic cross section. Filled circles show the results for Au+Au collisions at $\sqrt{s_{NN}}$=200 GeV (centrality region 20-70%) and open squares indicate the results for Au+Au collisions at $\sqrt{s_{NN}}$=62.4 GeV (centrality region 0-80%). Only statistical uncertainties are shown.

(Color online) Global polarization of $\overline{\Lambda}$–hyperons as a function of $\overline{\Lambda}$ transverse momentum $p^{\overline{\Lambda}}_{t}$. Filled circles show the results for Au+Au collisions at $\sqrt{s_{NN}}$=200 GeV (centrality region 20-70%) and open squares indicate the results for Au+Au collisions at $\sqrt{s_{NN}}$=62.4 GeV (centrality region 0-80%). Only statistical uncertainties are shown.

(Color online) Global polarization of $\overline{\Lambda}$–hyperons as a function of $\overline{\Lambda}$ pseudorapidity $\eta^{\overline{\Lambda}}$. Filled circles show the results for Au+Au collisions at $\sqrt{s_{NN}}$=200 GeV (centrality region 20-70%). A constant line fit to these data points yields $P_{\Lambda}=(1.8\pm 10.8)\times 10^{-3}$ with $\chi^{2}/ndf=5.5/10$. Open squares show the results for Au+Au collisions as $\sqrt{s_{NN}}$=62.4 GeV (centrality region 0-80%). A constant line fit gives $P_{\Lambda}=(-17.6\pm 11.1)\times 10^{-3}$ with $\chi^{2}/ndf=8.0/10$. Only statistical uncertainties are shown.

(Color online) Global polarization of $\overline{\Lambda}$–hyperons as a function of as a function of centrality. Filled circles show the results for Au+Au collisions at $\sqrt{s_{NN}}$=200 GeV (centrality region 20-70%) and open squares indicate the result for Au+Au collisions at $\sqrt{s_{NN}}$=62.4 GeV (centrality region 0-80%). Only statistical uncertainties are shown.

(Color online) Acceptance correction function $A_{0}(p^{H}_{t}, \eta^{H})$ defined in (10) as a function of $\Lambda$ (filled circles) and $\overline{\Lambda}$ (open squares) transverse momentum (top) and pseudorapidity (bottom). The deviation of this function from unity affects the overall scale of the measured global polarization according to Eq. 9. See the text for details and discussions on $A_{0}$ $p^{H}_{t}$ and $\eta^{H}$ dependence.

(Color online) Acceptance correction function $A_{0}(p^{H}_{t}, \eta^{H})$ defined in (10) as a function of $\Lambda$ (filled circles) and $\overline{\Lambda}$ (open squares) transverse momentum (top) and pseudorapidity (bottom). The deviation of this function from unity affects the overall scale of the measured global polarization according to Eq. 9. See the text for details and discussions on $A_{0}$ $p^{H}_{t}$ and $\eta^{H}$ dependence.

(Color online) Function $A_{2}(p^{H}_{t}, \eta^{H})$ defined in (11) as a function of $\Lambda$ (filled circles) and $\overline{\Lambda}$ (open squares) transverse momentum (top) and pseudorapidity (bottom). The deviation of this function from zero defines the contribution to the observable (3) from $P^{(2)}_{H} (p^{H}_{t}, \eta^{H})$ in the expansion (6).

(Color online) Function $A_{2}(p^{H}_{t}, \eta^{H})$ defined in (11) as a function of $\Lambda$ (filled circles) and $\overline{\Lambda}$ (open squares) transverse momentum (top) and pseudorapidity (bottom). The deviation of this function from zero defines the contribution to the observable (3) from $P^{(2)}_{H} (p^{H}_{t}, \eta^{H})$ in the expansion (6).

Measured and calculated asymmetries vs. di-jet pseudorapidity sum for $+\hat{z}$ (left) and $−\hat{z}$ (right) beams. (a,b): Fraction of the calculated di-jet cross section with a quark (gluon) from the $+\hat{z}$ $(−\hat{z})$ beam. (c,d): Unweighted asymmetries compared with pQCD calculations [20] (histograms) for two models of quark Sivers functions fitted to SIDIS results [8]. (e,f): Asymmetries for $|\sin\zeta|$-weighted yields, compared with calculations [20, 21] based on twist-3 quark-gluon correlations. Vertical (horizontal) bars on the data indicate statistical uncertainties (bin widths). The systematic error bands exclude a $\pm12\%$ beam polarization normalization uncertainty.

Centrality and pT dependence of RCP for pions and protons for Au+Au at sqrt(sNN) = 62.4 GeV. For studying the energy dependence, the corresponding RCP for central 0-12% Au+Au at sqrt(sNN) = 200 GeV are shown. RAA for pions at 62.4 GeV (0-10%) and 200 GeV (0-12%) are also shown.

The proton (anti-proton) over pion ratios versus pT at sqrt(sNN) = 62.4 and 200 GeV for central and peripheral collisions.

(color online) The $R_{CP}$ of mid-rapidity $\phi$-mesons produced in $\sqrt{s_{NN}}$ = 200 GeV Au+Au collisions: (top) 0-5% vs. 40-60% and (bottom) 0-5% vs. 60-80% The shaded bands represent the uncertainties in the Glauber model calculations for $\langle N_{bin}\rangle$ and $\langle N_{part}\rangle$ [32]. Also shown are results for $\Lambda$ and $K^{0}_{S}$ [11] and protons and $\pi^{+}$ [33].

(color online) The $R_{CP}$ of mid-rapidity $\phi$-mesons produced in $\sqrt{s_{NN}}$ = 200 GeV Au+Au collisions: (top) 0-5% vs. 40-60% and (bottom) 0-5% vs. 60-80% The shaded bands represent the uncertainties in the Glauber model calculations for $\langle N_{bin}\rangle$ and $\langle N_{part}\rangle$ [32]. Also shown are results for $\Lambda$ and $K^{0}_{S}$ [11] and protons and $\pi^{+}$ [33].

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

Minimum-bias (0–80% of the collision cross section) v2(pT ) for identified hadrons at |η| < 1 from Au+Au collisions at √sNN = 62.4 GeV. To facilitate comparisons between panels, v2 values for inclusive charged hadrons are displayed in each panel. The error bars on the data points represent statistical uncertainties. Systematic uncertainties for the identified particles are shown as shaded bands around v2 = 0.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The unidentified charged hadron, charged pion, $K^0_S$, charged kaon, proton and Λ+Λ v2 as a function of pT for 10%–40%, 0%–10% and 40%–80% of the Au+Au interaction cross section at √sNN = 62.4 GeV. Weak-decay feed-down errors are included in the error bars on the data points while non-flow and tracking error uncertainties are plotted as bands around v2 = 0, which apply to all identified particles. The errors are asymmetric and the portion of the error band above (below) zero represents the negative (positive) error.

The ratio of Λ v2 to Λbar v2. The data are from minimum bias Au+Au collisions at √sNN = 62.4 and 200 GeV. The bands show the average values of the ratios within the indicated pT ranges.

The ratio of Λ v2 to Λbar v2. The data are from minimum bias Au+Au collisions at √sNN = 62.4 and 200 GeV. The bands show the average values of the ratios within the indicated pT ranges.

The pT integrated ratio of Λ v2 to Λbar v2 for three centrality intervals 0%–10%, 10%–40%, and 40%–80%. The data are from Au+Au collisions at √sNN = 62.4 and 200 GeV.

The pT integrated ratio of Λ v2 to Λbar v2 for three centrality intervals 0%–10%, 10%–40%, and 40%–80%. The data are from Au+Au collisions at √sNN = 62.4 and 200 GeV.

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

Identified particle v2 from minimum bias collisions at √sNN = 62.4 GeV scaled by the number of valence quarks in the hadron (nq) and plotted versus pT/nq (a) and (mT − m0)/nq (b). In each case a polynomial curve is fit to all particles except pions. The ratio of v2/nq to the fit function is shown in the bottom panels (c) and (d).

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

v2/nq scaled by the mean eccentricity of the initial overlap region versus (mT − m0)/nq for 0%–10%, 10%–40%, and 40%–80% most central Au+Au collisions at √sNN = 62.4 GeV.

Top panels - minimum bias $v_4$ for pions, charged kaons, $K^0_S$, anti-protons and Λ + Λbar at √sNN = 62.4 GeV. In the left panel the solid (dashed) line shows thevalue for $v_2^2$ for pions (kaons). Intheright panel the dashedline is $v_2^2$ for Λ + Λbar. Bottom panels - $v_4$ scaled by $v_2^2$ (points where $v_4$ and $v_2$ fluctuate around zero are not plotted). Grey bands correspond to the fit results described in the text and Table II. The systematic errors on the $v_4/v_2^2$ ratio from nonflow are included in the error bars leading to asymmetric errors.

Top panels - minimum bias $v_4$ for pions, charged kaons, $K^0_S$, anti-protons and Λ + Λbar at √sNN = 62.4 GeV. In the left panel the solid (dashed) line shows thevalue for $v_2^2$ for pions (kaons). Intheright panel the dashedline is $v_2^2$ for Λ + Λbar. Bottom panels - $v_4$ scaled by $v_2^2$ (points where $v_4$ and $v_2$ fluctuate around zero are not plotted). Grey bands correspond to the fit results described in the text and Table II. The systematic errors on the $v_4/v_2^2$ ratio from nonflow are included in the error bars leading to asymmetric errors.

Top panels - minimum bias $v_4$ for pions, charged kaons, $K^0_S$, anti-protons and Λ + Λbar at √sNN = 62.4 GeV. In the left panel the solid (dashed) line shows thevalue for $v_2^2$ for pions (kaons). Intheright panel the dashedline is $v_2^2$ for Λ + Λbar. Bottom panels - $v_4$ scaled by $v_2^2$ (points where $v_4$ and $v_2$ fluctuate around zero are not plotted). Grey bands correspond to the fit results described in the text and Table II. The systematic errors on the $v_4/v_2^2$ ratio from nonflow are included in the error bars leading to asymmetric errors.

Top panels - minimum bias $v_4$ for pions, charged kaons, $K^0_S$, anti-protons and Λ + Λbar at √sNN = 62.4 GeV. In the left panel the solid (dashed) line shows thevalue for $v_2^2$ for pions (kaons). Intheright panel the dashedline is $v_2^2$ for Λ + Λbar. Bottom panels - $v_4$ scaled by $v_2^2$ (points where $v_4$ and $v_2$ fluctuate around zero are not plotted). Grey bands correspond to the fit results described in the text and Table II. The systematic errors on the $v_4/v_2^2$ ratio from nonflow are included in the error bars leading to asymmetric errors.

Top panels - minimum bias $v_4$ for pions, charged kaons, $K^0_S$, anti-protons and Λ + Λbar at √sNN = 62.4 GeV. In the left panel the solid (dashed) line shows thevalue for $v_2^2$ for pions (kaons). Intheright panel the dashedline is $v_2^2$ for Λ + Λbar. Bottom panels - $v_4$ scaled by $v_2^2$ (points where $v_4$ and $v_2$ fluctuate around zero are not plotted). Grey bands correspond to the fit results described in the text and Table II. The systematic errors on the $v_4/v_2^2$ ratio from nonflow are included in the error bars leading to asymmetric errors.

Top panels - minimum bias $v_4$ for pions, charged kaons, $K^0_S$, anti-protons and Λ + Λbar at √sNN = 62.4 GeV. In the left panel the solid (dashed) line shows thevalue for $v_2^2$ for pions (kaons). Intheright panel the dashedline is $v_2^2$ for Λ + Λbar. Bottom panels - $v_4$ scaled by $v_2^2$ (points where $v_4$ and $v_2$ fluctuate around zero are not plotted). Grey bands correspond to the fit results described in the text and Table II. The systematic errors on the $v_4/v_2^2$ ratio from nonflow are included in the error bars leading to asymmetric errors.

Top panels - minimum bias $v_4$ for pions, charged kaons, $K^0_S$, anti-protons and Λ + Λbar at √sNN = 62.4 GeV. In the left panel the solid (dashed) line shows thevalue for $v_2^2$ for pions (kaons). Intheright panel the dashedline is $v_2^2$ for Λ + Λbar. Bottom panels - $v_4$ scaled by $v_2^2$ (points where $v_4$ and $v_2$ fluctuate around zero are not plotted). Grey bands correspond to the fit results described in the text and Table II. The systematic errors on the $v_4/v_2^2$ ratio from nonflow are included in the error bars leading to asymmetric errors.

Top panels - minimum bias $v_4$ for pions, charged kaons, $K^0_S$, anti-protons and Λ + Λbar at √sNN = 62.4 GeV. In the left panel the solid (dashed) line shows thevalue for $v_2^2$ for pions (kaons). Intheright panel the dashedline is $v_2^2$ for Λ + Λbar. Bottom panels - $v_4$ scaled by $v_2^2$ (points where $v_4$ and $v_2$ fluctuate around zero are not plotted). Grey bands correspond to the fit results described in the text and Table II. The systematic errors on the $v_4/v_2^2$ ratio from nonflow are included in the error bars leading to asymmetric errors.

Top panel - $v_2$ for pions and protons at √sNN = 62.4 and 17.3 GeV. The 62.4 GeV data are from TOF and dE/dx measurements combined. Middle and bottom panel - ratios of $v_2$ for $\pi^++\pi^-$, $K^0_S$, p+p, Λ+Λbar and at different center-of-mass energies scaled by the values at 62.4 GeV. The grey and yellow bands represent systematic uncertainties in the v2 ratios arising from non-flow effects. The grey bands (above unity) are the uncertainties for the 200 GeV/62.4 GeV data and the yellow bands (below unity) are for the 17.3 GeV/62.4 GeV data.

Top panel - $v_2$ for pions and protons at √sNN = 62.4 and 17.3 GeV. The 62.4 GeV data are from TOF and dE/dx measurements combined. Middle and bottom panel - ratios of $v_2$ for $\pi^++\pi^-$, $K^0_S$, p+p, Λ+Λbar and at different center-of-mass energies scaled by the values at 62.4 GeV. The grey and yellow bands represent systematic uncertainties in the v2 ratios arising from non-flow effects. The grey bands (above unity) are the uncertainties for the 200 GeV/62.4 GeV data and the yellow bands (below unity) are for the 17.3 GeV/62.4 GeV data.

Top panel - $v_2$ for pions and protons at √sNN = 62.4 and 17.3 GeV. The 62.4 GeV data are from TOF and dE/dx measurements combined. Middle and bottom panel - ratios of $v_2$ for $\pi^++\pi^-$, $K^0_S$, p+p, Λ+Λbar and at different center-of-mass energies scaled by the values at 62.4 GeV. The grey and yellow bands represent systematic uncertainties in the v2 ratios arising from non-flow effects. The grey bands (above unity) are the uncertainties for the 200 GeV/62.4 GeV data and the yellow bands (below unity) are for the 17.3 GeV/62.4 GeV data.

Top panel - $v_2$ for pions and protons at √sNN = 62.4 and 17.3 GeV. The 62.4 GeV data are from TOF and dE/dx measurements combined. Middle and bottom panel - ratios of $v_2$ for $\pi^++\pi^-$, $K^0_S$, p+p, Λ+Λbar and at different center-of-mass energies scaled by the values at 62.4 GeV. The grey and yellow bands represent systematic uncertainties in the v2 ratios arising from non-flow effects. The grey bands (above unity) are the uncertainties for the 200 GeV/62.4 GeV data and the yellow bands (below unity) are for the 17.3 GeV/62.4 GeV data.

Top panel - $v_2$ for pions and protons at √sNN = 62.4 and 17.3 GeV. The 62.4 GeV data are from TOF and dE/dx measurements combined. Middle and bottom panel - ratios of $v_2$ for $\pi^++\pi^-$, $K^0_S$, p+p, Λ+Λbar and at different center-of-mass energies scaled by the values at 62.4 GeV. The grey and yellow bands represent systematic uncertainties in the v2 ratios arising from non-flow effects. The grey bands (above unity) are the uncertainties for the 200 GeV/62.4 GeV data and the yellow bands (below unity) are for the 17.3 GeV/62.4 GeV data.

Top panel - $v_2$ for pions and protons at √sNN = 62.4 and 17.3 GeV. The 62.4 GeV data are from TOF and dE/dx measurements combined. Middle and bottom panel - ratios of $v_2$ for $\pi^++\pi^-$, $K^0_S$, p+p, Λ+Λbar and at different center-of-mass energies scaled by the values at 62.4 GeV. The grey and yellow bands represent systematic uncertainties in the v2 ratios arising from non-flow effects. The grey bands (above unity) are the uncertainties for the 200 GeV/62.4 GeV data and the yellow bands (below unity) are for the 17.3 GeV/62.4 GeV data.

Top panel - $v_2$ for pions and protons at √sNN = 62.4 and 17.3 GeV. The 62.4 GeV data are from TOF and dE/dx measurements combined. Middle and bottom panel - ratios of $v_2$ for $\pi^++\pi^-$, $K^0_S$, p+p, Λ+Λbar and at different center-of-mass energies scaled by the values at 62.4 GeV. The grey and yellow bands represent systematic uncertainties in the v2 ratios arising from non-flow effects. The grey bands (above unity) are the uncertainties for the 200 GeV/62.4 GeV data and the yellow bands (below unity) are for the 17.3 GeV/62.4 GeV data.

Top panel - $v_2$ for pions and protons at √sNN = 62.4 and 17.3 GeV. The 62.4 GeV data are from TOF and dE/dx measurements combined. Middle and bottom panel - ratios of $v_2$ for $\pi^++\pi^-$, $K^0_S$, p+p, Λ+Λbar and at different center-of-mass energies scaled by the values at 62.4 GeV. The grey and yellow bands represent systematic uncertainties in the v2 ratios arising from non-flow effects. The grey bands (above unity) are the uncertainties for the 200 GeV/62.4 GeV data and the yellow bands (below unity) are for the 17.3 GeV/62.4 GeV data.

High transverse-momentum spectra ($p_{T} > 2.5$ GeV/c) of charged pions, protons, and antiprotons for the rapidity regions $|y| < 0.5$ (solid symbols) and $0.5 < |y| < 1.0$ (open symbols) for $d+Au$ collisions and various event centrality classes at $\sqrt{s_{NN}}=200$ GeV.

High transverse-momentum spectra ($p_{T} > 2.5$ GeV/c) of charged pions, protons, and antiprotons for the rapidity regions $|y| < 0.5$ (solid symbols) and $0.5 < |y| < 1.0$ (open symbols) for $d+Au$ collisions and various event centrality classes at $\sqrt{s_{NN}}=200$ GeV.

High transverse-momentum spectra ($p_{T} > 2.5$ GeV/c) of charged pions, protons, and antiprotons for the rapidity regions $|y| < 0.5$ (solid symbols) and $0.5 < |y| < 1.0$ (open symbols) for $d+Au$ collisions and various event centrality classes at $\sqrt{s_{NN}}=200$ GeV.

High transverse-momentum spectra ($p_{T} > 2.5$ GeV/c) of charged pions, protons, and antiprotons for the rapidity regions $|y| < 0.5$ (solid symbols) and $0.5 < |y| < 1.0$ (open symbols) for $d+Au$ collisions and various event centrality classes at $\sqrt{s_{NN}}=200$ GeV.

High transverse-momentum spectra ($p_{T} > 2.5$ GeV/c) of charged pions, protons, and antiprotons for the rapidity regions $|y| < 0.5$ (solid symbols) and $0.5 < |y| < 1.0$ (open symbols) for $d+Au$ collisions and various event centrality classes at $\sqrt{s_{NN}}=200$ GeV.

High transverse-momentum rapidity asymmetry factor ($Y_{\scriptsize{\mbox{Asym}}}$) for $\pi^{+}+\pi^{-}$ and $p+\bar{p}$ for $|y| < 0.5$ and $0.5 < |y| < 1.0$ for minimum bias $d+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV. For comparison the inclusive charged hadron results from STAR [5] are also shown by the curves.

Centrality dependence of high transverse-momentum rapidity asymmetry factor ($Y_{\scriptsize{\mbox{Asym}}}$) for $\pi^{+}+\pi^{-}$ and $p+\bar{p}$ for $|y| < 0.5$ (left panels) and $0.5 < |y| < 1.0$ (right panels) for $d+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV.

Centrality dependence of high transverse-momentum rapidity asymmetry factor ($Y_{\scriptsize{\mbox{Asym}}}$) for $\pi^{+}+\pi^{-}$ and $p+\bar{p}$ for $|y| < 0.5$ (left panels) and $0.5 < |y| < 1.0$ (right panels) for $d+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV.

Variation of nuclear modification factor ($R_{CP}$) for $\pi^{+}+\pi^{-}$ and $p+\bar{p}$ with rapidity for $p_{T} > 2.5$ GeV/c for $d+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV. Also shown for comparison are the $R_{CP}$ values for inclusive charged hadrons as a function of pseudorapidity [5]. The errors shown as boxes are the systematic errors. The error due to number of binary collisions is $\sim 14\%$ and is not shown in the figure.

Variation of $\pi^{-} / \pi^{+}$ and $\bar{p} / p$ with rapidity for $2.5 < p_{T} < 10$ GeV/c for peripheral ($40-100\%$) and n-tag events for $d+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV. Also shown for comparison are the $\bar{p} / p$ for minimum bias $p+p$ collisions. The boxes are the systematic errors. The ratios for n-tag events are shifted by 0.05 units in rapidity and those for $p+p$ collisions by -0.05 units in rapidity for clarity of presentation.

Variation of $\pi^{-} / \pi^{+}$ and $\bar{p} / p$ with rapidity for $2.5 < p_{T} < 10$ GeV/c for peripheral ($40-100\%$) and n-tag events for $d+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV. Also shown for comparison are the $\bar{p} / p$ for minimum bias $p+p$ collisions. The boxes are the systematic errors. The ratios for n-tag events are shifted by 0.05 units in rapidity and those for $p+p$ collisions by -0.05 units in rapidity for clarity of presentation.

Variation of double ratio $(\pi^{-} / \pi^{+})_{\scriptsize{\mbox{n-tag}}}/(\pi^{-} / \pi^{+})_{40-100\%}$ and $(\bar{p} / p)_{\scriptsize{\mbox{n-tag}}}/(\bar{p} / p)_{40-100\%}$ with rapidity for $2.5 < p_{T} < 10$ GeV/c for $d+Au$ collisions at $\sqrt{s_{NN}}=200$ GeV. The boxes are the systematic errors.

K0sK0s correlation function

K0sK0s theoretical correlation function

K0sK0s correlation function with fit functions

K0sK0s radius versus pair purity

K0sK0s lambda versus pair purity

K0sK0s radius versus Mean KT

K0sK0s radius versus MT

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

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

Corrected mid-rapidity (|y| < 0.5) pT spectra for $K^{+}$, $K^{−}$, $K^{0}_{S}$, Λ, Ξ, and Ω. Λ spectra that have been corrected for feed-down are shown as open symbols in the Λ panel. The dashed lines are fits using Equation 11 except for the $\Omega+\overline{\Omega}$ where the fit uses Equation 9. The error bars displayed include systematic errors while the fits were done using statistical errors only for all species except the charged kaons.

Mid-rapidity (|y| < 0.5) ratios of $\overline{\Lambda}$ to $\Lambda$ and $\overline{\Xi}^{+}$ to $\Xi^{-}$ vs $p_{T}$. The dashed lines in 10(b) and 10(c) are the errorweighted means over the measured $p_{T}$ range, 0.882±0.017 for $\overline{\Lambda}/\Lambda$, 0.921±0.062 for $\overline{\Xi}/\Xi$. Figure 10(a) shows the $p_{T}$ averaged ratio for our measured species compared with measurements from Au+Au. The error bars are statistical only. The dashed line in 10(a) is at unity for reference.

Mid-rapidity (|y| < 0.5) ratios of $\overline{\Lambda}$ to $\Lambda$ and $\overline{\Xi}^{+}$ to $\Xi^{-}$ vs $p_{T}$. The dashed lines in 10(b) and 10(c) are the errorweighted means over the measured $p_{T}$ range, 0.882±0.017 for $\overline{\Lambda}/\Lambda$, 0.921±0.062 for $\overline{\Xi}/\Xi$. Figure 10(a) shows the $p_{T}$ averaged ratio for our measured species compared with measurements from Au+Au. The error bars are statistical only. The dashed line in 10(a) is at unity for reference.

Mid-rapidity (|y| < 0.5) ratios of $\overline{\Lambda}$ to $\Lambda$ and $\overline{\Xi}^{+}$ to $\Xi^{-}$ vs $p_{T}$. The dashed lines in 10(b) and 10(c) are the errorweighted means over the measured $p_{T}$ range, 0.882±0.017 for $\overline{\Lambda}/\Lambda$, 0.921±0.062 for $\overline{\Xi}/\Xi$. Figure 10(a) shows the $p_{T}$ averaged ratio for our measured species compared with measurements from Au+Au. The error bars are statistical only. The dashed line in 10(a) is at unity for reference.

⟨$p_{T}$⟩ vs particle mass for different particles measured by STAR. Error bars include systematic errors. The ISR parameterization is given in reference [32]

⟨$p_{T}$⟩ vs charged multiplicity for $K^{+}$, $K^{−}$, $K^{0}_{S}$, $\Lambda+\overline{\Lambda}$, and $\Xi+\overline{\Xi}$. The points for $\Xi+\overline{\Xi}$ have been determined using only the measured region. The error bars are statistical only. See text for more details.

Ratios of multiplicity binned spectra to minimum bias spectra ($R_{pp}$) for $K^{0}_{S}$ and Λ and the ratio of the Λ spectrum to the $K^{0}_{S}$ spectrum in each multiplicity bin. See text for further details.

Ratios of multiplicity binned spectra to minimum bias spectra ($R_{pp}$) for $K^{0}_{S}$ and Λ and the ratio of the Λ spectrum to the $K^{0}_{S}$ spectrum in each multiplicity bin. See text for further details.

Ratios of multiplicity binned spectra to minimum bias spectra ($R_{pp}$) for $K^{0}_{S}$ and Λ and the ratio of the Λ spectrum to the $K^{0}_{S}$ spectrum in each multiplicity bin. See text for further details.

Parameters of ratio data to statistical model fit using THERMUS. Filled circles are ratios from $\sqrt{s}$=200 GeV collisions in STAR. Solid lines are the results from the statis- tical model fit. All ratios to the left of the vertical line were used in the fit. The $(\Omega^{-}+\overline{\Omega^{+}})/\Xi^{-}$ ratio was then predicted from the fit results. The dashed lines in the lower panel are guides for the eye at 1 σ.

Parameters of ratio data to statistical model fit using THERMUS. Filled circles are ratios from $\sqrt{s}$=200 GeV collisions in STAR. Solid lines are the results from the statis- tical model fit. All ratios to the left of the vertical line were used in the fit. The $(\Omega^{-}+\overline{\Omega^{+}})/\Xi^{-}$ ratio was then predicted from the fit results. The dashed lines in the lower panel are guides for the eye at 1 σ.

FIG. 7. Left: Parameters $n$ and $n_{h} / \tilde{n}_{s}$ vs $\hat{n}_{\mathrm{ch}}$ from free $\chi^{2}$ fits (solid symbols) in Table I. The open symbols represent similar free fits but with $\bar{y}_{t}$ held fixed at $2.65$ (see text for discussion). The bands represent the corresponding two-component fixed parametrizations with their stated errors also given in Table $I$. Right: Ratio $n_{h}\left(\hat{n}_{\mathrm{ch}}\right) / n_{s}\left(\hat{n}_{\mathrm{ch}}\right)$ derived from integrals in Fig. 4 (left panel) (open squares) and from Gaussian amplitudes in Fig. 5 (left panel) (solid dots). The dashed and dotted lines have slopes $0.105$ and $0.095$, respectively. The solid curve is described in the text. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

FIG. 8. $\left\langle p_{t}\right\rangle\left(\hat{n}_{\mathrm{ch}}\right)$ derived from the two-component $H_{0}$ Gaussian amplitudes (solid dots), from the running integrals (open squares), and from power-law fits to STAR and UA1 data (open FIG. 8. $\left\langle p_{t}\right\rangle\left(\hat{n}_{\mathrm{ch}}\right)$ derived from the two-component $H_{0}$ Gaussian amplitudes (solid dots), from the running integrals (open squares), and from power-law fits to STAR and UA1 data (open circles, triangles). The solid and dotted lines correspond to Eq. (6) with $\alpha=0.0095$ and $0.015$, respectively.circles, triangles). The solid and dotted lines correspond to Eq. (6) with $\alpha=0.0095$ and $0.015$, respectively. NOTE 1: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty. NOTE 2: The "ext" uncertainty corresponds to the systematic uncertainty of the extrapolation of the particle yield below $p_T < 0.2 GeV/c$ (hatched region).

(a) The corrected integrated yield $dN/dy$ at mid-rapidity for $\overline{\Xi}^{+}$ , $\overline{\Lambda}$ and $\overline{p}$ divided by $N_{part}$, normalized to the most peripheral centrality interval (60 − 80%), plotted as a function of $N_{part}$. The gray, black and dashed bands represent the errors on the normalization to the most peripheral bin for the $\overline{\Xi}^{+}$, $\overline{\Lambda}$ and $\overline{p}$. Other errors shown are statistical only. (b) $\gamma_{s}$ as a function of $N_{part}$ calculated from thermal model fits to the measured particle yields ($\pi$,K,p [24], $\Lambda$, $\Xi$, $\Omega$ and their anti-particles) at 200 GeV. Values for $e^{+}+e^{−}$ and $p+\overline{p}$ collisions at $\sqrt{s_{NN}}$ = 91 and 200 GeV respectively and for Pb+Pb SPS collisions at $\sqrt{s_{NN}}$ = 17.2 GeV are shown for comparison [26, 27, 28].

(a) The corrected integrated yield $dN/dy$ at mid-rapidity for $\overline{\Xi}^{+}$ , $\overline{\Lambda}$ and $\overline{p}$ divided by $N_{part}$, normalized to the most peripheral centrality interval (60 − 80%), plotted as a function of $N_{part}$. The gray, black and dashed bands represent the errors on the normalization to the most peripheral bin for the $\overline{\Xi}^{+}$, $\overline{\Lambda}$ and $\overline{p}$. Other errors shown are statistical only. (b) $\gamma_{s}$ as a function of $N_{part}$ calculated from thermal model fits to the measured particle yields ($\pi$,K,p [24], $\Lambda$, $\Xi$, $\Omega$ and their anti-particles) at 200 GeV. Values for $e^{+}+e^{−}$ and $p+\overline{p}$ collisions at $\sqrt{s_{NN}}$ = 91 and 200 GeV respectively and for Pb+Pb SPS collisions at $\sqrt{s_{NN}}$ = 17.2 GeV are shown for comparison [26, 27, 28].

$R_{CP}$ for $\Xi^{-}+\overline{\Xi}^{+}$ and $\Omega^{-}+\overline{\Omega}^{+}$ at mid-rapidity (centrality interval : 0 − 5% vs. 40 − 60%). A dashed line for charged hadrons and gray band for $\Lambda^{-}+\overline{\Lambda}$ are shown as comparison. The gray rectangles represent participant and binary scalings.

Nuclear modification factors Rcp for $\pi^{+}$ + $\pi^{-}$ and p + $\bar{p}$ in 200 GeV AuAu collisions.

Nuclear modification factors Rcp for $\pi^{+}$ + $\pi^{-}$ and p + $\bar{p}$ in 200 GeV AuAu collisions.

Nuclear modification factors Rcp for $\pi^{+}$ + $\pi^{-}$ and p + $\bar{p}$ in 200 GeV AuAu collisions.

The $\pi^{-}$/$\pi{+}$ and $\bar{p}$/p ratios in 12 $\%$ central, MB Au+Au and d+Au collisions at $\sqrt{s}_{NN}$ = 200 GeV. Errors are statistical and point-to-point systematic.

The p/$\pi^{+}$ and $\bar{p}$/$\pi^{-}$ ratios from d+Au and Au+Au collisions at $\sqrt{s}_{NN}$ = 200 GeV. The systematic uncertainties for 60-80$\%$ Au+Au collisions are similar.

$\sqrt{s_{NN}}$ dependence of $<p_t>$ fluctuations for central collisions and STAR full acceptance. CERES fluctuation data (12.3 and 17.3 GeV) were linearly extrapolated to the STAR $\eta$ acceptance

Hijing-1.37 $<p_t>$ fluctuations for central collisions and the STAR acceptance at three energies, for quench-on, quench-off and jets-off collisions.

Trigger-normalized fragment distribution $D(z_T)$ with $8 < p_T^{trig} < 15$ GeV/c for near-side ($|\Delta\phi-\pi|$<0.63) correlations in Au+Au collisions at 200 GeV.

Trigger-normalized fragment distribution $D(z_T)$ with $8 < p_T^{trig} < 15$ GeV/c for away-side ($|\Delta\phi-\pi|$<0.63) correlations in d+Au and Au+Au collisions at 200 GeV.

Coincidence probability versus azimuthal angle difference between the forward $\pi^{0}$ and a leading charged particle at midrapidity with $p_{T}$ > 0.5 GeV/c. The left (right) column is p+p (d+Au) data. The curves are fits described in the text, including the area of the back-to-back peak (S).

Coincidence probability versus azimuthal angle difference between the forward $\pi^{0}$ and a leading charged particle at midrapidity with $p_{T}$ > 0.5 GeV/c. The left (right) column is p+p (d+Au) data. The curves are fits described in the text, including the area of the back-to-back peak (S).

Coincidence probability versus azimuthal angle difference between the forward $\pi^{0}$ and a leading charged particle at midrapidity with $p_{T}$ > 0.5 GeV/c. The left (right) column is p+p (d+Au) data. The curves are fits described in the text, including the area of the back-to-back peak (S).

Coincidence probability versus azimuthal angle difference between the forward $\pi^{0}$ and a leading charged particle at midrapidity with $p_{T}$ > 0.5 GeV/c. The left (right) column is p+p (d+Au) data. The curves are fits described in the text, including the area of the back-to-back peak (S).

Transverse momentum distribution for $\pi^+$ production in d+Au collisions with centrality 20-40% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for $\pi^+$ production in d+Au collisions with centrality 40-100% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for $\pi^-$ production in d+Au minbias events in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for $\pi^-$ production in p+p NSD events in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for $\pi^-$ production in d+Au collisions with centrality 0-20% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for $\pi^-$ production in d+Au collisions with centrality 20-40% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for $\pi^-$ production in d+Au collisions with centrality 40-100% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for proton production in d+Au minbias events in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for proton production in p+p NSD events in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for proton production in d+Au collisions with centrality 0-20% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for proton production in d+Au collisions with centrality 20-40% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for proton production in d+Au collisions with centrality 40-100% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for pbar production in d+Au minbias events in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for pbar production in p+p NSD events in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for pbar production in d+Au collisions with centrality 0-20% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for pbar production in d+Au collisions with centrality 20-40% in the mid rapidity region, $|y|<0.5$.

Transverse momentum distribution for pbar production in d+Au collisions with centrality 40-100% in the mid rapidity region, $|y|<0.5$.

Nuclear modification factors $R_{dAu}$ for $\pi^+$ + $\pi^-$ production in the mid rapidity region, $|y|<0.5$. There is an overall additional normalization uncertainty of the order of 17%. For the measurement of $R_{dAu}$ there is an additional uncertainty of 5% due to uncertainty on $N_{bin}$ for minbias collisions.

Nuclear modification factors $R_{dAu}$ for $p + \bar{p}$ production in the mid rapidity region, $|y|<0.5$. There is an overall additional normalization uncertainty of the order of 17%. For the measurement of $R_{dAu}$ there is an additional uncertainty of 5% due to uncertainty on $N_{bin}$ for minbias collisions.

Ratios of $ \pi^-/\pi^+ $ in DEUT AU and P P minbias collisions in the mid rapidity region, $|y|<0.5$. In addition there are approximately 7% correlated errors for the ratios in P P collisions and about 10% in the DEUT AU collisions.

Ratios of $ \bar{p}/p $ in DEUT AU and P P minbias collisions in the mid rapidity region, $|y|<0.5$. In addition there are approximately 7% correlated errors for the ratios in P P collisions and about 10% in the DEUT AU collisions.

Ratios of $ p/\pi^+ $ in DEUT AU and P P minbias collisions in the mid rapidity region, $|y|<0.5$. In addition there are approximately 7% correlated errors for the ratios in P P collisions and about 10% in the DEUT AU collisions.

Ratios of $ \bar{p}/\pi^- $ in DEUT AU and P P minbias collisions in the mid rapidity region, $|y|<0.5$. In addition there are approximately 7% correlated errors for the ratios in P P collisions and about 10% in the DEUT AU collisions.

(Color Online) Variation of $N_{ch}$ normalized to the number of collisions in the FTPC coverage $(2.9 \leq \eta \leq 3.9)$ and $N_{\gamma}$ normalized to number of collisions, in the PMD coverage $(2.3 \leq \eta \leq 3.7)$ as a function of $N_{coll}$. The lower band shows the uncertainty in the ratio due to uncertainties in $N_{coll}$ calculations.

(Color Online) $dN/d\eta$ for charged particles and photons for Au + Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV for various event centrality classes.

(Color Online) $dN/d\eta$ for charged particles and photons for Au + Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV for various event centrality classes.

(Color Online) Half width at half maximum of the pseudorapidity distributions ($\eta_{h}$) of charged particles as a function of total charged particle multiplicity ($N_{T}$) normalized to the center of mass energy. The Au + Au collision data are from the PHOBOS [8] experiment and p + p collision data are from the ISR [31] experiments.

(Color Online) Half width at half maximum of the pseudorapidity distributions ($\eta_{h}$) of charged particles as a function of total charged particle multiplicity ($N_{T}$) normalized to the center of mass energy. The Au + Au collision data are from the PHOBOS [8] experiment and p + p collision data are from the ISR [31] experiments.

(Color Online) Variation of $dN_{ch}/d\eta$ normalized to $N_{part}$ with $\eta – y_{beam}$ for central and peripheral collisions for positively charged hadrons ($h^{+}$) and negatively charged hadrons ($h^{−}$).

(Color Online) The top panel shows the variation of pion rapidity density normalized to $N_{part}$ with $y – y_{beam}$ for central collisions at various collision energies. Also shown is the estimated $dN_{\pi^{0}}/dy$ obtained from $dN\_{\gamma}/dy$ normalized to $N_{part}$. The bottom panel shows the variation of net proton rapidity density normalized to $N_{part}$ with $y – y_{beam}$ for central collisions at various collision energies.

(Color Online) The top panel shows the variation of pion rapidity density normalized to $N_{part}$ with $y – y_{beam}$ for central collisions at various collision energies. Also shown is the estimated $dN_{\pi^{0}}/dy$ obtained from $dN\_{\gamma}/dy$ normalized to $N_{part}$. The bottom panel shows the variation of net proton rapidity density normalized to $N_{part}$ with $y – y_{beam}$ for central collisions at various collision energies.

The purity and momentum-resolution corrected correlation functions $C_{true}(k^{*})$ for $p-\Lambda$, $\bar{p}-\bar{\Lambda}$ (a), $\bar{p}-\Lambda$, $p-\bar{\Lambda}$ (b). Curves correspond to fits done using the Lednicky and Lyuboshitz analytical model [12].

$(p-\Lambda) \oplus (\bar{p}-\bar{\Lambda})$ and $(\bar{p}-\Lambda) \oplus (p-\bar{\Lambda})$ combined correlation functions. Correlation functions are corrected for purity and momentum resolution. Curves correspond to fits done using the Lednicky and Lyuboshitz analytical model [12].

The combined $(\bar{p}-\Lambda) \oplus (p-\bar{\Lambda})$ spin-averaged $s$-wave scattering length compared with the previous measurements for the $p-\bar{p}$ system [24-27]. The curves show the one standard deviation contours. Note that for $(\bar{p}-\Lambda) \oplus (p-\bar{\Lambda})$ only, one should read $0.1973 \times \mbox{Im}f_{0}$ instead of $\mbox{Im}f_{0}$ on the x-axis.

The combined $(\bar{p}-\Lambda) \oplus (p-\bar{\Lambda})$ spin-averaged $s$-wave scattering length compared with the previous measurements for the $p-\bar{p}$ system [24-27]. The curves show the one standard deviation contours. Note that for $(\bar{p}-\Lambda) \oplus (p-\bar{\Lambda})$ only, one should read $0.1973 \times \mbox{Im}f_{0}$ instead of $\mbox{Im}f_{0}$ on the x-axis.

Pion source size [5] compared with proton and lambda source sizes as a function of the mean transverse mass $(\langle m_{t}\rangle)$ The curve shows the $\langle m_{t}\rangle^{-1/2}$ dependence with an arbitrary normalization.