The measurement of the production of ${\rm f}_{0}(980)$ in inelastic pp collisions at $\sqrt{s} = 5.02$ TeV is presented. This is the first reported measurement of inclusive ${\rm f}_{0}(980)$ yield at LHC energies. The production is measured at midrapidity, $|y| < 0.5$, in a wide transverse momentum range, $0 < p_{\rm T} < 16$ GeV/$c$, by reconstructing the resonance in the ${\rm f}_{0}(980) \rightarrow \pi^{+}\pi^{-}$ hadronic decay channel using the ALICE detector. The $p_{\rm T}$-differential yields are compared to those of pions, protons and $\phi$ mesons as well as to predictions from the HERWIG 7.2 QCD-inspired Monte Carlo event generator and calculations from a coalescence model that uses the AMPT model as an input. The ratio of the $p_{\rm T}$-integrated yield of ${\rm f}_{0}(980)$ relative to pions is compared to measurements in ${\rm e}^{+}{\rm e}^{-}$ and pp collisions at lower energies and predictions from statistical hadronisation models and HERWIG 7.2. A mild collision energy dependence of the ${\rm f}_{0}(980)$ to pion production is observed in pp collisions from SPS to LHC energies. All considered models underpredict the $p_{\rm T}$-integrated $2{\rm f}_{0}(980)/(\pi^{+}+\pi^{-})$ ratio. The prediction from the canonical statistical hadronisation model assuming a zero total strangeness content of ${\rm f}_{0}(980)$ is consistent with the data within 1.9$\sigma$ and is the closest to the data. The results provide an essential reference for future measurements of the particle yield and nuclear modification in p$-$Pb and Pb$-$Pb collisions, which have been proposed to be instrumental to probe the elusive nature and quark composition of the ${\rm f}_{0}(980)$ scalar meson.

Modifications of the properties of jets in PbPb collisions, relative to those in pp collisions, are studied at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV via correlations of charged particles with the jet axis in relative pseudorapidity ($\Delta \eta$), relative azimuth ($\Delta \phi$), and relative angular distance from the jet axis $\Delta \mathrm{r} = \sqrt{{(\Delta\eta)^{2}+(\Delta\phi)^{2}}}$. This analysis uses data collected with the CMS detector at the LHC, corresponding to integrated luminosities of 404 $\mu$b$^{-1}$ and 27.4 pb$^{-1}$ for PbPb and pp collisions, respectively. Charged particle number densities, jet fragmentation functions, and jet shapes are presented as a function of PbPb collision centrality and charged-particle track transverse momentum, providing a differential description of jet modifications due to interactions with the quark-gluon plasma.

The first study of W boson production in pPb collisions is presented, for bosons decaying to a muon or electron, and a neutrino. The measurements are based on a data sample corresponding to an integrated luminosity of 34.6 inverse nanobarns at a nucleon-nucleon centre-of-mass energy of sqrt(s[NN]) = 5.02 TeV, collected by the CMS experiment. The W boson differential cross sections, lepton charge asymmetry, and forward-backward asymmetries are measured for leptons of transverse momentum exceeding 25 GeV, and as a function of the lepton pseudorapidity in the abs(eta[lab]) < 2.4 range. Deviations from the expectations based on currently available parton distribution functions are observed, showing the need for including W boson data in nuclear parton distribution global fits.

The production of $\Lambda$ baryons and ${\rm K}^{0}_{\rm S}$ mesons (${\rm V}^{0}$ particles) was measured in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV and pp collisions at $\sqrt{s} = 7$ TeV with ALICE at the LHC. The production of these strange particles is studied separately for particles associated with hard scatterings and the underlying event to shed light on the baryon-to-meson ratio enhancement observed at intermediate transverse momentum ($p_{\rm T}$) in high multiplicity pp and p-Pb collisions. Hard scatterings are selected on an event-by-event basis with jets reconstructed with the anti-$k_{\rm T}$ algorithm using charged particles. The production of strange particles associated with jets $p_{\rm T,\;jet}^{\rm ch}>10$ and $p_{\rm T,\;jet}^{\rm ch}>20$ GeV/$c$ in p-Pb collisions, and with jet $p_{\rm T,\;jet}^{\rm ch}>10$ GeV/$c$ in pp collisions is reported as a function of $p_{\rm T}$. Its dependence on angular distance from the jet axis, $R({\rm V}^{0},\;{\rm jet})$, for jets with $p_{\rm T,\;jet}^{\rm ch}>10$ GeV/$c$ in p-Pb collisions is reported as well. The $p_{\rm T}$-differential production spectra of strange particles associated with jets are found to be harder compared to that in the underlying event and both differ from the inclusive measurements. In events containing a jet, the density of the ${\rm V}^{0}$ particles in the underlying event is found to be larger than the density in the minimum bias events. The $\Lambda/{\rm K}^{0}_{\rm S}$ ratio associated with jets in p-Pb collisions is consistent with the ratio in pp collisions and follows the expectation of jets fragmenting in vacuum. On the other hand, this ratio within jets is consistently lower than the one obtained in the underlying event and it does not show the characteristic enhancement of baryons at intermediate $p_{\rm T}$ often referred to as "baryon anomaly" in the inclusive measurements.

