Measurement of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _c(1S)$$\end{document}ηc(1S) production cross-section in proton–proton collisions via the decay \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\eta _c(1S)\,{\rightarrow } \,{{{ p}}} \overline{{{{ p}}}}} $$\end{document}ηc(1S)→pp¯

The production of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _c (1S)$$\end{document}ηc(1S) state in proton-proton collisions is probed via its decay to the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p\overline{p}$$\end{document}pp¯ final state with the LHCb detector, in the rapidity range \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2.0 < y < 4.5$$\end{document}2.0 6.5 \mathrm{{\,GeV/}{ c}} $$\end{document}pT>6.5GeV/c. The cross-section for prompt production of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _c (1S)$$\end{document}ηc(1S) mesons relative to the prompt \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{ J}}/{\psi } $$\end{document}J/ψ cross-section is measured, for the first time, to be \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma _{\eta _c (1S)}/\sigma _{{{{ J}}/{\psi }}} = 1.74\, \pm \,0.29\, \pm \, 0.28\, \pm \,0.18 _{{\mathcal{B}}}$$\end{document}σηc(1S)/σJ/ψ=1.74±0.29±0.28±0.18B at a centre-of-mass energy \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\sqrt{s}} = 7 {~\mathrm{TeV}}$$\end{document}s=7TeV using data corresponding to an integrated luminosity of 0.7 fb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{-1}$$\end{document}-1, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma _{\eta _c (1S)}/\sigma _{{{{ J}}/{\psi }}} = 1.60 \pm 0.29 \pm 0.25 \pm 0.17 _{{\mathcal{B}}}$$\end{document}σηc(1S)/σJ/ψ=1.60±0.29±0.25±0.17B at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\sqrt{s}} = 8 {~\mathrm{TeV}}$$\end{document}s=8TeV using 2.0 fb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{-1}$$\end{document}-1. The uncertainties quoted are, in order, statistical, systematic, and that on the ratio of branching fractions of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _c (1S)$$\end{document}ηc(1S) and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{ J}}/{\psi } $$\end{document}J/ψ decays to the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p\overline{p}$$\end{document}pp¯ final state. In addition, the inclusive branching fraction of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${b} $$\end{document}b-hadron decays into \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _c (1S)$$\end{document}ηc(1S) mesons is measured, for the first time, to be \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathcal{B}}( b {\rightarrow } \eta _c X ) = (4.88\, \pm \,0.64\, \pm \,0.29\, \pm \, 0.67 _{{\mathcal{B}}}) \times 10^{-3}$$\end{document}B(b→ηcX)=(4.88±0.64±0.29±0.67B)×10-3, where the third uncertainty includes also the uncertainty on the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{ J}}/{\psi } $$\end{document}J/ψ inclusive branching fraction from \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${b} $$\end{document}b-hadron decays. The difference between the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{ J}}/{\psi } $$\end{document}J/ψ and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _c (1S)$$\end{document}ηc(1S) meson masses is determined to be \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$114.7 \pm 1.5 \pm 0.1 {\mathrm {\,MeV\!/}c^2} $$\end{document}114.7±1.5±0.1MeV/c2.


