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

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The $\pi^0$ invariant mass peak positions and widths as a function of $p_{T}$ in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Photon conversion point radius distribution from real data and MC simulation in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Conversion probability correction factor for $\pi^0$ as a function of $p_{T}$

The overall detection efficiency of $\pi^0$ from embedding studt in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The overall detection efficiency of $\pi^0$ from embedding studt in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The overall detection efficiency of $\pi^0$ from embedding studt in 0-80% Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Invariant yield of STAR $\pi^0$ as a function of $p_{T}$ at midrapidity for different collision centrality bins in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Invariant yield of STAR $\pi^0$ as a function of $p_{T}$ at midrapidity for different collision centrality bins in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

Invariant yield of STAR $\pi^0$ as a function of $p_{T}$ at midrapidity for different collision centrality bins in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The ratios of STAR $\pi^0$ spectra over $\pi^{\pm}$ from STAR and $\pi^0$ from PHENIX for different collision centrality bins in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The nuclear modification factor $R_{CP}$ as a functon of $p_{T}$ of STAR $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.

The nuclear modification factor $R_{AA}$ as a functon of $p_{T}$ of STAR $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV.