The production of the $\psi(2S)$ charmonium state was measured with ALICE in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV, in the dimuon decay channel. A significant signal was observed for the first time at LHC energies down to zero transverse momentum, at forward rapidity ($2.5<y<4$). The measurement of the ratio of the inclusive production cross sections of the $\psi(2S)$ and J/$\psi$ resonances is reported as a function of the centrality of the collisions and of transverse momentum, in the region $p_{\rm T}<12$ GeV/$c$. The results are compared with the corresponding measurements in pp collisions, by forming the double ratio $[\sigma^{\psi(2S)}/\sigma^{J/\psi}]_{\rm{Pb-Pb}}/[\sigma^{\psi(2S)}/\sigma^{J/\psi}]_{\rm{pp}}$. It is found that in Pb-Pb collisions the $\psi(2S)$ is suppressed by a factor of $\sim 2$ with respect to the J/$\psi$. The $\psi(2S)$ nuclear modification factor $R_{\rm AA}$ was also obtained as a function of both centrality and $p_{\rm T}$. The results show that the $\psi(2S)$ resonance yield is strongly suppressed in Pb-Pb collisions, by a factor up to $\sim 3$ with respect to pp. Comparisons of cross section ratios with previous SPS findings by the NA50 experiment and of $R_{\rm AA}$ with higher-$p_{\rm T}$ results at LHC energy are also reported. These results and the corresponding comparisons with calculations of transport and statistical models address questions on the presence and properties of charmonium states in the quark-gluon plasma formed in nuclear collisions at the LHC.

The cross section for coherent photonuclear production of J/$\psi$ is presented as a function of the electromagnetic dissociation (EMD) of Pb. The measurement is performed with the ALICE detector in ultra-peripheral Pb-Pb collisions at a centre-of-mass energy per nucleon pair of $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV. Cross sections are presented in five different J/$\psi$ rapidity ranges within $|y|<4$, with the J/$\psi$ reconstructed via its dilepton decay channels. In some events the J/$\psi$ is not accompanied by EMD, while other events do produce neutrons from EMD at beam rapidities either in one or the other beam direction, or in both. The cross sections in a given rapidity range and for different configurations of neutrons from EMD allow for the extraction of the energy dependence of this process in the range $17 < W_{\gamma\, \mathrm{Pb, n}} <920$ GeV, where $W_{\gamma\, \mathrm{Pb, n}}$ is the centre-of-mass energy per nucleon of the $\gamma\,\mathrm{Pb}$ system. This range corresponds to a Bjorken-$x$ interval spanning about three orders of magnitude: $ 1.1\times10^{-5}<x<3.3\times 10^{-2}$. In addition to the ultra-peripheral and photonuclear cross sections, the nuclear suppression factor is obtained. These measurements point to a strong depletion of the gluon distribution in Pb nuclei over a broad, previously unexplored, energy range. These results, together with previous ALICE measurements, provide unprecedented information to probe quantum chromodynamics at high energies.

The azimuthal ($\Delta\varphi$) correlation distributions between heavy-flavor decay electrons and associated charged particles are measured in pp and p$-$Pb collisions at $\sqrt{s_{\rm{NN}}} = 5.02$ TeV. Results are reported for electrons with transverse momentum $4<p_{\rm T}<16$ GeV/$c$ and pseudorapidity $|\eta|<0.6$. The associated charged particles are selected with transverse momentum $1<p_{\rm T}<7$ GeV/$c$, and relative pseudorapidity separation with the leading electron $|\Delta\eta| < 1$. The correlation measurements are performed to study and characterize the fragmentation and hadronization of heavy quarks. The correlation structures are fitted with a constant and two von Mises functions to obtain the baseline and the near- and away-side peaks, respectively. The results from p$-$Pb collisions are compared with those from pp collisions to study the effects of cold nuclear matter. In the measured trigger electron and associated particle kinematic regions, the two collision systems give consistent results. The $\Delta\varphi$ distribution and the peak observables in pp and p$-$Pb collisions are compared with calculations from various Monte Carlo event generators.

