The study of (anti-)deuteron production in pp collisions has proven to be a powerful tool to investigate the formation mechanism of loosely bound states in high energy hadronic collisions. In this paper the production of (anti-)deuterons is studied as a function of the charged particle multiplicity in inelastic pp collisions at $\sqrt{s}=13$ TeV using the ALICE experiment. Thanks to the large number of accumulated minimum bias events, it has been possible to measure (anti-)deuteron production in pp collisions up to the same charged particle multiplicity ($\rm{d} N_{ch}/\rm{d}\eta\sim26$) as measured in p-Pb collisions at similar centre-of-mass energies. Within the uncertainties, the deuteron yield in pp collisions resembles the one in p-Pb interactions, suggesting a common formation mechanism behind the production of light nuclei in hadronic interactions. In this context the measurements are compared with the expectations of coalescence and Statistical Hadronisation Models (SHM).

This article reports measurements characterizing the Underlying Event (UE) associated with hard scatterings at midrapidity in pp collisions at $\sqrt{s}=13$ TeV. The hard scatterings are identified by the leading particle, the charged particle with the highest transverse momentum ($p_{\rm T}^{\rm leading}$) in the event. Charged-particle number and summed transverse-momentum densities are measured in different azimuthal regions defined with respect to the leading particle direction: Toward, Transverse, and Away. The Toward and Away regions contain the fragmentation products of the hard scatterings in addition to the UE contribution, whereas particles in the Transverse region are expected to originate predominantly from the UE. The study is performed as a function of $p_{\rm T}^{\rm leading}$ with three different $p_{\rm T}$ thresholds for the associated particles, $p_{\rm T}^{\rm min} >$ 0.15, 0.5, and 1.0 GeV/$c$. The charged-particle density in the Transverse region rises steeply for low values of $p_{\rm T}^{\rm leading}$ and reaches a plateau. The results confirm the trend that the charged-particle density in the Transverse region shows a stronger increase with $\sqrt{s}$ than the inclusive charged-particle density at midrapidity. The UE activity is increased by approximately 20% when going from 7 to 13 TeV. The plateau in the Transverse region ($5 < p_{\rm T}^{\rm leading} < ~ 40$ GeV/$c$ ) is further characterized by the probability distribution of its charged-particle multiplicity normalized to its average value (relative transverse activity, $R_{T}$) and the mean transverse momentum as a function of $R_{T}$. Experimental results are compared to model calculations using PYTHIA 8 and EPOS LHC. The overall agreement between models and data is within 30%. These measurements provide new insights on the interplay between hard scatterings and the associated UE in pp collisions.

This Letter presents the first direct investigation of the p-$\Sigma^{0}$ interaction, using the femtoscopy technique in high-multiplicity pp collisions at $\sqrt{s}$ = 13 TeV measured by the ALICE detector. The $\Sigma^{0}$ is reconstructed via the decay channel to $\Lambda \gamma$, and the subsequent decay of $\Lambda$ to p$\pi^-$. The photon is detected via the conversion in material to e$^{+}$e$^{-}$ pairs exploiting the unique capability of the ALICE detector to measure electrons at low transverse momenta. The measured p-$\Sigma^{0}$ correlation indicates a shallow strong interaction. The comparison of the data to several theoretical predictions obtained employing the $Correlation~Analysis~Tool~using~the~Schr\"odinger~Equation$ (CATS) and the Lednick\'y-Lyuboshits approach shows that the current experimental precision does not yet allow to discriminate between different models, as it is the case for the available scattering and hypernuclei data. Nevertheless, the p-$\Sigma^{0}$ correlation function is found to be sensitive to the strong interaction, and driven by the interplay of the different spin and isospin channels. This pioneering study demonstrates the feasibility of a femtoscopic measurement in the p-$\Sigma^{0}$ channel and with the expected larger data samples in LHC Run 3 and Run 4, the p-$\Sigma^{0}$ interaction will be constrained with high precision.

Deuteron production in high-energy collisions is sensitive to the space-time evolution of the collision system, and is typically described by a coalescence mechanism. For the first time, we present results on jet-associated deuteron production in \pp collisions at $\sqrt{s}=13$ TeV, providing an opportunity to test the established picture for deuteron production in events with a hard scattering. Using a trigger particle with high transverse-momentum ($p_{\rm T}>5$ GeV/$c$) as a proxy for the presence of a jet at midrapidity, we observe a measurable population of deuterons being produced around the jet proxy. The associated deuteron yield measured in a narrow angular range around the trigger particle differs by 2.4-4.8 standard deviations from the uncorrelated background. The data are described by PYTHIA model calculations featuring baryon coalescence.

Understanding the production mechanism of light (anti)nuclei is one of the key challenges of nuclear physics and has important consequences for astrophysics, since it provides an input for indirect dark-matter searches in space. In this paper, the latest results about the production of light (anti)nuclei in pp collisions at $\sqrt{s} = 13$ TeV are presented, focusing on the comparison with the predictions of coalescence and thermal models. For the first time, the coalescence parameters $B_2$ for deuterons and $B_3$ for helions are compared with parameter-free theoretical predictions that are directly constrained by the femtoscopic measurement of the source radius in the same event class. A fair description of the data with a Gaussian wave function is observed for both deuteron and helion, supporting the coalescence mechanism for the production of light (anti)nuclei in pp collisions. This method paves the way for future investigations of the internal structure of more complex nuclear clusters, including the hypertriton.

The cross section of jets reconstructed from charged particles is measured in the transverse momentum range of $5<p_\mathrm{T}<100\ \mathrm{GeV}/c$ in pp collisions at the center-of-mass energy of $\sqrt{s} = 5.02\ \mathrm{TeV}$ with the ALICE detector. The jets are reconstructed using the anti-$k_\mathrm{T}$ algorithm with resolution parameters $R=0.2$, $0.3$, $0.4$, and $0.6$ in the pseudorapidity range $|\eta|< 0.9-R$. The charged jet cross sections are compared with the leading order (LO) and to next-to-leading order (NLO) perturbative Quantum ChromoDynamics (pQCD) calculations. It was found that the NLO calculations agree better with the measurements. The cross section ratios for different resolution parameters were also measured. These ratios increase from low $p_\mathrm{T}$ to high $p_\mathrm{T}$ and saturate at high $p_\mathrm{T}$, indicating that jet collimation is larger at high $p_\mathrm{T}$ than at low $p_\mathrm{T}$. These results provide a precision test of pQCD predictions and serve as a baseline for the measurement in Pb$-$Pb collisions at the same energy to quantify the effects of the hot and dense medium created in heavy-ion collisions at the LHC.