The production of beauty hadrons was measured via semi-leptonic decays at mid-rapidity with the ALICE detector at the LHC in the transverse momentum interval $1<p_{\rm T}<8$ GeV/$c$ in minimum-bias p-Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV and in $1.3<p_{\rm T}<8$ GeV/$c$ in the 20% most central Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76$ TeV. The pp reference spectra at $\sqrt{s}=5.02$ TeV and $\sqrt{s}=2.76$ TeV, needed for the calculation of the nuclear modification factors $R_{\rm pPb}$ and $R_{\rm PbPb}$, were obtained by a pQCD-driven scaling of the cross section of electrons from beauty-hadron decays measured at $\sqrt{s}=7$ TeV. In the $p_{\rm T}$ interval $3<p_{\rm T}<8$ GeV/$c$ a suppression of the yield of electrons from beauty-hadron decays is observed in Pb-Pb compared to pp collisions. Towards lower $p_{\rm T}$, the $R_{\rm PbPb}$ values increase with large systematic uncertainties. The $R_{\rm pPb}$ is consistent with unity within systematic uncertainties and is well described by theoretical calculations that include cold nuclear matter effects in p-Pb collisions. The measured $R_{\rm pPb}$ and these calculations indicate that cold nuclear matter effects are small at high transverse momentum also in Pb-Pb collisions. Therefore, the observed reduction of $R_{\rm PbPb}$ below unity at high $p_{\rm T}$ may be ascribed to an effect of the hot and dense medium formed in Pb-Pb collisions.

Bose-Einstein correlations between identical particles are measured in samples of proton-proton collisions at 0.9 and 7 TeV centre-of-mass energies, recorded by the CMS experiment at the LHC. The signal is observed in the form of an enhancement of number of pairs of same-sign charged particles with small relative momentum. The dependence of this enhancement on kinematic and topological features of the event is studied.

The $p_{\rm T}$-differential production cross sections of the prompt charmed mesons $D^0$, $D^+$, $D^{*+}$ and $D_{\rm s}^{+}$ and their charge conjugate in the rapidity interval $-0.96 < y_{\rm cms} < 0.04$ were measured in p-Pb collisions at a centre-of-mass energy $\sqrt{s_{\rm NN}} = 5.02$ TeV with the ALICE detector at the LHC. The nuclear modification factor $R_{\rm pPb}$, quantifying the D-meson yield in p-Pb collisions relative to the yield in pp collisions scaled by the number of binary nucleon-nucleon collisions, is compatible within the 15-20% uncertainties with unity in the transverse momentum interval $1 < p_{\rm T} < 24$ GeV/$c$. No significant difference among the $R_{\rm pPb}$ of the four D-meson species is observed. The results are described within uncertainties by theoretical calculations that include initial-state effects. The measurement adds experimental evidence that the modification of the momentum spectrum of D mesons observed in Pb-Pb collisions with respect to pp collisions is due to strong final-state effects induced by hot partonic matter.

Three- and four-pion Bose-Einstein correlations are presented in pp, p-Pb, and Pb-Pb collisions at the LHC. We compare our measured four-pion correlations to the expectation derived from two- and three-pion measurements. Such a comparison provides a method to search for coherent pion emission. We also present mixed-charge correlations in order to demonstrate the effectiveness of several analysis procedures such as Coulomb corrections. Same-charge four-pion correlations in pp and p-Pb appear consistent with the expectations from three-pion measurements. However, the presence of non-negligible background correlations in both systems prevent a conclusive statement. In Pb-Pb collisions, we observe a significant suppression of three- and four-pion Bose-Einstein correlations compared to expectations from two-pion measurements. There appears to be no centrality dependence of the suppression within the 0-50% centrality interval. The origin of the suppression is not clear. However, by postulating either coherent pion emission or large multibody Coulomb effects, the suppression may be explained.

