Away-side charged hadron yield per π 0 trigger as a function of xE, which is equivalent to zT in the collinear limit cos(∆φ) = 1 & Away-side isolated direct photon trigger as a function of xE, which is equivalent to zT in the collinear limit cos(∆φ) = 1.

Away-side charged hadron yield per π 0 trigger as a function of xE, which is equivalent to zT in the collinear limit cos(∆φ) = 1 & Away-side isolated direct photon trigger as a function of xE, which is equivalent to zT in the collinear limit cos(∆φ) = 1.

Away-side charged hadron yield per π 0 trigger as a function of xE, which is equivalent to zT in the collinear limit cos(∆φ) = 1.

Away-side isolated direct photon trigger as a function of xE, which is equivalent to zT in the collinear limit cos(∆φ) = 1.

Away-side charged hadron yield per π 0 trigger as a function of xE, which is equivalent to zT in the collinear limit cos(∆φ) = 1 & Away-side isolated direct photon trigger as a function of xE, which is equivalent to zT in the collinear limit cos(∆φ) = 1.

Away-side charged hadron yield per π 0 trigger as a function of xE, which is equivalent to zT in the collinear limit cos(∆φ) = 1.

Away-side isolated direct photon trigger as a function of xE, which is equivalent to zT in the collinear limit cos(∆φ) = 1.

Breakdown of syst. uncertainties (fully bin-by-bin correlated) for differential cross section for $e^+e^-\rightarrow\pi^+\pi^-\gamma$, with $50^o<\theta_\gamma<130^o$

Bare cross section for $e^+e^-\rightarrow\pi^+\pi^-$

Statistical covariance matrix for bare cross section for $e^+e^-\rightarrow\pi^+\pi^-$

Inverse statistical covariance matrix for bare cross section for $e^+e^-\rightarrow\pi^+\pi$

Breakdown of syst. uncertainties (fully bin-by-bin correlated) for the cross section for $e^+e^-\rightarrow\pi^+\pi^-$

Pion form factor $|F_\pi|^2$

Statistical covariance matrix for pion form factor $|F_\pi|^2$

Inverse statistical covariance matrix $|F_\pi|^2$

Breakdown of syst. uncertainties (fully bin-by-bin correlated) for the pion form factor $|F_\pi|^2$

Vacuum polarisation correction $\delta_{VP}(s)=(\frac{\alpha_{em}(s)}{\alpha_{em}(0)})^2$, with $\sigma^{bare} = \sigma^{dressed}/\delta$ - calculated using the routine from F.Jegerlehner: https://web.archive.org/web/20170828074416/http://www-com.physik.hu-berlin.de/~fjeger/alphaQEDn.uu (2003, archived on August 28, 2017)

Final State Radiation correction $(1.+\eta_{FSR}(s))$, evaluated using the formula in Eur. Phys. J. C 24, 51-69 (2002)

The differential radiator function cross section $H(M_{\pi\pi}^2,s)\cdot \frac{\pi \alpha^2 \beta_\pi^3}{3 M_{\pi\pi}^2 s}$, inclusive in $\theta_\gamma$, in bins of 0.01 GeV$^2$ in $M_{\pi\pi}^2$. The value used for $s$ in the Monte Carlo production is $s = 999.85\; GeV^2$, corresponding to the mean value of DA$\Phi$NE energy for data collected in 2006. Obtained with the PHOKHARA5 MC generator, see Eur. Phys. J. C 27, 563 (2003) and https://looptreeduality.csic.es/phokhara/phokhara5.0.tar.gz

The K$\pi$ pair invariant mass distribution integrated over the $K^{*0}$ $p_T$ for minimum bias Cu+Cu collisions at $\sqrt{s_{NN}}$ =62.4 GeV after mixed-event background subtraction.

The Kπ pair invariant mass distribution for various pT bins (top left) pT = 0.4–0.6 GeV/c in Au+Au collisions at √sNN = 200 GeV after the mixed-event background subtraction.

The Kπ pair invariant mass distribution for various pT bins (top right) pT = 0.6–0.8 GeV/c in Au+Au collisions at √sNN = 62.4 GeV after the mixed-event background subtraction.

