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.

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 production cross-sections of $\Upsilon(1S)$, $\Upsilon(2S)$ and $\Upsilon(3S)$ mesons in proton-proton collisions at $\sqrt{s}$= 13 TeV are measured with a data sample corresponding to an integrated luminosity of $277 \pm 11$ $\rm pb^{-1}$ recorded by the LHCb experiment in 2015. The $\Upsilon$ mesons are reconstructed in the decay mode $\Upsilon\to\mu^{+}\mu^{-}$. The differential production cross-sections times the dimuon branching fractions are measured as a function of the $\Upsilon$ transverse momentum, $p_{\rm T}$, and rapidity, $y$, over the range $0 < p_{\rm T}< 30$ GeV/c and $2.0 < y < 4.5$. The ratios of the cross-sections with respect to the LHCb measurement at $\sqrt{s}$= 8 TeV are also determined. The measurements are compared with theoretical predictions based on NRQCD.

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.

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

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 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 $J/\psi$ mesons with rapidity $1.5<y<4.0$ or $-5.0<y<-2.5$ and transverse momentum $p_\mathrm{T}<14 \mathrm{GeV}/c$ is studied with the LHCb detector in proton-lead collisions at a nucleon-nucleon centre-of-mass energy $\sqrt{s_{NN}}=5 \mathrm{TeV}$. The analysis is based on a data sample corresponding to an integrated luminosity of about $1.6 \mathrm{nb}^{-1}$. For the first time the nuclear modification factor and forward-backward production ratio are determined separately for prompt $J/\psi$ mesons and $J/\psi$ from $b$-hadron decays. Clear suppression of prompt $J/\psi$ production with respect to proton-proton collisions at large rapidity is observed, while the production of $J/\psi$ from $b$-hadron decays is less suppressed. These results show good agreement with available theoretical predictions. The measurement shows that cold nuclear matter effects are important for interpretations of the related quark-gluon plasma signatures in heavy-ion collisions.

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

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.

The production of J/psi and Upsilon mesons in pp collisions at sqrt(s) = 8 TeV is studied with the LHCb detector. The J/psi and Upsilon mesons are reconstructed in the mu+mu- decay mode and the signal yields are determined with a fit to the mu+mu- invariant mass distributions. The analysis is performed in the rapidity range 2.0<y<4.5 and transverse momentum range 0<p_T<14(15) GeV/c of the J/psi(Upsilon) mesons. The J/psi and Upsilon production cross-sections and the fraction of J/psi mesons from b-hadron decays are measured as a function of the meson p_T and y.

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

Measurements of the cross-section for producing \bquark quarks in the reaction $pp\to b\bar{b} X$ are reported in 7 and 13 TeV collisions at the LHC as a function of the pseudorapidity $\eta$ in the range $2<\eta<5$ covered by the acceptance of the LHCb experiment. The measurements are done using semileptonic decays of $b$-flavored hadrons decaying into a ground-state charmed hadron in association with a muon. The cross-sections in the covered $\eta$ range are $72.0\pm 0.3\pm6.8~\mu$b and $144\pm 1\pm 21~\mu$b for 7 and 13 TeV. The ratio is $2.00\pm0.02\pm0.26$, where the quoted uncertainties are statistical and systematic, respectively. The agreement with theoretical expectation is good at 7 TeV, but differs somewhat at 13 TeV. The measured ratio of cross-sections is larger at lower $\eta$ than the model prediction.

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.