Nuclear modification factors $Q_{\rm pPb}$ as a function of \pt in the 0-20\%, 20-40\%, 40-60\%, and 60-100\% multiplicity classes selected with the ZNA estimator in \pPb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV. The different panels of the figure are for different multiplicity classes. The vertical error bars and the empty boxes represent the statistical and systematic uncertainties, respectively. The solid boxes at high \pt at $Q_{\rm pPb}$ = 1 represent the normalisation uncertainties. The results are compared with the PHENIX results on electrons from heavy-flavour hadron decays  in d$+$Au collisions at $\sqrt{s_{\rm NN}}$ = 200 GeV and with ALICE charged particle results  in \pPb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV.