Inclusive J/$\psi$ nuclear modification factor $R_{\rm AA}$ at midrapidity, integrated over $p_{\rm T}$, as a function of mean number of participating nucleons in Pb$-$Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV. The correlated systematic uncertainty due to the pp reference is not included and amounts to 6.9%, shown as the red box around unity in the paper figure.

Inclusive J/$\psi$ nuclear modification factor $R_{\rm AA}$ at midrapidity for $p_{\rm T} > 0.15$ GeV as a function of mean number of participating nucleons in Pb$-$Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV. The low-$p_{\rm T}$ interval is removed to reject J/$\psi$ from photo-production. The correlated systematic uncertainty due to the pp reference cross section is not listed and amounts to 6.9%, shown as the red box around unity in the paper figure.

Inclusive J/$\psi$ nuclear modification factor, $R_{\rm AA}$, at midrapidity in Pb$-$Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV as a function of transverse momentum, $p_{\rm T}$, for different centrality intervals. The systematic uncertianty from the pp reference, which is correlated with centrality and shown as grey band in the paper figure, is listed as sys,Reference. The systematic uncertainties on the $T_{\rm AA}$ which are correlated in $p_{\rm T}$ for each centrality bin are not listed and amount additionally to 2\%, 3\%, and 5\% for 0-20\%, 20-40\%, and 40-90\% centrality, respectively. These uncertainties are shown as boxes around unity in the paper figure.

Nuclear modification factor $R_{\rm AA}$ of inclusive J/$\psi$ at mid-rapidity. For the midrapidity point all normalization uncertainties are included.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) with $16 < p_{\rm T} < 24$ GeV/$c$ and charged particles with $p_{\rm T} > 0.3$ GeV/$c$, in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, after baseline subtraction. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) with $3 < p_{\rm T} < 5$ GeV/$c$ and charged particles with $0.3 < p_{\rm T} < 1.0$ GeV/$c$, in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, after baseline subtraction. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) with $5 < p_{\rm T} < 8$ GeV/$c$ and charged particles with $0.3 < p_{\rm T} < 1.0$, in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, after baseline subtraction. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $0.3 < p_{\rm T} < 1.0$ GeV/$c$, in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, after baseline subtraction. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) with $16 < p_{\rm T} < 24$ GeV/$c$ and charged particles with $0.3 < p_{\rm T} < 1.0$ GeV/$c$, in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, after baseline subtraction. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) with $3 < p_{\rm T} < 5$ GeV/$c$ and charged particles with $1.0 < p_{\rm T} < 2.0$ GeV/$c$, in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, after baseline subtraction. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) with $5 < p_{\rm T} < 8$ GeV/$c$ and charged particles with $1.0 < p_{\rm T} < 2.0$ GeV/$c$, in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, after baseline subtraction. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $1.0 < p_{\rm T} < 2.0$ GeV/$c$, in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, after baseline subtraction. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) with $16 < p_{\rm T} < 24$ GeV/$c$ and charged particles with $1.0 < p_{\rm T} < 2.0$ GeV/$c$, in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, after baseline subtraction. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) with $3 < p_{\rm T} < 5$ GeV/$c$ and charged particles with $2.0 < p_{\rm T} < 3.0$ GeV/$c$, in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, after baseline subtraction. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) with $5 < p_{\rm T} < 8$ GeV/$c$ and charged particles with $2.0 < p_{\rm T} < 3.0$ GeV/$c$, in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, after baseline subtraction. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $2.0 < p_{\rm T} < 3.0$ GeV/$c$, in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, after baseline subtraction. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) with $16 < p_{\rm T} < 24$ GeV/$c$ and charged particles with $2.0 < p_{\rm T} < 3.0$ GeV/$c$, in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, after baseline subtraction. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the near-side associated peak yields measured in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $p_{\rm T} > 0.3$ GeV/$c$. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak yields measured in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1.0$ GeV/$c$. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak yields measured in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1.0 < p_{\rm T} < 2.0$. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak yields measured in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2.0 < p_{\rm T} < 3.0$. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak width measured in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $p_{\rm T} > 0.3$. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak width measured in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1.0$. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak width measured in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1.0 < p_{\rm T} < 2.0$. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak width measured in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2.0 < p_{\rm T} < 3.0$. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak yields measured in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $p_{\rm T} > 0.3$ GeV/$c$. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak yields measured in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1.0$ GeV/$c$. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak yields measured in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1.0 < p_{\rm T} < 2.0$. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak yields measured in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2.0 < p_{\rm T} < 3.0$. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak width measured in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $p_{\rm T} > 0.3$. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak width measured in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1.0$. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak width measured in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1.0 < p_{\rm T} < 2.0$. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak width measured in pp collisions at $\sqrt{s} = 5.02$ TeV and p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2.0 < p_{\rm T} < 3.0$. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp, $-0.96 < y^{\rm D}_{\rm cms} < 0.04$ in p-Pb. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Azimuthal-correlation distributions of D mesons and charged particles obtained for D$^0$, D$^{+}$ and D$^{*+}$ mesons for $3 < p_{\rm T} < 5$ GeV/$c$, charged particles $p_{\rm T} > 0.3$ GeV/$c$, in p-Pb collisions at $\sqrt{s} = 5.02$ TeV in 0-20%, 20-60% and 60-100% centrality classes. Rapidity range for the D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi< \pi$.