The first experimental information on the strong interaction between $\Lambda$ and $\Xi^-$ strange baryons is presented in this Letter. The correlation function of $\Lambda-\Xi^-$ and $\overline{\Lambda}-\overline{\Xi}^{+}$ pairs produced in high-multiplicity proton-proton (pp) collisions at $\sqrt{s}$ = 13 TeV at the LHC is measured as a function of the relative momentum of the pair. The femtoscopy method is used to calculate the correlation function, which is then compared with theoretical expectations obtained using a meson exchange model, chiral effective field theory, and Lattice QCD calculations close to the physical point. Data support predictions of small scattering parameters while discarding versions with large ones, thus suggesting a weak $\Lambda-\Xi^{-}$ interaction. The limited statistical significance of the data does not yet allow one to constrain the effects of coupled channels like $\Sigma-\Xi$ and N$-\Omega$.

We report on the inclusive J/$\psi$ production cross section measured at the CERN Large Hadron Collider in proton-proton collisions at a centre-of-mass energy $\sqrt{s}~=~13$ TeV. The J/$\psi$ mesons are reconstructed in the $\rm e^{+} e^{-}$ decay channel and the measurements are performed at midrapidity ($|y|<0.9$) in the transverse-momentum interval $0< p_{\rm T} <40$ GeV/$c$, using a minimum bias data sample corresponding to an integrated luminosity $L_{\text{int}} = 32.2~\text{nb}^{-1}$ and an Electromagnetic Calorimeter triggered data sample with $L_{\text{int}} = 8.3~\mathrm{pb}^{-1}$. The $p_{\rm T}$-integrated J/$\psi$ production cross section at midrapidity, computed using the minimum bias data sample, is $\text{d}\sigma/\text{d}y|_{y=0} = 8.97\pm0.24~(\text{stat})\pm0.48~(\text{syst})\pm0.15~(\text{lumi})~\mu\text{b}$. An approximate logarithmic dependence with the collision energy is suggested by these results and available world data, in agreement with model predictions. The integrated and $p_{\rm T}$-differential measurements are compared with measurements in pp collisions at lower energies and with several recent phenomenological calculations based on the non-relativistic QCD and Color Evaporation models.

The study of the strength and behaviour of the antikaon-nucleon ($\mathrm{\overline{K}N}$) interaction constitutes one of the key focuses of the strangeness sector in low-energy Quantum Chromodynamics (QCD). In this letter a unique high-precision measurement of the strong interaction between kaons and protons, close and above the kinematic threshold, is presented. The femtoscopic measurements of the correlation function at low pair-frame relative momentum of (K$^+$ p $\oplus$ K$^-$ $\overline{\mathrm{p}}$) and (K$^-$ p $\oplus$ K$^+$ $\overline{\mathrm{p}}$) pairs measured in pp collisions at $\sqrt{s}$ = 5, 7 and 13 TeV are reported. A structure observed around a relative momentum of 58 MeV/$c$ in the measured correlation function of (K$^-$ p $\oplus$ K$^+$ $\overline{\mathrm{p}}$) with a significance of 4.4. $\sigma$ constitutes the first experimental evidence for the opening of the $(\mathrm{\overline{K}^0 n} \oplus \mathrm{K^0 \overline{n}})$ isospin breaking channel due to the mass difference between charged and neutral kaons. The measured correlation functions have been compared to J\"{u}lich and Kyoto models in addition to the Coulomb potential. The high-precision data at low relative momenta presented in this work prove femtoscopy to be a powerful complementary tool to scattering experiments and provide new constraints above the $\mathrm{\overline{K}N}$ threshold for low-energy QCD chiral models.

The inclusive production of the charm-strange baryon $\rm \Omega_c^0$ is measured for the first time via its hadronic decay into $\rm \Omega^-\pi^+$ at midrapidity ($|y|<0.5$) in proton-proton (pp) collisions at the centre-of-mass energy $\sqrt{s}=13$ TeV with the ALICE detector at the LHC. The transverse momentum ($p_{\rm T}$) differential cross section multiplied by the branching ratio is presented in the interval $2<p_{\rm T}<12~{\rm GeV}/c$. The $p_{\rm T}$ dependence of the $\rm \Omega_c^0$-baryon production relative to the prompt $\rm D^0$-meson and to the prompt $\rm \Xi_c^0$-baryon production is compared to various models that take different hadronisation mechanisms into consideration. In the measured $p_{\rm T}$ interval, the ratio of the $p_{\rm T}$-integrated cross sections of $\rm \Omega_c^0$ and prompt $\Lambda_{\rm c}^{+}$ baryons multiplied by the $\rm \Omega^-\pi^+$ branching ratio is found to be larger by a factor of about 20 with a significance of about 4$\sigma$ when compared to $\rm e^+e^-$ collisions.

The production of $\phi$ mesons has been studied in pp collisions at LHC energies with the ALICE detector via the dimuon decay channel in the rapidity region $2.5 < y < 4$. Measurements of the differential cross section ${\rm d}^2\sigma/{\rm d}y {\rm d}p_{\rm T}$ are presented as a function of the transverse momentum ($p_{\rm T}$) at the center-of-mass energies $\sqrt{s}=5.02$, 8 and 13 TeV and compared with the ALICE results at midrapidity. The differential cross sections at $\sqrt{s}=5.02$ and 13 TeV are also studied in several rapidity intervals as a function of $p_{\rm T}$, and as a function of rapidity in three $p_{\rm T}$ intervals. A hardening of the $p_{\rm T}$-differential cross section with the collision energy is observed, while, for a given energy, $p_{\rm T}$ spectra soften with increasing rapidity and, conversely, rapidity distributions get slightly narrower at increasing $p_{\rm T}$. The new results, complementing the published measurements at $\sqrt{s}=2.76$ and 7 TeV, allow one to establish the energy dependence of $\phi$ meson production and to compare the measured cross sections with phenomenological models. None of the considered models manages to describe the evolution of the cross section with $p_{\rm T}$ and rapidity at all the energies.

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.

Measurements of the inclusive J/$\psi$ yield as a function of charged-particle pseudorapidity density ${\rm d}N_{\rm ch}/\rm{d}\eta$ in pp collisions at $\sqrt{s}$ = 13 TeV with ALICE at the LHC are reported. The J/$\psi$ meson yield is measured at midrapidity ($|y|<0.9$) in the dielectron channel, for events selected based on the charged-particle multiplicity at midrapidity ($|\eta|<1$) and at forward rapidity ($-3.7<\eta<-1.7$ and $2.8<\eta<5.1$); both observables are normalized to their corresponding averages in minimum bias events. The increase of the normalized J/$\psi$ yield with normalized ${\rm d}N_{\rm ch}/\rm{d}\eta$ is significantly stronger than linear and dependent on the transverse momentum. The data are compared to theoretical predictions, which describe the observed trends well, albeit not always quantitatively.

The production rates and the transverse momentum distribution of strange hadrons at mid-rapidity ($\ |y\ | < 0.5$) are measured in proton-proton collisions at $\sqrt{s}$ = 13 TeV as a function of the charged particle multiplicity, using the ALICE detector at the LHC. The production rates of $\rm{K}^{0}_{S}$, $\Lambda$, $\Xi$, and $\Omega$ increase with the multiplicity faster than what is reported for inclusive charged particles. The increase is found to be more pronounced for hadrons with a larger strangeness content. Possible auto-correlations between the charged particles and the strange hadrons are evaluated by measuring the event-activity with charged particle multiplicity estimators covering different pseudorapidity regions. When comparing to lower energy results, the yields of strange hadrons are found to depend only on the mid-rapidity charged particle multiplicity. Several features of the data are reproduced qualitatively by general purpose QCD Monte Carlo models that take into account the effect of densely-packed QCD strings in high multiplicity collisions. However, none of the tested models reproduce the data quantitatively. This work corroborates and extends the ALICE findings on strangeness production in proton-proton collisions at 7 TeV.

