Nuclear modification factor $R_{AA}$ in 10% most central $Au$+$Au$ collisions.

Nuclear modification factor $R_{dA}$ for ($h^+$+$h^-$)/2 in minimum bias $d$+$Au$.

$R_{AA}$ in the 10% most central $Au$+$Au$ collisions.

$R_{dA}$ for ($h^+$+$h^-$)/2 measurement in $d$+$Au$.

$R_{dA}$ for $\pi^0$ measurement in $d$+$Au$.

Ratio of $R_{dA}$ in minimum bias $d$+$Au$ and $R_{pA}$ from neutron tagged $d$+$Au$ collisions for ($h^+$+$h^-$)/2.

Ratio of $R_{dA}$ in minimum bias $d$+$Au$ and $R_{pA}$ from neutron tagged $d$+$Au$ collisions for $\pi^0$

FIG. 1. (Color online) Perspective views of $2 \mathrm{D}$ charge-independent angular correlations $\Delta \rho / \sqrt{\rho_{\mathrm{ref}}}$ on $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ for Au-Au collisions at $\sqrt{s_{N N}}=200$ and $62 \mathrm{GeV}$ (top and bottom rows, respectively). Centrality increases left to right from most peripheral to most central. Corrected total cross-section fractions are (left to right) $84 \%-93 \%, 55 \%-64 \%, 18 \%-28 \%,$ and $0 \%-5 \%$ for the $200-\mathrm{GeV}$ data and $84 \%-95 \%, 56 \%-65 \%$ $18 \%-28 \%,$ and $0 \%-5 \%$ for the $62 \mathrm{GeV}$ data (see Tables III and IV).

FIG. 1. (Color online) Perspective views of $2 \mathrm{D}$ charge-independent angular correlations $\Delta \rho / \sqrt{\rho_{\mathrm{ref}}}$ on $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ for Au-Au collisions at $\sqrt{s_{N N}}=200$ and $62 \mathrm{GeV}$ (top and bottom rows, respectively). Centrality increases left to right from most peripheral to most central. Corrected total cross-section fractions are (left to right) $84 \%-93 \%, 55 \%-64 \%, 18 \%-28 \%,$ and $0 \%-5 \%$ for the $200-\mathrm{GeV}$ data and $84 \%-95 \%, 56 \%-65 \%$ $18 \%-28 \%,$ and $0 \%-5 \%$ for the $62 \mathrm{GeV}$ data (see Tables III and IV).

FIG. 1. (Color online) Perspective views of $2 \mathrm{D}$ charge-independent angular correlations $\Delta \rho / \sqrt{\rho_{\mathrm{ref}}}$ on $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ for Au-Au collisions at $\sqrt{s_{N N}}=200$ and $62 \mathrm{GeV}$ (top and bottom rows, respectively). Centrality increases left to right from most peripheral to most central. Corrected total cross-section fractions are (left to right) $84 \%-93 \%, 55 \%-64 \%, 18 \%-28 \%,$ and $0 \%-5 \%$ for the $200-\mathrm{GeV}$ data and $84 \%-95 \%, 56 \%-65 \%$ $18 \%-28 \%,$ and $0 \%-5 \%$ for the $62 \mathrm{GeV}$ data (see Tables III and IV).

FIG. 1. (Color online) Perspective views of $2 \mathrm{D}$ charge-independent angular correlations $\Delta \rho / \sqrt{\rho_{\mathrm{ref}}}$ on $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ for Au-Au collisions at $\sqrt{s_{N N}}=200$ and $62 \mathrm{GeV}$ (top and bottom rows, respectively). Centrality increases left to right from most peripheral to most central. Corrected total cross-section fractions are (left to right) $84 \%-93 \%, 55 \%-64 \%, 18 \%-28 \%,$ and $0 \%-5 \%$ for the $200-\mathrm{GeV}$ data and $84 \%-95 \%, 56 \%-65 \%$ $18 \%-28 \%,$ and $0 \%-5 \%$ for the $62 \mathrm{GeV}$ data (see Tables III and IV).

FIG. 1. (Color online) Perspective views of $2 \mathrm{D}$ charge-independent angular correlations $\Delta \rho / \sqrt{\rho_{\mathrm{ref}}}$ on $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ for Au-Au collisions at $\sqrt{s_{N N}}=200$ and $62 \mathrm{GeV}$ (top and bottom rows, respectively). Centrality increases left to right from most peripheral to most central. Corrected total cross-section fractions are (left to right) $84 \%-93 \%, 55 \%-64 \%, 18 \%-28 \%,$ and $0 \%-5 \%$ for the $200-\mathrm{GeV}$ data and $84 \%-95 \%, 56 \%-65 \%$ $18 \%-28 \%,$ and $0 \%-5 \%$ for the $62 \mathrm{GeV}$ data (see Tables III and IV).

FIG. 2. (Color online) Fit decomposition of the $46 \%-56 \%$ centrality data for 62 -GeV Au-Au collisions. The top panels show from left to right the corrected data, model fit, fit residuals (data - model), and same-side 2D Gaussian. The bottom panels similarly show the away-side azimuth dipole, the nonjet azimuth quadrupole, the 1D $\eta_{\Delta}$ Gaussian, and the 2D exponential. This centrality is just below the sharp transition at $v_{\text {trans }}=3$. Fit residuals (c) are scaled up eightfold relative to the data.

FIG. 3. Fit parameters for $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ correlation data from Au-Au collisions at $\sqrt{s_{N N}}=62$ (open symbols) and 200 GeV (solid symbols) versus centrality measure $\nu$ computed at fixed energy $(200 \mathrm{GeV})$. The SS $2 \mathrm{D}$ Gaussian amplitudes, $\eta_{\Delta}$ widths, and $\phi_{\Delta}$ widths are shown in the left, center, and right panels, respectively of the top row. The bottom row shows from left to right the amplitudes for the dipole, quadrupole, and SS peak width aspect ratio $\sigma_{\eta_{\Delta}} / \sigma_{\phi_{\Delta}} .$ Fitting errors are indicated by error bars where larger than the symbols. Solid lines connect the points for clarity. The dotted and dashed curves indicate Glauber linear superposition estimates for 62 - and 200 -GeV peak amplitudes respectively, as discussed in the text. The quadrupole data are consistent with Ref. [60]. The hatched regions indicate the full range of systematic uncertainties listed in Appendix F. The vertical dark bands indicate estimated $v$ equivalents for $N-N$ collisions and $b=0$ Au-Au collisions.

FIG. 3. Fit parameters for $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ correlation data from Au-Au collisions at $\sqrt{s_{N N}}=62$ (open symbols) and 200 GeV (solid symbols) versus centrality measure $\nu$ computed at fixed energy $(200 \mathrm{GeV})$. The SS $2 \mathrm{D}$ Gaussian amplitudes, $\eta_{\Delta}$ widths, and $\phi_{\Delta}$ widths are shown in the left, center, and right panels, respectively of the top row. The bottom row shows from left to right the amplitudes for the dipole, quadrupole, and SS peak width aspect ratio $\sigma_{\eta_{\Delta}} / \sigma_{\phi_{\Delta}} .$ Fitting errors are indicated by error bars where larger than the symbols. Solid lines connect the points for clarity. The dotted and dashed curves indicate Glauber linear superposition estimates for 62 - and 200 -GeV peak amplitudes respectively, as discussed in the text. The quadrupole data are consistent with Ref. [60]. The hatched regions indicate the full range of systematic uncertainties listed in Appendix F. The vertical dark bands indicate estimated $v$ equivalents for $N-N$ collisions and $b=0$ Au-Au collisions.

FIG. 3. Fit parameters for $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ correlation data from Au-Au collisions at $\sqrt{s_{N N}}=62$ (open symbols) and 200 GeV (solid symbols) versus centrality measure $\nu$ computed at fixed energy $(200 \mathrm{GeV})$. The SS $2 \mathrm{D}$ Gaussian amplitudes, $\eta_{\Delta}$ widths, and $\phi_{\Delta}$ widths are shown in the left, center, and right panels, respectively of the top row. The bottom row shows from left to right the amplitudes for the dipole, quadrupole, and SS peak width aspect ratio $\sigma_{\eta_{\Delta}} / \sigma_{\phi_{\Delta}} .$ Fitting errors are indicated by error bars where larger than the symbols. Solid lines connect the points for clarity. The dotted and dashed curves indicate Glauber linear superposition estimates for 62 - and 200 -GeV peak amplitudes respectively, as discussed in the text. The quadrupole data are consistent with Ref. [60]. The hatched regions indicate the full range of systematic uncertainties listed in Appendix F. The vertical dark bands indicate estimated $v$ equivalents for $N-N$ collisions and $b=0$ Au-Au collisions.

