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.

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

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.

The production cross-sections of B mesons are measured in pp collisions at a centre-of-mass energy of 7 TeV, using data collected with the LHCb detector corresponding to a integrated luminosity of 0.36 fb-1. The B+, B0 and Bs0 mesons are reconstructed in the exclusive decays B+ -> J/psi K+, B0 -> J/psi K*0 and Bs0 -> J/psi phi, with J/psi -> mu+ mu-, K*0 -> K+ pi- and phi -> K+ K-. The differential cross-sections are measured as functions of B meson transverse momentum pT and rapidity y, in the range 0 < pT < 40 GeV/c and 2.0 < y < 4.5. The integrated cross-sections in the same pT and y ranges, including charge-conjugate states, are measured to be sigma(pp -> B+ + X) = 38.9 +- 0.3 (stat.) +- 2.5 (syst.) +- 1.3 (norm.) mub, sigma(pp -> B0 + X) = 38.1 +- 0.6 (stat.) +- 3.7 (syst.) +- 4.7 (norm.) mub, sigma(pp -> Bs0 + X) = 10.5 +- 0.2 (stat.) +- 0.8 (syst.) +- 1.0 (norm.) mub, where the third uncertainty arises from the pre-existing branching fraction measurements.

The production of J/psi mesons is studied with the LHCb detector using data from pp collisions at sqrt(s)=2.76 TeV corresponding to an integrated luminosity of 71 nb^-1. The differential cross-section for inclusive J/psi production is measured as a function of its transverse momentum pT. The cross-section in the fiducial region 0<pT<12 GeV/c and rapidity 2.0<y<4.5 is measured to be 5.6 +/- 0.1(stat) +/- 0.4 (syst) mub, with the assumption of unpolarised J/psi production. The fraction of J/psi production from b-hadron decays is measured to be (7.1 +/- 0.6 (stat) +/- 0.7 (syst)) %.

The production of $J/\psi$ mesons in proton-proton collisions at $\sqrt{s}$= 7 tev is studied with the LHCb detector at the LHC. The differential cross-section for prompt $J/\psi$ production is measured as a function of the $J/\psi$ transverse momentum p_{\rm T}\ and rapidity $y$ in the fiducial region $p_{\rm T}\in[0;14]$ GeV/c and $y\in[2.0;4.5]$. The differential cross-section and fraction of \j psi\ from $b$-hadron decays are also measured in the same p_{\rm T}\ and $y$ range s. The analysis is based on a data sample corresponding to an integrated luminos ity of 5.2 pb$^{-1}$. The measured cross-sections integrated over the fiducial region are $10.52\pm 0. 04\pm 1.40^{+1.64}_{-2.20}\mu b$ for $\mathrm{prompt}~J/\psi$ production and $1.14 \pm 0. 01\pm 0.16 \mu b$ for $J/\psi$ from $b$-hadron decays, where the first uncertainty is statistical and the second systematic. The $\mathrm{prompt}~J/\psi$ production cross-section is obtained assuming no $J/\psi$ polarisation and the third error indicates the acceptance uncertainty due to this assumption.

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 production of Upsilon(1S), Upsilon(2S) and Upsilon(3S) mesons in proton-proton collisions at the centre-of-mass energy of sqrt(s)=7 TeV is studied with the LHCb detector. The analysis is based on a data sample of 25 pb-1 collected at the Large Hadron Collider. The Upsilon mesons are reconstructed in the decay mode Upsilon -&gt; mu+ mu- and the signal yields are extracted from a fit to the mu+ mu- invariant mass distributions. The differential production cross-sections times dimuon branching fractions are measured as a function of the Upsilon transverse momentum pT and rapidity y, over the range pT &lt; 15 GeV/c and 2.0 &lt; y &lt; 4.5. The cross-sections times branching fractions, integrated over these kinematic ranges, are measured to be sigma(pp -&gt; Upsilon(1S) X) x B(Upsilon(1S)-&gt;mu+ mu-) = 2.29 {\pm} 0.01 {\pm} 0.10 -0.37 +0.19 nb, sigma(pp -&gt; Upsilon(2S) X) x B(Upsilon(2S)-&gt;mu+ mu-) = 0.562 {\pm} 0.007 {\pm} 0.023 -0.092 +0.048 nb, sigma(pp -&gt; Upsilon(3S) X) x B(Upsilon(3S)-&gt;mu+ mu-) = 0.283 {\pm} 0.005 {\pm} 0.012 -0.048 +0.025 nb, where the first uncertainty is statistical, the second systematic and the third is due to the unknown polarisation of the three Upsilon states.

