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.

$F_{p_T}$ (in percent) of non-random fluctuations as a function of the $p_T$ range over which $M_{p_T}$ is calculated, 0.2 GeV/$c$ < $p_T$ < $p^{max}_T$, for the 20-25% centrality class ($N_{part}$ = 181.6).

$\Delta\phi$ distributions for meson triggers with 2.5 < $p_T$ < 4.0 GeV/$c$ and associated charged hadrons with 1.7 < $p_T$ < 2.5 GeV/$c$ for 60-70% centrality in Au+Au collisions.

$\Delta\phi$ distributions for baryon triggers with 2.5 < $p_T$ < 4.0 GeV/$c$ and associated charged hadrons with 1.7 < $p_T$ < 2.5 GeV/$c$ for 60-70% centrality in Au+Au collisions.

Combinatorial Background Levels from Figure 1 Top Panels.

Combinatorial Background Levels from Figure 1 Bottom Panels.

Trigger Particle $v_2$ values.

Partner Particle $v_2$ values.

Yield per trigger for associated charged hadrons between 1.7 < $p_T$ < 2.5 GeV/$c$ for the near- (top) and away- (bottom) side jets.

Yield per trigger for associated charged hadrons between 1.7 < $p_T$ < 2.5 GeV/$c$ for the near- (top) and away- (bottom) side jets.

Uncertainties in Figure 2 as a function of $N_{part}$ and trigger type for AuAu and $d$+Au.

Uncertainties in Figure 2 as a function of $N_{part}$ and trigger type for AuAu and $d$+Au.

$p_T$ spectra of the near side associated charged hadrons corrected to the full jet yield for meson triggers at 2.5 < $p_T$ < 4.0 GeV/$c$ and $|\eta|$ < 0.35 for six centralities in Au+Au and $d$+Au collisions.

$p_T$ spectra of the near side associated charged hadrons corrected to the full jet yield for meson triggers at 2.5 < $p_T$ < 4.0 GeV/$c$ and $|\eta|$ < 0.35 for six centralities in Au+Au and $d$+Au collisions.

$p_T$ spectra of the near side associated charged hadrons corrected to the full jet yield for baryon triggers at 2.5 < $p_T$ < 4.0 GeV/$c$ and $|\eta|$ < 0.35 for six centralities in Au+Au, $d$+Au and $p$+$p$ collisions.

Inverse slopes from the fits in Figure 3.

Inverse slopes from the fits in Figure 3.

Inverse slopes from the fits in Figure 3.

Transverse momentum in GeV/$c$ for $K^{\pm}$.

Transverse momentum in GeV/$c$ for $K^{\pm}$.

Transverse momentum in GeV/$c$ for $p$ and $\bar{p}$.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The $p_{out}$ distributions are shown for dihadron and direct photon-hadron correlations, binned in $p^{trig}_T$ in $x_E$.

The $p_{out}$ distributions are shown for dihadron and direct photon-hadron correlations, binned in $p^{trig}_T$ in $x_E$.

The $p_{out}$ distributions are shown for dihadron and direct photon-hadron correlations, binned in $p^{trig}_T$ in $x_E$.

The $p_{out}$ distributions are shown for dihadron and direct photon-hadron correlations, binned in $p^{trig}_T$ in $x_E$.

The $p_{out}$ distributions are shown for dihadron and direct photon-hadron correlations, binned in $p^{trig}_T$ in $x_E$.

The $p_{out}$ distributions are shown for dihadron and direct photon-hadron correlations, binned in $p^{trig}_T$ in $x_E$.

The $p_{out}$ distributions are shown for dihadron and direct photon-hadron correlations, binned in $p^{trig}_T$ in $x_E$.

The $p_{out}$ distributions are shown for dihadron and direct photon-hadron correlations, binned in $p^{trig}_T$ in $x_E$.

The $p_{out}$ distributions are shown for dihadron and direct photon-hadron correlations, binned in $p^{trig}_T$ in $x_E$.

The nonperturbative away-side jet widths as a function of $p^{trig}_T$ at $\sqrt{s}$ = 200 GeV for both direct photon-hadron and dihadron correlations.

The nonperturbative away-side jet widths as a function of $p^{trig}_T$ at $\sqrt{s}$ = 200 GeV for both direct photon-hadron and dihadron correlations.

The $p_{out}$ distributions are shown in several bins of $x_E$, integrated over a range of $p_T^{trig}$.

The $p_{out}$ distributions are shown in several bins of $x_E$, integrated over a range of $p_T^{trig}$.