FIG. 3. Fit parameters for $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ correlation data from Au-Au collisions at $\sqrt{s_{N N}}=62$ (open symbols) and 200 GeV (solid symbols) versus centrality measure $\nu$ computed at fixed energy $(200 \mathrm{GeV})$. The SS $2 \mathrm{D}$ Gaussian amplitudes, $\eta_{\Delta}$ widths, and $\phi_{\Delta}$ widths are shown in the left, center, and right panels, respectively of the top row. The bottom row shows from left to right the amplitudes for the dipole, quadrupole, and SS peak width aspect ratio $\sigma_{\eta_{\Delta}} / \sigma_{\phi_{\Delta}} .$ Fitting errors are indicated by error bars where larger than the symbols. Solid lines connect the points for clarity. The dotted and dashed curves indicate Glauber linear superposition estimates for 62 - and 200 -GeV peak amplitudes respectively, as discussed in the text. The quadrupole data are consistent with Ref. [60]. The hatched regions indicate the full range of systematic uncertainties listed in Appendix F. The vertical dark bands indicate estimated $v$ equivalents for $N-N$ collisions and $b=0$ Au-Au collisions.

FIG. 3. Fit parameters for $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ correlation data from Au-Au collisions at $\sqrt{s_{N N}}=62$ (open symbols) and 200 GeV (solid symbols) versus centrality measure $\nu$ computed at fixed energy $(200 \mathrm{GeV})$. The SS $2 \mathrm{D}$ Gaussian amplitudes, $\eta_{\Delta}$ widths, and $\phi_{\Delta}$ widths are shown in the left, center, and right panels, respectively of the top row. The bottom row shows from left to right the amplitudes for the dipole, quadrupole, and SS peak width aspect ratio $\sigma_{\eta_{\Delta}} / \sigma_{\phi_{\Delta}} .$ Fitting errors are indicated by error bars where larger than the symbols. Solid lines connect the points for clarity. The dotted and dashed curves indicate Glauber linear superposition estimates for 62 - and 200 -GeV peak amplitudes respectively, as discussed in the text. The quadrupole data are consistent with Ref. [60]. The hatched regions indicate the full range of systematic uncertainties listed in Appendix F. The vertical dark bands indicate estimated $v$ equivalents for $N-N$ collisions and $b=0$ Au-Au collisions.

FIG. 3. Fit parameters for $\left(\eta_{\Delta}, \phi_{\Delta}\right)$ correlation data from Au-Au collisions at $\sqrt{s_{N N}}=62$ (open symbols) and 200 GeV (solid symbols) versus centrality measure $\nu$ computed at fixed energy $(200 \mathrm{GeV})$. The SS $2 \mathrm{D}$ Gaussian amplitudes, $\eta_{\Delta}$ widths, and $\phi_{\Delta}$ widths are shown in the left, center, and right panels, respectively of the top row. The bottom row shows from left to right the amplitudes for the dipole, quadrupole, and SS peak width aspect ratio $\sigma_{\eta_{\Delta}} / \sigma_{\phi_{\Delta}} .$ Fitting errors are indicated by error bars where larger than the symbols. Solid lines connect the points for clarity. The dotted and dashed curves indicate Glauber linear superposition estimates for 62 - and 200 -GeV peak amplitudes respectively, as discussed in the text. The quadrupole data are consistent with Ref. [60]. The hatched regions indicate the full range of systematic uncertainties listed in Appendix F. The vertical dark bands indicate estimated $v$ equivalents for $N-N$ collisions and $b=0$ Au-Au collisions.

FIG. 4. (Color online) (a) 2D data histogram from $0 \%-5 \%$ central 200 GeV Au-Au collisions $[(0,0)$ bin suppressed]. (b) Additional model component fitted to the $2 \mathrm{D}$ histogram and consisting of an AS azimuth dipole modulated by function $F\left(\eta_{\Delta}\right)$ described in the text. (c) Fit residuals including the new model component. (d) Data histogram in the first panel minus the additional model component and fitted offset. Remaining structure is described accurately by an AS dipole and SS 2D Gaussian.

