The production of Z bosons associated with jets is measured in pp collisions at $\sqrt{s}$ = 13 TeV with data recorded with the CMS experiment at the LHC corresponding to an integrated luminosity of 36.3 fb$^{-1}$. The multiplicity of jets with transverse momentum $p_\mathrm{T}$$\gt$ 30 GeV is measured for different regions of the Z boson's $p_\mathrm{T}$(Z), from lower than 10 GeV to higher than 100 GeV. The azimuthal correlation $\Delta \phi$ between the Z boson and the leading jet, as well as the correlations between the two leading jets are measured in three regions of $p_\mathrm{T}$(Z). The measurements are compared with several predictions at leading and next-to-leading orders, interfaced with parton showers. Predictions based on transverse-momentum dependent parton distributions and corresponding parton showers give a good description of the measurement in the regions where multiple parton interactions and higher jet multiplicities are not important. The effects of multiple parton interactions are shown to be important to correctly describe the measured spectra in the low $p_\mathrm{T}$(Z) regions.

The double differential cross sections of the Drell-Yan lepton pair ($\ell^+\ell^-$, dielectron or dimuon) production are measured as functions of the invariant mass $m_{\ell\ell}$, transverse momentum $p_\mathrm{T}(\ell\ell)$, and $\phi^*_\eta$. The $\phi^*_\eta$ observable, derived from angular measurements of the leptons and highly correlated with $p_\mathrm{T}(\ell\ell)$, is used to probe the low-$p_\mathrm{T}(\ell\ell)$ region in a complementary way. Dilepton masses up to 1 TeV are investigated. Additionally, a measurement is performed requiring at least one jet in the final state. To benefit from partial cancellation of the systematic uncertainty, the ratios of the differential cross sections for various $m_{\ell\ell}$ ranges to those in the Z mass peak interval are presented. The collected data correspond to an integrated luminosity of 36.3 fb$^{-1}$ of proton-proton collisions recorded with the CMS detector at the LHC at a centre-of-mass energy of 13 TeV. Measurements are compared with predictions based on perturbative quantum chromodynamics, including soft-gluon resummation.

The dependence of the ratio between the B$_\mathrm{s}^0$ and B$^+$ hadron production fractions, $f_\mathrm{s} / f_\mathrm{u}$, on the transverse momentum ($p_\mathrm{T}$) and rapidity of the B mesons is studied using the decay channels B$_\mathrm{s}^0$$\to$ J$/\psi\,\phi$ and B$^+$$\to$ J$/\psi$ K$^+$. The analysis uses a data sample of proton-proton collisions at a center-of-mass energy of 13 TeV, collected by the CMS experiment in 2018 and corresponding to an integrated luminosity of 61.6 fb$^{-1}$. The $f_\mathrm{s} / f_\mathrm{u}$ ratio is observed to depend on the B $p_\mathrm{T}$ and to be consistent with becoming asymptotically constant at large $p_\mathrm{T}$. No rapidity dependence is observed. The ratio of the B$^0$ to B$^+$ hadron production fractions, $f_\mathrm{d} / f_\mathrm{u}$, measured using the B$^0$$\to$ J$/\psi$ K$^{*0}$ decay channel, is found to be consistent with unity and independent of $p_\mathrm{T}$ and rapidity, as expected from isospin invariance.

The $B_\mathrm{c}^+$ meson is observed for the first time in heavy ion collisions. Data from the CMS detector are used to study the production of the $B_\mathrm{c}^+$ meson in lead-lead (PbPb) and proton-proton (pp) collisions at a center-of-mass energy per nucleon pair of $\sqrt{s_{\mathrm{NN}}} =$ 5.02 TeV, via the $B_\mathrm{c}^+ \to (J/\psi\to\mu^+\mu^-)\mu^+\nu_\mu$ decay. The $B_\mathrm{c}^+$ nuclear modification factor, derived from the PbPb-to-pp ratio of production cross sections, is measured in two bins of the trimuon transverse momentum and of the PbPb collision centrality. The B$_\mathrm{c}^+$ meson is shown to be less suppressed than quarkonia and most of the open heavy-flavor mesons, suggesting that effects of the hot and dense nuclear matter created in heavy ion collisions contribute to its production. This measurement sets forth a promising new probe of the interplay of suppression and enhancement mechanisms in the production of heavy-flavor mesons in the quark-gluon plasma.

Double-parton scattering is investigated using events with a Z boson and jets. The Z boson is reconstructed using only the dimuon channel. The measurements are performed with proton-proton collision data recorded by the CMS experiment at the LHC at $\sqrt{s} =$ 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$ collected in the year 2016. Differential cross sections of Z + $\geq$ 1 jet and Z + $\geq$ 2 jets are measured with transverse momentum of the jets above 20 GeV and pseudorapidity $|\eta|$$\lt$ 2.4. Several distributions with sensitivity to double-parton scattering effects are measured as functions of the angle and the transverse momentum imbalance between the Z boson and the jets. The measured distributions are compared with predictions from several event generators with different hadronization models and different parameter settings for multiparton interactions. The measured distributions show a dependence on the hadronization and multiparton interaction simulation parameters, and are important input for future improvements of the simulations.

The production of $\Upsilon (nS)$ mesons ($n=1,2,3$) in $p$Pb and Pb$p$ collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{NN}}=8.16$ TeV is measured by the LHCb experiment, using a data sample corresponding to an integrated luminosity of 31.8 nb$^{-1}$. The $\Upsilon (nS)$ mesons are reconstructed through their decays into two opposite-sign muons. The measurements comprise the differential production cross-sections of the $\Upsilon (1S)$ and $\Upsilon (2S)$ states, their forward-to-backward ratios and nuclear modification factors, performed as a function of the transverse momentum $p_{\mathrm{T}}$ and rapidity in the nucleon-nucleon centre-of-mass frame $y^*$ of the $\Upsilon (nS)$ states, in the kinematic range $p_{\rm{T}}<25$ GeV/$c$ and $1.5<y^*<4.0$ ($-5.0<y^*<-2.5$) for $p$Pb (Pb$p$) collisions. In addition, production cross-sections for $\Upsilon (3S)$ are measured integrated over phase space and the production ratios between all three $\Upsilon (nS)$ states are determined. The measurements are compared to theoretical predictions and suppressions for quarkonium in $p$Pb collisions are observed.

