Top: invariant mass distribution of $\mathrm{D^{0}}$-jet candidates for one jet-$p_{\rm T}$ interval of ${5< p_{\mathrm{T,ch\,jet}} <7}$ GeV/$c$ and $R=0.4$ in pp collisions at $\sqrt{s}$ = 5.02 TeV, in $\mathrm{D^{0}}$-meson transverse momentum intervals: $3< p_{\mathrm{T,D^0}} <4$ GeV/$c$ (left), $4< p_{\mathrm{T,D^0}} <5$ GeV/$c$ (centre), and $5< p_{\mathrm{T,D^0}} <6$ GeV/$c$ (right). The total fit function is represented by the blue solid line, while the red dashed line represents the background fit function. The red and green shaded areas correspond to the peak and sideband regions, respectively. Bottom: $\mathrm{D^{0}}$-jet raw yields as a function of $z_{||}^{\textrm{ch}}$ in the signal and sideband regions, and their subtracted yields.