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.

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 polarisation of prompt $\psi(2S)$ mesons is measured by performing an angular analysis of $\psi(2S)\rightarrow \mu^{+} \mu^{-}$ decays using proton-proton collision data, corresponding to an integrated luminosity of 1.0 fb$^{-1}$, collected by the LHCb detector at a centre-of-mass energy of 7 TeV. The polarisation is measured in bins of transverse momentum $p_\mathrm{T}$ and rapidity $y$ in the kinematic region $3.5<p_\mathrm{T}<15$ GeV$/c$ and $2.0<y<4.5$, and is compared to theoretical models. No significant polarisation is observed.

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.

A measurement of the inclusive $W \to \mu\nu$ production cross-section using data from $pp$ collisions at a centre-of-mass energy of $\sqrt{s} = 7 {\rm \, TeV}$ is presented. The analysis is based on an integrated luminosity of about $1.0 {\rm \, fb}^{-1}$ recorded with the LHCb detector. Results are reported for muons with a transverse momentum greater than $20 {\rm \, Ge\kern -0.1em V\!/}c$ and pseudorapidity between 2.0 and 4.5. The $W^{+}$ and $W^{-}$ production cross-sections are measured to be $$\sigma_{W^{+} \to \mu^{+}\nu} = 861.0 \pm 2.0 \pm 11.2 \pm 14.7 {\rm \, pb},$$ $$\sigma_{W^{-} \to \mu^{-}\overline{\nu}} = 675.8 \pm 1.9 \pm 8.8 \pm 11.6 {\rm \, pb},$$ where the first uncertainty is statistical, the second is systematic and the third is due to the luminosity determination. Cross-section ratios and differential distributions as functions of the muon pseudorapidity are also presented. The ratio of $W^{+}$ to $W^{-}$ cross-sections in the same fiducial kinematic region is determined to be $$\frac{\sigma_{W^{+} \to \mu^{+}\nu}}{\sigma_{W^{-} \to \mu^{-}\overline{\nu}}} = 1.274 \pm 0.005 \pm 0.009,$$ where the uncertainties are statistical and systematic, respectively. Results are in good agreement with theoretical predictions at next-to-next-to-leading order in perturbative quantum chromodynamics.

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.

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 of $\Upsilon$ mesons in proton-lead collisions at a nucleon-nucleon centre-of-mass energy $\sqrt{s_{NN}}=5 \mathrm{TeV}$ is studied with the LHCb detector. The analysis is based on a data sample corresponding to an integrated luminosity of $1.6 \mathrm{nb}^{-1}$. The $\Upsilon$ mesons of transverse momenta up to $15 \mathrm{GeV}/c$ are reconstructed in the dimuon decay mode. The rapidity coverage in the centre-of-mass system is $1.5 < y < 4.0$ (forward region) and $-5.0 < y < -2.5$ (backward region). The forward-backward production ratio and the nuclear modification factor for $\Upsilon(1S)$ mesons are determined. The data are compatible with the predictions for a suppression of $\Upsilon(1S)$ production with respect to proton-proton collisions in the forward region, and an enhancement in the backward region. The suppression is found to be smaller than in the case of prompt $J/\psi$ mesons.

The differential cross-section of prompt inclusive production of long-lived charged particles in proton-proton collisions is measured using a data sample recorded by the LHCb experiment at a centre-of-mass energy of ${\sqrt{s} = 13\,\mathrm{TeV}}$. The data sample, collected with an unbiased trigger, corresponds to an integrated luminosity of ${5.4\,\mathrm{nb}^{-1}}$. The differential cross-section is measured as a function of transverse momentum and pseudorapidity in the ranges ${p_\mathrm{T} \in [0.08, 10)\,\mathrm{GeV}\,c^{-1}}$ and ${\eta \in [2.0, 4.8)}$ and is determined separately for positively and negatively charged particles. The results are compared with predictions from various hadronic-interaction models.

The differential cross-section as a function of rapidity has been measured for the exclusive production of $J/\psi$ and $\psi(2S)$ mesons in proton-proton collisions at $\sqrt{s}=7$ TeV, using data collected by the LHCb experiment, corresponding to an integrated luminosity of 930 pb$^{-1}$. The cross-sections times branching fractions to two muons having pseudorapidities between 2.0 and 4.5 are measured to be $$\begin{array}{rl} \sigma_{pp\rightarrow J/\psi\rightarrow{\mu^+}{\mu^-}}(2.0<\eta_{\mu^\pm }<4.5)=&291\pm 7\pm19 {\rm \ pb},\\ \sigma_{pp\rightarrow\psi(2S)\rightarrow{\mu^+}{\mu^-}}(2.0<\eta_{\mu^\pm}<4.5)=&6.5\pm 0.9\pm 0.4 {\rm \ pb},\end{array}$$ where the first uncertainty is statistical and the second is systematic. The measurements agree with next-to-leading order QCD predictions as well as with models that include saturation effects.

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 study of $\chi_{b}$ meson production at LHCb is performed on proton-proton collision data, corresponding to 3.0fb$^{-1}$ of integrated luminosity collected at centre-of-mass energies $\sqrt{s}=7$ and 8TeV. The fraction of $\Upsilon(nS)$ mesons originating from $\chi_{b}$ decays is measured as a function of the $\Upsilon$ transverse momentum in the rapidity range $2.0 < y^{\Upsilon} < 4.5$. The radiative transition of the $\chi_{b}(3P)$ meson to $\Upsilon(3S)$ is observed for the first time. The $\chi_{b1}(3P)$ mass is determined to be $$m(\chi_{b1}(3P)) = 10\,511.3 \pm 1.7 \pm 2.5 MeV/c^2,$$ where the first uncertainty is statistical and the second is systematic.

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