Introduction
High centre-of-mass energies available in proton-proton collisions at the LHC allow models describing charmonium production to be tested. We distinguish promptly produced charmonia from those originating from b-hadron decays. Promptly produced charmonia include charmonia directly produced in parton interactions and those originating from the decays of heavier quarkonium states, which are in turn produced in parton interactions. While measurements of J/ψ and ψ(2S) meson production rates at the LHC [1][2][3][4][5][6] are successfully described by next-to-leading order (NLO) calculations in non-relativistic quantum chromodynamics (QCD) [7], the observation of small or no polarization in J/ψ meson prompt production [2] remains unexplained within e-mail: sergey.barsuk@cern.ch the available theoretical framework [8]. The investigation of the lowest state, the η c (1S) meson, can provide important additional information on the long-distance matrix elements [9,10]. In particular, the heavy-quark spin-symmetry relation between the η c (1S) and J/ψ matrix elements can be tested, with the NLO calculations predicting a different dependence of the production rates on charmonium transverse momentum, p T , for spin singlet (η c (1S)) and triplet (J/ψ, χ cJ ) states [11][12][13]. Thus, a measurement of the p T dependence of the η c (1S) production rate, in particular in the low p T region, can have important implications. Recent LHCb results on prompt production of χ c states [14] provide information on the production of the P-wave states χ c0 and χ c2 at low p T , using the well-understood χ c1 production as a reference. A measurement of the cross-section of prompt η c (1S) production may allow an important comparison with the χ c0 results and yields indirect information on the production of heavier states.
At LHC energies, all b-hadron species are produced, including weakly decaying B − , B 0 , B 0 s , B − c mesons, bbaryons, and their charge-conjugate states. A previous study of inclusive η c (1S) meson production in b-hadron decays by the CLEO experiment, based on a sample of B − and B 0 mesons, placed an upper limit on the combined inclusive branching fraction of B − and B 0 meson decays into final states containing an η c (1S) meson of B(B − , B 0 →η c (1S)X ) < 9 × 10 −3 at 90 % confidence level [15]. Exclusive analyses of η c (1S) and J/ψ meson production in b-hadron decays using the B→K (pp) decay mode have been performed by the BaBar experiment [16], by the Belle experiment [17] and recently by the LHCb experiment [18].
In the present paper we report the first measurement of the cross-section for the prompt production of η c (1S) mesons in pp collisions at √ s = 7 TeV and √ s = 8 TeV centreof-mass energies, as well as the b-hadron inclusive branching fraction into η c (1S) final states. This paper extends the scope of previous charmonium production studies reported by LHCb, which were restricted to the use of J/ψ or ψ(2S) decays to dimuon final states [1,2,14,19]. In order to explore states that do not have J PC = 1 −− quantum numbers, while avoiding reconstruction of radiative decays with low-energy photons, the authors of Ref. [20] suggested to investigate hadronic final states. In the present analysis, we reconstruct η c (1S) mesons decaying into the p p final state. All well-established charmonium states decay to p p final states [20,21]. With its powerful charged-hadron identification and high charmonium production rate, the LHCb experiment is well positioned for these studies. The measurements are performed relative to the topologically and kinematically similar J/ψ→pp channel, which allows partial cancellation of systematic uncertainties in the ratio. This is the first such inclusive analysis using decays to hadronic final states performed at a hadron collider.
In addition, a departure in excess of two standard deviations between the recent BES III results [22,23] and earlier measurements [21] motivates the determination of the difference between J/ψ and η c (1S) meson masses using a different technique and final state. In the present analysis, the low-background sample of charmonia produced in b-hadron decays is used to determine M J/ψ, η c (1S) and the η c (1S) natural width, η c (1S) . In Sect. 2 we present the LHCb detector and data sample used for the analysis. Section 3 describes the analysis details, while the systematic uncertainties are discussed in Sect. 4. The results are given in Sect. 5 and summarized in Sect. 6.

