Ratio of $p_{\mathrm{T,ch\ jet}}$-differential cross section of charm jets tagged with $\mathrm{D^{0}}$ mesons for different $R$: $\sigma(R=0.2)/\sigma(R=0.4)$ and $\sigma(R=0.2)/\sigma(R=0.6)$ in pp collisions at $\sqrt{s}=13$ TeV.

Ratio of $p_{\mathrm{T,ch\ jet}}$-differential cross section of charm jets tagged with $\mathrm{D^{0}}$ mesons for different $R$: $\sigma(R=0.2)/\sigma(R=0.4)$ and $\sigma(R=0.2)/\sigma(R=0.6)$ in pp collisions at $\sqrt{s}=5.02$ TeV.

The fraction of $\mathrm{D^{0}}$ jets over inclusive charged-particle jets in pp collisions at $\sqrt{s}=5.02$ TeV for $R=0.2$, $0.4$, and $0.6$.

Distributions of $z^{\mathrm{ch}}_{||}$-differential yield of charm jets tagged with $\mathrm{D^{0}}$ mesons normalised by the number of $\mathrm{D^{0}}$ jets within each distribution in pp collisions at $\sqrt{s}=13$ TeV in four jet-$p_{\mathrm{T}}$ intervals $5<p_{\mathrm{T,ch\ jet}}<7$ GeV/$c$, $7<p_{\mathrm{T,ch\ jet}}<10$ GeV/$c$, $10<p_{\mathrm{T,ch\ jet}}<15$ GeV/$c$, and $15<p_{\mathrm{T,ch\ jet}}<50$ GeV/$c$ for jet radii $R=0.2$, $0.4$, and $0.6$.

Distributions of $z^{\mathrm{ch}}_{||}$-differential yield of charm jets tagged with $\mathrm{D^{0}}$ mesons normalised by the number of $\mathrm{D^{0}}$ jets within each distribution in pp collisions at $\sqrt{s}=5.02$ TeV in four jet-$p_{\mathrm{T}}$ intervals $5<p_{\mathrm{T,ch\ jet}}<7$ GeV/$c$, $7<p_{\mathrm{T,ch\ jet}}<10$ GeV/$c$, $10<p_{\mathrm{T,ch\ jet}}<15$ GeV/$c$, and $15<p_{\mathrm{T,ch\ jet}}<50$ GeV/$c$ for jet radii $R=0.2$, $0.4$, and $0.6$.

$p_{\mathrm{T,ch\ jet}}$-differential cross section of charm jets tagged with $\mathrm{D^{0}}$ mesons for $R=0.3$ in pp collisions at $\sqrt{s}=5.02$ TeV.

The fraction of $\mathrm{D^{0}}$ jets over inclusive charged-particle jets in pp collisions at $\sqrt{s}=5.02$ TeV for $R=0.3$.

Distributions of $z^{\mathrm{ch}}_{||}$-differential yield of charm jets tagged with $\mathrm{D^{0}}$ mesons normalised by the number of $\mathrm{D^{0}}$ jets within each distribution in pp collisions at $\sqrt{s}=5.02$ TeV in four jet-$p_{\mathrm{T}}$ intervals $5<p_{\mathrm{T,ch\ jet}}<7$ GeV/$c$, $7<p_{\mathrm{T,ch\ jet}}<10$ GeV/$c$, $10<p_{\mathrm{T,ch\ jet}}<15$ GeV/$c$, and $15<p_{\mathrm{T,ch\ jet}}<50$ GeV/$c$ for jet radius $R=0.3$.

