Average transverse momentum per event for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV for the 5% most central collisions.

$<\Delta p_{t,i}\Delta p_{t,j}>$ as a function of centrality and incident energy for Au+Au collisions compared with HIJING results.

(d$N/\textrm{d}\eta)<\Delta p_{t,i}\Delta p_{t,j}>$ as a function of centrality and incident energy for Au+Au collisions compared with HIJING results.

$(<\Delta p_{t,i}\Delta p_{t,j}>)^{1/2}/<<p_{t}>>$ as a function of centrality and incident energy for Au+Au collisions compared with HIJING results.

$(<\Delta p_{t,i}\Delta p_{t,j}>)^{1/2}/<<p_{t}>>$ as a function of incident energy for 0-5% most central Au+Au collisions compared with CERES results.

$v_2(p_T)$ of $\Xi^{-} + \Xi^{+}$ and $\Omega^{-} + \Omega^{+}$ from 200 GeV Au+Au minimum bias collisions.

Number of quark ($n_q$) scaled $v_2$ as a function of scaled $p_T$ for $\Xi^{-} + \Xi^{+}$ and $\Omega^{-} + \Omega^{+}$ from 200 GeV Au+Au minimum bias collisions.

$\frac{\mathrm{d}N}{\mathrm{d}y}\frac{2}{N_{\mathrm{part}}}$ versus $N_{\mathrm{part}}$ for $\mathrm{\pi}^{+}$,$\mathrm{\pi}^{-}$,$\mathrm{K}^{+}$,$\mathrm{K}^{-}$,$\mathrm{p}$,$\overline{\mathrm{p}}$,$\mathrm{p}-\overline{\mathrm{p}}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ for $0-1$% central

$\mathrm{K}^{+}/\mathrm{\pi}^{+}$ versus $p_{\mathrm{T}}$ for $\mathrm{K}^{+}$,$\mathrm{\pi}^{+}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=0$ for $0-10$% central, $40-60$% central

$\mathrm{K}^{-}/\mathrm{\pi}^{-}$ versus $p_{\mathrm{T}}$ for $\mathrm{K}^{-}$,$\mathrm{\pi}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=0$ for $0-10$% central, $40-60$% central

$\mathrm{K}^{+}/\mathrm{\pi}^{+}$ versus $p_{\mathrm{T}}$ for $\mathrm{K}^{+}$,$\mathrm{\pi}^{+}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=1$ for $0-10$% central, $40-60$% central

$\mathrm{K}^{-}/\mathrm{\pi}^{-}$ versus $p_{\mathrm{T}}$ for $\mathrm{K}^{-}$,$\mathrm{\pi}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=1$ for $0-10$% central, $40-60$% central

$\mathrm{p}/\mathrm{\pi}^{+}$ versus $p_{\mathrm{T}}$ for $\mathrm{p}$,$\mathrm{\pi}^{+}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=0$ for $0-10$% central, $40-60$% central

$\overline{\mathrm{p}}/\mathrm{\pi}^{-}$ versus $p_{\mathrm{T}}$ for $\overline{\mathrm{p}}$,$\mathrm{\pi}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=0$ for $0-10$% central, $40-60$% central

$\mathrm{p}/\mathrm{\pi}^{+}$ versus $p_{\mathrm{T}}$ for $\mathrm{p}$,$\mathrm{\pi}^{+}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=1$ for $0-10$% central, $40-60$% central

$\overline{\mathrm{p}}/\mathrm{\pi}^{-}$ versus $p_{\mathrm{T}}$ for $\overline{\mathrm{p}}$,$\mathrm{\pi}^{-}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=1$ for $0-10$% central, $40-60$% central

$\mathrm{K}^{+}/\mathrm{\pi}^{+}$,$\mathrm{K}^{-}/\mathrm{\pi}^{-}$,$\mathrm{p}/\mathrm{\pi}^{+}$,$\overline{\mathrm{p}}/\mathrm{\pi}^{-}$ versus $N_{\mathrm{part}}$ for $\mathrm{K}^{+}$,$\mathrm{K}^{-}$,$\mathrm{p}$,$\overline{\mathrm{p}}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=0$