This paper presents the measurements of $\pi^{\pm}$, $\rm{K}^{\pm}$, $\rm{p}$ and $\bar{\rm{p}}$ transverse momentum ($p_{\rm{T}}$) spectra as a function of charged-particle multiplicity density in proton-proton (pp) collisions at $\sqrt{s}$ = 13 TeV with the ALICE detector at the LHC. Such study allows us to isolate the center-of-mass energy dependence of light-flavour particle production. The measurements reported here cover a $p_{\rm{T}}$ range from 0.1 GeV/$c$ to 20 GeV/$c$ and are done in the rapidity interval $|y|<0.5$. The $p_{\rm{T}}$-differential particle ratios exhibit an evolution with multiplicity, similar to that observed in pp collisions at $\sqrt{s}$ = 7 TeV, which is qualitatively described by some of the hydrodynamical and pQCD-inspired models discussed in this paper. Furthermore, the $p_{\rm{T}}$-integrated hadron-to-pion yield ratios measured in pp collisions at two different center-of-mass energies are consistent when compared at similar multiplicities. This also extends to strange and multistrange hadrons, suggesting that, at LHC energies, particle hadrochemistry scales with particle multiplicity the same way under different collision energies and colliding systems.

We present results on transverse momentum ($p_{\rm T}$) and rapidity ($y$) differential production cross sections, mean transverse momentum and mean transverse momentum square of inclusive J/$\psi$ and $\psi(2S)$ at forward rapidity ($2.5<y<4$) as well as $\psi(2S)$-to-J/$\psi$ cross section ratios. These quantities are measured in pp collisions at center of mass energies $\sqrt{s}=5.02$ and 13 TeV with the ALICE detector. Both charmonium states are reconstructed in the dimuon decay channel, using the muon spectrometer. A comprehensive comparison to inclusive charmonium cross sections measured at $\sqrt{s}=2.76$, 7 and 8 TeV is performed. A comparison to non-relativistic quantum chromodynamics and fixed-order next-to-leading logarithm calculations, which describe prompt and non-prompt charmonium production respectively, is also presented. A good description of the data is obtained over the full $p_{\rm T}$ range, provided that both contributions are summed. In particular, it is found that for $p_{\rm T}>15$ GeV/$c$ the non-prompt contribution reaches up to 50% of the total charmonium yield.

The production of $\pi^{\pm}$, $\rm{K}^{\pm}$, $\rm{K}^{0}_{S}$, $\rm{K}^{*}(892)^{0}$, $\rm{p}$, $\phi(1020)$, $\Lambda$, $\Xi^{-}$, $\Omega^{-}$, and their antiparticles was measured in inelastic proton-proton (pp) collisions at a center-of-mass energy of $\sqrt{s}$ = 13 TeV at midrapidity ($|y|<0.5$) as a function of transverse momentum ($p_{\rm{T}}$) using the ALICE detector at the CERN LHC. Furthermore, the single-particle $p_{\rm{T}}$ distributions of $\rm{K}^{0}_{S}$, $\Lambda$, and $\overline{\Lambda}$ in inelastic pp collisions at $\sqrt{s}$ = 7 TeV are reported here for the first time. The $p_{\rm{T}}$ distributions are studied at midrapidity within the transverse momentum range $0\leq p_{\rm{T}}\leq20$ GeV/$c$, depending on the particle species. The $p_{\rm{T}}$ spectra, integrated yields, and particle yield ratios are discussed as a function of collision energy and compared with measurements at lower $\sqrt{s}$ and with results from various general-purpose QCD-inspired Monte Carlo models. A hardening of the spectra at high $p_{\rm{T}}$ with increasing collision energy is observed, which is similar for all particle species under study. The transverse mass and $x_{\rm{T}}\equiv2p_{\rm{T}}/\sqrt{s}$ scaling properties of hadron production are also studied. As the collision energy increases from $\sqrt{s}$ = 7 to 13 TeV, the yields of non- and single-strange hadrons normalized to the pion yields remain approximately constant as a function of $\sqrt{s}$, while ratios for multi-strange hadrons indicate enhancements. The $p_{\rm{T}}$-differential cross sections of $\pi^{\pm}$, $\rm{K}^{\pm}$ and $\rm{p}$ ($\overline{\rm{p}}$) are compared with next-to-leading order perturbative QCD calculations, which are found to overestimate the cross sections for $\pi^{\pm}$ and $\rm{p}$ ($\overline{\rm{p}}$) at high $p_{\rm{T}}$.

The production cross section of inclusive J/$\psi$ pairs in pp collisions at a centre-of-mass energy $\sqrt{s} = 13$ TeV is measured with ALICE. The measurement is performed for J/$\psi$ in the rapidity interval $2.5 < y < 4.0$ and for transverse momentum $p_{\rm T} > 0$. The production cross section of inclusive J/$\psi$ pairs is reported to be $10.3 \pm 2.3 {\rm (stat.)} \pm 1.3 {\rm (syst.)}$ nb in this kinematic interval. The contribution from non-prompt J/$\psi$ (i.e. originated from beauty-hadron decays) to the inclusive sample is evaluated. The results are discussed and compared with data.

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.

This work presents new constraints on the existence and the binding energy of a possible $\Lambda$-$\Lambda$ bound state, the H-dibaryon, derived from $\Lambda$-$\Lambda$ femtoscopic measurements by the ALICE collaboration. The results are obtained from a new measurement using the femtoscopy technique in pp collisions at $\sqrt{s}=13$ TeV and p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV, combined with previously published results from p-Pb collisions at $\sqrt{s}=7$ TeV. The $\Lambda$-$\Lambda$ scattering parameter space, spanned by the inverse scattering length $f_0^{-1}$ and the effective range $d_0$, is constrained by comparing the measured $\Lambda$-$\Lambda$ correlation function with calculations obtained within the Lednicky model. The data are compatible with hypernuclei results and lattice computations, both predicting a shallow attractive interaction, and permit to test different theoretical approaches describing the $\Lambda$-$\Lambda$ interaction. The region in the $(f_0^{-1},d_0)$ plane which would accommodate a $\Lambda$-$\Lambda$ bound state is substantially restricted compared to previous studies. The binding energy of the possible $\Lambda$-$\Lambda$ bound state is estimated within an effective-range expansion approach and is found to be $B_{\Lambda\Lambda}=3.2^{+1.6}_{-2.4}\mathrm{(stat)}^{+1.8}_{-1.0}\mathrm{(syst)}$ MeV.