FIG. 9. (a) Amplitude of SS 2D peak (jet-correlated pairs per final-state hadron) for data (solid points), GLS extrapolation from $p$ - $p$ data (dashed curve) and from the HIJING Monte Carlo (open points, dotted curve). (b) Single-particle production extrapolated from $p-p$ data (dashed curve, GLS), and from more-central Au-Au data (solid line) and HIJING (dash-dotted line, open points). (c) Single-particle hard-component production per $N-N$ binary collision from $p-p$ data, from more-central Au-Au data and from HIJING (open points). (d) Jet-correlated pair production (SS 2D peak) per $N-N$ binary collision extrapolated from $p-p$ data (GLS), and from Au-Au data (solid points) and HIJING (open points).

FIG. 9. (a) Amplitude of SS 2D peak (jet-correlated pairs per final-state hadron) for data (solid points), GLS extrapolation from $p$ - $p$ data (dashed curve) and from the HIJING Monte Carlo (open points, dotted curve). (b) Single-particle production extrapolated from $p-p$ data (dashed curve, GLS), and from more-central Au-Au data (solid line) and HIJING (dash-dotted line, open points). (c) Single-particle hard-component production per $N-N$ binary collision from $p-p$ data, from more-central Au-Au data and from HIJING (open points). (d) Jet-correlated pair production (SS 2D peak) per $N-N$ binary collision extrapolated from $p-p$ data (GLS), and from Au-Au data (solid points) and HIJING (open points).

FIG. 9. (a) Amplitude of SS 2D peak (jet-correlated pairs per final-state hadron) for data (solid points), GLS extrapolation from $p$ - $p$ data (dashed curve) and from the HIJING Monte Carlo (open points, dotted curve). (b) Single-particle production extrapolated from $p-p$ data (dashed curve, GLS), and from more-central Au-Au data (solid line) and HIJING (dash-dotted line, open points). (c) Single-particle hard-component production per $N-N$ binary collision from $p-p$ data, from more-central Au-Au data and from HIJING (open points). (d) Jet-correlated pair production (SS 2D peak) per $N-N$ binary collision extrapolated from $p-p$ data (GLS), and from Au-Au data (solid points) and HIJING (open points).

FIG. 9. (a) Amplitude of SS 2D peak (jet-correlated pairs per final-state hadron) for data (solid points), GLS extrapolation from $p$ - $p$ data (dashed curve) and from the HIJING Monte Carlo (open points, dotted curve). (b) Single-particle production extrapolated from $p-p$ data (dashed curve, GLS), and from more-central Au-Au data (solid line) and HIJING (dash-dotted line, open points). (c) Single-particle hard-component production per $N-N$ binary collision from $p-p$ data, from more-central Au-Au data and from HIJING (open points). (d) Jet-correlated pair production (SS 2D peak) per $N-N$ binary collision extrapolated from $p-p$ data (GLS), and from Au-Au data (solid points) and HIJING (open points).

FIG. 11. (Color online) (Left) Uncorrected angular correlations from 62-GeV $37 \%-46 \%$ central Au-Au collisions showing pileup distortions, especially evident as the W-shaped nonuniformity of the AS ridge on $\eta_{\Delta}$. (Right) The same data with pileup correction applied.

$I_{dAu}$ vs. $p_T^{assoc}$ for different centrality, $p_T^{trig}$ and $\eta^{trig}$ bins.

The ratio of correlation peak widths between $d + Au$ and $p + p$ collisions. Only statistic errors are shown.

The ratio of correlation peak widths between $d + Au$ and $p + p$ collisions. Only statistic errors are shown.

Upper panel, Azimuthal distributions of associated particles for trigger particles in-plane (squares) and out-of-plane (triangles) for Au+Au collisions at centrality 20-60%. Open symbols are reflections of solid symbols around $\Delta \phi$ = 0 and $\Delta \phi$ = $\pi$. Elliptic flow contribution is shown by dashed lines. Lower panel, Distributions after substracting elliptic flow, and the corresponding measurement in p + p collisions (histogram).

$v_{2}$ at 3 ≤ $p_{t}$ ≤ 6 GeV/c versus impact parameter, b, compared to models of particle emission by a static source (see text).

Charged hadron v2 for the centrality bin 0-50%.

v2(4) vs. p_T for pions for 20 - 60% centrality

v2(4) vs. p_T for kaons for 20 - 60% centrality

v2(4) vs. p_T for anti-protons for 20 - 60% centrality

v2(2) vs. centrality for pions

v2(2) vs. centrality for kaons

v2(2) vs. centrality for anti-protons

v2 vs. p_T for particles identified in the RICH (min. bias) for pions+kaons and protons/antiprotons

v2 for Kshort for min. bias

v2 for kaons from kinks for min. bias

v2 for kaons from dE/dx for min. bias

v2 for Lambda+Lambdabar for min. bias [PRL/92/052302/2004]. The PHENIX data are from [PRL/91/182301/2003]. The hydro calculations are from [P. Huovinen/private communication].

v2 for charged pions for min. bias. kaon data points are from Figure 9.

v2 for anti-protons for min. bias. kaon data points are from Figure 9.

Azimuthal correlations (\sum_i cos(2(\phi_pT -\phi_i))) in Au+Au collisions as a function of p_T along with results from p+p collisions.

Azimuthal correlations (\sum_i cos(2(\phi_pT -\phi_i))) from d+Au as a function of p_T for min bias collisions.

Azimuthal correlations (\sum_i cos(2(\phi_pT -\phi_i))) from d+Au as a function of p_T for three centralities.

Charged hadron v2 and v4{EP_2} as a function of pseudorapidity for minimum bias.

v2(2) and v2(4) vs. pseudorapidity

v2(2) and v2(4) vs. pseudorapidity,FTPC only

v2(4) 20-70% in TPC+FTPC and v2(2) 20-70% in FTPC.

MC v2 with FTPC momentum resolution.

Charged hadron v2 for the centrality bin 0-5% in standard and modified event plane methods.

Charged hadron v2 for the centrality bin 5-10% in standard and modified event plane methods.

Charged hadron v2 for the centrality bin 10-20% in standard and modified event plane methods.

Charged hadron v2 for the centrality bin 20-30% in standard and modified event plane methods.

Charged hadron v2 for the centrality bin 30-40% in standard and modified event plane methods.

Charged hadron v2 for the centrality bin 40-60% in standard and modified event plane methods.

Charged hadron v2{2} from modified reaction plane method and two-particle cumulant results (after subtracting the correlations measured in p+p collisions) for centrality bin 20-60% .

v4(EP2) for charged hadrons min. bias

v4(EP2) vs. centrality for charged hadrons

v4(EP2) for pions in min. bias

v4(EP2) for protons in min. bias

v4(EP2) for kaons in min. bias

v4(EP2) for Lambda in min. bias

Charged hadron v4{EP_2} as a function of pseudorapidity for minimum bias (TPC only).

charged hadrons high p_T (3 - 6 GeV/c) v4(3)

charged hadron integrated v2 as a function of centrality for the different methods

charged hadron v2(6) vs. centrality bin. v2(2) and v2(4) are from Figure 29(a).

g2 as a function of centrality from 200 GeV, 130 GeV,NA49 and q-dist methods. The solid points are from the cumulant method, while open circles are from the q distribution method. Please note that centrality percentage differs from different method. For 200 GeV-cumulant (0, 64.1, 53.9,44.7,35.2,25.4,15.1,7.7,0), for 200 GeV-q-dist (73.8,64.1,53.9,44.7,35.2,25.4,15.1,7.7,2.3), for 130 GeV (65,47,36,27.5,20,13,7.5,2.5,0) and for 17 GeV (0,0,0,61.2,38.2,29.1,18.2,8.5,3.8).

wounded nucleons as a function of centrality from 200 GeV, 130 GeV,NA49 and q-dist methods.Please note that centrality percentage differs from different method. For 200 GeV-cumulant (0, 64.1, 53.9,44.7,35.2,25.4,15.1,7.7,0), for 200 GeV-q-dist (73.8,64.1,53.9,44.7,35.2,25.4,15.1,7.7,2.3), for 130 GeV (65,47,36,27.5,20,13,7.5,2.5,0) and for 17 GeV (0,0,0,61.2,38.2,29.1,18.2,8.5,3.8).

v2/n for pions for three centrality bins

v2/n for protons for three centrality bins

v2/n for Kshort for three centrality bins

v2/n for Lambda + Lambda bar for three centrality bins

Blast Wave parameters rho_2 from pion v2(2),pion v2 and hadron v2

Blast Wave parameters rho_4 from hadron v2

Blast Wave parameters epsilon from pion v2(2),pion v2, hadron v2 and initial

Blast Wave parameters rho_0 from pion v2(2),pion v2 and hadron v2

Blast Wave parameters T from pion v2(2),pion v2 and hadron v2