The Kπ pair invariant mass distribution for various pT bins (bottom left) pT = 0.8–1.0 GeV/c in Au+Au collisions at √sNN = 200 GeV after the mixed-event background subtraction.

The Kπ pair invariant mass distribution for various pT bins (bottom right) pT = 1.0–1.2 GeV/c in Au+Au collisions at √sNN = 62.4 GeV after the mixed-event background subtraction.

The signal-to-background ratio for $K^{*0}$ measurements as a function of $p_T$ for different collision centrality bins (0-10%, 10-40%, 40-60%, 60-80%) in Au+Au collisions at 200 GeV.

$K^{*0}$ mass as a function of $p_T$ for minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

$K^{*0}$ mass as a function of $p_T$ for minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

$K^{*0}$ mass as a function of $p_T$ for minimum bias Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

$K^{*0}$ mass as a function of $p_T$ for minimum bias Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

$K^{*0}$ width as a function of $p_T$ for minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

$K^{*0}$ width as a function of $p_T$ for minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

$K^{*0}$ width as a function of $p_T$ for minimum bias Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

$K^{*0}$ width as a function of $p_T$ for minimum bias Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

The $K^{*0}$ reconstruction efficiency multiplied by the detector acceptance as a function of $p_T$ in Au+Au (|$\eta$| < 0.8) collisions at 200 GeV for different collision centrality bins (0-20% ,20-40% , 40-60%)

The $K^{*0}$ reconstruction efficiency multiplied by the detector acceptance as a function of $p_T$ in Cu+Cu (|$\eta$| < 1.0) collisions at 200 GeV for different collision centrality bins (0-20% ,20-40% , 40-60%)

Mid-rapidity $K^{*0}$ $p_T$ spectra for various collision centrality bins (0-20%, 20-40%, 40-60%, 60-80%) in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

Mid-rapidity $K^{*0}$ $p_T$ spectra for various collision centrality bins (0-20%, 20-40%, 40-60%) in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

Mid-rapidity $K^{*0}$ $p_T$ spectra for various collision centrality bins (0-20%, 20-40%, 40-60%, 60-80%) in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

Mid-rapidity $K^{*0}$ $p_T$ spectra for various collision centrality bins (0-20%, 20-40%, 40-60%) in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

The mid-rapidity yields dN/dy of $K^{*0}$ as a function of the average number of participating nucleons, $⟨N_{part}⟩$, for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

The mid-rapidity yields dN/dy of $K^{*0}$ as a function of the average number of participating nucleons, $⟨N_{part}⟩$, for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

The mid-rapidity yields dN/dy of $K^{*0}$ as a function of the average number of participating nucleons, $⟨N_{part}⟩$, for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

The mid-rapidity yields dN/dy of $K^{*0}$ as a function of the average number of participating nucleons, $⟨N_{part}⟩$, for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

The mid-rapidity $K^{*0}$ $⟨p_T⟩$ as a function $⟨N_{part}⟩$ for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

The mid-rapidity $K^{*0}$ $⟨p_T⟩$ as a function $⟨N_{part}⟩$ for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV

The mid-rapidity $K^{*0}$ $⟨p_T⟩$ as a function $⟨N_{part}⟩$ for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV

The mid-rapidity $K^{*0}$ $⟨p_T⟩$ as a function $⟨N_{part}⟩$ for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV

The mid-rapidity $⟨p_T⟩$ of $\pi$, K, p and $K^{*0}$ as a function of $⟨N_{part}⟩$ for Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio for Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio for Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio for Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})N(K^-)$ in Au+Au collisions divided by $N(K^{*0})N(K^-)$ ratio in p+p collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$.

