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.

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.

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

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.

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.

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

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

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.

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.

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.

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.

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

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.

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.

Non photonic electron yield in minimum bias D+AU collisions versus PT 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 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 10-40%.

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

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.

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