Nuclear modification factor of inclusive UPSI(1S) as a function of rapidity in p-Pb collisions at sqrt(s_NN) = 5.02 TeV. The data was collected in 2013 with two beam configurations, p-Pb and Pb-p with integrated luminosities of 5.0 nb-1 and 5.8 nb-1, respectively.

Nuclear modification factor of inclusive UPSI(1S) as a function of rapidity in p-Pb collisions at sqrt(sNN) = 5.02 TeV. The data was collected in 2013 with two beam configurations, p-Pb and Pb-p with integrated luminosities of 5.0 nb-1 and 5.8 nb-1, respectively.

Forward to backward ratio RFB of inclusive UPSI(1S) yields in p-Pb collisions at sqrt(sNN) = 5.02 TeV. The data was collected in 2013 with two beam configurations, p-Pb and Pb-p with integrated luminosities of 5.0 nb-1 and 5.8 nb-1, respectively.

Ratios of the cross sections of UPSI(2S)->MU+MU- to UPSI(1S)->MU+MU- in p-Pb collisions at sqrt(sNN) = 5.02 TeV. The data was collected in 2013 with two beam configurations, p-Pb and Pb-p with integrated luminosities of 5.0 nb-1 and 5.8 nb-1, respectively.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the $N_{offline}^{trk}$ < 35 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K_{S}^{0}$ as a function of $p_{T}$ from the correlation in the 35 $\leq N_{offline}^{trk}$ < 60 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K_{S}^{0}$ as a function of $p_{T}$ from the correlation in the 60 $\leq N_{offline}^{trk}$ < 120 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the $N_{offline}^{trk}$ < 35 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 35 $\leq N_{offline}^{trk}$ < 60 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 60 $\leq N_{offline}^{trk}$ < 120 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the $N_{offline}^{trk}$ < 35 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 35 $\leq N_{offline}^{trk}$ < 60 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 60 $\leq N_{offline}^{trk}$ < 120 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K_{S}^{0}$ as a function of $p_{T}$ from the correlation in the $N_{offline}^{trk}$ < 35 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K_{S}^{0}$ as a function of $p_{T}$ from the correlation in the 35 $\leq N_{offline}^{trk}$ < 60 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K_{S}^{0}$ as a function of $p_{T}$ from the correlation in the 60 $\leq N_{offline}^{trk}$ < 120 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the $N_{offline}^{trk}$ < 35 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 35 $\leq N_{offline}^{trk}$ < 60 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 60 $\leq N_{offline}^{trk}$ < 120 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 120 $\leq N_{offline}^{trk}$ < 150 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 150 $\leq N_{offline}^{trk}$ < 185 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 220 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 220 $\leq N_{offline}^{trk}$ < 260 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the 120 $\leq N_{offline}^{trk}$ < 150 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the 150 $\leq N_{offline}^{trk}$ < 185 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 220 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the 220 $\leq N_{offline}^{trk}$ < 260 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 120 $\leq N_{offline}^{trk}$ < 150 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 150 $\leq N_{offline}^{trk}$ < 185 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 220 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 220 $\leq N_{offline}^{trk}$ < 260 multiplicity class in pPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 120 $\leq N_{offline}^{trk}$ < 150 multiplicity class in pPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 150 $\leq N_{offline}^{trk}$ < 185 multiplicity class in pPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 220 multiplicity class in pPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 220 $\leq N_{offline}^{trk}$ < 260 multiplicity class in pPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $\Lambda/\bar{\Lambda}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 120 $\leq N_{offline}^{trk}$ < 150 multiplicity class in pPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $\Lambda/\bar{\Lambda}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 150 $\leq N_{offline}^{trk}$ < 185 multiplicity class in pPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $\Lambda/\bar{\Lambda}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 220 multiplicity class in pPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $\Lambda/\bar{\Lambda}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 220 $\leq N_{offline}^{trk}$ < 260 multiplicity class in pPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 120 $\leq N_{offline}^{trk}$ < 150 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 150 $\leq N_{offline}^{trk}$ < 185 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 220 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 220 $\leq N_{offline}^{trk}$ < 260 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the 120 $\leq N_{offline}^{trk}$ < 150 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the 150 $\leq N_{offline}^{trk}$ < 185 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 220 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the 220 $\leq N_{offline}^{trk}$ < 260 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 120 $\leq N_{offline}^{trk}$ < 150 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 150 $\leq N_{offline}^{trk}$ < 185 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 220 multiplicity class in PbPb.