The interpretation of cosmic antiproton flux measurements from space-borne experiments is currently limited by the knowledge of the antiproton production cross-section in collisions between primary cosmic rays and the interstellar medium. Using collisions of protons with an energy of 6.5 TeV incident on helium nuclei at rest in the proximity of the interaction region of the LHCb experiment, the ratio of antiprotons originating from antihyperon decays to prompt production is measured for antiproton momenta between 12 and 110 GeV/c. The dominant antihyperon contribution, namely $\bar{\Lambda} \to \bar{p} \pi^+$ decays from promptly produced $\bar{\Lambda}$ particles, is also exclusively measured. The results complement the measurement of prompt antiproton production obtained from the same data sample. At the energy scale of this measurement, the antihyperon contributions to antiproton production are observed to be significantly larger than predictions of commonly used hadronic production models.

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.

The charged particle production in proton-proton collisions is studied with the LHCb detector at a centre-of-mass energy of ${\sqrt{s} =7}$TeV in different intervals of pseudorapidity $\eta$. The charged particles are reconstructed close to the interaction region in the vertex detector, which provides high reconstruction efficiency in the $\eta$ ranges $-2.5&lt;\eta&lt;-2.0$ and $2.0&lt;\eta&lt;4.5$. The data were taken with a minimum bias trigger, only requiring one or more reconstructed tracks in the vertex detector. By selecting an event sample with at least one track with a transverse momentum greater than 1 GeV/c a hard QCD subsample is investigated. Several event generators are compared with the data; none are able to describe fully the multiplicity distributions or the charged particle density distribution as a function of $\eta$. In general, the models underestimate the charged particle production.

The charged-particle production ratios $\bar{p}/p$, $K^-/K^+$, $\pi^-/\pi^+$, $(p + \bar{p})/(\pi^+ + \pi^-)$, $(K^+ + K^-)/(\pi^+ + \pi^-)$ and $(p + \bar{p})/(K^+ + K^-)$ are measured with the LHCb detector using $0.3 {\rm nb^{-1}}$ of $pp$ collisions delivered by the LHC at $\sqrt{s} = 0.9$ TeV and $1.8 {\rm nb^{-1}}$ at $\sqrt{s} = 7$ TeV. The measurements are performed as a function of transverse momentum $p_{\rm T}$ and pseudorapidity $\eta$. The production ratios are compared to the predictions of several Monte Carlo generator settings, none of which are able to describe adequately all observables. The ratio $\bar{p}/p$ is also considered as a function of rapidity loss, $\Delta y \equiv y_{\rm beam} - y$, and is used to constrain models of baryon transport.

The production of B+- mesons in proton-proton collisions at sqrt(s)=7 TeV is studied using 35 pb-1 of data collected by the LHCb detector. The B+- mesons are reconstructed exclusively in the B+- -&gt; J/psi K+- mode, with J/psi -&gt; mu+ mu-. The differential production cross-section is measured as a function of the B+- transverse momentum in the fiducial region 0 &lt; pT &lt; 40 GeV/c and with rapidity 2.0 &lt; y &lt; 4.5. The total cross-section, summing up B+ and B-, is measured to be sigma(pp -&gt; B+- X, 0 &lt; pT &lt; 40 GeV/c, 2.0 &lt; y &lt; 4.5) = 41.4 +- 1.5 (stat.) +- 3.1 (syst.) mub.

A measurement of the cross-section for pp$ \rightarrow$Z$ \rightarrow$e$^+$e$^-$ is presented using data at $\sqrt{s}=7$ TeV corresponding to an integrated luminosity of 0.94 fb$^{-1}$. The process is measured within the kinematic acceptance $p_{\mathrm{T}}&gt;20$GeV/$c$ and $2&lt;\eta&lt;4.5$ for the daughter electrons and dielectron invariant mass in the range 60-120 GeV/$c^2$. The cross-section is determined to be $$\sigma(pp \rightarrow Z \rightarrow e^+ e^- )=76.0\pm0.8\pm2.0\pm2.6{\rm pb}$$ where the first uncertainty is statistical, the second is systematic and the third is the uncertainty in the luminosity. The measurement is performed as a function of Z rapidity and as a function of an angular variable which is closely related to the Z transverse momentum. The results are compared with previous LHCb measurements and with theoretical predictions from QCD.

