The production cross section of inclusive J/$\psi$ pairs in pp collisions at a centre-of-mass energy $\sqrt{s} = 13$ TeV is measured with ALICE. The measurement is performed for J/$\psi$ in the rapidity interval $2.5 < y < 4.0$ and for transverse momentum $p_{\rm T} > 0$. The production cross section of inclusive J/$\psi$ pairs is reported to be $10.3 \pm 2.3 {\rm (stat.)} \pm 1.3 {\rm (syst.)}$ nb in this kinematic interval. The contribution from non-prompt J/$\psi$ (i.e. originated from beauty-hadron decays) to the inclusive sample is evaluated. The results are discussed and compared with data.

The transverse momentum ($p_{\rm T}$) differential cross section of the charm-strange baryon $\Xi^0_{\rm c}$ is measured at midrapidity ($|y|<$ 0.5) via its semileptonic decay into ${\rm e^{+}}\Xi^{-}\nu_{\rm e}$ in pp collisions at $\sqrt{s}$ = 5.02 TeV with the ALICE detector at the LHC. The ratio of the $p_{\rm T}$-differential $\Xi^0_{\rm c}$-baryon and ${\rm D^0}$-meson production cross sections is also reported. The measurements are compared with simulations with different tunes of the PYTHIA 8 event generator, with predictions from a statistical hadronisation model (SHM) with a largely augmented set of charm-baryon states beyond the current lists of the Particle Data Group, and with models including hadronisation via quark coalescence. The $p_{\rm T}$-integrated cross section of prompt $\Xi^0_{\rm c}$-baryon production at midrapidity is also reported, which is used to calculate the baryon-to-meson ratio $\Xi^0_{\rm c}/{\rm D^0} = 0.20 \pm 0.04~{\rm (stat.)} ^{+0.08}_{-0.07}~{\rm (syst.)}$. These results provide an additional indication of a modification of the charm fragmentation from $\rm e^+e^-$ and $\rm e^{-}p$ collisions to pp collisions.

Measurements of the production cross sections of prompt ${\rm D^0}$, ${\rm D^+}$, ${\rm D^{\ast +}}$, ${\rm D_s^+}$, ${\rm \Lambda_{c}^{+}}$, and ${\rm \Xi_{c}^{+}}$ charm hadrons at midrapidity in proton$-$proton collisions at $\sqrt{s}=13$ TeV with the ALICE detector are presented. The D-meson cross sections as a function of transverse momentum ($p_{\rm T}$) are provided with improved precision and granularity. The ratios of $p_{\rm T}$-differential meson production cross sections based on this publication and on measurements at different rapidity and collision energy provide a constraint on gluon parton distribution functions at low values of Bjorken-$x$ ($10^{-5}-10^{-4}$). The measurements of ${\rm \Lambda_{c}^{+}}$ (${\rm \Xi_{c}^{+}}$) baryon production extend the measured $p_{\rm T}$ intervals down to $p_{\rm T}=0(3)$~GeV$/c$. These measurements are used to determine the charm-quark fragmentation fractions and the ${\rm c\overline{c}}$ production cross section at midrapidity ($|y|<0.5$) based on the sum of the cross sections of the weakly-decaying ground-state charm hadrons ${\rm D^0}$, ${\rm D^+}$, ${\rm D_s^+}$, ${\rm \Lambda_{c}^{+}}$, ${\rm \Xi_{c}^{0}}$ and, for the first time, ${\rm \Xi_{c}^{+}}$, and of the strongly-decaying J/$psi$ mesons. The first measurements of ${\rm \Xi_{c}^{+}}$ and ${\rm \Sigma_{c}^{0,++}}$ fragmentation fractions at midrapidity are also reported. A significantly larger fraction of charm quarks hadronising to baryons is found compared to e$^+$e$^-$ and ep collisions. The ${\rm c\overline{c}}$ production cross section at midrapidity is found to be at the upper bound of state-of-the-art perturbative QCD calculations.

The first measurement of the production of pions, kaons, (anti-)protons and $\phi$ mesons at midrapidity in Xe-Xe collisions at $\sqrt{s_{\rm NN}} = 5.44$ TeV is presented. Transverse momentum ($p_{\rm T}$) spectra and $p_{\rm T}$-integrated yields are extracted in several centrality intervals bridging from p-Pb to mid-central Pb-Pb collisions in terms of final-state multiplicity. The study of Xe-Xe and Pb-Pb collisions allows systems at similar charged-particle multiplicities but with different initial geometrical eccentricities to be investigated. A detailed comparison of the spectral shapes in the two systems reveals an opposite behaviour for radial and elliptic flow. In particular, this study shows that the radial flow does not depend on the colliding system when compared at similar charged-particle multiplicity. In terms of hadron chemistry, the previously observed smooth evolution of particle ratios with multiplicity from small to large collision systems is also found to hold in Xe-Xe. In addition, our results confirm that two remarkable features of particle production at LHC energies are also valid in the collision of medium-sized nuclei: the lower proton-to-pion ratio with respect to the thermal model expectations and the increase of the $\phi$-to-pion ratio with increasing final-state multiplicity.

Neutral pion and $\eta$ meson production in the transverse momentum range 1 < $p_{T}$ < 20 GeV/$c$ have been measured at mid-rapidity by the ALICE experiment at the Large Hadron Collider (LHC) in central and semi-central Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV. These results were obtained using the photon conversion method as well as the PHOS and EMCal detectors. The results extend the upper $p_{T}$ reach of the previous ALICE $\pi^{0}$ measurements from 12 GeV/$c$ to 20 GeV/$c$ and present the first measurement of $\eta$ meson production in heavy-ion collisions at the LHC. The $\eta/\pi^{0}$ ratio is similar for the two centralities and reaches at high $p_{T}$ a plateau value of 0.457 $\pm$ 0.013$^{stat}$ $\pm$ 0.018$^{syst}$. A suppression of similar magnitude for $\pi^{0}$ and $\eta$ meson production is observed in Pb-Pb collisions with respect to their production in pp collisions scaled by the number of binary nucleon-nucleon collisions. We discuss the results in terms of NLO pQCD predictions and hydrodynamic models. The measurements show a stronger suppression with respect to what was observed at lower center-of-mass energies in the $p_{T}$ range 6 < $p_{T}$ < 10 GeV/$c$. At $p_{T}$ < 3 GeV/$c$, hadronization models describe the $\pi^{0}$ results while for the $\eta$ some tension is observed.

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.

The production of muons from heavy-flavour hadron decays in p-Pb collisions at $\sqrt{{\textit s}_{\rm NN}}=5.02$ TeV was studied for $2 < p_{\rm T} < 16$ GeV/$c$ with the ALICE detector at the CERN LHC. The measurement was performed at forward (p-going direction) and backward (Pb-going direction) rapidity, in the ranges of rapidity in the center-of-mass system (cms) $2.03<y_{\rm cms}<3.53$ and $-4.46<y_{\rm cms}<-2.96$, respectively. The production cross sections and nuclear modification factors are presented as a function of transverse momentum ($p_{\rm T}$). At forward rapidity, the nuclear modification factor is compatible with unity while at backward rapidity, in the interval $2.5<p_{\rm T}<3.5$ GeV/$c$, it is above unity by more than 2$\sigma$. The ratio of the forward-to-backward production cross sections is also measured in the overlapping interval $2.96 < \vert y_{\rm cms} \vert < 3.53$ and is smaller than unity by 3.7$\sigma$ in $2.5<p_{\rm T}<3.5$ GeV/$c$. The data are described by model calculations including cold nuclear matter effects.

The inclusive production of the charm-strange baryon $\rm \Omega_c^0$ is measured for the first time via its hadronic decay into $\rm \Omega^-\pi^+$ at midrapidity ($|y|<0.5$) in proton-proton (pp) collisions at the centre-of-mass energy $\sqrt{s}=13$ TeV with the ALICE detector at the LHC. The transverse momentum ($p_{\rm T}$) differential cross section multiplied by the branching ratio is presented in the interval $2<p_{\rm T}<12~{\rm GeV}/c$. The $p_{\rm T}$ dependence of the $\rm \Omega_c^0$-baryon production relative to the prompt $\rm D^0$-meson and to the prompt $\rm \Xi_c^0$-baryon production is compared to various models that take different hadronisation mechanisms into consideration. In the measured $p_{\rm T}$ interval, the ratio of the $p_{\rm T}$-integrated cross sections of $\rm \Omega_c^0$ and prompt $\Lambda_{\rm c}^{+}$ baryons multiplied by the $\rm \Omega^-\pi^+$ branching ratio is found to be larger by a factor of about 20 with a significance of about 4$\sigma$ when compared to $\rm e^+e^-$ collisions.

This paper presents the measurements of $\pi^{\pm}$, $\rm{K}^{\pm}$, $\rm{p}$ and $\bar{\rm{p}}$ transverse momentum ($p_{\rm{T}}$) spectra as a function of charged-particle multiplicity density in proton-proton (pp) collisions at $\sqrt{s}$ = 13 TeV with the ALICE detector at the LHC. Such study allows us to isolate the center-of-mass energy dependence of light-flavour particle production. The measurements reported here cover a $p_{\rm{T}}$ range from 0.1 GeV/$c$ to 20 GeV/$c$ and are done in the rapidity interval $|y|<0.5$. The $p_{\rm{T}}$-differential particle ratios exhibit an evolution with multiplicity, similar to that observed in pp collisions at $\sqrt{s}$ = 7 TeV, which is qualitatively described by some of the hydrodynamical and pQCD-inspired models discussed in this paper. Furthermore, the $p_{\rm{T}}$-integrated hadron-to-pion yield ratios measured in pp collisions at two different center-of-mass energies are consistent when compared at similar multiplicities. This also extends to strange and multistrange hadrons, suggesting that, at LHC energies, particle hadrochemistry scales with particle multiplicity the same way under different collision energies and colliding 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.