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

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.

A study is presented of central exclusive production of $\Upsilon(nS)$ states, where the $\Upsilon(nS)$ resonances decay to the $\mu^+\mu^-$ final state, using $pp$ collision data recorded by the LHCb experiment. The cross-section is measured in the rapidity range $2<y(\Upsilon)<4.5$ where the muons are reconstructed in the pseudorapidity range $2<\eta(\mu^\pm)<4.5$. The data sample corresponds to an integrated luminosity of 2.9 fb$^{-1}$ and was collected at centre-of-mass energies of $7$ TeV and $8$ TeV. The measured $\Upsilon(1S)$ and $\Upsilon(2S)$ production cross-sections are \begin{eqnarray} \sigma(pp \to p\Upsilon(1S)p) &=& 9.0 \pm 2.1 \pm 1.7\textrm{ pb and}\nonumber\\ \sigma(pp \to p\Upsilon(2S)p) &=& 1.3 \pm 0.8 \pm 0.3\textrm{ pb},\nonumber \end{eqnarray} where the first uncertainties are statistical and the second are systematic. The $\Upsilon(1S)$ cross-section is also measured as a function of rapidity and is found to be in good agreement with Standard Model predictions. An upper limit is set at 3.4 pb at the 95\% confidence level for the exclusive $\Upsilon(3S)$ production cross-section, including possible contamination from $\chi_b(3P)\to\Upsilon(3S)\gamma$ decays.

Associated production of bottomonia and open charm hadrons in $pp$ collisions at $\sqrt{s}=7$ and $8$TeV is observed using data corresponding to an integrated luminosity of 3$fb^{-1}$ accumulated with the LHCb detector. The observation of five combinations, $\Upsilon(1S)D^0$, $\Upsilon(2S)D^0$, $\Upsilon(1S)D^+$, $\Upsilon(2S)D^+$ and $\Upsilon(1S)D^+_{s}$, is reported. Production cross-sections are measured for $\Upsilon(1S)D^0$ and $\Upsilon(1S)D^+$ pairs in the forward region. The measured cross-sections and the differential distributions indicate the dominance of double parton scattering as the main production mechanism. This allows a precise measurement of the effective cross-section for double parton scattering.

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.

Using 1.8 fb-1 of pp collisions at a center-of-mass energy of 7 TeV recorded by the ATLAS detector at the Large Hadron Collider, we present measurements of the production cross sections of Upsilon(1S,2S,3S) mesons. Upsilon mesons are reconstructed using the di-muon decay mode. Total production cross sections for p_T<70 GeV and in the rapidity interval |Upsilon|<2.25 are measured to be 8.01+-0.02+-0.36+-0.31 nb, 2.05+-0.01+-0.12+-0.08 nb, 0.92+-0.01+-0.07+-0.04 nb respectively, with uncertainties separated into statistical, systematic, and luminosity measurement effects. In addition, differential cross section times di-muon branching fractions for Upsilon(1S), Upsilon(2S), and Upsilon(3S) as a function of Upsilon transverse momentum p_T and rapidity are presented. These cross sections are obtained assuming unpolarized production. If the production polarization is fully transverse or longitudinal with no azimuthal dependence in the helicity frame the cross section may vary by approximately +-20%. If a non-trivial azimuthal dependence is considered, integrated cross sections may be significantly enhanced by a factor of two or more. We compare our results to several theoretical models of Upsilon meson production, finding that none provide an accurate description of our data over the full range of Upsilon transverse momenta accessible with this dataset.

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.