Nuclear modicifation factor vs. N_coll $J\psi$ with|y|<1 and pT < 3 GeV/c in d Au collisions.

The invariant cross section versus transverse momentum for J/$\psi$ with $|y|\lt 1$ in $p+p$ collisions.

The invariant yield versus transverse momentum for $J/\Psi$ with |y| < 1 in 0 - 100 percent central d+Au collisions.

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.

The three-point correlator, $\gamma$, as a function of centrality for Au+Au collisions at 19.6 GeV.

The three-point correlator, $\gamma$, as a function of centrality for Au+Au collisions at 11.5 GeV.

The three-point correlator, $\gamma$, as a function of centrality for Au+Au collisions at 7.7.

The two-particle correlation as a function of centrality for Au+Au collisions at 62.4 GeV.

The two-particle correlation as a function of centrality for Au+Au collisions at 39 GeV.

The two-particle correlation as a function of centrality for Au+Au collisions at 27 GeV.

The two-particle correlation as a function of centrality for Au+Au collisions at 19.6 GeV.

The two-particle correlation as a function of centrality for Au+Au collisions at 11.5 GeV.

The two-particle correlation as a function of centrality for Au+Au collisions at 7.7 GeV.

$H_{SS}-H{OS}$, as a function of beam energy for 60-80% centrality in Au+Au collisions.

$H_{SS}-H{OS}$, as a function of beam energy for 30-60% centrality in Au+Au collisions.

$H_{SS}-H{OS}$, as a function of beam energy for 10-30% centrality in Au+Au collisions.

The $v_{dyn}$ and the three terms of $v_{dyn}$ vs $\sqrt{\langle N_{ch}\rangle \langle N_{\gamma}\rangle }$ for mixed events. $Corr_{\gamma\,-ch}^{real}$ is plotted.

The $v_{dyn}$ and the three terms of $v_{dyn}$ vs $\sqrt{\langle N_{ch}\rangle \langle N_{\gamma}\rangle }$ for mixed events. $Corr_{\gamma\,-ch}^{mixed}$ is plotted.

The $v_{dyn}$ and the three terms of $v_{dyn}$ vs $\sqrt{\langle N_{ch}\rangle \langle N_{\gamma}\rangle }$ for mixed events. $\omega_{\gamma}^{real}$ is plotted.

The $v_{dyn}$ and the three terms of $v_{dyn}$ vs $\sqrt{\langle N_{ch}\rangle \langle N_{\gamma}\rangle }$ for mixed events. $\omega_{\gamma}^{mixed}$ is plotted.

The $v_{dyn}$ and the three terms of $v_{dyn}$ vs $\sqrt{\langle N_{ch}\rangle \langle N_{\gamma}\rangle }$ for mixed events. $v_{dyn}^{real}$ is plotted.

The $v_{dyn}$ and the three terms of $v_{dyn}$ vs $\sqrt{\langle N_{ch}\rangle \langle N_{\gamma}\rangle }$ for mixed events. $v_{dyn}^{mixed}$ is plotted.

The values of the observable $v_{dyn}^{ch-\gamma}$ for the same side.

The values of the observable $v_{dyn}^{ch-\gamma}$ for the same side.

The values of the observable $v_{dyn}^{ch-\gamma}$ for the away-side.

The values of the observable $v_{dyn}^{ch-\gamma}$ for the away-side.

The correlation between positive and negative charged particles measured by the FTPC and photons measured by the PMD using $v_{dyn}$ for real events.

$r_{m,1}$ vs multiplicity for first order of m.

$r_{m,1}$ vs multiplicity for first order of m.

$r_{m,1}$ vs multiplicity for second order of m.

$r_{m,1}$ vs multiplicity for second order of m.

$r_{m,1}$ vs multiplicity for third order of m.

$r_{m,1}$ vs multiplicity for third order of m.