The ratios of the production cross sections between the excited $\Upsilon$(2S) and $\Upsilon$(3S) mesons and the $\Upsilon$(1S) ground state, detected via their decay into two muons, are studied as a function of the number of charged particles in the event. The data are from proton-proton collisions at $\sqrt{s} =$ 7 TeV, corresponding to an integrated luminosity of 4.8 fb$^{-1}$, collected with the CMS detector at the LHC. Evidence of a decrease in these ratios as a function of the particle multiplicity is observed, more pronounced at low transverse momentum $p_\mathrm{T}^{\mu\mu}$. For $\Upsilon$(nS) mesons with $p_\mathrm{T}^{\mu\mu}$ $\gt$ 7 GeV, where most of the data were collected, the correlation with multiplicity is studied as a function of the underlying event transverse sphericity and the number of particles in a cone around the $\Upsilon$(nS) direction. The ratios are found to be multiplicity independent for jet-like events. The mean $p_\mathrm{T}^{\mu\mu}$ values for the $\Upsilon$(nS) states as a function of particle multiplicity are also measured and found to grow more steeply as their mass increases.

The dependence of inclusive jet production in proton-proton collisions with a center-of-mass energy of 13 TeV on the distance parameter $R$ of the anti-$k_\mathrm{T}$ algorithm is studied using data corresponding to integrated luminosities up to 35.9 fb$^{-1}$ collected by the CMS experiment in 2016. The ratios of the inclusive cross sections as functions of transverse momentum $p_\mathrm{T}$ and rapidity $y$, for $R$ in the range 0.1 to 1.2 to those using $R = $ 0.4 are presented in the region 84 $\lt$ $p_\mathrm{T}$ $\lt$ 1588 GeV and $|y|$ $\lt$ 2.0. The results are compared to calculations at leading and next-to-leading order in the strong coupling constant using different parton shower models. The variation of the ratio of cross sections with $R$ is well described by calculations including a parton shower model, but not by a leading-order quantum chromodynamics calculation including nonperturbative effects. The agreement between the data and the theoretical predictions for the ratios of cross sections is significantly improved when next-to-leading order calculations with nonperturbative effects are used.

The production of K_short mesons in pp collisions at a centre-of-mass energy of 0.9 TeV is studied with the LHCb detector at the Large Hadron Collider. The luminosity of the analysed sample is determined using a novel technique, involving measurements of the beam currents, sizes and positions, and is found to be 6.8 +/- 1.0 microbarn^-1. The differential prompt K_short production cross-section is measured as a function of the K_short transverse momentum and rapidity in the region 0 &lt; pT &lt; 1.6 GeV/c and 2.5 &lt; y &lt; 4.0. The data are found to be in reasonable agreement with previous measurements and generator expectations.

Charm production at the LHC in pp collisions at sqrt(s)=7 TeV is studied with the LHCb detector. The decays D0 -> K- pi+, D+ -> K- pi+ pi+, D*+ -> D0(K- pi+) pi+, D_s+ -> phi(K- K+) pi+, Lambda_c+ -> p K- pi+, and their charge conjugates are analysed in a data set corresponding to an integrated luminosity of 15 nb^{-1}. Differential cross-sections dsigma/dp_T are measured for prompt production of the five charmed hadron species in bins of transverse momentum and rapidity in the region 0 < p_T < 8 GeV/c and 2.0 < y < 4.5. Theoretical predictions are compared to the measured differential cross-sections. The integrated cross-sections of the charm hadrons are computed in the above p_T-y range, and their ratios are reported. A combination of the five integrated cross-section measurements gives sigma(c\bar{c})_{p_T < 8 GeV/c, 2.0 < y < 4.5} = 1419 +/- 12 (stat) +/- 116 (syst) +/- 65 (frag) microbarn, where the uncertainties are statistical, systematic, and due to the fragmentation functions.

The production of the $\Lambda$(1520) baryonic resonance has been measured at midrapidity in inelastic pp collisions at $\sqrt{s}$ = 7 TeV and in p-Pb collisions at $\sqrt{s_{\rm{NN}}}$ = 5.02 TeV for non-single diffractive events and in multiplicity classes. The resonance is reconstructed through its hadronic decay channel $\Lambda$(1520) $\rightarrow$ pK$^{-}$ and the charge conjugate with the ALICE detector. The integrated yields and mean transverse momenta are calculated from the measured transverse momentum distributions in pp and p-Pb collisions. The mean transverse momenta follow mass ordering as previously observed for other hyperons in the same collision systems. A Blast-Wave function constrained by other light hadrons ($\pi$, K, K$_{\rm{S}}^0$, p, $\Lambda$) describes the shape of the $\Lambda$(1520) transverse momentum distribution up to 3.5 GeV/$c$ in p-Pb collisions. In the framework of this model, this observation suggests that the $\Lambda(1520)$ resonance participates in the same collective radial flow as other light hadrons. The ratio of the yield of $\Lambda(1520)$ to the yield of the ground state particle $\Lambda$ remains constant as a function of charged-particle multiplicity, suggesting that there is no net effect of the hadronic phase in p-Pb collisions on the $\Lambda$(1520) yield.

The production of Upsilon(1S), Upsilon(2S) and Upsilon(3S) mesons in proton-proton collisions at the centre-of-mass energy of sqrt(s)=7 TeV is studied with the LHCb detector. The analysis is based on a data sample of 25 pb-1 collected at the Large Hadron Collider. The Upsilon mesons are reconstructed in the decay mode Upsilon -&gt; mu+ mu- and the signal yields are extracted from a fit to the mu+ mu- invariant mass distributions. The differential production cross-sections times dimuon branching fractions are measured as a function of the Upsilon transverse momentum pT and rapidity y, over the range pT &lt; 15 GeV/c and 2.0 &lt; y &lt; 4.5. The cross-sections times branching fractions, integrated over these kinematic ranges, are measured to be sigma(pp -&gt; Upsilon(1S) X) x B(Upsilon(1S)-&gt;mu+ mu-) = 2.29 {\pm} 0.01 {\pm} 0.10 -0.37 +0.19 nb, sigma(pp -&gt; Upsilon(2S) X) x B(Upsilon(2S)-&gt;mu+ mu-) = 0.562 {\pm} 0.007 {\pm} 0.023 -0.092 +0.048 nb, sigma(pp -&gt; Upsilon(3S) X) x B(Upsilon(3S)-&gt;mu+ mu-) = 0.283 {\pm} 0.005 {\pm} 0.012 -0.048 +0.025 nb, where the first uncertainty is statistical, the second systematic and the third is due to the unknown polarisation of the three Upsilon states.

