The average of the azimuthal-correlation functions of $\Dzero$, $\Dplus$, and $\Dstar$ mesons with associated particles, after the subtraction of the baseline, in \pp collisions at $\s = 5.02$~\TeV, compared to predictions from the PYTHIA, POWHEG+PYTHIA6, POWHEG LO+PYTHIA6, HERWIG, and EPOS 3 event generators with various configurations (see text for details). The functions are shown for $3 < \ptD < 5$~\GeVc, $5 < \ptD < 8$~\GeVc, $8 < \ptD < 16$~\GeVc, and $16 < \ptD < 24$~\GeVc (from left to right) and $\ptass > 0.3$~\GeVc, $0.3 < \ptass < 1$~\GeVc, $1 < \ptass < 2$~\GeVc, and $2 < \ptass < 3$~\GeVc (from top to bottom). Statistical and $\Dphi$-dependent systematic uncertainties are shown as vertical error bars and boxes, respectively, while the $\Dphi$-independent uncertainties are written as text. The uncertainties from the subtraction of the baseline are displayed as boxes at $\Dphi > \pi$.