The $p_{out}$ distributions are shown in several bins of $x_E$, integrated over a range of $p_T^{trig}$.

The $p_{out}$ distributions are shown in several bins of $x_E$, integrated over a range of $p_T^{trig}$.

The $p_{out}$ distributions are shown in several bins of $x_E$, integrated over a range of $p_T^{trig}$.

The $p_{out}$ distributions are shown in several bins of $x_E$, integrated over a range of $p_T^{trig}$.

The $p_{out}$ distributions are shown in several bins of $x_E$, integrated over a range of $p_T^{trig}$.

The $p_{out}$ distributions are shown in several bins of $x_E$, integrated over a range of $p_T^{trig}$.

The Gaussian widths of $p_{out}$ as a function of $x_E$ are shown for both $\pi^0$ and direct photon triggered correlations.

The nonperturbative Gaussian widths are shown as a function of $x_E$ for $\sqrt{s}$ = 510 GeV. The points for $\sqrt{s}$ = 200 GeV are given in figure 8 and its associated table.

The nonperturbative Gaussian widths are shown as a function of $x_E$ for $\sqrt{s}$ = 510 GeV. The points for $\sqrt{s}$ = 200 GeV are given in figure 8 and its associated table.

Away-side jet widths from a Gaussian fit by $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta in Au+Au collisions.

Away-side $I_{AA}$ for the entire away-side $|\Delta \phi - \pi| < \pi /2$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

Away-side $I_{AA}$ for the entire away-side $|\Delta \phi - \pi| < \pi /2$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

Away-side $I_{AA}$ for a narrow “head” $|\Delta \phi - \pi| < \pi /6$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

Away-side $I_{AA}$ for a narrow “head” $|\Delta \phi - \pi| < \pi /6$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

Supplemental near-side jet widths from a Gaussian fit by $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta in $p+p$ collisions.

Supplemental near-side jet widths from a Gaussian fit by $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta in Au+Au collisions.

Supplemental near-side jet widths from a Gaussian fit by $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta in Au+Au collisions.

Supplemental away-side $I_{AA}$ for $|\Delta \phi - \pi| < \pi /3$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

Supplemental away-side $I_{AA}$ for $|\Delta \phi - \pi| < \pi /3$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

$p+p$ yield vs $p_T^a$ supplemental tables.

Au+Au yield vs $p_T^a$ supplemental tables.

Au+Au yield vs $p_T^a$ supplemental tables.

$J_{dA}$ versus $x^{frag}_{Au}$ for $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for different centrality classes.

$J_{dA}$ versus $x^{frag}_{Au}$ for $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for different centrality classes.

$J_{dA}$ versus $x^{frag}_{Au}$ for $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for different centrality classes.

$p_T$ versus $R_{dA}$ for $\pi^0$ and range of 3.0 < $\eta$ < 3.8.

PHENIX $\pi^0$ invariant cross-section for $p$+$p$ collisions for range 3.0 < $\eta$ < 3.8.

Charged hadron $dN_{ch}/d\eta$ as a function of pseudorapidity in high-multiplicity 0%-5% central $p$+Al collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Charged hadron $dN_{ch}/d\eta$ as a function of pseudorapidity in various mutiplicity classes of $^3$He+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Charged hadron $dN_{ch}/d\eta$ as a function of pseudorapidity in various mutiplicity classes of $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Charged hadron $dN_{ch}/d\eta$ as a function of pseudorapidity in various mutiplicity classes of $p$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Charged hadron $dN_{ch}/d\eta$ as a function of pseudorapidity in various mutiplicity classes of $p$+Al collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Midrapidity charged hadron $dN_{ch}/d\eta$ per participating quark pair ($N_{qp}$/2) as a function of the number of participating quarks ($N_{qp}$).

Midrapidity charged hadron $dN_{ch}/d\eta$ per participating quark pair ($N_{qp}$/2) as a function of the number of participating quarks ($N_{qp}$).

Midrapidity charged hadron $dN_{ch}/d\eta$ per participating quark pair ($N_{qp}$/2) as a function of the number of participating quarks ($N_{qp}$).

Midrapidity charged hadron $dN_{ch}/d\eta$ per participating quark pair ($N_{qp}$/2) as a function of the number of participating quarks ($N_{qp}$).

Midrapidity charged hadron $dN_{ch}/d\eta$ per participating quark pair ($N_{qp}$/2) as a function of the number of participating quarks ($N_{qp}$).

Elliptic flow $v_2$ as a function of pseudorapidity in high-multiplicity 0%-5% central $p$+Al, $p$+Au, $d$+Au, and $^3$He+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

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