Differential cross sections for the Drell-Yan process, including Z boson production, using the dimuon decay channel are measured in proton-lead (pPb) collisions at a nucleon-nucleon centre-of-mass energy of 8.16 TeV. A data sample recorded with the CMS detector at the LHC is used, corresponding to an integrated luminosity of 173 nb$^{-1}$. The differential cross section as a function of the dimuon mass is measured in the range 15-600 GeV, for the first time in proton-nucleus collisions. It is also reported as a function of dimuon rapidity over the mass ranges 15-60 GeV and 60-120 GeV, and ratios for the p-going over the Pb-going beam directions are built. In both mass ranges, the differential cross sections as functions of the dimuon transverse momentum $p_\mathrm{T}$ and of a geometric variable $\phi^*$ are measured, where $\phi^*$ highly correlates with $p_\mathrm{T}$ but is determined with higher precision. In the Z mass region, the rapidity dependence of the data indicate a modification of the distribution of partons within a lead nucleus as compared to the proton case. The data are more precise than predictions based upon current models of parton distributions.

The inclusive $D_s^{\pm}$ production asymmetry is measured in $pp$ collisions collected by the LHCb experiment at centre-of-mass energies of $\sqrt{s} =7$ and 8 TeV. Promptly produced $D_s^{\pm}$ mesons are used, which decay as $D_s^{\pm}\to\phi\pi^{\pm}$, with $\phi\to K^+K^-$. The measurement is performed in bins of transverse momentum, $p_{\rm T}$, and rapidity, $y$, covering the range $2.5

Using a data sample corresponding to an integrated luminosity of $2.0\,fb^{-1}$, collected by the LHCb experiment, the production of the $\eta_c(1S)$ state in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=13 \text{ TeV}$ is studied in the rapidity range ${2.0 < y < 4.5}$ and in the transverse momentum range ${6.5 < p_{T} < 14.0\text{ GeV}}$. The cross-section for prompt production of $\eta_c(1S)$ mesons relative to that of the $J/\psi$ meson is measured using the ${p\bar{p}}$ decay mode and is found to be ${\sigma_{\eta_c(1S)}/\sigma_{J/\psi} = 1.69 \pm 0.15 \pm 0.10 \pm 0.18}$. The quoted uncertainties are, in order, statistical, systematic and due to uncertainties on the branching fractions of the ${J/\psi\to p \bar{p}}$ and ${\eta_c\to p \bar{p}}$ decays. The prompt $\eta_c(1S)$ production cross-section is determined to be ${\sigma_{\eta_c(1S)} = 1.26 \pm 0.11\pm 0.08 \pm 0.14 \,\mu b}$, where the last uncertainty includes that on the ${J/\psi}$ meson cross-section. The ratio of the branching fractions of $b$-hadron decays to the $\eta_c(1S)$ and ${J/\psi}$ states is measured to be ${\mathcal{B}_{b\to\eta_c X}/\mathcal{B}_{b\to J/\psi X} = 0.48 \pm 0.03 \pm 0.03 \pm 0.05}$, where the last uncertainty is due to those on the branching fractions of the ${J/\psi \to p \bar{p}}$ and ${\eta_c\to p \bar{p}}$ decays. The difference between the ${J/\psi}$ and $\eta_c(1S)$ masses is also determined to be ${113.0 \pm 0.7 \pm 0.1\text{ MeV}}$, which is the most precise single measurement of this quantity to date.

The production of Upsilon(1S), Upsilon(2S) and Upsilon(3S) mesons decaying into the dimuon final state is studied with the LHCb detector using a data sample corresponding to an integrated luminosity of 3.3 pb^{-1} collected in proton-proton collisions at a centre-of-mass energy of sqrt{s}=2.76 TeV. The differential production cross-sections times dimuon branching fractions are measured as functions of the Upsilon transverse momentum and rapidity, over the ranges p_T<15 GeV/c and 2.0<y<4.5. The total cross-sections in this kinematic region, assuming unpolarised production, are measured to be sigma(pp -> Upsilon(1S) X) x B(Upsilon(1S) -> mu+mu-) = 1.111 +/- 0.043 +/- 0.044 nb, sigma(pp -> Upsilon(2S) X) x B(Upsilon(2S) -> mu+mu-) = 0.264 +/- 0.023 +/- 0.011 nb, sigma(pp -> Upsilon(3S) X) x B(Upsilon(3S) -> mu+mu-) = 0.159 +/- 0.020 +/- 0.007 nb, where the first uncertainty is statistical and the second systematic.

A measurement of the transverse momentum spectra of jets in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76$ TeV is reported. Jets are reconstructed from charged particles using the anti-$k_{\rm T}$ jet algorithm with jet resolution parameters $R$ of $0.2$ and $0.3$ in pseudo-rapidity $|\eta|<0.5$. The transverse momentum $p_{\rm T}$ of charged particles is measured down to $0.15$ GeV/$c$ which gives access to the low $p_{\rm T}$ fragments of the jet. Jets found in heavy-ion collisions are corrected event-by-event for average background density and on an inclusive basis (via unfolding) for residual background fluctuations and detector effects. A strong suppression of jet production in central events with respect to peripheral events is observed. The suppression is found to be similar to the suppression of charged hadrons, which suggests that substantial energy is radiated at angles larger than the jet resolution parameter $R=0.3$ considered in the analysis. The fragmentation bias introduced by selecting jets with a high $p_{\rm T}$ leading particle, which rejects jets with a soft fragmentation pattern, has a similar effect on the jet yield for central and peripheral events. The ratio of jet spectra with $R=0.2$ and $R=0.3$ is found to be similar in Pb-Pb and simulated PYTHIA pp events, indicating no strong broadening of the radial jet structure in the reconstructed jets with $R<0.3$.

