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 crosssections and the fraction of J/psi mesons from bhadron decays are measured as a function of the meson p_T and y.
Projections of the fit result for a selected bin in $p_{\rm T}$ and $y$ for (a) the $ J/\psi$ invariant dimuon mass and (b) $t_z$. For the former, the total fitted function is shown (blue solid line) together with the signal distribution (red dotted line). In the $t_z$ projection the total fitted function is shown together with the $ J/\psi \mathrm{from} b $ component, the prompt signal, the background and the tail component due to the association of a $ J/\psi$ candidate with a wrong PV. 
Fig1a.pdf [19 KiB] HiDef png [215 KiB] Thumbnail [59 KiB] *.C file 

Fig1b.pdf [27 KiB] HiDef png [807 KiB] Thumbnail [121 KiB] *.C file 

Invariant mass distribution of the selected $\varUpsilon\rightarrow \mu^+\mu^$ candidates in the range $ p_{\rm T} <15 {\mathrm{ Ge V /}c} $ and $2.0<y<4.5$. The three peaks correspond to the $\varUpsilon(1S)$, $\varUpsilon(2S)$ and $\varUpsilon(3S)$ meson signals (from left to right). The superimposed curve and the signal yields (dotted) are the result of the fit described in the text. 
Fig2.pdf [91 KiB] HiDef png [219 KiB] Thumbnail [58 KiB] *.C file 

Differential production crosssection for (a) prompt $ J/\psi$ mesons and (b) $ J/\psi \mathrm{from} b $ as a function of $p_{\rm T}$ in bins of $y$. It is assumed that prompt $ J/\psi$ mesons are produced unpolarised. The errors are the quadratic sums of the statistical and systematic uncertainties. 
Fig3a.pdf [7 KiB] HiDef png [139 KiB] Thumbnail [58 KiB] *.C file 

Fig3b.pdf [8 KiB] HiDef png [144 KiB] Thumbnail [61 KiB] *.C file 

Differential production crosssection as a function of $y$ integrated over $p_{\rm T}$ , for (a) unpolarised $\mathrm{prompt} J/\psi$ mesons and (b) $ J/\psi \mathrm{from} b $ . The errors are the quadratic sums of the statistical and systematic uncertainties. 
Fig4a_[..].pdf [36 KiB] HiDef png [85 KiB] Thumbnail [38 KiB] *.C file 

Fig4b.pdf [4 KiB] HiDef png [83 KiB] Thumbnail [37 KiB] *.C file 

Fraction of $ J/\psi \mathrm{from} b $ as a function of $p_{\rm T}$ , in bins of $y$. 
Fig5.pdf [18 KiB] HiDef png [188 KiB] Thumbnail [77 KiB] *.C file 

Doubledifferential crosssections times dimuon branching fractions as a function of $p_{\rm T}$ in bins of $y$ for (a) the $\varUpsilon(1S)$, (b) $\varUpsilon(2S)$ and (c) $\varUpsilon(3S)$ mesons. 
Fig6a.pdf [54 KiB] HiDef png [230 KiB] Thumbnail [85 KiB] *.C file 

Fig6b.pdf [54 KiB] HiDef png [222 KiB] Thumbnail [85 KiB] *.C file 

Fig6c.pdf [54 KiB] HiDef png [224 KiB] Thumbnail [83 KiB] *.C file 

Differential production crosssections for $\varUpsilon(1S)$, $\varUpsilon(2S)$ and $\varUpsilon(3S)$ mesons times dimuon branching fraction (a) as a function of $ p_{\rm T} $ integrated over $y$, and (b) as a function of $y$ integrated over $ p_{\rm T} $. 
Fig7a.pdf [17 KiB] HiDef png [124 KiB] Thumbnail [56 KiB] *.C file 

Fig7b.pdf [14 KiB] HiDef png [92 KiB] Thumbnail [43 KiB] *.C file 

{\it } Ratio of the $\varUpsilon(2S)$ to $\varUpsilon(1S)$ and $\varUpsilon(3S)$ to $\varUpsilon(1S)$ crosssections times dimuon branching fractions (a) as a function of $ p_{\rm T} $ integrated over $y$, and (b) as a function of $y$ integrated over $ p_{\rm T} $. 
Fig8a.pdf [16 KiB] HiDef png [132 KiB] Thumbnail [56 KiB] *.C file 

