$D^0\,R_{AA}$ for different centralities.

Centrality dependence of the $D^0$ $p_{\rm T}$ differential invariant yield in Au+Au collisions (solid symbols). The curves are number-of-binary-collision-scaled Levy functions from fitting to the p+p result (open circles), which has been updated with the latest global analysis of charm fragmentation ratios from Ref and also taking into account the $p_{\rm T}$ dependence of the fragmentation ratio between $D^0$ and $D^{*{\pm}}$ from PYTHIA 6.4. The arrow denotes the upper limit with 90% confidence level of the last data point for 10$-$40% collisions. The systematic uncertainties are shown as square brackets.

Centrality dependence of the $D^0$ $p_{\rm T}$ differential invariant yield in Au+Au collisions (solid symbols). The curves are number-of-binary-collision-scaled Levy functions from fitting to the p+p result (open circles), which has been updated with the latest global analysis of charm fragmentation ratios from Ref and also taking into account the $p_{\rm T}$ dependence of the fragmentation ratio between $D^0$ and $D^{*{\pm}}$ from PYTHIA 6.4. The arrow denotes the upper limit with 90% confidence level of the last data point for 10$-$40% collisions. The systematic uncertainties are shown as square brackets.

Integrated $D^0$ $R_{AA}$ as a function of $N_{part}$ in different $p_T$ regions.

Panels (ab), $D^0$ $R_{\rm AA}$ for peripheral 40$-$80% and semi a central 10$-$40% collisions; Panel (c), $D^0$ $R_{\rm AA}$ for 0$-$10% most central events (blue circles) compared with model calculations from the TAMU (solid curve), SUBATECH (dashed curve), Torino (dot-dashed curve), Duke (long-dashed and long-dot-dashed curves), and LANL groups (filled band). The open symbol indicates the result with the extrapolated p+p reference. The vertical lines and brackets around the data points denote the statistical and systematic uncertainties respectively. The vertical bars around unity denote the overall normalization uncertainties in the Au+Au and p+p data, respectively. The $R_{\rm AA}$ probability distribution for the 0$-$0.7 GeV/$c$ data point is largely skewed. The uncertainty we report is the 68.3% probability range with respect to the measured central value assuming Gaussian distribution.

Panels (ab), $D^0$ $R_{\rm AA}$ for peripheral 40$-$80% and semi a central 10$-$40% collisions; Panel (c), $D^0$ $R_{\rm AA}$ for 0$-$10% most central events (blue circles) compared with model calculations from the TAMU (solid curve), SUBATECH (dashed curve), Torino (dot-dashed curve), Duke (long-dashed and long-dot-dashed curves), and LANL groups (filled band). The open symbol indicates the result with the extrapolated p+p reference. The vertical lines and brackets around the data points denote the statistical and systematic uncertainties respectively. The vertical bars around unity denote the overall normalization uncertainties in the Au+Au and p+p data, respectively. The $R_{\rm AA}$ probability distribution for the 0$-$0.7 GeV/$c$ data point is largely skewed. The uncertainty we report is the 68.3% probability range with respect to the measured central value assuming Gaussian distribution.

Integrated $D^0$ $R_{\rm AA}$ as a function of $N_{\mathrm{part}}$ in different $p_{\rm T}$ regions, 0$-$8 ${\mathrm{GeV/}}c$ (squares), 0.7$-$2.2 ${\mathrm{GeV/}}c$ (diamonds) and 3$-$8 ${\mathrm{GeV/}}c$ (circles). Open symbols are for the 0$-$80% minimum bias events. The vertical bar around unity denotes the overall normalization uncertainty from p+p reference.

Integrated $D^0$ $R_{\rm AA}$ as a function of $N_{\mathrm{part}}$ in different $p_{\rm T}$ regions, 0$-$8 ${\mathrm{GeV/}}c$ (squares), 0.7$-$2.2 ${\mathrm{GeV/}}c$ (diamonds) and 3$-$8 ${\mathrm{GeV/}}c$ (circles). Open symbols are for the 0$-$80% minimum bias events. The vertical bar around unity denotes the overall normalization uncertainty from p+p reference.

