Differential cross section (multiplied by the dimuon branching fraction) of prompt J/$\psi$ mesons in pPb collisions at backward $y_{\mathrm{CM}}$.

Rapidity dependence of the cross section (multiplied by the dimuon branching fraction) for prompt J/$\psi$ mesons in the $p_{\mathrm{T}}$ intervals 6.5<$p_{\mathrm{T}}$<10 GeV and 10<$p_{\mathrm{T}}$<30 GeV in pp collisions.

Rapidity dependence of the cross section (multiplied by the dimuon branching fraction) for prompt J/$\psi$ mesons in the $p_{\mathrm{T}}$ intervals 6.5<$p_{\mathrm{T}}$<10 GeV and 10<$p_{\mathrm{T}}$<30 GeV in pPb collisions.

Transverse momentum dependence of $R_{\mathrm{pPb}}$ for prompt J/$\psi$ mesons in seven $y_{\mathrm{CM}}$ ranges.

Rapidity dependence of $R_{\mathrm{pPb}}$ for prompt J/$\psi$ mesons in $p_{\mathrm{T}}$ range 6.5<$p_{\mathrm{T}}$<10 GeV.

Rapidity dependence of $R_{\mathrm{pPb}}$ for prompt J/$\psi$ mesons in $p_{\mathrm{T}}$ range 10<$p_{\mathrm{T}}$<30 GeV.

Transverse momentum dependence of $R_{\mathrm{FB}}$ for prompt J/$\psi$ mesons in the 0<$y_{\mathrm{CM}}$<0.9 range.

Transverse momentum dependence of $R_{\mathrm{FB}}$ for prompt J/$\psi$ mesons in the 0.9<$y_{\mathrm{CM}}$<1.5 range.

Transverse momentum dependence of $R_{\mathrm{FB}}$ for prompt J/$\psi$ mesons in the 1.5<$y_{\mathrm{CM}}$<1.93 range.

Dependence of $R_{\mathrm{FB}}$ for prompt J/$\psi$ mesons on the hadronic activity in the event, given by the transverse energy deposited in the CMS detector at large pseudorapidities $E_{\mathrm{T}}^{\mathrm{HF}|\eta|>4}$.

Differential cross section (multiplied by the dimuon branching fraction) of nonprompt J/$\psi$ mesons in pp collisions at forward $y_{\mathrm{CM}}$.

Differential cross section (multiplied by the dimuon branching fraction) of nonprompt J/$\psi$ mesons in pp collisions at backward $y_{\mathrm{CM}}$.

Differential cross section (multiplied by the dimuon branching fraction) of nonprompt J/$\psi$ mesons in pPb collisions at forward $y_{\mathrm{CM}}$.

Differential cross section (multiplied by the dimuon branching fraction) of nonprompt J/$\psi$ mesons in pPb collisions at backward $y_{\mathrm{CM}}$.

Rapidity dependence of the cross section (multiplied by the dimuon branching fraction) for nonprompt J/$\psi$ mesons in the $p_{\mathrm{T}}$ intervals 6.5<$p_{\mathrm{T}}$<10 GeV and 10<$p_{\mathrm{T}}$<30 GeV in pp collisions.

Rapidity dependence of the cross section (multiplied by the dimuon branching fraction) for nonprompt J/$\psi$ mesons in the $p_{\mathrm{T}}$ intervals 6.5<$p_{\mathrm{T}}$<10 GeV and 10<$p_{\mathrm{T}}$<30 GeV in pPb collisions.

Transverse momentum dependence of $R_{\mathrm{pPb}}$ for nonprompt J/$\psi$ mesons in seven $y_{\mathrm{CM}}$ ranges.

Rapidity dependence of $R_{\mathrm{pPb}}$ for nonprompt J/$\psi$ mesons in $p_{\mathrm{T}}$ range 6.5<$p_{\mathrm{T}}$<10 GeV.

Rapidity dependence of $R_{\mathrm{pPb}}$ for nonprompt J/$\psi$ mesons in $p_{\mathrm{T}}$ range 10<$p_{\mathrm{T}}$<30 GeV.

Transverse momentum dependence of $R_{\mathrm{FB}}$ for nonprompt J/$\psi$ mesons in the 0<$y_{\mathrm{CM}}$<0.9 range.

Transverse momentum dependence of $R_{\mathrm{FB}}$ for nonprompt J/$\psi$ mesons in the 0.9<$y_{\mathrm{CM}}$<1.5 range.

Transverse momentum dependence of $R_{\mathrm{FB}}$ for nonprompt J/$\psi$ mesons in the 1.5<$y_{\mathrm{CM}}$<1.93 range.

Dependence of $R_{\mathrm{FB}}$ for nonprompt J/$\psi$ mesons on the hadronic activity in the event, given by the transverse energy deposited in the CMS detector at large pseudorapidities $E_{\mathrm{T}}^{\mathrm{HF}|\eta|>4}$.