Measurements of identified hadrons as a function of the charged-particle multiplicity in pp collisions enable a search for the onset of collective effects in small collision systems. With such measurements, it is possible to study the mechanisms that determine the shapes of hadron transverse momentum ($p_{\rm{T}}$) spectra, to search for possible modifications of the yields of short-lived hadronic resonances due to scattering effects in the hadron-gas phase, and to investigate different explanations for the multiplicity evolution of strangeness production provided by phenomenological models. In this paper, these topics are addressed through measurements of the $\rm{K}^{*}(892)^{0}$ and $\phi(1020)$ mesons at midrapidity in pp collisions at $\sqrt{s}$ = 13 TeV as a function of the charged-particle multiplicity. The results include the $p_{\rm{T}}$ spectra, $p_{\rm{T}}$-integrated yields, mean transverse momenta, and the ratios of the yields of these resonances to those of longer-lived hadrons. Comparisons with results from other collision systems and energies, as well as predictions from phenomenological models, are also discussed.

The pseudorapidity ($\eta$) and transverse-momentum ($p_{\rm T}$) distributions of charged particles produced in proton-proton collisions are measured at the centre-of-mass energy $\sqrt{s}$ = 13 TeV. The pseudorapidity distribution in $|\eta|<$ 1.8 is reported for inelastic events and for events with at least one charged particle in $|\eta|<$ 1. The pseudorapidity density of charged particles produced in the pseudorapidity region $|\eta|<$ 0.5 is 5.31 $\pm$ 0.18 and 6.46 $\pm$ 0.19 for the two event classes, respectively. The transverse-momentum distribution of charged particles is measured in the range 0.15 $<$ $p_{\rm T}$ $<$ 20 GeV/c and $|\eta|<$ 0.8 for events with at least one charged particle in $|\eta|<$ 1. The correlation between transverse momentum and particle multiplicity is also investigated by studying the evolution of the spectra with event multiplicity. The results are compared with calculations from PYTHIA and EPOS Monte Carlo generators.

The pseudorapidity density of charged particles with minimum transverse momentum ($p_{\rm T}$) thresholds of 0.15, 0.5, 1, and 2 GeV$/c$ is measured in pp collisions at the centre of mass energies of $\sqrt{s} =$ 5.02 and 13 TeV with the ALICE detector. The study is carried out for inelastic collisions with at least one primary charged particle having a pseudorapidity ($\eta$) within $\pm0.8$ and $p_{\rm T}$ larger than the corresponding threshold. In addition, measurements without $p_{\rm T}$-thresholds are performed for inelastic and non-single-diffractive events as well as for inelastic events with at least one charged particle having $|\eta|<1$ in pp collisions at $\sqrt{s} =$ 5.02 TeV for the first time at the LHC. These measurements are compared to the PYTHIA 6, PYTHIA 8, and EPOS-LHC models. In general, the models describe the $\eta$ dependence of particle production well. However, discrepancies are observed for the highest transverse momentum threshold ($p_{\rm T}>2 {\rm\ GeV}/c$), highlighting the importance of such measurements for tuning event generators. The new measurements agree within uncertainties with results from the ATLAS and CMS experiments obtained at $\sqrt{s} = 13$ TeV.

The energy deposited at very forward rapidities (very forward energy) is a powerful tool for characterising proton fragmentation in pp and p$-$Pb collisions. The correlation of very forward energy with particle production at midrapidity provides direct insights into the initial stages and the subsequent evolution of the collision. Furthermore, the correlation with the production of particles with large transverse momenta at midrapidity provides information complementary to the measurements of the underlying event, which are usually interpreted in the framework of models implementing centrality-dependent multiple parton interactions. Results about very forward energy, measured by the ALICE zero degree calorimeters (ZDCs), and its dependence on the activity measured at midrapidity in pp collisions at $\sqrt{s}=13$ TeV and in p$-$Pb collisions at $\sqrt{s_{\rm{NN}}}=8.16$ TeV are discussed. The measurements performed in pp collisions are compared with the expectations of three hadronic interaction event generators: PYTHIA 6 (Perugia 2011 tune), PYTHIA 8 (Monash tune), and EPOS LHC. These results provide new constraints on the validity of models in describing the beam remnants at very forward rapidities, where perturbative QCD cannot be used.

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.

Multiplicity ($N_{\rm ch}$) distributions and transverse momentum ($p_{\rm T}$) spectra of inclusive primary charged particles in the kinematic range of $|\eta| < 0.8$ and 0.15 GeV/$c$$< p_{T} <$ 10 GeV/$c$ are reported for pp, p-Pb, Xe-Xe and Pb-Pb collisions at centre-of-mass energies per nucleon pair ranging from $\sqrt{s_{\rm NN}} = 2.76$ TeV up to $13$ TeV. A sequential two-dimensional unfolding procedure is used to extract the correlation between the transverse momentum of primary charged particles and the charged-particle multiplicity of the corresponding collision. This correlation sharply characterises important features of the final state of a collision and, therefore, can be used as a stringent test of theoretical models. The multiplicity distributions as well as the mean and standard deviation derived from the $p_{\rm T}$ spectra are compared to state-of-the-art model predictions. Providing these fundamental observables of bulk particle production consistently across a wide range of collision energies and system sizes can serve as an important input for tuning Monte Carlo event generators.

The measurement of dielectron production is presented as a function of invariant mass and transverse momentum ($p_{\rm T}$) at midrapidity ($|y_{\rm e}|<0.8$) in proton-proton (pp) collisions at a centre-of-mass energy of $\sqrt{s}=13$ TeV. The contributions from light-hadron decays are calculated from their measured cross sections in pp collisions at $\sqrt{s}=7$ TeV or 13 TeV. The remaining continuum stems from correlated semileptonic decays of heavy-flavour hadrons. Fitting the data with templates from two different MC event generators, PYTHIA and POWHEG, the charm and beauty cross sections at midrapidity are extracted for the first time at this collision energy: ${\rm d}\sigma_{\rm c\bar{c}}/{\rm d}y|_{y=0}=974\pm138(\rm{stat.})\pm140(\rm{syst.})~\mu{\rm b}$ and ${\rm d}\sigma_{\rm b\bar{b}}/{\rm d}y|_{y=0}=79\pm14(\rm{stat.})\pm11(\rm{syst.})~\mu{\rm b}$ using PYTHIA simulations and ${\rm d}\sigma_{\rm c\bar{c}}/{\rm d}y|_{y=0}=1417\pm184(\rm{stat.})\pm204(\rm{syst.})~\mu{\rm b}$ and ${\rm d}\sigma_{\rm b\bar{b}}/{\rm d}y|_{y=0}=48\pm14(\rm{stat.})\pm7(\rm{syst.})~\mu{\rm b}$ for POWHEG. These values, whose uncertainties are fully correlated between the two generators, are consistent with extrapolations from lower energies. The different results obtained with POWHEG and PYTHIA imply different kinematic correlations of the heavy-quark pairs in these two generators. Furthermore, comparisons of dielectron spectra in inelastic events and in events collected with a trigger on high charged-particle multiplicities are presented in various $p_{\rm T}$ intervals. The differences are consistent with the already measured scaling of light-hadron and open-charm production at high charged-particle multiplicity as a function of $p_{\rm T}$. Upper limits for the contribution of virtual direct photons are extracted at 90% confidence level and found to be in agreement with pQCD calculations.

