Measurement of relative branching fractions of B decays to $\psi(2S)$ and $J/\psi$ mesons

The relative rates of B-meson decays into $J/\psi$ and $\psi(2S)$ mesons are measured for the three decay modes in pp collisions recorded with the LHCb detector. The ratios of branching fractions ($\mathcal{B}$) are measured to be $\frac{\mathcal{B}(B^+ \to \psi(2S) K^+)}{\mathcal{B}(B^+ \to J/\psi K^+)} = 0.594 \pm 0.006 (stat) \pm 0.016 (syst) \pm 0.015 (R_{\psi})$, $\frac{\mathcal{B}(B^0 \to \psi(2S) K^{*0})}{\mathcal{B}(B^0 \to J/\psi K^{*0})} = 0.476 \pm 0.014 (stat) \pm 0.010 (syst) \pm 0.012\,(R_{\psi})$, $\frac{\mathcal{B}^{0}_{s}(B^0_s \to \psi(2S)\phi)}{\mathcal{B}(B^0_s \to J/\psi\phi)} = 0.489 \pm 0.026 (stat) \pm 0.021 (syst) \pm 0.012\,(R_{\psi})$ where the third uncertainty is from the ratio of the $\psi(2S)$ and $J/\psi$ branching fractions to $\mu\mu$.


arXiv:1205.0918v1 [hep-ex] 4 May 2012 1 Introduction
Decays of B mesons to two-body final states containing a charmonium resonance such as a J/ψ or ψ(2S) offer a powerful way of studying electroweak transitions.Such decays probe charmonium properties and play a role in the study of CP violation and mixing in the neutral B system [1].
The relative branching fractions of B + , B 0 and B 0 s mesons into J/ψ and ψ(2S) mesons have previously been studied by both the CDF and D0 collaborations [2][3][4].Since the current experimental results for the study of CP violation in B 0 s mixing using the B 0 s → J/ψ φ decay [5][6][7] are statistically limited, it is important to establish other channels where this analysis can be done.One such channel is the B 0 s → ψ(2S)φ decay.In this paper, measurements of the ratios of the branching fractions of B mesons decaying to ψ(2S)X and J/ψ X are reported, where B denotes a B + , B 0 or B 0 s meson (charge conjugate decays are implicitly included) and X denotes a K + , K * 0 or φ meson.The data were collected by the LHCb experiment in pp collisions at the centre-of-mass energy √ s = 7 TeV during 2011 and correspond to an integrated luminosity of 0.37 fb −1 .

Detector description
The LHCb detector [8] is a single-arm forward spectrometer covering the pseudorapidity range 2 < η < 5, designed for the study of b-and c-hadrons.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.The combined tracking system has a momentum resolution ∆p/p that varies from 0.4% at 5 GeV/c to 0.6% at 100 GeV/c, and an impact parameter resolution of 20 µm for tracks with high transverse momentum.Data were taken with both magnet polarities to reduce systematic effects due to detector asymmetries.Charged hadrons are identified using two ring-imaging Cherenkov (RICH) detectors.Photon, electron and hadron candidates are identified by a calorimeter system consisting of scintillating-pad and pre-shower detectors, and electromagnetic and hadronic calorimeters.Muons are identified by a muon 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 with a J/ψ → µ + µ − final state are triggered using two hardware trigger decisions: the single-muon decision, which requires one muon candidate with a transverse momentum p T larger than 1.5 GeV/c, and the di-muon decision, which requires two muon candidates with transverse momenta p T 1 and p T 2 satisfying the relation 3 GeV/c.The di-muon trigger decision in the software trigger requires muon pairs of opposite charge with p T > 500 MeV/c, forming a common vertex and with an invariant mass in excess of 2.9 GeV/c 2 .

