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$.