Summary of results for cross-section of non-prompt $\psi(2S)$ decaying to a muon pair for 13 TeV data in fb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $J/\psi$ decaying to a muon pair as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $\psi(2S)$ decaying to a muon pair as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Correction factors for an assumed angular dependence $\propto 1+\lambda_{\theta} \cos^2\theta $ in the helicity frame, for several values of the parameter $\lambda_{\theta}$. The correction factors were found to be the same (within 1% - 2%) for $J/\psi$ and $\psi(2S)$ mesons, for prompt and non-prompt production mechanisms, and for the three rapidity intervals considered in this paper.

pT-differential cross section of inclusive Dzero mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV in the rapidity interval -0.96<y<0.04. Branching ratio of D0->Kpi : 0.0388.

pT-differential cross section of prompt Dzero mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV in the rapidity interval -0.96<y<0.04. Branching ratio of D0->Kpi : 0.0388. Data points for pt<2 GeV/c from analysis "without vertexing". Data points for pt>2 GeV/c from the analysis "with vertexing" taken from PRL 113 (2014) 232301 (http://hepdata.cedar.ac.uk/view/ins1296081).

First column: production cross sections per unit of rapidity for prompt D0 mesons, inclusive D0 mesons (no feed-down subtraction) and charm quarks in -0.96 < ycm < 0.04 in p-Pb collisions at 5.02 TeV. For D0 mesons, the second (sys) error is from the luminosity uncertainty, the third (sys) error is from the branching-ratio uncertainties. For charm quarks, the second (sys) error is from the luminosity uncertainty, the third (sys) error is from the Fragmentation Function uncertainties, the fourth (sys) error is from the rapidity shapes of D0 mesons and single charm quarks. Second column: value of <pT> of prompt D0 mesons in p-Pb collisions at 5.02 TeV in the rapidity interval -0.96 < y < 0.04. The first uncertainty is statistical, the second is the systematic uncertainty.

Ratio of D+ to D0 meson yield in p-Pb collisions at sqrt{sNN}=5.02 TeV in -0.96<y< 0.04 as a function of pT. Branching ratio of D+->Kpipi : 0.0913. Branching ratio of D0->Kpi : 0.0388.

Ratio of D*+ to D0 meson yield in p-Pb collisions at sqrt{sNN}=5.02 TeV in -0.96<y< 0.04 as a function of pT. Branching ratio of D*+->D0pi->Kpipi : 0.0388*0.677. Branching ratio of D0->Kpi : 0.0388.

Ratio of Ds to D0 meson yield in p-Pb collisions at sqrt{sNN}=5.02 TeV in -0.96<y< 0.04 as a function of pT. Branching ratio of Ds->phipi->KKpi : 0.0224. Branching ratio of D0->Kpi : 0.0388.

Ratio of Ds to D+ meson yield in p-Pb collisions at sqrt{sNN}=5.02 TeV in -0.96<y< 0.04 as a function of pT. Branching ratio of Ds->phipi->KKpi : 0.0224. Branching ratio of D+->Kpipi : 0.0913.

cross section of prompt Dzero mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 2<pt<5 GeV/c. Branching ratio of D0->Kpi : 0.0388. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dzero mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 5<pt<8 GeV/c. Branching ratio of D0->Kpi : 0.0388. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dzero mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 8<pt<16 GeV/c. Branching ratio of D0->Kpi : 0.0388. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dplus mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 2<pt<5 GeV/c. Branching ratio of D+->Kpipi : 0.0913. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dplus mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 5<pt<8 GeV/c. Branching ratio of D+->Kpipi : 0.0913. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dplus mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 8<pt<16 GeV/c. Branching ratio of D+->Kpipi : 0.0913. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dstar mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 2<pt<5 GeV/c. Branching ratio of D*+->D0pi->Kpipi : 0.0388*0.677. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dstar mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 5<pt<8 GeV/c. Branching ratio of D*+->D0pi->Kpipi : 0.0388*0.677. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dstar mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 8<pt<16 GeV/c. Branching ratio of D*+->D0pi->Kpipi : 0.0388*0.677. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

Nuclear modification factor RpPb of D0 mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV in -0.96<y<0.04 as a function of pT. Data point for 0<pt<1 GeV/c from analysis "without vertexing". Data points for pt>1 GeV/c from the analysis "with vertexing" taken from PRL 113 (2014) 232301 (http://hepdata.cedar.ac.uk/view/ins1296081).

Pt-integrated nuclear modification factor RpA of prompt D0 mesons in -0.96 < ycm < 0.04. The first error is the statistical uncertainty, the second is the systematic uncertainty.

$R_{pPb}$ as a function of $p_{T}$ integrated over rapidity range $-1.8<y*<1.3$ for the 0--90% centrality interval for the three geometric models: (a) Glauber, (b) Glauber-Gribov with $\omega_{\sigma}=0.11$ and (c) Glauber-Gribov with $\omega_{\sigma}=0.2$.

$R_{pPb}$ values as a function of $p_{T}$ in the 0--90% centrality interval averaged over $|y*|<0.5$.The $\langle T_{Pb} \rangle$ value for the ATLAS centrality correction is calculated with the Glauber model. The total systematic uncertainties, which include the uncertainty in $\langle T_{Pb} \rangle$, are indicated by lines of the same colour. Strict quantitative agreement is not expected as each measurement uses different rapidity intervals for the centrality determination and apply different event selection criteria to reject diffractive collisions.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $-1.8<y*<1.3$ for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $-1.8<y*<1.3$ for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $-1.8<y*<1.3$ for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{CP}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $|\eta|<2.3$ for seven centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{CP}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $|\eta|<2.3$ for seven centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{CP}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $|\eta|<2.3$ for seven centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

Centrality dependence of the Psi(2S) Q_pPb in the backward rapidity range (-4.46<y_cms<-2.96). The first uncertainty is statistical, the second is a systematic one. The third systematic uncertainty is fully correlated over centrality, while the fourth one is fully correlated over centrality and correlated between the JPsi and the Psi(2S) charmonium states.

Centrality dependence of the Psi(2S) Q_pPb in the forward rapidity range (2.03<y_cms<3.53). The first uncertainty is statistical, the second is the systematic one. The third systematic uncertainty is fully correlated over centrality, while the fourth one is fully correlated over centrality and correlated between the JPsi and the Psi(2S) charmonium states.

Psi(2S)/J/Psi)_pA/(Psi(2S)/J/Psi)_pp double ratio in the backward rapidity range (-4.46<y_cms<-2.96) plotted as a function of the crossing time. The first uncertainty is statistical, the second is the systematic one.

Psi(2S)/J/Psi)_pA/(Psi(2S)/J/Psi)_pp double ratio in the forward rapidity range (2.03<y_cms<3.53) plotted as a function of the crossing time. The first uncertainty is statistical, the second is the systematic one.