Mid-rapidity $N(K^{*0})N(K^-)$ in Cu+Cu collisions divided by $N(K^{*0})N(K^-)$ ratio in p+p collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})N(K^-)$ in d+Au collisions divided by $N(K^{*0})N(K^-)$ ratio in d+Au collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias Au+Au collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias Cu+Cu collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias p+p collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias Au+Au collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias Cu+Cu collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(K^{*0})/N(K^-)$ ratio in minimum bias p+p collisions as a function of $\sqrt{s_{NN}}.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio for Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio for Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio for Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $[N(\phi)/N(K^{*0})]$ in Au+Au collisions divided by $[N(\phi)/N(K^{*0})]$ ratio in p+p collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $[N(\phi)/N(K^{*0})]$ in Cu+Cu collisions divided by $[N(\phi)/N(K^{*0})]$ ratio in p+p collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $[N(\phi)/N(K^{*0})]$ in d+Au collisions divided by $[N(\phi)/N(K^{*0})]$ ratio in p+p collisions at $\sqrt{s_{NN}}$=200 GeV as a function of $⟨N_{part}⟩$

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias Au+Au collisions as a function of $\sqrt{s_{NN}}$.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias Cu+Cu collisions as a function of $\sqrt{s_{NN}}$.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias p+p collisions as a function of $\sqrt{s_{NN}}$.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias Au+Au collisions as a function of $\sqrt{s_{NN}}$.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias Cu+Cu collisions as a function of $\sqrt{s_{NN}}$.

Mid-rapidity $N(\phi)/N(K^{*0})$ ratio in minimum bias p+p collisions as a function of $\sqrt{s_{NN}}$.

The $K^{*0}$ $v_2$ (Run IV) as a function of $p_T$ in minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $K^{*0}$ $v_2$ (Run II) as a function of $p_T$ in minimum bias Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $K^{*0}$ $R_{CP}$ as a function of $p_T$ in Au+Au collisions at 62.4 and 200 GeV compared to the $R_{CP}$ of $K^0_S$ and $\Lambda$ at 200 GeV.

The $K^{*0}$ $R_{CP}$ as a function of $p_T$ in Au+Au collisions at 62.4 and 200 GeV compared to the $R_{CP}$ of $K^0_S$ and $\Lambda$ at 200 GeV.

The $K^{*0}$ $R_{CP}$ as a function of $p_T$ in Au+Au collisions at 62.4 and 200 GeV compared to the $R_{CP}$ of $K^0_S$ and $\Lambda$ at 200 GeV.

The $K^{*0}$ ~$R_{CP}$~ as a function of $p_T$ in Au+Au collisions at 62.4 and 200 GeV compared to the $R_{CP}$ of $K^0_S$ and $\Lambda$ at 200 GeV.

$v_2$ vs $N_{part}$ in two $p_T$ ranges and $R_{AA}$ vs $N_{part}$ in the same $p_T$ ranges.

$v_2$ vs $N_{part}$ in two $p_T$ ranges and $R_{AA}$ vs $N_{part}$ in the same $p_T$ ranges.

$v_2$ vs $N_{part}$ and $R_{AA}$ vs $N_{part}$ in two $p_T$ ranges compared with WHDG model.

$v_2$ vs $N_{part}$ in two $p_T$ ranges compared with various models.

$R_{AA}$ vs $N_{part}$ in two $p_T$ ranges compared with various models.

Measured differential yield of charged hadrons as a function oftransverse momentum for the pseudorapidity range -2.4 to +2.4 for centre-of-mass energy 7000 GeV. Errors are statistical and systematic added in quadrature.

Reconstructed differential charged hadron yields as a function of pseudorapidity averaged over the three methods methodsat a centre-of-mass energy of 7000 GeV.Errors shown are systematic (4 pct).Systematic errors are less than 0.5 pct of the values.

The COLLINS asymmetry for negatively charged hadrons as a function of X.

The COLLINS asymmetry for negatively charged hadrons as a function of Z.

The COLLINS asymmetry for negatively charged hadrons as a function of PT.

The SIVERS asymmetry for positively charged hadrons as a function of X.

The SIVERS asymmetry for positively charged hadrons as a function of Z.

The SIVERS asymmetry for positively charged hadrons as a function of PT.

The SIVERS asymmetry for negatively charged hadrons as a function of X.

The SIVERS asymmetry for negatively charged hadrons as a function of Z.

The SIVERS asymmetry for negatively charged hadrons as a function of PT.

