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.

The multi-strange baryon yields in Pb--Pb collisions have been shown to exhibit an enhancement relative to pp reactions. In this work, $\Xi$ and $\Omega$ production rates have been measured with the ALICE experiment as a function of transverse momentum, ${p_{\rm T}}$, in p-Pb collisions at a centre-of-mass energy of ${\sqrt{s_{\rm NN}}}$ = 5.02 TeV. The results cover the kinematic ranges 0.6 GeV/$c<{p_{\rm T}} <$7.2 GeV/$c$ and 0.8 GeV/$c<{p_{\rm T}}<$ 5 GeV/$c$, for $\Xi$ and $\Omega$ respectively, in the common rapidity interval -0.5 $<{y_{\rm CMS}}<$ 0. Multi-strange baryons have been identified by reconstructing their weak decays into charged particles. The ${p_{\rm T}}$ spectra are analysed as a function of event charged-particle multiplicity, which in p-Pb collisions ranges over one order of magnitude and lies between those observed in pp and Pb-Pb collisions. The measured ${p_{\rm T}}$ distributions are compared to the expectations from a Blast-Wave model. The parameters which describe the production of lighter hadron species also describe the hyperon spectra in high multiplicity p-Pb. The yield of hyperons relative to charged pions is studied and compared with results from pp and Pb-Pb collisions. A statistical model is employed, which describes the change in the ratios with volume using a canonical suppression mechanism, in which the small volume causes a species-dependent relative reduction of hadron production. The calculations, in which the magnitude of the effect depends on the strangeness content, show good qualitative agreement with the data.

We report on two-particle charge-dependent correlations in pp, p-Pb, and Pb-Pb collisions as a function of the pseudorapidity and azimuthal angle difference, $\mathrm{\Delta}\eta$ and $\mathrm{\Delta}\varphi$ respectively. These correlations are studied using the balance function that probes the charge creation time and the development of collectivity in the produced system. The dependence of the balance function on the event multiplicity as well as on the trigger and associated particle transverse momentum ($p_{\mathrm{T}}$) in pp, p-Pb, and Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 7$, 5.02, and 2.76 TeV, respectively, are presented. In the low transverse momentum region, for $0.2 < p_{\mathrm{T}} < 2.0$ GeV/$c$, the balance function becomes narrower in both $\mathrm{\Delta}\eta$ and $\mathrm{\Delta}\varphi$ directions in all three systems for events with higher multiplicity. The experimental findings favor models that either incorporate some collective behavior (e.g. AMPT) or different mechanisms that lead to effects that resemble collective behavior (e.g. PYTHIA8 with color reconnection). For higher values of transverse momenta the balance function becomes even narrower but exhibits no multiplicity dependence, indicating that the observed narrowing with increasing multiplicity at low $p_{\mathrm{T}}$ is a feature of bulk particle production.

Measurements of charged jet production as a function of centrality are presented for p-Pb collisions recorded at $\sqrt{s_{\rm NN}} = 5.02$ TeV with the ALICE detector. Centrality classes are determined via the energy deposit in neutron calorimeters at zero degree, close to the beam direction, to minimise dynamical biases of the selection. The corresponding number of participants or binary nucleon-nucleon collisions is determined based on the particle production in the Pb-going rapidity region. Jets have been reconstructed in the central rapidity region from charged particles with the anti-$k_{\rm T}$ algorithm for resolution parameters $R = 0.2$ and $R = 0.4$ in the transverse momentum range 20 to 120 GeV/$c$. The reconstructed jet momentum and yields have been corrected for detector effects and underlying-event background. In the five centrality bins considered, the charged jet production in p-Pb collisions is consistent with the production expected from binary scaling from pp collisions. The ratio of jet yields reconstructed with the two different resolution parameters is also independent of the centrality selection, demonstrating the absence of major modifications of the radial jet structure in the reported centrality classes.

We present the first azimuthally differential measurements of the pion source size relative to the second harmonic event plane in Pb-Pb collisions at a center-of-mass energy per nucleon-nucleon pair of $\sqrt{s_{\rm NN}}=2.76$ TeV. The measurements have been performed in the centrality range 0-50% and for pion pair transverse momenta $0.2 < k_{\rm T} < 0.7$ GeV/$c$. We find that the $R_{\rm side}$ and $R_{\rm out}$ radii, which characterize the pion source size in the directions perpendicular and parallel to the pion transverse momentum, oscillate out of phase, similar to what was observed at the Relativistic Heavy Ion Collider (RHIC). The final-state source eccentricity, estimated via $R_{\rm side}$ oscillations, is found to be significantly smaller than the initial-state source eccentricity, but remains positive; indicating that even after a stronger expansion in the in-plane direction, the pion source at the freeze-out is still elongated in the out-of-plane direction. The 3+1D hydrodynamic calculations are in qualitative agreement with observed centrality and transverse momentum $R_{\rm side}$ oscillations, but systematically underestimate the oscillation magnitude.

