The production of $\Upsilon (nS)$ mesons ($n=1,2,3$) in $p$Pb and Pb$p$ collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{NN}}=8.16$ TeV is measured by the LHCb experiment, using a data sample corresponding to an integrated luminosity of 31.8 nb$^{-1}$. The $\Upsilon (nS)$ mesons are reconstructed through their decays into two opposite-sign muons. The measurements comprise the differential production cross-sections of the $\Upsilon (1S)$ and $\Upsilon (2S)$ states, their forward-to-backward ratios and nuclear modification factors, performed as a function of the transverse momentum $p_{\mathrm{T}}$ and rapidity in the nucleon-nucleon centre-of-mass frame $y^*$ of the $\Upsilon (nS)$ states, in the kinematic range $p_{\rm{T}}<25$ GeV/$c$ and $1.5<y^*<4.0$ ($-5.0<y^*<-2.5$) for $p$Pb (Pb$p$) collisions. In addition, production cross-sections for $\Upsilon (3S)$ are measured integrated over phase space and the production ratios between all three $\Upsilon (nS)$ states are determined. The measurements are compared to theoretical predictions and suppressions for quarkonium in $p$Pb collisions are observed.

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 the W$^\pm$ bosons measured in p$-$Pb collisions at a centre-of-mass energy per nucleon$-$nucleon collision $\sqrt{s_{NN}} = 8.16$ TeV and Pb$-$Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV with ALICE at the LHC is presented. The W$^\pm$ bosons are measured via their muonic decay channel, with the muon reconstructed in the pseudorapidity region $-4 < \eta^\mu_{\rm lab} < -2.5$ with transverse momentum $p_{\rm T}^\mu > 10$ GeV/$c$. While in Pb$-$Pb collisions the measurements are performed in the forward ($2.5 < y^\mu_{\rm cms} < 4$) rapidity region, in p$-$Pb collisions, where the centre-of-mass frame is boosted with respect to the laboratory frame, the measurements are performed in the backward ($-4.46 < y^\mu_{\rm cms} < -2.96$) and forward ($2.03 < y^\mu_{\rm cms} < 3.53$) rapidity regions. The W$^{-}$ and W$^{+}$ production cross sections, lepton-charge asymmetry, and nuclear modification factors are evaluated as a function of the muon rapidity. In order to study the production as a function of the p$-$Pb collision centrality, the production cross sections of the W$^{-}$ and W$^{+}$ bosons are combined and normalised to the average number of binary nucleon$-$nucleon collision $\langle N_\mathrm{coll} \rangle$. In Pb$-$Pb collisions, the same measurements are presented as a function of the collision centrality. Study of the binary scaling of the W$^\pm$-boson cross sections in p$-$Pb and Pb$-$Pb collisions is also reported. The results are compared with perturbative QCD (pQCD) calculations, with and without nuclear modifications of the Parton Distribution Functions (PDFs), as well as with available data at the LHC. Significant deviations from the theory expectations are found in the two collision systems, indicating that the measurements can provide additional constraints for the determination of nuclear PDF (nPDFs) and in particular of the light-quark distributions.

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.

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.

We report on the inclusive production cross sections of J/$\psi$, $\psi$(2S), $\Upsilon$(1S), $\Upsilon$(2S) and $\Upsilon$(3S), measured at forward rapidity with the ALICE detector in pp collisions at a center-of-mass energy $\sqrt{s}=8$ TeV. The analysis is based on data collected at the LHC and corresponds to an integrated luminosity of 1.28 pb$^{-1}$. Quarkonia are reconstructed in the dimuon-decay channel. The differential production cross sections are measured as a function of the transverse momentum $p_{\rm T}$ and rapidity $y$, over the $p_{\rm T}$ ranges $0<p_{\rm T}<20$ GeV/$c$ for J/$\psi$, $0<p_{\rm T}<12$ GeV/$c$ for all other resonances, and for $2.5<y<4$. The cross sections, integrated over $p_{\rm T}$ and $y$, and assuming unpolarized quarkonia, are $\sigma_{{\rm J}/\psi} = 8.98\pm0.04\pm0.82$ $\mu$b, $\sigma_{\psi{\rm (2S)}} = 1.23\pm0.08\pm0.22$ $\mu$b, $\sigma_{\Upsilon{\rm(1S)}} = 71\pm6\pm7$ nb, $\sigma_{\Upsilon{\rm(2S)}} = 26\pm5\pm4$ nb and $\sigma_{\Upsilon{\rm(3S)}} = 9\pm4\pm1$ nb, where the first uncertainty is statistical and the second one is systematic. These values agree, within at most $1.4\sigma$, with measurements performed by the LHCb collaboration in the same rapidity range.

A detailed study of pseudorapidity densities and multiplicity distributions of primary charged particles produced in proton-proton collisions, at $\sqrt{s} =$ 0.9, 2.36, 2.76, 7 and 8 TeV, in the pseudorapidity range $|\eta|<2$, was carried out using the ALICE detector. Measurements were obtained for three event classes: inelastic, non-single diffractive and events with at least one charged particle in the pseudorapidity interval $|\eta|<1$. The use of an improved track-counting algorithm combined with ALICE's measurements of diffractive processes allows a higher precision compared to our previous publications. A KNO scaling study was performed in the pseudorapidity intervals $|\eta|<$ 0.5, 1.0 and 1.5. The data are compared to other experimental results and to models as implemented in Monte Carlo event generators PHOJET and recent tunes of PYTHIA6, PYTHIA8 and EPOS.

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