The collective behavior of K$^0_\mathrm{S}$ and $\Lambda/\bar{\Lambda}$ strange hadrons is studied by measuring the elliptic azimuthal anisotropy ($v_2$) using the scalar-product and multiparticle correlation methods. Proton-lead (pPb) collisions at a nucleon-nucleon center-of-mass energy $\sqrt{s_\mathrm{NN}}$ = 8.16 TeV and lead-lead (PbPb) collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV collected by the CMS experiment at the LHC are investigated. Nonflow effects in the pPb collisions are studied by using a subevent cumulant analysis and by excluding events where a jet with transverse momentum greater than 20\GeV is present. The strange hadron $v_2$ values extracted in \pPb collisions via the four- and six-particle correlation method are found to be nearly identical, suggesting the collective behavior. Comparisons of the pPb and PbPb results for both strange hadrons and charged particles illustrate how event-by-event flow fluctuations depend on the system size.

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 production of inclusive, prompt and non-prompt J/$\psi$ was studied for the first time at midrapidity ($ -1.37 < y_{\rm cms} < 0.43$) in p$-$Pb collisions at $\sqrt{s_{\rm NN}} = 8.16$ TeV with the ALICE detector at the LHC. The inclusive J/$\psi$ mesons were reconstructed in the dielectron decay channel in the transverse momentum ($p_{\rm T}$) interval $0 < p_{\rm T} < 14$ GeV/$c$ and the prompt and non-prompt contributions were separated on a statistical basis for $p_{\rm T} > 2$ GeV/$c$. The study of the J/$\psi$ mesons in the dielectron channel used for the first time in ALICE online single-electron triggers from the Transition Radiation Detector, providing a data sample corresponding to an integrated luminosity of $689 \pm 13 \mu{\rm b}^{-1}$. The proton$-$proton reference cross section for inclusive J/$\psi$ was obtained based on interpolations of measured data at different centre-of-mass energies and a universal function describing the $p_{\rm T}$-differential J/$\psi$ production cross sections. The $p_{\rm T}$-differential nuclear modification factors $R_{\rm pPb}$ of inclusive, prompt, and non-prompt J/$\psi$ are consistent with unity and described by theoretical models implementing only nuclear shadowing.

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

We present the charged-particle multiplicity distributions over a wide pseudorapidity range ($-3.4<\eta<5.0$) for pp collisions at $\sqrt{s}=$ 0.9, 7, and 8 TeV at the LHC. Results are based on information from the Silicon Pixel Detector and the Forward Multiplicity Detector of ALICE, extending the pseudorapidity coverage of the earlier publications and the high-multiplicity reach. The measurements are compared to results from the CMS experiment and to PYTHIA, PHOJET and EPOS LHC event generators, as well as IP-Glasma calculations.

The second- and third-order azimuthal anisotropy Fourier harmonics of charged particles produced in pPb collisions, at $\sqrt{s_\mathrm{NN}} =$ 8.16 TeV, are studied over a wide range of event multiplicities. Multiparticle correlations are used to isolate global properties stemming from the collision overlap geometry. The second-order "elliptic" harmonic moment is obtained with high precision through four-, six-, and eight-particle correlations and, for the first time, the third-order "triangular" harmonic moment is studied using four-particle correlations. A sample of peripheral PbPb collisions at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV that covers a similar range of event multiplicities as the pPb results is also analyzed. Model calculations of initial-state fluctuations in pPb and PbPb collisions can be directly compared to the high precision experimental results. This work provides new insight into the fluctuation-driven origin of the $v_3$ coefficients in pPb and PbPb collisions, and into the dominating overall collision geometry in PbPb collisions at the earliest stages of heavy ion interactions.

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.

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.