The measured differential cross sections $\rm{d}\sigma/\rm{d}|\eta|$ of $D^{*\pm}$ and $D^\pm$ production with $20<p_T<100\,$GeV. The first and second errors are the statistical and systematic uncertainties, respectively. The systematic uncertainty fractions corresponding to the tracking ($\delta_2$) uncertainties (Table 2 of the paper) are strongly correlated. The fully correlated uncertainties linked with the luminosity measurement ($3.5\%$) and branching fractions ($1.5\%$ and $2.1\%$ for $D^{*\pm}$ and $D^\pm$, respectively) are not shown.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} p_T(D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} y(\Upsilon(1S))}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} y(\Upsilon(1S))}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} y(D^0)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} y(D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{\pi}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} \left| \Delta \phi \right| }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{\pi}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} \left| \Delta \phi \right| }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} \left| \Delta y \right| }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} \left| \Delta y \right| }$ for $2<y(\Upsilon(1S))<4.5$, for $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} p_T(\Upsilon(1S)D^0) }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} p_T(\Upsilon(1S)D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} y(\Upsilon(1S)D^0) }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} y(\Upsilon(1S)D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} \mathcal{A}_T }$, where $\mathcal{A}_T = \frac{p_T(\Upsilon(1S))-p_T(D^0)}{p_T(\Upsilon(1S))-p_T(D^0)}$, for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} \mathcal{A}_T }$, where $\mathcal{A}_T = \frac{p_T(\Upsilon(1S))-p_T(D^0)}{p_T(\Upsilon(1S))-p_T(D^0)}$, for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} m(\Upsilon(1S)D^0) }$, for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} m(\Upsilon(1S)D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

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