Slopes and intercepts of $<p_x>/<p_T>(\eta)$ in Cu+Au and Au+Au collisions.

Slopes and intercepts of $v_1(\eta)$ in Cu+Au and Au+Au collisions.

Centrality dependence of $v_1^{fluc(even)}$ and $<p_x>/<p_T>$ in Cu+Au and Au+Au collisions.

$p_T$ dependence of $v_1^{fluc(even)}$ and $v_1^{conv(odd)}$ in Cu+Au and Au+Au collisions.

$p_T$ dependence of $v_1$ measured with three-point correlator in Cu+Au collisions

$p_T$ dependence of $v_1$ measured with three-point correlator in Cu+Au collisions

$\eta$ dependence of $v_1$ measured with three-point correlator in Cu+Au collisions

Directed flow $v_1(p_T)$ of $\pi^+ + \pi^-$ in Cu+Au collisions.

Directed flow $v_1(p_T)$ of $K^++K^-$ in Cu+Au collisions.

Directed flow $v_1(p_T)$ of $p+\bar{p}$ in Cu+Au collisions.

Centrality dpendence of $<p_x>$ of positively and negatively charged particles and its differece in Cu+Au and Au+Au collisions.

Higher harmonic flow coefficients $v_n{\Psi_n}$ of charged particles for different centrality bins in Cu+Au collisions.

Higher harmonic flow coefficients of charged particles for different centrality bins in Cu+Au collisions.

Higher harmonic flow coefficients $v_n{\Psi_n}$ of charged particles for different centrality bins in Cu+Au collisions.

Higher harmonic flow coefficients $v_n$ of charged particles for two selected $p_T$ bins as a function of the number of participants calculated with a Monte-Carlo Glauber simulationin Cu+Au collisions.

Higher harmonic flow coefficients $v_n$ of charged particles for two selected $p_T$ bins as a function of the number of participants calculated with a Monte-Carlo Glauber simulationin Cu+Au collisions.

The second and third harmonic flow coefficients $v_n(p_T)$ of charged particles measured with the scalar product method for 0-5% and 20-30% centrality bins in Cu+Au collisions.

The second $v_2(p_T)$ harmonic flow coefficients of $\pi$ for different centrality bins in Cu+Au collisions.

The second $v_2(p_T)$ harmonic flow coefficients of p for different centrality bins in Cu+Au collisions.

The third $v_3(p_T)$ harmonic flow coefficients of $\pi$ for different centrality bins in Cu+Au collisions.

The third $v_3(p_T)$ harmonic flow coefficients of p for different centrality bins in Cu+Au collisions.

The second $v_2(p_T)$ and third $v_3(p_T)$ harmonic flow coefficients of K for different centrality bins in Cu+Au collisions.

Higher harmonic flow coefficients $v_2(p_T)$, $v_3(p_T)$, and $v_4(p_T)$ of $\pi$, K, and p for 0-40% centrality bin in Cu+Au collisions.

Higher harmonic flow coefficients $v_2(p_T)$, $v_3(p_T)$, and $v_4(p_T)$ of $\pi$, K, and p for 0-40% centrality bin in Cu+Au collisions.

Higher harmonic flow coefficients $v_2(p_T)$, $v_3(p_T)$, and $v_4(p_T)$ of $\pi$, K, and p for 0-40% centrality bin in Cu+Au collisions.

Higher harmonic flow coefficients $v_2(p_T)$, $v_3(p_T)$, and $v_4(p_T)$ of $\pi$, K, and p for 0-40% centrality bin in Cu+Au collisions.

Higher harmonic flow coefficients $v_2(p_T)$, $v_3(p_T)$, and $v_4(p_T)$ of $\pi$, K, and p for 0-40% centrality bin in Cu+Au collisions.

Higher harmonic flow coefficients $v_2(p_T)$, $v_3(p_T)$, and $v_4(p_T)$ of $\pi$, K, and p for 0-40% centrality bin in Cu+Au collisions.

Higher harmonic flow coefficients $v_2(p_T)$, $v_3(p_T)$, and $v_4(p_T)$ of $\pi$, K, and p for 0-40% centrality bin in Cu+Au collisions.

Higher harmonic flow coefficients $v_2(p_T)$, $v_3(p_T)$, and $v_4(p_T)$ of $\pi$, K, and p for 0-40% centrality bin in Cu+Au collisions.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

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