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.

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.

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.

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

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.

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.

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

The production of $\chi_b$ mesons in proton-proton collisions is studied using a data sample collected by the LHCb detector, at centre-of-mass energies of $\sqrt{s}=7$ and $8$ TeV and corresponding to an integrated luminosity of 3.0 fb$^{-1}$. The $\chi_b$ mesons are identified through their decays to $\Upsilon(1S)\gamma$ and $\Upsilon(2S)\gamma$ using photons that converted to $e^+e^-$ pairs in the detector. The relative prompt production rate of $\chi_{b1}(1P)$ and $\chi_{b2}(1P)$ mesons is measured as a function of the $\Upsilon(1S)$ transverse momentum in the $\chi_b$ rapidity range $2.0<y<4.5$. A precise measurement of the $\chi_b(3P)$ mass is also performed. Assuming a mass splitting between the $\chi_{b1}(3P)$ and the $\chi_{b2}(3P)$ states of 10.5 MeV/$c^2$, the measured mass of the $\chi_{b1}(3P)$ meson is $ m(\chi_{b1}(3P))= 10515.7^{+2.2}_{-3.9}(stat) ^{+1.5}_{-2.1}(syst) MeV/c^2.$

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.