Production cross-sections of prompt charm mesons are measured using data from $pp$ collisions at the LHC at a centre-of-mass energy of $5\,$TeV. The data sample corresponds to an integrated luminosity of $8.60\pm0.33\,$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_{\text{T}}$, and rapidity, $y$. They cover the rapidity range $2.0<y<4.5$ and transverse momentum ranges $0 < p_{\text{T}} < 10\, \text{GeV}/c$ for $D^0$ and $D^+$ and $1 < p_{\text{T}} < 10\, \text{GeV}/c$ for $D_s^+$ and $D^{*+}$ mesons. The inclusive cross-sections for the four mesons, including charge-conjugate states, within the range of $1 < p_{\text{T}} < 8\, \text{GeV}/c$ are determined to be \sigma(pp\rightarrow D^0 X) = 1004 \pm 3 \pm 54\,\mu\text{b} \sigma(pp\rightarrow D^+ X) = 402 \pm 2 \pm 30\,\mu\text{b} \sigma(pp\rightarrow D_s^+ X) = 170 \pm 4 \pm 16\,\mu\text{b} \sigma(pp\rightarrow D^{*+} X)= 421 \pm 5 \pm 36\,\mu\text{b} where the uncertainties are statistical and systematic, respectively.

This letter presents the first measurement of jet mass in Pb-Pb and p-Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV and 5.02 TeV, respectively. Both the jet energy and the jet mass are expected to be sensitive to jet quenching in the hot Quantum Chromodynamics (QCD) matter created in nuclear collisions at collider energies. Jets are reconstructed from charged particles using the anti-$k_{\rm T}$ jet algorithm and resolution parameter $R = 0.4$. The jets are measured in the pseudorapidity range $|\eta_{\rm{jet}}|<0.5$ and in three intervals of transverse momentum between 60 GeV/$c$ and 120 GeV/$c$. The measurement of the jet mass in central Pb-Pb collisions is compared to the jet mass as measured in p-Pb reference collisions, to vacuum event generators, and to models including jet quenching. It is observed that the jet mass in central Pb-Pb collisions is consistent within uncertainties with p-Pb reference measurements. Furthermore, the measured jet mass in Pb-Pb collisions is not reproduced by the quenching models considered in this letter and is found to be consistent with PYTHIA expectations within systematic uncertainties.

Comprehensive results on the production of unidentified charged particles, $\pi^{\pm}$, $\rm{K}^{\pm}$, $\rm{K}^{0}_{S}$, $\rm{K}$*(892)$^{0}$, $\rm{p}$, $\overline{\rm{p}}$, $\phi$(1020), $\Lambda$, $\overline{\Lambda}$, $\Xi^{-}$, $\overline{\Xi}^{+}$, $\Omega^{-}$ and $\overline{\Omega}^{+}$ hadrons in proton-proton (pp) collisions at $\sqrt{s}$ = 7 TeV at midrapidity ($|y| < 0.5$) as a function of charged-particle multiplicity density are presented. In order to avoid auto-correlation biases, the actual transverse momentum ($p_{\rm{T}}$) spectra of the particles under study and the event activity are measured in different rapidity windows. In the highest multiplicity class, the charged-particle density reaches about 3.5 times the value measured in inelastic collisions. While the yield of protons normalized to pions remains approximately constant as a function of multiplicity, the corresponding ratios of strange hadrons to pions show a significant enhancement that increases with increasing strangeness content. Furthermore, all identified particle to pion ratios are shown to depend solely on charged-particle multiplicity density, regardless of system type and collision energy. The evolution of the spectral shapes with multiplicity and hadron mass shows patterns that are similar to those observed in p-Pb and Pb-Pb collisions at LHC energies. The obtained $p_{\rm{T}}$ distributions and yields are compared to expectations from QCD-based pp event generators as well as to predictions from thermal and hydrodynamic models. These comparisons indicate that traces of a collective, equilibrated system are already present in high-multiplicity pp collisions.

