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

We report on measurements of the ϒ(1S), ϒ(2S), and ϒ(3S) differential, (d2σdPtdy)y=0, and integrated cross sections in pp¯ collisions at s=1.8 TeV using a sample of 16.6 ± 0.6 pb−1 collected by the Collider Detector at Fermilab. The three resonances were reconstructed through the decay ϒ→μ+μ−. Comparison is made to a leading order QCD prediction.

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

Differential production cross sections of J/$\psi$ and $\psi$(2S) charmonium and $\Upsilon$(nS) (n = 1, 2, 3) bottomonium states are measured in proton-proton collisions at $\sqrt{s} =$ 13 TeV, with data collected by the CMS detector at the LHC, corresponding to an integrated luminosity of 2.3 fb$^{-1}$ for the J/$\psi$ and 2.7 fb$^{-1}$ for the other mesons. The five quarkonium states are reconstructed in the dimuon decay channel, for dimuon rapidity $|y| <$ 1.2. The double-differential cross sections for each state are measured as a function of $y$ and transverse momentum, and compared to theoretical expectations. In addition, ratios are presented of cross sections for prompt $\psi$(2S) to J/$\psi$, $\Upsilon$(2S) to $\Upsilon$(1S), and $\Upsilon$(3S) to $\Upsilon$(1S) production.

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

We report on measurements of the ϒ(1S), ϒ(2S), and ϒ(3S) differential cross sections (d2σ/dpTdy)|y|<0.4, as well as on the ϒ(1S) polarization in pp¯ collisions at s=1.8TeV using a sample of 77±3pb−1 collected by the collider detector at Fermilab. The three resonances were reconstructed through the decay ϒ→μ+μ−. The measured angular distribution of the muons in the ϒ(1S) rest frame is consistent with unpolarized meson production.

The proton-nucleon cross section ratio $R=Br(\Upsilon\to l^+l^-) d\sigma(\Upsilon)/dy|_{y=0} / {\sigma(J/\psi)}$ has been measured with the HERA-B spectrometer in fixed-target proton-nucleus collisions at 920 GeV proton beam energy corresponding to a proton-nucleon cms energy of sqrt{s}=41.6 GeV. The combined results for the Upsilon decay channels Upsilon $\to e^+e^-$ and Upsilon $\to\mu^+\mu^-$ yield a ratio $R=(9.0 \pm 2.1) 10^{-6}$. The corresponding Upsilon production cross section per nucleon at mid-rapidity (y=0) has been determined to be $Br(\Upsilon\to{}l^+l^-) {d\sigma(\Upsilon)/dy}|_{y=0}= 4.5 \pm 1.1 $ pb/nucleon.

We report measurements of Upsilon meson production in p+p, d+Au, and Au+Au collisions using the STAR detector at RHIC. We compare the Upsilon yield to the measured cross section in p+p collisions in order to quantify any modifications of the yield in cold nuclear matter using d+Au data and in hot nuclear matter using Au+Au data separated into three centrality classes. Our p+p measurement is based on three times the statistics of our previous result. We obtain a nuclear modification factor for Upsilon(1S+2S+3S) in the rapidity range |y|<1 in d+Au collisions of R_dAu = 0.79 +/- 0.24 (stat.) +/- 0.03 (sys.) +/- 0.10 (pp sys.). A comparison with models including shadowing and initial state parton energy loss indicates the presence of additional cold-nuclear matter suppression. Similarly, in the top 10% most-central Au+Au collisions, we measure a nuclear modification factor of R_AA=0.49 +/- 0.1 (stat.) +/- 0.02 (sys.) +/- 0.06 (pp sys.), which is a larger suppression factor than that seen in cold nuclear matter. Our results are consistent with complete suppression of excited-state Upsilon mesons in Au+Au collisions. The additional suppression in Au+Au is consistent with the level expected in model calculations that include the presence of a hot, deconfined Quark-Gluon Plasma. However, understanding the suppression seen in d+Au is still needed before any definitive statements about the nature of the suppression in Au+Au can be made.

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.

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.

The production yields of Upsilon(1S), Upsilon(2S), and Upsilon(3S) quarkonium states are measured through their decays into muon pairs in the CMS detector, in PbPb and pp collisions at the centre-of-mass energy per nucleon pair of 2.76 TeV. The data correspond to integrated luminosities of 166 inverse microbarns and 5.4 inverse picobarns for PbPb and pp collisions, respectively. Differential production cross sections are reported as functions of Upsilon rapidity y up to 2.4, and transverse momentum pT up to 20 GeV/c. A strong centrality-dependent suppression is observed in PbPb relative to pp collisions, by factors of up to approximately 2 and 8, for the Upsilon(1S) and Upsilon(2S) states, respectively. No significant dependence of this suppression is observed as a function of y or pT. The Upsilon(3S) state is not observed in PbPb collisions, which corresponds to a suppression for the centrality-integrated data by at least a factor of approximately 7 at a 95% confidence level. The observed suppression is in agreement with theoretical scenarios modeling the sequential melting of quarkonium states in a quark gluon plasma.