The production of $J/\psi$ mesons in proton-proton collisions at $\sqrt{s}$= 7 tev is studied with the LHCb detector at the LHC. The differential cross-section for prompt $J/\psi$ production is measured as a function of the $J/\psi$ transverse momentum p_{\rm T}\ and rapidity $y$ in the fiducial region $p_{\rm T}\in[0;14]$ GeV/c and $y\in[2.0;4.5]$. The differential cross-section and fraction of \j psi\ from $b$-hadron decays are also measured in the same p_{\rm T}\ and $y$ range s. The analysis is based on a data sample corresponding to an integrated luminos ity of 5.2 pb$^{-1}$. The measured cross-sections integrated over the fiducial region are $10.52\pm 0. 04\pm 1.40^{+1.64}_{-2.20}\mu b$ for $\mathrm{prompt}~J/\psi$ production and $1.14 \pm 0. 01\pm 0.16 \mu b$ for $J/\psi$ from $b$-hadron decays, where the first uncertainty is statistical and the second systematic. The $\mathrm{prompt}~J/\psi$ production cross-section is obtained assuming no $J/\psi$ polarisation and the third error indicates the acceptance uncertainty due to this assumption.

The production of J/psi mesons is studied with the LHCb detector using data from pp collisions at sqrt(s)=2.76 TeV corresponding to an integrated luminosity of 71 nb^-1. The differential cross-section for inclusive J/psi production is measured as a function of its transverse momentum pT. The cross-section in the fiducial region 0<pT<12 GeV/c and rapidity 2.0<y<4.5 is measured to be 5.6 +/- 0.1(stat) +/- 0.4 (syst) mub, with the assumption of unpolarised J/psi production. The fraction of J/psi production from b-hadron decays is measured to be (7.1 +/- 0.6 (stat) +/- 0.7 (syst)) %.

Production of $B_c^+$ mesons in proton-proton collisions at a center-of-mass energy of 8 TeV is studied with data corresponding to an integrated luminosity of $2.0~{\rm fb}^{-1}$ recorded by the LHCb experiment. The ratio of production cross-sections times branching fractions between the $B_c^+\to J/\psi \pi^+$ and $B^+\to J/\psi K^+$ decays is measured as a function of transverse momentum and rapidity in the regions $0 < p_{\rm T} < 20~{\rm GeV}/c$ and $2.0 < y < 4.5$. The ratio in this kinematic range is measured to be $(0.683\pm0.018\pm0.009)\%$, where the first uncertainty is statistical and the second systematic.

The production of $J/\psi$ mesons with rapidity $1.5<y<4.0$ or $-5.0<y<-2.5$ and transverse momentum $p_\mathrm{T}<14 \mathrm{GeV}/c$ is studied with the LHCb detector in proton-lead collisions at a nucleon-nucleon centre-of-mass energy $\sqrt{s_{NN}}=5 \mathrm{TeV}$. The analysis is based on a data sample corresponding to an integrated luminosity of about $1.6 \mathrm{nb}^{-1}$. For the first time the nuclear modification factor and forward-backward production ratio are determined separately for prompt $J/\psi$ mesons and $J/\psi$ from $b$-hadron decays. Clear suppression of prompt $J/\psi$ production with respect to proton-proton collisions at large rapidity is observed, while the production of $J/\psi$ from $b$-hadron decays is less suppressed. These results show good agreement with available theoretical predictions. The measurement shows that cold nuclear matter effects are important for interpretations of the related quark-gluon plasma signatures in heavy-ion collisions.

The differential cross-section for the inclusive production of $\psi(2S)$ mesons in $pp$ collisions at $\sqrt{s}$ =7 TeV has been measured with the LHCb detector. The data sample corresponds to an integrated luminosity of 36 pb$^{-1}$ collected during 2010. The $\psi(2S)$ mesons are reconstructed in the decay channels $\psi(2S) \rightarrow \mu^{+} \mu^{-}$ and $\psi(2S)$ $\rightarrow$ $J$/$\psi \pi^{+}\pi^{-}$, with the $J$/$\psi$ meson decaying into two muons. Results are presented both for promptly produced $\psi(2S)$ mesons and for those originating from $b$-hadron decays. In the kinematic range $p_{\rm T}(\psi(2S))$ $\leq$ 16 GeV/$c$ and 2 < $y$($\psi(2S))$ $\leq$ 4.5 we measure \begin{eqnarray*} \sigma_{\rm prompt}(\psi(2S)) &=& 1.44 \pm 0.01~(\text{stat})\pm 0.12~(\text{syst})^{+0.20}_{-0.40}~(\text{pol})~{\rm \mu b}, \\ \sigma_{b}(\psi(2S)) &=& 0.25 \pm 0.01~(\text{stat}) \pm 0.02~(\text{syst})~{\rm \mu b}, \end{eqnarray*} where the last uncertainty on the prompt cross-section is due to the unknown $\psi(2S)$ polarization. Recent QCD calculations are found to be in good agreement with our measurements. Combining the present result with the LHCb $J$/$\psi$ measurements we determine the inclusive branching fraction \begin{equation*} \mathcal{B}(b \rightarrow \psi(2S) X) = (2.73 \pm 0.06~(\text{stat}) \pm 0.16~(\text{syst}) \pm 0.24~(\text{BF})) \times 10^{-3}, \end{equation*} where the last uncertainty is due to the $\mathcal{B}(b \rightarrow J/\psi X)$, $\mathcal{B}(J/\psi \rightarrow \mu^+ \mu^-)$ and $\mathcal{B}(\psi(2S) \rightarrow e^+ e^-)$ branching fraction uncertainties. All above results are corrected by an erratum included as an appendix.

The production of $\Upsilon$ mesons in $pp$ collisions at $\sqrt=7$ and $8\,\mathrm{TeV}$ is studied with the LHCb detector using data samples corresponding to an integrated luminosity of $1\,\mathrm{fb}^{-1}$ and $2\,\mathrm{fb}^{-1}$ respectively. The production cross-sections and ratios of cross-sections are measured as functions of the meson transverse momentum $p_T$ and rapidity $y$, for $p_T<30\,\mathrm{GeV}/c$} and $2.0<y<4.5$.

A measurement of the $Z(\rightarrow\mu^+\mu^-)$+jet production cross-section in $pp$ collisions at a centre-of-mass energy $\sqrt{s} = 7$ TeV is presented. The analysis is based on an integrated luminosity of $1.0\,\text{fb}^{-1}$ recorded by the LHCb experiment. Results are shown with two jet transverse momentum thresholds, 10 and 20 GeV, for both the overall cross-section within the fiducial volume, and for six differential cross-section measurements. The fiducial volume requires that both the jet and the muons from the Z boson decay are produced in the forward direction ($2.0<\eta<4.5$). The results show good agreement with theoretical predictions at the second-order expansion in the coupling of the strong interaction.