The moments $r_{m,1}$(m=1-3) for real and mixed events as a function of order m in a fixed multiplicity bin of 47$\langle \sqrt{\langle N_{ch}\rangle \langle N_{\gamma}\rangle }\langle 54$.

K$^+$p (K$^+$p $\oplus$ K$^-\overline{\mathrm p}$) correlation function in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV (40-100%).

K$^+$p (K$^+$p $\oplus$ K$^-\overline{\mathrm p}$) correlation function in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV (60-70%).

K$^+$p (K$^+$p $\oplus$ K$^-\overline{\mathrm p}$) correlation function in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV (70-80%).

K$^+$p (K$^+$p $\oplus$ K$^-\overline{\mathrm p}$) correlation function in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV (80-90%).

K$^-$p (K$^-$p $\oplus$ K$^+\overline{\mathrm p}$) correlation function in pp collisions at $\sqrt{s}=13$ TeV.

K$^-$p (K$^-$p $\oplus$ K$^+\overline{\mathrm p}$) correlation function in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV (0-20%).

K$^-$p (K$^-$p $\oplus$ K$^+\overline{\mathrm p}$) correlation function in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV (20-40%).

K$^-$p (K$^-$p $\oplus$ K$^+\overline{\mathrm p}$) correlation function in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV (40-100%).

K$^-$p (K$^-$p $\oplus$ K$^+\overline{\mathrm p}$) correlation function in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV (60-70%).

K$^-$p (K$^-$p $\oplus$ K$^+\overline{\mathrm p}$) correlation function in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV (70-80%).

K$^-$p (K$^-$p $\oplus$ K$^+\overline{\mathrm p}$) correlation function in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV (80-90%).

Scaling factor ($\alpha_j $) for $\mathrm{\overline{K}^0 n}$ extracted from the different fits of the K$^-$ p correlation function as a function of the core radius $r_{\mathrm{core}}$.

Scaling factor ($\alpha_j $) for $\pi \Sigma$ extracted from the different fits of the K$^-$ p correlation function as a function of the core radius $r_{\mathrm{core}}$.

This table shows the correlation coefficient between the $\alpha_{\mathrm{\overline{K}^0 n}}$ and $\alpha_{\pi \Sigma}$ parameters extracted from the different fits to the K$^-$ p correlation function with free weights. The results are shown as a function of the core radius $r_{\mathrm{core}}$ shown in Table 7. The error on the correlation coefficient, shown as statistical errors in the table, are obteined by using the Pearson's formula to evaluate the probable error of correlation coefficient.

Correlation Asymmetry data from the upgraded detector.

Combined correlation data.

Combined correlation asymmetry data.

