The production cross sections of the prompt charmed mesons D$^0$, D$^+$, D$^{*+}$ and D$_s$ were measured at mid-rapidity in p-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}}=5.02$ TeV with the ALICE detector at the LHC. D mesons were reconstructed from their decays D$^0\rightarrow{\rm K}^-\pi^+$, D$^+\rightarrow{\rm K}^-\pi^+\pi^+$, D$^{*+}\rightarrow D^0\pi^+$, D$_s^+\rightarrow\phi\pi^+\rightarrow{\rm K}^-{\rm K}^+\pi^+$, and their charge conjugates. The $p_{\rm T}$-differential production cross sections were measured at mid-rapidity in the interval $1<p_{\rm T}<24$ GeV/$c$ for D$^0$, D$^+$ and D$^{*+}$ mesons and in $2<p_{\rm T}<12$ GeV/$c$ for D$_s$ mesons, using an analysis method based on the selection of decay topologies displaced from the interaction vertex. The production cross sections of the D$^0$, D$^+$ and D$^{*+}$ mesons were also measured in three $p_{\rm T}$ intervals as a function of the rapidity $y_{\rm cms}$ in the centre-of-mass system in $-1.26<y_{\rm cms}<0.34$. In addition, the prompt D$^0$ cross section was measured in pp collisions at $\sqrt{s}=7$ TeV and p-Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV down to $p_{\rm T}=0$ using an analysis technique that is based on the estimation and subtraction of the combinatorial background, without reconstruction of the D$^0$ decay vertex. The nuclear modification factor $R_{\rm pPb}(p_{\rm T})$, defined as the ratio of the $p_{\rm T}$-differential D-meson cross section in p-Pb collisions and that in pp collisions scaled by the mass number of the Pb nucleus, was calculated for the four D-meson species and found to be compatible with unity within experimental uncertainties. The results are compared to theoretical calculations that include cold-nuclear-matter effects and to transport model calculations incorporating the interactions of charm quarks with an expanding deconfined medium.

The first study of $\phi$-meson production in p-Pb collisions at forward and backward rapidity, at a nucleon-nucleon centre-of-mass energy $\sqrt{s_{\rm NN}} = 5.02$~TeV, has been performed with the ALICE apparatus at the LHC. The $\phi$-mesons have been identified in the dimuon decay channel in the transverse momentum ($p_{\rm T}$) range $1 < p_{\rm T} < 7$ GeV/$c$, both in the p-going ($2.03 < y < 3.53$) and the Pb-going ($-4.46 < y < -2.96$) directions, where $y$ stands for the rapidity in the nucleon-nucleon centre-of-mass, the integrated luminosity amounting to $5.01 \pm 0.19$~nb$^{-1}$ and $5.81 \pm 0.20$~nb$^{-1}$, respectively, for the two data samples. Differential cross sections as a function of transverse momentum and rapidity are presented. The forward-backward ratio for $\phi$-meson production is measured for $2.96<|y|<3.53$, resulting in a ratio $\sim 0.5$ with no significant $p_{\rm T}$ dependence within the uncertainties. The $p_{\rm T}$ dependence of the $\phi$ nuclear modification factor $R_{\rm pPb}$ exhibits an enhancement up to a factor 1.6 at $p_{\rm T}$ = 3-4 GeV/$c$ in the Pb-going direction. The $p_{\rm T}$ dependence of the $\phi$-meson cross section in pp collisions at $\sqrt{s}$ = 2.76 TeV, which is used to determine a reference for the p-Pb results, is also presented here for $1 < p_{\rm T} < 5$ GeV/$c$ and $2.5 <y < 4$ for a $78 \pm 3$~nb$^{-1}$ integrated luminosity sample.

The charged-particle pseudorapidity density measured over 4 units of pseudorapidity in non-single-diffractive (NSD) p-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}} = 5.02$ TeV is presented. The average value at midrapidity is measured to be $16.81 \pm 0.71$ (syst.), which corresponds to $2.14 \pm 0.17$ (syst.) per participating nucleon. This is 16% lower than in NSD pp collisions interpolated to the same collision energy, and 84% higher than in d-Au collisions at $\sqrt{s_{\rm NN}} = 0.2$ TeV. The measured pseudorapidity density in p-Pb collisions is compared to model predictions, and provides new constraints on the description of particle production in high-energy nuclear collisions.

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.

Measurements of multi-particle azimuthal correlations (cumulants) for charged particles in p-Pb and Pb-Pb collisions are presented. They help address the question of whether there is evidence for global, flow-like, azimuthal correlations in the p-Pb system. Comparisons are made to measurements from the larger Pb-Pb system, where such evidence is established. In particular, the second harmonic two-particle cumulants are found to decrease with multiplicity, characteristic of a dominance of few-particle correlations in p-Pb collisions. However, when a $|\Delta \eta|$ gap is placed to suppress such correlations, the two-particle cumulants begin to rise at high-multiplicity, indicating the presence of global azimuthal correlations. The Pb-Pb values are higher than the p-Pb values at similar multiplicities. In both systems, the second harmonic four-particle cumulants exhibit a transition from positive to negative values when the multiplicity increases. The negative values allow for a measurement of $v_{2}\{4\}$ to be made, which is found to be higher in Pb-Pb collisions at similar multiplicities. The second harmonic six-particle cumulants are also found to be higher in Pb-Pb collisions. In Pb-Pb collisions, we generally find $v_{2}\{4\} \simeq v_{2}\{6\}\neq 0$ which is indicative of a Bessel-Gaussian function for the $v_{2}$ distribution. For very high-multiplicity Pb-Pb collisions, we observe that the four- and six-particle cumulants become consistent with 0. Finally, third harmonic two-particle cumulants in p-Pb and Pb-Pb are measured. These are found to be similar for overlapping multiplicities, when a $|\Delta\eta| > 1.4$ gap is placed.

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.

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.

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.

The ALICE Collaboration has studied the inclusive production of the charmonium state $\psi(2S)$ in proton-lead (p-Pb) collisions at the nucleon-nucleon centre of mass energy $\sqrt{s_{NN}}$ = 5.02 TeV at the CERN LHC. The measurement was performed at forward ($2.03<y_{cms}<3.53$) and backward ($-4.46<y_{cms}<-2.96$) centre of mass rapidities, studying the decays into muon pairs. In this paper, we present the inclusive production cross sections $\sigma_{\psi(2S)}$, both integrated and as a function of the transverse momentum $p_{T}$, for the two $y_{cms}$ domains. The results are compared to those obtained for the 1S vector state (J/$\psi$), by showing the ratios between the production cross sections, as well as the double ratios $[\sigma_{\psi(2S)}/\sigma_{J/\psi}]_{pPb}/[\sigma_{\psi(2S)}/\sigma_{J/\psi}]_{pp}$ between p-Pb and proton-proton collisions. Finally, the nuclear modification factor for inclusive $\psi(2S)$ is evaluated and compared to the measurement of the same quantity for J/$\psi$ and to theoretical models including parton shadowing and coherent energy loss mechanisms. The results show a significantly larger suppression of the $\psi(2S)$ compared to that measured for J/$\psi$ and to models. These observations represent a clear indication for sizeable final state effects on $\psi(2S)$ production.