Fig. 6b: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Away and Toward regions after the subtraction of Number density $N_{\rm ch}$ and $\Sigma p_{\rm T}$ distributions in the transverse region for p-Pb collisions for $p_{\rm T} >$ 0.5 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. 7a: $<p_{\rm T}>$ as a function of $p_{\rm T}^{\rm trig}$ in Toward (left) and Away (right) regions after the subtraction of transverse region for pp collisions for $p_{\rm T} >$ 0.5 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. 7b: $<p_{\rm T}>$ as a function of $p_{\rm T}^{\rm trig}$ in Toward (left) and Away (right) regions after the subtraction of transverse region for p-Pb collisions for $p_{\rm T} >$ 0.5 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A1: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 0.15 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A2: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 1.0 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A3: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 0.15 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A4: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 1.0 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Groomed jet radius (scaled) $\theta_{{\mathrm{g}}}$ $60<p_{\mathrm{T}}^{\mathrm{ch\;jet}}<80$ GeV/$c$, soft drop $z_{\mathrm{cut}}=0.1, \beta=0$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding") no correlation information is specified ($\pm$ is always used).

Groomed jet radius (scaled) $\theta_{{\mathrm{g}}}$ $60<p_{\mathrm{T}}^{\mathrm{ch\;jet}}<80$ GeV/$c$, soft drop $z_{\mathrm{cut}}=0.1, \beta=1$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding") no correlation information is specified ($\pm$ is always used).

Groomed jet radius (scaled) $\theta_{{\mathrm{g}}}$ $60<p_{\mathrm{T}}^{\mathrm{ch\;jet}}<80$ GeV/$c$, soft drop $z_{\mathrm{cut}}=0.1, \beta=2$. Note: The first bin corresponds to the Soft Drop untagged fraction. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding") no correlation information is specified ($\pm$ is always used).

Groomed jet momentum splitting fraction $z_{{\mathrm{g}}}$ $60<p_{\mathrm{T}}^{\mathrm{ch\;jet}}<80$ GeV/$c$, dynamical grooming $a=0.1$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding") no correlation information is specified ($\pm$ is always used).

Groomed jet momentum splitting fraction $z_{{\mathrm{g}}}$ $60<p_{\mathrm{T}}^{\mathrm{ch\;jet}}<80$ GeV/$c$, dynamical grooming $a=1.0$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding") no correlation information is specified ($\pm$ is always used).

Groomed jet momentum splitting fraction $z_{{\mathrm{g}}}$ $60<p_{\mathrm{T}}^{\mathrm{ch\;jet}}<80$ GeV/$c$, dynamical grooming $a=2.0$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding") no correlation information is specified ($\pm$ is always used).

Groomed jet radius (scaled) $\theta_{{\mathrm{g}}}$ $60<p_{\mathrm{T}}^{\mathrm{ch\;jet}}<80$ GeV/$c$, dynamical grooming $a=0.1$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding") no correlation information is specified ($\pm$ is always used).

Groomed jet radius (scaled) $\theta_{{\mathrm{g}}}$ $60<p_{\mathrm{T}}^{\mathrm{ch\;jet}}<80$ GeV/$c$, dynamical grooming $a=1.0$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding") no correlation information is specified ($\pm$ is always used).

Groomed jet radius (scaled) $\theta_{{\mathrm{g}}}$ $60<p_{\mathrm{T}}^{\mathrm{ch\;jet}}<80$ GeV/$c$, dynamical grooming $a=2.0$. For the "trkeff" and "generator" systematic uncertainty sources, the signed systematic uncertainty breakdowns ($\pm$ vs. $\mp$), denote correlation across bins (both within this table, and across tables). For the remaining sources ("unfolding") no correlation information is specified ($\pm$ is always used).

(Anti)proton spectrum in V0M multiplicity class IV - V

(Anti)proton spectrum in V0M multiplicity class VI

(Anti)proton spectrum in V0M multiplicity class VII

(Anti)proton spectrum in V0M multiplicity class VIII

(Anti)proton spectrum in V0M multiplicity class IX

(Anti)proton spectrum in V0M multiplicity class X

(Anti)proton spectrum in V0M multiplicity class INEL > 0

(Anti)deuteron spectrum in V0M multiplicity class I