Average of relative D$^{0}$, D$^{+}$ and D*$^{+}$ meson yields for the sum of particles and antiparticles in several multiplicity and PT(D) intervals for p-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV as a function of the relative charged-particle multiplicity at central rapidity. The values are reported together with their uncertainties, which are quoted in the order: statistical, systematic and feed-down contribution uncertainties. The yields reported here are per inelastic event. The global normalisation uncertainty of $3.1\%$ is not shown.

Average of relative D$^{0}$, D$^{+}$ and D*$^{+}$ meson yields for the sum of particles and antiparticles in several multiplicity and PT(D) intervals for p-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV as a function of the relative average multiplicity in the VZERO detector, $N_{V0A} \big/ \langle N_{V0A} \rangle$. The uncertainties are shown in the following order: statistical, systematic, and feed-down contribution uncertainties. The global normalisation uncertainty of $3.1\%$ is not shown. The yields reported here are normalised to the inelastic cross section. For $1 < PT(D) < 2$ GeV/c and $12 < PT(D) < 24$ GeV/c, the final two multiplicity intervals are merged.

$p_{\rm T}$-differential yield of (K$^{*0}$ + $\overline{K^{*0}}$)/2 in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (20-40% multiplicity class).

$p_{\rm T}$-differential yield of (K$^{*0}$ + $\overline{K^{*0}}$)/2 in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (40-60% multiplicity class).

$p_{\rm T}$-differential yield of (K$^{*0}$ + $\overline{K^{*0}}$)/2 in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (60-80% multiplicity class).

$p_{\rm T}$-differential yield of (K$^{*0}$ + $\overline{K^{*0}}$)/2 in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (80-100multiplicity class).

$p_{\rm T}$-differential invariant yield of $phi$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (NSD). Additional systematic error: +- 3.1% (normalization).

$p_{\rm T}$-differential yield of $\phi$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (0-5% multiplicity class).

$p_{\rm T}$-differential yield of $\phi$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (5-10% multiplicity class).

$p_{\rm T}$-differential yield of $\phi$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (10-20% multiplicity class).

$p_{\rm T}$-differential yield of $\phi$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (20-40% multiplicity class).

$p_{\rm T}$-differential yield of $\phi$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (40-60% multiplicity class).

$p_{\rm T}$-differential yield of $\phi$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (60-80% multiplicity class).

$p_{\rm T}$-differential yield of $\phi$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (80-100% multiplicity class).

$p_{\rm T}$-integrated yield d$N$/d$y$ of (K$^{*0}$ + $\overline{K^{*0}}$)/2 in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV. Yields are normalised to the visible cross section for each V0A multiplicity event class and to NSD events for the minimum bias case (additional normalization systematic error: +- 3.1%).

Mean $p_{\rm T}$ of (K$^{*0}$ + $\overline{K^{*0}}$)/2 in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV for each V0A multiplicity event class.

$p_{\rm T}$-integrated yield d$N$/d$y$ of $phi$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV. Yields are normalised to the visible cross section for each V0A multiplicity event class and to NSD events for the minimum bias case (additional normalization systematic error: +- 3.1%).

Mean $p_{\rm T}$ of $phi$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV for each V0A multiplicity event class.

Mean $p_{\rm T}$ as function of the mass for particles measured with the ALICE detector in INEL pp collisions at $\sqrt{s}$ = 7 TeV in |y| < 0.5.

Mean $p_{\rm T}$ as function of the mass for particles measured with the ALICE detector in 0-20% p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV in -0.5 < y < 0.

Ratio of $p_{\rm T}$-differential yields of (p + $\bar{p}$) and $\phi$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (0-5$\%$ multiplicity class).

Ratio of $p_{\rm T}$-differential yields of (p + $\bar{p}$) and $\phi$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (5-10% multiplicity class).

Ratio of $p_{\rm T}$-differential yields of (p + $\bar{p}$) and $\phi$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (10-20% multiplicity class).

Ratio of $p_{\rm T}$-differential yields of (p + $\bar{p}$) and $\phi$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (20-40% multiplicity class).

Ratio of $p_{\rm T}$-differential yields of (p + $\bar{p}$) and $\phi$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (40-60% multiplicity class).

Ratio of $p_{\rm T}$-differential yields of (p + $\bar{p}$) and $\phi$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (60-80% multiplicity class).

Ratio of $p_{\rm T}$-differential yields of (p + $\bar{p}$) and $\phi$ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV (80-100% multiplicity class).

Ratio of $p_{\rm T}$-integrated yields of $(K^{*0}+\overline{K^{*0}})$ to long lived hadrons for different V0A multiplicity classes in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV.

Ratio of $p_{\rm T}$-integrated yields of 2$\phi$ to long lived hadrons for different V0A multiplicity classes in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV.

pT-differential invariant yield of charged kaons in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV, measured for NSD events. The first uncertainty is statistical, the second one is the total systematic uncertainty, while the third one is the uncorrelated systematic uncertainty which is multiplicity dependent.

pT-differential invariant yield of protons+antiprotons in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV, measured for different V0A multiplicity classes. The first uncertainty is statistical, the second one is the total systematic uncertainty, while the third one is the uncorrelated systematic uncertainty which is multiplicity dependent.

pT-differential invariant yield of protons+antiprotons in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV, measured for NSD events. The first uncertainty is statistical, the second one is the total systematic uncertainty, while the third one is the uncorrelated systematic uncertainty which is multiplicity dependent.

pT-differential invariant yield of charged pions in INEL p-p collisions with centre-of-mass energies 5.02 TeV and 7.0 TeV. The first uncertainty is statistical, the second is the systematic one.

pT-differential invariant yield of charged kaons in INEL p-p collisions with centre-of-mass energies 5.02 TeV and 7.0 TeV. The first uncertainty is statistical, the second is the systematic one.

pT-differential invariant yield of protons+antiprotons in INEL p-p collisions with centre-of-mass energies 5.02 TeV and 7.0 TeV. The first uncertainty is statistical, the second is the systematic one.

The nuclear modification factor RPPB of charged pions, charged kaons and protons+antiprotons as a function of pT in NSD p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV. The first uncertainty is statistical, the second is the systematic one.

pT-dependent charged kaon to charged pion production ratios in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV, measured for different V0A multiplicity classes. The first uncertainty is statistical, the second one is the total systematic uncertainty, while the third one is the uncorrelated systematic uncertainty which is multiplicity dependent.

pT-dependent proton+antiproton to charged pion production ratios in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV, measured for different V0A multiplicity classes. The first uncertainty is statistical, the second one is the total systematic uncertainty, while the third one is the uncorrelated systematic uncertainty which is multiplicity dependent.

pT-integrated charged kaon to charged pion and proton+antiproton to charged pion production ratios as a function of average multiplicities in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV, measured for different V0A multiplicity classes. The first uncertainty is statistical, the second one is the total systematic uncertainty, while the third one is the uncorrelated systematic uncertainty which is multiplicity dependent.

Xi pt spectrum in 20-40% multiplicity class.

Xi pt spectrum in 40-60% multiplicity class.

Xi pt spectrum in 60-80% multiplicity class.

Xi pt spectrum in 80-100% multiplicity class.

Omega pt spectrum in 0-5% multiplicity class.

Omega pt spectrum in 5-10% multiplicity class.