The polarization of the $\Upsilon(1S)$, $\Upsilon(2S)$ and $\Upsilon(3S) $mesons, produced in $pp$ collisions at centre-of-mass energies $\sqrt{s}$=7 and 8TeV, is measured using data samples collected by the LHCb experiment, corresponding to integrated luminosities of 1 and 2fb$^{-1}$, respectively. The measurements are performed in three polarization frames, using $\Upsilon\to\mu^+\mu^-$ decays in the kinematic region of the transverse momentum $p_{T}(\Upsilon)<30GeV/c$, and rapidity $2.2<y(\Upsilon)<4.5$. No large polarization is observed.

The production of J/psi and Upsilon mesons in pp collisions at sqrt(s) = 8 TeV is studied with the LHCb detector. The J/psi and Upsilon mesons are reconstructed in the mu+mu- decay mode and the signal yields are determined with a fit to the mu+mu- invariant mass distributions. The analysis is performed in the rapidity range 2.0<y<4.5 and transverse momentum range 0<p_T<14(15) GeV/c of the J/psi(Upsilon) mesons. The J/psi and Upsilon production cross-sections and the fraction of J/psi mesons from b-hadron decays are measured as a function of the meson p_T and y.

The cross-section for inclusive phi meson production in pp collisions at a centre-of-mass energy of sqrt(s) = 7 TeV has been measured with the LHCb detector at the Large Hadron Collider. The differential cross-section is measured as a function of the phi transverse momentum p_T and rapidity y in the region 0.6 < p_T < 5.0 GeV/c and 2.44 < y < 4.06. The cross-section for inclusive phi production in this kinematic range is sigma(pp -> phi X) = 1758 pm 19(stat) ^{+43}_{-14}(syst) pm 182(scale) microbarn, where the first systematic uncertainty depends on the p_T and y region and the second is related to the overall scale. Predictions based on the Pythia 6.4 generator underestimate the cross-section.

The production of B+- mesons in proton-proton collisions at sqrt(s)=7 TeV is studied using 35 pb-1 of data collected by the LHCb detector. The B+- mesons are reconstructed exclusively in the B+- -&gt; J/psi K+- mode, with J/psi -&gt; mu+ mu-. The differential production cross-section is measured as a function of the B+- transverse momentum in the fiducial region 0 &lt; pT &lt; 40 GeV/c and with rapidity 2.0 &lt; y &lt; 4.5. The total cross-section, summing up B+ and B-, is measured to be sigma(pp -&gt; B+- X, 0 &lt; pT &lt; 40 GeV/c, 2.0 &lt; y &lt; 4.5) = 41.4 +- 1.5 (stat.) +- 3.1 (syst.) mub.

The production cross-sections of $J/\psi$ mesons in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=5$ TeV are measured using a data sample corresponding to an integrated luminosity of $9.13\pm0.18~\text{pb}^{-1}$, collected by the LHCb experiment. The cross-sections are measured differentially as a function of transverse momentum, $p_{\text{T}}$, and rapidity, $y$, and separately for $J/\psi$ mesons produced promptly and from beauty hadron decays (nonprompt). With the assumption of unpolarised $J/\psi$ mesons, the production cross-sections integrated over the kinematic range $0<p_{\text{T}}<20~\text{GeV}/c$ and $2.0<y<4.5$ are $8.154\pm0.010\pm0.283~\mu\text{b}$ for prompt $J/\psi$ mesons and $0.820\pm0.003\pm0.034~\mu\text{b}$ for nonprompt $J/\psi$ mesons, where the first uncertainties are statistical and the second systematic. These cross-sections are compared with those at $\sqrt{s}=8$ TeV and $13$ TeV, and are used to update the measurement of the nuclear modification factor in proton-lead collisions for $J/\psi$ mesons at a centre-of-mass energy per nucleon pair of $\sqrt{s_{\text{NN}}}=5$ TeV. The results are compared with theoretical predictions.

The ALICE collaboration at the CERN Large Hadron Collider reports the first measurement of the inclusive differential jet cross section at mid-rapidity in pp collisions at $\sqrt{s} = 2.76$ TeV, with integrated luminosity of 13.6 nb$^{-1}$. Jets are measured over the transverse momentum range 20 to 125 GeV/c and are corrected to the particle level. Calculations based on Next-to-Leading Order perturbative QCD are in good agreement with the measurements. The ratio of inclusive jet cross sections for jet radii $R = 0.2$ and $R = 0.4$ is reported, and is also well reproduced by a Next-to-Leading Order perturbative QCD calculation when hadronization effects are included.

The production of $J/\psi$ mesons accompanied by open charm, and of pairs of open charm hadrons are observed in pp collisions at a centre-of-mass energy of 7 TeV using an integrated luminosity of $355pb^{-1}$ collected with the LHCb detector. Model independent measurements of absolute cross-sections are given together with ratios to the measured $J/\psi$ and open charm cross-sections. The properties of these events are studied and compared to theoretical predictions.

The production cross-sections of B mesons are measured in pp collisions at a centre-of-mass energy of 7 TeV, using data collected with the LHCb detector corresponding to a integrated luminosity of 0.36 fb-1. The B+, B0 and Bs0 mesons are reconstructed in the exclusive decays B+ -> J/psi K+, B0 -> J/psi K*0 and Bs0 -> J/psi phi, with J/psi -> mu+ mu-, K*0 -> K+ pi- and phi -> K+ K-. The differential cross-sections are measured as functions of B meson transverse momentum pT and rapidity y, in the range 0 < pT < 40 GeV/c and 2.0 < y < 4.5. The integrated cross-sections in the same pT and y ranges, including charge-conjugate states, are measured to be sigma(pp -> B+ + X) = 38.9 +- 0.3 (stat.) +- 2.5 (syst.) +- 1.3 (norm.) mub, sigma(pp -> B0 + X) = 38.1 +- 0.6 (stat.) +- 3.7 (syst.) +- 4.7 (norm.) mub, sigma(pp -> Bs0 + X) = 10.5 +- 0.2 (stat.) +- 0.8 (syst.) +- 1.0 (norm.) mub, where the third uncertainty arises from the pre-existing branching fraction measurements.

