$p_{T}^{Part}$ vs. $p_{T}^{Det}$ for SubLeading jets with $R = 0.4$.

$p_{T}^{Part}$ vs. $p_{T}^{Det}$ for SubLeading jets with $R = 0.2$.

$A_{J}$ distributions for Au+Au data (filled symbols) and p+p HT $\oplus$ Au+Au MB (open symbols) for low constituent $p_{T}^{Cut}$ di-jets from Fig. 1 compared to $A_{J}$ distributions calculated assuming the RC and EC null hypotheses, respectively, shown as colored bands; see the text for details.

$dN/d\delta p_{T}$ vs. $\delta p_{T}$ for $p_{T}^{Cut} > 2$ GeV/$c$ jets with $R = 0.4$.

$dN/d\delta p_{T}$ vs. $\delta p_{T}$ for $p_{T}^{Cut} > 0.2$ GeV/$c$ matched jets with $R = 0.4$.

$dN/d\delta p_{T}$ vs. $\delta p_{T}$ for $p_{T}^{Cut} > 2$ GeV/$c$ jets with $R = 0.2$.

$dN/d\delta p_{T}$ vs. $\delta p_{T}$ for $p_{T}^{Cut} > 0.2$ GeV/$c$ matched jets with $R = 0.2$.

Repetition of the analysis shown in Fig. 1 with a smaller resolution parameter $R = 0.2$. Normalized $A_{J}$ distributions for Au+Au HT data (filled symbols) and p+p HT $\oplus$ Au+Au MB (open symbols). The red circles are for jets found using only constituents with $p_{T}^{Cut} > 2$ GeV/$c$ and the black squares are for matched jets found using constituents with $p_{T}^{Cut} > 0.2$ GeV/$c$.

$R_{dAu}$ vs. $y$.

Comparison of our d+Au measurements to the pA measurements from E772. Ratio of $\Upsilon$ production in pA to pp scaled by mass number as a function of mass number. Shown are the 1S and 2S+3S $\Upsilon$ measurements from E772 and our 1S measurement.

Comparison of our d+Au measurements to the pA measurements from E772. Exponent $\alpha$ as a function of $x_{F}$.

Invariant mass distributions of electron pairs in the region $|y_{ee}| < 1.0$ for the centrality selections $30-60\%$.

Invariant mass distributions of electron pairs in the region $|y_{ee}| < 1.0$ for the centrality selections $10-30\%$.

Invariant mass distributions of electron pairs in the region $|y_{ee}| < 1.0$ for the centrality selections $0-10\%$.

Nuclear modification factor for $\Upsilon$(1S+2S+3S), in $|y| < 1.0$ in d+Au and Au+Au collisions as a function of $N_{part}$.

Nuclear modification factor for $\Upsilon$(1S+2S+3S), in $|y| < 0.5$, in d+Au and Au+Au collisions as a function of $N_{part}$.

Nuclear modification factor for $\Upsilon$(1S) in $|y| < 1.0$, in d+Au and Au+Au collisions as a function of $N_{part}$.

Nuclear modification factor for $\Upsilon$(1S) in $|y| < 0.5$, in d+Au and Au+Au collisions as a function of $N_{part}$.

Nuclear modification factor of quarkonium states as a function of binding energy as measured by STAR.

The same as the left but for Cu+Cu collisions. The systematic errors are shown as thin lines with wide caps at the ends and statistical errors are shown as thick lines with small caps at the end. Statistical and systematic errors are very small.

The difference of charge-independent (CI) v2{2} and like-sign (LS) $v_2\{2\}$ for Au+Au and Cu+Cu collisions at 200 (top panel) and 62.4 (bottom panel) GeV vs. the log of $\langle dN_{ch}/d\eta\rangle$.The statistical errors are smaller than the marker size and not visible for most of the data.

The difference of charge-independent (CI) v2{2} and like-sign (LS) $v_2\{2\}$ for Au+Au and Cu+Cu collisions at 200 (top panel) and 62.4 (bottom panel) GeV vs. the log of $\langle dN_{ch}/d\eta\rangle$.The statistical errors are smaller than the marker size and not visible for most of the data.

The difference of charge-independent (CI) v2{2} and like-sign (LS) $v_2\{2\}$ for Au+Au and Cu+Cu collisions at 200 (top panel) and 62.4 (bottom panel) GeV vs. the log of $\langle dN_{ch}/d\eta\rangle$.The statistical errors are smaller than the marker size and not visible for most of the data.