Direct $\gamma$ invariant yields as a function of transverse momentum for 9 centrality selections and minimum bias Au+AU collisions at $\sqrt{s_{NN}}$ = 200 GeV. Data with no errors represents 90% confidence level upper limit. The bin range is not an uncertainty in the x-axis because the actual uncertainty by having the finite bin width is corrected for by the bin-shift correction. These bins were constructed using the corrected finite values as centers.

Ratio of Au+Au yield to $p$+$p$ yield normalized by the number of binary nucleon collisions as a function of centrality given by $N_{part}$ for direct $\gamma$ and $\pi^0$ yields integrated above 6 GeV/$c$. Data with no errors represents 90% confidence level upper limit.

Dphi correlation functions for 2 < pT < 4 GEV/c and 4 < p_T^trig < 6 GEV/c.

Deta correlation functions for 0.15 < pT < 4 GEV/c and 4 < p_T^trig < 6 GEV/c.

Deta correlation functions for 0.15<pT<4 GEV/c and 4 < p_T^trig < 6 GEV/c.

Deta awayside correlation functions for 0.15 < pT < 4 GEV/c and 4 < p_T^trig < 6 GEV/c.

Deta correlation functions for 2<pT<4 GEV/c and 4 < p_T^trig < 6 GEV/c.

Dphi correlation functions for 2<pT<4 GEV/c and 4 < p_T^trig < 6 GEV/c.

Associated multiplicity vs centrality for 4 < p_T^trig < 6 GEV/c.

Associated multiplicity vs centrality for 6 < p_T^trig < 10 GEV/c.

Associated total p_T vs centrality for 4 < p_T^trig < 6 GEV/c.

Associated total p_T vs centrality for 6 < p_T^trig < 10 GEV/c.

Near-side associated charged hadron p_T spectra for 4 < p_T^trig < 6 GEV/c.

Near-side associated charged hadron p_T spectra for 4 < p_T^trig < 6 GEV/c.

Near-side associated charged hadron p_T spectra for 4 < p_T^trig < 6 GEV/c.

Away-side associated charged hadron p_T spectra for 4 < p_T^trig < 6 GEV/c.

Away-side associated charged hadron p_T spectra for 4 < p_T^trig < 6 GEV/c.

Away-side associated charged hadron p_T spectra for 4 < p_T^trig < 6 GEV/c.

Near-side AA/pp ratio for 4 < p_T^trig < 6 GEV/c.

Near-side AA/pp ratio for 4 < p_T^trig < 6 GEV/c.

Away-side AA/pp ratio for 4 < p_T^trig < 6 GEV/c.

Away-side AA/pp ratio for 4 < p_T^trig < 6 GEV/c.

Away-side associated <p_T> vs centrality for 4 < p_T^trig < 6 GEV/c.

Away-side associated <p_T> vs centrality for 6 < p_T^trig < 10 GEV/c.

The $K^{*0}$ mass (upper panel) and width (lower panel) as a function of $p_T$ for minimum bias $p + p$ interactions and for central Au+Au collisions. The solid straight lines are the standard $K^{*0}$ mass (896.1 MeV/$c^2$ ) and width (50.7 MeV/$c^2$), respectively. The dashed and dotted curves are the MC results in minimum bias $p+p$ and for central Au+Au collisions, respectively, after considering detector effects and kinematic cuts. The grey shadows (caps) indicate the systematic uncertainties for the measurement in minimum bias $p + p$ interactions (central Au+Au collisions).

The $K^{*0}$ mass (upper panel) and width (lower panel) as a function of $p_T$ for minimum bias $p + p$ interactions and for central Au+Au collisions. The solid straight lines are the standard $K^{*0}$ mass (896.1 MeV/$c^2$ ) and width (50.7 MeV/$c^2$), respectively. The dashed and dotted curves are the MC results in minimum bias $p+p$ and for central Au+Au collisions, respectively, after considering detector effects and kinematic cuts. The grey shadows (caps) indicate the systematic uncertainties for the measurement in minimum bias $p + p$ interactions (central Au+Au collisions).

The $K^{*0}$ mass (upper panel) and width (lower panel) as a function of $p_T$ for minimum bias $p + p$ interactions and for central Au+Au collisions. The solid straight lines are the standard $K^{*0}$ mass (896.1 MeV/$c^2$ ) and width (50.7 MeV/$c^2$), respectively. The dashed and dotted curves are the MC results in minimum bias $p+p$ and for central Au+Au collisions, respectively, after considering detector effects and kinematic cuts. The grey shadows (caps) indicate the systematic uncertainties for the measurement in minimum bias $p + p$ interactions (central Au+Au collisions).

The $K^{*0}$ mass (upper panel) and width (lower panel) as a function of $p_T$ for minimum bias $p + p$ interactions and for central Au+Au collisions. The solid straight lines are the standard $K^{*0}$ mass (896.1 MeV/$c^2$ ) and width (50.7 MeV/$c^2$), respectively. The dashed and dotted curves are the MC results in minimum bias $p+p$ and for central Au+Au collisions, respectively, after considering detector effects and kinematic cuts. The grey shadows (caps) indicate the systematic uncertainties for the measurement in minimum bias $p + p$ interactions (central Au+Au collisions).

The $K_S^0\pi^\pm$ invariant mass distribution integrated over the $K^{*\pm}$ $p_T$ for minimum bias $p + p$ collisions (upper panel) and for the $50-80\%$ of the inelastic hadronic Au+Au crosssection (lower panel) after the mixed-event background subtraction. The solid curves are fits to Eq. 10 and the dashed lines are the linear function representing the residual background.

The $K^{*0}$ and $K^{*\pm}$ reconstruction efficiency multiplied by the detector acceptance as a function of $p_T$ in minimum bias $p + p$ and different centralities in Au+Au collisions.

The $K^{*0}$ and $K^{*\pm}$ reconstruction efficiency multiplied by the detector acceptance as a function of $p_T$ in minimum bias $p + p$ and different centralities in Au+Au collisions.

The $K^{*0}$ and $K^{*\pm}$ reconstruction efficiency multiplied by the detector acceptance as a function of $p_T$ in minimum bias $p + p$ and different centralities in Au+Au collisions.

The $K^{*0}$ and $K^{*\pm}$ reconstruction efficiency multiplied by the detector acceptance as a function of $p_T$ in minimum bias $p + p$ and different centralities in Au+Au collisions.

The $(K^* + \bar{K^*})/2$ invariant yields as a function of $m_T − m_0$ for |y| < 0.5 from minimum bias $p + p$ and different centralities in Au+Au collisions. The top $10\% $central data have been multiplied by two for clarity. The lines are fits to Eq. 12. The errors shown are the quadratic sum of statistical and systematic (in the level of $10\% $) uncertainties.

The $(K^* + \bar{K^*})/2$ invariant yields as a function of $m_T − m_0$ for |y| < 0.5 from minimum bias $p + p$ and different centralities in Au+Au collisions. The top $10\% $central data have been multiplied by two for clarity. The lines are fits to Eq. 12. The errors shown are the quadratic sum of statistical and systematic (in the level of $10\% $) uncertainties.

The $(K^* + \bar{K^*})/2$ invariant yields as a function of $m_T − m_0$ for |y| < 0.5 from minimum bias $p + p$ and different centralities in Au+Au collisions. The top $10\% $central data have been multiplied by two for clarity. The lines are fits to Eq. 12. The errors shown are the quadratic sum of statistical and systematic (in the level of $10\% $) uncertainties.

(Color online) The invariant yields for both $(K^{*0} + \bar{K^{*0}})/2$ and $(K^{*+} + K^{*-})/2$ as a function of $p_T$ for |y| < 0.5 in minimum bias $p + p$ interactions. The dotted curve is the fit to the power-law function from Equation 13 for $p_T$ > 0.5 GeV/$c$ and extended to lower values of $p_T$ . The dashed curve is the $K^{*0}$ spectrum fit to the exponential function from Equation 12 and extended to higher values of $p_T$ . The dashed-dotted curve is the fit to the Levy function from Equation 14 for $p_T$ < 4 GeV/$c$. Errors are statistical only.

(Color online) The invariant yields for both $(K^{*0} + \bar{K^{*0}})/2$ and $(K^{*+} + K^{*-})/2$ as a function of $p_T$ for |y| < 0.5 in minimum bias $p + p$ interactions. The dotted curve is the fit to the power-law function from Equation 13 for $p_T$ > 0.5 GeV/$c$ and extended to lower values of $p_T$ . The dashed curve is the $K^{*0}$ spectrum fit to the exponential function from Equation 12 and extended to higher values of $p_T$ . The dashed-dotted curve is the fit to the Levy function from Equation 14 for $p_T$ < 4 GeV/$c$. Errors are statistical only.