Event selection
In this analysis, the decays B + → ψK + (B 0 → ψK * 0 , B 0 s → ψφ) are reconstructed, where ψ represents ψ(2S) or J/ψ , reconstructed in the ψ → µ + µ − decay modes.A K + (K * 0 , φ) candidate is added to the di-muon pair to form a B + (B 0 , B 0 s ) candidate.The starting point of the analysis is the reconstruction of either a J/ψ or ψ(2S) meson decaying into a di-muon pair.Candidates are formed from pairs of opposite sign tracks that both have a transverse momentum larger than 500 MeV/c.Good reconstruction quality is assured by requiring the χ 2 per degree of freedom of the track fit to satisfy χ 2 /ndf < 5.Both tracks must be identified as muons.This is achieved by requiring the muon identification variable, the difference in logarithm of the likelihood of the muon and hadron hypotheses [9] provided by the muon detection system, to satisfy ∆ log L µ−h > −5.The muons are required to form a common vertex of good quality (χ 2 vtx < 20).The resulting di-muon candidate is required to have decay length significance from its associated primary vertex greater than 5 and have an invariant mass between 3020 and 3135 MeV/c 2 in the case of a J/ψ candidate or between 3597 and 3730 MeV/c 2 for a ψ(2S) candidate.These correspond to [−5σ; 3σ] windows around the nominal mass.The asymmetric window allows for the QED radiative tail.
The selected J/ψ and ψ(2S) candidates are then combined with a K + , K * 0 or φ to create B meson candidates.Only the K * 0 → K + π − and φ → K + K − decay modes are considered.Pion-kaon separation is provided by the ring-imaging Cherenkov detectors.To identify kaons the difference in logarithm of the likelihood of the kaon and pion hypotheses [9] is required to satisfy ∆ log L K−π > −5.In the case of pions the difference in logarithm of the likelihood of the pion and kaon hypotheses [9] is required to satisfy ∆ log L π−K > −5.As in the case of muons, a cut is applied on the track χ 2 /ndf provided by the track fit at 5. The kaons and pions are required to have a transverse momentum larger than 250 MeV/c and to have an impact parameter significance with respect to any primary vertex larger than 2. In the B 0 channel, the mass of the kaon and pion system is required to be 842 < M K + π − < 942 MeV/c 2 and in the B 0 s channel the mass of the kaon pair is required to be 1010 In addition, we require the decay time of the B candidate (cτ ) to be larger than 100 µm to reduce the large combinatorial background from particles produced in the primary pp interaction.A global refit of the three-prong (four-prong) combination is performed with a primary vertex constraint and with the di-muon pair mass constrained to the nominal value [10] using the Decay Tree Fit (DTF) procedure [11].The reduced χ 2 of this fit (χ 2 DTF /ndf) is required to be less than 5, where the DTF algorithm takes into account the number of decay products to determine the number of degrees of freedom.The B + candidates, where a muon from the ψ(2S) → µ + µ − decay is reconstructed as both muon and kaon, are removed by requiring the angle between the same sign muon and kaon to be greater than 3 mrad.4 Measurement of N ψ(2S)X /N J/ψ X The mass distributions for selected candidates are shown in Fig. 1.The number of the B + → ψK + candidates is estimated by performing an unbinned maximum likelihood fit.The same procedure is used to determine the number of the B 0 → ψK + π − candidates in a 842 < M K + π − < 942 MeV/c 2 mass window and the number of the B 0 s → ψK + K − candidates in a 1010 < M K + K − < 1030 MeV/c 2 mass window.The number of signal candidates is determined by fitting a double-sided Crystal Ball function [12,13] for signal together with an exponential function to model the background.The tail parameters of the Crystal Ball function are fixed to values determined from simulation.
The B 0 mass distributions include the contributions from resonant decays (B 0 → ψK * 0 ), non-resonant decays (B 0 → ψK + π − ) and combinatorial background.The contributions from resonant and non-resonant modes are separated with the sPlot technique [14].The K + π − invariant mass is used as a discriminating variable to unfold the B 0 mass dis-Table 1: Summary of the signal yields for the six B modes considered and the ratios of the number of J/ψ and ψ(2S) decays: N total is the summed signal yield for resonant and nonresonant modes, N non−res is the signal yield for non-resonant modes only and N res ψX is the signal yield for resonant decays (through K * 0 or φ).The uncertainties are statistical only.

B decay modes
A fit is then performed to the unfolded B 0 distribution, which contains both non-resonant B 0 → ψK + π − decays and background, to determine the number of non-resonant decays.The final number of resonant decays is calculated by subtracting the number of non-resonant decays from the total number of decays.For the B 0 s modes the number of non-resonant decays are obtained by a similar procedure using the K + K − invariant mass as the discriminating variable.The signal yields and their ratios are summarized in Table 1.