The difference of charge-independent (CI) v2{2} and like-sign (LS) $v_2\{2\}$ for Au+Au and Cu+Cu collisions at 200 (top panel) and 62.4 (bottom panel) GeV vs. the log of $\langle dN_{ch}/d\eta\rangle$.The statistical errors are smaller than the marker size and not visible for most of the data.

The LS and CI four-particle cumulant $v_2\{4\}^4$ for Au+Au collisions at 200 and 62.4 GeV for $0.15 < pT < 2.0$ GeV/$c$. The systematic errors are shown as narrow lines with wide caps at the end and statistical errors are shown as thick lines with narrow caps at the end. Statistical errors are not visible for most of the points.

The LS and CI four-particle cumulant $v_2\{4\}^4$ for Au+Au collisions at 200 and 62.4 GeV for $0.15 < pT < 2.0$ GeV/$c$. The systematic errors are shown as narrow lines with wide caps at the end and statistical errors are shown as thick lines with narrow caps at the end. Statistical errors are not visible for most of the points.

The LS and CI four-particle cumulant $v_2\{4\}^4$ for Cu+Cu collisions at 200 and 62.4 GeV for $0.15 < p_T < 2.0$ GeV/c. The most central points (two points for Cu+Cu 62.4 GeV) gives $v_2\{4\}^4 < 0$ for all the data sets. The negative values are probably due to large fluctuations in agreement with Eq. (1). These may include contributions from impact parameter spread and finite multiplicity bin width.

The LS and CI four-particle cumulant $v_2\{4\}^4$ for Cu+Cu collisions at 200 and 62.4 GeV for $0.15 < p_T < 2.0$ GeV/c. The most central points (two points for Cu+Cu 62.4 GeV) gives $v_2\{4\}^4 < 0$ for all the data sets. The negative values are probably due to large fluctuations in agreement with Eq. (1). These may include contributions from impact parameter spread and finite multiplicity bin width.

The difference of charge-independent (CI) $v_2\{4\}$ and like-sign (LS) $v_2\{4\}$ for Au+Au collisions at 200 and 62.4 GeV vs. the log of $\langle dN_{ch}/d\eta\rangle$.

The difference of charge-independent (CI) $v_2\{4\}$ and like-sign (LS) $v_2\{4\}$ for Au+Au collisions at 200 and 62.4 GeV vs. the log of $\langle dN_{ch}/d\eta\rangle$.

(Left) The difference between $v_2\{2\}^2$ and $v_2\{4\}^2$ for 200 GeV Au+Au and Cu+Cu collisions for both LS and CI combinations.

(Left) The difference between $v_2\{2\}^2$ and $v_2\{4\}^2$ for 200 GeV Au+Au and Cu+Cu collisions for both LS and CI combinations.

(Right) The difference between $v_2\{2\}^2$ and $v_2\{4\}^2$ for 62.4 GeV Au+Au and Cu+Cu collisions for both LS and CI combinations. The statistical and systematic errors are shown as in previous figures.

(Right) The difference between $v_2\{2\}^2$ and $v_2\{4\}^2$ for 62.4 GeV Au+Au and Cu+Cu collisions for both LS and CI combinations. The statistical and systematic errors are shown as in previous figures.

The upper limit on $\sigma_{v_2}/\langle v_2 \rangle$ for 200 GeV (left) and 62.4 GeV (right) Au+Au collisions from Eq. (9) compared to $\sigma_\varepsilon/\varepsilon$ from Eq. (10) for three different models. The upper limit is found using the LS results for $v_2\{2\}$. Data are from the range $0.15 < p_T < 2.0$ GeV/$c$. The shaded bands reflect the uncertainties on the models which are dominated by uncertainty on the distribution of nucleons inside the nucleus. The uncertainty is only shown for the MCG-N and fKLN-CGC models. The uncertainty on the MCG-Q model is the same as for the MCG-N model but is not shown for the visual clarity.

The upper limit on $\sigma_{v_2}/\langle v_2 \rangle$ for 200 GeV (left) and 62.4 GeV (right) Au+Au collisions from Eq. (9) compared to $\sigma_\varepsilon/\varepsilon$ from Eq. (10) for three different models. The upper limit is found using the LS results for $v_2\{2\}$. Data are from the range $0.15 < p_T < 2.0$ GeV/$c$. The shaded bands reflect the uncertainties on the models which are dominated by uncertainty on the distribution of nucleons inside the nucleus. The uncertainty is only shown for the MCG-N and fKLN-CGC models. The uncertainty on the MCG-Q model is the same as for the MCG-N model but is not shown for the visual clarity.