LHCb detector and data sample
The LHCb detector [24] is a single-arm forward spectrometer covering the pseudorapidity range 2 < η < 5, designed for the study of particles containing b or c quarks. The detector includes a high-precision tracking system consisting of a silicon-strip vertex detector surrounding the pp interaction region, a large-area silicon-strip detector located upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip detectors and straw drift tubes placed downstream of the magnet. The combined tracking system provides a momentum measurement with a relative uncertainty that varies from 0.4 % at low momentum to 0.6 % at 100 GeV/c, and an impact parameter measurement with a resolution of 20 μm for charged particles with large transverse momentum. Different types of charged hadrons are distinguished using information from two ring-imaging Cherenkov detectors. Photon, electron, and hadron candidates are identified by a system consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter, and a hadronic calorimeter. Muons are identified by a system composed of alternating layers of iron and multiwire proportional chambers.
The trigger consists of a hardware stage, based on information from the calorimeter and muon systems, followed by a software stage, which applies a full event reconstruction.
Events enriched in signal decays are selected by the hardware trigger, based on the presence of a single highenergy deposit in the calorimeter. The subsequent software trigger specifically rejects high-multiplicity events and selects events with two oppositely charged particles having good track-fit quality and transverse momentum larger than 1.9 GeV/c. Proton and antiproton candidates are identified using the information from Cherenkov and tracking detectors [25]. Selected p and p candidates are required to form a good quality vertex. In order to further suppress the dominant background from accidental combinations of random tracks (combinatorial background), charmonium candidates are required to have high transverse momentum, p T > 6.5 GeV/c.
The present analysis uses pp collision data recorded by the LHCb experiment at √ s = 7 TeV, corresponding to an integrated luminosity of 0.7 fb −1 , and at √ s = 8 TeV, corresponding to an integrated luminosity of 2.0 fb −1 .
Simulated samples of η c (1S) and J/ψ mesons decaying to the pp final state, and J/ψ decaying to the ppπ 0 final state, are used to estimate efficiency ratios, the contribution from the decay J/ψ→ppπ 0 , and to evaluate systematic uncertainties. In the simulation, pp collisions are generated using Pythia [26] with a specific LHCb configuration [27]. Decays of hadronic particles are described by EvtGen [28], in which final-state radiation is generated using Photos [29]. The interaction of the generated particles with the detector and its response are implemented using the Geant4 toolkit [30,31] as described in Ref. [32].

