The study of the azimuthal anisotropy of inclusive muons produced in p-Pb collisions at $\sqrt{s_{\rm NN}} = 8.16$ TeV, using the ALICE detector at the LHC is reported. The measurement of the second-order Fourier coefficient of the particle azimuthal distribution, $v_2$, is performed as a function of transverse momentum $p_{\rm T}$ in the 0-20% high-multiplicity interval at both forward ($2.03 < y_{\rm CMS} < 3.53$) and backward ($-4.46 < y_{\rm CMS} < -2.96$) rapidities over a wide $p_{\rm T}$ range, $0.5 < p_{\rm T} < 10$ GeV/$c$, in which a dominant contribution of muons from heavy-flavour hadron decays is expected at $p_{\rm T} > 2$ GeV/$c$. The $v_2$ coefficient of inclusive muons is extracted using two different techniques, namely two-particle cumulants, used for the first time for heavy-flavour measurements, and forward-central two-particle correlations. Both techniques give compatible results. A positive $v_2$ is measured at both forward and backward rapidities with a significance larger than $4.7\sigma$ and $7.6\sigma$, respectively, in the interval $2 < p_{\rm T} < 6$ GeV/$c$. Comparisons with previous measurements in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV, and with AMPT and CGC-based theoretical calculations are discussed. The findings impose new constraints on the theoretical interpretations of the origin of the collective behaviour in small collision systems.

Inclusive J/$\psi$ yields and average transverse momenta in p-Pb collisions at a center-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}}$ = 8.16 TeV are measured as a function of the charged-particle pseudorapidity density with ALICE. The J/$\psi$ mesons are reconstructed at forward $(2.03 < y_{\rm cms} < 3.53)$ and backward ($-4.46 < y_{\rm cms} < -2.96$) center-of-mass rapidity in their dimuon decay channel while the charged-particle pseudorapidity density is measured around midrapidity. The J/$\psi$ yields at forward and backward rapidity normalized to their respective average values increase with the normalized charged-particle pseudorapidity density, the former showing a weaker increase than the latter. The normalized average transverse momenta at forward and backward rapidity manifest a steady increase from low to high charged-particle pseudorapidity density with a saturation beyond the average value.

The production of J/$\psi$ is measured as a function of charged-particle multiplicity at forward rapidity in proton$-$proton (pp) collisions at center-of-mass energies $\sqrt{s} =$ 5.02 and 13 TeV. The J/$\psi$ mesons are reconstructed via their decay into dimuons in the rapidity interval (2.5 $< y <$ 4.0), whereas the charged-particle multiplicity density (${\rm d}N_{\rm{ch}}/{\rm d}\eta$) is measured at midrapidity $(|\eta| < 1)$. The production rate as a function of multiplicity is reported as the ratio of the yield in a given multiplicity interval to the multiplicity-integrated one. This observable shows a linear increase with charged-particle multiplicity normalized to the corresponding average value for inelastic events (${{\rm d}N_{\rm ch}/{\rm d}\eta}/{\langle {\rm d}N_{\rm ch}/{\rm d}\eta \rangle}$), at both the colliding energies. Measurements are compared with available ALICE results at midrapidity and theoretical model calculations. First measurement of the mean transverse momentum ($\langle p_{\mathrm{T}}\rangle$) of J/$\psi$ in pp collisions exhibits an increasing trend as a function of ${{\rm d}N_{\rm ch}/{\rm d}\eta}/{\langle {\rm d}N_{\rm ch}/{\rm d}\eta \rangle}$ showing a saturation towards high charged-particle multiplicities.

