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

Measurements of the charge-dependent two-particle angular correlation function in proton-lead (pPb) collisions at a nucleon-nucleon center-of-mass energy of $\sqrt{s_\mathrm{NN}}$ = 8.16 TeV and lead-lead (PbPb) collisions at$\sqrt{s_\mathrm{NN}}$ = 5.02 TeV are reported. The pPb and PbPb datasets correspond to integrated luminosities of 186\nbinv and 0.607 nb$^{-1}$, respectively, and were collected using the CMS detector at the CERN LHC. The charge-dependent correlations are characterized by balance functions of same- and opposite-sign particle pairs. The balance functions, which contain information about the creation time of charged particle pairs and the development of collectivity, are studied as functions of relative pseudorapidity ($\Delta \eta$) and relative azimuthal angle ($\Delta \phi$), for various multiplicity and transverse momentum ($p_\mathrm{T}$) intervals. A multiplicity dependence of the balance function is observed in $\Delta \eta$ and $\Delta \phi$ for both systems. The width of the balance functions decreases towards high-multiplicity collisions in the momentum region $\lt$2 GeV, for pPb and PbPb results. No multiplicity dependence is observed at higher transverse momentum. The data are compared with HYDJET, HIJING and AMPT generator predictions, none of which capture completely the multiplicity dependence seen in the data.

The correlations between different moments of two flow amplitudes, extracted with the recently developed asymmetric cumulants, are measured in Pb$-$Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV recorded by the ALICE detector at the CERN Large Hadron Collider. The magnitudes of the measured observables show a dependence on the different moments as well as on the collision centrality, indicating the presence of non-linear response in all even moments up to the eighth. Furthermore, the higher-order asymmetric cumulants show different signatures than the symmetric and lower-order asymmetric cumulants. Comparisons with state-of-the-art event generators using two different parametrizations obtained from Bayesian optimization show differences between data and simulations in many of the studied observables, indicating a need for further tuning of the models behind those event generators. These results provide new and independent constraints on the initial conditions and transport properties of the system created in heavy-ion collisions.

The production of $\pi^\pm$, ${\rm K}^\pm$, and $(\overline{\rm p})$p is measured in pp collisions at $\sqrt{s}=13$ TeV in different topological regions. Particle transverse momentum ($p_{\rm T}$) spectra are measured in the ``toward'', ``transverse'', and ``away'' angular regions defined with respect to the direction of the leading particle in the event. While the toward and away regions contain the fragmentation products of the near-side and away-side jets, respectively, the transverse region is dominated by particles from the Underlying Event (UE). The relative transverse activity classifier, $R_{\rm T}=N_{\rm T}/\langle N_{\rm T}\rangle$, is used to group events according to their UE activity, where $N_{\rm T}$ is the measured charged-particle multiplicity per event in the transverse region and $\langle N_{\rm T}\rangle$ is the mean value over all the analysed events. The first measurements of identified particle $p_{\rm T}$ spectra as a function of $R_{\rm T}$ in the three topological regions are reported. The yield of high transverse momentum particles relative to the $R_{\rm T}$-integrated measurement decreases with increasing $R_{\rm T}$ in both the toward and away regions, indicating that the softer UE dominates particle production as $R_{\rm T}$ increases and validating that $R_{\rm T}$ can be used to control the magnitude of the UE. Conversely, the spectral shapes in the transverse region harden significantly with increasing $R_{\rm T}$. This hardening follows a mass ordering, being more significant for heavier particles. The $p_{\rm T}$-differential particle ratios $({\rm p+\overline{p}})/(\pi^+ +\pi^-)$ and $({\rm K^+ +K^-})/(\pi^+ +\pi^-)$ in the low UE limit $(R_{\rm T}\rightarrow 0)$ approach expectations from Monte Carlo generators such as PYTHIA 8 with Monash 2013 tune and EPOS LHC, where the jet-fragmentation models have been tuned to reproduce ${\rm e^+ e^-}$ results.