The COLLINS asymmetry for positively charged hadrons as a function of Z for X > 0.05.

The COLLINS asymmetry for positively charged hadrons as a function of PT for X > 0.05.

The COLLINS asymmetry for negatively charged hadrons as a function of Z for X > 0.05.

The COLLINS asymmetry for negatively charged hadrons as a function of PT for X > 0.05.

The SIVERS asymmetry for positively charged hadrons as a function of Z for X > 0.032.

The SIVERS asymmetry for positively charged hadrons as a function of PT for X > 0.032.

The SIVERS asymmetry for negatively charged hadrons as a function of Z for X > 0.032.

The SIVERS asymmetry for negatively charged hadrons as a function of PT for X > 0.032.

The COLLINS asymmetry for positively charged hadrons as a function of W for X > 0.032.

The COLLINS asymmetry for negatively charged hadrons as a function of W for X > 0.032.

The SIVERS asymmetry for positively charged hadrons as a function of W for X > 0.032.

The SIVERS asymmetry for negatively charged hadrons as a function of W for X > 0.032.

The COLLINS asymmetry for positively charged hadrons as a function of W for X < 0.032.

The COLLINS asymmetry for negatively charged hadrons as a function of W for X < 0.032.

The SIVERS asymmetry for positively charged hadrons as a function of W for X < 0.032.

The SIVERS asymmetry for negatively charged hadrons as a function of W for X < 0.032.

The COLLINS asymmetry for positively charged hadrons as a function of X for W > 7.5 GeV.

The COLLINS asymmetry for negatively charged hadrons as a function of X for W > 7.5 GeV.

The SIVERS asymmetry for positively charged hadrons as a function of X for W > 7.5 GeV.

The SIVERS asymmetry for negatively charged hadrons as a function of X for W > 7.5 GeV.

The COLLINS asymmetry for positively charged hadrons as a function of X for W < 7.5 GeV.

The COLLINS asymmetry for negatively charged hadrons as a function of X for W < 7.5 GeV.

The SIVERS asymmetry for positively charged hadrons as a function of X for W < 7.5 GeV.

The SIVERS asymmetry for negatively charged hadrons as a function of X for W < 7.5 GeV.

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=1$ for $0-20$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{d}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=1$ for $0-20$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=2$ for $0-20$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{d}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=2$ for $0-20$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=3$ for $0-20$% central

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\mathrm{d}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=3$ for $0-20$% central

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

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\overline{\mathrm{d}}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=0$ for $0-20$% central

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

$\frac{1}{2\pi p_{\mathrm{T}}}\frac{\mathrm{d}^2N}{\mathrm{d}p_{\mathrm{T}}\mathrm{d}y}$ versus $p_{\mathrm{T}}$ for $\overline{\mathrm{d}}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=1$ for $0-20$% central

$B_{2}$ versus $p_{\mathrm{T}}/A$ for $\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=0$ for $0-20$% central

$B_{2}$ versus $p_{\mathrm{T}}/A$ for $\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=1$ for $0-20$% central

$B_{2}$ versus $p_{\mathrm{T}}/A$ for $\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=2$ for $0-20$% central

$B_{2}$ versus $p_{\mathrm{T}}/A$ for $\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=3$ for $0-20$% central

$B_{2}$ versus $p_{\mathrm{T}}/A$ for $\overline{\mathrm{p}}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=0$ for $0-20$% central

$B_{2}$ versus $p_{\mathrm{T}}/A$ for $\overline{\mathrm{p}}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=1$ for $0-20$% central

$\langle f\rangle$ versus $m_{\mathrm{T}}/A$ for $\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=0$ for $0-20$% central

$\langle f\rangle$ versus $m_{\mathrm{T}}/A$ for $\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=1$ for $0-20$% central

$\langle f\rangle$ versus $m_{\mathrm{T}}/A$ for $\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=2$ for $0-20$% central

$\langle f\rangle$ versus $m_{\mathrm{T}}/A$ for $\mathrm{p}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=3$ for $0-20$% central

