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.

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.

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.

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 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 first observation of $Z$ boson production in proton-lead collisions at a centre-of-mass energy per proton-nucleon pair of $\sqrt{s_{NN}}=5~\text{TeV}$ is presented. The data sample corresponds to an integrated luminosity of $1.6~\text{nb}^{-1}$ collected with the LHCb detector. The $Z$ candidates are reconstructed from pairs of oppositely charged muons with pseudorapidities between 2.0 and 4.5 and transverse momenta above $20~\text{GeV}/c$. The invariant dimuon mass is restricted to the range $60-120~\text{GeV}/c^2$. The $Z$ production cross-section is measured to be \begin{eqnarray*} \sigma_{Z\to\mu^+\mu^-}(\text{fwd})&=&13.5^{+5.4}_{-4.0}\text{(stat.)}\pm1.2\text{(syst.)}~\text{nb} \end{eqnarray*} in the direction of the proton beam and \begin{eqnarray*} \sigma_{Z\to\mu^+\mu^-}(\text{bwd}) & =&10.7^{+8.4}_{-5.1}\text{(stat.)}\pm1.0\text{(syst.)}~\text{nb} \end{eqnarray*} in the direction of the lead beam, where the first uncertainty is statistical and the second systematic.

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