Top panels: zch || -differential yield of R = 0.3 charm jets tagged with D0 mesons normalised by the number of D0 jets within each distribution in pp collisions at √s = 5.02 TeV in four pT,ch jet intervals (top left) 5 < pT,ch jet < 7 GeV/c, (top right) 7 < pT,ch jet < 10 GeV/c, (bottom left) 10 < pT,ch jet < 15 GeV/c and (bottom right) 15 < pT,ch jet < 50 GeV/c GeV/c. They are compared to PYTHIA 8 Monash 2013 (dashed-dotted lines), PYTHIA 8 Monash 2013 SoftQCD Mode 2 (dashed lines) and POWHEG hvq + PYTHIA 8 (open circles) predictions. The shaded bands indicate the systematic uncertainty on the distributions while open boxes represent the theoretical uncertainties on the POWHEG predictions. Bottom panels present ratios of MC predictions to the data