(Color online) The $K^*$ $<p_T>$ as a function of $dN_{ch}/d\eta$ compared to that of $\pi^−$, $K^−$, and $p$ for minimum bias $p+p$ (solid symbols) and Au+Au (open symbols) collisions. The errors shown are the quadratic sum of the statistical and systematic uncertainties.

The $K^*/K$ (upper panel) and $\phi/K^*$ (lower panel) yield ratios as a function of the c.m. system energies. The yield ratios for central Au+Au collisions at $\sqrt{s_{NN}}$=130 [35] and 200 GeV and minimum bias $p+p$ interactions at $\sqrt{s_{NN}}$=200 GeV are compared to measurements from e+e- at sqrt(s) of 10.45 GeV [48], 29 GeV [49] and 91 GeV [50,51], pbar-p at sqrt(s) of 5.6 GeV [52] and pp at sqrt(s) of 27.5 GeV [53], 52.5 GeV [54] and 63 GeV [55]. The errors at $\sqrt{s_{NN}}$=130 and 200 GeV correspond to the quadratic sum of the statistical and systematic errors.

The $K^*/K$ (upper panel) and $\phi/K^*$ (lower panel) yield ratios as a function of the c.m. system energies. The yield ratios for central Au+Au collisions at $\sqrt{s_{NN}}$=130 [35] and 200 GeV and minimum bias $p+p$ interactions at $\sqrt{s_{NN}}$=200 GeV are compared to measurements from e+e- at sqrt(s) of 10.45 GeV [48], 29 GeV [49] and 91 GeV [50,51], pbar-p at sqrt(s) of 5.6 GeV [52] and pp at sqrt(s) of 27.5 GeV [53], 52.5 GeV [54] and 63 GeV [55]. The errors at $\sqrt{s_{NN}}$=130 and 200 GeV correspond to the quadratic sum of the statistical and systematic errors.

The $K^{*0}$ v2 (filled stars) as a function of $p_T$ for minimum bias Au+Au collisions compared to the $K^{0}_S$ (open triangles), $\Lambda$ (open circles), and charged hadron (open diamonds) $v_2$. The errors shown are statistical only.

The $K^*$ $R_{AA}$ (filled triangles) and $R_{CP}$ (filled circles) as a function of $p_T$ compared to the $K^{0}_S$ (open circles) and $\Lambda$ (open triangles) $R_{CP}$. The errors shown are statistical only. The dashed line represents the number of binary collisions scaling. The widths of the grey bands represent the systematic uncertainties of $R_{AA}$ (left) and $R_{CP}$ (right) due to the model calculations of $N_{bin}$.

Two-particle CI joint autocorrelations $\widehat{N}(\widehat{r}-1)$ on $(\eta_{\Delta}, \phi_{\Delta})$ for peripheral collisions.

HBT parameters for the three possible fitting procedures to the correlation functions described in this paper depending on how Coulomb interaction is taken into account from the 0-5% most central events.

HBT parameters for the three possible fitting procedures to the correlation functions described in this paper depending on how Coulomb interaction is taken into account from the 0-5% most central events.

HBT parameters for the three possible fitting procedures to the correlation functions described in this paper depending on how Coulomb interaction is taken into account from the 0-5% most central events.

HBT parameters for the three possible fitting procedures to the correlation functions described in this paper depending on how Coulomb interaction is taken into account from the 0-5% most central events.

1D correlation function for pi+pi- compared to standard, Bowler-Sinyukov, and a thoretical calculation that includes Coulomb and strong interactions.

Momentum resolution for pions at midrapidity expressed by the widths dp_T/p_T, dphi and dtheta as a function of p.

Squared HBT radii relative to the reaction plane angle, without and with the reaction plane resolution applied, for two different centrality-kT ranges. 20-30%, 0.15 < kT < 0.25, RP uncorrected.

Squared HBT radii relative to the reaction plane angle, without and with the reaction plane resolution applied, for two different centrality-kT ranges. 20-30%, 0.15 < kT < 0.25, RP res. corrected.

Squared HBT radii relative to the reaction plane angle, without and with the reaction plane resolution applied, for two different centrality-kT ranges. 10-20%, 0.35 < kT < 0.45, RP uncorrected.