This paper presents measurements of charged-hadron spectra obtained in $pp$, $p$+Pb, and Pb+Pb collisions at $\sqrt{s}$ or $\sqrt{s_{_\text{NN}}}=5.02$ TeV, and in Xe+Xe collisions at $\sqrt{s_{_\text{NN}}}=5.44$ TeV. The data recorded by the ATLAS detector at the LHC have total integrated luminosities of 25 pb${}^{-1}$, 28 nb${}^{-1}$, 0.50 nb${}^{-1}$, and 3 $\mu$b${}^{-1}$, respectively. The nuclear modification factors $R_{p\text{Pb}}$ and $R_\text{AA}$ are obtained by comparing the spectra in heavy-ion and $pp$ collisions in a wide range of charged-particle transverse momenta and pseudorapidity. The nuclear modification factor $R_{p\text{Pb}}$ shows a moderate enhancement above unity with a maximum at $p_{\mathrm{T}} \approx 3$ GeV; the enhancement is stronger in the Pb-going direction. The nuclear modification factors in both Pb+Pb and Xe+Xe collisions feature a significant, centrality-dependent suppression. They show a similar distinct $p_{\mathrm{T}}$-dependence with a local maximum at $p_{\mathrm{T}} \approx 2$ GeV and a local minimum at $p_{\mathrm{T}} \approx 7$ GeV. This dependence is more distinguishable in more central collisions. No significant $|\eta|$-dependence is found. A comprehensive comparison with several theoretical predictions is also provided. They typically describe $R_\text{AA}$ better in central collisions and in the $p_{\mathrm{T}}$ range from about 10 to 100 GeV.

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.

We report the measurement of $K^{*0}$ meson at midrapidity ($|y|<$ 1.0) in Au+Au collisions at $\sqrt{s_{\rm NN}}$~=~7.7, 11.5, 14.5, 19.6, 27 and 39 GeV collected by the STAR experiment during the RHIC beam energy scan (BES) program. The transverse momentum spectra, yield, and average transverse momentum of $K^{*0}$ are presented as functions of collision centrality and beam energy. The $K^{*0}/K$ yield ratios are presented for different collision centrality intervals and beam energies. The $K^{*0}/K$ ratio in heavy-ion collisions are observed to be smaller than that in small system collisions (e+e and p+p). The $K^{*0}/K$ ratio follows a similar centrality dependence to that observed in previous RHIC and LHC measurements. The data favor the scenario of the dominance of hadronic re-scattering over regeneration for $K^{*0}$ production in the hadronic phase of the medium.

We report the triton ($t$) production in mid-rapidity ($|y| <$ 0.5) Au+Au collisions at $\sqrt{s_\mathrm{NN}}$= 7.7--200 GeV measured by the STAR experiment from the first phase of the beam energy scan at the Relativistic Heavy Ion Collider (RHIC). The nuclear compound yield ratio ($\mathrm{N}_t \times \mathrm{N}_p/\mathrm{N}_d^2$), which is predicted to be sensitive to the fluctuation of local neutron density, is observed to decrease monotonically with increasing charged-particle multiplicity ($dN_{ch}/d\eta$) and follows a scaling behavior. The $dN_{ch}/d\eta$ dependence of the yield ratio is compared to calculations from coalescence and thermal models. Enhancements in the yield ratios relative to the coalescence baseline are observed in the 0%-10% most central collisions at 19.6 and 27 GeV, with a significance of 2.3$\sigma$ and 3.4$\sigma$, respectively, giving a combined significance of 4.1$\sigma$. The enhancements are not observed in peripheral collisions or model calculations without critical fluctuation, and decreases with a smaller $p_{T}$ acceptance. The physics implications of these results on the QCD phase structure and the production mechanism of light nuclei in heavy-ion collisions are discussed.