Inclusive J/$\psi$ production is studied in minimum-bias proton-proton collisions at a centre-of-mass energy of $\sqrt{s}$ = 5.02 TeV by ALICE at the CERN LHC. The measurement is performed at mid-rapidity ($|y| < 0.9$) in the dielectron decay channel down to zero transverse momentum $p_{\rm T}$, using a data sample corresponding to an integrated luminosity of $L_{\rm int} = 19.4 \pm$ 0.4 nb$^{-1}$. The measured $p_{\rm T}$-integrated inclusive J/$\psi$ production cross section is d$\sigma$/d$y$ = 5.64 $\pm$ 0.22 (stat.) $\pm 0.33$ (syst.) $\pm 0.12$ (lumi.) $\mu$b. The $p_{\rm T}$-differential cross section d$^{2} \sigma$/d$p_{\rm T}$d$y$ is measured in the $p_{\rm T}$ range 0$-$10 GeV/$c$ and compared with state-of-the-art perturbative QCD calculations. The J/$\psi$ $\langle p_{\rm T} \rangle$ and $\langle p_{\rm T}^{2} \rangle$ are extracted and compared with results obtained at other collision energies.

The global polarization of the $\Lambda$ and $\overline\Lambda$ hyperons is measured for Pb-Pb collisions at $\sqrt{s_{\rm{NN}}}$ = 2.76 and 5.02 TeV recorded with the ALICE at the LHC. The results are reported differentially as a function of collision centrality and hyperon's transverse momentum ($p_{\rm{T}}$) for the range of centrality 5-50%, $0.5 < p_{\rm{T}} <5$ GeV/$c$, and rapidity $|y|<0.5$. The hyperon global polarization averaged for Pb-Pb collisions at $\sqrt{s_{\rm{NN}}}$ = 2.76 and 5.02 TeV is found to be consistent with zero, $\langle P_{\rm{H}}\rangle$ (%) $\approx$ - 0.01 $\pm$ 0.05 (stat.) $\pm$ 0.03 (syst.) in the collision centrality range 15-50%, where the largest signal is expected. The results are compatible with expectations based on an extrapolation from measurements at lower collision energies at RHIC, hydrodynamical model calculations, and empirical estimates based on collision energy dependence of directed flow, all of which predict the global polarization values at LHC energies of the order of 0.01%.

Modifications to the distribution of charged particles with respect to high transverse momentum ($p_\mathrm{T}$) jets passing through a quark-gluon plasma are explored using the CMS detector. Back-to-back dijets are analyzed in lead-lead and proton-proton collisions at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV via correlations of charged particles in bins of relative pseudorapidity and angular distance from the leading and subleading jet axes. In comparing the lead-lead and proton-proton collision results, modifications to the charged-particle relative distance distribution and to the momentum distributions around the jet axis are found to depend on the dijet momentum balance $x_j$, which is the ratio between the subleading and leading jet $p_\mathrm{T}$. For events with $x_j$$\approx$ 1, these modifications are observed for both the leading and subleading jets. However, while subleading jets show significant modifications for events with a larger dijet momentum imbalance, much smaller modifications are found for the leading jets in these events.

A systematic study of the factorization of long-range azimuthal two-particle correlations into a product of single-particle anisotropies is presented as a function of pt and eta of both particles, and as a function of the particle multiplicity in PbPb and pPb collisions. The data were taken with the CMS detector for PbPb collisions at sqrt(s[NN]) = 2.76 TeV and pPb collisions at sqrt(s[NN]) = 5.02 TeV, covering a very wide range of multiplicity. Factorization is observed to be broken as a function of both particle pt and eta. When measured with particles of different pt, the magnitude of the factorization breakdown for the second Fourier harmonic reaches 20% for very central PbPb collisions but decreases rapidly as the multiplicity decreases. The data are consistent with viscous hydrodynamic predictions, which suggest that the effect of factorization breaking is mainly sensitive to the initial-state conditions rather than to the transport properties (e.g., shear viscosity) of the medium. The factorization breakdown is also computed with particles of different eta. The effect is found to be weakest for mid-central PbPb events but becomes larger for more central or peripheral PbPb collisions, and also for very high-multiplicity pPb collisions. The eta-dependent factorization data provide new insights to the longitudinal evolution of the medium formed in heavy ion collisions.

