Top: invariant mass distributions of $\mathrm{D^{0}}$-jet candidates for $2< p_{\mathrm{T,ch\,jet}} <50$ GeV/$c$ and $R=0.4$ in pp collisions at $\sqrt{s}$ = 5.02 TeV, in the $\mathrm{D^{0}}$-meson transverse momentum intervals: $3< p_{\mathrm{T,D^0}} <4$ GeV/$c$ (left), $7< p_{\mathrm{T,D^0}} <8$ GeV/$c$ (centre), and $10< p_{\mathrm{T,D^0}} <12$ GeV/$c$ (right). The total fit function is represented by the blue solid line, while the red dashed line represents the sum of the background and reflection fit functions. 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 $p_{\mathrm{T,ch\,jet}}$ in the signal and sideband regions, and their subtracted yields.