$\langle f\rangle$ versus $m_{\mathrm{T}}/A$ for $\overline{\mathrm{p}}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=0$ for $0-20$% central

$\langle f\rangle$ versus $m_{\mathrm{T}}/A$ for $\overline{\mathrm{p}}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ near $y=1$ for $0-20$% central

$E\frac{dN^3}{dp^3}$ vs. $p_T$, 0% to 20% centrality $Au+Au$. 90% Limit on 18-20 and 20-22 GeV/c bins.

$E\frac{dN^3}{dp^3}$ vs. $p_T$, 20% to 40% centrality $Au+Au$. 90% Limit for 20-22 GeV/c bin.

$E\frac{dN^3}{dp^3}$ vs. $p_T$, 40% to 60% centrality $Au+Au$. 90% Limit for 18-20 GeV/c bin.

$E\frac{dN^3}{dp^3}$ vs. $p_T$, 20% to 60% centrality $Au+Au$. 90% Limit for 18-20 GeV/c bin.

$E\frac{dN^3}{dp^3}$ vs. $p_T$, 60% to 92% centrality $Au+Au$. 90% Limit for 18-20 GeV/c bin.

$E\frac{dN^3}{dp^3}$ vs. $p_T$, 0% to 92% centrality $Au+Au$.

$E\frac{dN^3}{dp^3}$ vs. $p_T$, $p+p$.

${\eta}$ $R_{AA}$, 0% to 5% centrality $Au+Au$.

${\eta}$ $R_{AA}$, 0% to 10% centrality $Au+Au$.

${\eta}$ $R_{AA}$, 10% to 20% centrality $Au+Au$.

${\eta}$ $R_{AA}$, 0% to 20% centrality $Au+Au$.

${\eta}$ $R_{AA}$, 20% to 40% centrality $Au+Au$.

${\eta}$ $R_{AA}$, 40% to 60% centrality $Au+Au$.

${\eta}$ $R_{AA}$, 20% to 60% centrality $Au+Au$.

${\eta}$ $R_{AA}$, 60% to 92% centrality $Au+Au$.

${\eta}$ $R_{AA}$, 0% to 92% centrality $Au+Au$.

${\pi^{0}}$ $R_{AA}$, 0% to 92% centrality $Au+Au$.

${\eta}$ ${R_{AA}}$, 0% to 92% centrality $Au+Au$.

${\eta}$ $R_{AA}$ vs centrality ($p_T>$20 GeV/c).

${\eta}$ $R_{AA}$ vs centrality ($p_T>$5 GeV/c).

Slope ${\eta}$ $R_{AA}$ ($5<p_T<$ 18 GeV/c) vs centrality.

$\chi^{2}/{NDF}$, Prob%, ${\eta}$ $R_{AA}$ ($5<p_T<$ 18 GeV/c) Prob, $R_{AA}$ ($p_T>$ 5 GeV/c) Prob vs centrality

The balance function in terms of $\Delta y$ for identified charged pion pairs from Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for nine centrality bins.

The balance function in terms of $\Delta y$ for identified charged kaon pairs from Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for nine centrality bins.

The balance function for p+p collisions at $\sqrt{s_{NN}}$ = 200 GeV. The top panel shows the balance function for all charged particles in terms of $\Delta \eta$. The bottom panel gives the balance function for charged pion pairs and charged kaon pairs in terms of $\Delta y$.

The balance function for p+p collisions at $\sqrt{s_{NN}}$ = 200 GeV. The top panel shows the balance function for all charged particles in terms of $\Delta \eta$. The bottom panel gives the balance function for charged pion pairs and charged kaon pairs in terms of $\Delta y$.

The balance function in terms of $\Delta \eta$ for all charged particles from d+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for three centrality bins.

The balance function in terms of $q_{inv}$ for charged pion pairs from Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV in nine centrality bins. Curves correspond to a thermal distribution (Eq. 3) plus $K^{0}_{S}$ decay.

The balance function in terms of $q_{inv}$ for charged-pion pairs in two centrality bins over the range 0 < $q_{inv}$ < 0.2 GeV/c.

The balance function in terms of $q_{inv}$ for charged kaon pairs from Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV in nine centrality bins. Curves correspond to a thermal (Eq. 3) distribution plus $\phi$ decay.

The balance function in terms of $q_{inv}$ for charged pion pairs [part a)] and charged kaon pairs [part b)] from p+p collisions at $\sqrt{s_{NN}}$ = 200 GeV integrated over all multiplicities. Solid curves correspond to a thermal distribution (Eq. 3) plus $K^{0}_{S}$ and $\rho^{0}$ decay for pions and φ decay for kaons. The dashed curve for pions represents a fit to a thermal distribution (Eq. 3) plus $K^{0}_{S}$ decay and $\rho^{0}$ decay, with the $\rho^{0}$ mass shifted down by 0.04 $GeV/c^{2}$.

The balance function in terms of $q_{long}$ for charged pion pairs from Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV in nine centrality bins.

The balance function in terms of $q_{out}$ for charged pion pairs from Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV in nine centrality bins.

(Color online) The balance function in terms of $q_{side}$ for charged pion pairs from Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for nine centrality bins.

(Color online) The balance function in terms of $\Delta\phi$ for all charged particles with 0.2< pt<2.0 GeV/c from Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV in nine centrality bins.The closed circles represent the real data minus the mixed events.

(Color online) The balance function in terms of $\Delta\phi$ for all charged particles with 1.0< pt<10.0 GeV/c from Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV in nine centrality bins.The closed circles represent the real data minus the mixed events.

(Color online) The balance function in terms of $\Delta y$ for charged pions from central collisions of Au+Au at $\sqrt{s_{NN}}$ = 200 GeV compared with predictions from the blast-wave model from Ref.[30] and filtered HIJING calculations taking into account acceptance and efficiency.

(Color online) The balance function in terms of $\Delta y$ for charged pions from central collisions of Au+Au at $\sqrt{s_{NN}}$ = 200 GeV compared with predictions from the blast-wave model from Ref.[30] and filtered HIJING calculations taking into account acceptance and efficiency.

(Color online) The balance function in terms of $q_{inv}$ for charged pions from central collisions of Au+Au at $\sqrt{s_{NN}}$ = 200 GeV compared with predictions from the blast-wave model from Ref.[30] and predictions from filtered HIJING calculations including acceptance and efficiency. For the blast-wave calculations, HBT is not included and the decays of the $K^0_S$ and $\rho^0$ are not shown.

(Color online) The balance function in terms of $q_{inv}$ for charged pions from central collisions of Au+Au at $\sqrt{s_{NN}}$ = 200 GeV compared with predictions from the blast-wave model from Ref.[30] and predictions from filtered HIJING calculations including acceptance and efficiency. For the blast-wave calculations, HBT is not included and the decays of the $K^0_S$ and $\rho^0$ are not shown.

(Color online) The balance function width $<\Delta\eta>$ for all charged particles from Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV compared with the widths of balance functions calculated using shuffled events. Also shown are the balance function widths for p+p and d+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.Filtered HIJING calculations are also shown for the widths of the balance function from p+p and Au+Au collisions. Filtered UrQMD calculations are shown for the widths of the balance function from Au+Au collisions.

(Color online) The balance function width $<\Delta\eta>$ for all charged particles from Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV compared with the widths of balance functions calculated using shuffled events. Also shown are the balance function widths for p+p and d+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.Filtered HIJING calculations are also shown for the widths of the balance function from p+p and Au+Au collisions. Filtered UrQMD calculations are shown for the widths of the balance function from Au+Au collisions.

(Color online) The balance function widths for identified charged pions and charged kaons from Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV and p+p collisions at $\sqrt{s}$ = 200 GeV. Filtered HIJING calculations are shown for the same systems. Filtered UrQMD calculations are shown for Au+Au. Also shown is the width of the balance function for pions predicted by the blast-wave model of Ref.[30].

