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.
Two-particle transverse momentum differential correlators, recently measured in Pb--Pb collisions at energies available at the CERN Large Hadron Collider (LHC), provide an additional tool to gain insights into particle production mechanisms and infer transport properties, such as the ratio of shear viscosity to entropy density, of the medium created in Pb-Pb collisions. The longitudinal long-range correlations and the large azimuthal anisotropy measured at low transverse momenta in small collision systems, namely pp and p-Pb, at LHC energies resemble manifestations of collective behaviour. This suggests that locally equilibrated matter may be produced in these small collision systems, similar to what is observed in Pb-Pb collisions. In this work, the same two-particle transverse momentum differential correlators are exploited in pp and p-Pb collisions at $\sqrt{s} = 7$ TeV and $\sqrt{s_{\rm NN}} = 5.02$ TeV, respectively, to seek evidence for viscous effects. Specifically, the strength and shape of the correlators are studied as a function of the produced particle multiplicity to identify evidence for longitudinal broadening that might reveal the presence of viscous effects in these smaller systems. The measured correlators and their evolution from pp and p--Pb to Pb--Pb collisions are additionally compared to predictions from Monte Carlo event generators, and the potential presence of viscous effects is discussed.
Event-by-event fluctuations of the event-wise mean transverse momentum, $\langle p_{\mathrm{T}}\rangle$, of charged particles produced in proton-proton (pp) collisions at $\sqrt{s}$ = 5.02 TeV, Xe-Xe collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.44 TeV, and Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV are studied using the ALICE detector based on the integral correlator $\langle\langle \Delta p_{\rm T}\Delta p_{\rm T}\rangle\rangle $. The correlator strength is found to decrease monotonically with increasing produced charged-particle multiplicity measured at midrapidity in all three systems. In Xe-Xe and Pb-Pb collisions, the multiplicity dependence of the correlator deviates significantly from a simple power-law scaling as well as from the predictions of the HIJING and AMPT models. The observed deviation from power-law scaling is expected from transverse radial flow in semicentral to central Xe-Xe and Pb-Pb collisions. In pp collisions, the correlation strength is also studied by classifying the events based on the transverse spherocity, $S_0$, of the particle production at midrapidity, used as a proxy for the presence of a pronounced back-to-back jet topology. Low-spherocity (jetty) events feature a larger correlation strength than those with high spherocity (isotropic). The strength and multiplicity dependence of jetty and isotropic events are well reproduced by calculations with the PYTHIA 8 and EPOS LHC models.
Two-particle azimuthal correlations are measured with the ALICE apparatus in pp collisions at $\sqrt{s} = 13$ TeV to explore strangeness- and multiplicity-related effects in the fragmentation of jets and the transition regime between bulk and hard production, probed with the condition that a strange meson (K$^{0}_{\rm S}$) or baryon ($\Lambda$) with transverse momentum $p_{\rm T} > 3$ GeV/c is produced. Azimuthal correlations between kaons or $\Lambda$ hyperons with other hadrons are presented at midrapidity for a broad range of the trigger ($3 < p_{\rm T}^{\rm trigg} < 20$ GeV/$c$) and associated particle $p_{\rm T}$ (1 GeV/$c$$< p_{\rm T}^{\rm assoc} < p_{\rm T}^{\rm trigg}$), for minimum-bias events and as a function of the event multiplicity. The near- and away-side peak yields are compared for the case of either K$^{0}_{\rm S}$ or $\Lambda$($\overline{\Lambda}$) being the trigger particle with that of inclusive hadrons (a sample dominated by pions). In addition, the measurements are compared with predictions from PYTHIA 8 and EPOS LHC event generators.
Recent measurements of charm-baryon production in hadronic collisions have questioned the universality of charm-quark fragmentation across different collision systems. In this work the fragmentation of charm quarks into charm baryons is probed, by presenting the first measurement of the longitudinal jet momentum fraction carried by $\Lambda_{\rm c}^{+}$ baryons, $z^{\mathrm{ch}}_\mathrm{||}$, in hadronic collisions. The results are obtained in proton-proton (pp) collisions at $\sqrt{s}=13$ TeV at the LHC, with $\Lambda_{\rm c}^{+}$ baryons and charged (track-based) jets reconstructed in the transverse momentum intervals of $3 \leq p_{\rm T}^{\Lambda_{\rm c}^{+}} < 15$ GeV/$c$ and $7 \leq p_{\rm T}^{\rm jet\;ch} < 15$ GeV/$c$, respectively. The $z^{\mathrm{ch}}_\mathrm{||}$ distribution is compared to a measurement of ${\rm D}^0$-tagged charged jets in pp collisions as well as to PYTHIA 8 simulations. The data hints that the fragmentation of charm quarks into charm baryons is softer with respect to charm mesons, in the measured kinematic interval, as predicted by hadronisation models which include colour correlations beyond leading-colour in the string formation.
The interactions of kaons (K) and antikaons ($\mathrm{\overline{K}}$) with few nucleons (N) were studied so far using kaonic atom data and measurements of kaon production and interaction yields in nuclei. Some details of the three-body KNN and $\mathrm{\overline{K}}$NN dynamics are still not well understood, mainly due to the overlap with multi-nucleon interactions in nuclei. An alternative method to probe the dynamics of three-body systems with kaons is to study the final state interaction within triplet of particles emitted in pp collisions at the Large Hadron Collider, which are free from effects due to the presence of bound nucleons. This Letter reports the first femtoscopic study of p$-$p$-$K$^+$ and p$-$p$-$K$^-$ correlations measured in high-multiplicity pp collisions at $\sqrt{s}$ = 13 TeV by the ALICE Collaboration. The analysis shows that the measured p$-$p$-$K$^+$ and p$-$p$-$K$^-$ correlation functions can be interpreted in terms of pairwise interactions in the triplets, indicating that the dynamics of such systems is dominated by the two-body interactions without significant contributions from three-body effects or bound states.
The production of electrons from beauty-hadron decays was measured at midrapidity in proton-proton (pp) and central Pb-Pb collisions at center-of-mass energy per nucleon-nucleon pair $\sqrt{s_{\rm NN}}$ = 5.02 TeV, using the ALICE detector at the LHC. The cross section measured in pp collisions in the transverse momentum interval $2 < p_{\rm T} < 8$ GeV/$c$ was compared with models based on perturbative quantum chromodynamics calculations. The yield in the 10% most central Pb-Pb collisions, measured in the interval $2 < p_{\rm T} < 26$ GeV/$c$, was used to compute the nuclear modification factor $R_{\rm AA}$, extrapolating the pp reference cross section to $p_{\rm T}$ larger than 8 GeV/$c$. The measured $R_{\rm AA}$ shows significant suppression of the yield of electrons from beauty-hadron decays at high $p_{\rm T}$ and does not show a significant dependence on $p_{\rm T}$ above 8 GeV/$c$ within uncertainties. The results are described by several theoretical models based on different implementations of the interaction of heavy quarks with a quark-gluon plasma, which predict a smaller energy loss for beauty quarks compared to light and charm quarks.
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 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 interaction of $\rm{K}^{-}$ with protons is characterised by the presence of several coupled channels, systems like ${\rm \overline{K}^0}$n and $\pi\Sigma$ with a similar mass and the same quantum numbers as the $\rm{K}^{-}$p state. The strengths of these couplings to the $\rm{K}^{-}$p system are of crucial importance for the understanding of the nature of the $\Lambda(1405)$ resonance and of the attractive $\rm{K}^{-}$p strong interaction. In this article, we present measurements of the $\rm{K}^{-}$p correlation functions in relative momentum space obtained in pp collisions at $\sqrt{s}~=~13$ TeV, in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~5.02$ TeV, and (semi)peripheral Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~5.02$ TeV. The emitting source size, composed of a core radius anchored to the $\rm{K}^{+}$p correlation and of a resonance halo specific to each particle pair, varies between 1 and 2 fm in these collision systems. The strength and the effects of the ${\rm \overline{K}^0}$n and $\pi\Sigma$ inelastic channels on the measured $\rm{K}^{-}$p correlation function are investigated in the different colliding systems by comparing the data with state-of-the-art models of chiral potentials. A novel approach to determine the conversion weights $\omega$, necessary to quantify the amount of produced inelastic channels in the correlation function, is presented. In this method, particle yields are estimated from thermal model predictions, and their kinematic distribution from blast-wave fits to measured data. The comparison of chiral potentials to the measured $\rm{K}^{-}$p interaction indicates that, while the $\pi\Sigma-\rm{K}^{-}$p dynamics is well reproduced by the model, the coupling to the ${\rm \overline{K}^0}$n channel in the model is currently underestimated.
K$^+$p (K$^+$p $\oplus$ K$^-\overline{\mathrm p}$) correlation function in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV (60-70%).