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 production of $W$ and $Z$ bosons in association with jets is studied in the forward region of proton-proton collisions collected at a centre-of-mass energy of 8 TeV by the LHCb experiment, corresponding to an integrated luminosity of 1.98 $\pm$ 0.02 fb$^{-1}$. The $W$ boson is identified using its decay to a muon and a neutrino, while the $Z$ boson is identified through its decay to a muon pair. Total cross-sections are measured and combined into charge ratios, asymmetries, and ratios of $W+$jet and $Z$+jet production cross-sections. Differential measurements are also performed as a function of both boson and jet kinematic variables. All results are in agreement with Standard Model predictions.

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

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$ mesons is studied in proton-lead collisions at the centre-of-mass energy per nucleon pair $\sqrt{s_{\text{NN}}}=8.16$ TeV with the LHCb detector at the LHC. The double differential cross-sections of prompt and nonprompt J/$\psi$ production are measured as functions of the J/$\psi$ transverse momentum and rapidity in the nucleon-nucleon centre-of-mass frame. Forward-to-backward ratios and nuclear modification factors are determined. The results are compared with theoretical calculations based on collinear factorisation using nuclear parton distribution functions, on the colour glass condensate or on coherent energy loss models.

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

The product of the $\Lambda_b^0$ ($\overline{B}^0$) differential production cross-section and the branching fraction of the decay $\Lambda_b^0\rightarrow J/\psi pK^-$ ($\overline{B}^0\rightarrow J/\psi\overline{K}^*(892)^0$) is measured as a function of the beauty hadron transverse momentum, $p_{\rm T}$, and rapidity, $y$. The kinematic region of the measurements is $p_{\rm T}<20~{\rm GeV}/c$ and $2.0<y<4.5$. The measurements use a data sample corresponding to an integrated luminosity of $3~{\rm fb}^{-1}$ collected by the LHCb detector in $pp$ collisions at centre-of-mass energies $\sqrt{s}=7~{\rm TeV}$ in 2011 and $\sqrt{s}=8~{\rm TeV}$ in 2012. Based on previous LHCb results of the fragmentation fraction ratio, $f_{\Lambda_B^0}/f_d$, the branching fraction of the decay $\Lambda_b^0\rightarrow J/\psi pK^-$ is measured to be \begin{equation*} \mathcal{B}(\Lambda_b^0\rightarrow J/\psi pK^-)= (3.17\pm0.04\pm0.07\pm0.34^{+0.45}_{-0.28})\times10^{-4}, \end{equation*} where the first uncertainty is statistical, the second is systematic, the third is due to the uncertainty on the branching fraction of the decay $\overline{B}^0\rightarrow J/\psi\overline{K}^*(892)^0$, and the fourth is due to the knowledge of $f_{\Lambda_b^0}/f_d$. The sum of the asymmetries in the production and decay between $\Lambda_b^0$ and $\overline{\Lambda}_b^0$ is also measured as a function of $p_{\rm T}$ and $y$. The previously published branching fraction of $\Lambda_b^0\rightarrow J/\psi p\pi^-$, relative to that of $\Lambda_b^0\rightarrow J/\psi pK^-$, is updated. The branching fractions of $\Lambda_b^0\rightarrow P_c^+(\rightarrow J/\psi p)K^-$ are determined.

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