A precision measurement of the $Z$ boson production cross-section at $\sqrt{s} = 13$ TeV in the forward region is presented, using $pp$ collision data collected by the LHCb detector, corresponding to an integrated luminosity of 5.1 fb$^{-1}$. The production cross-section is measured using $Z\rightarrow\mu^+\mu^-$ events within the fiducial region defined as pseudorapidity $2.0<\eta<4.5$ and transverse momentum $p_{T}>20$ GeV/$c$ for both muons and dimuon invariant mass $60<M_{\mu\mu}<120$ GeV/$c^2$. The integrated cross-section is determined to be $\sigma (Z \rightarrow \mu^+ \mu^-)$ = 196.4 $\pm$ 0.2 $\pm$ 1.6 $\pm$ 3.9~pb, where the first uncertainty is statistical, the second is systematic, and the third is due to the luminosity determination. The measured results are in agreement with theoretical predictions within uncertainties.

The production of $J/\psi$ mesons in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=13$ TeV is studied with the \lhcb detector. Cross-section measurements are performed as a function of the transverse momentum $p_\mathrm{T}$ and the rapidity $y$ of the $J/\psi$ meson in the region $p_\mathrm{T}<14\mathrm{GeV}/c$ and $2.0<y<4.5$, for both prompt $J/\psi$ mesons and $J/\psi$ mesons from $b$-hadron decays. The production cross-sections integrated over the kinematic coverage are $15.03\pm 0.03\pm 0.94$ $\mu$b for prompt $J/\psi$ and $2.25\pm 0.01\pm 0.14$ $\mu$b for $J/\psi$ from $b$-hadron decays, assuming zero polarization of the $J/\psi$ meson. The first uncertainties are statistical and the second systematic. The cross-section reported for $J/\psi$ mesons from $b$-hadron decays is used to extrapolate to a total $b\bar{b}$ cross-section. The ratios of the cross-sections with respect to $\sqrt{s}=8$ TeV are also determined.

Production cross-sections of prompt charm mesons are measured with the first data from $pp$ collisions at the LHC at a centre-of-mass energy of $13\,\mathrm{TeV}$. The data sample corresponds to an integrated luminosity of $4.98 \pm 0.19\,\mathrm{pb}^{-1}$ collected by the LHCb experiment. The production cross-sections of $D^{0}$, $D^{+}$, $D_{s}^{+}$, and $D^{*+}$ mesons are measured in bins of charm meson transverse momentum, $p_{\mathrm{T}}$, and rapidity, $y$, and cover the range $0 < p_{\mathrm{T}} < 15\,\mathrm{GeV}/c$ and $2.0 < y < 4.5$. The inclusive cross-sections for the four mesons, including charge conjugation, within the range of $1 < p_{\mathrm{T}} < 8\,\mathrm{GeV}/c$ are found to be \begin{equation} \sigma(pp \to D^{0} X) = 2072 \pm 2 \pm 124\,\mu\mathrm{b}\\ \sigma(pp \to D^{+} X) = 834 \pm 2 \pm \phantom{1}78\,\mu\mathrm{b}\\ \sigma(pp \to D_{s}^{+} X) = 353 \pm 9 \pm \phantom{1}76\,\mu\mathrm{b}\\ \sigma(pp \to D^{*+} X) = 784 \pm 4 \pm \phantom{1}87\,\mu\mathrm{b} \end{equation} where the uncertainties are due to statistical and systematic uncertainties, respectively.

Charged-hadron transverse-momentum and pseudorapidity distributions in proton-proton collisions at sqrt(s) = 7 TeV are measured with the inner tracking system of the CMS detector at the LHC. The charged-hadron yield is obtained by counting the number of reconstructed hits, hit-pairs, and fully reconstructed charged-particle tracks. The combination of the three methods gives a charged-particle multiplicity per unit of pseudorapidity, dN(charged)/d(eta), for |eta| < 0.5, of 5.78 +/- 0.01 (stat) +/- 0.23 (syst) for non-single-diffractive events, higher than predicted by commonly used models. The relative increase in charged-particle multiplicity from sqrt(s) = 0.9 to 7 TeV is 66.1% +/- 1.0% (stat) +/- 4.2% (syst). The mean transverse momentum is measured to be 0.545 +/- 0.005 (stat) +/- 0.015 (syst) GeV/c. The results are compared with similar measurements at lower energies.

A measurement of the production cross-section for $Z$ bosons that decay to muons is presented. The data were recorded by the LHCb detector during $pp$ collisions at a centre-of-mass energy of 7 TeV, and correspond to an integrated luminosity of 1.0 fb$^{-1}$. The cross-section is measured for muons in the pseudorapidity range $2.0 < \eta < 4.5$ with transverse momenta $p_{T} > 20$ GeV/c. The dimuon mass is restricted to $60 < M_{\mu^{+}\mu^{-}} < 120$ GeV/c$^{2}$. The measured cross-section is $$\sigma_{Z\rightarrow\mu^{+}\mu^{-}} = (76.0 \pm 0.3 \pm 0.5 \pm 1.0 \pm 1.3) \, \text{pb}$$ where the uncertainties are due to the sample size, systematic effects, the beam energy and the luminosity. This result is in good agreement with theoretical predictions at next-to-next-to-leading order in perturbative quantum chromodynamics. The cross-section is also measured differentially as a function of kinematic variables of the $Z$ boson. Ratios of the production cross-sections of electroweak bosons are presented using updated LHCb measurements of $W$ boson production. A precise test of the Standard Model is provided by the measurement of the ratio $$\frac{\sigma_{W^{+}\rightarrow\mu^{+}\nu_{\mu}} + \sigma_{W^{-}\rightarrow\mu^{-}\bar{\nu}_{\mu}}}{\sigma_{Z\rightarrow\mu^{+}\mu^{-}}} = 20.63\pm0.09\pm0.12\pm0.05,$$ where the uncertainty due to luminosity cancels.

The production of prompt charged particles in proton-lead collisions and in proton-proton collisions at the nucleon-nucleon centre-of-mass energy ${\sqrt{s_{\scriptscriptstyle\mathrm{NN}}}=5\,\mathrm{TeV}}$ is studied at LHCb as a function of pseudorapidity ($\eta$) and transverse momentum ($p_{\mathrm{T}}$) with respect to the proton beam direction. The nuclear modification factor for charged particles is determined as a function of $\eta$ between ${-4.8<\eta<-2.5}$ (backward region) and ${2.0<\eta<4.8}$ (forward region), and $p_{\mathrm{T}}$ between ${0.2<p_{\mathrm{T}}<8.0\,\mathrm{GeV}/c}$. The results show a suppression of charged particle production in proton-lead collisions relative to proton-proton collisions in the forward region and an enhancement in the backward region for $p_{\mathrm{T}}$ larger than $1.5\,\mathrm{GeV}/c$. This measurement constrains nuclear PDFs and saturation models at previously unexplored values of the parton momentum fraction down to $10^{-6}$.

Differential production cross sections of J/$\psi$ and $\psi$(2S) charmonium and $\Upsilon$(nS) (n = 1, 2, 3) bottomonium states are measured in proton-proton collisions at $\sqrt{s} =$ 13 TeV, with data collected by the CMS detector at the LHC, corresponding to an integrated luminosity of 2.3 fb$^{-1}$ for the J/$\psi$ and 2.7 fb$^{-1}$ for the other mesons. The five quarkonium states are reconstructed in the dimuon decay channel, for dimuon rapidity $|y| <$ 1.2. The double-differential cross sections for each state are measured as a function of $y$ and transverse momentum, and compared to theoretical expectations. In addition, ratios are presented of cross sections for prompt $\psi$(2S) to J/$\psi$, $\Upsilon$(2S) to $\Upsilon$(1S), and $\Upsilon$(3S) to $\Upsilon$(1S) production.

Charged particle multiplicities are studied in proton-proton collisions in the forward region at a centre-of-mass energy of $\sqrt{s} = 7\;$TeV with data collected by the LHCb detector. The forward spectrometer allows access to a kinematic range of $2.0<\eta<4.8$ in pseudorapidity, momenta down to $2\;$GeV/$c$ and transverse momenta down to $0.2\;$GeV/$c$. The measurements are performed using minimum-bias events with at least one charged particle in the kinematic acceptance. The results are presented as functions of pseudorapidity and transverse momentum and are compared to predictions from several Monte Carlo event generators.

Production cross-sections of prompt charm mesons are measured using data from $pp$ collisions at the LHC at a centre-of-mass energy of $5\,$TeV. The data sample corresponds to an integrated luminosity of $8.60\pm0.33\,$pb$^{-1}$ collected by the LHCb experiment. The production cross-sections of $D^0$, $D^+$, $D_s^+$, and $D^{*+}$ mesons are measured in bins of charm meson transverse momentum, $p_{\text{T}}$, and rapidity, $y$. They cover the rapidity range $2.0<y<4.5$ and transverse momentum ranges $0 < p_{\text{T}} < 10\, \text{GeV}/c$ for $D^0$ and $D^+$ and $1 < p_{\text{T}} < 10\, \text{GeV}/c$ for $D_s^+$ and $D^{*+}$ mesons. The inclusive cross-sections for the four mesons, including charge-conjugate states, within the range of $1 < p_{\text{T}} < 8\, \text{GeV}/c$ are determined to be \sigma(pp\rightarrow D^0 X) = 1004 \pm 3 \pm 54\,\mu\text{b} \sigma(pp\rightarrow D^+ X) = 402 \pm 2 \pm 30\,\mu\text{b} \sigma(pp\rightarrow D_s^+ X) = 170 \pm 4 \pm 16\,\mu\text{b} \sigma(pp\rightarrow D^{*+} X)= 421 \pm 5 \pm 36\,\mu\text{b} where the uncertainties are statistical and systematic, respectively.

Measurements are presented of electroweak boson production using data from $pp$ collisions at a centre-of-mass energy of $\sqrt{s} = 8\mathrm{\,Te\kern -0.1em V}$. The analysis is based on an integrated luminosity of $2.0\mathrm{\,fb}^{-1}$ recorded with the LHCb detector. The bosons are identified in the $W\rightarrow\mu\nu$ and $Z\rightarrow\mu^{+}\mu^{-}$ decay channels. The cross-sections are measured for muons in the pseudorapidity range $2.0 < \eta < 4.5$, with transverse momenta $p_{\rm T} > 20{\mathrm{\,Ge\kern -0.1em V\!/}c}$ and, in the case of the $Z$ boson, a dimuon mass within $60 < M_{\mu^{+}\mu^{-}} < 120{\mathrm{\,Ge\kern -0.1em V\!/}c^{2}}$. The results are \begin{align*} \sigma_{W^{+}\rightarrow\mu^{+}\nu} &= 1093.6 \pm 2.1 \pm 7.2 \pm 10.9 \pm 12.7{\rm \,pb} \, , \sigma_{W^{-}\rightarrow\mu^{-}\bar{\nu}} &= \phantom{0}818.4 \pm 1.9 \pm 5.0 \pm \phantom{0}7.0 \pm \phantom{0}9.5{\rm \,pb} \, , \sigma_{Z\rightarrow\mu^{+}\mu^{-}} &= \phantom{00}95.0 \pm 0.3 \pm 0.7 \pm \phantom{0}1.1 \pm \phantom{0}1.1{\rm \,pb} \, , \end{align*} where the first uncertainties are statistical, the second are systematic, the third are due to the knowledge of the LHC beam energy and the fourth are due to the luminosity determination. The evolution of the $W$ and $Z$ boson cross-sections with centre-of-mass energy is studied using previously reported measurements with $1.0\mathrm{\,fb}^{-1}$ of data at $7\mathrm{\,Te\kern -0.1em V}$. Differential distributions are also presented. Results are in good agreement with theoretical predictions at next-to-next-to-leading order in perturbative quantum chromodynamics.

The inclusive b-jet production cross section in pp collisions at a center-of-mass energy of 7 TeV is measured using data collected by the CMS experiment at the LHC. The cross section is presented as a function of the jet transverse momentum in the range 18 < pT < 200 GeV for several rapidity intervals. The results are also given as the ratio of the b-jet production cross section to the inclusive jet production cross section. The measurement is performed with two different analyses, which differ in their trigger selection and b-jet identification: a jet analysis that selects events with a b jet using a sample corresponding to an integrated luminosity of 34 inverse picobarns, and a muon analysis requiring a b jet with a muon based on an integrated luminosity of 3 inverse picobarns. In both approaches the b jets are identified by requiring a secondary vertex. The results from the two methods are in agreement with each other and with next-to-leading order calculations, as well as with predictions based on the PYTHIA event generator.