Squared HBT radii relative to the reaction plane angle, without and with the reaction plane resolution applied, for two different centrality-kT ranges. 10-20%, 0.35 < kT < 0.45, RP res. uncorrected.

Squared HBT radii relative to the reaction plane angle, for the case where the lambda parameter is averaged and fixed in the fit, and the case where lambda is a free fit parameter for each Phi. 20-30%, 0.15 < kT < 0.25, fixed lambda.

Squared HBT radii relative to the reaction plane angle, for the case where the lambda parameter is averaged and fixed in the fit, and the case where lambda is a free fit parameter for each Phi. 20-30%, 0.15 < kT < 0.25, lambda varied.

Squared HBT radii relative to the reaction plane angle, for the case where the lambda parameter is averaged and fixed in the fit, and the case where lambda is a free fit parameter for each Phi. 10-20%, 0.35 < kT < 0.45, lambda fixed.

Squared HBT radii relative to the reaction plane angle, for the case where the lambda parameter is averaged and fixed in the fit, and the case where lambda is a free fit parameter for each Phi. 10-20%, 0.35 < kT < 0.45, lambda varied.

Fourier coefficients as a function of the maximum fraction of merged hits for the 5% most central events and kT between 150 and 250 MeV/c.

Projections of the 3 dimensional correlation function and fits to Eq. (10) and with the Edgeworth expansion to Eq. (17) to 6th order.

HBT parameters for the 0-5% most central events for fits to Eq. (10) and to Eq. (17) with no expansion.

HBT parameters for the 0-5% most central events for fits to Eq. (10) and to Eq. (17) upto 4th order expansion.

HBT parameters for the 0-5% most central events for fits to Eq. (10) and to Eq. (17) upto 6th order expansion.

HBT parameters for the 0-5% most central events for pi+pi+ and pi-pi- correlation functions.

HBT parameters from STAR and PHENIX at the same beam energy for the 0-30% most central events.

Energy dependence of the pi- HBT parameters for central Au+Au, Pb+Pb, and Pb+Au collisions at midrapidity and k_T ~ 0.2 GeV/c.

HBT parameters vs. m_T for 6 different centralities.

Fourier coefficients of azimuthal oscillations of HBT radii vs. numer of participating nucleons, for pi+ and pi- pairs separately (0.25 < kT < 0.35 GeV/c).

Comparison between the HBT radii obtained from the azimuthally integrated (traditional) HBT analysis and the 0th-order Fourier coefficients from the azimuthally-sensitive HBT analysis.

HBT radius parametres for 6 different centralities.

Extracted parameters R' and alpha, from the power-law fits to the HBT radius parameters.

Extracted freeze-out source radius extracted from a blast wave fit; source radius Rgeom from fits to Rside; and 2*Rside for the lowest k_T bin as a function of the number of participants.

HBT parameter Rside.

R - Rinitial (top panel) and R/Rinitial (bottom panel) for the azimuthally integrated analysis for the azimuthally-sensitive case vs. number of participants.

R - Rinitial (top panel) and R/Rinitial (bottom panel) for in the x (in-plane) and y (out-of-plane) directions for the azimuthally-sensitive case vs. number of participants.

R - Rinitial for the azimuthally integrated analysis for the azimuthally-sensitive case vs. (dN/dy)/Rinitial.

R - Rinitial in the x (in-plane) direction for the azimuthally-sensitive case vs. (dN/dy)/Rinitial.

R - Rinitial in the y (out-of-plane) direction for the azimuthally-sensitive case vs. (dN/dy)/Rinitial.

Longitudinal HBT radius Rl.

Evolution time tau vs. number of participants as extracted from a fit to Rl, lines in Fig. 26 (triangles), and from a blast wave fit to HBT parameters and spectra (circles).

R - Rinital for the azimuthally integrated analysis and for the in-plane and out-of-plane directions vs. beta_T,max*tau.

R - Rinital for the azimuthally integrated analysis and for the in-plane and out-of-plane directions vs. beta_T,max*tau.

R - Rinital for the azimuthally integrated analysis and for the in-plane and out-of-plane directions vs. beta_T,max*tau.

Emission duration time dTau vs. number of participants as extracted using a blast wave fit to HBT parameters and spectra.

Final source eccentricity as calculated from the Fourier coefficients (2R_{s,2}^2/R_{s,0}^2) and from the final in-plane and out-of-plane radii vs. initial eccentricity.