Efficiencies and systematic uncertainties
The branching fraction ratio is calculated using where N res is the number of signal candidates and ε is the overall efficiency.The overall efficiency is the product of the geometrical acceptance of the detector, the combined reconstruction and selection efficiency, and the trigger efficiency.The efficiency ratio is estimated using simulation for all six modes.The simulation samples used are based on the Pythia 6.4 generator [15] configured with the parameters detailed in Ref. [16].Final state QED radiative corrections are included using the Photos package [17].The EvtGen [18] and Geant4 [19] packages are used to generate hadron decays and simulate interactions in the detector, respectively.The digitized output is passed through a full simulation of both the hardware and software trigger and then reconstructed in the same way as the data.
The overall efficiency ratio is 0.901 ± 0.016, 1.011 ± 0.014 and 0.994 ± 0.014 for the B + , the B 0 and the B 0 s channels respectively.Since the selection criteria for B → J/ψ X and B → ψ(2S)X decays are identical, the ratio of efficiencies is expected to be close to unity.The deviation of the overall efficiency ratio from unity in the case of the B + → ψK + decays is due to the difference between the p T spectra of muons for the J/ψ and ψ(2S) decays.For the B 0 and B 0 s channels this difference is small.It has been checked that the behaviour of the efficiencies of all selection criteria is consistent in the data and simulation.
Since the decay products in each of the pairs of channels considered have similar kinematics, most uncertainties cancel in the ratio.The different contributions to the systematic uncertainties affecting this analysis are discussed in the following and summarized in Table 2.
The dominant source of systematic uncertainty arises from the subtraction of the nonresonant components in the B 0 and the B 0 s decays.The non-resonant background is studied with two alternative methods.First, determining the number of B 0 (s) → ψK * 0 (φ) decays directly using the sPlot technique by unfolding and fitting the B 0 (s) mass distribution of candidates containing genuine K * 0 (φ) resonances.Second, using the B 0 (s) mass distribution as the discriminating variable to unfold the K + π − (K + K − ) mass distribution of genuine B 0 (s) candidates and fitting this distribution to determine the number of nonresonant decays.The corresponding uncertainties are found to be 1.5% in the B 0 channel and 3.4% in the B 0 s channel.The other important source of uncertainty arises from the estimation of the efficiencies due to the potential disagreement between data and simulation.This is studied by varying the selection criteria in data and simulation.The corresponding uncertainties are found to be 1.7% in the B + channel, 0.5% in the B 0 channel and 2.0% in the B 0 s channel.The observed difference in the efficiency ratios for the two magnet polarities is conservatively taken as an estimate of the systematic uncertainty.This is 1.4% in the B + channel, 0.6% in the B 0 channel and 0.7% in the B 0 s channel.The trigger is highly efficient in selecting B meson decays with two muons in the final state.For this analysis the di-muon pair is required to trigger the event.Differences in the trigger efficiency between data and simulation are studied in the data using events which were triggered independently on the di-muon pair [20].Based on these studies, an uncertainty of 1.1% is assigned.
A further uncertainty arises from the imperfect knowledge of the shape of the signal and background in the B meson mass distribution.To estimate this effect, a linear and a quadratic function are considered as alternative models for the background mass distribution.In addition, a double Gaussian shape and a sum of double-sided Crystal Ball and Gaussian shapes are used as alternative models for the signal shape.The maximum observed change in the ratio of yields in the ψ(2S) and J/ψ modes is taken as systematic uncertainty.
The central value of the relative efficiency is determined by assuming that the angular distribution of the B → ψ(2S)X decay is the same as that of the B → J/ψ X.The systematic uncertainty due to the unknown polarization of the ψ(2S) in the B meson decays is estimated as follows.The simulation samples were re-weighted to match the angular distributions found from the data and the relative efficiency was recalculated.The difference between the baseline analysis and the re-weighted simulation is taken as the systematic uncertainty, as shown in Table 2.
Finally, the uncertainty due to potential contribution from the Cabibbo-suppressed mode with a π misidentified as K is found to be 0.4% in the B + channel and negligible in the B 0 and B 0 s channels.The uncertainty due to the cross-feed between B 0 and B 0 s channels with a π misidentified as K (or a K misidentified as π) is negligible.
The results are combined using Eq.where the first uncertainty is statistical, the second is systematic and the third is the uncertainty on the R ψ value [10].

Table 2 :
Systematic uncertainties (in %) on the relative branching fractions.