A study of multiplicity and pseudorapidity distributions of inclusive photons measured in pp and p$-$Pb collisions at a center-of-mass energy per nucleon$-$nucleon collision of $\sqrt{s_{\rm NN}} = 5.02$ TeV using the ALICE detector in the forward pseudorapidity region $2.3 < \eta_{\rm lab} < 3.9$ is presented. Measurements in p$-$Pb collisions are reported for two beam configurations in which the directions of the proton and lead ion beam were reversed. The pseudorapidity distributions in p$-$Pb collisions are obtained for seven centrality classes which are defined based on different event activity estimators, i.e., the charged-particle multiplicity measured at midrapidity as well as the energy deposited in a calorimeter at beam rapidity. The inclusive photon multiplicity distributions for both pp and p$-$Pb collisions are described by double negative binomial distributions. The pseudorapidity distributions of inclusive photons are compared to those of charged particles at midrapidity in \pp collisions and for different centrality classes in p$-$Pb collisions. The results are compared to predictions from various Monte Carlo event generators. None of the generators considered in this paper reproduces the inclusive photon multiplicity distributions in the reported multiplicity range. The pseudorapidity distributions are, however, better described by the same generators.

The measurement of the production of deuterons, tritons and $^{3}\mathrm{He}$ and their antiparticles in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV is presented in this article. The measurements are carried out at midrapidity ($|y| < $ 0.5) as a function of collision centrality using the ALICE detector. The $p_{\rm T}$-integrated yields, the coalescence parameters and the ratios to protons and antiprotons are reported and compared with nucleosynthesis models. The comparison of these results in different collision systems at different centre-of-mass collision energies reveals a suppression of nucleus production in small systems. In the Statistical Hadronisation Model framework, this can be explained by a small correlation volume where the baryon number is conserved, as already shown in previous fluctuation analyses. However, a different size of the correlation volume is required to describe the proton yields in the same data sets. The coalescence model can describe this suppression by the fact that the wave functions of the nuclei are large and the fireball size starts to become comparable and even much smaller than the actual nucleus at low multiplicities.

The fractions of non-prompt (i.e. originating from beauty-hadron decays) D$^0$ and D$^+$ mesons with respect to the inclusive yield are measured as a function of the charged-particle multiplicity in proton$-$proton collisions at a centre-of-mass energy of $\sqrt{s} = 13$ TeV with the ALICE detector at the LHC. The results are reported in intervals of transverse momentum ($p_{\rm T}$) and integrated in the range $1 < p_{\rm T} < 24$ GeV/$c$. The fraction of non-prompt D$^0$ and D$^+$ mesons is found to increase slightly as a function of $p_{\rm T}$ in all the measured multiplicity intervals, while no significant dependence on the charged-particle multiplicity is observed. In order to investigate the production and hadronisation mechanisms of charm and beauty quarks, the results are compared to PYTHIA 8 as well as EPOS 3 and EPOS 4 Monte Carlo simulations, and to calculations based on the colour glass condensate including three-pomeron fusion.

A newly developed observable for correlations between symmetry planes, which characterize the direction of the anisotropic emission of produced particles, is measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV with ALICE. This so-called Gaussian Estimator allows for the first time the study of these quantities without the influence of correlations between different flow amplitudes. The centrality dependence of various correlations between two, three and four symmetry planes is presented. The ordering of magnitude between these symmetry plane correlations is discussed and the results of the Gaussian Estimator are compared with measurements of previously used estimators. The results utilizing the new estimator lead to significantly smaller correlations than reported by studies using the Scalar Product method. Furthermore, the obtained symmetry plane correlations are compared to state-of-the-art hydrodynamic model calculations for the evolution of heavy-ion collisions. While the model predictions provide a qualitative description of the data, quantitative agreement is not always observed, particularly for correlators with significant non-linear response of the medium to initial state anisotropies of the collision system. As these results provide unique and independent information, their usage in future Bayesian analysis can further constrain our knowledge on the properties of the QCD matter produced in ultrarelativistic heavy-ion collisions.

Two-particle correlations with ${\rm K}^{0}_{\rm{S}}$, $\Lambda$/$\bar{\Lambda}$, and charged hadrons as trigger particles in the transverse momentum range $8 < p_\mathrm{T,trig}<16$ GeV/$c$, and associated charged particles within $1 < p_\mathrm{T,assoc}<8$ GeV/$c$, are studied at mid-rapidity in pp and central Pb-Pb collisions at a centre-of-mass energy per nucleon-nucleon collision $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV with the ALICE detector at the LHC. After subtracting the contributions of the flow background, the per-trigger yields are extracted on both the near and away sides, and the ratio in Pb-Pb collisions with respect to pp collisions ($I_{\mathrm {AA}}$) is computed. The per-trigger yield in Pb-Pb collisions on the away side is strongly suppressed to the level of $I_{\mathrm {AA}} \approx 0.6$ for $p_\mathrm{T,assoc}>3$ GeV/$c$ as expected from strong in-medium energy loss, while an enhancement develops at low $p_\mathrm{T,assoc}$ on both the near and away sides, reaching $I_{\mathrm {AA}} \approx 1.8$ and $2.7$ respectively. These findings are in good agreement with previous ALICE measurements from two-particle correlations triggered by neutral pions ($\pi^{0}$-h) and charged hadrons (h-h) in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 2.76$ TeV. Moreover, the correlations with ${\rm K}^{0}_{\rm{S}}$ mesons and $\Lambda$/$\bar{\Lambda}$ baryons as trigger particles are compared to those of inclusive charged hadrons. The results are compared with the predictions of Monte Carlo models.

