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.

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.

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.

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 production of $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $\mathrm{\phi(1020)}$ resonances has been measured in p-Pb collisions at $\sqrt{s_{\rm NN}}$ = 8.16 TeV using the ALICE detector. Resonances are reconstructed via their hadronic decay channels in the rapidity interval $-$0.5 $<$$y$$<$ 0 and the transverse momentum spectra are measured for various multiplicity classes up to $p_{\rm T}$ = 20 GeV/$c$ for $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $p_{\rm T}$ = 16 GeV/$c$ for $\mathrm{\phi(1020)}$. The $p_{\rm T}$ -integrated yields and mean transverse momenta are reported and compared with previous results in pp, p-Pb and Pb-Pb collisions. The $x_{\mathrm{T}}$ scaling for $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $\mathrm{\phi(1020)}$ resonance production is newly tested in p-Pb collisions and found to hold in the high-$p_{\rm T}$ region at LHC energies. The nuclear modification factors ($R_{\rm pPb}$) as a function of $p_{\rm T}$ for $\mathrm{K}^{*0}$ and $\mathrm{\phi}$ at $\sqrt{s_{NN}}$ = 8.16 TeV are presented along with the new $R_{\rm pPb}$ measurements of $\mathrm{K}^{*0}$, $\mathrm{\phi}$ , $\Xi$, and $\Omega$ at $\sqrt{s_{\rm NN}}$ = 5.02 TeV. At intermediate $p_{\rm T}$ (2-8 GeV/$c$), $R_{\rm pPb}$ of $\Xi$, $\Omega$ show a Cronin-like enhancement, while $\mathrm{K}^{*0}$ and $\mathrm{\phi}$ show no or little nuclear modification. At high $p_{\rm T}$ ($>$ 8 GeV/$c$), the $R_{\rm pPb}$ values of all hadrons are consistent with unity within uncertainties. The $R_{\rm pPb}$ of $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $\mathrm{\phi(1020)}$ at $\sqrt{s_{\rm NN}}$ = 8.16 and 5.02 TeV show no significant energy dependence.

The ALICE Collaboration reports three measurements in ultra-peripheral proton$-$lead collisions at forward rapidity. The exclusive two-photon process \ggmm and the exclusive photoproduction of J/$\psi$ are studied. J/$\psi$ photoproduction with proton dissociation is measured for the first time at a hadron collider. The cross section for the two-photon process of dimuons in the invariant mass range from 1 to 2.5 GeV/$c^2$ agrees with leading order quantum electrodynamics calculations. The exclusive and dissociative cross sections for J/$\psi$ photoproductions are measured for photon$-$proton centre-of-mass energies from 27 to 57 GeV. They are in good agreement with HERA results.

Inclusive $\psi$(2S) production is measured in p-Pb collisions at the centre-of-mass energy per nucleon-nucleon pair $\sqrt{s_{\rm{NN}}}=8.16$ TeV, using the ALICE detector at the CERN LHC. The production of $\psi$(2S) is studied at forward ($2.03 < y_{\rm cms} < 3.53$) and backward ($-4.46 < y_{\rm cms} < -2.96$) centre-of-mass rapidity and for transverse momentum $p_{\rm{T}}$ $<$ 12 GeV/$c$ via the decay to muon pairs. In this paper, we report the integrated as well as the $y_{\rm cms}$- and $p_{\rm T}$-differential inclusive production cross sections. Nuclear effects on $\psi$(2S) production are studied via the determination of the nuclear modification factor that shows a strong suppression at both forward and backward centre-of-mass rapidities. Comparisons with corresponding results for inclusive J/$\psi$ show a similar suppression for the two states at forward rapidity (p-going direction), but a stronger suppression for $\psi$(2S) at backward rapidity (Pb-going direction). As a function of $p_{\rm T}$, no clear dependence of the nuclear modification factor is found. The relative size of nuclear effects on $\psi$(2S) production compared to J/$\psi$ is also studied via the double ratio of production cross sections $[\sigma_{\psi(2S)}/\sigma_{\rm{J/}\psi}]_{\rm pPb}/[\sigma_{\psi(2S)}/\sigma_{\rm{J/}\psi}]_{\rm pp}$ between p-Pb and pp collisions. The results are compared with theoretical models that include various effects related to the initial and final state of the collision system and also with previous measurements at $\sqrt{s_{\rm{NN}}}$ = 5.02 TeV.

The energy deposited at very forward rapidities (very forward energy) is a powerful tool for characterising proton fragmentation in pp and p$-$Pb collisions. The correlation of very forward energy with particle production at midrapidity provides direct insights into the initial stages and the subsequent evolution of the collision. Furthermore, the correlation with the production of particles with large transverse momenta at midrapidity provides information complementary to the measurements of the underlying event, which are usually interpreted in the framework of models implementing centrality-dependent multiple parton interactions. Results about very forward energy, measured by the ALICE zero degree calorimeters (ZDCs), and its dependence on the activity measured at midrapidity in pp collisions at $\sqrt{s}=13$ TeV and in p$-$Pb collisions at $\sqrt{s_{\rm{NN}}}=8.16$ TeV are discussed. The measurements performed in pp collisions are compared with the expectations of three hadronic interaction event generators: PYTHIA 6 (Perugia 2011 tune), PYTHIA 8 (Monash tune), and EPOS LHC. These results provide new constraints on the validity of models in describing the beam remnants at very forward rapidities, where perturbative QCD cannot be used.

The production of $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $\mathrm{\phi(1020)}$ in pp collisions at $\sqrt{s}$ = 8 TeV was measured using Run 1 data collected by the ALICE collaboration at the LHC. The $p_{\rm{T}}$-differential yields d$^{\rm 2}N$/d$y$d$p_{\rm{T}}$ in the range $0 < p_{\rm{T}} < 20$ GeV/$c$ for $\rm{K}^{*0}$ and $0.4 < p_{\rm{T}} < 16$ GeV/$c$ for $\phi$ have been measured at midrapidity, $|y| < 0.5$. Moreover, improved measurements of the $\rm{K}^{*}(892)^{0}$ and $\phi(1020)$ at $\sqrt{s} = 7$TeV are presented. The collision energy dependence of $p_{\rm{T}}$ distributions, $p_{\rm{T}}$-integrated yields and particle ratios in inelastic pp collisions are examined. The results are also compared with different collision systems. The values of the particle ratios are found to be similar to those measured at other LHC energies. In pp collisions a hardening of the particle spectra is observed with increasing energy, but at the same time it is also observed that the relative particle abundances are independent of the collision energy. The $p_{\rm{T}}$-differential yields of $\rm{K}^{*0}$ and $\phi$ in pp collisions at $\sqrt{s} = 8$ TeV are compared with the expectations of different Monte Carlo event generators.