Production rate for PI0 production in the rapidity range -9.4 to -9.6.

Production rate for PI0 production in the rapidity range -9.6 to -10.0.

Production rate for PI0 production in the rapidity range -10.0 to -11.0.

The differential cross section as a function of the transverse momentum of the top quark/antiquark, PT(TOP/TOPBAR).

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.

The K0S cross section for 5 jet ET ranges.

The LAMBDA cross section for 5 jet ET ranges.

The triple differential GAMMA+JET cross section for |y_gamma| < 1.0, |y_jet| 2.4 TO 3.2 and y_gamma*y_jet > 0 A common 6.8% nomalization is included in the (sys) error.

The triple differential GAMMA+JET cross section for |y_gamma| < 1.0, |y_jet| <= 0.8 and y_gamma*y_jet <= 0 A common 6.8% nomalization is included in the (sys) error.

The triple differential GAMMA+JET cross section for |y_gamma| < 1.0, |y_jet| 0.8 TO 1.6 and y_gamma*y_jet <= 0 A common 6.8% nomalization is included in the (sys) error.

The triple differential GAMMA+JET cross section for |y_gamma| < 1.0, |y_jet| 1.6 TO 2.4 and y_gamma*y_jet <= 0 A common 6.8% nomalization is included in the (sys) error.

The triple differential GAMMA+JET cross section for |y_gamma| < 1.0, |y_jet| 2.4 TO 3.2 and y_gamma*y_jet <= 0 A common 6.8% nomalization is included in the (sys) error.

The triple differential GAMMA+JET cross section for |y_gamma| 1.5 TO 2.5, |y_jet| <= 0.8 and y_gamma*y_jet > 0 A common 6.8% nomalization is included in the (sys) error.

The triple differential GAMMA+JET cross section for |y_gamma| 1.5 TO 2.5, |y_jet| 0.8 TO 1.6 and y_gamma*y_jet > 0 A common 6.8% nomalization is included in the (sys) error.

The triple differential GAMMA+JET cross section for |y_gamma| 1.5 TO 2.5, |y_jet| 1.6 TO 2.4 and y_gamma*y_jet > 0 A common 6.8% nomalization is included in the (sys) error.

The triple differential GAMMA+JET cross section for |y_gamma| 1.5 TO 2.5, |y_jet| 2.4 TO 3.2 and y_gamma*y_jet > 0 A common 6.8% nomalization is included in the (sys) error.

The triple differential GAMMA+JET cross section for |y_gamma| 1.5 TO 2.5, |y_jet| <= 0.8 and y_gamma*y_jet <= 0 A common 6.8% nomalization is included in the (sys) error.

The triple differential GAMMA+JET cross section for |y_gamma| 1.5 TO 2.5, |y_jet| 0.8 TO 1.6 and y_gamma*y_jet <= 0 A common 6.8% nomalization is included in the (sys) error.

The triple differential GAMMA+JET cross section for |y_gamma| 1.5 TO 2.5, |y_jet| 1.6 TO 2.4 and y_gamma*y_jet <= 0 A common 6.8% nomalization is included in the (sys) error.

The triple differential GAMMA+JET cross section for |y_gamma| 1.5 TO 2.5, |y_jet| 2.4 TO 3.2 and y_gamma*y_jet <= 0 A common 6.8% nomalization is included in the (sys) error.

Double Differential distributions for B/S0 production.

PT distributions for B0 production.

PT distributions for B+ production.

PT distributions for B/S0 production.

Rapidity distributions for B0 production.

Rapidity distributions for B+ production.

Rapidity distributions for B/S0 production.

The double-differential cross sections for prompt J/PSI production (assuming no polarisation) and production of J/PSI from b-hadron decays as a function of transverse momentum for the rapidity range 2.5-3.0. Also shown in the final column is the fraction (in %) of J/PSIs from the latter.

The double-differential cross sections for prompt J/PSI production (assuming no polarisation) and production of J/PSI from b-hadron decays as a function of transverse momentum for the rapidity range 3.0-3.5. Also shown in the final column is the fraction (in %) of J/PSIs from the latter.