The production of $\pi^\pm$, ${\rm K}^\pm$, and $(\overline{\rm p})$p is measured in pp collisions at $\sqrt{s}=13$ TeV in different topological regions. Particle transverse momentum ($p_{\rm T}$) spectra are measured in the ``toward'', ``transverse'', and ``away'' angular regions defined with respect to the direction of the leading particle in the event. While the toward and away regions contain the fragmentation products of the near-side and away-side jets, respectively, the transverse region is dominated by particles from the Underlying Event (UE). The relative transverse activity classifier, $R_{\rm T}=N_{\rm T}/\langle N_{\rm T}\rangle$, is used to group events according to their UE activity, where $N_{\rm T}$ is the measured charged-particle multiplicity per event in the transverse region and $\langle N_{\rm T}\rangle$ is the mean value over all the analysed events. The first measurements of identified particle $p_{\rm T}$ spectra as a function of $R_{\rm T}$ in the three topological regions are reported. The yield of high transverse momentum particles relative to the $R_{\rm T}$-integrated measurement decreases with increasing $R_{\rm T}$ in both the toward and away regions, indicating that the softer UE dominates particle production as $R_{\rm T}$ increases and validating that $R_{\rm T}$ can be used to control the magnitude of the UE. Conversely, the spectral shapes in the transverse region harden significantly with increasing $R_{\rm T}$. This hardening follows a mass ordering, being more significant for heavier particles. The $p_{\rm T}$-differential particle ratios $({\rm p+\overline{p}})/(\pi^+ +\pi^-)$ and $({\rm K^+ +K^-})/(\pi^+ +\pi^-)$ in the low UE limit $(R_{\rm T}\rightarrow 0)$ approach expectations from Monte Carlo generators such as PYTHIA 8 with Monash 2013 tune and EPOS LHC, where the jet-fragmentation models have been tuned to reproduce ${\rm e^+ e^-}$ results.

The production of strange hadrons (K$^{0}_{\rm S}$, $\Lambda$, $\Xi^{\pm}$, and $\Omega^{\pm}$), baryon-to-meson ratios ($\Lambda/{\rm K}^0_{\rm S}$, $\Xi/{\rm K}^0_{\rm S }$, and $\Omega/{\rm K}^0_{\rm S}$), and baryon-to-baryon ratios ($\Xi/\Lambda$, $\Omega/\Lambda$, and $\Omega/\Xi$) associated with jets and the underlying event were measured as a function of transverse momentum ($p_{\rm T}$) in pp collisions at $\sqrt{s} = 13$ TeV and p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV with the ALICE detector at the LHC. The inclusive production of the same particle species and the corresponding ratios are also reported. The production of multi-strange hadrons, $\Xi^{\pm}$ and $\Omega^{\pm}$, and their associated particle ratios in jets and in the underlying event are measured for the first time. In both pp and p-Pb collisions, the baryon-to-meson and baryon-to-baryon yield ratios measured in jets differ from the inclusive particle production for low and intermediate hadron $p_{\rm T}$ (0.6$-$6 GeV/$c$). Ratios measured in the underlying event are in turn similar to those measured for inclusive particle production. In pp collisions, the particle production in jets is compared with PYTHIA 8 predictions with three colour-reconnection implementation modes. None of them fully reproduces the data in the measured hadron $p_{\rm T}$ region. The maximum deviation is observed for $\Xi^{\pm}$ and $\Omega^{\pm}$, which reaches a factor of about six. In p-Pb collisions, there is no significant event-multiplicity dependence for particle production in jets, in contrast to what is observed in the underlying event. The presented measurements provide novel constraints on hadronisation and its Monte Carlo description. In particular, they demonstrate that the fragmentation of jets alone is insufficient to describe the strange and multi-strange particle production in hadronic collisions at LHC energies.