The energy flow created in pp collisions at is studied within the pseudorapidity range 1.9<η<4.9 with data collected by the LHCb experiment. The measurements are performed for inclusive minimum-bias interactions, hard scattering processes and events with an enhanced or suppressed diffractive contribution. The results are compared to predictions given by Pythia-based and cosmic-ray event generators, which provide different models of soft hadronic interactions.

The production of \chib1P mesons in $pp$ collisions at a centre-of-mass energy of $7\tev$ is studied using $32\invpb$ of data collected with the \lhcb detector. The $\chib1P$ mesons are reconstructed in the decay mode $\chib1P \to \Y1S\g \to \mumu\g$. The fraction of \Y1S originating from \chib1P decays in the \Y1S transverse momentum range $6 < \pt^{\Y1S} < 15\gevc$ and rapidity range $2.0 < y^{\Y1S} < 4.5$ is measured to be $(20.7\pm 5.7\pm 2.1^{+2.7}_{-5.4})%$, where the first uncertainty is statistical, the second is systematic and the last gives the range of the result due to the unknown \Y1S and \chib1P polarizations.

The cross-section for inclusive phi meson production in pp collisions at a centre-of-mass energy of sqrt(s) = 7 TeV has been measured with the LHCb detector at the Large Hadron Collider. The differential cross-section is measured as a function of the phi transverse momentum p_T and rapidity y in the region 0.6 < p_T < 5.0 GeV/c and 2.44 < y < 4.06. The cross-section for inclusive phi production in this kinematic range is sigma(pp -> phi X) = 1758 pm 19(stat) ^{+43}_{-14}(syst) pm 182(scale) microbarn, where the first systematic uncertainty depends on the p_T and y region and the second is related to the overall scale. Predictions based on the Pythia 6.4 generator underestimate the cross-section.

The prompt production of charmonium $\chi_{c}$ and $J/\psi$ states is studied in proton-proton collisions at a centre-of-mass energy of $\sqrt{s}=7$ TeV at the Large Hadron Collider. The $\chi_{c}$ and $J/\psi$ mesons are identified through their decays $\chi_{c}\rightarrow J/\psi \gamma$ and $J/\psi\rightarrow \mu^+\mu^-$ using 36 pb$^{-1}$ of data collected by the LHCb detector in 2010. The ratio of the prompt production cross-sections for $\chi_{c}$ and $J/\psi$, $\sigma (\chi_{c}\rightarrow J/\psi \gamma)/ \sigma (J/\psi)$, is determined as a function of the $J/\psi$ transverse momentum in the range $2 &lt; p_{\mathrm T}^{J/\psi} &lt; 15$ GeV/$c$. The results are in excellent agreement with next-to-leading order non-relativistic expectations and show a significant discrepancy compared with the colour singlet model prediction at leading order, especially in the low $p_{\mathrm T}^{J/\psi}$ region.

The production of $J/\psi$ pairs in proton-proton collisions at a centre-of-mass energy of 7 TeV has been observed using an integrated luminosity of $37.5 pb^{-1}$ collected with the LHCb detector. The production cross-section for pairs with both \jpsi in the rapidity range $2&lt;y^{J/\psi}&lt;4.5$ and transverse momentum $p_{T}^{J/\psi}&lt;10 GeV/c$ is $$ \sigma^{J/\psi J/\psi} = 5.1\pm1.0\pm1.1 nb,$$ where the first uncertainty is statistical and the second systematic.

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.

The production of $J/\psi$ mesons accompanied by open charm, and of pairs of open charm hadrons are observed in pp collisions at a centre-of-mass energy of 7 TeV using an integrated luminosity of $355pb^{-1}$ collected with the LHCb detector. Model independent measurements of absolute cross-sections are given together with ratios to the measured $J/\psi$ and open charm cross-sections. The properties of these events are studied and compared to theoretical predictions.