We have reconstructed the radiative decays $\chi_{b}(1P) \to \Upsilon(1S) \gamma $ and $\chi_{b}(2P) \to \Upsilon(1S) \gamma $ in $p \bar{p}$ collisions at $\sqrt{s} = 1.8$ TeV, and measured the fraction of $\Upsilon(1S)$ mesons that originate from these decays. For $\Upsilon(1S)$ mesons with $p^{\Upsilon}_{T}>8.0$ GeV/$c$, the fractions that come from $\chi_{b}(1P)$ and $\chi_{b}(2P)$ decays are $(27.1\pm6.9(stat)\pm4.4(sys))%$ and $(10.5\pm4.4(stat)\pm1.4(sys))%$, respectively. We have derived the fraction of directly produced $\Upsilon(1S)$ mesons to be $(50.9\pm8.2(stat)\pm9.0(sys))%$.

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.

Differential cross sections as a function of transverse momentum pt are presented for the production of Y(nS) (n = 1, 2, 3) states decaying into a pair of muons. Data corresponding to an integrated luminosity of 4.9 inverse femtobarns in pp collisions at sqrt(s) = 7 TeV were collected with the CMS detector at the LHC. The analysis selects events with dimuon rapidity abs(y) < 1.2 and dimuon transverse momentum in the range 10 < pt < 100 GeV. The measurements show a transition from an exponential to a power-law behavior at pt ~ 20 GeV for the three Y states. Above that transition, the Y spectrum is significantly harder than that of the Y(1S) and Y(2S). The ratios of the Y(3S) and Y(2S) differential cross sections to the Y(1S) cross section show a rise as pt increases at low pt, then become flatter at higher pt.

We present measurements of the inclusive production cross sections of the Upsilon(1S) bottomonium state in ppbar collisions at sqrt(s) = 1.96 TeV. Using the Upsilon(1S) to mu+mu- decay mode for a data sample of 159 +- 10 pb^-1 collected by the D0 detector at the Fermilab Tevatron collider, we determine the differential cross sections as a function of the Upsilon(1S) transverse momentum for three ranges of the Upsilon(1S) rapidity: 0 < |y| < 0.6, 0.6 < |y| < 1.2, and 1.2 < |y| < 1.8.

The production of $\Upsilon$ mesons in Pb-Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{\rm NN}}$ = 5 TeV is measured with the muon spectrometer of the ALICE detector at the LHC. The yields as well as the nuclear modification factors are determined in the forward rapidity region $2.5<y<4.0$, as a function of rapidity, transverse momentum and collision centrality. The results show that the production of the $\Upsilon$(1S) meson is suppressed by a factor of about three with respect to the production in proton-proton collisions. For the first time, a significant signal for the $\Upsilon$(2S) meson is observed at forward rapidity, indicating a suppression stronger by about a factor 2-3 with respect to the ground state. The measurements are compared with transport, hydrodynamic, comover and statistical hadronisation model calculations.

The production of $\Upsilon$(2S) and $\Upsilon$(3S) mesons in lead-lead (PbPb) and proton-proton (pp) collisions is studied in their dimuon decay channel using the CMS detector at the LHC. The $\Upsilon$(3S) meson is observed for the first time in PbPb collisions, with a significance above five standard deviations. The ratios of yields measured in PbPb and pp collisions are reported for both the $\Upsilon$(2S) and $\Upsilon$(3S) mesons, as functions of transverse momentum and PbPb collision centrality. These ratios, when appropriately scaled, are significantly less than unity, indicating a suppression of $\Upsilon$ yields in PbPb collisions. This suppression increases from peripheral to central PbPb collisions. Furthermore, the suppression is stronger for $\Upsilon$(3S) mesons compared to $\Upsilon$(2S) mesons, extending the pattern of sequential suppression of quarkonium states in nuclear collisions previously seen for the $\psi$/J, $\psi$(2S), $\Upsilon$(1S), and $\Upsilon$(2S) mesons.

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.

Measurements of the cross section for the production of electron pairs with invariant masses between 4 and 8.7 GeV are presented as a function of the centre-of-mass energy ( s = 28 to s = 62 GeV ) of the colliding proton beams. A significant excess of events is observed in the region 8.7 to 10.3 GeV; these are ascribed to the ϒ(9.5 GeV) resonances and estimates of the production cross sections are given.

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.