Production cross sections of $\Upsilon$(1S), $\Upsilon$(2S), and $\Upsilon$(3S) states decaying into $\mu^+\mu^-$ in proton-lead (pPb) collisions are reported using data collected by the CMS experiment at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV. A comparison is made with corresponding cross sections obtained with pp data measured at the same collision energy and scaled by the Pb nucleus mass number. The nuclear modification factor for $\Upsilon$(1S) is found to be $R_\mathrm{pPb}(\Upsilon(1S))$ = 0.806 $\pm$ 0.024 (stat) $\pm$ 0.059 (syst). Similar results for the excited states indicate a sequential suppression pattern, such that $R_\mathrm{pPb}(\Upsilon(1S))$$\gt$$R_\mathrm{pPb}(\Upsilon(2S))$$\gt$$R_\mathrm{pPb}(\Upsilon(3S))$. The suppression is much less pronounced in pPb than in PbPb collisions, and independent of transverse momentum $p_\mathrm{T}^\Upsilon$ and center-of-mass rapidity $y_\mathrm{CM}^\Upsilon$ of the individual $\Upsilon$ state in the studied range $p_\mathrm{T}^\Upsilon$$\lt$ 30 GeV$/c$ and $\vert y_\mathrm{CM}^\Upsilon\vert$$\lt$ 1.93. Models that incorporate sequential suppression of bottomonia in pPb collisions are in better agreement with the data than those which only assume initial-state modifications.

The relative modification of the prompt psi(2S) and J/psi yields from pp to PbPb collisions, at the center of mass energy of 5.02 TeV per nucleon pair, is presented. The analysis is based on pp and PbPb data samples collected by the CMS experiment at the LHC in 2015, corresponding to integrated luminosities of 28.0 inverse picobarns and 464 inverse microbarns, respectively. The double ratio of measured yields of prompt charmonia reconstructed through their decays into muon pairs, (N[psi(2S)]/N[J/psi])[PbPb] / (N[psi(2S)]/N[J/psi])[pp], is determined as a function of PbPb collision centrality and charmonium transverse momentum pt, in two kinematic intervals: abs(y) < 1.6 covering 6.5 < pt < 30 GeV/c and 1.6 < abs(y) < 2.4 covering 3 < pt < 30 GeV/c. The centrality-integrated double ratios are 0.36 +/- 0.08 (stat) +/-0.05 (syst) in the first interval and 0.24 +/- 0.22 (stat) +/- 0.09 (syst) in the second. The double ratio is lower than unity in all the measured bins, suggesting that the psi(2S) yield is more suppressed than the J/psi yield in the explored phase space.

A search for Z bosons in the mu^+mu^- decay channel has been performed in PbPb collisions at a nucleon-nucleon centre of mass energy = 2.76 TeV with the CMS detector at the LHC, in a 7.2 inverse microbarn data sample. The number of opposite-sign muon pairs observed in the 60--120 GeV/c2 invariant mass range is 39, corresponding to a yield per unit of rapidity (y) and per minimum bias event of (33.8 ± 5.5 (stat) ± 4.4 (syst)) 10^{-8}, in the |y|<2.0 range. Rapidity, transverse momentum, and centrality dependencies are also measured. The results agree with next-to-leading order QCD calculations, scaled by the number of incoherent nucleon-nucleon collisions.

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.

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.

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 cross sections for $\Upsilon$(1S), $\Upsilon$(2S), and $\Upsilon$(3S) production in lead-lead (PbPb) and proton-proton (pp) collisions at $\sqrt{s_{_\mathrm{NN}}} =$ 5.02 TeV have been measured using the CMS detector at the LHC. The nuclear modification factors, R$_\mathrm{AA}$, derived from the PbPb-to-pp ratio of yields for each state, are studied as functions of meson rapidity and transverse momentum, as well as PbPb collision centrality. The yields of all three states are found to be significantly suppressed, and compatible with a sequential ordering of the suppression, R$_\mathrm{AA}$($\Upsilon$(1S)) $>$ R$_\mathrm{AA}$($\Upsilon$(2S)) $>$ R$_\mathrm{AA}$($\Upsilon$(3S)) . The suppression of $\Upsilon$(1S) is larger than that seen at $\sqrt{s_{_\mathrm{NN}}} =$ 2.76 TeV, although the two are compatible within uncertainties. The upper limit on the R$_\mathrm{AA}$ of $\Upsilon$(3S) integrated over $p_\mathrm{T}$ and rapidity is 0.094 at 95% confidence level, which is the strongest suppression observed for any hadron species in heavy ion collisions to date.

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 second-order Fourier coefficients ($v_2$) characterizing the azimuthal distribution of $\Upsilon$(1S) and $\Upsilon$(2S) mesons arising from PbPb collisions at $\sqrt{s_\mathrm{NN}} =$ 5.02 TeV are studied. The $\Upsilon$ mesons are reconstructed in their dimuon decay channel, as measured by the CMS detector. The data set corresponds to an integrated luminosity of 1.7 nb$^{-1}$. The scalar product method is used to extract the $v_2$ coefficients of the azimuthal distribution. Results are reported for the rapidity range $|y|$$\lt$ 2.4, with the transverse momentum 0 $\lt$$p_\mathrm{T}$$\lt$ 50 GeV/$c$, and in three centrality ranges of 10-30%, 30-50% and 50-90%. In contrast to the J/$\psi$ mesons, the measured $v_2$ values for the $\Upsilon$ mesons are found to be consistent with zero.

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.