Collective behavior has been observed in high-energy heavy-ion collisions for several decades. Collectivity is driven by the high particle multiplicities that are produced in these collisions. At the Large Hadron Collider (LHC), features of collectivity have also been seen in high-multiplicity proton-proton collisions that can attain particle multiplicities comparable to peripheral Pb-Pb collisions. One of the possible signatures of collective behavior is the decrease of femtoscopic radii extracted from pion and kaon pairs emitted from high-multiplicity collisions with increasing pair transverse momentum. This decrease can be described in terms of an approximate transverse mass scaling. In the present work, femtoscopic analyses are carried out by the ALICE collaboration on charged pion and kaon pairs produced in pp collisions at $\sqrt{s}=13$ TeV from the LHC to study possible collectivity in pp collisions. The event-shape analysis method based on transverse sphericity is used to select for spherical versus jet-like events, and the effects of this selection on the femtoscopic radii for both charged pion and kaon pairs are studied. This is the first time this selection method has been applied to charged kaon pairs. An approximate transverse-mass scaling of the radii is found in all multiplicity ranges studied when the difference in the Lorentz boost for pions and kaons is taken into account. This observation does not support the hypothesis of collective expansion of hot and dense matter that should only occur in high-multiplicity events. A possible alternate explanation of the present results is based on a scenario of common emission conditions for pions and kaons in pp collisions for the multiplicity ranges studied.

Measurements of anisotropic flow coefficients ($v_n$) and their cross-correlations using two- and multi-particle cumulant methods are reported in collisions of pp at $\sqrt{s} = 13$ TeV, p-Pb at $\sqrt{s_{_{\rm NN}}} = 5.02$ TeV, Xe-Xe at $\sqrt{s_{_{\rm NN}}} = 5.44$ TeV, and Pb-Pb at $\sqrt{s_{_{\rm NN}}} = 5.02$ TeV recorded with the ALICE detector. The multiplicity dependence of $v_n$ is studied in a very wide range from 20 to 3000 particles produced in the mid-rapidity region $|\eta|<0.8$ for the transverse momentum range $0.2 < p_{\rm T} < 3.0$ GeV/$c$. An ordering of the coefficients $v_2 > v_3 > v_4$ is found in pp and p-Pb collisions, similar to that seen in large collision systems, while a weak $v_2$ multiplicity dependence is observed relative to nucleus-nucleus collisions in the same multiplicity range. Using a novel subevent method, $v_{2}$ measured with four-particle cumulants is found to be compatible with that from six-particle cumulants in pp and p-Pb collisions. The magnitude of the correlation between $v_n^2$ and $v_m^2$, evaluated with the symmetric cumulants SC$(m,n)$ is observed to be positive at all multiplicities for $v_2$ and $v_4$, while for $v_2$ and $v_3$ it is negative and changes sign for multiplicities below 100, which may indicate a different $v_n$ fluctuation pattern in this multiplicity range. The observed long-range multi-particle azimuthal correlations in high multiplicity pp and p-Pb collisions can neither be described by PYTHIA 8 nor by IP-Glasma+MUSIC+UrQMD model calculations, and hence provide new insights into the understanding of collective effects in small collision systems.

The production of J/$\psi$ is measured at midrapidity ($|y|<0.9$) in proton-proton collisions at $\sqrt{s}$ = 5.02 and 13 TeV, through the dielectron decay channel, using the ALICE detector at the Large Hadron Collider. The data sets used for the analyses correspond to integrated luminosities of $\mathcal{L}_{\rm int}$ = 19.4 $\pm$ 0.4 nb$^{-1}$ and $\mathcal{L}_{\rm int}$ = 32.2 $\pm$ 0.5 nb$^{-1}$ at $\sqrt{s}$ = 5.02 and 13 TeV, respectively. The fraction of non-prompt J/$\psi$ mesons, i.e. those originating from the decay of beauty hadrons, is measured down to a transverse momentum $p_{\rm T}$ = 2 GeV/$c$ (1 GeV/$c$) at $\sqrt{s}$ = 5.02 TeV (13 TeV). The $p_{\rm T}$ and rapidity ($y$) differential cross sections, as well as the corresponding values integrated over $p_{\rm T}$ and $y$, are carried out separately for prompt and non-prompt J/$\psi$ mesons. The results are compared with measurements from other experiments and theoretical calculations based on quantum chromodynamics (QCD). The shapes of the $p_{\rm T}$ and $y$ distributions of beauty quarks predicted by state-of-the-art perturbative QCD models are used to extrapolate an estimate of the $\rm b\overline{b}$ pair cross section at midrapidity and in the total phase space. The total $\rm b\overline{b}$ cross sections are found to be $\sigma_{\rm b \overline{\rm b}} = 541 \pm 45 (\rm stat.) \pm 69 (\rm syst.)_{-12}^{+10} (\rm extr.)~{\rm \mu b}$ and $\sigma_{\rm b \overline{\rm b}}~=~218 \pm 37 (\rm stat.) \pm 31 (\rm syst.)_{-9.1}^{+8.2} (\rm extr.)~{\rm \mu b}$ at $\sqrt{s}$ = 13 and 5.02 TeV, respectively. The value obtained from the combination of ALICE and LHCb measurements in pp collisions at $\sqrt{s}$ = 13 TeV is also provided.

The production of J/$\psi$ is measured as a function of charged-particle multiplicity at forward rapidity in proton$-$proton (pp) collisions at center-of-mass energies $\sqrt{s} =$ 5.02 and 13 TeV. The J/$\psi$ mesons are reconstructed via their decay into dimuons in the rapidity interval (2.5 $< y <$ 4.0), whereas the charged-particle multiplicity density (${\rm d}N_{\rm{ch}}/{\rm d}\eta$) is measured at midrapidity $(|\eta| < 1)$. The production rate as a function of multiplicity is reported as the ratio of the yield in a given multiplicity interval to the multiplicity-integrated one. This observable shows a linear increase with charged-particle multiplicity normalized to the corresponding average value for inelastic events (${{\rm d}N_{\rm ch}/{\rm d}\eta}/{\langle {\rm d}N_{\rm ch}/{\rm d}\eta \rangle}$), at both the colliding energies. Measurements are compared with available ALICE results at midrapidity and theoretical model calculations. First measurement of the mean transverse momentum ($\langle p_{\mathrm{T}}\rangle$) of J/$\psi$ in pp collisions exhibits an increasing trend as a function of ${{\rm d}N_{\rm ch}/{\rm d}\eta}/{\langle {\rm d}N_{\rm ch}/{\rm d}\eta \rangle}$ showing a saturation towards high charged-particle multiplicities.

