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.

The azimuthal anisotropy of particles associated with jets (jet particles) at midrapidity is measured for the first time in p-Pb and Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV down to transverse momentum ($p_{\rm T}$) of 0.5 GeV/$c$ and 2 GeV/$c$, respectively, with ALICE. The results obtained in p-Pb collisions are based on a novel three-particle correlation technique. The azimuthal anisotropy coefficient $v_2$ in high-multiplicity p-Pb collisions is positive, with a significance reaching 6.8$\sigma$ at low $p_{\rm T}$, and its magnitude is smaller than in semicentral Pb-Pb collisions. In contrast to the measurements in Pb-Pb collisions, the $v_2$ coefficient is also found independent of $p_{\rm T}$ within uncertainties. Comparisons with the inclusive charged-particle $v_2$ and with AMPT calculations are discussed. The predictions suggest that parton interactions play an important role in generating a non-zero jet-particle $v_2$ in p-Pb collisions, even though they overestimate the reported measurement. These observations shed new insights on the understanding of the origin of the collective behaviour of jet particles in small systems such as p-Pb collisions, and provide significant stringent new constraints to models.

The production yield and angular anisotropy of prompt ${D_s^+}$ mesons were measured as a function of transverse momentum ($p_{ T}$) in Pb-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{ NN}} = 5.02$ TeV collected with the ALICE detector at the LHC. ${D_s^+}$ mesons and their charge conjugates were reconstructed at midrapidity ($|y|<0.5$) from their hadronic decay channel ${D_s^+ \to \phi \pi^+}$, with ${\phi \to K^-K^+}$, in the $p_{ T}$ intervals $2<p_{ T}<50$ GeV/$c$ and $2<p_{ T}<36$ GeV/$c$ for the 0-10% and 30-50% centrality intervals. For $p_{ T}>10$ GeV/$c$, the measured ${D_s^+}$-meson nuclear modification factor $R_{ AA}$ is consistent with the one of non-strange D mesons within uncertainties, while at lower $p_{ T}$ a hint for a ${D_s^+}$-meson $R_{ AA}$ larger than that of non-strange D mesons is seen. The enhanced production of ${D_s^+}$ relative to non-strange D mesons is also studied by comparing the $p_{ T}$-dependent ${D_s^+/D^0}$ production yield ratios in Pb-Pb and in pp collisions. The ratio measured in Pb-Pb collisions is found to be on average higher than that in pp collisions in the interval $2<p_{ T} <8$ GeV/$c$ with a significance of 2.3$\sigma$ and 2.4$\sigma$ for the 0-10% and 30-50% centrality intervals. The azimuthal anisotropy coefficient $v_2$ of prompt ${D_s^+}$ mesons was measured in Pb-Pb collisions in the 30-50% centrality interval and is found to be compatible with that of non-strange D mesons. The main features of the measured $R_{ AA}$, ${D_s^+/D^0}$ ratio, and $v_2$ as a function of $p_{ T}$ are described by theoretical calculations of charm-quark transport in a hydrodynamically expanding quark-gluon plasma including hadronisation via charm-quark recombination with light quarks from the medium. The $p_{ T}$-integrated production yield of ${D_s^+}$ mesons is compatible with the prediction of the statistical hadronisation model.

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.

Measurements of charged-particle production in pp, p$-$Pb, and Pb$-$Pb collisions in the toward, away, and transverse regions with the ALICE detector are discussed. These regions are defined event-by-event relative to the azimuthal direction of the charged trigger particle, which is the reconstructed particle with the largest transverse momentum ($p_{\mathrm{T}}^{\rm trig}$) in the range $8<p_{\mathrm{T}}^{\rm trig}<15$ GeV$/c$. The toward and away regions contain the primary and recoil jets, respectively; both regions are accompanied by the underlying event (UE). In contrast, the transverse region perpendicular to the direction of the trigger particle is dominated by the so-called UE dynamics, and includes also contributions from initial- and final-state radiation. The relative transverse activity classifier, $R_{\mathrm{T}}=N_{\mathrm{ch}}^{\mathrm{T}}/\langle N_{\mathrm{ch}}^{\mathrm{T}}\rangle$, is used to group events according to their UE activity, where $N_{\mathrm{ch}}^{\mathrm{T}}$ is the charged-particle multiplicity per event in the transverse region and $\langle N_{\mathrm{ch}}^{\mathrm{T}}\rangle$ is the mean value over the whole analysed sample. The energy dependence of the $R_{\mathrm{T}}$ distributions in pp collisions at $\sqrt{s}=2.76$, 5.02, 7, and 13 TeV is reported, exploring the Koba-Nielsen-Olesen (KNO) scaling properties of the multiplicity distributions. The first measurements of charged-particle $p_{\rm T}$ spectra as a function of $R_{\mathrm{T}}$ in the three azimuthal regions in pp, p$-$Pb, and Pb$-$Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV are also reported. Data are compared with predictions obtained from the event generators PYTHIA 8 and EPOS LHC. This set of measurements is expected to contribute to the understanding of the origin of collective-like effects in small collision systems (pp and p$-$Pb).

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.

