The measurements of the inclusive J/$\psi$ yield at midrapidity ($\left | y \right | < 0.9$) and forward rapidity (2.5 $< y <$ 4) in Pb$-$Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV with the ALICE detector at the LHC are reported. The inclusive J/$\psi$ production yields and nuclear modification factors, $R_{\rm AA}$, are measured as a function of the collision centrality, J/$\psi$ transverse momentum ($p_{\rm T}$), and rapidity. The J/$\psi$ average transverse momentum and squared transverse momentum ($\langle p_{\mathrm{T}}\rangle$ and $\langle p_{\mathrm{T}}^{\mathrm{2}}\rangle$) are evaluated as a function of the centrality at midrapidity. Compared to the previous ALICE publications, here the entire Pb$-$Pb collisions dataset collected during the LHC Run 2 is used, which improves the precision of the measurements and extends the $p_{\rm T}$ coverage. The $p_{\rm T}$-integrated $R_{\rm AA}$ shows a hint of an increasing trend towards unity from semicentral to central collisions at midrapidity, while it is flat at forward rapidity. The $p_{\rm T}$-differential $R_{\rm AA}$ shows a strong suppression at high $p_{\rm T}$ with less suppression at low $p_{\rm T}$ where it reaches a larger value at midrapidity compared to forward rapidity. The ratio of the $p_{\rm T}$-integrated yields of J/$\psi$ to those of D$^{0}$ mesons is reported for the first time for the central and semicentral event classes at midrapidity. Model calculations implementing charmonium production via the coalescence of charm quarks and antiquarks during the fireball evolution (transport models) or in a statistical approach with thermal weights are in good agreement with the data at low $p_{\rm T}$. At higher $p_{\rm T}$, the data are well described by transport models and a model based on energy loss in the strongly-interacting medium produced in nuclear collisions at the LHC.

The elliptic flow ($v_2$) of $\rm D^0$ mesons from beauty-hadron decays (non-prompt $\rm D^0$) was measured in midcentral (30-50%) Pb-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}}$ = 5.02 TeV with the ALICE detector at the LHC. The $\rm D^0$ mesons were reconstructed at midrapidity ($|y|<0.8$) from their hadronic decay $\mathrm{D^0 \to K^-\pi^+}$, in the transverse momentum interval $2 < p_{\rm T} < 12$ GeV/$c$. The result indicates a positive $v_2$ for non-prompt $\rm D^0$ mesons with a significance of 2.7$\sigma$. The non-prompt $\rm D^0$-meson $v_2$ is lower than that of prompt non-strange D mesons with 3.2$\sigma$ significance in $2 < p_{\rm T} < 8$ GeV/$c$, and compatible with the $v_2$ of beauty-decay electrons. Theoretical calculations of beauty-quark transport in a hydrodynamically expanding medium describe the measurement within uncertainties.

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 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.

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 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.

The $p_{\rm T}$-differential production cross sections of non-prompt ${\rm D^0}$, ${\rm D^+}$, and ${\rm D_s^+}$ mesons originating from beauty-hadron decays are measured in proton$-$proton collisions at a centre-of-mass energy $\sqrt{s}=13$ TeV. The measurements are performed at midrapidity, $|y| < 0.5$, with the data sample collected by ALICE from 2016 to 2018. The results are in agreement with predictions from several perturbative QCD calculations. The fragmentation fraction of beauty quarks to strange mesons divided by the one to non-strange mesons, $f_{\rm{s}}/(f_{\rm{u}} + f_{\rm{d}})$, is found to be $0.114 \pm 0.016~{\rm (stat.)} \pm 0.006~{\rm (syst.)} \pm 0.003~{\rm (BR)} \pm 0.003~{\rm (extrap.)}$. This value is compatible with previous measurements at lower centre-of-mass energies and in different collision systems in agreement with the assumption of universality of fragmentation functions. In addition, the dependence of the non-prompt D meson production on the centre-of-mass energy is investigated by comparing the results obtained at $\sqrt{s} = 5.02$ and 13 TeV, showing a hardening of the non-prompt D-meson $p_{\rm T}$-differential production cross section at higher $\sqrt{s}$. Finally, the ${\rm b\overline{b}}$ production cross section per unit of rapidity at midrapidity is calculated from the non-prompt ${\rm D^0}$, ${\rm D^+}$, ${\rm D_s^+}$, and $\Lambda_{\rm c}^+$ hadron measurements, obtaining ${\rm d}\sigma/{\rm d}y = 75.2\pm 3.2~(\mathrm{stat.}) \pm 5.2~(\mathrm{syst.})^{+12.3}_{-3.2} ~(\mathrm{extrap.})\text{ } \rm \mu b \;.$