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.

We present a measurement of inclusive J/$\psi$ production in p-Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV as a function of the centrality of the collision, as estimated from the energy deposited in the Zero Degree Calorimeters. The measurement is performed with the ALICE detector down to zero transverse momentum, $p_{\rm T}$, in the backward ($-4.46 < y_{\rm cms} < -2.96$) and forward ($2.03 < y_{\rm cms} < 3.53$) rapidity intervals in the dimuon decay channel and in the mid-rapidity region ($-1.37 < y_{\rm cms} < 0.43$) in the dielectron decay channel. The backward and forward rapidity intervals correspond to the Pb-going and p-going direction, respectively. The $p_{\rm T}$-differential J/$\psi$ production cross section at backward and forward rapidity is measured for several centrality classes, together with the corresponding average $p_{\rm T}$ and $p^2_{\rm T}$ values. The nuclear modification factor, $Q_{\rm pPb}$, is presented as a function of centrality for the three rapidity intervals, and, additionally, at backward and forward rapidity, as a function of $p_{\rm T}$ for several centrality classes. At mid- and forward rapidity, the J/$\psi$ yield is suppressed up to 40% compared to that in pp interactions scaled by the number of binary collisions. The degree of suppression increases towards central p-Pb collisions at forward rapidity, and with decreasing $p_{\rm T}$ of the J/$\psi$. At backward rapidity, the $Q_{\rm pPb}$ is compatible with unity within the total uncertainties, with an increasing trend from peripheral to central p-Pb collisions.

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 elliptic flow, $v_{2}$, of muons from heavy-flavour hadron decays at forward rapidity ($2.5 < y < 4$) is measured in Pb--Pb collisions at $\sqrt{s_{\rm NN}}$~=~2.76 TeV with the ALICE detector at the LHC. The scalar product, two- and four-particle $Q$ cumulants and Lee-Yang zeros methods are used. The dependence of the $v_2$ of muons from heavy-flavour hadron decays on the collision centrality, in the range 0--40\%, and on transverse momentum, $p_{\rm T}$, is studied in the interval $3 < p_{\rm T} < 10$~GeV/$c$. A positive $v_2$ is observed with the scalar product and two-particle $Q$ cumulants in semi-central collisions (10--20\% and 20--40\% centrality classes) for the $p_{\rm T}$ interval from 3 to about 5 GeV/$c$. The $v_2$ magnitude tends to decrease towards more central collisions and with increasing $p_{\rm T}$. It becomes compatible with zero in the interval $6<p_{\rm T}<10~{\rm GeV/}c$. The results are compared to models describing the interaction of heavy quarks and open heavy-flavour hadrons with the high-density medium formed in high-energy heavy-ion collisions. The model calculations describe the measured $v_2$ within uncertainties.

We have studied the transverse-momentum ($p_{\rm T}$) dependence of the inclusive J/$\psi$ production in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV, in three center-of-mass rapidity ($y_{\rm cms}$) regions, down to zero $p_{\rm T}$. Results in the forward and backward rapidity ranges ($2.03 < y_{\rm cms} < 3.53$ and $-4.46 <y_{\rm cms}< -2.96$) are obtained by studying the J/$\psi$ decay to $\mu^+\mu^-$, while the mid-rapidity region ($-1.37 < y_{\rm cms} < 0.43$) is investigated by measuring the ${\rm e}^+{\rm e}^-$ decay channel. The $p_{\rm T}$ dependence of the J/$\psi$ production cross section and nuclear modification factor are presented for each of the rapidity intervals, as well as the J/$\psi$ mean $p_{\rm T}$ values. Forward and mid-rapidity results show a suppression of the J/$\psi$ yield, with respect to pp collisions, which decreases with increasing $p_{\rm T}$. At backward rapidity no significant J/$\psi$ suppression is observed. Theoretical models including a combination of cold nuclear matter effects such as shadowing and partonic energy loss, are in fair agreement with the data, except at forward rapidity and low transverse momentum. The implications of the p-Pb results for the evaluation of cold nuclear matter effects on J/$\psi$ production in Pb-Pb collisions are also discussed.

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.

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

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.

A measurement of continuum dimuon production in proton-copper collisions at 800-GeV incident energy is presented. The dimuons observed in this experiment cover the mass range from 6.5 to 18 GeV near y=0 in the proton-nucleon center-of-momentum frame. Scaling forms of the cross section for the continuum are compared with the results of other experiments in the context of the parton model and quantum chromodynamics. The present limitations of such scaling comparisons are discussed.