The ALICE Collaboration reports measurements of the large relative transverse momentum ($k_{\text{T}}$) component of jet substructure in pp and Pb$-$Pb collisions at center-of-mass energy per nucleon pair $\sqrt{s_{\text{NN}}}=5.02$ TeV. Enhancement in the yield of such large-$k_{\text{T}}$ emissions in head-on Pb$-$Pb collisions is predicted to arise from partonic scattering with quasi-particles of the quark-gluon plasma. The analysis utilizes charged-particle jets reconstructed by the anti-$k_{\text{T}}$ algorithm with resolution parameter $R=0.2$ in the transverse-momentum interval $60 < p_{\text{T,ch jet}} < 80$ GeV/$c$. The soft drop and dynamical grooming algorithms are used to identify high transverse momentum splittings in the jet shower. Comparison of measurements in Pb$-$Pb and pp collisions shows medium-induced narrowing, corresponding to yield suppression of high-$k_{\text{T}}$ splittings, in contrast to the expectation of yield enhancement due to quasi-particle scattering. The measurements are compared to theoretical model calculations incorporating jet modification due to jet-medium interactions ("jet quenching"), both with and without quasi-particle scattering effects. These measurements provide new insight into the underlying mechanisms and theoretical modeling of jet quenching.

The ALICE Collaboration at the CERN LHC has measured the inclusive production cross section of isolated photons at midrapidity as a function of the photon transverse momentum ($p_{\rm T}^{\gamma}$), in Pb$-$Pb collisions in different centrality intervals, and in pp collisions, at centre-of-momentum energy per nucleon pair of $\sqrt{s_{\rm NN}}=5.02$ TeV. The photon transverse momentum range is between 10-14 and 40-140 GeV/$c$, depending on the collision system and on the Pb$-$Pb centrality class. The result extends to lower $p_{\rm T}^{\gamma}$ than previously published results by the ATLAS and CMS experiments at the same collision energy. The covered pseudorapidity range is $|\eta^{\gamma}| <0.67$. The isolation selection is based on a charged particle isolation momentum threshold $p_{\rm T}^{\rm iso,~ch} = 1.5$ GeV/$c$ within a cone of radii $R=0.2$ and $0.4$. The nuclear modification factor is calculated and found to be consistent with unity in all centrality classes, and also consistent with the HG-PYTHIA model, which describes the event selection and geometry biases that affect the centrality determination in peripheral Pb$-$Pb collisions. The measurement is compared to next-to-leading order perturbative QCD calculations and to the measurements of isolated photons and Z$^0$ bosons from the CMS experiment, which are all found to be in agreement.

ALICE is a large experiment at the CERN Large Hadron Collider. Located 52 meters underground, its detectors are suitable to measure muons produced by cosmic-ray interactions in the atmosphere. In this paper, the studies of the cosmic muons registered by ALICE during Run 2 (2015--2018) are described. The analysis is limited to multimuon events defined as events with more than four detected muons ($N_\mu>4$) and in the zenith angle range $0^{\circ}<\theta<50^{\circ}$. The results are compared with Monte Carlo simulations using three of the main hadronic interaction models describing the air shower development in the atmosphere: QGSJET-II-04, EPOS-LHC, and SIBYLL 2.3d. The interval of the primary cosmic-ray energy involved in the measured muon multiplicity distribution is about $ 4 \times 10^{15}<E_\mathrm{prim}< 6 \times 10^{16}$~eV. In this interval none of the three models is able to describe precisely the trend of the composition of cosmic rays as the energy increases. However, QGSJET-II-04 is found to be the only model capable of reproducing reasonably well the muon multiplicity distribution, assuming a heavy composition of the primary cosmic rays over the whole energy range, while SIBYLL 2.3d and EPOS-LHC underpredict the number of muons in a large interval of multiplicity by more than $20\%$ and $30\%$, respectively. The rate of high muon multiplicity events ($N_\mu>100$) obtained with QGSJET-II-04 and SIBYLL 2.3d is compatible with the data, while EPOS-LHC produces a significantly lower rate ($55\%$ of the measured rate). For both QGSJET-II-04 and SIBYLL 2.3d, the rate is close to the data when the composition is assumed to be dominated by heavy elements, an outcome compatible with the average energy $E_\mathrm{prim} \sim 10^{17}$ eV of these events. This result places significant constraints on more exotic production mechanisms.

The production yields of the orbitally excited charm-strange mesons $\mathrm{D_{s1}(1^{+})(2536)^+}$ and $\mathrm{D_{s2}^{*}(2^{+})(2573)^+}$were measured for the first time in proton-proton (pp) collisions at a center-of-mass energy of $\sqrt{s}=13$ TeV with the ALICE experiment at the LHC. The $\mathrm{D_{s1}^+}$ and $\mathrm{D_{s2}^{*+}}$ mesons were measured at midrapidity ($|y|<0.5$) in minimum-bias and high-multiplicity pp collisions in the transverse-momentum interval $2<p_{\rm T}<24$ GeV/$c$. Their production yields relative to the $\mathrm{D_{s}^{+}}$ ground-state yield were found to be compatible between minimum-bias and high-multiplicity collisions, as well as with previous measurements in ${\rm e}^{\pm}$p and ${\rm e^+ e^-}$ collisions. The measured $\mathrm{D_{s1}^+/D_{s}^{+}}$ and $\mathrm{D_{s2}^{*+}/D_{s}^{+}}$ yield ratios are described by statistical hadronization models and can be used to tune the parameters governing the production of excited charm-strange hadrons in Monte Carlo generators, such as PYTHIA 8.