Azimuthal-correlation distributions of D mesons and charged particles obtained for D$^0$, D$^{+}$ and D$^{*+}$ mesons for $5 < p_{\rm T} < 8$ GeV/$c$, charged particles $p_{\rm T} > 0.3$ GeV/$c$, in p-Pb collisions at $\sqrt{s} = 5.02$ TeV in 0-20%, 20-60% and 60-100% centrality classes. Rapidity range for the D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi< \pi$.

Azimuthal-correlation distributions of D mesons and charged particles obtained for D$^0$, D$^{+}$ and D$^{*+}$ mesons for $8 < p_{\rm T} < 16$ GeV/$c$, charged particles $p_{\rm T} > 0.3$ GeV/$c$, in p-Pb collisions at $\sqrt{s} = 5.02$ TeV in 0-20%, 20-60% and 60-100% centrality classes. Rapidity range for the D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi< \pi$.

Azimuthal-correlation distributions of D mesons and charged particles obtained for D$^0$, D$^{+}$ and D$^{*+}$ mesons for $3 < p_{\rm T} < 5$ GeV/$c$, charged particles $0.3 < p_{\rm T}< 1.0$ GeV/$c$, in p-Pb collisions at $\sqrt{s} = 5.02$ TeV in 0-20%, 20-60% and 60-100% centrality classes. Rapidity range for the D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi< \pi$.

Azimuthal-correlation distributions of D mesons and charged particles obtained for D$^0$, D$^{+}$ and D$^{*+}$ mesons for $5 < p_{\rm T} < 8$ GeV/$c$, charged particles $0.3 < p_{\rm T}< 1.0$ GeV/$c$, in p-Pb collisions at $\sqrt{s} = 5.02$ TeV in 0-20%, 20-60% and 60-100% centrality classes. Rapidity range for the D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi< \pi$.

Azimuthal-correlation distributions of D mesons and charged particles obtained for D$^0$, D$^{+}$ and D$^{*+}$ mesons for $8 < p_{\rm T} < 16$ GeV/$c$, charged particles $0.3 < p_{\rm T}< 1.0$ GeV/$c$, in p-Pb collisions at $\sqrt{s} = 5.02$ TeV in 0-20%, 20-60% and 60-100% centrality classes. Rapidity range for the D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi< \pi$.

Azimuthal-correlation distributions of D mesons and charged particles obtained for D$^0$, D$^{+}$ and D$^{*+}$ mesons for $3 < p_{\rm T} < 5$ GeV/$c$, charged particles $p_{\rm T}> 1.0$ GeV/$c$, in p-Pb collisions at $\sqrt{s} = 5.02$ TeV in 0-20% and 60-100% centrality classes. Rapidity range for the D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi< \pi$.