(Anti)deuteron spectrum in V0M multiplicity class II

(Anti)deuteron spectrum in V0M multiplicity class III

(Anti)deuteron spectrum in V0M multiplicity class IV - V

(Anti)deuteron spectrum in V0M multiplicity class VI

(Anti)deuteron spectrum in V0M multiplicity class VII

(Anti)deuteron spectrum in V0M multiplicity class VIII

(Anti)deuteron spectrum in V0M multiplicity class IX

(Anti)deuteron spectrum in V0M multiplicity class X

(Anti)deuteron spectrum in V0M multiplicity class INEL > 0

(Anti)helion spectrum in V0M multiplicity class I - III

(Anti)helion spectrum in V0M multiplicity class IV - X

(Anti)helion spectrum in V0M multiplicity class INEL > 0

Mean transverse momentum $<p_\mathrm{T}>$ of protons as a function of the charged particle multiplicity density at mid-rapidity

Mean transverse momentum $<p_\mathrm{T}>$ of deuterons as a function of the charged particle multiplicity density at mid-rapidity

Mean transverse momentum $<p_\mathrm{T}>$ of helions as a function of the charged particle multiplicity density at mid-rapidity

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class I

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class II

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class III

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class IV - V

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class VI

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class VII

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class VIII

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class IX

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class X

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class INEL > 0

Coalescence parameter $B_3$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class I - III

Coalescence parameter $B_3$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class IV - X

Coalescence parameter $B_3$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class INEL > 0

Coalescence parameter $B_2$ as a function of the charged particle multiplicity density at mid-rapidity

Coalescence parameter $B_3$ as a function of the charged particle multiplicity density at mid-rapidity

Deuteron over proton ratio as a function of the charged particle multiplicity density at mid-rapidity

Helion over proton ratio as a function of the charged particle multiplicity density at mid-rapidity

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) with $3 < p_{\rm T} < 5$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} =$ 5 TeV and 13 TeV, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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$^{+}$, and D$^{*+}$) with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} =$ 5 TeV and 13 TeV, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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$^{+}$, and D$^{*+}$) with $16 < p_{\rm T} < 24$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} =$ 5 TeV and 13 TeV, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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, 7, and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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 TeV and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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, 7, and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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 TeV and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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 TeV and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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 TeV and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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 TeV and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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 TeV and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the azimuthal-correlation distributions of of D$^{0}$ mesons with $3 < p_{\rm T} < 5$ GeV/$c$ and charged particles with $p_{\rm T} > 0.3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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$^{0}$ mesons with $5 < p_{\rm T} < 8$ GeV/$c$ and charged particles with $p_{\rm T} > 0.3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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$^{0}$ mesons with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $p_{\rm T} > 0.3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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$^{0}$ mesons 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} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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$^{0}$ mesons with $3 < p_{\rm T} < 5$ GeV/$c$ and charged particles with $0.3 < p_{\rm T} < 1$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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$^{0}$ mesons with $5 < p_{\rm T} < 8$ GeV/$c$ and charged particles with $0.3 < p_{\rm T} < 1$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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$^{0}$ mesons with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $0.3 < p_{\rm T} < 1$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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$^{0}$ mesons with $16 < p_{\rm T} < 24$ GeV/$c$ and charged particles with $0.3 < p_{\rm T} < 1$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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$^{0}$ mesons with $3 < p_{\rm T} < 5$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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$^{0}$ mesons with $5 < p_{\rm T} < 8$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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$^{0}$ mesons with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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$^{0}$ mesons with $16 < p_{\rm T} < 24$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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$^{0}$ mesons with $3 < p_{\rm T} < 5$ GeV/$c$ and charged particles with $2 < p_{\rm T} < 3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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$^{0}$ mesons with $5 < p_{\rm T} < 8$ GeV/$c$ and charged particles with $2 < p_{\rm T} < 3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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$^{0}$ mesons with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $2 < p_{\rm T} < 3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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$^{0}$ mesons with $16 < p_{\rm T} < 24$ GeV/$c$ and charged particles with $2 < p_{\rm T} < 3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. 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 among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons 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 is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak yields among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak yields among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak yields among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2 < p_{\rm T} < 3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak width among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons 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 is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak width among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak width among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak width among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2 < p_{\rm T} < 3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak yields among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons 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 is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak width among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons 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 is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with 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 is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2 < p_{\rm T} < 3$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with 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 is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2 < p_{\rm T} < 3$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with 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 is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2 < p_{\rm T} < 3$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with 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 is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2 < p_{\rm T} < 3$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Transverse momentum distributions d$N$/(d$y$d$p_{\mathrm{T}}$) of prompt D$^{+}$ (middle) meson in the 30-50% (diamonds) centrality classes in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in |y|<0.5. The uncertainties on the BRs are quoted separately and the horizontal bars represent bin widths. The points in the figure of the paper are scaled for visibility by a factor 10$^{-1}$. Branching ratio of D$^+$ --> K$\pi\pi$ = 0.0938