The upper limit on $\sigma_{v_2}/\langle v_2 \rangle$ for 200 GeV (left) and 62.4 GeV (right) Au+Au collisions from Eq. (9) compared to $\sigma_\varepsilon/\varepsilon$ from Eq. (10) for three different models. The upper limit is found using the LS results for $v_2\{2\}$. Data are from the range $0.15 < p_T < 2.0$ GeV/$c$. The shaded bands reflect the uncertainties on the models which are dominated by uncertainty on the distribution of nucleons inside the nucleus. The uncertainty is only shown for the MCG-N and fKLN-CGC models. The uncertainty on the MCG-Q model is the same as for the MCG-N model but is not shown for the visual clarity.

The upper limit on $\sigma_{v_2}/\langle v_2 \rangle$ for 200 GeV (left) and 62.4 GeV (right) Au+Au collisions from Eq. (9) compared to $\sigma_\varepsilon/\varepsilon$ from Eq. (10) for three different models. The upper limit is found using the LS results for $v_2\{2\}$. Data are from the range $0.15 < p_T < 2.0$ GeV/$c$. The shaded bands reflect the uncertainties on the models which are dominated by uncertainty on the distribution of nucleons inside the nucleus. The uncertainty is only shown for the MCG-N and fKLN-CGC models. The uncertainty on the MCG-Q model is the same as for the MCG-N model but is not shown for the visual clarity.

The STAR data compared to PHOBOS data [34] on $\sigma_{v_2}/\langle v_2 \rangle$ with $\delta_2$ for $\Delta\eta > 2$ taken to be zero (see Fig. 6 from Ref. [34]). The shaded band shows the errors quoted from Ref. [34].

The STAR data compared to PHOBOS data [34] on $\sigma_{v_2}/\langle v_2 \rangle$ with $\delta_2$ for $\Delta\eta > 2$ taken to be zero (see Fig. 6 from Ref. [34]). The shaded band shows the errors quoted from Ref. [34].

The upper limit on $\sigma_{v_2}/\langle v_2 \rangle$ for 200 GeV (left) and 62.4 GeV (right) Cu+Cu collisions from Eq. (9) compared to $\sigma_\varepsilon/\varepsilon$ from Eq. (10) for three different models.

The upper limit on $\sigma_{v_2}/\langle v_2 \rangle$ for 200 GeV (left) and 62.4 GeV (right) Cu+Cu collisions from Eq. (9) compared to $\sigma_\varepsilon/\varepsilon$ from Eq. (10) for three different models.

The upper limit on $\sigma_{v_2}/\langle v_2 \rangle$ for 200 GeV (left) and 62.4 GeV (right) Cu+Cu collisions from Eq. (9) compared to $\sigma_\varepsilon/\varepsilon$ from Eq. (10) for three different models.

The upper limit on $\sigma_{v_2}/\langle v_2 \rangle$ for 200 GeV (left) and 62.4 GeV (right) Cu+Cu collisions from Eq. (9) compared to $\sigma_\varepsilon/\varepsilon$ from Eq. (10) for three different models.

$\rho^0$ production cross section determined with the data set collected with trigger B, as a function of momentum transfer $t$, fitted with a double exponential fit function in the range $t < 0.1$. The statistical errors are shown. The ratio of incoherent to coherent cross sections is measured to be $0.20 \pm 0.8$. The fit parameters are shown in Table I.

Comparison of theoretical predictions to the measured total cross section for coherent $\rho^0$ production as a functionof $\sqrt{s_{NN}}$ . The measured cross section is based on the dataset collected with trigger B. The error bars show the sum of the statistical and systematic uncertainties. See text for details.

The correlation function C, C is plotted in units of (GeV/c)$^2$ and the relative azimuthal angle ∆φ in radians for 30-40% centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Relative statistical errors range from 0.8% in peripheral collisions to 0.9% in the most central collisions at the peak of the distribution.

The correlation function C, C is plotted in units of (GeV/c)$^2$ and the relative azimuthal angle ∆φ in radians for 0-5% centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV..Relative statistical errors range from 0.8% in peripheral collisions to 0.9% in the most central collisions at the peak of the distribution.

The correlation function C, C is plotted in units of (GeV/c)$^2$ and the relative azimuthal angle ∆φ in radians for 0-5% centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV..Relative statistical errors range from 0.8% in peripheral collisions to 0.9% in the most central collisions at the peak of the distribution.

Projection of the correlation function C, for|∆φ|<1.0 radians on the ∆η axis for 70-80% centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The correlation function C is plottedin units of (GeV/c)$^2$.The correlation function C is plotted in units of (GeV/c)$^2$. The solid line shows the fit obtained with Eq.2. The dotted line corresponds to the baseline, b, obtained in the fit and shaded band shows uncertainty in determining b.

