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

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

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

Charged hadron azimuthal anisotropy v2 as a function of pT in 50-60% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 as a function of pT in 40-50% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 as a function of pT in 30-40% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 as a function of pT in 20-30% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 as a function of pT in 10-20% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 as a function of pT in 0-10% Cu+Cu collisions at 200 GeV using FTPC flow vectors, and that with subtracting the azimuthal correlations in p+p collisions.

Charged hadron azimuthal anisotropy v2 integrated over pT and eta vs centrality in Cu+Cu collisions at 200 GeV for the various methods.

Charged hadron azimuthal anisotropy v2 integrated over pT and eta vs centrality in Cu+Cu collisions at 62.4 GeV for the various methods.

Charged hadron azimuthal anisotropy v2 as a function of pT in 50-60% Cu+Cu collisions at 62.4 GeV. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 40-50% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 30-40% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 20-30% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 10-20% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 0-10% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.4.

Charged hadron azimuthal anisotropy v2 as a function of pT in 0-60% Cu+Cu collisions at 62.4 GeV using FTPC flow vectors. Data tables for 200 GeV are shown in Fig.3.

K0S azimuthal anisotropy v2 vs pT in 0-60% central Cu+Cu collisions at 200 GeV.

K0S azimuthal anisotropy v2 vs pT in 0-20% central Cu+Cu collisions at 200 GeV.

K0S azimuthal anisotropy v2 vs pT in 20-60% central Cu+Cu collisions at 200 GeV.

Lambda azimuthal anisotropy v2 vs pT in 0-60% central Cu+Cu collisions at 200 GeV.

Lambda azimuthal anisotropy v2 vs pT in 0-20% central Cu+Cu collisions at 200 GeV.

Lambda azimuthal anisotropy v2 vs pT in 20-60% central Cu+Cu collisions at 200 GeV.

Xi azimuthal anisotropy v2 vs pT in 0-60% central Cu+Cu collisions at 200 GeV.

phi azimuthal anisotropy v2 vs pT in 0-60% central Cu+Cu collisions at 200 GeV.

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.

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.

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

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

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