Transverse momentum distributions d$N$/(d$y$d$p_{\mathrm{T}}$) of prompt D*$^{+}$ (right) meson in the 0-10% (crosses) centrality classes in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in |y|<0.5. The uncertainties on the BRs are quoted separately and the horizontal bars represent bin widths. Branching ratio of D*$^+$ --> D$^0\pi$ = 0.677

Transverse momentum distributions d$N$/(d$y$d$p_{\mathrm{T}}$) of prompt D*$^{+}$ (right) meson in the 30-50% (diamonds) centrality classes in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in |y|<0.5. The uncertainties on the BRs are quoted separately and the horizontal bars represent bin widths. The points in the figure of the paper are scaled for visibility by a factor 10$^{-1}$. Branching ratio of D*$^+$ --> D$^0\pi$ = 0.677

Ratio D$^{+}$/D$^{0}$ in central Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in |y|<0.5 and as a function of $p_\mathrm{T}$. Statistical (bars) and systematic (boxes) uncertainties are shown. Branching ratio of D$^0$ --> K$\pi$ = 0.0395 Branching ratio of D$^+$ --> K$\pi\pi$ = 0.0938

Ratio D$^{+}$/D$^{0}$ in semicentral Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in |y|<0.5 and as a function of $p_\mathrm{T}$. Statistical (bars) and systematic (boxes) uncertainties are shown. Branching ratio of D$^0$ --> K$\pi$ = 0.0395 Branching ratio of D$^+$ --> K$\pi\pi$ = 0.0938

Ratio D*$^{+}$/D$^{0}$ in central Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in |y|<0.5 and as a function of $p_\mathrm{T}$. Statistical (bars) and systematic (boxes) uncertainties are shown. Branching ratio of D$^0$ --> K$\pi$ = 0.0395 Branching ratio of D*$^+$ --> D$^0\pi$ = 0.677

Ratio D*$^{+}$/D$^{0}$ in semicentral Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in |y|<0.5 and as a function of $p_\mathrm{T}$. Statistical (bars) and systematic (boxes) uncertainties are shown. Branching ratio of D$^0$ --> K$\pi$ = 0.0395 Branching ratio of D*$^+$ --> D$^0\pi$ = 0.677

Nuclear modification factor of D$^{0}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 0-10% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Nuclear modification factor of D$^{+}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 0-10% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Nuclear modification factor of D*$^{+}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 0-10% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Nuclear modification factor of D$^{0}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 30-50% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Nuclear modification factor of D$^{+}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 30-50% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Nuclear modification factor of D*$^{+}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 30-50% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Average nuclear modification factor of D$^{0}$, D$^{+}$ and D*$^{+}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 0-10% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Average nuclear modification factor of D$^{0}$, D$^{+}$ and D*$^{+}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 30-50% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

