Fourier fit coefficients for CNT-MPCN ($d$-going) correlations, as a function of collision system and $\pi^0$ $p_T$: (a) the negative of the dipole coefficient, $-c_1$; (b) the quadrupole coefficient $c_2$; (c) the ratio ${-c_2}/{c_1}$.

Fourier fit coefficients for CNT-MPCN ($d$-going) correlations, as a function of collision system and $\pi^0$ $p_T$: Fractional systematic uncertainty on the quadrupole coefficient $c_2$ for $d$+Au.

Fourier fit coefficients for CNT-MPCN ($d$-going) correlations, as a function of collision system and $\pi^0$ $p_T$: Fractional systematic uncertainty on the quadrupole coefficient $c_2$ for $d$+Au and $p$+$p$.

System centrality dependence of correlation coefficients $-c_1$, $c_2$, and the ratio ${-c_2}/{c_1}$ as a function of $N_{part}$.

System centrality dependence of correlation coefficients $-c_1$, $c_2$, and the ratio ${-c_2}/{c_1}$ as a function of $N_{part}$.

$J/\psi$ $R_{CP}$ for 0-20%, 20-40%, and 40-60% (central) relative to 60-86% (peripheral) Au+Au collisions at 62.4 GeV. Type A errors are point-to-point uncorrelated, Type B errors are point-to-point correlated, and Type C are global systematic from the uncertainty in the peripheral $<N_{coll}>$ value.

$J/\psi$ $R_{CP}$ for 0-40% (central) relative to 40-86% (peripheral) Au+Au collisions at 39 GeV. Type A errors are point-to-point uncorrelated, Type B errors are point-to-point correlated, and Type C are global systematic from the uncertainty in the peripheral $<N_{coll}>$ value.

$J/\psi$ $R_{AA}$ at $\sqrt{s_{NN}}$ = 39 and 62.4 GeV. Type A errors are point-to-point uncorrelated, Type B errors are point-to-point correlated, and Type C are global systematic from the uncertainty in the peripheral $<N_{coll}>$ value.

$J/\psi$ $R_{AA}$ at $\sqrt{s_{NN}}$ = 39 and 62.4 GeV.

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.

$R_{AA}$ for $\phi$ integrated at $p_T$ > 2.2 GeV/$c$ in Cu+Cu and Au+Au collisions vs. $N_{part}$. The global uncertainty related to the $p$+$p$ reference normalization is shown as a box on the right in the figure from the paper.

$R_{AA}$ for $\phi$ integrated at $p_T$ > 2.2 GeV/$c$ in Cu+Cu and Au+Au collisions vs. $N_{part}$. The global uncertainty related to the $p$+$p$ reference normalization is shown as a box on the right in the figure from the paper.

$R_{dA}$ vs. $p_T$ for $\phi$, $\pi^0$, and ($p$+$\bar{p}$) for 0-20% centrality $d$+Au collisions, and $R_{dA}$ vs. $p_T$ for $\phi$, $\pi^0$, and ($p$+$\bar{p}$) for 60-88% peripheral $d$+Au collisions. The global uncertainty of ~ 10% related to the $p$+$p$ reference normalization is not shown.

Additional information containing number of events which were used to reconstruct the numvers matching to Figure 1 and 2.

Additional information containing number of events which were used to reconstruct the numvers matching to Figure 1 and 2.

Additional information containing number of events which were used to reconstruct the numvers matching to Figure 1 and 2.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 22.5 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 22.5 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV $p$+$p$ collisions.