Distribution of the transverse momentum ratio between the Z boson and the jet in PbPb collisions. The distribution is normalized by the number of Z events.

The mean value of the X_JZ distribution in pp collisions. The resolutions of the measured jet energy and azimuthal angle are smeared to match those of the PbPb sample.

The mean value of the X_JZ distribution in PbPb collisions.

The average number of jet partners per Z boson in pp collisions. The resolutions of the measured jet energy and azimuthal angle are smeared to match those of the PbPb sample.

The average number of jet partners per Z boson in PbPb collisions.

Distribution of the transverse momentum ratio between the Z boson and the jet in pp collisions. The distribution is normalized by the number of Z events. No resolution smearing is applied to the measured jet energy and azimuthal angle.

The $v_{2}$ result from SP method as a function of $p_{T}$ in 20-30\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{2}$ result from SP method as a function of $p_{T}$ in 30-40\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{2}$ result from SP method as a function of $p_{T}$ in 40-50\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{2}$ result from SP method as a function of $p_{T}$ in 50-60\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{3}$ result from SP method as a function of $p_{T}$ in 0-5\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{3}$ result from SP method as a function of $p_{T}$ in 5-10\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{3}$ result from SP method as a function of $p_{T}$ in 10-20\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{3}$ result from SP method as a function of $p_{T}$ in 20-30\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{3}$ result from SP method as a function of $p_{T}$ in 30-40\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{3}$ result from SP method as a function of $p_{T}$ in 40-50\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{3}$ result from SP method as a function of $p_{T}$ in 50-60\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{2}$ result from 4-, 6- and 8-particle cumulant methods as a function of $p_{T}$ in 5-10\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{2}$ result from 4-, 6- and 8-particle cumulant methods as a function of $p_{T}$ in 10-20\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{2}$ result from 4-, 6- and 8-particle cumulant methods as a function of $p_{T}$ in 20-30\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{2}$ result from 4-, 6- and 8-particle cumulant methods as a function of $p_{T}$ in 30-40\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{2}$ result from 4-, 6- and 8-particle cumulant methods as a function of $p_{T}$ in 40-50\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{2}$ result from 4-, 6- and 8-particle cumulant methods as a function of $p_{T}$ in 50-60\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{2}^{high}$ as a function of $v_{2}^{low}$ results from SP method in PbPb collisions at $sqrt{s_{NN}}$ = 5.02 TeV. Only statistical uncertainties are shown.

The $v_{2}^{high}$ as a function of $v_{2}^{low}$ results from 4-particle cumulant method in PbPb collisions at $sqrt{s_{NN}}$ = 5.02 TeV. Only statistical uncertainties are shown.

95% CL intervals on the double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of pT, for the forward rapidity analysis bin.

Double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of centrality, for the midrapidity analysis bin.

95% CL intervals on the double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of centrality, for the midrapidity analysis bin.

Double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of centrality, for the forward rapidity analysis bin.

95% CL intervals on the double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of centrality, for the forward rapidity analysis bin.

Double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, integrated over centrality, for the midrapidity and forward rapidity analysis bins.

95% CL interval on the double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, integrated over centrality, for the midrapidity and forward rapidity analysis bins.

Three-particle correlation with respect to the 2nd order event plane from Pb-going side in pPb collisions as a function of multiplicity with individual track pT between 0.3 to 3.0 GeV/c.

Three-particle correlation with respect to the 2nd order event plane in PbPb collisions as a function of multiplicity with individual track pT between 0.3 to 3.0 GeV/c.

Three-particle correlation with respect to the 2nd order event plane in PbPb collisions as a function of centrality with individual track pT between 0.3 to 3.0 GeV/c, where the average multiplicity value is plotted on figure in order to compare with pPb collisions.

Three-particle correlation with respect to the 2nd order event plane from Pb-going side, p-going side in pPb collisions as a function of multiplicity with individual track pT between 0.3 to 3.0 GeV/c.

Three-particle correlation with respect to the 2nd order event plane in PbPb collisions as a function of multiplicity with individual track pT between 0.3 to 3.0 GeV/c.

Three-particle correlation with respect to the 2nd order event plane in PbPb collisions as a function of centrality with individual track pT between 0.3 to 3.0 GeV/c, where the average multiplicity value is plotted on figure in order to compare with pPb collisions.

Three-particle correlation with respect to the 2nd order event plane from Pb-going and p-going direction in pPb collisions and in PbPb collisions as a function DeltaEta for multiplicity [185,220) with individual track pT between 0.3 to 3.0 GeV/c. Data points are plotted at the bin center.

The second-order Fourier coefficients, $V_{3\Delta}(2, |\Delta\eta| > 2)$, as a function of $N_{offline}^{trk}$ for charged particles, after correcting for back-to-back jet correlations, estimated from the 10 $\leq$ $N_{offline}^{trk}$ < 20 range.

