Measurements of two- and multi-particle angular correlations in pp collisions at sqrt(s) = 5, 7, and 13 TeV are presented as a function of charged-particle multiplicity. The data, corresponding to integrated luminosities of 1.0 inverse picobarn (5 TeV), 6.2 inverse picobarns (7 TeV), and 0.7 inverse picobarns (13 TeV), were collected using the CMS detector at the LHC. The second-order (v[2]) and third-order (v[3]) azimuthal anisotropy harmonics of unidentified charged particles, as well as v[2] of K0 short and Lambda/anti-Lambda particles, are extracted from long-range two-particle correlations as functions of particle multiplicity and transverse momentum. For high-multiplicity pp events, a mass ordering is observed for the v[2] values of charged hadrons (mostly pions), K0 short, and Lambda/anti-Lambda, with lighter particle species exhibiting a stronger azimuthal anisotropy signal below pt of about 2 GeV/c. For 13 TeV data, the v[2] signals are also extracted from four- and six-particle correlations for the first time in pp collisions, with comparable magnitude to those from two-particle correlations. These observations are similar to those seen in pPb and PbPb collisions, and support the interpretation of a collective origin for the observed long-range correlations in high-multiplicity pp collisions.

The second-order azimuthal anisotropy Fourier harmonics, v2, are obtained in pPb and PbPb collisions over a wide pseudorapidity (eta) range based on correlations among six or more charged particles. The pPb data, corresponding to an integrated luminosity of 35 inverse nanobarns, were collected during the 2013 LHC pPb run at a nucleon-nucleon center-of-mass energy of 5.02 TeV by the CMS experiment. A sample of semi-peripheral PbPb collision data at sqrt(s[NN])= 2.76 TeV, corresponding to an integrated luminosity of 2.5 inverse microbarns and covering a similar range of particle multiplicities as the pPb data, is also analyzed for comparison. The six- and eight-particle cumulant and the Lee-Yang zeros methods are used to extract the v2 coefficients, extending previous studies of two- and four-particle correlations. For both the pPb and PbPb systems, the v2 values obtained with correlations among more than four particles are consistent with previously published four-particle results. These data support the interpretation of a collective origin for the previously observed long-range (large Delta[eta]) correlations in both systems. The ratios of v2 values corresponding to correlations including different numbers of particles are compared to theoretical predictions that assume a hydrodynamic behavior of a pPb system dominated by fluctuations in the positions of participant nucleons. These results provide new insights into the multi-particle dynamics of collision systems with a very small overlapping region.

The nuclear modification factors of J/$\psi$ and $\psi$(2S) mesons are measured in PbPb collisions at a centre-of-mass energy per nucleon pair of $\sqrt{s_{\mathrm{NN}}} =$ 5.02 TeV. The analysis is based on PbPb and pp data samples collected by CMS at the LHC in 2015, corresponding to integrated luminosities of 464 $\mu$b$^{-1}$ and 28 pb$^{-1}$, respectively. The measurements are performed in the dimuon rapidity range of $|y| <$ 2.4 as a function of centrality, rapidity, and transverse momentum (p$_\mathrm{T}$) from p$_\mathrm{T}=$ 3 GeV/$c$ in the most forward region and up to 50 GeV/$c$. Both prompt and nonprompt (coming from b hadron decays) mesons are observed to be increasingly suppressed with centrality, with a magnitude similar to the one observed at $\sqrt{s_{\mathrm{NN}}}=$ 2.76 TeV for the two J/$\psi$ meson components. No dependence on rapidity is observed for either prompt or nonprompt J/$\psi$ mesons. An indication of a lower prompt J/$\psi$ meson suppression at p$_\mathrm{T} >$ 25 GeV/$c$ is seen with respect to that observed at intermediate p$_\mathrm{T}$. The prompt $\psi$(2S) meson yield is found to be more suppressed than that of the prompt J/$\psi$ mesons in the entire p$_\mathrm{T}$ range.

The production cross sections of the B+, B0, and B0s mesons, and of their charge conjugates, are measured via exclusive hadronic decays in pPb collisions at the center-of-mass energy sqrt(s_NN) = 5.02 TeV with the CMS detector at the CERN LHC. The data set used for this analysis corresponds to an integrated luminosity of 34.6 inverse nanobarns. The production cross sections are measured in the transverse momentum range between 10 and 60 GeV/c. No significant modification is observed compared to proton-proton perturbative QCD calculations scaled by the number of incoherent nucleon-nucleon collisions. These results provide a baseline for the study of in-medium b quark energy loss in PbPb collisions.

The first study of W boson production in pPb collisions is presented, for bosons decaying to a muon or electron, and a neutrino. The measurements are based on a data sample corresponding to an integrated luminosity of 34.6 inverse nanobarns at a nucleon-nucleon centre-of-mass energy of sqrt(s[NN]) = 5.02 TeV, collected by the CMS experiment. The W boson differential cross sections, lepton charge asymmetry, and forward-backward asymmetries are measured for leptons of transverse momentum exceeding 25 GeV, and as a function of the lepton pseudorapidity in the abs(eta[lab]) < 2.4 range. Deviations from the expectations based on currently available parton distribution functions are observed, showing the need for including W boson data in nuclear parton distribution global fits.