An invariant differential cross section measurement of inclusive $\pi^{0}$ and $\eta$ meson production at mid-rapidity in pp collisions at $\sqrt{s}=8$ TeV was carried out by the ALICE experiment at the LHC. The spectra of $\pi^{0}$ and $\eta$ mesons were measured in transverse momentum ranges of $0.3<p_{\rm T}<35$ GeV/$c$ and $0.5<p_{\rm T}<35$ GeV/$c$, respectively. Next-to-leading order perturbative QCD calculations using fragmentation functions DSS14 for the $\pi^{0}$ and AESSS for the $\eta$ overestimate the cross sections of both neutral mesons, although such calculations agree with the measured $\eta/\pi^{0}$ ratio within uncertainties. The results were also compared with PYTHIA~8.2 predictions for which the Monash~2013 tune yields the best agreement with the measured neutral meson spectra. The measurements confirm a universal behavior of the $\eta/\pi^{0}$ ratio seen for NA27, PHENIX and ALICE data for pp collisions from $\sqrt{s}=27.5$ GeV to $\sqrt{s}=8$ TeV within experimental uncertainties. A relation between the $\pi^{0}$ and $\eta$ production cross sections for pp collisions at $\sqrt{s}=8$ TeV is given by $m_{\rm T}$ scaling for $p_{\rm T}>3.5$ GeV/$c$. However, a deviation from this empirical scaling rule is observed for transverse momenta below $p_{\rm T}<3.5$ GeV/$c$ in the $\eta/\pi^0$ ratio with a significance of $6.2\sigma$.

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.

We present a measurement of b jet transverse momentum (pt) spectra in proton-lead (pPb) collisions using a dataset corresponding to about 35 inverse nanobarns collected with the CMS detector at the LHC. Jets from b quark fragmentation are found by exploiting the long lifetime of hadrons containing a b quark through tagging methods using distributions of the secondary vertex mass and displacement. Extracted cross sections for b jets are scaled by the effective number of nucleon-nucleon collisions and are compared to a reference obtained from PYTHIA simulations of pp collisions. The PYTHIA-based estimate of the nuclear modification factor is found to be 1.22 +/- 0.15 (stat + syst pPb) +/- 0.27 (syst PYTHIA) averaged over all jets with pt between 55 and 400 GeV/c and with abs(eta[lab]) < 2. We also compare this result to predictions from models using perturbative calculations in quantum chromodynamics.

For the first time a principle-component analysis is used to separate out different orthogonal modes of the two-particle correlation matrix from heavy ion collisions. The analysis uses data from sqrt(s[NN]) = 2.76 TeV PbPb and sqrt(s[NN]) = 5.02 TeV pPb collisions collected by the CMS experiment at the LHC. Two-particle azimuthal correlations have been extensively used to study hydrodynamic flow in heavy ion collisions. Recently it has been shown that the expected factorization of two-particle results into a product of the constituent single-particle anisotropies is broken. The new information provided by these modes may shed light on the breakdown of flow factorization in heavy ion collisions. The first two modes ("leading" and "subleading") of two-particle correlations are presented for elliptical and triangular anisotropies in PbPb and pPb collisions as a function of pt over a wide range of event activity. The leading mode is found to be essentially equivalent to the anisotropy harmonic previously extracted from two-particle correlation methods. The subleading mode represents a new experimental observable and is shown to account for a large fraction of the factorization breaking recently observed at high transverse momentum. The principle-component analysis technique has also been applied to multiplicity fluctuations. These also show a subleading mode. The connection of these new results to previous studies of factorization is discussed.

The measurement of the production of deuterons, tritons and $^{3}\mathrm{He}$ and their antiparticles in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV is presented in this article. The measurements are carried out at midrapidity ($|y| < $ 0.5) as a function of collision centrality using the ALICE detector. The $p_{\rm T}$-integrated yields, the coalescence parameters and the ratios to protons and antiprotons are reported and compared with nucleosynthesis models. The comparison of these results in different collision systems at different centre-of-mass collision energies reveals a suppression of nucleus production in small systems. In the Statistical Hadronisation Model framework, this can be explained by a small correlation volume where the baryon number is conserved, as already shown in previous fluctuation analyses. However, a different size of the correlation volume is required to describe the proton yields in the same data sets. The coalescence model can describe this suppression by the fact that the wave functions of the nuclei are large and the fireball size starts to become comparable and even much smaller than the actual nucleus at low multiplicities.