The elliptic flow after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of $N_{offline}^{trk}$ for charged particles.

The elliptic flow after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of $N_{offline}^{trk}$ for charged particles.

The elliptic flow after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of $N_{offline}^{trk}$ for charged particles.

The triangular flow after correcting for back-to-back jet correlations, $v_{3}^{sub}(2, |\Delta\eta| > 2)$, as a function of $N_{offline}^{trk}$ for charged particles.

The triangular flow after correcting for back-to-back jet correlations, $v_{3}^{sub}(2, |\Delta\eta| > 2)$, as a function of $N_{offline}^{trk}$ for charged particles.

The triangular flow after correcting for back-to-back jet correlations, $v_{3}^{sub}(2, |\Delta\eta| > 2)$, as a function of $N_{offline}^{trk}$ for charged particles.

The elliptic flow, $v_{2}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for charged particles.

The elliptic flow, $v_{2}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for charged particles.

The elliptic flow, $v_{2}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for charged particles.

The elliptic flow after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for charged particles.

The elliptic flow after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for charged particles.

The elliptic flow after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for charged particles.

The elliptic flow, $v_{2}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for charged particles.

The elliptic flow, $v_{2}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for $K^{0}_{S}$.

The elliptic flow, $v_{2}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for $\Lambda/\bar{\Lambda}$.

The elliptic flow, $v_{2}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for charged particles.

The elliptic flow, $v_{2}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for $K^{0}_{S}$.

The elliptic flow, $v_{2}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for $\Lambda/\bar{\Lambda}$.

The elliptic flow after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for charged particles.

The elliptic flow after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for $K^{0}_{S}$.

The elliptic flow after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for $\Lambda/\bar{\Lambda}$.

The elliptic flow per constituent quark after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)/n_{q}$, as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ for $K^{0}_{S}$.

The elliptic flow per constituent quark after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ for $\Lambda/\bar{\Lambda}$.

The four-particle cumulant, $c_{2}(4)$, as a function of $N_{offline}^{trk}$ for charged particles.

The four-particle cumulant, $c_{2}(4)$, as a function of $N_{offline}^{trk}$ for charged particles.

The four-particle cumulant, $c_{2}(4)$, as a function of $N_{offline}^{trk}$ for charged particles.

The six-particle cumulant, $c_{2}(6)$, as a function of $N_{offline}^{trk}$ for charged particles.

The elliptic flow after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of $N_{offline}^{trk}$ for charged particles.

The elliptic flow, $v_{2}(4)$, as a function of $N_{offline}^{trk}$ for charged particles.

The elliptic flow, $v_{2}(6)$, as a function of $N_{offline}^{trk}$ for charged particles.

Distributions of the b-tagging efficiency as a function of the mistag rate of charm jets for pp collisions in a PYTHIA simulation.

Distributions of the b-tagging efficiency as a function of the mistag rate of light jets for pPb collisions in a PYTHIA+HIJING simulation.

Distributions of the b-tagging efficiency as a function of the mistag rate of charm jets for pPb collisions in a PYTHIA+HIJING simulation.

Distributions of the Secondary Vertex Mass after applying the SSV tagger selection.

b-jet cross-sections in pPb, scaled by the Glauber nuclear overlap function, as well as a b-jet cross-section from a PYTHIA simulation from -2.5 < eta_CM < -1.5.

b-jet cross-sections in pPb, scaled by the Glauber nuclear overlap function, as well as a b-jet cross-section from a PYTHIA simulation from -1.5 < eta_CM < -0.5 (x 10^2).

b-jet cross-sections in pPb, scaled by the Glauber nuclear overlap function, as well as a b-jet cross-section from a PYTHIA simulation from -0.5 < eta_CM < 0.5 (x10^4).

b-jet cross-sections in pPb, scaled by the Glauber nuclear overlap function, as well as a b-jet cross-section from a PYTHIA simulation from 0.5 < eta_CM < 1.5 (x10^6).

b-jet RpA (PYTHIA) from 0.5 < eta_CM < 1.5.

b-jet RpA (PYTHIA) from -0.5 < eta_CM < 0.5.

b-jet RpA (PYTHIA) from -1.5 < eta_CM < -0.5.

b-jet RpA (PYTHIA) from -2.5 < eta_CM < -1.5.

b-jet Fraction.