Signal selection and data analysis
The signal selection is largely performed at the trigger level. The offline analysis, in addition, requires the transverse momentum of p and p to be p T > 2.0 GeV/c, and restricts charmonium candidates to the rapidity range 2.0 < y < 4.5.
Discrimination between promptly produced charmonium candidates and those from b-hadron decays is achieved using the pseudo-decay time t z = z · M/ p z , where z is the distance along the beam axis between the corresponding pp collision vertex (primary vertex) and the candidate decay vertex, M is the candidate mass, and p z is the longitudinal component of its momentum. Candidates with t z < 80 fs are classified as prompt, while those with t z > 80 fs are classified as having originated from b-hadron decays. For charmonium candidates from b-hadron decays, a significant displacement of the proton tracks with respect to the primary vertex is also required.
The selected samples of prompt charmonium candidates and charmonia from b-hadron decays have some candidates wrongly classified (cross-feed). The cross-feed probability is estimated using simulated samples and is scaled using the observed signal candidate yields in data. The crossfeed component is subtracted to obtain the ratio of produced η c (1S) and J/ψ mesons decaying into the p p final state. Corrections range from 2 % to 3 % for the ratio of promptly produced η c (1S) and J/ψ mesons, and from 8 % to 10 % for the ratio of charmonia produced in b-hadron decays.
The ratios of signal yields are expressed in terms of ratios of cross-sections multiplied by the decay branching fractions as where N P and N b are the numbers of charmonia from prompt production and b-hadron decays, respectively. The simulation describes the kinematic-related differences between the η c (1S) and J/ψ decay modes reasonably well and predicts that the relative efficiencies for selecting and reconstructing η c (1S) and J/ψ mesons differ by less than 0.5%. Equal efficiencies are assumed for the η c (1S) and J/ψ meson reconstruction and selection criteria. The efficiency for selecting and reconstructing prompt J/ψ mesons is corrected for polarization effects, as a function of rapidity and p T , according to Ref. [2]. The numbers of reconstructed η c (1S) and J/ψ candidates are extracted from an extended maximum likelihood fit to the unbinned p p invariant mass distribution. The J/ψ peak position M J/ψ and the mass difference M J/ψ,η c (1S) are fitted in the sample of charmonia from b-hadron decays, where the signal is more prominent because of the reduced background level due to charmonium decay-vertex displacement requirements. The results are then used to apply Gaussian constraints in the fit to the p p invariant mass spectrum in the prompt production analysis, where the signalto-background ratio is smaller, due to large combinatorial backgrounds.
The signal shape is defined by the detector response, combined with the natural width in the case of the η c (1S) resonance. The detector response is described using two Gaussian functions with a common mean value. In the description of each resonance, the ratio of narrow to wide Gaussian widths, , the fraction of the narrow Gaussian component, and the ratio of the η c (1S) and J/ψ narrow Gaussian widths, σ a η c (1S) /σ a J/ψ , are fixed in the fit to the values observed in simulation. The only resolution parameter left free in the fit to the low-background sample from b-hadron decays, σ a J/ψ , is fixed to its central value in the fit to the prompt sample. The natural width η c (1S) of the η c (1S) resonance is also extracted from the fit to the bhadron decays sample, and is fixed to that value in the prompt production analysis. Gaussian constraints on the J/ψ meson mass and the M J/ψ, η c (1S) mass difference from the fit to the b-hadron decays sample are applied in the prompt production analysis. The fit with free mass values gives consistent results.
The combinatorial background is parametrized by an exponential function in the fit of the sample from b-hadron decays, and by a third-order polynomial in the fit to the prompt sample.
Combinations of p p from the decay J/ψ→ppπ 0 potentially affect the region close to the η c (1S) signal; hence, this contribution is specifically included in the background description. It produces a non-peaking contribution, and its mass distribution is described by a square-root shape to account for the phase space available to the p p system from the J/ψ→ppπ 0 decay, convolved with two Gaussian functions to account for the detector mass resolution. In the fit to the p p invariant mass spectrum, the normalization of this contribution is fixed using the number of candidates found in the J/ψ signal peak and the ratios of branching fractions and efficiencies for the J/ψ→ppπ 0 and J/ψ→pp decay modes.
The p p invariant mass spectra for charmonium candidates from b-hadron decays in the 7 TeV and 8 TeV data are observed to be consistent. The two data samples are therefore combined and the resulting spectrum is shown in Fig. 1 with the fit overlaid.
The J/ψ meson signal is modelled using a double-Gaussian function. The η c (1S) signal is modelled using a relativistic Breit-Wigner function convolved with a double-Gaussian function. The background contribution from the J/ψ→ppπ 0 decay with an unreconstructed pion, is small. The fit yields 2020±230 η c (1S) signal decays and 6110±116 J/ψ signal decays.
The results of the fit to the p p invariant mass spectrum of the prompt sample are shown in Fig. 2a and b for data collected at √ s = 7 TeV and √ s = 8 TeV, respectively. The fits yield 13 370 ± 2260 η c (1S) and 11 052 ± 1004 J/ψ signal decays for the data taken at √ s = 7 TeV, and 22 416± 4072 η c (1S) and 20 217 ± 1403 J/ψ signal decays for the √ s = 8 TeV data. In order to assess the quality of these unbinned fits to the invariant pp mass spectra, the chisquare per degree of freedom was calculated for the binning schemes shown in Figs. 1, and 2a, b. The values are 1.3, 1.7 and 1.8, respectively.
From the observed η c (1S) and J/ψ yields, and taking into account cross-feed between the samples, the yield ratios are obtained as for the prompt production and charmonium production in bhadron decays. Only statistical uncertainties are given in the above ratios.