The first measurements of elliptic flow of $\pi^\pm$, ${\rm K}^\pm$, p+$\overline{\rm p}$, ${\rm K_{S}^0}$, $\Lambda$+$\overline{\Lambda}$, $\phi$, $\Xi^-$+$\Xi^+$, and $\Omega^-$+$\Omega^+$ using multiparticle cumulants in Pb$-$Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV are presented. Results obtained with two- ($v_2\{2\}$) and four-particle cumulants ($v_2\{4\}$) are shown as a function of transverse momentum, $p_{\rm T}$, for various collision centrality intervals. Combining the data for both $v_2\{2\}$ and $v_2\{4\}$ also allows us to report the first measurements of the mean elliptic flow, elliptic flow fluctuations, and relative elliptic flow fluctuations for various hadron species. These observables probe the event-by-event eccentricity fluctuations in the initial state and the contributions from the dynamic evolution of the expanding quark-gluon plasma. The characteristic features observed in previous $p_{\rm T}$-differential anisotropic flow measurements for identified hadrons with two-particle correlations, namely the mass ordering at low $p_{\rm T}$ and the approximate scaling with the number of constituent quarks at intermediate $p_{\rm T}$, are similarly present in the four-particle correlations and the combinations of $v_2\{2\}$ and $v_2\{4\}$. In addition, a particle species dependence of flow fluctuations is observed that could indicate a significant contribution from final state hadronic interactions. The comparison between experimental measurements and CoLBT model calculations, which combine the various physics processes of hydrodynamics, quark coalescence, and jet fragmentation, illustrates their importance over a wide $p_{\rm T}$ range.

Azimuthal anisotropy of produced particles is one of the most important observables used to access the collective properties of the expanding medium created in relativistic heavy-ion collisions. In this paper, we present second ($v_{2}$) and third ($v_{3}$) order azimuthal anisotropies of $K_{S}^{0}$, $\phi$, $\Lambda$, $\Xi$ and $\Omega$ at mid-rapidity ($|y|<$1) in Au+Au collisions at $\sqrt{s_{\text{NN}}}$ = 54.4 GeV measured by the STAR detector. The $v_{2}$ and $v_{3}$ are measured as a function of transverse momentum and centrality. Their energy dependence is also studied. $v_{3}$ is found to be more sensitive to the change in the center-of-mass energy than $v_{2}$. Scaling by constituent quark number is found to hold for $v_{2}$ within 10%. This observation could be evidence for the development of partonic collectivity in 54.4 GeV Au+Au collisions. Differences in $v_{2}$ and $v_{3}$ between baryons and anti-baryons are presented, and ratios of $v_{3}$/$v_{2}^{3/2}$ are studied and motivated by hydrodynamical calculations. The ratio of $v_{2}$ of $\phi$ mesons to that of anti-protons ($v_{2}(\phi)/v_{2}(\bar{p})$) shows centrality dependence at low transverse momentum, presumably resulting from the larger effects from hadronic interactions on anti-proton $v_{2}$.

The correlations between flow harmonics $v_n$ for $n=2$, 3 and 4 and mean transverse momentum $[p_\mathrm{T}]$ in $^{129}$Xe+$^{129}$Xe and $^{208}$Pb+$^{208}$Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.44$ TeV and 5.02 TeV, respectively, are measured using charged particles with the ATLAS detector. The correlations are sensitive to the shape and size of the initial geometry, nuclear deformation, and initial momentum anisotropy. The effects from non-flow and centrality fluctuations are minimized, respectively, via a subevent cumulant method and event activity selection based on particle production in the very forward rapidity. The results show strong dependences on centrality, harmonic number $n$, $p_{\mathrm{T}}$ and pseudorapidity range. Current models describe qualitatively the overall centrality- and system-dependent trends but fail to quantitatively reproduce all the data. In the central collisions, where models generally show good agreement, the $v_2$-$[p_\mathrm{T}]$ correlations are sensitive to the triaxiality of the quadruple deformation. The comparison of model to the Pb+Pb and Xe+Xe data suggests that the $^{129}$Xe nucleus is a highly deformed triaxial ellipsoid that is neither a prolate nor an oblate shape. This provides strong evidence for a triaxial deformation of $^{129}$Xe nucleus using high-energy heavy-ion collision.

The transverse-momentum ($p_{\rm T}$) spectra of ${\rm K}^*(892)^0$ and $\phi(1020)$ measured with the ALICE detector up to $p_{\rm T}$ = 16 GeV/$c$ in the rapidity range $-1.2 < y < 0.3$, in p-Pb collisions at the center-of-mass energy per nucleon-nucleon collision $\sqrt{s_{\rm NN}}$ = 5.02 TeV are presented as a function of charged particle multiplicity and rapidity. The measured $p_{\rm T}$ distributions show a dependence on both multiplicity and rapidity at low $p_{\rm T}$ whereas no significant dependence is observed at high $p_{\rm T}$. A rapidity dependence is observed in the $p_{\rm T}$-integrated yield (d$N$/d$y$), whereas the mean transverse momentum ($\left< p_{\rm T} \right>$) shows a flat behavior as a function of rapidity. The rapidity asymmetry ($Y_{\rm asym}$) at low $p_{\rm T}$ ( < 5 GeV/$c$) is more significant for higher multiplicity classes. At high $p_{\rm T}$, no significant rapidity asymmetry is observed in any of the multiplicity classes. Both ${\rm K}^*(892)^0$ and $\phi(1020)$ show similar $Y_{\rm asym}$. The nuclear modification factor ($Q_{\rm CP}$) as a function of $p_{\rm T}$ shows a Cronin-like enhancement at intermediate $p_{\rm T}$, which is more prominent at higher rapidities (Pb-going direction) and in higher multiplicity classes. At high $p_{\rm T}$ (> 5 GeV/$c$), the $Q_{\rm CP}$ values are greater than unity and no significant rapidity dependence is observed.

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.

The production of $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $\mathrm{\phi(1020)}$ resonances has been measured in p-Pb collisions at $\sqrt{s_{\rm NN}}$ = 8.16 TeV using the ALICE detector. Resonances are reconstructed via their hadronic decay channels in the rapidity interval $-$0.5 $<$$y$$<$ 0 and the transverse momentum spectra are measured for various multiplicity classes up to $p_{\rm T}$ = 20 GeV/$c$ for $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $p_{\rm T}$ = 16 GeV/$c$ for $\mathrm{\phi(1020)}$. The $p_{\rm T}$ -integrated yields and mean transverse momenta are reported and compared with previous results in pp, p-Pb and Pb-Pb collisions. The $x_{\mathrm{T}}$ scaling for $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $\mathrm{\phi(1020)}$ resonance production is newly tested in p-Pb collisions and found to hold in the high-$p_{\rm T}$ region at LHC energies. The nuclear modification factors ($R_{\rm pPb}$) as a function of $p_{\rm T}$ for $\mathrm{K}^{*0}$ and $\mathrm{\phi}$ at $\sqrt{s_{NN}}$ = 8.16 TeV are presented along with the new $R_{\rm pPb}$ measurements of $\mathrm{K}^{*0}$, $\mathrm{\phi}$ , $\Xi$, and $\Omega$ at $\sqrt{s_{\rm NN}}$ = 5.02 TeV. At intermediate $p_{\rm T}$ (2-8 GeV/$c$), $R_{\rm pPb}$ of $\Xi$, $\Omega$ show a Cronin-like enhancement, while $\mathrm{K}^{*0}$ and $\mathrm{\phi}$ show no or little nuclear modification. At high $p_{\rm T}$ ($>$ 8 GeV/$c$), the $R_{\rm pPb}$ values of all hadrons are consistent with unity within uncertainties. The $R_{\rm pPb}$ of $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $\mathrm{\phi(1020)}$ at $\sqrt{s_{\rm NN}}$ = 8.16 and 5.02 TeV show no significant energy dependence.

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 chiral magnetic effect (CME) is predicted to occur as a consequence of a local violation of $\cal P$ and $\cal CP$ symmetries of the strong interaction amidst a strong electro-magnetic field generated in relativistic heavy-ion collisions. Experimental manifestation of the CME involves a separation of positively and negatively charged hadrons along the direction of the magnetic field. Previous measurements of the CME-sensitive charge-separation observables remain inconclusive because of large background contributions. In order to better control the influence of signal and backgrounds, the STAR Collaboration performed a blind analysis of a large data sample of approximately 3.8 billion isobar collisions of $^{96}_{44}$Ru+$^{96}_{44}$Ru and $^{96}_{40}$Zr+$^{96}_{40}$Zr at $\sqrt{s_{\rm NN}}=200$ GeV. Prior to the blind analysis, the CME signatures are predefined as a significant excess of the CME-sensitive observables in Ru+Ru collisions over those in Zr+Zr collisions, owing to a larger magnetic field in the former. A precision down to 0.4% is achieved, as anticipated, in the relative magnitudes of the pertinent observables between the two isobar systems. Observed differences in the multiplicity and flow harmonics at the matching centrality indicate that the magnitude of the CME background is different between the two species. No CME signature that satisfies the predefined criteria has been observed in isobar collisions in this blind analysis.

The jet angularities are a class of jet substructure observables which characterize the angular and momentum distribution of particles within jets. These observables are sensitive to momentum scales ranging from perturbative hard scatterings to nonperturbative fragmentation into final-state hadrons. We report measurements of several groomed and ungroomed jet angularities in pp collisions at $\sqrt{s}=5.02$ TeV with the ALICE detector. Jets are reconstructed using charged particle tracks at midrapidity ($|\eta| < 0.9$). The anti-$k_{\rm T}$ algorithm is used with jet resolution parameters $R=0.2$ and $R=0.4$ for several transverse momentum $p_{\rm T}^{\text{ch jet}}$ intervals in the 20$-$100 GeV/$c$ range. Using the jet grooming algorithm Soft Drop, the sensitivity to softer, wide-angle processes, as well as the underlying event, can be reduced in a way which is well-controlled in theoretical calculations. We report the ungroomed jet angularities, $\lambda_{\alpha}$, and groomed jet angularities, $\lambda_{\alpha\text{,g}}$, to investigate the interplay between perturbative and nonperturbative effects at low jet momenta. Various angular exponent parameters $\alpha = 1$, 1.5, 2, and 3 are used to systematically vary the sensitivity of the observable to collinear and soft radiation. Results are compared to analytical predictions at next-to-leading-logarithmic accuracy, which provide a generally good description of the data in the perturbative regime but exhibit discrepancies in the nonperturbative regime. Moreover, these measurements serve as a baseline for future ones in heavy-ion collisions by providing new insight into the interplay between perturbative and nonperturbative effects in the angular and momentum substructure of jets. They supply crucial guidance on the selection of jet resolution parameter, jet transverse momentum, and angular scaling variable for jet quenching studies.

Measurements of elliptic ($v_2$) and triangular ($v_3$) flow coefficients of $\pi^{\pm}$, K$^{\pm}$, p+$\rm \overline{p}$, K$^0_{\rm S}$, and $\Lambda + \overline{\Lambda}$ obtained with the scalar product method in Xe-Xe collisions at $\sqrt{s_{\rm NN}}$ = 5.44 TeV are presented. The results are obtained in the rapidity range $\left | y \right |<0.5$ and reported as a function of transverse momentum, $p_{\rm T}$, for several collision centrality classes. The flow coefficients exhibit a particle mass dependence for $p_{\rm T}<3$ GeV/$c$, while a grouping according to particle type (i.e., meson and baryon) is found at intermediate transverse momenta (3< $p_{\rm T}$ <8 GeV/$c$). The magnitude of the baryon $v_{2}$ is larger than that of mesons up to $p_{\rm T}$ = 6 GeV/$c$. The centrality dependence of the shape evolution of the $p_{\rm T}$-differential $v_2$ is studied for the various hadron species. The $v_2$ coefficients of $\pi^{\pm}$, K$^{\pm}$, and p+$\rm \overline{p}$ are reproduced by MUSIC hydrodynamic calculations coupled to a hadronic cascade model (UrQMD) for $p_{\rm T} <1$ GeV/$c$. A comparison with $v_{\rm n}$ measurements in the corresponding centrality intervals in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV yields an enhanced $v_2$ in central collisions and diminished value in semicentral collisions.

The production of K$^{*}(892)^{0}$ and $\phi(1020)$ mesons in proton-proton (pp) and lead-lead (Pb-Pb) collisions at $\sqrt{s_\mathrm{NN}} = 5.02$ TeV has been measured using the ALICE detector at the Large Hadron Collider (LHC). The transverse momentum ($p_{\mathrm{T}}$) distributions of K$^{*}(892)^{0}$ and $\phi(1020)$ mesons have been measured at midrapidity $(|y|<0.5)$ up to $p_{\mathrm{T}} = 20$ GeV$/c$ in inelastic pp collisions and for several Pb-Pb collision centralities. The collision centrality and collision energy dependence of the average transverse momenta agree with the radial flow scenario observed with stable hadrons, showing that the effect is stronger for more central collisions and higher collision energies. The $\mathrm{K^{*0}/K}$ ratio is found to be suppressed in Pb-Pb collisions relative to pp collisions: this indicates a loss of the measured K$^{*}(892)^{0}$ signal due to rescattering of its decay products in the hadronic phase. In contrast, for the longer-lived $\phi(1020)$ mesons, no such suppression is observed. The nuclear modification factors ($R_{\rm AA}$) of K$^{*}(892)^{0}$ and $\phi(1020)$ mesons are calculated using pp reference spectra at the same collision energy. In central Pb-Pb collisions for $p_{\rm T} > 8$ GeV$/c$, the $R_{\rm AA}$ values of K$^{*}(892)^{0}$ and $\phi(1020)$ are below unity and observed to be similar to those of pions, kaons, and (anti)protons. The $R_{\rm AA}$ values at high $p_{\mathrm T}$ ($>$~8 GeV$/c$) for K$^{*}(892)^{0}$ and $\phi(1020)$ mesons are in agreement within uncertainties for $\sqrt{s_\mathrm{NN}} = 5.02$ and 2.76 TeV.

The spectra measurements of charged hadrons produced in the shower of a parton originating in the same hard scattering with a leptonically decaying Z boson, are reported in lead-lead (PbPb) and proton-proton (pp) collisions at a nucleon-nucleon center-of-mass energy of 5.02 TeV. Both PbPb and pp data sets are recorded by the CMS experiment at the LHC, and correspond to an integrated luminosity of 1.7 nb$^{-1}$ and 320 pb$^{-1}$, respectively. Hadronic collision data with one reconstructed Z boson candidate with the transverse momentum $p_\mathrm{T}$$\gt$ 30 GeV/$c$ are analyzed. The Z boson constrains the initial energy and direction of the associated parton. In heavy ion events, azimuthal angular distributions of charged hadrons with respect to the direction of a Z boson are sensitive to modifications of the in-medium parton shower and medium response. Compared to reference data from pp interactions, the results for central PbPb collisions indicate a modification of the angular correlations. The measurements of the fragmentation functions and $p_\mathrm{T}$ spectra of charged particles in Z boson events, which are sensitive to medium modifications of the parton shower longitudinal structure, are also reported. Significant modifications in central PbPb events compared to pp reference data are also found for these observables.

Differential cross sections for the Drell-Yan process, including Z boson production, using the dimuon decay channel are measured in proton-lead (pPb) collisions at a nucleon-nucleon centre-of-mass energy of 8.16 TeV. A data sample recorded with the CMS detector at the LHC is used, corresponding to an integrated luminosity of 173 nb$^{-1}$. The differential cross section as a function of the dimuon mass is measured in the range 15-600 GeV, for the first time in proton-nucleus collisions. It is also reported as a function of dimuon rapidity over the mass ranges 15-60 GeV and 60-120 GeV, and ratios for the p-going over the Pb-going beam directions are built. In both mass ranges, the differential cross sections as functions of the dimuon transverse momentum $p_\mathrm{T}$ and of a geometric variable $\phi^*$ are measured, where $\phi^*$ highly correlates with $p_\mathrm{T}$ but is determined with higher precision. In the Z mass region, the rapidity dependence of the data indicate a modification of the distribution of partons within a lead nucleus as compared to the proton case. The data are more precise than predictions based upon current models of parton distributions.

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.

The first measurement of the production of pions, kaons, (anti-)protons and $\phi$ mesons at midrapidity in Xe-Xe collisions at $\sqrt{s_{\rm NN}} = 5.44$ TeV is presented. Transverse momentum ($p_{\rm T}$) spectra and $p_{\rm T}$-integrated yields are extracted in several centrality intervals bridging from p-Pb to mid-central Pb-Pb collisions in terms of final-state multiplicity. The study of Xe-Xe and Pb-Pb collisions allows systems at similar charged-particle multiplicities but with different initial geometrical eccentricities to be investigated. A detailed comparison of the spectral shapes in the two systems reveals an opposite behaviour for radial and elliptic flow. In particular, this study shows that the radial flow does not depend on the colliding system when compared at similar charged-particle multiplicity. In terms of hadron chemistry, the previously observed smooth evolution of particle ratios with multiplicity from small to large collision systems is also found to hold in Xe-Xe. In addition, our results confirm that two remarkable features of particle production at LHC energies are also valid in the collision of medium-sized nuclei: the lower proton-to-pion ratio with respect to the thermal model expectations and the increase of the $\phi$-to-pion ratio with increasing final-state multiplicity.

The NA61/SHINE experiment at the CERN Super Proton Synchrotron (SPS) studies the onset of deconfinement in hadron matter by a scan of particle production in collisions of nuclei with various sizes at a set of energies covering the SPS energy range. This paper presents results on inclusive double-differential spectra, transverse momentum and rapidity distributions and mean multiplicities of $\pi^\pm$, $K^\pm$, $p$ and $\bar{p}$ produced in the 20$\%$ most $central$$^7$Be+$^9$Be collisions at beam momenta of 19$A$, 30$A$, 40$A$, 75$A$ and 150$A$ GeV/$c$. The energy dependence of the $K^\pm$/$\pi^\pm$ ratios as well as of inverse slope parameters of the $K^\pm$ transverse mass distributions are close to those found in inelastic $p$+$p$ reactions. The new results are compared to the world data on $p$+$p$ and Pb+Pb collisions as well as to predictions of the EPOS, UrQMD, AMPT, PHSD and SMASH models.

The inclusive production of the J/$\psi$ and $\psi$(2S) charmonium states is studied as a function of centrality in p-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}} = 8.16$ TeV at the LHC. The measurement is performed in the dimuon decay channel with the ALICE apparatus in the centre-of-mass rapidity intervals $-4.46 < y_{\rm cms} < -2.96$ (Pb-going direction) and $2.03 < y_{\rm cms} < 3.53$ (p-going direction), down to zero transverse momentum ($p_{\rm T}$). The J/$\psi$ and $\psi$(2S) production cross sections are evaluated as a function of the collision centrality, estimated through the energy deposited in the zero degree calorimeter located in the Pb-going direction. The $p_{\rm T}$-differential J/$\psi$ production cross section is measured at backward and forward rapidity for several centrality classes, together with the corresponding average $\langle p_{\rm T} \rangle$ and $\langle p^{2}_{\rm T} \rangle$ values. The nuclear effects affecting the production of both charmonium states are studied using the nuclear modification factor. In the p-going direction, a suppression of the production of both charmonium states is observed, which seems to increase from peripheral to central collisions. In the Pb-going direction, however, the centrality dependence is different for the two states: the nuclear modification factor of the J/$\psi$ increases from below unity in peripheral collisions to above unity in central collisions, while for the $\psi$(2S) it stays below or consistent with unity for all centralities with no significant centrality dependence. The results are compared with measurements in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV and no significant dependence on the energy of the collision is observed. Finally, the results are compared with theoretical models implementing various nuclear matter effects.