(Color online) The balance function width $\sigma$ extracted from $B(q_{inv})$ for identified charged pions and kaons from Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV and p+p collisions at $\sqrt{s}$ = 200 GeV using a thermal fit (Eq. 3) where $\sigma$ is the width. Filtered HIJING and UrQMD calculations are shown for pions and kaons from Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Values are shown for $\sqrt{2mT_{kin}}$ from Au+Au collisions, where m is the mass of a pion or a kaon, and $T_{kin}$ is calculated from identified particle spectra [45]. The width predicted by the blast-wave model of Ref. [30] is also shown for pions.

(Color online) The balance function width $\sigma$ extracted from $B(q_{inv})$ for identified charged pions and kaons from Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV and p+p collisions at $\sqrt{s}$ = 200 GeV using a thermal fit (Eq. 3) where $\sigma$ is the width. Filtered HIJING and UrQMD calculations are shown for pions and kaons from Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Values are shown for $\sqrt{2mT_{kin}}$ from Au+Au collisions, where m is the mass of a pion or a kaon, and $T_{kin}$ is calculated from identified particle spectra [45]. The width predicted by the blast-wave model of Ref. [30] is also shown for pions.

(Color online) The widths for the balance functions for pions in terms of $q_{long}$, $q_{out}$, and $q_{side}$ compared with UrQMD calculations.

(Color online) The weighted average cosine of the relative azimuthal angle,$<\cos(\Delta\phi)$, extracted from $B(\Delta\phi)$ for all charged particles with $0.2< p_t <2.0$ from Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV and from all charged particleswith $1.0< p_t <10.0$ GeV/c, compared with predictions using filtered UrQMD calculations.

Parameters of the Tsallis fit with parameter $n$ constrained to a fixed linear dependence on mass (for mesons). Cross sections are in $\mu$b for $J/\psi$ and $\psi^{\prime}$ and in mb for all other particles.

Constant and linear fits to the Tsallis parameter $T$ of mesons with fixed parameter $n$. The last column (Prob.) gives the probability estimated by the $\chi^2$/$n.d.f.$ of the fit.

Parameters of the Tsallis fit with Eq. 8 in the paper with parameters $n$ and $T$ constrained to have a fixed linear dependence on mass (for mesons). Cross sections are in $\mu$b for $J/\psi$ and $\psi^{\prime}$ and in mb for all other particles.

Cross sections in mb and $<p_T>$ in GeV/$c$ of different particles in $p$+$p$ collisions at $\sqrt{s}$ = 200 GeV. PHENIX columns show the values obtained by fits to the experiment spectra with the Tsallis functional form as described in the text. One should state explicitly that these values do not supersede values given in the ”Published” column by the experiments, in their publications listed in the last column, or elsewhere. An additional 9.7% systematic uncertainty should be added to all $d\sigma$/$dy$ values listed in the column “PHENIX” to account for the trigger uncertainties. Values in the column “Published” are also given without these systematic uncertainties. The column “S.M.” is the prediction of the statistical model discussed in the text. The characteristic widths of the particle spectra are $\langle p_T \rangle$ for all species except for $J/\psi$ and $\psi^{\prime}$ for which the values given in the Table are $\langle p_T^2 \rangle$. For $\psi^{\prime}$ the integration is done in the $p_T$ regionbelow 5 GeV/$c$. All errors are the combined statistical and systematic uncertainties.

Cross sections in mb and $<p_T>$ in GeV/$c$ of different particles in $p$+$p$ collisions at $\sqrt{s}$ = 200 GeV. PHENIX columns show the values obtained by fits to the experiment spectra with the Tsallis functional form as described in the text. One should state explicitly that these values do not supersede values given in the ”Published” column by the experiments, in their publications listed in the last column, or elsewhere. An additional 9.7% systematic uncertainty should be added to all $d\sigma$/$dy$ values listed in the column “PHENIX” to account for the trigger uncertainties. Values in the column “Published” are also given without these systematic uncertainties. The column “S.M.” is the prediction of the statistical model discussed in the text. The characteristic widths of the particle spectra are $\langle p_T \rangle$ for all species except for $J/\psi$ and $\psi^{\prime}$ for which the values given in the Table are $\langle p_T^2 \rangle$. For $\psi^{\prime}$ the integration is done in the $p_T$ regionbelow 5 GeV/$c$. All errors are the combined statistical and systematic uncertainties.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.