V1Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-2% having the associated particle PT in the range 0.25-0.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V1Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-2% having the associated particle PT in the range 1-1.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V1Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-10% having the associated particle PT in the range 0.25-0.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V1Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-10% having the associated particle PT in the range 1-1.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V1Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-10% having the associated particle PT in the range 2-2.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V1Delta coefficients as a function of the trigger particle PT for events in the centrality class 40-50% having the associated particle PT in the range 0.25-0.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V1Delta coefficients as a function of the trigger particle PT for events in the centrality class 40-50% having the associated particle PT in the range 1-1.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V1Delta coefficients as a function of the trigger particle PT for events in the centrality class 40-50% having the associated particle PT in the range 2-2.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V2Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-2% having the associated particle PT in the range 0.25-0.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V2Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-2% having the associated particle PT in the range 1-1.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V2Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-10% having the associated particle PT in the range 0.25-0.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V2Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-10% having the associated particle PT in the range 1-1.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V2Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-10% having the associated particle PT in the range 2-2.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V2Delta coefficients as a function of the trigger particle PT for events in the centrality class 40-50% having the associated particle PT in the range 0.25-0.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V2Delta coefficients as a function of the trigger particle PT for events in the centrality class 40-50% having the associated particle PT in the range 1-1.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V2Delta coefficients as a function of the trigger particle PT for events in the centrality class 40-50% having the associated particle PT in the range 2-2.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V3Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-2% having the associated particle PT in the range 0.25-0.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V3Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-2% having the associated particle PT in the range 1-1.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V3Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-10% having the associated particle PT in the range 0.25-0.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V3Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-10% having the associated particle PT in the range 1-1.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V3Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-10% having the associated particle PT in the range 2-2.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V3Delta coefficients as a function of the trigger particle PT for events in the centrality class 40-50% having the associated particle PT in the range 0.25-0.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V3Delta coefficients as a function of the trigger particle PT for events in the centrality class 40-50% having the associated particle PT in the range 1-1.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V3Delta coefficients as a function of the trigger particle PT for events in the centrality class 40-50% having the associated particle PT in the range 2-2.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V4Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-2% having the associated particle PT in the range 0.25-0.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V4Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-2% having the associated particle PT in the range 1-1.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V4Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-10% having the associated particle PT in the range 0.25-0.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V4Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-10% having the associated particle PT in the range 1-1.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V4Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-10% having the associated particle PT in the range 2-2.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V4Delta coefficients as a function of the trigger particle PT for events in the centrality class 40-50% having the associated particle PT in the range 0.25-0.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V4Delta coefficients as a function of the trigger particle PT for events in the centrality class 40-50% having the associated particle PT in the range 1-1.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V4Delta coefficients as a function of the trigger particle PT for events in the centrality class 40-50% having the associated particle PT in the range 2-2.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V5Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-2% having the associated particle PT in the range 0.25-0.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V5Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-2% having the associated particle PT in the range 1-1.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V5Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-10% having the associated particle PT in the range 0.25-0.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V5Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-10% having the associated particle PT in the range 1-1.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V5Delta coefficients as a function of the trigger particle PT for events in the centrality class 0-10% having the associated particle PT in the range 2-2.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V5Delta coefficients as a function of the trigger particle PT for events in the centrality class 40-50% having the associated particle PT in the range 0.25-0.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V5Delta coefficients as a function of the trigger particle PT for events in the centrality class 40-50% having the associated particle PT in the range 1-1.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V5Delta coefficients as a function of the trigger particle PT for events in the centrality class 40-50% having the associated particle PT in the range 2-2.5 GeV. Note that in the paper the data are plotted multiplied by 100.

V1Delta for all pt combinations for the centrality classes 0-10% and 40-50%.

V2Delta for all pt combinations for the centrality classes 0-10% and 40-50%.

V3Delta for all pt combinations for the centrality classes 0-10% and 40-50%.

V4Delta for all pt combinations for the centrality classes 0-10% and 40-50%.

V5Delta for all pt combinations for the centrality classes 0-10% and 40-50%.