Azimuthal-correlation distributions of D mesons and charged particles obtained for D$^0$, D$^{+}$ and D$^{*+}$ mesons for $5 < p_{\rm T} < 8$ GeV/$c$, charged particles $p_{\rm T}> 1.0$ GeV/$c$, in p-Pb collisions at $\sqrt{s} = 5.02$ TeV in 0-20%, 20-60% and 60-100% centrality classes. Rapidity range for the D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi< \pi$.

Azimuthal-correlation distributions of D mesons and charged particles obtained for D$^0$, D$^{+}$ and D$^{*+}$ mesons for $8 < p_{\rm T} < 16$ GeV/$c$, charged particles $p_{\rm T}> 1.0$ GeV/$c$, in p-Pb collisions at $\sqrt{s} = 5.02$ TeV in 0-20%, 20-60% and 60-100% centrality classes. Rapidity range for the D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi< \pi$.

Comparison of the near-side associated peak yields measured in p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV in 0-20%, 20-60% and 60-100% centrality classes, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $p_{\rm T} > 0.3$ GeV/$c$. The rapidity range of D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak yields measured in p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV in 0-20%, 20-60% and 60-100% centrality classes, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1.0$ GeV/$c$. The rapidity range of D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak yields measured in p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV in 0-20%, 20-60% and 60-100% centrality classes, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $p_{\rm T} > 1.0$ GeV/$c$. The rapidity range of D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side peak widths measured in p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV in 0-20%, 20-60% and 60-100% centrality classes, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $p_{\rm T} > 0.3$ GeV/$c$. The rapidity range of D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side peak widths measured in p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV in 0-20%, 20-60% and 60-100% centrality classes, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1.0$ GeV/$c$. The rapidity range of D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side peak widths measured in p-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV in 0-20%, 20-60% and 60-100% centrality classes, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $p_{\rm T} > 1.0$ GeV/$c$. The rapidity range of D mesons is $-0.96 < y^{\rm D}_{\rm cms} < 0.04$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Baseline height measured in pp collisions at $\sqrt{s} = 5.02$ TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $p_{\rm T} > 0.3$ GeV/$c$. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Baseline height measured in pp collisions at $\sqrt{s} = 5.02$ TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1.0$ GeV/$c$. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Baseline height measured in pp collisions at $\sqrt{s} = 5.02$ TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1.0 < p_{\rm T} < 2.0$ GeV/$c$. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Baseline height measured in pp collisions at $\sqrt{s} = 5.02$ TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2.0 < p_{\rm T} < 3.0$ GeV/$c$. Rapidity range for the D mesons are $|y^{\rm D}_{\rm cms}| < 0.5$ in pp. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Fig. A1: Number density $N_{ch}$ (left) and $\\Sigma p_{T}$ (right) distributions as a function of pleadingT and the comparisons to MC predictions in Toward (top), Transverse (middle), and Away (bottom) regions for $p_{T}^{track} >$ 0.5 GeV/$c$. The shaded areas in the upper panels represent the systematic uncertainties and vertical error bars indicate statistical uncertainties. In the lower panels, the shaded areas are the sum in quadrature of statistical and systematic uncertainties from the upper panels. No uncertainties are given for the MC calculations.

Fig. A2: Number density $N_{ch}$ (left) and $\\Sigma p_{T}$ (right) distributions as a function of pleadingT and the comparisons to MC predictions in Toward (top), Transverse (middle), and Away (bottom) regions for $p_{T}^{track} >$ 1.0 GeV/$c$. The shaded areas in the upper panels represent the systematic uncertainties and vertical error bars indicate statistical uncertainties. In the lower panels, the shaded areas are the sum in quadrature of statistical and systematic uncertainties from the upper panels. No uncertainties are given for the MC calculations.

Measured second cumulants of anti-protons.

Measured first cumulants of protons.

Measured first cumulants of anti-protons.

Measured R1 values.

Measured R2 values.

R2 values from Model.

R1 values vs. delaY.

Global baryon number conservation.

Local baryon number conservation, deltaY_corr = 5.

Local baryon number conservation, deltaY_corr = 2.

HIJING results.

$^3\overline{\mathrm{H}} /\,^3\mathrm{H}$ ratio in $\mathrm{INEL}>0$ p--Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Transverse momentum spectra obtained from the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for $\mathrm{INEL}>0$ events in p$-$Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Transverse momentum spectra obtained from the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for events in 0$-$10% multiplicity class in p$-$Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Transverse momentum spectra obtained from the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for events in 10$-$20% multiplicity class in p$-$Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Transverse momentum spectra obtained from the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for events in 20$-$40% multiplicity class in p--Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Transverse momentum spectra obtained from the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for events in 40$-$100% multiplicity class in p--Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Mean transverse momentum obtained from the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ as a function of the mean charged-particle multiplicity density in p--Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$ TeV

$(^3\mathrm{He} +\,^3\overline{\mathrm{He}}) / (\mathrm{p} + \overline{\mathrm{p}})$ ratio in p$-$Pb as a function of the mean charged-particle multiplicity density

$(^3\mathrm{H} +\,^3\overline{\mathrm{H}}) / (\mathrm{p} + \overline{\mathrm{p}})$ ratio in p$-$Pb as a function of the mean charged-particle multiplicity density

Coalescence parameter $B_3$ calculated using the average of $\mathrm{INEL}>0$ $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ yields

Coalescence parameter $B_3$ calculated using the average of $\mathrm{INEL}>0$ $^3\mathrm{H}$ and $^3\overline{\mathrm{H}}$ yields

Ratio of (anti-)$^3$H and (anti-)$^3$He yields

Coalescence parameter $B_3$ calculated using the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for $\mathrm{INEL}>0$ events in p$-$Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Coalescence parameter $B_3$ calculated using the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for events in 0$-$10% multiplicity class in p$-$Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Coalescence parameter $B_3$ calculated using the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for events in 10$-$20% multiplicity class in p$-$Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Coalescence parameter $B_3$ calculated using the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for events in 20$-$40% multiplicity class in p--Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Coalescence parameter $B_3$ calculated using the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for events in 40$-$100% multiplicity class in p--Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Integrated production yield obtained from the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ corrected for the spin degeneracy factor $(2J +1)$ for minimum-bias p--Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Upper limit on the integrated production yield of $^4\overline{\mathrm{He}}$ at 90$\%$ confidence level corrected for the spin degeneracy factor $(2J +1)$ for minimum-bias p--Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Integrated production yield obtained from the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for the different multiplicity classes and $\mathrm{INEL} > 0$ p--Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Integrated production yield obtained from the average of $^3\mathrm{H}$ and $^3\overline{\mathrm{H}}$ for $\mathrm{INEL} > 0$ p--Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

$cos\theta^{*}$ distribution of $\phi$ meson w.r.t. Production Plane in transverse momentum range 0.5-0.8 GeV/$c$ measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$cos\theta^{*}$ distribution of $\rm{K}^{*0}$ (average of particle and anti-particle) meson w.r.t. Random Plane in transverse momentum range 0.8-1.2 GeV/$c$ measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$cos\theta^{*}$ distribution of $\phi$ meson w.r.t. Random Plane in transverse momentum range 0.5-0.7 GeV/$c$ measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$cos\theta^{*}$ distribution of K$^{0}_{S}$ meson w.r.t. Event Plane in transverse momentum range 0.6-0.8 GeV/$c$ measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$cos\theta^{*}$ distribution of K$^{0}_{S}$ meson w.r.t. Production Plane in transverse momentum range 0.6-0.8 GeV/$c$ measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$cos\theta^{*}$ distribution of $\rm{K}^{*0}$ (average of particle and anti-particle) meson w.r.t. Production Plane in transverse momentum range 0-0.6 GeV/$c$ measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$cos\theta^{*}$ distribution of $\phi$ meson w.r.t. Production Plane in transverse momentum range 0.5-0.8 GeV/$c$ measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$\rho_{00}$ as a function of $p_{\rm T}$ w.r.t. Event Plane for $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $p_{\rm T}$ w.r.t. Event Plane for $\rm{K}^{0}_{S}$ meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $p_{\rm T}$ w.r.t. Event Plane for $\phi$ meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $p_{\rm T}$ w.r.t. Production Plane for $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $p_{\rm T}$ w.r.t. Production Plane for $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$\rho_{00}$ as a function of $p_{\rm T}$ w.r.t. Production Plane for $\rm{K}^{0}_{S}$ meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $p_{\rm T}$ w.r.t. Production Plane for $\phi$ meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $p_{\rm T}$ w.r.t. Production Plane for $\phi$ meson measured in pp collisions at $\sqrt{s}$ = 13 TeV.

$\rho_{00}$ as a function of $p_{\rm T}$ w.r.t. Random Plane for $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $p_{\rm T}$ w.r.t. Random Plane for $\phi$ meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $\langle N_{\rm{part}} \rangle$ w.r.t. Production Plane in transverse momentum range 0.4-1.2 GeV/$c$ for $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $\langle N_{\rm{part}} \rangle$ w.r.t. Production Plane in transverse momentum range 3.0-5.0 GeV/$c$ for $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $\langle N_{\rm{part}} \rangle$ w.r.t. Production Plane in transverse momentum range 0.5-0.8 GeV/$c$ for $\phi$ meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $\langle N_{\rm{part}} \rangle$ w.r.t. Production Plane in transverse momentum range 3.0-5.0 GeV/$c$ for $\phi$ meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $\langle N_{\rm{part}} \rangle$ w.r.t. Event Plane in transverse momentum range 0.8-1.2 GeV/$c$ for $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $\langle N_{\rm{part}} \rangle$ w.r.t. Event Plane in transverse momentum range 3.0-5.0 GeV/$c$ for $\rm{K}^{*0}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $\langle N_{\rm{part}} \rangle$ w.r.t. Event Plane in transverse momentum range 0.5-0.7 GeV/$c$ for $\phi$ meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\rho_{00}$ as a function of $\langle N_{\rm{part}} \rangle$ w.r.t. Event Plane in transverse momentum range 3.0-5.0 GeV/$c$ for $\phi$ meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV.

$\Upsilon$(1S) differential cross section as a function of $p_{\rm T}$, at forward rapidity at $\sqrt{s_{\rm NN}}$ = 8.16 TeV .The first uncertainty is statistical, while the second is the systematic.

$\Upsilon$(1S) differential cross section as a function of $p_{\rm T}$, at backward rapidity at $\sqrt{s_{\rm NN}}$ = 8.16 TeV .The first uncertainty is statistical, while the second is the systematic.

Ratio of $\Upsilon$(2S) over $\Upsilon$(1S) yields in p--Pb collisions at $\sqrt{s_{\rm NN}}$ = 8.16 TeV .The first uncertainty is statistical, while the second is the systematic.

Ratio of $\Upsilon$(3S) over $\Upsilon$(1S) yields in p--Pb collisions at $\sqrt{s_{\rm NN}}$ = 8.16 TeV .The first uncertainty is statistical, while the second is the systematic.

$\Upsilon$(1S) $R_{\rm pPb}$ values at $\sqrt{s_{\rm NN}}$ = 8.16 TeV. The first uncertainty is statistical, the second is the uncorrelated systematic, while the third one is a correlated systematic uncertainty.

$\Upsilon$(1S) $R_{\rm pPb}$ as a function of $y_{\rm {cms}}$ at $\sqrt{s_{\rm NN}}$ = 8.16 TeV. The first uncertainty is statistical, the second is the uncorrelated systematic, while the third one is a correlated systematic uncertainty.

$\Upsilon$(1S) $R_{\rm pPb}$ as a function of $p_{\rm T}$ at backward rapidity at $\sqrt{s_{\rm NN}}$ = 8.16 TeV. The first uncertainty is statistical, the second is the uncorrelated systematic, while the third one is a correlated systematic uncertainty.

$\Upsilon$(1S) $R_{\rm pPb}$ as a function of $p_{\rm T}$ at forward rapidity at $\sqrt{s_{\rm NN}}$ = 8.16 TeV. The first uncertainty is statistical, the second is the uncorrelated systematic, while the third one is a correlated systematic uncertainty.

$\Upsilon$(1S) $Q_{\rm pPb}$ as a function of $<N_{\rm coll}>$ at backward rapidity at $\sqrt{s_{\rm NN}}$ = 8.16 TeV. The first uncertainty is statistical, the second is the uncorrelated systematic, while the third one is a correlated systematic uncertainty.

$\Upsilon$(1S) $Q_{\rm pPb}$ as a function of $<N_{\rm coll}>$ at forward rapidity at $\sqrt{s_{\rm NN}}$ = 8.16 TeV. The first uncertainty is statistical, the second is the uncorrelated systematic, while the third one is a correlated systematic uncertainty.

$\Upsilon$(1S) $R_{\rm pPb}$ values at $\sqrt{s_{\rm NN}}$ = 8.16 TeV. The first uncertainty is statistical, the second is the uncorrelated systematic, while the third one is a correlated systematic uncertainty.

$\Upsilon$(2S) $R_{\rm pPb}$ values at $\sqrt{s_{\rm NN}}$ = 8.16 TeV. The first uncertainty is statistical, the second is the uncorrelated systematic, while the third one is a correlated systematic uncertainty.

$\Upsilon$(3S) $R_{\rm pPb}$ values at $\sqrt{s_{\rm NN}}$ = 8.16 TeV. The first uncertainty is statistical, the second is the uncorrelated systematic, while the third one is a correlated systematic uncertainty.

Ratio of $\Upsilon$(2S) to $\Upsilon$(1S) $R_{\rm pPb}$ at $\sqrt{s_{\rm NN}}$ = 8.16 TeV as a function of $y_{\rm cms}$ .The first uncertainty is statistical, while the second is the systematic.

Ratio of $\Upsilon$(3S) to $\Upsilon$(1S) $R_{\rm pPb}$ at $\sqrt{s_{\rm NN}}$ = 8.16 TeV as a function of $y_{\rm cms}$ .The first uncertainty is statistical, while the second is the systematic.

Directed flow D0 meson 10-40%

Directed flow D0-bar meson 10-40%

Delta Directed flow D0 meson 10-40%

Slope Dv1 charged hadrons 540%

Slope Dv1 D0 meson 540%

K$^{*0}$ transverse momentum spectrum - V0M multiplicity class IV+V, average of particle and antiparticle

K$^{*0}$ transverse momentum spectrum - V0M multiplicity class VI, average of particle and antiparticle

K$^{*0}$ transverse momentum spectrum - V0M multiplicity class VII, average of particle and antiparticle

K$^{*0}$ transverse momentum spectrum - V0M multiplicity class VIII, average of particle and antiparticle

K$^{*0}$ transverse momentum spectrum - V0M multiplicity class IX, average of particle and antiparticle

K$^{*0}$ transverse momentum spectrum - V0M multiplicity class X, average of particle and antiparticle

K$^{*0}$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class I

K$^{*0}$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class II

K$^{*0}$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class III

K$^{*0}$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class IV+V

K$^{*0}$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class VI

K$^{*0}$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class VII

K$^{*0}$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class VIII

K$^{*0}$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class IX

K$^{*0}$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class X

$\phi$ transverse momentum spectrum - V0M multiplicity class I

$\phi$ transverse momentum spectrum - V0M multiplicity class II

$\phi$ transverse momentum spectrum - V0M multiplicity class III

$\phi$ transverse momentum spectrum - V0M multiplicity class IV

$\phi$ transverse momentum spectrum - V0M multiplicity class V