The azimuthal correlations of D mesons and charged particles were measured with the ALICE detector in pp collisions at $\sqrt{s}=7$ TeV and p-Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV at the Large Hadron Collider. D$^0$, D$^+$, and D$^{*+}$ mesons and their charge conjugates with transverse momentum $3<p_{\rm T}<16$ GeV/$c$ and rapidity in the nucleon-nucleon centre-of-mass system $|y_{\rm cms}|<0.5$ (pp collisions) and $-0.96<y_{\rm cms}<0.04$ (p-Pb collisions) were correlated to charged particles with $p_{\rm T}>0.3$ Gev/$c$. The properties of the correlation peak induced by the jet containing the D meson, described in terms of the yield of charged particles in the peak and peak width, are compatible within uncertainties between the two collision systems, and described by Monte-Carlo simulations based on the PYTHIA, POWHEG and EPOS 3 event generators.

The elliptic, triangular, quadrangular and pentagonal anisotropic flow coefficients for $\pi^{\pm}$, $\mathrm{K}^{\pm}$ and p+$\overline{\mathrm{p}}$ in Pb-Pb collisions at $\sqrt{s_\mathrm{{NN}}} = 2.76$ TeV were measured with the ALICE detector at the Large Hadron Collider. The results were obtained with the Scalar Product method, correlating the identified hadrons with reference particles from a different pseudorapidity region. Effects not related to the common event symmetry planes (non-flow) were estimated using correlations in pp collisions and were subtracted from the measurement. The obtained flow coefficients exhibit a clear mass ordering for transverse momentum ($p_{\mathrm{T}}$) values below $\approx$ 3 GeV/$c$. In the intermediate $p_{\mathrm{T}}$ region ($3 < p_{\mathrm{T}} < 6$ GeV/$c$), particles group at an approximate level according to the number of constituent quarks, suggesting that coalescence might be the relevant particle production mechanism in this region. The results for $p_{\mathrm{T}} < 3$ GeV/$c$ are described fairly well by a hydrodynamical model (iEBE-VISHNU) that uses initial conditions generated by A Multi-Phase Transport model (AMPT) and describes the expansion of the fireball using a value of 0.08 for the ratio of shear viscosity to entropy density ($\eta/s$), coupled to a hadronic cascade model (UrQMD). Finally, expectations from AMPT alone fail to quantitatively describe the measurements for all harmonics throughout the measured transverse momentum region. However, the comparison to the AMPT model highlights the importance of the late hadronic rescattering stage to the development of the observed mass ordering at low values of $p_{\mathrm{T}}$ and of coalescence as a particle production mechanism for the particle type grouping at intermediate values of $p_{\mathrm{T}}$ for all harmonics.

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.

In ultrarelativistic heavy-ion collisions, the event-by-event variation of the elliptic flow $v_2$ reflects fluctuations in the shape of the initial state of the system. This allows to select events with the same centrality but different initial geometry. This selection technique, Event Shape Engineering, has been used in the analysis of charge-dependent two- and three-particle correlations in Pb-Pb collisions at $\sqrt{s_{_{\rm NN}}} =2.76$ TeV. The two-particle correlator $\langle \cos(\varphi_\alpha - \varphi_\beta) \rangle$, calculated for different combinations of charges $\alpha$ and $\beta$, is almost independent of $v_2$ (for a given centrality), while the three-particle correlator $\langle \cos(\varphi_\alpha + \varphi_\beta - 2\Psi_2) \rangle$ scales almost linearly both with the event $v_2$ and charged-particle pseudorapidity density. The charge dependence of the three-particle correlator is often interpreted as evidence for the Chiral Magnetic Effect (CME), a parity violating effect of the strong interaction. However, its measured dependence on $v_2$ points to a large non-CME contribution to the correlator. Comparing the results with Monte Carlo calculations including a magnetic field due to the spectators, the upper limit of the CME signal contribution to the three-particle correlator in the 10-50% centrality interval is found to be 26-33% at 95% confidence level.

The measurement of azimuthal correlations of charged particles is presented for Pb-Pb collisions at $\sqrt{s_{\rm NN}}=$ 2.76 TeV and p-Pb collisions at $\sqrt{s_{\rm NN}}=$ 5.02 TeV with the ALICE detector at the CERN Large Hadron Collider. These correlations are measured for the second, third and fourth order flow vector in the pseudorapidity region $|\eta|<0.8$ as a function of centrality and transverse momentum $p_{\rm T}$ using two observables, to search for evidence of $p_{\rm T}$-dependent flow vector fluctuations. For Pb-Pb collisions at 2.76 TeV, the measurements indicate that $p_{\rm T}$-dependent fluctuations are only present for the second order flow vector. Similar results have been found for p-Pb collisions at 5.02 TeV. These measurements are compared to hydrodynamic model calculations with event-by-event geometry fluctuations in the initial state to constrain the initial conditions and transport properties of the matter created in Pb-Pb and p-Pb collisions.