Systematic uncertainties
We consider systematic uncertainties due to limited knowledge of the detector mass resolution, the J/ψ polarization, the η c (1S) natural width, possible differences of the prompt charmonium production spectra in data and simulation, cross-feed between the prompt charmonium sample and the charmonium sample from b-hadron decays, background description and feed-down from J/ψ→ppπ 0 decays. Uncertainties due to limited knowledge of the detector mass resolution are estimated by assigning the same σ a value to the η c (1S) and J/ψ signal description for the bhadron sample, and by varying the σ a parameters in the prompt production analysis within their uncertainties. Uncertainties associated with the J/ψ polarization in the prompt production reflect those of the polarization measurement in Ref. [2]. We evaluate a potential contribution from J/ψ polarization in b-hadron decays using a BaBar study [32] of the J/ψ polarization in inclusive decays of B mesons. Simulations are used to estimate the effective polarization parameter for the LHCb kinematic region where the b-hadrons have a high boost and the longitudinal polarization is significantly reduced. A conservative value for the polarization parameter of −0.2 is used to estimate the corresponding systematic uncertainty. In order to estimate the systematic uncertainty associated with the η c (1S) natural width, which enters the results for the prompt production analysis, the world average η c (1S) value of 32.0 MeV from Ref. [21] is used. Possible differences of the prompt charmonium production spectra in data and simulation are estimated by correcting the efficiency derived from simulation according to the observed p T distribution. The uncertainty related to the cross-feed is estimated by varying the signal yields in each sample according to their uncertainties. Uncertainties associated with the background description are estimated by using an alternative parametrization and varying the fit range. The uncertainty due to the contribution from the J/ψ→ppπ 0 decay is dominated by the modelling of the pp invariant mass shape, and is estimated by using an alternative parametrization, which is linear instead of the square root. Possible systematic effect related to separation between prompt and b-decays samples, was checked by varying the t z discriminant value from 80 to 120 fs. The results are found to be stable under variation of the value of the t z discriminant, and no related systematic uncertainty is assigned. Table 1 lists the systematic uncertainties for the production yield ratio. The total systematic uncertainty is estimated as the quadratic sum of the uncertainties from the sources listed in Table 1 and, in the case of the prompt production measurement, is dominated by the uncertainty associated with the η c (1S) natural width. For the measurement with b-hadron decays the uncertainties associated with the background model, the J/ψ polarization and the cross-feed provide significant contributions.

Results
The yield ratio for charmonium production in b-hadron decays is obtained as .039 ± 0.015. In all quoted results, the first uncertainty refers to the statistical contribution and the second to the systematic contribution. By correcting the yield ratio with the ratio of branching fractions B(J/ψ→pp)/B(η c (1S)→pp) = 1.39 ± 0.15 [21], the ratio of the inclusive b-hadron branching fractions into η c (1S) and J/ψ final states for charmonium transverse momentum p T > 6.5 GeV/c is measured to be where the third uncertainty is due to that on the J/ψ→pp and η c (1S)→pp branching fractions [21]. Assuming that the p T > 6.5 GeV/c requirement does not bias the distribution of charmonium momentum in the b-hadron rest frame, and After correcting with the ratio of branching fractions B(J/ψ →pp)/B(η c (1S)→pp) [21], the relative η c (1S) to J/ψ prompt production rates in the kinematic regime 2.0 < y < 4.5 and p T > 6.5 GeV/c are found to be for the data sample collected at √ s = 7 TeV, and for the data sample collected at √ s = 8 TeV. The third contribution to the uncertainty is due to that on the J/ψ→pp and η c (1S)→pp branching fractions.
The absolute η c (1S) prompt cross-section is calculated using the J/ψ prompt cross-section measured in Refs. [2] and [3] and integrated in the kinematic range of the present analysis, 2.0 < y < 4.5 and p T > 6.5 GeV/c. The corresponding J/ψ prompt cross-sections were determined to be (σ J/ψ ) √ s=7 TeV = 296.9±1.8±16.9 nb for √ s = 7 TeV [2], and (σ J/ψ ) √ s=8 TeV = 371.4 ± 1.4 ± 27.1 nb for √ s = 8 TeV [3]. The J/ψ meson is assumed to be produced unpolarized. The prompt η c (1S) cross-sections in this kinematic region are determined to be (σ η c (1S) ) √ s=7 TeV = 0.52 ± 0.09 ± 0.08 ± 0.06 σ J/ψ , B μb, for √ s = 7 TeV, and (σ η c (1S) ) √ s=8 TeV = 0.59 ± 0.11 ± 0.09 ± 0.08 σ J/ψ , B μb, for √ s = 8 TeV. Uncertainties associated with the J/ψ→pp and η c (1S)→pp branching fractions, and with the J/ψ crosssection measurement, are combined into the last uncertainty component, dominated by the knowledge of the branching fractions. This is the first measurement of prompt η c (1S) production in pp collisions. The cross-section for the η c (1S) prompt production is in agreement with the colour-singlet leading order (LO) calculations, while the predicted crosssection exceeds the observed value by two orders of magnitude when the colour-octet LO contribution is taken into account [33]. However, the NLO contribution is expected to significantly modify the LO result [11]. Future measurements at the LHC design energy of √ s = 14 TeV may allow a study of the energy dependence of the η c (1S) prompt production.
The η c (1S) differential cross-section as a function of p T is obtained by fitting the p p invariant mass spectrum in three or four bins of p T . The same procedure as used to extract the η c (1S) cross-section is followed. The J/ψ p T spectrum measured in Refs. [1][2][3] is used to obtain the η c (1S) p T spectrum for both prompt production and inclusive η c (1S) production in b-hadron decays (Fig. 3). The p T dependence of the η c (1S) production rate exhibits similar behaviour to the J/ψ meson rate in the kinematic region studied.
The performance of the LHCb tracking system and the use of a final state common to J/ψ and η c (1S) decays allows a precise measurement of the mass difference between the two mesons. In order to measure the η c (1S) mass relative to the well-reconstructed and well-known J/ψ mass, a momentum scale calibration [34] is applied on data, and validated with the J/ψ mass measurement. The M J/ψ and M J/ψ, η c (1S) values are extracted from the fit to the p p invariant mass in the low-background sample of charmonium candidates produced in b-hadron decays (Fig. 1). The J/ψ mass measurement, M J/ψ = 3096.66 ± 0.19 ± 0.02 MeV/c 2 , agrees well with the average from Ref. [21]. The mass difference is measured to be The systematic uncertainty is dominated by the parametrization of the J/ψ→ppπ 0 contribution. The mass difference agrees with the average from Ref. [21]. In addition, the η c (1S) natural width is obtained from the fit to the p p invariant mass (Fig. 1), η c (1S) = 25.8 ± 5.2 ± 1.9 MeV. The systematic uncertainty is dominated by knowledge of the detector mass resolution. The value of η c (1S) obtained is in good agreement with the average from Ref. [21], but it is less precise than previous measurements.

