Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 0.15 GeV, for centrality 30-50%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 0.15 GeV, for centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 5 GeV, for centrality 0-10%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 5 GeV, for centrality 10-30%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 5 GeV, for centrality 30-50%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 5 GeV, for centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 10 GeV, for centrality 0-10%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 10 GeV, for centrality 10-30%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 10 GeV, for centrality 30-50%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Charged jet spectra using two cone radius parameters R = 0.2 and 0.3 and a leading track selection of pT > 10 GeV, for centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Ratio of charged jet spectra with leading track selection of pT > 0.15 GeV to pT > 5 GeV using two cone radius parameters R = 0.2 and 0.3, for centrality 0-10%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Ratio of charged jet spectra with leading track selection of pT > 0.15 GeV to pT > 5 GeV using two cone radius parameters R = 0.2 and 0.3, for centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Ratio of charged jet spectra with leading track selection of pT > 10 GeV to pT > 5 GeV using two cone radius parameters R = 0.2 and 0.3, for centrality 0-10%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Ratio of charged jet spectra with leading track selection of pT > 10 GeV to pT > 5 GeV using two cone radius parameters R = 0.2 and 0.3, for centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 0.15 GeV. Central collisions are defined to have centrality 0-10% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 0.15 GeV. Central collisions are defined to have centrality 10-30% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 0.15 GeV. Central collisions are defined to have centrality 30-50% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 5 GeV. Central collisions are defined to have centrality 0-10% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 5 GeV. Central collisions are defined to have centrality 10-30% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 5 GeV. Central collisions are defined to have centrality 30-50% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 10 GeV. Central collisions are defined to have centrality 0-10% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 10 GeV. Central collisions are defined to have centrality 10-30% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, for charged jets with either R = 0.2 or R = 0.3 and a leading charged particle with pT > 10 GeV. Central collisions are defined to have centrality 30-50% and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, as a function of the average number of participants in the collision (<Npart>), for charged jets with pT = 60-70 GeV and either R = 0.2 or R = 0.3 and a leading charged particle with pT > 0.15 GeV. The correspondence between <Npart> and the centrality classes is given in Table 1 of the paper, and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, as a function of the average number of participants in the collision (<Npart>), for charged jets with pT = 60-70 GeV and either R = 0.2 or R = 0.3 and a leading charged particle with pT > 5 GeV. The correspondence between <Npart> and the centrality classes is given in Table 1 of the paper, and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Nuclear modification factor, constructed as the ratio of jet pT spectra in central and peripheral collisions normalized by the nuclear overlap functions, as a function of the average number of participants in the collision (<Npart>), for charged jets with pT = 60-70 GeV and either R = 0.2 or R = 0.3 and a leading charged particle with pT > 10 GeV. The correspondence between <Npart> and the centrality classes is given in Table 1 of the paper, and peripheral collisions are defined to have centrality 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Ratio of charged jet pT-spectra with radius parameter R = 0.2 and 0.3 and a leading charged particle with pT > 5 GeV, for two centrality classes 0-10% and 50-80%. The two systematic uncertainties correspond to the shape uncertainty and the correlated uncertainty.

Differential cross-section for $Z$+jet in the leading jet pseudorapidity, $\eta^{\mathrm{jet}}$.

Differential cross-section for $Z$+jet in the boson $Z$ rapidity, $y^Z$.

Cross-section for $Z$+jet production, differential in the $Z$ boson transverse momentum, $p_{T}^{Z}$.

Cross-section for $Z$+jet production, differential in the difference in $\phi$ between the $Z$ boson and the leading jet, $\Delta \phi$.

Cross-section for $Z$+jet production, differential in the difference in rapidity between the $Z$ boson and the leading jet, $\Delta y$.

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.

$\Delta N_p$ multiplicity distributions in Au+Au collisions at $\sqrt{S_{NN}}=27$ GeV for 0-5 percent, 30-40 percent and 70-80 percent collision centralities at midrapidity.

$\Delta N_p$ multiplicity distributions in Au+Au collisions at $\sqrt{S_{NN}}=39$ GeV for 0-5 percent, 30-40 percent and 70-80 percent collision centralities at midrapidity.

$\Delta N_p$ multiplicity distributions in Au+Au collisions at $\sqrt{S_{NN}}=62.4$ GeV for 0-5 percent, 30-40 percent and 70-80 percent collision centralities at midrapidity.

$\Delta N_p$ multiplicity distributions in Au+Au collisions at $\sqrt{S_{NN}}=200$ GeV for 0-5 percent, 30-40 percent and 70-80 percent collision centralities at midrapidity.