Omega pt spectrum in 10-20% multiplicity class.

Omega pt spectrum in 20-40% multiplicity class.

Omega pt spectrum in 40-60% multiplicity class.

Omega pt spectrum in 60-80% multiplicity class.

Omega pt spectrum in 80-100% multiplicity class.

Mean pT of XI- in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV for each V0A multiplicity event class.

Mean pT of XIBAR+ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV for each V0A multiplicity event class.

Mean pT of OMEGA- in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV for each V0A multiplicity event class.

Mean pT of OMEGABAR+ in p-Pb collisions with centre-of-mass energy/nucleon=5.02 TeV for each V0A multiplicity event class.

XiMinus yield as a function of charged particle multiplicity per unit pseudorapidity.

XiPlus yield as a function of charged particle multiplicity per unit pseudorapidity.

OmegaMinus yield as a function of charged particle multiplicity per unit pseudorapidity.

OmegaPlus yield as a function of charged particle multiplicity per unit pseudorapidity.

Xi/pi ratio as a function of charged particle multiplicity per unit pseudorapidity.

Omega/pi ratio as a function of charged particle multiplicity per unit pseudorapidity.

Summary of results for cross-section of non-prompt $J/\psi$ decaying to a muon pair for 8 TeV data in nb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for cross-section of prompt $\psi(2S)$ decaying to a muon pair for 7 TeV data in nb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for cross-section of prompt $\psi(2S)$ decaying to a muon pair for 8 TeV data in nb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for cross-section of non-prompt $\psi(2S)$ decaying to a muon pair for 7 TeV data in nb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for cross-section of non-prompt $\psi(2S)$ decaying to a muon pair for 8 TeV data in nb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $J/\psi$ decaying to a muon pair as a percentage for 7 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $J/\psi$ decaying to a muon pair as a percentage for 8 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $\psi(2S)$ decaying to a muon pair as a percentage for 7 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $\psi(2S)$ decaying to a muon pair as a percentage for 8 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 7 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 8 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 7 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 8 TeV data. Uncertainties are statistical and systematic, respectively.

Mean weight correction factor for $J/\psi$ under the "longitudinal" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $J/\psi$ under the "transverse zero" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $J/\psi$ under the "transverse positive" spin-alignment transverse positive hypothesis for 7 TeV.

Mean weight correction factor for $J/\psi$ under the "transverse negative" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $J/\psi$ under the "off-($\lambda_{\theta}$--$\lambda_{\phi}$)-plane positive" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $J/\psi$ under the "off-($\lambda_{\theta}$--$\lambda_{\phi}$)-plane negative" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $\psi(2S)$ under the "longitudinal" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $\psi(2S)$ under the "transverse zero" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $\psi(2S)$ under the "transverse positive" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $\psi(2S)$ under the "transverse negative" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $\psi(2S)$ under the "off-($\lambda_{\theta}$--$\lambda_{\phi}$)-plane positive" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $\psi(2S)$ under the "off-($\lambda_{\theta}$--$\lambda_{\phi}$)-plane negative" spin-alignment hypothesis for 7 TeV.

Mean weight correction factor for $J/\psi$ under the "longitudinal" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $J/\psi$ under the "transverse zero" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $J/\psi$ under the "transverse positive" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $J/\psi$ under the "transverse negative" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $J/\psi$ under the "off-($\lambda_{\theta}$--$\lambda_{\phi}$)-plane positive" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $J/\psi$ under the "off-($\lambda_{\theta}$--$\lambda_{\phi}$)-plane negative" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $\psi(2S)$ under the "longitudinal" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $\psi(2S)$ under the "transverse zero" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $\psi(2S)$ under the "transverse positive" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $\psi(2S)$ under the "transverse negative" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $\psi(2S)$ under the "off-($\lambda_{\theta}$--$\lambda_{\phi}$)-plane positive" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Mean weight correction factor for $\psi(2S)$ under the "off-($\lambda_{\theta}$--$\lambda_{\phi}$)-plane negative" spin-alignment hypothesis for 8 TeV. Those intervals not measured in the analysis at low $p_{T}$, high rapidity are also excluded here.

Raw number of J/$\psi$ (N$_{\rm AA}^{\rm J/\psi}$), expected raw number of hadronic J/$\psi$ (N$_{\rm AA}^{\rm{h J}/\psi}$) and measured excess in the number of J/$\psi$ (N$_{\rm AA}^{\rm{exess J}/\psi}$), all three numbers in the $p_{\rm T}$ range (0-0.3) GeV/$c$, and J/$\psi$ coherent photoproduction cross section, as a function of centrality, with their statistical and systematic uncertainties. An additional correlated systematic uncertainty also applies to the cross section. In the most central classes, an upper limit (95% CL) on the J/$\psi$ yield excess and on the cross section is given.