$\phi$ transverse momentum spectrum - V0M multiplicity class VI

$\phi$ transverse momentum spectrum - V0M multiplicity class VII

$\phi$ transverse momentum spectrum - V0M multiplicity class VIII

$\phi$ transverse momentum spectrum - V0M multiplicity class IX

$\phi$ transverse momentum spectrum - V0M multiplicity class X

$\phi$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class I

$\phi$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class II

$\phi$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class III

$\phi$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class IV

$\phi$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class V

$\phi$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class VI

$\phi$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class VII

$\phi$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class VIII

$\phi$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class IX

$\phi$ transverse momentum spectrum ratio to INEL>0 - V0M multiplicity class X

K$^{*0}$ average transverse momentum vs multiplicity - V0M multiplicity classes

$\phi$ average transverse momentum vs multiplicity - V0M multiplicity classes

K$^{*0}$ transverse momentum integrated yield vs multiplicity - V0M multiplicity classes, average of particle and antiparticle

$\phi$ transverse momentum integrated yield vs multiplicity - V0M multiplicity classes

K$^{*0}$/K$_{\mathrm{S}}^{0}$ transverse momentum integrated yield ratio vs multiplicity - V0M multiplicity classes, average of particle and antiparticle

$\phi$/K$_{\mathrm{S}}^{0}$ transverse momentum integrated yield ratio vs multiplicity - V0M multiplicity classes

$(\Xi^{-} + \bar{\Xi}^{+})$/$\phi$ transverse momentum integrated yield ratio vs multiplicity - V0M multiplicity classes

$\phi$/$(\pi^{+} + \pi^{-})$ transverse momentum integrated yield ratio vs multiplicity - V0M multiplicity classes

K$^{*0}$/K$_{\mathrm{S}}^{0}$ yield ratio vs transverse momentum - V0M multiplicity class II, average of particle and antiparticle

K$^{*0}$/K$_{\mathrm{S}}^{0}$ yield ratio vs transverse momentum - V0M multiplicity class X, average of particle and antiparticle

$\phi$/K$_{\mathrm{S}}^{0}$ yield ratio vs transverse momentum - V0M multiplicity class II

$\phi$/K$_{\mathrm{S}}^{0}$ yield ratio vs transverse momentum - V0M multiplicity class X

K$^{*0}$/K$_{\mathrm{S}}^{0}$ double yield ratio vs transverse momentum - V0M multiplicity class II / V0M multiplicity class X, average of particle and antiparticle

$\phi$/K$_{\mathrm{S}}^{0}$ double yield ratio vs transverse momentum - V0M multiplicity class II / V0M multiplicity class X

Significance of K$^{*0}$/K$_{\mathrm{S}}^{0}$ double yield ratio vs transverse momentum - V0M multiplicity class II / V0M multiplicity class X

Significance of $\phi$/K$_{\mathrm{S}}^{0}$ double yield ratio vs transverse momentum - V0M multiplicity class II / V0M multiplicity class X

($\rm{K}^{*0}$+$\bar{\rm{K}}^{*0}$)/(K$^{+}$+K$^{-}$) transverse momentum integrated yield ratio vs multiplicity - V0M multiplicity classes

($\rm{K}^{*0}$+$\bar{\rm{K}}^{*0}$)/($\pi^{+}$+$\pi^{-}$) transverse momentum integrated yield ratio vs multiplicity - V0M multiplicity classes

2$\phi$/(K$^{+}$+K$^{-}$) transverse momentum integrated yield ratio vs multiplicity - V0M multiplicity classes

$\phi$/($\rm{P}$ + $\bar{\rm{P}}$) transverse momentum integrated yield ratio vs multiplicity - V0M multiplicity classes

$\phi$/$(\Lambda~+~\bar{\Lambda})$ transverse momentum integrated yield ratio vs multiplicity - V0M multiplicity classes

$(\Omega^{-}~+~\bar{\Omega}^{+})$/$\phi$ transverse momentum integrated yield ratio vs multiplicity - V0M multiplicity classes