The spectra of strange hadrons are measured in proton-proton collisions, recorded by the CMS experiment at the CERN LHC, at centre-of-mass energies of 0.9 and 7 TeV. The K^0_s, Lambda, and Xi^- particles and their antiparticles are reconstructed from their decay topologies and the production rates are measured as functions of rapidity and transverse momentum. The results are compared to other experiments and to predictions of the PYTHIA Monte Carlo program. The transverse momentum distributions are found to differ substantially from the PYTHIA results and the production rates exceed the predictions by up to a factor of three.

The production of the $\eta_c (1S)$ state in proton-proton collisions is probed via its decay to the $p \bar{p}$ final state with the LHCb detector, in the rapidity range $2.0 < y < 4.5$ and in the meson transverse-momentum range $p_T > 6.5$ GeV/c. The cross-section for prompt production of $\eta_c (1S)$ mesons relative to the prompt $J/\psi$ cross-section is measured, for the first time, to be $\sigma_{\eta_c (1S)}/\sigma_{J/\psi} = 1.74 \pm 0.29 \pm 0.28 \pm 0.18 _{B}$ at a centre-of-mass energy $\sqrt{s} = 7$ TeV using data corresponding to an integrated luminosity of 0.7 fb$^{-1}$, and $\sigma_{\eta_c (1S)}/\sigma_{J/\psi} = 1.60 \pm 0.29 \pm 0.25 \pm 0.17 _{B}$ at $\sqrt{s} = 8$ TeV using 2.0 fb$^{-1}$. The uncertainties quoted are, in order, statistical, systematic, and that on the ratio of branching fractions of the $\eta_c (1S)$ and $J/\psi$ decays to the $p \bar{p}$ final state. In addition, the inclusive branching fraction of $b$-hadron decays into $\eta_c (1S)$ mesons is measured, for the first time, to be $B ( b \rightarrow \eta_c X ) = (4.88 \pm 0.64 \pm 0.29 \pm 0.67 _{B}) \times 10^{-3}$, where the third uncertainty includes also the uncertainty on the $J/\psi$ inclusive branching fraction from $b$-hadron decays. The difference between the $J/\psi$ and $\eta_c (1S)$ meson masses is determined to be $114.7 \pm 1.5 \pm 0.1$ MeV/c$^2$.

A measurement of the production cross-section of Z bosons in pp collisions at $\sqrt{s} = 13$ TeV is presented using dimuon and dielectron final states in LHCb data. The cross-section is measured for leptons with pseudorapidities in the range $2.0 < \eta < 4.5$, transverse momenta $p_\text{T} > 20$ GeV and dilepton invariant mass in the range $60<m(\ell\ell)<120$ GeV. The integrated cross-section from averaging the two final states is \begin{equation*}\sigma_{\text{Z}}^{\ell\ell} = 194.3 \pm 0.9 \pm 3.3 \pm 7.6\text{ pb,}\end{equation*} where the first uncertainty is statistical, the second is due to systematic effects, and the third is due to the luminosity determination. In addition, differential cross-sections are measured as functions of the Z boson rapidity, transverse momentum and the angular variable $\phi^*_\eta$.

Associated production of bottomonia and open charm hadrons in $pp$ collisions at $\sqrt{s}=7$ and $8$TeV is observed using data corresponding to an integrated luminosity of 3$fb^{-1}$ accumulated with the LHCb detector. The observation of five combinations, $\Upsilon(1S)D^0$, $\Upsilon(2S)D^0$, $\Upsilon(1S)D^+$, $\Upsilon(2S)D^+$ and $\Upsilon(1S)D^+_{s}$, is reported. Production cross-sections are measured for $\Upsilon(1S)D^0$ and $\Upsilon(1S)D^+$ pairs in the forward region. The measured cross-sections and the differential distributions indicate the dominance of double parton scattering as the main production mechanism. This allows a precise measurement of the effective cross-section for double parton scattering.

The inclusive jet cross section is measured in pp collisions with a center-of-mass energy of 7 TeV at the LHC using the CMS experiment. The data sample corresponds to an integrated luminosity of 34 inverse picobarns. The measurement is made for jet transverse momenta in the range 18-1100 GeV and for absolute values of rapidity less than 3. The measured cross section extends to the highest values of jet pT ever observed and, within the experimental and theoretical uncertainties, is generally in agreement with next-to-leading-order perturbative QCD predictions.

Measurements of the differential production cross sections in transverse momentum and rapidity for B0 mesons produced in pp collisions at sqrt(s) = 7 TeV are presented. The dataset used was collected by the CMS experiment at the LHC and corresponds to an integrated luminosity of 40 inverse picobarns. The production cross section is measured from B0 meson decays reconstructed in the exclusive final state J/Psi K-short, with the subsequent decays J/Psi to mu^+ mu^- and K-short to pi^+ pi^-. The total cross section for pt(B0) > 5 GeV and y(B0) < 2.2 is measured to be 33.2 +/- 2.5 +/- 3.5 microbarns, where the first uncertainty is statistical and the second is systematic.

Production of prompt J/$\psi$ meson pairs in proton-proton collisions at $\sqrt{s}$ = 7 TeV is measured with the CMS experiment at the LHC in a data sample corresponding to an integrated luminosity of about 4.7 inverse-femtobarns. The two J/$\psi$ mesons are fully reconstructed via their decays into $\mu^+\mu^-$ pairs. This observation provides for the first time access to the high-transverse-momentum region of J/$\psi$ pair production where model predictions are not yet established. The total and differential cross sections are measured in a phase space defined by the individual J/$\psi$ transverse momentum ($p_T^{J/\psi}$) and rapidity (|$y^{J/\psi}$|): |$y^{J/\psi}$| lower than 1.2 for $p_T^{J/\psi}$ greater than 6.5 GeV/c; |$y^{J/\psi}$| in [1.2,1.43] for a $p_T$ threshold that scales linearly with |$y^{J/\psi}$| from 6.5 to 4.5 GeV/c; and |$y^{J/\psi}$| in [1.43,2.2] for $p_T^{J/\psi}$ greater than 4.5 GeV/c. The total cross section, assuming unpolarized prompt J/$\psi$ pair production is 1.49 $\pm$ 0.07 (stat.) $\pm$ 0.13 (syst.) nb. Different assumptions about the J/$\psi$ polarization imply modifications to the cross section ranging from -31% to +27%.