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 \;.$

The first measurement of $_{\Lambda}^{3}\mathrm{H}$ and $^3_ {\overline{\Lambda}}\overline{\mathrm{H}}$ differential production with respect to transverse momentum and centrality in Pb$-$Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02$~TeV is presented. The $_{\Lambda}^{3}\mathrm{H}$ has been reconstructed via its two-charged-body decay channel, i.e., $_{\Lambda}^{3}\mathrm{H} \rightarrow {}^{3}\mathrm{He} + \pi^{-}$. A Blast-Wave model fit of the $p_{\rm T}$-differential spectra of all nuclear species measured by the ALICE collaboration suggests that the $_{\Lambda}^{3}\mathrm{H}$ kinetic freeze-out surface is consistent with that of other nuclei. The ratio between the integrated yields of $_{\Lambda}^{3}\mathrm{H}$ and $^3\mathrm{He}$ is compared to predictions from the statistical hadronisation model and the coalescence model, with the latter being favoured by the presented measurements.

The production cross sections of ${\rm D^0}$ and $\Lambda^+_{\rm c}$ hadrons originating from beauty-hadron decays (i.e. non-prompt) were measured for the first time at midrapidity ($|y|<0.5$) by the ALICE Collaboration in proton-proton collisions at a center-of-mass energy $\sqrt{s}=13$ TeV. They are described within uncertainties by perturbative QCD calculations employing the fragmentation fractions of beauty quarks to baryons measured at forward rapidity by the LHCb Collaboration. The ${\rm b\overline{b}}$ production cross section per unit of rapidity at midrapidity, estimated from these measurements, is ${\rm d}\sigma_{\rm b\overline{b}}/{\rm d}y|_{|y|<0.5} = 83.1 \pm 3.5 (\mathrm{stat.}) \pm 5.4(\mathrm{syst.}) ^{+12.3}_{-3.2} (\mathrm{extrap.})\,\mu$b. The baryon-to-meson ratios are computed to investigate the hadronization mechanism of beauty quarks. The non-prompt $\Lambda^+_{\rm c}/{\rm D^0}$ production ratio has a similar trend to the one measured for the promptly produced charmed particles and to the p$/\pi^+$ and $\Lambda/{\rm K^0_S}$ ratios, suggesting a similar baryon-formation mechanism among light, strange, charm, and beauty hadrons. The $p_{\rm T}$-integrated non-prompt $\Lambda_{\rm c}/{\rm D^0}$ ratio is found to be significantly higher than the one measured in e$^+$e$^-$ collisions.

We report on the measurement of inclusive, non-prompt, and prompt J/$ψ$-hadron correlations by the ALICE Collaboration at the CERN Large Hadron Collider in pp collisions at a center-of-mass energy of $13~\text{TeV}$. The correlations are studied at midrapidity ($|y| < 0.9$) in the transverse momentum ranges $p_{\rm T} < 40~\text{GeV}/c$ for the J/$ψ$ and $0.15 < p_{\rm T} < 10~\text{GeV}/c$ and $|η|<0.9$ for the associated hadrons. The measurement is based on minimum bias and high multiplicity data samples corresponding to integrated luminosities of $L_{\text{int}} = 34~\text{nb}^{-1}$ and $L_{\text{int}} = 6.9~\text{pb}^{-1}$, respectively. In addition, two more data samples are employed, requiring, on top of the minimum bias condition, a threshold on the tower energy of $E = 4$ and $9~\text{GeV}$ in the ALICE electromagnetic calorimeters, which correspond to integrated luminosities of $L_{\text{int}} = 0.9~\text{pb}^{-1}$ and $L_{\text{int}} = 8.4~\text{pb}^{-1}$, respectively. The azimuthally integrated near and away side yields of associated charged hadrons per J/$ψ$ trigger are presented as a function of the J/$ψ$ and associated hadron transverse momentum. The measurements are discussed in comparison to PYTHIA calculations.

In this letter, measurements of (anti)alpha production in central (0$-$10%) Pb$-$Pb collisions at a center-of-mass energy per nucleon$-$nucleon pair of $\sqrt{s_{\rm NN}}$ = 5.02 TeV are presented, including the first measurement of an antialpha transverse-momentum spectrum. Owing to its large mass, the production of (anti)alpha is expected to be sensitive to different particle production models. The production yields and transverse-momentum spectra of nuclei are of particular interest because they provide a stringent test of these models. The averaged antialpha and alpha spectrum is compared to the spectra of lighter particles, by including it into a common blast-wave fit capturing the hydrodynamic-like flow of all particles. This fit is indicating that the (anti)alpha also participates in the collective expansion of the medium created in the collision. A blast-wave fit including only protons, (anti)alpha, and other light nuclei results in a similar flow velocity as the fit that includes all particles. A similar flow velocity, but a significantly larger kinetic freeze-out temperature is obtained when only protons and light nuclei are included in the fit. The coalescence parameter $B_4$ is well described by calculations from a statistical hadronization model but significantly underestimated by calculations assuming nucleus formation via coalescence of nucleons. Similarly, the (anti)alpha-to-proton ratio is well described by the statistical hadronization model. On the other hand, coalescence calculations including approaches with different implementations of the (anti)alpha substructure tend to underestimate the data.

The ALICE Collaboration reports a search for jet quenching effects in high-multiplicity (HM) proton$-$proton collisions at $\sqrt{s}$ = 13 TeV, using the semi-inclusive azimuthal-difference distribution $\Delta\varphi$ of charged-particle jets recoiling from a high transverse momentum (high-$p_{\mathrm{T,trig}}$) trigger hadron. Jet quenching may broaden the $\Delta\varphi$ distribution measured in HM events compared to that in minimum bias (MB) events. The measurement employs a $p_{\mathrm{T,trig}}$-differential observable for data-driven suppression of the contribution of multiple partonic interactions, which is the dominant background. While azimuthal broadening is indeed observed in HM compared to MB events, similar broadening for HM events is observed for simulations based on the PYTHIA 8 Monte Carlo generator, which does not incorporate jet quenching. Detailed analysis of these data and simulations show that the azimuthal broadening is due to bias of the HM selection towards events with multiple jets in the final state. The identification of this bias has implications for all jet quenching searches where selection is made on the event activity.

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 multiplicity dependence of jet production in pp collisions at the centre-of-mass energy of $\sqrt{s} = 13\ \mathrm{TeV}$ is studied for the first time. Jets are reconstructed from charged particles using the anti-$k_\mathrm{T}$ algorithm with resolution parameters $R$ varying from $0.2$ to $0.7$. The jets are measured in the pseudorapidity range $|\eta_{\rm jet}|< 0.9-R$ and in the transverse momentum range $5<p_\mathrm{T,jet}^{\rm ch}<140\ \mathrm{GeV}/c$. The multiplicity intervals are categorised by the ALICE forward detector V0. The $p_{\mathrm{T}}$ differential cross section of charged-particle jets are compared to leading order (LO) and next-to-leading order (NLO) perturbative quantum chromodynamics (pQCD) calculations. It is found that the data are better described by the NLO calculation, although the NLO prediction overestimates the jet cross section below $20\ \mathrm{GeV}/c$. The cross section ratios for different $R$ are also measured and compared to model calculations. These measurements provide insights into the angular dependence of jet fragmentation. The jet yield increases with increasing self-normalised charged-particle multiplicity. This increase shows only a weak dependence on jet transverse momentum and resolution parameter at the highest multiplicity. While such behaviour is qualitatively described by the present version of PYTHIA, quantitative description may require implementing new mechanisms for multi-particle production in hadronic collisions.