The global fit paramter v_2{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 0-2%.

The global fit paramter v_2{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 2-10%.

The global fit paramter v_2{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 10-20%.

The global fit paramter v_2{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 20-30%.

The global fit paramter v_2{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 30-40%.

The global fit paramter v_2{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 40-50%.

The global fit paramter v_3{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 0-2%.

The global fit paramter v_3{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 2-10%.

The global fit paramter v_3{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 10-20%.

The global fit paramter v_3{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 20-30%.

The global fit paramter v_3{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 30-40%.

The global fit paramter v_3{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 40-50%.

The global fit paramter v_4{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 0-2%.

The global fit paramter v_4{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 2-10%.

The global fit paramter v_4{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 10-20%.

The global fit paramter v_4{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 20-30%.

The global fit paramter v_4{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 30-40%.

The global fit paramter v_4{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 40-50%.

The global fit paramter v_5{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 0-2%.

The global fit paramter v_5{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 2-10%.

The global fit paramter v_5{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 10-20%.

The global fit paramter v_5{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 20-30%.

The global fit paramter v_5{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 30-40%.

The global fit paramter v_5{GF} vs. the trigger particle PT for fit range of 0.25-15 GeV/c for the centrality class 40-50%.

The global fit paramter v_1{GF} vs. the trigger particle PT for fit range of 0.25-5 GeV/c for the centrality class 0-20%.

The global fit paramter v_1{GF} vs. the trigger particle PT for fit range of 0.25-5 GeV/c for the centrality class 40-50%.

The global fit paramter v_2{GF} vs. the trigger particle PT for fit range of 0.25-5 GeV/c for the centrality class 0-20%.

The global fit paramter v_2{GF} vs. the trigger particle PT for fit range of 0.25-5 GeV/c for the centrality class 40-50%.

The global fit paramter v_3{GF} vs. the trigger particle PT for fit range of 0.25-5 GeV/c for the centrality class 0-20%.

The global fit paramter v_3{GF} vs. the trigger particle PT for fit range of 0.25-5 GeV/c for the centrality class 40-50%.

The global fit paramter v_4{GF} vs. the trigger particle PT for fit range of 0.25-5 GeV/c for the centrality class 0-20%.

The global fit paramter v_4{GF} vs. the trigger particle PT for fit range of 0.25-5 GeV/c for the centrality class 40-50%.

The global fit paramter v_5{GF} vs. the trigger particle PT for fit range of 0.25-5 GeV/c for the centrality class 0-20%.

The global fit paramter v_5{GF} vs. the trigger particle PT for fit range of 0.25-5 GeV/c for the centrality class 40-50%.

The global fit paramter v_1{GF} vs. the trigger particle PT for fit range of 5-15 GeV/c for the centrality class 0-20%.

The global fit paramter v_1{GF} vs. the trigger particle PT for fit range of 5-15 GeV/c for the centrality class 40-50%.

The global fit paramter v_2{GF} vs. the trigger particle PT for fit range of 5-15 GeV/c for the centrality class 0-20%.

The global fit paramter v_2{GF} vs. the trigger particle PT for fit range of 5-15 GeV/c for the centrality class 40-50%.

The global fit paramter v_3{GF} vs. the trigger particle PT for fit range of 5-15 GeV/c for the centrality class 0-20%.

The global fit paramter v_3{GF} vs. the trigger particle PT for fit range of 5-15 GeV/c for the centrality class 40-50%.

The global fit paramter v_4{GF} vs. the trigger particle PT for fit range of 5-15 GeV/c for the centrality class 0-20%.

The global fit paramter v_4{GF} vs. the trigger particle PT for fit range of 5-15 GeV/c for the centrality class 40-50%.

The global fit paramter v_5{GF} vs. the trigger particle PT for fit range of 5-15 GeV/c for the centrality class 0-20%.

The global fit paramter v_5{GF} vs. the trigger particle PT for fit range of 5-15 GeV/c for the centrality class 40-50%.

The global fit paramter v_1{GF} vs. the trigger particle PT for fit range of 0.25-1 GeV/c for the centrality class 0-10%.

The global fit paramter v_1{GF} vs. the trigger particle PT for fit range of 0.25-1 GeV/c for the centrality class 40-50%.

The global fit paramter v_2{GF} vs. the trigger particle PT for fit range of 0.25-1 GeV/c for the centrality class 0-10%.

The global fit paramter v_2{GF} vs. the trigger particle PT for fit range of 0.25-1 GeV/c for the centrality class 40-50%.

The global fit paramter v_1{GF} vs. the trigger particle PT for fit range of 2-4 GeV/c for the centrality class 0-10%.

The global fit paramter v_1{GF} vs. the trigger particle PT for fit range of 2-4 GeV/c for the centrality class 40-50%.

The global fit paramter v_2{GF} vs. the trigger particle PT for fit range of 2-4 GeV/c for the centrality class 0-10%.

The global fit paramter v_2{GF} vs. the trigger particle PT for fit range of 2-4 GeV/c for the centrality class 40-50%.

FIG. 1: (a) Raw two-particle correlation signal $Y_2$ (red), background $aB_{inc}F_2$ (solid histogram), and background systematic uncertainty from a (dashed histograms). (b) Background-subtracted two-particle correlation $\hat{Y}_2$ (red), and systematic uncertainties due to a (dashed histograms) and flow (blue histograms). (c) Raw three-particle correlation $Y_3$. (d) $ba^2Y_{inc}^2$ . (e) Sum of trig-corr-bkgd and trigger flow. Data are from 12% central Au+Au collisions. Statistical errors in (a,b) are smaller than the point size.

FIG. 1: (a) Raw two-particle correlation signal $Y_2$ (red), background $aB_{inc}F_2$ (solid histogram), and background systematic uncertainty from a (dashed histograms). (b) Background-subtracted two-particle correlation $\hat{Y}_2$ (red), and systematic uncertainties due to a (dashed histograms) and flow (blue histograms). (c) Raw three-particle correlation $Y_3$. (d) $ba^2Y_{inc}^2$ . (e) Sum of trig-corr-bkgd and trigger flow. Data are from 12% central Au+Au collisions. Statistical errors in (a,b) are smaller than the point size.

FIG. 2: Background subtracted three-particle correlations, $\hat{Y}_3$, for (a) pp, (b) d+Au, (c) 80-50% Au+Au, (d) 50-30% Au+Au, (e) 30-10% Au+Au, and (f) 12% central Au+Au. Statistical errors per bin are approximately $\pm$0.012 in (a) and $\pm$0.006 in (b), both at ($\pi$, $\pi$), and are $\pm$0.022, $\pm$0.049, $\pm$0.099 and $\pm$0.077 from (c) to (f), similar for all bins.

FIG. 2: Background subtracted three-particle correlations, $\hat{Y}_3$, for (a) pp, (b) d+Au, (c) 80-50% Au+Au, (d) 50-30% Au+Au, (e) 30-10% Au+Au, and (f) 12% central Au+Au. Statistical errors per bin are approximately $\pm$0.012 in (a) and $\pm$0.006 in (b), both at ($\pi$, $\pi$), and are $\pm$0.022, $\pm$0.049, $\pm$0.099 and $\pm$0.077 from (c) to (f), similar for all bins.

FIG. 2: Background subtracted three-particle correlations, $\hat{Y}_3$, for (a) pp, (b) d+Au, (c) 80-50% Au+Au, (d) 50-30% Au+Au, (e) 30-10% Au+Au, and (f) 12% central Au+Au. Statistical errors per bin are approximately $\pm$0.012 in (a) and $\pm$0.006 in (b), both at ($\pi$, $\pi$), and are $\pm$0.022, $\pm$0.049, $\pm$0.099 and $\pm$0.077 from (c) to (f), similar for all bins.

FIG. 2: Background subtracted three-particle correlations, $\hat{Y}_3$, for (a) pp, (b) d+Au, (c) 80-50% Au+Au, (d) 50-30% Au+Au, (e) 30-10% Au+Au, and (f) 12% central Au+Au. Statistical errors per bin are approximately $\pm$0.012 in (a) and $\pm$0.006 in (b), both at ($\pi$, $\pi$), and are $\pm$0.022, $\pm$0.049, $\pm$0.099 and $\pm$0.077 from (c) to (f), similar for all bins.

FIG. 2: Background subtracted three-particle correlations, $\hat{Y}_3$, for (a) pp, (b) d+Au, (c) 80-50% Au+Au, (d) 50-30% Au+Au, (e) 30-10% Au+Au, and (f) 12% central Au+Au. Statistical errors per bin are approximately $\pm$0.012 in (a) and $\pm$0.006 in (b), both at ($\pi$, $\pi$), and are $\pm$0.022, $\pm$0.049, $\pm$0.099 and $\pm$0.077 from (c) to (f), similar for all bins.

FIG. 2: Background subtracted three-particle correlations, $\hat{Y}_3$, for (a) pp, (b) d+Au, (c) 80-50% Au+Au, (d) 50-30% Au+Au, (e) 30-10% Au+Au, and (f) 12% central Au+Au. Statistical errors per bin are approximately $\pm$0.012 in (a) and $\pm$0.006 in (b), both at ($\pi$, $\pi$), and are $\pm$0.022, $\pm$0.049, $\pm$0.099 and $\pm$0.077 from (c) to (f), similar for all bins.

Projections of away-side three-particle correlations along the diagonal $\Sigma$ within 0 < $\Delta$ < 0.35 (squares) and along the off-diagonal $\Delta$ within |$\Sigma$| < 0.35 (points) in (a) d+Au and (b) 12% central Au+Au collisions. The shaded areas indicate systematic uncertainties on the off-diagonal projections. The histogram in (a) is the near-side off-diagonal projection. The histogram in (b) is the away-side off-diagonal projection of our result with a = b = 1. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Projections of away-side three-particle correlations along the diagonal $\Sigma$ within 0 < $\Delta$ < 0.35 (squares) and along the off-diagonal $\Delta$ within |$\Sigma$| < 0.35 (points) in (a) d+Au and (b) 12% central Au+Au collisions. The shaded areas indicate systematic uncertainties on the off-diagonal projections. The histogram in (a) is the near-side off-diagonal projection. The histogram in (b) is the away-side off-diagonal projection of our result with a = b = 1. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Three-particle correlation strength per radian$^2$ versus ($N_{part}$/2)$^{1/3}$ where $N_{part}$ is the number of participants. Some of the data points have been displaced in $N_{part}$ for clarity.

Parton inclusive-jet $p_T$ and $A_{LL}$ values with associated uncertainties for jet-$\eta$ region $0.5<|\eta|<1$.

Parton inclusive-jet $p_T$ and $A_{LL}$ values with associated uncertainties for jet-$\eta$ region $|\eta|<0.5$.

Parton dijet invariant mass $M_{inv}$ and $A_{LL}$ values with associated uncertainties for the $\textrm{sign}(\eta_1) = \textrm{sign}(\eta_2)$ topology.

Parton dijet invariant mass $M_{inv}$ and $A_{LL}$ values with associated uncertainties for the $\textrm{sign}(\eta_1) \neq \textrm{sign}(\eta_2)$ topology.

Parton inclusive-jet $p_T$, $x_T$, and $A_{LL}$ values with associated uncertainties for jet-$\eta$ region $|\eta|<1$.

The correlation matrix for the point-to-point uncertainties for inclusive jet measurements with jets in the $0.5<|\eta|<1$ region. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties for inclusive jet measurements with jets in the $|\eta|<0.5$ region. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties for inclusive jet measurements with jets in the $|\eta|<1$ region. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $|\eta|<0.5$ region and jets in the $0.5<|\eta|<1$ region. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $0.5<|\eta|<1$ region with dijet measurements with the the $\textrm{sign}(\eta_1) = \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $0.5<|\eta|<1$ region with dijet measurements with the the $\textrm{sign}(\eta_1) \neq \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $|\eta|<0.5$ region with dijet measurements with the the $\textrm{sign}(\eta_1) = \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $|\eta|<0.5$ region with dijet measurements with the the $\textrm{sign}(\eta_1) \neq \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $|\eta|<1$ region with dijet measurements with the the $\textrm{sign}(\eta_1) = \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling inclusive jet measurements with jets in the $|\eta|<1$ region with dijet measurements with the the $\textrm{sign}(\eta_1) \neq \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties for dijet measurements with the $\textrm{sign}(\eta_1) = \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties for dijet measurements with the $\textrm{sign}(\eta_1) \neq \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.

The correlation matrix for the point-to-point uncertainties coupling dijet measurements with the $\textrm{sign}(\eta_1) \neq \textrm{sign}(\eta_2)$ topology and the $\textrm{sign}(\eta_1) = \textrm{sign}(\eta_2)$ topology. The $A_{LL}$ uncertainty contribution of $0.0007$ from uncertainty in the relative luminosity measurement and $6.1\%$ from the beam polarization uncertainty, which are common to all the data points, are separated from the listed values.