b-jet RpA(PYTHIA).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and $N_{offline}^{trk}<$ 35 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and $N_{offline}^{trk}<$ 35 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and $N_{offline}^{trk}<$ 35 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and $N_{offline}^{trk}<$ 35 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and $N_{offline}^{trk}<$ 35 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 0.1 $<p_{T}<$ 1.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 80 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 0.1 $<p_{T}<$ 1.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 90 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 80 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 90 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 80 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 90 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 80 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 90 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 0.1 $<p_{T}<$ 1.0 $GeV/c$ and 80 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 0.1 $<p_{T}<$ 1.0 $GeV/c$ and 90 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and 80 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and 90 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and 80 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and 90 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and 80 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and 90 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 0.1 $<p_{T}<$ 1.0 $GeV/c$ and $N_{offline}^{trk}>$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 0.1 $<p_{T}<$ 1.0 $GeV/c$ and $N_{offline}^{trk}>$ 110 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and $N_{offline}^{trk}>$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and $N_{offline}^{trk}>$ 110 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and $N_{offline}^{trk}>$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and $N_{offline}^{trk}>$ 110 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and $N_{offline}^{trk}>$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and $N_{offline}^{trk}>$ 110 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Associated yield for the near side of the correlation function averaged over 2 $<|\Delta\eta|<$ 4 and integrated over the region $|\Delta\Phi| < \Delta\Phi_{ZYAM}$ for pp data at $\sqrt{s} =$ 13 $TeV$. The associated yield as a function of $p_{T}$ for events with $N^{offline}_{trk} \geq$ 105. The $p_{T}$ value for each $p_{T}$ bin is the average $p_{T}$ value.

Associated yield for the near side of the correlation function averaged over 2 $<|\Delta\eta|<$ 4 for pp data at $\sqrt{s} =$ 7 $TeV$. The associated yield as a function of $p_{T}$ for events with $N^{offline}_{trk} \geq$ 110. The $p_{T}$ value for each $p_{T}$ bin is the average $p_{T}$ value.

Associated yield for the near side of the correlation function averaged over 2 $<|\Delta\eta|<$ 4 and integrated over the region $|\Delta\Phi| < \Delta\Phi_{ZYAM}$ for pp data at $\sqrt{s} =$ 13 $TeV$. The associated yield as a function of $N_{trk}^{offline}$ for events with 1.0 $< p_{T} <$ 2.0 GeV/c. The $N_{trk}^{offline}$ value for each $N_{trk}^{offline}$ bin is the average $N_{trk}^{offline}$ value.

Associated yield for the near side of the correlation function averaged over 2 $<|\Delta\eta|<$ 4 for pp data at $\sqrt{s} =$ 7 $TeV$. The associated yield as a function of $N_{trk}^{offline}$ for events with 1.0 $< p_{T} <$ 2.0 GeV/c. The $N_{trk}^{offline}$ value for each $N_{trk}^{offline}$ bin is the average $N_{trk}^{offline}$ value.

Associated yield for the near side of the correlation function averaged over 2 $<|\Delta\eta|<$ 4 pPb data at $\sqrt{s} =$ 5.02 $TeV$. The associated yield as a function of $N_{trk}^{offline}$ for events with 1.0 $< p_{T} <$ 2.0 GeV/c. The $N_{trk}^{offline}$ value for each $N_{trk}^{offline}$ bin is the average $N_{trk}^{offline}$ value.

Associated yield for the near side of the correlation function averaged over 2 $<|\Delta\eta|<$ 4 PbPb data at $\sqrt{s} =$ 2.76 $TeV$. The associated yield as a function of $N_{trk}^{offline}$ for events with 1.0 $< p_{T} <$ 2.0 GeV/c. The $N_{trk}^{offline}$ value for each $N_{trk}^{offline}$ bin is the average $N_{trk}^{offline}$ value.

The nuclear modification factors $R_{p+A}^{\rm{FONLL}}(p_{\rm{T}})$ of $B^{+}$ measured in $p$ + Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. The statistical and systematic uncertainties on the $p$ + Pb data are shown, followed by upper and lower systematic uncertainties from the FONLL predictions.

The nuclear modification factors $R_{p+A}^{\rm{FONLL}}(p_{\rm{T}})$ of $B^{0}$ measured in $p$ + Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. The statistical and systematic uncertainties on the $p$ + Pb data are shown, followed by upper and lower systematic uncertainties from the FONLL predictions.

The nuclear modification factors $R_{p+A}^{\rm{FONLL}}(p_{\rm{T}})$ of $B_{s}^{0}$ measured in $p$ + Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. The statistical and systematic uncertainties on the $p$ + Pb data are shown, followed by upper and lower systematic uncertainties from the FONLL predictions.

The measured $y_{\rm{CM}}$ (rapidity in the center-of-mass frame) differential production cross section of $B^{+}$ in $p$ + Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, together with the cross section calculated by the FONLL model.

The nuclear modification factors $R_{p+A}^{\rm{FONLL}}(y_{\rm{CM}})$ of $B^{+}$ measured in $p$ + Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. The statistical and systematic uncertainties on the $p$ + Pb data are shown, followed by upper and lower systematic uncertainties from the FONLL predictions.