Fig8b.pdf [14 KiB] HiDef png [114 KiB] Thumbnail [49 KiB] *.C file 

Comparison of the differential crosssection for the production of prompt $ J/\psi$ meson (under the assumption of zero polarisation) as a function of $ p_{\rm T} $ with direct production in an NLO NRQCD model [47,48] (orange diagonal shading), an NNLO* CSM [10] (solid yellow) and an NLO CSM [7] (blue vertical shading). The points show the measurements reported in this analysis. 
Fig9_g[..].pdf [14 KiB] HiDef png [318 KiB] Thumbnail [83 KiB] *.C file 

Differential production crosssection for $ J/\psi \mathrm{from} b $ (a) as a function of $p_{\rm T}$ in the fiducial range $2.0<y<4.5$, and (b) as a function of $y$ in the fiducial range $ p_{\rm T} <14 {\mathrm{ Ge V /}c} $. The FONLL prediction [46,50] is shown in yellow. The points show the measurements reported in this analysis. 
Fig10a.pdf [5 KiB] HiDef png [141 KiB] Thumbnail [53 KiB] *.C file 

Fig10b.pdf [4 KiB] HiDef png [129 KiB] Thumbnail [48 KiB] *.C file 

Predictions based on the FONLL formalism [46,50] for the production crosssection for $ J/\psi \mathrm{from} b $ in the fiducial range $0< p_{\rm T} <14 {\mathrm{ Ge V /}c} $ and $2.0<y<4.5$ (yellow band). The uncertainty includes contributions from the renormalisation scale, quark masses and the choice of PDF set. The black dotted line shows the central value of the prediction. The points show the LHCb measurements at $\sqrt{s}=2.76$ [14], 7 [12], and 8$ \mathrm{ Te V} $. 
Fig11.pdf [47 KiB] HiDef png [137 KiB] Thumbnail [50 KiB] *.C file 

Comparison of the differential production crosssections times dimuon branching fractions for (a) $\varUpsilon(1S)$, (b) $\varUpsilon(2S)$ and (c) $\varUpsilon(3S)$ mesons as a function of $ p_{\rm T} $ with direct production in an NNLO$^*$ CSM [9] (solid yellow) and an NLO CSM [7] (blue vertical shading) model. The points show the measurements reported in this analysis. 
Fig12a.pdf [15 KiB] HiDef png [743 KiB] Thumbnail [86 KiB] *.C file 

Fig12b.pdf [15 KiB] HiDef png [749 KiB] Thumbnail [86 KiB] *.C file 

Fig12c.pdf [15 KiB] HiDef png [741 KiB] Thumbnail [85 KiB] *.C file 

Animated gif made out of all figures. 
PAPER2013016.gif Thumbnail 
Selection criteria for the $ J/\psi$ and $\varUpsilon$ meson analyses. Criteria common to both analyses are displayed between the two columns. 
Table_1.pdf [49 KiB] HiDef png [56 KiB] Thumbnail [8 KiB] tex code 

Relative systematic uncertainties (in %) on the $ J/\psi$ and $\varUpsilon$ crosssection results and on the fraction of $ J/\psi \mathrm{from} b $ . 
Table_2.pdf [65 KiB] HiDef png [162 KiB] Thumbnail [23 KiB] tex code 

Differential production crosssection ${\rm d}\sigma/{\rm d}y$ in nb for $\mathrm{prompt} J/\psi$ mesons (assumed unpolarised) and for $ J/\psi \mathrm{from} b $ , integrated over $p_{\rm T}$ . The first uncertainty is statistical, the second (third) is the part of the systematic uncertainty that is uncorrelated (correlated) between bins. 
Table_3.pdf [28 KiB] HiDef png [59 KiB] Thumbnail [10 KiB] tex code 

Doubledifferential crosssection $\frac{{\rm d}^2\sigma}{{\rm d}p_{\rm T}{\rm d}y}$ in nb/( $ {\mathrm{ Ge V /}c}$ ) for prompt $ J/\psi$ meson production in bins of of $p_{\rm T}$ and $y$, with the assumption of no polarisation. The first error is statistical, the second is the component of the systematic uncertainty that is uncorrelated between bins and the third is the correlated component. 
Table_4.pdf [36 KiB] HiDef png [268 KiB] Thumbnail [41 KiB] tex code 

