Total $J/\psi$ efficiency as a function of cos$\theta$ for $2 < p_{T}^{J/\psi} < 3$ GeV/c

Total $J/\psi$ efficiency as a function of cos$\theta$ for $3 < p_{T}^{J/\psi} < 4$ GeV/c

Total $J/\psi$ efficiency as a function of cos$\theta$ for $4 < p_{T}^{J/\psi} < 6$ GeV/c

Different contributions to total $J/\psi$ efficiency as a function of cos$\theta$ for $2 < p_{T}^{J/\psi} < 3$ GeV/c

Different contributions to total $J/\psi$ efficiency as a function of cos$\theta$ for $3 < p_{T}^{J/\psi} < 4$ GeV/c

Different contributions to total $J/\psi$ efficiency as a function of cos$\theta$ for $4 < p_{T}^{J/\psi} < 6$ GeV/c

Corrected cos$\theta$ distribution for $2 < p_{T}^{J/\psi} < 3$ GeV/c

Corrected cos$\theta$ distribution for $3 < p_{T}^{J/\psi} < 4$ GeV/c

Corrected cos$\theta$ distribution for $4 < p_{T}^{J/\psi} < 6$ GeV/c

$J/\psi$ polarization parameter $\lambda_{\theta}$ as a function of $J/\psi$ $p_{T}$ for rapidity |y|<1

The invariant yield versus transverse momentum for |y| < 1 in 0-60% centrality in Au+Au collisions (solid circles). The results are compared to high-$p_T$ (3 < $p_T$ < 10 GeV/c) results from STAR [9] (solid squares) and PHENIX data [8] (open squares). Also shown is the yield in 0-60% centrality Cu+Cu collisions for low $p_T$ (solid triangles) and high $p_T$ [48] (open triangles). The models are described in the text [49, 50]. The $J/\psi$ cross section in p+p collisions is also shown at STAR (stars) and PHENIX [51] (diamonds), and the scale is indicated on the right axis.

The $J/\psi$ nuclear modification factor versus transverse momentum for |y| < 1 and $p_T$ < 5 GeV/c for 0-60% centrality Au+Au collisions (solid circles). The data are compared with STAR high-$p_T$ (5 < $p_T$ < 10 GeV/c) results [9] and PHENIX results [8] in |y| < 0.35 (open squares). The statistical and systematic uncertainties on the baseline $J/\psi$ cross section in p+p collisions are indicated by the hatched and solid bands, respectively. Boxes on the vertical axes represent the uncertainty on $N_{coll}$ combined with the uncertainty on the inelastic cross section in p+p collisions

The $J/\psi$ nuclear modification factor versus transverse momentum for |y| < 1 and $p_T$ < 5 GeV/c for 0-20% centrality Au+Au collisions (solid circles). The data are compared with STAR high-$p_T$ (5 < $p_T$ < 10 GeV/c) results [9] and PHENIX results [8] in |y| < 0.35 (open squares). The statistical and systematic uncertainties on the baseline $J/\psi$ cross section in p+p collisions are indicated by the hatched and solid bands, respectively. Boxes on the vertical axes represent the uncertainty on $N_{coll}$ combined with the uncertainty on the inelastic cross section in p+p collisions

The $J/\psi$ nuclear modification factor versus transverse momentum for |y| < 1 and $p_T$ < 5 GeV/c for 20-40% centrality Au+Au collisions (solid circles). The data are compared with STAR high-$p_T$ (5 < $p_T$ < 10 GeV/c) results [9] and PHENIX results [8] in |y| < 0.35 (open squares). The statistical and systematic uncertainties on the baseline $J/\psi$ cross section in p+p collisions are indicated by the hatched and solid bands, respectively. Boxes on the vertical axes represent the uncertainty on $N_{coll}$ combined with the uncertainty on the inelastic cross section in p+p collisions

The $J/\psi$ nuclear modification factor versus transverse momentum for |y| < 1 and $p_T$ < 5 GeV/c for 40-60% centrality Au+Au collisions (solid circles). The data are compared with STAR high-$p_T$ (5 < $p_T$ < 10 GeV/c) results [9] and PHENIX results [8] in |y| < 0.35 (open squares). The statistical and systematic uncertainties on the baseline $J/\psi$ cross section in p+p collisions are indicated by the hatched and solid bands, respectively. Boxes on the vertical axes represent the uncertainty on $N_{coll}$ combined with the uncertainty on the inelastic cross section in p+p collisions

The $J/\psi$ invariant yield and nuclear modification factor as a function of centrality for |y| < 1 in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $J/\psi$ invariant yield $Bd^2N/(2\pi p_Tdydp_T) [(GeV/c)^{-2}]$ and nuclear modification factor as a function of transverse momentum for |y| < 1 in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The $J/\psi$ invariant yield and nuclear modification factor for |y| < 1 in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The differential cross section of $Z$ boson production scaled by 1/$p_\mathrm{T}^{Z}$, (1/$p_\mathrm{T}^{Z}$) $d\sigma /dp_\mathrm{T}^{Z}$, for $-3<y^\mathrm{*}_{Z}<2$.

The differential cross section of $Z$ boson production scaled by 1/$p_\mathrm{T}^{Z}$, (1/$p_\mathrm{T}^{Z}$) $d\sigma /dp_\mathrm{T}^{Z}$, for $-2<y^\mathrm{*}_{Z}<0$ and $0<y^\mathrm{*}_{Z}<2.

Centrality bias corrected $Z$ boson yields per event for $-3<y^\mathrm{*}_{Z}<2$ scaled by by $\langle N_{coll}\rangle$. (To remove the centrality bias correction each value may be multiplied by the approriate correction value found in arXiv:1412.0976.).

The rapidity differential Z boson yields per event scaled by $\langle N_{coll}\rangle$ for three centrality ranges.

$\rho^0$ production cross section determined with the data set collected with trigger B, as a function of momentum transfer $t$, fitted with a double exponential fit function in the range $t < 0.1$. The statistical errors are shown. The ratio of incoherent to coherent cross sections is measured to be $0.20 \pm 0.8$. The fit parameters are shown in Table I.

Comparison of theoretical predictions to the measured total cross section for coherent $\rho^0$ production as a functionof $\sqrt{s_{NN}}$ . The measured cross section is based on the dataset collected with trigger B. The error bars show the sum of the statistical and systematic uncertainties. See text for details.

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.