Two-particle angular correlations are measured in high-multiplicity proton-proton collisions at $\sqrt{s} =13$ TeV by the ALICE Collaboration. The yields of particle pairs at short-($\Delta\eta$$\sim$ 0) and long-range ($1.6 < |\Delta\eta| < 1.8$) in pseudorapidity are extracted on the near-side ($\Delta\varphi$$\sim$ 0). They are reported as a function of transverse momentum ($p_{\mathrm T}$) in the range $1 <p_{\mathrm T}< 4$ GeV/$c$. Furthermore, the event-scale dependence is studied for the first time by requiring the presence of high-$p_{\rm T}$ leading particles and jets for varying $p_{\rm T}$ thresholds. The results demonstrate that the long-range "ridge" yield, possibly related to the collective behavior of the system, is present in events with high-$p_{\mathrm T}$ processes. The magnitudes of the short- and long-range yields are found to grow with the event scale. The results are compared to EPOS LHC and PYTHIA 8 calculations, with and without string-shoving interactions. It is found that while both models describe the qualitative trends in the data, calculations from EPOS LHC show a better quantitative agreement, in particular for the $p_{\rm T}$ and event-scale dependencies.

Angular correlations of heavy-flavour and charged particles in high-energy proton-proton collisions are sensitive to the production mechanisms of heavy quarks and to their fragmentation as well as hadronisation processes. The measurement of the azimuthal-correlation function of prompt D mesons with charged particles in proton-proton collisions at a centre-of-mass energy of $\sqrt{s} = 13$ TeV with the ALICE detector is reported, considering D$^{0}$, D$^{+}$, and D$^{*+}$ mesons in the transverse-momentum interval $3 < p_{\rm T} < 36$ GeV/$c$ at midrapidity ($|y| < 0.5$), and charged particles with $p_{\rm T} > 0.3$ GeV/$c$ and pseudorapidity $|\eta| < 0.8$. This measurement has an improved precision and provides an extended transverse-momentum coverage compared to previous ALICE measurements at lower energies. The study is also performed as a function of the charged-particle multiplicity, showing no modifications of the correlation function with multiplicity within uncertainties. The properties and the transverse-momentum evolution of the near- and away-side correlation peaks are studied and compared with predictions from various Monte Carlo event generators. Among those considered, PYTHIA8 and POWHEG+PYTHIA8 provide the best description of the measured observables. The obtained results can provide guidance on tuning the generators.

The production of prompt $D^{0}$, $D^{+}_{\rm s}$, and $\Lambda_{\rm c}^{+}$ hadrons, and their ratios, $D^{+}_{\rm s}$/$D^{0}$ and $\Lambda_{\rm c}^{+}$/$D^{0}$, are measured in proton-proton collisions at $\sqrt{s}$ = 13 TeV at midrapidity ($|y| <0.5$) with the ALICE detector at the LHC. The measurements are performed as a function of the charm-hadron transverse momentum ($p_{\rm T}$) in intervals of charged-particle multiplicity, measured with two multiplicity estimators covering different pseudorapidity regions. While the strange to non-strange $D^{+}_{\rm s}$/$D^{0}$ ratio indicates no significant multiplicity dependence, the baryon-to-meson $p_{\rm T}$-differential $\Lambda_{\rm c}^{+}$/$D^{0}$ ratio shows a multiplicity-dependent enhancement, with a significance of 5.3$\sigma$ for $1< p_{\rm T} < 12$ GeV/$c$, comparing the highest multiplicity interval with respect to the lowest one. The measurements are compared with a theoretical model that explains the multiplicity dependence by a canonical treatment of quantum charges in the statistical hadronisation approach, and with predictions from event generators that implement colour reconnection mechanisms beyond the leading colour approximation to model the hadronisation process. The $\Lambda_{\rm c}^{+}$/$D^{0}$ ratios as a function of $p_{\rm T}$ present a similar shape and magnitude as the $\Lambda/K^{0}_{s}$ ratios in comparable multiplicity intervals, suggesting a potential common mechanism for light- and charm-hadron formation, with analogous multiplicity dependence. The $p_{\rm T}$-integrated ratios, extrapolated down to $p_{\rm T}$=0, do not show a significant dependence on multiplicity within the uncertainties.

The $p_{\rm T}$-differential cross sections of prompt charm-strange baryons $\Xi^0_{\rm c}$ and $\Xi^+_{\rm c}$ were measured at midrapidity ($|y| < 0.5$) in proton$-$proton (pp) collisions at a centre-of-mass energy $\sqrt{s}=13$ TeV with the ALICE detector at the LHC. The $\Xi^0_{\rm c}$ baryon was reconstructed via both the semileptonic decay ($\Xi^{-}{\rm e^{+}}\nu_{\rm e}$) and the hadronic decay ($\Xi^{-}{\rm \pi^{+}}$) channels. The $\Xi^+_{\rm c}$ baryon was reconstructed via the hadronic decay ($\Xi^{-}\pi^{+}\pi^{+}$) channel. The branching-fraction ratio $\rm {\rm BR}(\Xi_c^0\rightarrow \Xi^-e^+\nu_e)/\rm {\rm BR}(\Xi_c^0\rightarrow \Xi^{-}\pi^+)=$ 1.38 $\pm$ 0.14 (stat) $\pm$ 0.22 (syst) was measured with a total uncertainty reduced by a factor of about 3 with respect to the current world average reported by the Particle Data Group. The transverse momentum ($p_{\rm T}$) dependence of the $\Xi^0_{\rm c}$- and $\Xi^+_{\rm c}$-baryon production relative to the ${\rm D^0}$-meson and to the $\Sigma^{0,+,++}_{\rm c}$- and $\Lambda^+_{\rm c}$-baryon production are reported. The baryon-to-meson ratio increases towards low $p_{\rm T}$ up to a value of approximately 0.3. The measurements are compared with various models that take different hadronisation mechanisms into consideration. The results provide stringent constraints to these theoretical calculations and additional evidence that different processes are involved in charm hadronisation in electron$-$positron ($\rm e^+e^-$) and hadronic collisions.

A measurement of dielectron production in proton-proton (pp) collisions at $\sqrt{s} = 13$ TeV, recorded with the ALICE detector at the CERN LHC, is presented in this Letter. The data set was recorded with a reduced magnetic solenoid field. This enables the investigation of a kinematic domain at low dielectron invariant mass $m_{\rm ee}$ and pair transverse momentum $p_{\rm T,ee}$ that was previously inaccessible at the LHC. The cross section for dielectron production is studied as a function of $m_{\rm ee}$, $p_{\rm T,ee}$, and event multiplicity ${\rm d} N_{\rm ch}/{\rm d} \eta$. The expected dielectron rate from hadron decays, called hadronic cocktail, utilizes a parametrization of the measured $\eta/\pi^0$ ratio in pp and proton-nucleus (p-A) collisions, assuming that this ratio shows no strong dependence on collision energy at low transverse momentum. Comparison of the measured dielectron yield to the hadronic cocktail at $0.15<m_{\rm ee}<0.6$ GeV/$c^2$ and for $p_{\rm T,ee} < 0.4$ GeV/$c$ indicates an enhancement of soft dielectrons, reminiscent of the 'anomalous' soft-photon and -dilepton excess in hadron-hadron collisions reported by several experiments under different experimental conditions. The enhancement factor over the hadronic cocktail amounts to $1.61\pm 0.13\,(\rm{stat.})\pm 0.17\,(\rm{syst., data}) \pm 0.34\,(\rm{syst., cocktail})$ in the ALICE acceptance. Acceptance-corrected excess spectra in $m_{\rm ee}$ and $p_{\rm T,ee}$ are extracted and compared with calculations of dielectron production from hadronic bremsstrahlung and thermal radiation within a hadronic many-body approach.

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.