The mean from the NBD fit as a function of $N_{part}$ for 200 GeV Au+Au collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The mean from the NBD fit as a function of $N_{part}$ for 62.4 GeV Au+Au collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The mean from the NBD fit as a function of $N_{part}$ for 200 GeV Cu+Cu collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The mean from the NBD fit as a function of $N_{part}$ for 62.4 GeV Cu+Cu collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The mean from the NBD fit as a function of $N_{part}$ for 22.5 GeV Cu+Cu collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of centrality for 200 GeV Au+Au collisions in the range 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 200 GeV Au+Au collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 62.4 GeV Au+Au collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 200 GeV Cu+Cu collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 62.4 GeV Cu+Cu collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 22.5 GeV Cu+Cu collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Scaled variance for 200 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 200 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 200 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 200 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 200 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The scaled variance as a funciton of $N_{part}$ for 200 GeV Au+Au collisions in the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The fraction of the direct photon component as a function of $p_T$.

The fraction of the direct photon component as a function of $p_T$.

Invariant yield (Au+Au) of direct photons as a function of $p_T$.

Invariant yield (Au+Au) of direct photons as a function of $p_T$.

Invariant yield (Au+Au) of direct photons as a function of $p_T$.

Invariant cross section ($p$+$p$) of direct photons as a function of $p_T$.

Measured widths of the near- and away-angle $\pi^0$ - $h^{\pm}$ correlation peaks for various trigger momenta.

Measured widths of the near- and away-angle $\pi^0$ - $h^{\pm}$ correlation peaks for various trigger momenta.

Measured widths of the near- and away-angle $\pi^0$ - $h^{\pm}$ correlation peaks for various trigger momenta.

The $\sigma_N$ and $\sigma_A$ values in GeV/$c$.

The $\sqrt{<p^2_{out}>}$ values in GeV/$c$.

The $\sqrt{<p^2_{out}>}$ values in GeV/$c$.

The $\sqrt{<p^2_{out}>}$ values in GeV/$c$.

The near and away side conditional yield per number of triggers for 1.4 < $p_{Ta}$ < 5.0 GeV/$c$.

The $\sqrt{<j^2_T>}$ values in GeV/$c$.

The $\sqrt{<j^2_T>}$ values in GeV/$c$.

The $\sqrt{<j^2_T>}$ values in GeV/$c$.

The $\sqrt{<j^2_T>}$ values in GeV/$c$.

Measured $x_E$ distributions associated with various transverse momenta of the trigger $\pi^0$.

Measured $x_E$ distributions associated with various transverse momenta of the trigger $\pi^0$.

Measured $x_E$ distributions associated with various transverse momenta of the trigger $\pi^0$.

Measured $x_E$ distributions associated with various transverse momenta of the trigger $\pi^0$.

Measured $x_E$ distributions associated with various transverse momenta of the trigger $\pi^0$.

Extracted values of $D(x)$ parameters according from the fit to the LEP data and power $n$ of the unmeasured final state parton spectra $\Sigma_q(\bar{p_T})$ extracted from the fit to the single inclusive $\pi^0$ invariant cross section for corresponding fragmentation and fixed values of $\sqrt{<k^2_T>}$ = 2.5 GeV/$c$.

The $\bar{x^-1_h}$ $<z_t>$ $\sqrt{<k^2_T>}$ $\sqrt{<k^2_T>}$ values as a function of $p_{Ta}$ for two different trigger $\pi^0$ transverse momentum bins.

The $\bar{x^-1_h}$ $<z_t>$ $\sqrt{<k^2_T>}$ $\sqrt{<k^2_T>}$ values as a function of $p_{Ta}$ for two different trigger $\pi^0$ transverse momentum bins.

Values of $\bar{x^-1_h}$ $<z_t>$ $\sqrt{<k^2_T>}$ and $\sqrt{<k^2_T>}$ for various trigger particle $p_{Tt}$ and associated momenta in 1.4 < $p_{Ta}$ < 5.0 GeV/$c$ region.

The $<z_t>$ and $\bar{x_h}$ values with $p_{Tt}$

The $<z_t>$ and $\bar{x_h}$ values with $p_{Ta}$ for two trigger $\pi^0$ momenta bins.

The $<z_t>$ and $\bar{x_h}$ values with $p_{Ta}$ for two trigger $\pi^0$ momenta bins.