Summary
In summary, η c (1S) production is studied using pp collision data corresponding to integrated luminosities of 0.7 fb −1 and 2.0 fb −1 , collected at centre-of-mass energies √ s = 7 TeV and √ s = 8 TeV, respectively. The inclusive branching fraction of b-hadron decays into η c (1S) mesons with p T > 6.5 GeV/c, relative to the corresponding fraction into J/ψ mesons, is measured, for the first time, to be B(b→η c (1S)X )/B(b→J/ψ X ) = 0.421 ± 0.055 ± 0.025 ± 0.045 B .
The first uncertainty is statistical, the second is systematic, and the third is due to uncertainties in the branching fractions of η c (1S) and J/ψ meson decays to the p p final state. The inclusive branching fraction of b-hadrons into η c (1S) mesons is derived as where the third uncertainty component includes also the uncertainty on the inclusive branching fraction of b-hadrons into J/ψ mesons. The η c (1S) prompt production crosssection in the kinematic region 2.0 < y < 4.5 and p T > 6.5 GeV/c, relative to the corresponding J/ψ meson crosssection, is measured, for the first time, to be where the third uncertainty component is due to uncertainties in the branching fractions of η c (1S) and J/ψ meson decays to the p p final state. From these measurements, absolute η c (1S) prompt cross-sections are derived, yielding (σ η c (1S) ) √ s=7 TeV = 0.52 ± 0.09 ± 0.08 ± 0.06 σ J/ψ , B μb, (σ η c (1S) ) √ s=8 TeV = 0.59 ± 0.11 ± 0.09 ± 0.08 σ J/ψ , B μb, where the third uncertainty includes an additional contribution from the J/ψ meson cross-section. The above results assume that the J/ψ is unpolarized. The η c (1S) prompt cross-section is in agreement with the colour-singlet LO calculations, whereas the colour-octet LO contribution predicts a cross-section that exceeds the observed value by two orders of magnitude [33]. The p T dependences of the η c (1S) and J/ψ production rates exhibit similar behaviour in the kinematic region studied. The difference between the J/ψ and η c (1S) meson masses is also measured, yielding M J/ψ, η c (1S) = 114.7 ± 1.5 ± 0.1 MeV/c 2 , where the first uncertainty is statistical and the second is systematic. The result is consistent with the average from Ref. [21].