Top panels: \ptchjet-differential cross section of charm jets tagged with \Dzero mesons for $R=0.2$ (circles, scaled by 0.1), $0.4$ (squares) and $0.6$ (crosses, scaled by 10) in \pp collisions at \thirteen (left) and \five (right) compared to PYTHIA~8 HardQCD Monash~2013 (dash-dotted lines), PYTHIA~8 SoftQCD Mode 2 (dashed lines), and POWHEG hvq + PYTHIA~8 (open circles) predictions. The shaded bands indicate the systematic uncertainty on the data cross section while open boxes represent the theoretical uncertainties on the POWHEG predictions. Bottom panels: ratios of MC predictions to the data for $R=0.2$, $0.4$ and $0.6$.