Centrality dependence of the cumulants of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{S_{NN}}=7.7$ GeV.

Centrality dependence of the cumulants of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{S_{NN}}=11.5$ GeV.

Centrality dependence of the cumulants of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{S_{NN}}=19.6$ GeV.

Centrality dependence of the cumulants of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{S_{NN}}=27$ GeV.

Centrality dependence of the cumulants of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{S_{NN}}=39$ GeV.

Centrality dependence of the cumulants of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{S_{NN}}=62.4$ GeV.

Centrality dependence of the cumulants of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{S_{NN}}=200$ GeV.

Centrality dependence of $S\sigma$/Skellam and $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{S_{NN}}=7.7$ GeV.

Centrality dependence of $S\sigma$/Skellam and $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{S_{NN}}=11.5$ GeV.

Centrality dependence of $S\sigma$/Skellam and $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{S_{NN}}=19.6$ GeV.

Centrality dependence of $S\sigma$/Skellam and $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{S_{NN}}=27$ GeV.

Centrality dependence of $S\sigma$/Skellam and $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{S_{NN}}=39$ GeV.

Centrality dependence of $S\sigma$/Skellam and $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{S_{NN}}=62.4$ GeV.

Centrality dependence of $S\sigma$/Skellam and $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{S_{NN}}=200$ GeV.

Collision energy and centrality dependence of the net-proton $S\sigma$ and $\kappa\sigma^2$ from Au+Au and p+p collisions at RHIC.

Collision energy and centrality dependence of the net-proton $S\sigma$ and $\kappa\sigma^2$ from Au+Au and p+p collisions at RHIC.

Collision energy and centrality dependence of the net-proton $S\sigma$ and $\kappa\sigma^2$ from Au+Au and p+p collisions at RHIC.

Collision energy and centrality dependence of the net-proton $S\sigma$ and $\kappa\sigma^2$ from Au+Au and p+p collisions at RHIC.

Cumulants of net-proton distribution at $\sqrt{S_{NN}}=7.7$ GeV. (efficiency corrected).

Cumulants of net-proton distribution at $\sqrt{S_{NN}}=11.5$ GeV. (efficiency corrected).

Cumulants of net-proton distribution at $\sqrt{S_{NN}}=19.6$ GeV. (efficiency corrected).

Cumulants of net-proton distribution at $\sqrt{S_{NN}}=27$ GeV. (efficiency corrected).

Cumulants of net-proton distribution at $\sqrt{S_{NN}}=39$ GeV. (efficiency corrected).

Cumulants of net-proton distribution at $\sqrt{S_{NN}}=62.4$ GeV. (efficiency corrected).

Cumulants of net-proton distribution at $\sqrt{S_{NN}}=200$ GeV. (efficiency corrected).

Cumulants of proton distribution at $\sqrt{S_{NN}}=7.7$ GeV. (efficiency corrected).

Cumulants of proton distribution at $\sqrt{S_{NN}}=11.5$ GeV. (efficiency corrected).

Cumulants of proton distribution at $\sqrt{S_{NN}}=19.6$ GeV. (efficiency corrected).

Cumulants of proton distribution at $\sqrt{S_{NN}}=27$ GeV. (efficiency corrected).

Cumulants of proton distribution at $\sqrt{S_{NN}}=39$ GeV. (efficiency corrected).

Cumulants of proton distribution at $\sqrt{S_{NN}}=62.4$ GeV. (efficiency corrected).

Cumulants of proton distribution at $\sqrt{S_{NN}}=200$ GeV. (efficiency corrected).

Cumulants of anti-proton distribution at $\sqrt{S_{NN}}=7.7$ GeV. (efficiency corrected).

Cumulants of anti-proton distribution at $\sqrt{S_{NN}}=11.5$ GeV. (efficiency corrected).

Cumulants of anti-proton distribution at $\sqrt{S_{NN}}=19.6$ GeV. (efficiency corrected).

Cumulants of anti-proton distribution at $\sqrt{S_{NN}}=27$ GeV. (efficiency corrected).

Cumulants of anti-proton distribution at $\sqrt{S_{NN}}=39$ GeV. (efficiency corrected).

Cumulants of anti-proton distribution at $\sqrt{S_{NN}}=62.4$ GeV. (efficiency corrected).

Cumulants of anti-proton distribution at $\sqrt{S_{NN}}=200$ GeV. (efficiency corrected).

$\pi^0$ meson $v_4, 20-40\%$ centrality

$\pi^0$ meson $v_4, 40-60\%$ centrality

$\pi^0$ meson $v_4, 20-60\%$ centrality

$\pi^0$ meson $v_4/v_2^2, 0-20\%$ centrality

$\pi^0$ meson $v_4/v_2^2, 20-40\%$ centrality

$\pi^0$ meson $v_4/v_2^2, 40-60\%$ centrality

$\pi^0$ meson $v_4/v_2^2, 20-60\%$ centrality

Comparison of data to Monte Carlo for the distributions of two-photon invariant mass (left) and energy for the higher (center) and lower (right) energy photon.

Invariant mass distribution for the two-photon system with 7 < pT < 8 GeV/c.

Signal fractions calculated within the “peak region” of 0.1 < $M_{\gamma\gamma}$ < 0.2 $GeV/c^{2}$. Fractions for the full dataset as well as the subsets of longitudinal and transverse polarizations are shown as a function of pT (left), and the fractions for the transversely polarized data are shown as a function of $x_{F}$ integrated over 5 < pT < 12 GeV/c (right).

Signal fractions calculated within the “peak region” of 0.1 < $M_{\gamma\gamma}$ < 0.2 $GeV/c^{2}$. Fractions for the full dataset as well as the subsets of longitudinal and transverse polarizations are shown as a function of pT (left), and the fractions for the transversely polarized data are shown as a function of $x_{F}$ integrated over 5 < pT < 12 GeV/c (right).

Upper panel, the $\pi^{0}$ cross section (blue markers) is shown compared with an NLO pQCD cal- culation with three options for the scale parameter. The lower panel presents the ratio of the data to the pT-scale theory curve, as well as the ratio of the 2pT-scale and pT/2-scale theory curves to the pT-scale curve.

Upper panel, the $\pi^{0}$ cross section (blue markers) is shown compared with an NLO pQCD cal- culation with three options for the scale parameter. The lower panel presents the ratio of the data to the pT-scale theory curve, as well as the ratio of the 2pT-scale and pT/2-scale theory curves to the pT-scale curve.

The ALL results (blue markers) are presented with the DSSV prediction and the GRSV prediction using the best fit to polarized DIS ($\Delta$ g = std) and the maximum and minimum allowed values for gluon polarization.

The ALL results (blue markers) are presented with the DSSV prediction and the GRSV prediction using the best fit to polarized DIS ($\Delta$ g = std) and the maximum and minimum allowed values for gluon polarization.

The ALL results (blue markers) are presented with the DSSV prediction and the GRSV prediction using the best fit to polarized DIS ($\Delta$ g = std) and the maximum and minimum allowed values for gluon polarization.

The AN results are plotted versus $x_{F}$ integrated over 5 < pT < 12 GeV/c (left panel) and versus pT integrated over 0.06 < |$x_{F}$| < 0.27 (right panel).

The AN results are plotted versus $x_{F}$ integrated over 5 < pT < 12 GeV/c (left panel) and versus pT integrated over 0.06 < |$x_{F}$| < 0.27 (right panel).

The AN results are plotted versus $x_{F}$ integrated over 5 < pT < 12 GeV/c (left panel) and versus pT integrated over 0.06 < |x_{F}| < 0.27 (right panel).

The AN results are plotted versus $x_{F}$ integrated over 5 < pT < 12 GeV/c (left panel) and versus pT integrated over 0.06 < |$x_{F}$| < 0.27 (right panel).

Single differential production cross sections of prompt J/PSI mesons and of J/PSI from B decay as a function of rapidity in the BACKWARD region. The errors shown are statistical and the uncorrelated and correlated components of the systematic uncertainties.

Double differential production cross sections of prompt J/PSI mesons as a function of transverse momentum and rapidity in the FORWARD region. The errors shown are statistical and the uncorrelated and correlated components of the systematic uncertainties.

Double differential production cross sections of J/PSI mesons from B decays as a function of transverse momentum and rapidity in the FORWARD region. The errors shown are statistical and the uncorrelated and correlated components of the systematic uncertainties.

The nuclear modification factor RpPB for PT < 14 GeV/c in the BACKWARD and FORWARD regions.

The FORWARD/BACKWARD production ratio R(FB) as a function of rapidity for PT < 14 GeV/c.

The FORWARD/BACKWARD production ratioo (R(FB) as a function of PT with |yrap|=2.5-4.0.

The K0S cross section for 5 jet ET ranges.

The LAMBDA cross section for 5 jet ET ranges.

Double-differential cross-section $d^2\sigma/dp_T/dy$ for prompt $J/\psi$ production in bins of $p_T$ and $y$, with statistical, systematic and polarization uncertainties.