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 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 Upsilon production cross section in proton-proton collisions at sqrt(s) = 7 TeV is measured using a data sample collected with the CMS detector at the LHC, corresponding to an integrated luminosity of 3.1 +/- 0.3 inverse picobarns. Integrated over the rapidity range |y|<2, we find the product of the Upsilon(1S) production cross section and branching fraction to dimuons to be sigma(pp to Upsilon(1S) X) B(Upsilon(1S) to mu+ mu-) = 7.37 +/- 0.13^{+0.61}_{-0.42}\pm 0.81 nb, where the first uncertainty is statistical, the second is systematic, and the third is associated with the estimation of the integrated luminosity of the data sample. This cross section is obtained assuming unpolarized Upsilon(1S) production. If the Upsilon(1S) production polarization is fully transverse or fully longitudinal the cross section changes by about 20%. We also report the measurement of the Upsilon(1S), Upsilon(2S), and Upsilon(3S) differential cross sections as a function of transverse momentum and rapidity.

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.

The ratios of the production cross sections between the excited $\Upsilon$(2S) and $\Upsilon$(3S) mesons and the $\Upsilon$(1S) ground state, detected via their decay into two muons, are studied as a function of the number of charged particles in the event. The data are from proton-proton collisions at $\sqrt{s} =$ 7 TeV, corresponding to an integrated luminosity of 4.8 fb$^{-1}$, collected with the CMS detector at the LHC. Evidence of a decrease in these ratios as a function of the particle multiplicity is observed, more pronounced at low transverse momentum $p_\mathrm{T}^{\mu\mu}$. For $\Upsilon$(nS) mesons with $p_\mathrm{T}^{\mu\mu}$ $\gt$ 7 GeV, where most of the data were collected, the correlation with multiplicity is studied as a function of the underlying event transverse sphericity and the number of particles in a cone around the $\Upsilon$(nS) direction. The ratios are found to be multiplicity independent for jet-like events. The mean $p_\mathrm{T}^{\mu\mu}$ values for the $\Upsilon$(nS) states as a function of particle multiplicity are also measured and found to grow more steeply as their mass increases.