$p_\mathrm{T}$-integrated nuclear modification factor of D$^{0}$ meson in central (0-10%) and semicentral (30-50%) Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in |y|<0.5.

Fig. 11 right: The $p_\mathrm{T}$ differential production cross section of charged-particle anti-$k_{T}$ $R=0.4$ b jets in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}~=~5.02$ TeV using the IP method.

Fig. 12 left: The combined differential production cross section of charged-particle anti-$k_{T}$ $R=0.4$ b jets in pp collisions at $\sqrt{s}~=~5.02$ TeV.

Fig. 12 right: The combined differential production cross section of charged-particle anti-$k_{T}$ $R=0.4$ b jets in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}~=~5.02$ TeV.

Fig. 13 left: The b-jet fraction in pp collisions at $\sqrt{s}~=~5.02$ TeV.

Fig. 13 right: The b-jet fraction in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}~=~5.02$ TeV.

Fig. 14 left: The nuclear modification factor $R_{\mathrm{pPb}}$ of b jets as a function of $P_T$ from the SV method.

Fig. 14 left: The nuclear modification factor $R_{\mathrm{pPb}}$ of b jets as a function of $P_T$ from the IP method.

Fig. 14 right: The combined nuclear modification factor $R_{\mathrm{pPb}}$ of b jets as a function of $P_T$.

Rapidity Differential production cross sections of J/$\psi$ in the range 2 $\leq$ p$_{\rm T}$ < 4 GeV/c

Rapidity Differential production cross sections of J/$\psi$ in the range 4 $\leq$ p$_{\rm T}$ < 6 GeV/c

Rapidity Differential production cross sections of J/$\psi$ in the range 6 $\leq$ p$_{\rm T}$ < 12 GeV/c

Differential production cross sections of $\psi$(2S) as a function of $p_{\rm T}$.

Differential production cross sections of $\psi$(2S) as a function of rapidity.

$\psi$(2S) over J/$\psi$ ratio as a function of $p_{\rm T}$.

$\psi$(2S) over J/$\psi$ ratio as a function of rapidity.

Differential production cross sections of $\Upsilon(1S)$ as a function of $p_{\rm T}$.

Differential production cross sections of $\Upsilon(1S)$ as a function of rapidity.

J/$\psi$ 5-to-13 TeV cross section ratio as a function of $p_{\rm T}$.

J/$\psi$ 7-to-13 TeV cross section ratio as a function of $p_{\rm T}$.

J/$\psi$ 8-to-13 TeV cross section ratio as a function of $p_{\rm T}$.

J/$\psi$ 5-to-13 TeV cross section ratio as a function of rapidity.

J/$\psi$ 7-to-13 TeV cross section ratio as a function of rapidity.

J/$\psi$ 8-to-13 TeV cross section ratio as a function of rapidity.

$\psi$(2S) 5-to-13 TeV cross section ratio as a function of rapidity.

$\psi$(2S) 7-to-13 TeV cross section ratio as a function of rapidity.

$\psi$(2S) 8-to-13 TeV cross section ratio as a function of rapidity.

$\psi$(2S) 5-to-13 TeV cross section ratio as a function of $p_{\rm T}$.

$\psi$(2S) 7-to-13 TeV cross section ratio as a function of $p_{\rm T}$.

$\psi$(2S) 8-to-13 TeV cross section ratio as a function of $p_{\rm T}$.

$\psi$(2S)-to-J/$\psi$ cross section ratio as a function of $\sqrt{s}$.

Integrated production cross section of J/$\psi$.

Integrated production cross section of $\psi$(2S).

Integrated production cross section of $\Upsilon$(1S).

Integrated production cross section of $\Upsilon$(2S).

Integrated production cross section of $\Upsilon$(3S).

$\Upsilon$(2S)-to-$\Upsilon$(1S) cross section ratio.

$\Upsilon$(3S)-to-$\Upsilon$(1S) cross section ratio.