This article reports measurements of the angle between differently defined jet axes in pp collisions at $\sqrt{s} = 5.02$ TeV carried out by the ALICE Collaboration. Charged particles at midrapidity are clustered into jets with resolution parameters $R=0.2$ and 0.4. The jet axis, before and after Soft Drop grooming, is compared to the jet axis from the Winner-Takes-All (WTA) recombination scheme. The angle between these axes, $\Delta R_{\mathrm{axis}}$, probes a wide phase space of the jet formation and evolution, ranging from the initial high-momentum-transfer scattering to the hadronization process. The $\Delta R_{\mathrm{axis}}$ observable is presented for $20 < {p_{\mathrm{T}}^{\mathrm{ch\; jet}}}< 100$ GeV/$c$, and compared to predictions from the PYTHIA 8 and Herwig 7 event generators. The distributions can also be calculated analytically with a leading hadronization correction related to the non-perturbative component of the Collins$-$Soper$-$Sterman (CSS) evolution kernel. Comparisons to analytical predictions at next-to-leading-logarithmic accuracy with leading hadronization correction implemented from experimental extractions of the CSS kernel in Drell$-$Yan measurements are presented. The analytical predictions describe the measured data within 20% in the perturbative regime, with surprising agreement in the non-perturbative regime as well. These results are compatible with the universality of the CSS kernel in the context of jet substructure.

The transverse momentum ($p_\mathrm{t}$) spectrum of prompt D$^0$ mesons and their antiparticles has been measured via the hadronic decay channels D$^0 \to \mathrm{K}^- \pi^+$ and $\overline{\mathrm{D}}^0 \to \mathrm{K}^+ \pi^-$ in pp and PbPb collisions at a centre-of-mass energy of 5.02 TeV per nucleon pair with the CMS detector at the LHC. The measurement is performed in the D$^0$ meson $p_\mathrm{t}$ range of 2-100 GeV and in the rapidity range of $|y| < $1. The pp (PbPb) dataset used for this analysis corresponds to an integrated luminosity of 27.4 pb$^{-1}$ (530 $\mu$b$^{-1}$). The measured D$^0$ meson $p_\mathrm{t}$ spectrum in pp collisions is well described by perturbative QCD calculations. The nuclear modification factor, comparing D$^0$ meson yields in PbPb and pp collisions, was extracted for both minimum-bias and the 10% most central PbPb interactions. For central events, the D$^0$ meson yield in the PbPb collisions is suppressed by a factor of 5-6 compared to the pp reference in the $p_\mathrm{t}$ range of 6-10 GeV. For D$^0$ mesons in the high-$p_\mathrm{t}$ range of 60-100 GeV, a significantly smaller suppression is observed. The results are also compared to theoretical calculations.

The azimuthal anisotropy Fourier coefficients ($v_n$) in 8.16 TeV pPb data are extracted via long-range two-particle correlations as a function of event multiplicity and compared to corresponding results in pp and PbPb collisions. Using a four-particle cumulant technique, $v_n$ correlations are measured for the first time in pp and pPb collisions. The $v_2$ and $v_4$ coefficients are found to be positively correlated in all collision systems. For high multiplicity pPb collisions an anticorrelation of $v_2$ and $v_3$ is observed, with a similar correlation strength as in PbPb data at the same multiplicity. The new correlation results strengthen the case for a common origin of the collectivity seen in pPb and PbPb collisions in the measured multiplicity range.

Production cross sections of $\Upsilon$(1S), $\Upsilon$(2S), and $\Upsilon$(3S) states decaying into $\mu^+\mu^-$ in proton-lead (pPb) collisions are reported using data collected by the CMS experiment at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV. A comparison is made with corresponding cross sections obtained with pp data measured at the same collision energy and scaled by the Pb nucleus mass number. The nuclear modification factor for $\Upsilon$(1S) is found to be $R_\mathrm{pPb}(\Upsilon(1S))$ = 0.806 $\pm$ 0.024 (stat) $\pm$ 0.059 (syst). Similar results for the excited states indicate a sequential suppression pattern, such that $R_\mathrm{pPb}(\Upsilon(1S))$$\gt$$R_\mathrm{pPb}(\Upsilon(2S))$$\gt$$R_\mathrm{pPb}(\Upsilon(3S))$. The suppression is much less pronounced in pPb than in PbPb collisions, and independent of transverse momentum $p_\mathrm{T}^\Upsilon$ and center-of-mass rapidity $y_\mathrm{CM}^\Upsilon$ of the individual $\Upsilon$ state in the studied range $p_\mathrm{T}^\Upsilon$$\lt$ 30 GeV$/c$ and $\vert y_\mathrm{CM}^\Upsilon\vert$$\lt$ 1.93. Models that incorporate sequential suppression of bottomonia in pPb collisions are in better agreement with the data than those which only assume initial-state modifications.

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.