Projection of the correlation function C, for|∆φ|<1.0 radians on the ∆η axis for 70-80% centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The correlation function C is plottedin units of (GeV/c)$^2$.The correlation function C is plotted in units of (GeV/c)$^2$. The solid line shows the fit obtained with Eq.2. The dotted line corresponds to the baseline, b, obtained in the fit and shaded band shows uncertainty in determining b.

Projection of the correlation function C, for|∆φ|<1.0 radians on the ∆η axis for 30-40% centrality. Statistical errors at(\Deltaeta_1,Deltaeta_2~(0,0) are approximately 0.084 for Au+Au. The correlation function C is plotted in units of (GeV/c)$^2$. The dotted line corresponds to the baseline,b, obtained in the fit and shaded band shows uncertainty in determining b.

Projection of the correlation function C, for|∆φ|<1.0 radians on the ∆η axis for 30-40% centrality. Statistical errors at(\Deltaeta_1,Deltaeta_2~(0,0) are approximately 0.084 for Au+Au. The correlation function C is plotted in units of (GeV/c)$^2$. The dotted line corresponds to the baseline,b, obtained in the fit and shaded band shows uncertainty in determining b.

Projection of the correlation function C, for|∆φ|<1.0 radians on the ∆η axis for 0-5% centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The correlation functionCis plotted in units of (GeV/c)$^2$. The correlation function C is plotted in units of (GeV/c)$^2$. The solid line shows the fit obtained with Eq.2. The dotted line corresponds to the baseline,b, obtained in the fit and shaded band shows uncertainty in determining b.

Projection of the correlation function C, for|∆φ|<1.0 radians on the ∆η axis for 0-5% centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The correlation functionCis plotted in units of (GeV/c)$^2$. The correlation function C is plotted in units of (GeV/c)$^2$. The solid line shows the fit obtained with Eq.2. The dotted line corresponds to the baseline,b, obtained in the fit and shaded band shows uncertainty in determining b.

RMS as function of the number of participating nucleons for the correlation function C, for nine centrality classes in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The dotted line represents a lower limit estimate of the RMS explained in the text and the shaded band represents systematic uncertainties on the RMS.

RMS as function of the number of participating nucleons for the correlation function C, for nine centrality classes in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The dotted line represents a lower limit estimate of the RMS explained in the text and the shaded band represents systematic uncertainties on the RMS.

The simulated electron pair invariant mass distributions for electrons at $8 < p_{T} < 10$ GeV/c

The distribution of the minimum distance $\Delta$Z (cm) between an electron track projection point at the Barrel Electromagnetic Calorimeter (BEMC) and all BEMC clusters along the beam direction from unlike-sign electron candidate pairs, like-sign electron candidate pairs and unlike-minus-like.

The momentum over energy $p/E_0$ distribution ($p$ - momentum, $E_0$ - the energy of the most energetic tower in a Barrel Electromagnetic Calorimeter cluster) from unlike-sign electron candidate pairs, like-sign electron candidate pairs and unlike-minus-like.

Electron identification efficiency of the cuts on number of TPC points, n$\sigma_e$ and Barrel Electromagnetic Calorimeter cuts in the Run2008 analysis. The total efficiency is the product of all individual ones.

Electron identification efficiency of the cuts on number of TPC points, n$\sigma_e$ and Barrel Electromagnetic Calorimeter cuts in the Run2005 analysis. The total efficiency is the product of all individual ones.

n$\sigma_e$ distribution in the Run2008 analysis for unlike-sign, like-sign and unlike-minus-like pairs at $2.5 < p_T < 3.0$ GeV after applying all the electron identification cuts except the n$\sigma_e$ cut.

n$\sigma_e$ distribution in the Run2008 analysis for unlike-sign, like-sign and unlike-minus-like pairs at $8 < p_T < 10$ GeV after applying all the electron identification cuts except the n$\sigma_e$ cut.

The mean and width of the Gaussian fitting functions for pure photonic electron (unlike - minus-like) n$\sigma_e$ distribution as shown in Fig. 4(a) and Fig 4(b) for each $p_T$ bin.

n$\sigma_e$ distribution for inclusive electrons at (a) $2.5 < p_T < 3.0$ GeV/c in the Run2008 analysis, (b) $2.5< p_T < 3.5$ GeV/c in the Run2005 analysis and (c) $8.0 < p_T < 10.0$ GeV/c in the Run2008 analysis after applying all electron identification cuts except the n$\sigma_e$ cut.

Purity of the inclusive electron sample as a function of $p_T$ in data from Run 2008.

Purity of the inclusive electron sample as a function of $p_T$ in data from Run 2005 with High Tower HT1 trigger.

Purity of the inclusive electron sample as a function of $p_T$ in data from Run 2005 with High Tower HT2 trigger.

Derived $p_T$ spectrum for inclusive photons and the uncertainty represented by the region between the spectra of $\pi^0$ and $\eta$ decay photons and inclusive photon with doubled direct photon yield.

Ratio of photon $p_T$ spectra from Fig. 7(a) to inclusive photons $p_T$ spectrum.

Photonic electron reconstruction efficiency as a function of $p_T$ for $\gamma$ conversion, $\pi^0$ and $\eta$ Dalitz decay for the Run2008 analysis, and for their combination for the Run2008 analysis.

Photonic electron reconstruction efficiency as a function of $p_T$ for $\gamma$ conversion for the Run2005 analysis.

The adc0 distribution for high-tower trigger events. The adc0 is the offline ADC value of a BEMC cluster’s most energetic tower which is one of the high-towers responsible for firing a high-tower trigger.

The $\chi^2$ as a function of the adc0 cut.

Raw inclusive electron $p_T$ spectrum from VPD trigger in Run2008 d+Au collisions before (open squares) and after applying the adc0 > 193 cut.

The adc0 distribution for data at $p_T$ = 4.0 − 5.0 GeV/c

The adc0 distribution for simulations at $p_T$ = 4.0 − 5.0 GeV/c

The $p_T$ dependence of high-tower trigger efficiency from VPD data for Run2008 analysis.

The $p_T$ dependence of high-tower trigger efficiency from simulation for Run2008 analysis.

The $p_T$ dependence of high-tower trigger efficiency for combined VPD data and simulation for Run2008 analysis.

Raw inclusive electron $p_T$ spectrum for minimum-bias (MB) and two high-tower triggers, i.e. HT1 and HT2 for Run2005 analysis.

Raw inclusive electron $p_T$ spectrum for high-tower trigge HT1 for Run2005 analysis.

Raw inclusive electron $p_T$ spectrum for high-tower trigge HT2 for Run2005 analysis.

The $p_T$ dependence of trigger efficiency for the high-tower trigger HT1 data for Run2005.

The $p_T$ dependence of trigger efficiency for the high-tower trigger HT2 data for Run2005.

Distribution of the calculated BBC cross section $\sigma_{BBC}$ before and after removing events at the beginning and the end of Run2008.

Invariant cross section of the electron from decays of $J/\psi$.

Invariant cross section of the electron from decays of $\varUpsilon$.

Invariant cross section of the electron from Drell-Yan decays.

Invariant cross section of the electron from decays of light vector mesons ($\rho,\omega, \phi$).

The uncertainty of the $J/\psi$ feeddown.

Ratio of non-photonic to photonic electron yield from the Run2008 analysis.

Ratio of non-photonic to photonic electron yield from the Run2005 analysis.

Invariant cross section for non-photonic electron production $(e^+ + e^-)/2$ in p+p collisions from the Run2008 analysis.

Invariant cross section for non-photonic electron production $(e^+ + e^-)/2$ in p+p collisions from the Run2005 analysis.

Invariant cross section for non-photonic electron production $(e^+ + e^-)/2$ in p+p collisions from the combined Run2005 and 2008 analyses.

Invariant cross section for non-photonic electron production $(e^+ + e^-)/2$ in p+p collisions (update on PRL 98 (2007) 192301).

Ratio of invariant cross section for non-photonic electron production $(e^+ + e^-)/2$ in p+p collisions from the combined Run2005 and 2008 analyses over FONLL.

Ratio of invariant cross section for non-photonic electron production $(e^+ + e^-)/2$ in p+p collisions (update on PRL 98 (2007) 192301) over FONLL.

Invariant cross section for electrons $(e^+ + e^-)/2$ from bottom meson decay, together with the ratio of the measurements to the FONLL predictions.

Invariant cross section for electrons $(e^+ + e^-)/2$ from charm meson decay, together with the ratio of the measurements to the FONLL predictions.

The 3ΛH (solid squares) and Λ (open circles) yield distributions versus cτ. The solid lines represent the cτ fits. The inset depicts the $\chi^2$ distribution of the best 3ΛH cτ fit. The error bars represent the statistical uncertaintiesonly.

Particle ratios as a function of center-of-mass energy per nucleon-nucleon collision. The 4ΛH/(4He × Λ/p) ratio is corrected for the spin degeneracy factor (38). The error bars represent statistical uncertainties only.

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.