Measurements of the production of electrons from heavy-flavour hadron decays in pp collisions at $\sqrt{s} = 13$ TeV at midrapidity with the ALICE detector are presented down to a transverse momentum ($p_{\rm T}$) of 0.2 GeV$/c$ and up to $p_{\rm T} = 35$ GeV$/c$, which is the largest momentum range probed for inclusive electron measurements in ALICE. In p$-$Pb collisions, the production cross section and the nuclear modification factor of electrons from heavy-flavour hadron decays are measured in the $p_{\rm T}$ range $0.5 < p_{\rm T} < 26$ GeV$/c$ at $\sqrt{s_{\rm NN}} = 8.16$ TeV. The nuclear modification factor is found to be consistent with unity within the statistical and systematic uncertainties. In both collision systems, first measurements of the yields of electrons from heavy-flavour hadron decays in different multiplicity intervals normalised to the multiplicity-integrated yield (self-normalised yield) at midrapidity are reported as a function of the self-normalised charged-particle multiplicity estimated at midrapidity. The self-normalised yields in pp and p$-$Pb collisions grow faster than linear with the self-normalised multiplicity. A strong $p_{\rm T}$ dependence is observed in pp collisions, where the yield of high-$p_{\rm T}$ electrons increases faster as a function of multiplicity than the one of low-$p_{\rm T}$ electrons. The measurement in p$-$Pb collisions shows no $p_{\rm T}$ dependence within uncertainties. The self-normalised yields in pp and p$-$Pb collisions are compared with measurements of other heavy-flavour, light-flavour, and strange particles, and with Monte Carlo simulations.

The multiplicity dependence of the pseudorapidity density of charged particles in proton-proton (pp) collisions at centre-of-mass energies $\sqrt{s}$ = 5.02, 7 and 13 TeV measured by ALICE is reported. The analysis relies on track segments measured in the midrapidity range ($|\eta| < 1.5$). Results are presented for inelastic events having at least one charged particle produced in the pseudorapidity interval $|\eta|<1$ ($\mathrm{INEL}_{>0}$). The multiplicity dependence of the pseudorapidy density of charged particles is measured with mid and forward rapidity multiplicity estimators, the latter being less affected by autocorrelations. A detailed comparison with predictions from the PYTHIA 8 and EPOS LHC event generators is also presented. Both generators provide a good description of the data.

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 production cross-sections of $\Upsilon(1S)$, $\Upsilon(2S)$ and $\Upsilon(3S)$ mesons in proton-proton collisions at $\sqrt{s}$= 13 TeV are measured with a data sample corresponding to an integrated luminosity of $277 \pm 11$ $\rm pb^{-1}$ recorded by the LHCb experiment in 2015. The $\Upsilon$ mesons are reconstructed in the decay mode $\Upsilon\to\mu^{+}\mu^{-}$. The differential production cross-sections times the dimuon branching fractions are measured as a function of the $\Upsilon$ transverse momentum, $p_{\rm T}$, and rapidity, $y$, over the range $0 < p_{\rm T}< 30$ GeV/c and $2.0 < y < 4.5$. The ratios of the cross-sections with respect to the LHCb measurement at $\sqrt{s}$= 8 TeV are also determined. The measurements are compared with theoretical predictions based on NRQCD.

The production of $J/\psi$ mesons in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=13$ TeV is studied with the \lhcb detector. Cross-section measurements are performed as a function of the transverse momentum $p_\mathrm{T}$ and the rapidity $y$ of the $J/\psi$ meson in the region $p_\mathrm{T}<14\mathrm{GeV}/c$ and $2.0<y<4.5$, for both prompt $J/\psi$ mesons and $J/\psi$ mesons from $b$-hadron decays. The production cross-sections integrated over the kinematic coverage are $15.03\pm 0.03\pm 0.94$ $\mu$b for prompt $J/\psi$ and $2.25\pm 0.01\pm 0.14$ $\mu$b for $J/\psi$ from $b$-hadron decays, assuming zero polarization of the $J/\psi$ meson. The first uncertainties are statistical and the second systematic. The cross-section reported for $J/\psi$ mesons from $b$-hadron decays is used to extrapolate to a total $b\bar{b}$ cross-section. The ratios of the cross-sections with respect to $\sqrt{s}=8$ TeV are also determined.

Production cross-sections of prompt charm mesons are measured with the first data from $pp$ collisions at the LHC at a centre-of-mass energy of $13\,\mathrm{TeV}$. The data sample corresponds to an integrated luminosity of $4.98 \pm 0.19\,\mathrm{pb}^{-1}$ collected by the LHCb experiment. The production cross-sections of $D^{0}$, $D^{+}$, $D_{s}^{+}$, and $D^{*+}$ mesons are measured in bins of charm meson transverse momentum, $p_{\mathrm{T}}$, and rapidity, $y$, and cover the range $0 < p_{\mathrm{T}} < 15\,\mathrm{GeV}/c$ and $2.0 < y < 4.5$. The inclusive cross-sections for the four mesons, including charge conjugation, within the range of $1 < p_{\mathrm{T}} < 8\,\mathrm{GeV}/c$ are found to be \begin{equation} \sigma(pp \to D^{0} X) = 2072 \pm 2 \pm 124\,\mu\mathrm{b}\\ \sigma(pp \to D^{+} X) = 834 \pm 2 \pm \phantom{1}78\,\mu\mathrm{b}\\ \sigma(pp \to D_{s}^{+} X) = 353 \pm 9 \pm \phantom{1}76\,\mu\mathrm{b}\\ \sigma(pp \to D^{*+} X) = 784 \pm 4 \pm \phantom{1}87\,\mu\mathrm{b} \end{equation} where the uncertainties are due to statistical and systematic uncertainties, respectively.