A systematic study of the factorization of long-range azimuthal two-particle correlations into a product of single-particle anisotropies is presented as a function of pt and eta of both particles, and as a function of the particle multiplicity in PbPb and pPb collisions. The data were taken with the CMS detector for PbPb collisions at sqrt(s[NN]) = 2.76 TeV and pPb collisions at sqrt(s[NN]) = 5.02 TeV, covering a very wide range of multiplicity. Factorization is observed to be broken as a function of both particle pt and eta. When measured with particles of different pt, the magnitude of the factorization breakdown for the second Fourier harmonic reaches 20% for very central PbPb collisions but decreases rapidly as the multiplicity decreases. The data are consistent with viscous hydrodynamic predictions, which suggest that the effect of factorization breaking is mainly sensitive to the initial-state conditions rather than to the transport properties (e.g., shear viscosity) of the medium. The factorization breakdown is also computed with particles of different eta. The effect is found to be weakest for mid-central PbPb events but becomes larger for more central or peripheral PbPb collisions, and also for very high-multiplicity pPb collisions. The eta-dependent factorization data provide new insights to the longitudinal evolution of the medium formed in heavy ion collisions.

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.

Results on two-particle angular correlations for charged particles produced in pp collisions at a center-of-mass energy of 13 TeV are presented. The data were taken with the CMS detector at the LHC and correspond to an integrated luminosity of about 270 inverse nanobarns. The correlations are studied over a broad range of pseudorapidity (abs(eta) < 2.4) and over the full azimuth (phi) as a function of charged particle multiplicity and transverse momentum (pt). In high-multiplicity events, a long-range (abs(Delta eta) > 2.0), near-side (Delta phi approximately 0) structure emerges in the two-particle Delta eta-Delta phi correlation functions. The magnitude of the correlation exhibits a pronounced maximum in the range 1.0 < pt < 2.0 GeV/c and an approximately linear increase with the charged particle multiplicity, with an overall correlation strength similar to that found in earlier pp data at sqrt(s) = 7 TeV. The present measurement extends the study of near-side long-range correlations up to charged particle multiplicities of N[ch] approximately 180, a region so far unexplored in pp collisions. The observed long-range correlations are compared to those seen in pp, pPb, and PbPb collisions at lower collision energies.

Measurements of two-particle angular correlations between an identified strange hadron (K0S or Lambda/anti-Lambda) and a charged particle, emitted in pPb collisions, are presented over a wide range in pseudorapidity and full azimuth. The data, corresponding to an integrated luminosity of approximately 35 inverse nanobarns, were collected at a nucleon-nucleon center-of-mass energy (sqrt(s[NN])) of 5.02 TeV with the CMS detector at the LHC. The results are compared to semi-peripheral PbPb collision data at sqrt(s[NN]) = 2.76 TeV, covering similar charged-particle multiplicities in the events. The observed azimuthal correlations at large relative pseudorapidity are used to extract the second-order (v[2]) and third-order (v[3]) anisotropy harmonics of K0S and Lambda/anti-Lambda particles. These quantities are studied as a function of the charged-particle multiplicity in the event and the transverse momentum of the particles. For high-multiplicity pPb events, a clear particle species dependence of v[2] and v[3] is observed. For pt < 2 GeV, the v[2] and v[3] values of K0S particles are larger than those of Lambda/anti-Lambda particles at the same pt. This splitting effect between two particle species is found to be stronger in pPb than in PbPb collisions in the same multiplicity range. When divided by the number of constituent quarks and compared at the same transverse kinetic energy per quark, both v[2] and v[3] for K0S particles are observed to be consistent with those for Lambda/anti-Lambda particles at the 10% level in pPb collisions. This consistency extends over a wide range of particle transverse kinetic energy and event multiplicities.

The production yields of Upsilon(1S), Upsilon(2S), and Upsilon(3S) quarkonium states are measured through their decays into muon pairs in the CMS detector, in PbPb and pp collisions at the centre-of-mass energy per nucleon pair of 2.76 TeV. The data correspond to integrated luminosities of 166 inverse microbarns and 5.4 inverse picobarns for PbPb and pp collisions, respectively. Differential production cross sections are reported as functions of Upsilon rapidity y up to 2.4, and transverse momentum pT up to 20 GeV/c. A strong centrality-dependent suppression is observed in PbPb relative to pp collisions, by factors of up to approximately 2 and 8, for the Upsilon(1S) and Upsilon(2S) states, respectively. No significant dependence of this suppression is observed as a function of y or pT. The Upsilon(3S) state is not observed in PbPb collisions, which corresponds to a suppression for the centrality-integrated data by at least a factor of approximately 7 at a 95% confidence level. The observed suppression is in agreement with theoretical scenarios modeling the sequential melting of quarkonium states in a quark gluon plasma.