The elliptic flow v2(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 220 $\leq N_{offline}^{trk}$ < 260 multiplicity class in PbPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 120 $\leq N_{offline}^{trk}$ < 150 multiplicity class in PbPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 150 $\leq N_{offline}^{trk}$ < 185 multiplicity class in PbPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 220 multiplicity class in PbPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 220 $\leq N_{offline}^{trk}$ < 260 multiplicity class in PbPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $\Lambda/\bar{\Lambda}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 120 $\leq N_{offline}^{trk}$ < 150 multiplicity class in PbPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $\Lambda/\bar{\Lambda}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 150 $\leq N_{offline}^{trk}$ < 185 multiplicity class in PbPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $\Lambda/\bar{\Lambda}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 220 multiplicity class in PbPb.

The elliptic flow per constituent quark v2(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $\Lambda/\bar{\Lambda}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 220 $\leq N_{offline}^{trk}$ < 260 multiplicity class in PbPb.

The triangular flow v3(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 350 multiplicity class in PbPb.

The triangular flow v3(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 350 multiplicity class in PbPb.

The triangular flow v3(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 350 multiplicity class in PbPb.

The triangular flow per constituent quark v3(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 350 multiplicity class in PbPb.

The triangular flow per constituent quark v3(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $\Lambda/\bar{\Lambda}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 350 multiplicity class in PbPb.

The triangular flow v3(2, $|\Delta\eta| > 2$) extracted for all charged particles as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 350 multiplicity class in pPb.

The triangular flow v3(2, $|\Delta\eta| > 2$) extracted for $K^{0}_{S}$ as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 350 multiplicity class in pPb.

The triangular flow v3(2, $|\Delta\eta| > 2$) extracted for $\Lambda/\bar{\Lambda}$ as a function of $p_{T}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 350 multiplicity class in pPb.

The triangular flow per constituent quark v3(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 350 multiplicity class in pPb.

The triangular flow per constituent quark v3(2, $|\Delta\eta| > 2$)/$n_{q}$ extracted for $\Lambda/\bar{\Lambda}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ from the correlation in the 185 $\leq N_{offline}^{trk}$ < 350 multiplicity class in pPb.

$v_{2,2}^{unsub}$ and $v_{2,2}$ as a function of $\Delta\eta$ calculated from the 2-D per-trigger yields in figure 4(a) and 4(b), respectively.

$v_{3,3}^{unsub}$ and $v_{3,3}$ as a function of $\Delta\eta$ calculated from the 2-D per-trigger yields in figure 4(a) and 4(b), respectively.

$v_{4,4}^{unsub}$ and $v_{4,4}$ as a function of $\Delta\eta$ calculated from the 2-D per-trigger yields in figure 4(a) and 4(b), respectively.

The per-trigger yield distributions $Y^{corr}(\Delta\phi)$ and $Y^{recoil}(\Delta\phi)$ for events with $N_{ch}^{rec} \geq$ 220 in the long-range region $|\Delta\eta| >$ 2.

The per-trigger yield distributions $Y^{corr}(\Delta\phi)$ and $Y^{recoil}(\Delta\phi)$ for events with $N_{ch}^{rec} \geq$ 220 in the long-range region $|\Delta\eta| >$ 2.

The per-trigger yield distributions $Y^{corr}(\Delta\phi)$ and $Y^{recoil}(\Delta\phi)$ for events with $N_{ch}^{rec} \geq$ 220 in the long-range region $|\Delta\eta| >$ 2.

The per-trigger yield distributions $Y^{corr}(\Delta\phi)$ and $Y^{recoil}(\Delta\phi)$ for events with $N_{ch}^{rec} \geq$ 220 in the long-range region $|\Delta\eta| >$ 2.

The per-trigger yield distributions $Y^{corr}(\Delta\phi)$ and $Y^{recoil}(\Delta\phi)$ for events with $N_{ch}^{rec} \geq$ 220 in the long-range region $|\Delta\eta| >$ 2.

The per-trigger yield distributions $Y^{corr}(\Delta\phi)$ and $Y^{recoil}(\Delta\phi)$ for events with $N_{ch}^{rec} \geq$ 220 in the long-range region $|\Delta\eta| >$ 2.

The per-trigger yield distributions $Y^{corr}(\Delta\phi)$ and $Y^{recoil}(\Delta\phi)$ for events with $N_{ch}^{rec} \geq$ 220 in the long-range region $|\Delta\eta| >$ 2.

Integrated per-trigger yield, $Y_{int}$, on the near-side as a function of $p_{T}^{a}$ for 1 $< p_{T}^{b} <$ 3 GeV.

Integrated per-trigger yield, $Y_{int}$, on the near-side as a function of $p_{T}^{a}$ for 1 $< p_{T}^{b} <$ 3 GeV.

Integrated per-trigger yield, $Y_{int}$, on the near-side as a function of $p_{T}^{a}$ for 1 $< p_{T}^{b} <$ 3 GeV.

Integrated per-trigger yield, $Y_{int}$, on the near-side as a function of $p_{T}^{a}$ for 1 $< p_{T}^{b} <$ 3 GeV.

Integrated per-trigger yield, $Y_{int}$, on the near-side as a function of $p_{T}^{a}$ for 1 $< p_{T}^{b} <$ 3 GeV.

Integrated per-trigger yield, $Y_{int}$, on the away-side as a function of $p_{T}^{a}$ for 1 $< p_{T}^{b} <$ 3 GeV.

Integrated per-trigger yield, $Y_{int}$, on the away-side as a function of $p_{T}^{a}$ for 1 $< p_{T}^{b} <$ 3 GeV.

Integrated per-trigger yield, $Y_{int}$, on the away-side as a function of $p_{T}^{a}$ for 1 $< p_{T}^{b} <$ 3 GeV.

Integrated per-trigger yield, $Y_{int}$, on the away-side as a function of $p_{T}^{a}$ for 1 $< p_{T}^{b} <$ 3 GeV.

Integrated per-trigger yield, $Y_{int}$, on the away-side as a function of $p_{T}^{a}$ for 1 $< p_{T}^{b} <$ 3 GeV.

The integrated per-trigger yield, Y_{int}, on the near-side, the away-side and their difference and Y_{int} from the recoil as a function of event activity. Errors are statistical.

The integrated per-trigger yield, Y_{int}, on the near-side, the away-side and their difference and Y_{int} from the recoil as a function of event activity. Errors are statistical.

The Fourier coefficients $v_{n}$ as a function of $p_{T}^{a}$ extracted from the correlation functions, before and after the subtraction of the recoil component.

The Fourier coefficients $v_{n}$ as a function of $p_{T}^{a}$ extracted from the correlation functions, before and after the subtraction of the recoil component.

The Fourier coefficients $v_{n}$ as a function of $p_{T}^{a}$ extracted from the correlation functions, before and after the subtraction of the recoil component.

$v_{2}$, $v_{3}$, $v_{4}$ and $v_{5}$ as a function of $p_T^a$ for 1 $< p_{T}^{b} <$ 3 GeV for different $N_{ch}^{rec}$ intervals.

$v_{2}$, $v_{3}$, $v_{4}$ and $v_{5}$ as a function of $p_T^a$ for 1 $< p_{T}^{b} <$ 3 GeV for different $N_{ch}^{rec}$ intervals.

$v_{2}$, $v_{3}$, $v_{4}$ and $v_{5}$ as a function of $p_T^a$ for 1 $< p_{T}^{b} <$ 3 GeV for different $N_{ch}^{rec}$ intervals.

$v_{2}$, $v_{3}$, $v_{4}$ and $v_{5}$ as a function of $p_T^a$ for 1 $< p_{T}^{b} <$ 3 GeV for different $N_{ch}^{rec}$ intervals.

$v_{2}$, $v_{3}$, $v_{4}$ and $v_{5}$ as a function of $p_T^a$ for 1 $< p_{T}^{b} <$ 3 GeV for different $N_{ch}^{rec}$ intervals.

$v_{2}$, $v_{3}$, $v_{4}$ and $v_{5}$ as a function of $p_T^a$ for 1 $< p_{T}^{b} <$ 3 GeV for different $N_{ch}^{rec}$ intervals.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{2}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The values of factorization variable $r_{3}$ defined by Eq.(11) before and after the subtraction of the recoil component. Errors are total experimental uncertainties.

The centrality dependence of $v_{2}$ as a function of $N_{ch}^{rec}$. Values from before and after the recoil subtraction are included.

The centrality dependence of $v_{3}$ as a function of $N_{ch}^{rec}$. Values from before and after the recoil subtraction are included.

The centrality dependence of $v_{4}$ as a function of $N_{ch}^{rec}$. Values from before and after the recoil subtraction are included.

The centrality dependence of $v_{2}$ as a function of $E_{T}^{Pb}$. Values from before and after the recoil subtraction are included.

The centrality dependence of $v_{3}$ as a function of $E_{T}^{Pb}$. Values from before and after the recoil subtraction are included.

The centrality dependence of $v_{4}$ as a function of $E_{T}^{Pb}$. Values from before and after the recoil subtraction are included.

The $v_{2}$ as a function of $E_{T}^{Pb}$ obtained indirectly by mapping from the $N_{ch}^{rec}-dependence of $v_{2}$ using the correlation data shown in Fig. 2(b).

The $v_{3}$ as a function of $E_{T}^{Pb}$ obtained indirectly by mapping from the $N_{ch}^{rec}-dependence of $v_{3}$ using the correlation data shown in Fig. 2(b).

The first-order harmonic of 2PC before recoil subtraction, $v_{1,1}^{unsub}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC before recoil subtraction, $v_{1,1}^{unsub}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC before recoil subtraction, $v_{1,1}^{unsub}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC before recoil subtraction, $v_{1,1}^{unsub}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC before recoil subtraction, $v_{1,1}^{unsub}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC before recoil subtraction, $v_{1,1}^{unsub}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC after recoil subtraction, $v_{1,1}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC after recoil subtraction, $v_{1,1}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC after recoil subtraction, $v_{1,1}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC after recoil subtraction, $v_{1,1}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC after recoil subtraction, $v_{1,1}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic of 2PC after recoil subtraction, $v_{1,1}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic $v_1$ obtained using factorization from $v_{1,1}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic $v_1$ obtained using factorization from $v_{1,1}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

The first-order harmonic $v_1$ obtained using factorization from $v_{1,1}$, as a function of $p_T^a$ for different $p_T^b$ ranges for events with $N_{ch}^{rec} \geq$ 220.

$v_{2}$ for Pb+Pb collisions in 55-60% centrality interval obtained using an EP method.

$v_{2}$ for Pb+Pb collisions in 55-60% centrality interval obtained using an EP method, after the scaling.

$v_{3}$ for Pb+Pb collisions in 55-60% centrality interval obtained using an EP method.

$v_{3}$ for Pb+Pb collisions in 55-60% centrality interval obtained using an EP method, after the scaling.

$v_{4}$ for Pb+Pb collisions in 55-60% centrality interval obtained using an EP method.

$v_{4}$ for Pb+Pb collisions in 55-60% centrality interval obtained using an EP method, after the scaling.

Correlation between $E_{T}^{FCal}$ and $N_{ch}^{rec}$ for MB events (without weighting) and MB+HMT events (with weighting), errors are statistical.

CL1 per-trigger yield vs midrapidity charged-particle multiplicity.

ZNA per-trigger yield vs midrapidity charged-particle multiplicity.

near-side per-trigger yield vs V0A multiplicity with and without long-range correlations subtraction.

away-side per-trigger yield vs V0A multiplicity with and without long-range correlations subtraction.

number of uncorrelated seeds vs midrapidity V0A multiplicity.

ratio between number of uncorrelated seeds and Glauber Ncoll vs V0A multiplicity classes.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

The Psi(2S) pT-differential cross sections in p-Pb, in the centre of mass rapidity range 2.03<y_cms<3.53. First uncertainty is statistical, the second one is the pT-correlated systematic, the third is the uncorrelated systematic. The fourth is a global uncertainty common to p-Pb and Pb-p results.

The Psi(2S) pT-differential cross sections in Pb-p, in the centre of mass rapidity range -4.46<y_cms<-2.96. First uncertainty is statistical, the second one is the pT-correlated systematic, the third the uncorrelated systematic. The fourth is a global uncertainty common to p-Pb and Pb-p results.

Double ratio [sigma(Psi(2S))/sigma(J/Psi)]_pA/[sigma(Psi(2S))/sigma(J/Psi)]_pp in p-Pb collisions as a function of pT. Results are obtained in the centre of mass rapidity range 2.03<y_cms<3.53. First uncertainty is statistical, the second one is the correlated systematic, the third the uncorrelated systematic.

Double ratio [sigma(Psi(2S))/sigma(J/Psi)]_pA/[sigma(Psi(2S))/sigma(J/Psi)]_pp in Pb-p collisions as a function of pT. Results are obtained in the centre of mass rapidity range -4.46<y_cms<-2.96. First uncertainty is statistical, the second one is the correlated systematic, the third the uncorrelated systematic.

The Psi(2S) RpA versus pT in p-Pb collisions. Results are obtained in the centre of mass rapidity range 2.03<y_cms<3.53. First uncertainty is statistical, the second one is the correlated systematic, the third the uncorrelated systematic. The fourth is a global uncertainty common to the JPsi and PsiP RpA.

The Psi(2S) RpA versus pT in Pb-p collisions. Results are obtained in the centre of mass rapidity range -4.46<y_cms<-2.96. First uncertainty is statistical, the second one is the correlated systematic, the third the uncorrelated systematic. The fourth is a global uncertainty common to the JPsi and PsiP RpA.

pT-differential cross section of prompt Ds+ mesons measured in p-Pb collisions in the rapidity interval -0.96<y_cms<0.04.

Nuclear modification factor R(pPb) of prompt D0 mesons vs. pT measured in p-Pb collisions in the rapidity interval -0.96<y_cms<0.04.

Nuclear modification factor R(pPb) of D+ mesons vs. pT measured in p-Pb collisions in the rapidity interval -0.96<y_cms<0.04.

Nuclear modification factor R(pPb) of D*+ mesons vs. pT measured in p-Pb collisions in the rapidity interval -0.96<y_cms<0.04.

Nuclear modification factor R(pPb) of Ds+ mesons vs. pT measured in p-Pb collisions in the rapidity interval -0.96<y_cms<0.04.

Average of D0, D+, D*+ nuclear modification factors R(pPb) vs. pT.

The nuclear modification factor of charged particles as a function of transverse momentum, measured in minimum-bias (NSD) p-Pb collisions at sqrt(sNN) = 5.02 TeV in two pseudorapidity ranges, |eta(cms)| < 0.3 and -1.3 < eta(cms) < 0.3.