The double-differential cross sections for prompt J/PSI production (assuming no polarisation) and production of J/PSI from b-hadron decays as a function of transverse momentum for the rapidity range 3.5-4.0. Also shown in the final column is the fraction (in %) of J/PSIs from the latter.

The double-differential cross sections for prompt J/PSI production (assuming no polarisation) and production of J/PSI from b-hadron decays as a function of transverse momentum for the rapidity range 4.0-4.5. Also shown in the final column is the fraction (in %) of J/PSIs from the latter.

The total integrated cross sections for UPSILON 1S 2S and 3S times their dimuon branching ratio in the rapidity range 2.0-4.5 and transverse momentum 0-15 GeV/c.

Differential production cross section for UPSILON 1S 2s and 3S mesons times their dimuon branching fraction as a function of transverse momentum ;.

Differential production cross section for UPSILON 1S 2s and 3S mesons times their dimuon branching fraction as a function of rapidity ;.

The double-differential cross section (times the dimuon decay probability) for the production of UPSILON 1S, 2S and 3S mesons as a function of their transverse momentum for the rapidity range 2.0-2.5, assuming no polarisation.

The double-differential cross section (times the dimuon decay probability) for the production of UPSILON 1S, 2S and 3S mesons as a function of their transverse momentum for the rapidity range 2.5-3.0, assuming no polarisation.

The double-differential cross section (times the dimuon decay probability) for the production of UPSILON 1S, 2S and 3S mesons as a function of their transverse momentum for the rapidity range 3.0-3.5, assuming no polarisation.

The double-differential cross section (times the dimuon decay probability) for the production of UPSILON 1S, 2S and 3S mesons as a function of their transverse momentum for the rapidity range 3.5-4.0, assuming no polarisation.

The double-differential cross section (times the dimuon decay probability) for the production of UPSILON 1S, 2S and 3S mesons as a function of their transverse momentum for the rapidity range 4.0-4.5, assuming no polarisation.

Ratios of the cross-sections for UPSILON(2S)->MU+MU- and UPSILON(3S)->MU+MU- with respect toe UPSILON(1S)->MU+MU- as a function of PT in the rapidity range 2.0-4.0 ;.

Ratios of the cross-sections for UPSILON(2S)->MU+MU- and UPSILON(3S)->MU+MU- with respect to UPSILON(1S)->MU+MU- as a function of rapidity in the rapidity range 2.0-4.0 ;.

D(SIG)/DQ**2 IN NB/GEV**2 as a function of Q**2 IN GEV**2.

D(SIG)/DY IN NB as a function of Y.

D(SIG)/DX IN NB as a function of X.

SIG(C=VIS) IN PB as a function of Y for the Q^2 range 5-9 GeV^2.

SIG(C=VIS) IN PB as a function of Y for the Q^2 range 9-14 GeV^2.

SIG(C=VIS) IN PB as a function of Y for the Q^2 range 14-23 GeV^2.

SIG(C=VIS) IN PB as a function of Y for the Q^2 range 23-45 GeV^2.

SIG(C=VIS) IN PB as a function of Y for the Q^2 range 45-100 GeV^2.

SIG(C=VIS) IN PB as a function of Y for the Q^2 range 100-158 GeV^2.

SIG(C=VIS) IN PB as a function of Y for the Q^2 range 158-251 GeV^2.

SIG(C=VIS) IN PB as a function of Y for the Q^2 range 251-1000 GeV^2.

Reduced cross section as a function of X for Q^2 = 7 GeV^2.

Reduced cross section as a function of X for Q^2 = 12 GeV^2.

Reduced cross section as a function of X for Q^2 = 18 GeV^2.

Reduced cross section as a function of X for Q^2 = 32 GeV^2.

Reduced cross section as a function of X for Q^2 = 60 GeV^2.

Reduced cross section as a function of X for Q^2 = 120 GeV^2.

Reduced cross section as a function of X for Q^2 = 200 GeV^2.

Reduced cross section as a function of X for Q^2 = 350 GeV^2.