Using a data sample corresponding to an integrated luminosity of $2.0\,fb^{-1}$, collected by the LHCb experiment, the production of the $\eta_c(1S)$ state in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=13 \text{ TeV}$ is studied in the rapidity range ${2.0 < y < 4.5}$ and in the transverse momentum range ${6.5 < p_{T} < 14.0\text{ GeV}}$. The cross-section for prompt production of $\eta_c(1S)$ mesons relative to that of the $J/\psi$ meson is measured using the ${p\bar{p}}$ decay mode and is found to be ${\sigma_{\eta_c(1S)}/\sigma_{J/\psi} = 1.69 \pm 0.15 \pm 0.10 \pm 0.18}$. The quoted uncertainties are, in order, statistical, systematic and due to uncertainties on the branching fractions of the ${J/\psi\to p \bar{p}}$ and ${\eta_c\to p \bar{p}}$ decays. The prompt $\eta_c(1S)$ production cross-section is determined to be ${\sigma_{\eta_c(1S)} = 1.26 \pm 0.11\pm 0.08 \pm 0.14 \,\mu b}$, where the last uncertainty includes that on the ${J/\psi}$ meson cross-section. The ratio of the branching fractions of $b$-hadron decays to the $\eta_c(1S)$ and ${J/\psi}$ states is measured to be ${\mathcal{B}_{b\to\eta_c X}/\mathcal{B}_{b\to J/\psi X} = 0.48 \pm 0.03 \pm 0.03 \pm 0.05}$, where the last uncertainty is due to those on the branching fractions of the ${J/\psi \to p \bar{p}}$ and ${\eta_c\to p \bar{p}}$ decays. The difference between the ${J/\psi}$ and $\eta_c(1S)$ masses is also determined to be ${113.0 \pm 0.7 \pm 0.1\text{ MeV}}$, which is the most precise single measurement of this quantity to date.

A measurement of the production cross-section of Z bosons in pp collisions at $\sqrt{s} = 13$ TeV is presented using dimuon and dielectron final states in LHCb data. The cross-section is measured for leptons with pseudorapidities in the range $2.0 < \eta < 4.5$, transverse momenta $p_\text{T} > 20$ GeV and dilepton invariant mass in the range $60<m(\ell\ell)<120$ GeV. The integrated cross-section from averaging the two final states is \begin{equation*}\sigma_{\text{Z}}^{\ell\ell} = 194.3 \pm 0.9 \pm 3.3 \pm 7.6\text{ pb,}\end{equation*} where the first uncertainty is statistical, the second is due to systematic effects, and the third is due to the luminosity determination. In addition, differential cross-sections are measured as functions of the Z boson rapidity, transverse momentum and the angular variable $\phi^*_\eta$.

Forward top quark pair production is studied in $pp$ collisions in the $\mu eb$ final state using a data sample corresponding to an integrated luminosity of 1.93 fb$^{-1}$ collected with the LHCb experiment at a centre-of-mass energy of 13 TeV. The cross-section is measured in a fiducial region where both leptons have a transverse momentum greater than 20 GeV and a pseudorapidity between 2.0 and 4.5. The quadrature sum of the azimuthal separation and the difference in pseudorapidities, denoted $\Delta R$, between the two leptons must be larger than 0.1. The $b$-jet axis is required to be separated from both leptons by a $\Delta R$ of 0.5, and to have a transverse momentum in excess of 20 GeV and a pseudorapidity between 2.2 and 4.2. The cross-section is measured to be $$\sigma_{t\bar{t}}= 126\pm19\,(\mathrm{stat})\pm16\,(\mathrm{syst})\pm5\,(\mathrm{lumi})\,\,\mathrm{ fb}$$ where the first uncertainty is statistical, the second is systematic, and the third is due to the luminosity determination. The measurement is compatible with the Standard Model prediction.

The production of $B^{\pm}$ mesons is studied in $pp$ collisions at centre-of-mass energies of 7 and 13 TeV, using $B^{\pm}\rightarrow J/\psi K^{\pm}$ decays and data samples corresponding to 1.0 fb$^{-1}$ and 0.3 fb$^{-1}$, respectively. The production cross-sections summed over both charges and integrated over the transverse momentum range $0<p_{\text{T}}< 40$ GeV/$c$ and the rapidity range $2.0<y<4.5$ are measured to be $\sigma(pp \rightarrow B^{\pm} X, \sqrt{s} = \text{7 TeV}) = 43.0 \pm 0.2 \pm 2.5 \pm 1.7\mu b,$ $\sigma(pp \rightarrow B^{\pm} X, \sqrt{s} = \text{13 TeV}) = 86.6 \pm 0.5 \pm 5.4 \pm 3.4\mu b,$ where the first uncertainties are statistical, the second are systematic, and the third are due to the limited knowledge of the $B^{\pm}\rightarrow J/\psi K^{\pm}$ branching fraction. The ratio of the cross-section at 13 TeV to that at 7 TeV is determined to be $2.02\pm0.02\text{(stat)}\pm0.12\text{(syst)}$. Differential cross-sections are also reported as functions of $p_{\text{T}}$ and $y$. All results are in agreement with theoretical calculations based on the state-of-art fixed next-to-leading order quantum chromodynamics.

The inclusive $b \bar{b}$- and $c \bar{c}$-dijet production cross-sections in the forward region of $pp$ collisions are measured using a data sample collected with the LHCb detector at a centre-of-mass energy of 13 TeV in 2016. The data sample corresponds to an integrated luminosity of 1.6 fb$^{-1}$. Differential cross-sections are measured as a function of the transverse momentum and of the pseudorapidity of the leading jet, of the rapidity difference between the jets, and of the dijet invariant mass. A fiducial region for the measurement is defined by requiring that the two jets originating from the two $b$ or $c$ quarks are emitted with transverse momentum greater than 20 GeV$/c$, pseudorapidity in the range $2.2 < \eta < 4.2$, and with a difference in the azimuthal angle between the two jets greater than 1.5. The integrated $b \bar{b}$-dijet cross-section is measured to be $53.0 \pm 9.7$ nb, and the total $c \bar{c}$-dijet cross-section is measured to be $73 \pm 16$ nb. The ratio between $c \bar{c}$- and $b \bar{b}$-dijet cross-sections is also measured and found to be $1.37 \pm 0.27$. The results are in agreement with theoretical predictions at next-to-leading order.

Production of inclusive charmonia in pp collisions at center-of-mass energy of $\sqrt{s}$ = 13 TeV and p-Pb collisions at center-of-mass energy per nucleon pair of $\sqrt{s_{\rm NN}}$ = 8.16 TeV is studied as a function of charged-particle pseudorapidity density with ALICE. Ground and excited charmonium states (J/$\psi$, $\psi$(2S)) are measured from their dimuon decays in the interval of rapidity in the center-of-mass frame $2.5 < y_{\rm cms} < 4.0$ for pp collisions, and $2.03 < y_{\rm cms} < 3.53$ and $-4.46 < y_{\rm cms} < -2.96$ for p-Pb collisions. The charged-particle pseudorapidity density is measured around midrapidity ($|\eta|<1.0$). In pp collisions, the measured charged-particle multiplicity extends to about six times the average value, while in p-Pb collisions at forward (backward) rapidity a multiplicity corresponding to about three (four) times the average is reached. The $\psi$(2S) yield increases with the charged-particle pseudorapidity density. The ratio of $\psi$(2S) over J/$\psi$ yield does not show a significant multiplicity dependence in either colliding system, suggesting a similar behavior of J/$\psi$ and $\psi$(2S) yields with respect to charged-particle pseudorapidity density. Results for the $\psi$(2S) yield and its ratio with respect to J/$\psi$ agree with available model calculations.

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.