Exclusive $J/\psi$ and $\psi(2S)$ vector meson production has been observed in the dimuon channel using the LHCb detector. The cross-section times branching fractions to two muons with pseudorapidities between 2.0 and 4.5 are measured to be \begin{equation*} \sigma_{pp\rightarrow J/\psi (\rightarrow \mu^{+} \mu^{-})} (2.0 <\eta_{\mu^{\pm}}< 4.5) = 307 \pm 21 \pm 36~\text{pb}, \end{equation*} \begin{equation*} \sigma_{pp\rightarrow \psi(2S) (\rightarrow \mu^{+} \mu^{-})} (2.0 <\eta_{\mu^{\pm}}< 4.5) = 7.8 \pm 1.3 \pm 1.0~\text{pb}, \end{equation*} where the first uncertainties are statistical and the second are systematic. The measurements are found to be in good agreement with results from previous experiments and theoretical predictions. The $J/\psi$ photoproduction cross-section has been measured as a function of the photon-proton centre-of-mass energy. The results are consistent with measurements obtained at HERA and confirm a similar power law behaviour for the photoproduction cross-section.

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

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

An angular analysis of the decay $J/\psi \rightarrow \mu^{+} \mu^{-}$ is performed to measure the polarization of prompt $J/\psi$ mesons produced in $pp$ collisions at $\sqrt{s}=7\; \mathrm{TeV}$. The dataset corresponds to an integrated luminosity of 0.37 fb$^{-1}$ collected with the LHCb detector. The measurement is presented as a function of transverse momentum, $p_\mathrm{T}$, and rapidity, $y$, of the $J/\psi$ meson, in the kinematic region $2 < p_\mathrm{T} < 15\; \mathrm{GeV}/c$ and $2.0<y<4.5$.

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 $\bar{\Lambda} / \Lambda$ and $\bar{\Lambda} / K^0_\mathrm{S}$ production ratios are measured by the LHCb detector from $0.3\,\mathrm{nb}^{-1}$ of $pp$ collisions delivered by the LHC at $\sqrt{s} = 0.9$\,TeV and $1.8\,\mathrm{nb}^{-1}$ at $\sqrt{s} = 7$\,TeV. Both ratios are presented as a function of transverse momentum, $p_\mathrm{T}$, and rapidity, $y$, in the ranges {$0.15 < p_\mathrm{T} < 2.50\,\mathrm{GeV}/c$} and {$2.0<y<4.5$}. Results at the two energies are in good agreement as a function of rapidity loss, $\Delta y = y_\mathrm{beam} - y$, and are consistent with previous measurements. The ratio $\bar{\Lambda} / \Lambda$, measuring the transport of baryon number from the collision into the detector, is smaller in data than predicted in simulation, particularly at high rapidity. The ratio $\bar{\Lambda} / K^0_\mathrm{S}$, measuring the baryon-to-meson suppression in strange quark hadronisation, is significantly larger than expected.

The first measurement of the $Z$ boson production cross-section at centre-of-mass energy $\sqrt{s} = 5.02\,$TeV in the forward region is reported, using $pp$ collision data collected by the LHCb experiment in year 2017, corresponding to an integrated luminosity of $100 \pm 2\,\rm{pb^{-1}}$. The production cross-section is measured for final-state muons in the pseudorapidity range $2.0<\eta<4.5$ with transverse momentum $p_{\rm{T}}> 20\,\rm{GeV/}\it{c}$. The integrated cross-section is determined to be \[ \sigma_{Z \rightarrow \mu^{+}\mu^{-}} = 39.6 \pm 0.7\,(\rm{stat}) \pm 0.6\,(\rm{syst}) \pm 0.8\,(\rm{lumi}) \ \rm{pb} \] for the di-muon invariant mass in the range $60<M_{\mu\mu}<120\,\rm{GeV/}\it{c^{2}}$. This result and the differential cross-section results are in good agreement with theoretical predictions at next-to-next-to-leading order in the strong coupling. Based on a previous LHCb measurement of the $Z$ boson production cross-section in $p$Pb collisions at $\sqrt{s_{NN}}=5.02$ TeV, the nuclear modification factor $R_{p\rm{Pb}}$ is measured for the first time at this energy. The measured values are $1.2^{+0.5}_{-0.3}\,(\rm{stat}) \pm 0.1\,(\rm{syst})$ in the forward region ($1.53<y^*_{\mu}<4.03$) and $3.6^{+1.6}_{-0.9}\,(\rm{stat}) \pm 0.2\,(\rm{syst})$ in the backward region ($-4.97<y^*_{\mu}<-2.47$), where $y^*_{\mu}$ represents the muon rapidity in the centre-of-mass frame.

The production cross-sections of $\Upsilon(1S)$, $\Upsilon(2S)$ and $\Upsilon(3S)$ mesons in proton-proton collisions at $\sqrt{s}$= 13 TeV are measured with a data sample corresponding to an integrated luminosity of $277 \pm 11$ $\rm pb^{-1}$ recorded by the LHCb experiment in 2015. The $\Upsilon$ mesons are reconstructed in the decay mode $\Upsilon\to\mu^{+}\mu^{-}$. The differential production cross-sections times the dimuon branching fractions are measured as a function of the $\Upsilon$ transverse momentum, $p_{\rm T}$, and rapidity, $y$, over the range $0 < p_{\rm T}< 30$ GeV/c and $2.0 < y < 4.5$. The ratios of the cross-sections with respect to the LHCb measurement at $\sqrt{s}$= 8 TeV are also determined. The measurements are compared with theoretical predictions based on NRQCD.

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.