No description provided.

No description provided.

No description provided.

No description provided.

Data requested from authors.

No description provided.

No description provided.

No description provided.

Data supplied by S.Banerjee.

Data supplied by S.Banerjee.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

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

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 0.15 GeV, for centrality 30-50%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 0.15 GeV, for centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 5 GeV, for centrality 0-10%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 5 GeV, for centrality 10-30%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 5 GeV, for centrality 30-50%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 5 GeV, for centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 10 GeV, for centrality 0-10%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 10 GeV, for centrality 10-30%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 10 GeV, for centrality 30-50%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 10 GeV, for centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Ratio of charged jet spectra with leading track selection of pT > 0.15 GeV to pT > 5 GeV using two cone radius parameters R = 0.2 and 0.3, for centrality 0-10%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Ratio of charged jet spectra with leading track selection of pT > 0.15 GeV to pT > 5 GeV using two cone radius parameters R = 0.2 and 0.3, for centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Ratio of charged jet spectra with leading track selection of pT > 10 GeV to pT > 5 GeV using two cone radius parameters R = 0.2 and 0.3, for centrality 0-10%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Ratio of charged jet spectra with leading track selection of pT > 10 GeV to pT > 5 GeV using two cone radius parameters R = 0.2 and 0.3, for centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 0.15 GeV. Central collisions are defined to have centrality 0-10% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 0.15 GeV. Central collisions are defined to have centrality 10-30% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 0.15 GeV. Central collisions are defined to have centrality 30-50% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 5 GeV. Central collisions are defined to have centrality 0-10% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 5 GeV. Central collisions are defined to have centrality 10-30% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 5 GeV. Central collisions are defined to have centrality 30-50% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 10 GeV. Central collisions are defined to have centrality 0-10% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 10 GeV. Central collisions are defined to have centrality 10-30% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 10 GeV. Central collisions are defined to have centrality 30-50% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, as a function of the average number of participants in the collision (<Npart>), for charged jets with pT = 60-70 GeV and either R = 0.2 or R = 0.3 and a leading charged particle with pT > 0.15 GeV. The correspondence between <Npart> and the centrality classes is given in Table 1 of the paper, and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, as a function of the average number of participants in the collision (<Npart>), for charged jets with pT = 60-70 GeV and either R = 0.2 or R = 0.3 and a leading charged particle with pT > 5 GeV. The correspondence between <Npart> and the centrality classes is given in Table 1 of the paper, and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, as a function of the average number of participants in the collision (<Npart>), for charged jets with pT = 60-70 GeV and either R = 0.2 or R = 0.3 and a leading charged particle with pT > 10 GeV. The correspondence between <Npart> and the centrality classes is given in Table 1 of the paper, and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Ratio of charged jet pT-spectra with radius parameter R = 0.2 and 0.3 and a leading charged particle with pT > 5 GeV, for two centrality classes 0-10% and 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

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.

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.

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

Directed flow slope at mid-rapidity, $dv_{1}/dy|_{y=0}$, as a function of beam energy in 10-40\% Au+Au collisions. Dots represent deuterons. Open squares are the published results in~\cite{Adamczyk_2018} for protons. The band denotes the results for deuterons from AMPT transport model. Statistical uncertainties (bars) and systematic uncertainties (horizontal lines) are shown separately. For visibility, the data points are staggered horizontally.

The $p_{T}$ dependence of $v_1/A$ in $|y|<0.6$ for protons (open squares) in 10-40\% Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 14.5, 19.6, 27, and 39 GeV. Statistical uncertainties (bars) and systematic uncertainties (horizontal lines) are shown separately.

The $d_{T}$ dependence of $v_1/A$ in $|y|<0.6$ for deuterons (solid circles) in 10-40\% Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 14.5, 19.6, 27, and 39 GeV. Statistical uncertainties (bars) and systematic uncertainties (horizontal lines) are shown separately.

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.

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.

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.

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