The product of the $\Lambda_b^0$ ($\overline{B}^0$) differential production cross-section and the branching fraction of the decay $\Lambda_b^0\rightarrow J/\psi pK^-$ ($\overline{B}^0\rightarrow J/\psi\overline{K}^*(892)^0$) is measured as a function of the beauty hadron transverse momentum, $p_{\rm T}$, and rapidity, $y$. The kinematic region of the measurements is $p_{\rm T}<20~{\rm GeV}/c$ and $2.0<y<4.5$. The measurements use a data sample corresponding to an integrated luminosity of $3~{\rm fb}^{-1}$ collected by the LHCb detector in $pp$ collisions at centre-of-mass energies $\sqrt{s}=7~{\rm TeV}$ in 2011 and $\sqrt{s}=8~{\rm TeV}$ in 2012. Based on previous LHCb results of the fragmentation fraction ratio, $f_{\Lambda_B^0}/f_d$, the branching fraction of the decay $\Lambda_b^0\rightarrow J/\psi pK^-$ is measured to be \begin{equation*} \mathcal{B}(\Lambda_b^0\rightarrow J/\psi pK^-)= (3.17\pm0.04\pm0.07\pm0.34^{+0.45}_{-0.28})\times10^{-4}, \end{equation*} where the first uncertainty is statistical, the second is systematic, the third is due to the uncertainty on the branching fraction of the decay $\overline{B}^0\rightarrow J/\psi\overline{K}^*(892)^0$, and the fourth is due to the knowledge of $f_{\Lambda_b^0}/f_d$. The sum of the asymmetries in the production and decay between $\Lambda_b^0$ and $\overline{\Lambda}_b^0$ is also measured as a function of $p_{\rm T}$ and $y$. The previously published branching fraction of $\Lambda_b^0\rightarrow J/\psi p\pi^-$, relative to that of $\Lambda_b^0\rightarrow J/\psi pK^-$, is updated. The branching fractions of $\Lambda_b^0\rightarrow P_c^+(\rightarrow J/\psi p)K^-$ are determined.

The inclusive production of the J/$\psi$ and $\psi$(2S) charmonium states is studied as a function of centrality in p-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}} = 8.16$ TeV at the LHC. The measurement is performed in the dimuon decay channel with the ALICE apparatus in the centre-of-mass rapidity intervals $-4.46 < y_{\rm cms} < -2.96$ (Pb-going direction) and $2.03 < y_{\rm cms} < 3.53$ (p-going direction), down to zero transverse momentum ($p_{\rm T}$). The J/$\psi$ and $\psi$(2S) production cross sections are evaluated as a function of the collision centrality, estimated through the energy deposited in the zero degree calorimeter located in the Pb-going direction. The $p_{\rm T}$-differential J/$\psi$ production cross section is measured at backward and forward rapidity for several centrality classes, together with the corresponding average $\langle p_{\rm T} \rangle$ and $\langle p^{2}_{\rm T} \rangle$ values. The nuclear effects affecting the production of both charmonium states are studied using the nuclear modification factor. In the p-going direction, a suppression of the production of both charmonium states is observed, which seems to increase from peripheral to central collisions. In the Pb-going direction, however, the centrality dependence is different for the two states: the nuclear modification factor of the J/$\psi$ increases from below unity in peripheral collisions to above unity in central collisions, while for the $\psi$(2S) it stays below or consistent with unity for all centralities with no significant centrality dependence. The results are compared with measurements in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV and no significant dependence on the energy of the collision is observed. Finally, the results are compared with theoretical models implementing various nuclear matter effects.

Neutral pion ($\pi^{0}$) and $\eta$ meson production cross sections were measured up to unprecedentedly high transverse momenta ($p_{\rm T}$) in p-Pb collisions at $\sqrt{s_{\rm NN}}$ = 8.16 TeV. The mesons were reconstructed via their two-photon decay channel in the rapidity interval $-1.3< y <0.3$ in the ranges of $0.4<p_{\rm T}<200$ GeV/$c$ and $1.0<p_{\rm T}<50$ GeV/$c$, respectively. The respective nuclear modification factor ($R_{\rm pPb}$) is presented for $p_{\rm T}$ up to of 200 and 30 GeV/$c$, where the former was achieved by extending the $\pi^{0}$ measurement in pp collisions at $\sqrt{s}$ = 8 TeV using the merged cluster technique. The values of $R_{\rm pPb}$ are below unity for $p_{\rm T}<10$ GeV/$c$, while they are consistent with unity for $p_{\rm T}>10$ GeV/$c$, leaving essentially no room for final state energy loss. The new data provide strong constraints for nuclear parton distribution and fragmentation functions over a broad kinematic range and are compared to model predictions as well as previous results at $\sqrt{s_{\rm NN}}$ = 5.02 TeV.