Doubledifferential crosssection $\frac{{\rm d}^2\sigma}{{\rm d}p_{\rm T}{\rm d}y}$ in nb/( $ {\mathrm{ Ge V /}c}$ ) for the production of $ J/\psi \mathrm{from} b $ in bins of $p_{\rm T}$ and $y$. The first error is statistical, the second is the component of the systematic uncertainty that is uncorrelated between bins and the third is the correlated component. 
Table_5.pdf [36 KiB] HiDef png [282 KiB] Thumbnail [46 KiB] tex code 

Fraction of $ J/\psi \mathrm{from} b $ (in %) in bins of of $p_{\rm T}$ and $y$. The first uncertainty is statistical and the second systematic (uncorrelated between bins). 
Table_6.pdf [35 KiB] HiDef png [278 KiB] Thumbnail [58 KiB] tex code 

Doubledifferential production crosssections $\frac{{\rm d}^2\sigma}{{\rm d}p_{\rm T}{\rm d}y}\times \mathcal{B}^{1S}$ in pb/( $ {\mathrm{ Ge V /}c}$ ) for the $\varUpsilon(1S)$ meson in bins of transverse momentum and rapidity, assuming no polarisation. The first error is statistical, the second is the component of the systematic uncertainty that is uncorrelated between bins and the third is the correlated component. 
Table_7.pdf [36 KiB] HiDef png [261 KiB] Thumbnail [46 KiB] tex code 

Doubledifferential production crosssections $\frac{{\rm d}^2\sigma}{{\rm d}p_{\rm T}{\rm d}y}\times \mathcal{B}^{2S}$ in pb/( $ {\mathrm{ Ge V /}c}$ ) for the $\varUpsilon(2S)$ meson in bins of transverse momentum and rapidity, assuming no polarisation. The first error is statistical, the second is the component of the systematic uncertainty that is uncorrelated between bins and the third is the correlated component. Regions where the number of events was not large enough to perform a measurement are indicated with a dash. 
Table_8.pdf [36 KiB] HiDef png [240 KiB] Thumbnail [39 KiB] tex code 

Doubledifferential production crosssections $\frac{{\rm d}^2\sigma}{{\rm d}p_{\rm T}{\rm d}y}\times \mathcal{B}^{3S}$ in pb/( $ {\mathrm{ Ge V /}c}$ ) for the $\varUpsilon(3S)$ meson in bins of transverse momentum and rapidity, assuming no polarisation. The first error is statistical, the second is the component of the systematic uncertainty that is uncorrelated between bins and the third is the correlated component. Regions where the number of events was not large enough to perform a measurement are indicated with a dash. 
Table_9.pdf [36 KiB] HiDef png [234 KiB] Thumbnail [37 KiB] tex code 

Ratios of crosssections $\varUpsilon(2S)\rightarrow\mu^+\mu^$ and $\varUpsilon(3S)\rightarrow\mu^+\mu^$ with respect to {$\varUpsilon(1S)\rightarrow\mu^+\mu^$} as a function of $ p_{\rm T} $ in the range $2.0<y<4.0$, assuming no polarisation. The first error is statistical, the second is the component of the systematic uncertainty that is uncorrelated between bins and the third is the correlated component. 
Table_10.pdf [43 KiB] HiDef png [180 KiB] Thumbnail [29 KiB] tex code 

Ratios of crosssections $\varUpsilon(2S)\rightarrow\mu^+\mu^$ and $\varUpsilon(3S)\rightarrow\mu^+\mu^$ with respect to {$\varUpsilon(1S)\rightarrow\mu^+\mu^$} as a function of $y$ in the range $2.0<y<4.0$, assuming no polarisation. The first error is statistical, the second is the component of the systematic uncertainty that is uncorrelated between bins and the third is the correlated component. 
Table_11.pdf [41 KiB] HiDef png [63 KiB] Thumbnail [10 KiB] tex code 
This ZIP file contains supplemetary material for the publication LHCbPAPER2013016. The files are the three formats of the Figure 13, whose description is below. Fig.13: Comparison of the differential crosssection for the production of prompt J/psi as a function of y with direct production in a NLO NRQCD model 
Fig13.pdf [5 KiB] HiDef png [457 KiB] Thumbnail [97 KiB] *C file 
Created on 14 June 2019.Citation count from INSPIRE on 14 June 2019.