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 -> 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 < 15 GeV/c and 2.0 < y < 4.5. The cross-sections times branching fractions, integrated over these kinematic ranges, are measured to be sigma(pp -> Upsilon(1S) X) x B(Upsilon(1S)->mu+ mu-) = 2.29 {\pm} 0.01 {\pm} 0.10 -0.37 +0.19 nb, sigma(pp -> Upsilon(2S) X) x B(Upsilon(2S)->mu+ mu-) = 0.562 {\pm} 0.007 {\pm} 0.023 -0.092 +0.048 nb, sigma(pp -> Upsilon(3S) X) x B(Upsilon(3S)->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.

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.

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 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 in proton-lead collisions at the centre-of-mass energy per nucleon pair $\sqrt{s_{\text{NN}}}=8.16$ TeV with the LHCb detector at the LHC. The double differential cross-sections of prompt and nonprompt J/$\psi$ production are measured as functions of the J/$\psi$ transverse momentum and rapidity in the nucleon-nucleon centre-of-mass frame. Forward-to-backward ratios and nuclear modification factors are determined. The results are compared with theoretical calculations based on collinear factorisation using nuclear parton distribution functions, on the colour glass condensate or on coherent energy loss models.

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.

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.

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