Measurement of the branching fraction and $C\!P$ asymmetry in $B^{+}\rightarrow J/\psi \rho^{+}$ decays

The branching fraction and direct $C\!P$ asymmetry of the decay $B^{+}\rightarrow J/\psi \rho^{+}$ are measured using proton-proton collision data collected with the LHCb detector at centre-of-mass energies of 7 and 8 TeV, corresponding to a total integrated luminosity of 3$\mbox{fb}^{-1}$. The following results are obtained: \begin{align} \mathcal{B}(B^{+}\rightarrow J/\psi \rho^{+})&= (3.81 ^{+0.25}_{-0.24} \pm 0.35) \times 10^{-5}, \nonumber \\ \mathcal{A}^{C\!P} (B^{+}\rightarrow J/\psi \rho^{+})&= -0.045^{+0.056}_{-0.057} \pm 0.008, \nonumber \end{align} where the first uncertainties are statistical and the second systematic. Both measurements are the most precise to date.


Introduction
In the Standard Model of particle physics, the decay B + → J/ψ ρ + proceeds predominantly via a b → ccd transition involving tree and penguin amplitudes, 1 as shown in Fig. 1. Interference between these two amplitudes can lead to direct CP violation that is measured through an asymmetry defined as No precise prediction for A CP exists, though it is expected to have an absolute value 0.35 [1] assuming isospin symmetry between the B 0 → J/ψ ρ 0 and the B + → J/ψ ρ + decays. Measurements of A CP provide an estimate of the imaginary part of the penguinto-tree amplitude ratio for the b → ccd transition. Similarly to the B 0 → J/ψ ρ 0 decay [2], the CP asymmetry is expected to be enhanced in this decay compared to the decay B 0 s → J/ψ φ [3,4]. Therefore its value can be used to place constraints on penguin effects in measurements of the CP -violating phase φ s from the decay B 0 s → J/ψ φ, assuming approximate SU(3) flavour symmetry and neglecting exchange and annihilation diagrams. The branching fraction and the value of A CP for B + → J/ψ ρ + decays were measured previously by the BaBar collaboration to be (5.0±0.7±0.3)×10 −5 and −0.11±0.12±0.08, respectively [5].
In this paper, the branching fraction and the direct CP asymmetry of the decay B + → J/ψ ρ + are measured using proton-proton (pp) collision data collected with the LHCb detector at centre-of-mass energies of 7 TeV (in 2011) and 8 TeV (in 2012), corresponding to a total integrated luminosity of 3 fb −1 . The B + → J/ψ ρ + decay is analysed using the J/ψ → µ + µ − , ρ + → π + π 0 and π 0 → γγ decays. Its branching fraction is measured relative to that of the abundant decay B + → J/ψ K + , which has the same number of charged final-state particles and contains a J/ψ meson as the decay of interest.

Detector and simulation
The LHCb detector [6,7] 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

Selection
The same selection criteria are placed on the J/ψ candidates for both decays, B + → J/ψ ρ + and B + → J/ψ K + . Each J/ψ candidate is composed of two tracks compatible with being muons that form a good quality common vertex significantly displaced from any reconstructed PV in the event. The invariant mass of this two-track combination must be within ±100 MeV/c 2 of the known J/ψ mass [14].
For the B + → J/ψ ρ + decay, each ρ + candidate is formed from a charged and a neutral pion. The charged pion is required to have a χ 2 IP value that is inconsistent with originating from any PV in the event, where χ 2 IP is defined as the difference between the vertexfit χ 2 of a given PV reconstructed with and without the particle under consideration. The charged pion also needs to have a momentum larger than 3 GeV/c and a transverse momentum larger than 250 MeV/c. To reject background from B + → J/ψ K * + decays 2 , with K * + → K + π 0 , where a kaon is misidentified as a pion, a stringent criterion is placed on the pion-identification quality, which is mainly derived using information from the two RICH detectors. The neutral pion is reconstructed from two well-separated clusters in the electromagnetic calorimeter. It is required to have a p T larger than 800 MeV/c. The ρ + candidate must have a transverse momentum larger than 800 MeV/c and an invariant mass m π + γγ between 400 MeV/c 2 and 1100 MeV/c 2 .
For the B + → J/ψ K + decay, the K + candidate is selected among particles that fulfil the same requirements applied to the charged pion in the B + → J/ψ ρ + decay, have a transverse momentum larger than 800 MeV/c and are identified as a kaons.
The B + candidate is formed from the J/ψ candidate and the ρ + or K + candidate by requiring that the three reconstructed tracks form a good-quality vertex with significant displacement from the PV. In events with multiple PVs, that with which the B + candidate has the smallest χ 2 IP is chosen. In addition, the B + candidate is constrained to originate from the PV using a kinematic fit [21]. In the B + → J/ψ ρ + channel, the dimuon invariant mass is constrained to the known J/ψ mass, and when considering the invariant mass of the B + candidate, the π 0 candidate is additionally constrained to its known mass. This χ 2 IP value is required to be small enough to be consistent with the B + meson having originated from the PV.
Decays of b hadrons with an additional charged particle are rejected by ensuring that the quality of the B + decay vertex significantly degrades when the closest additional track is included in the vertex fit. Finally, the mass of the B + candidate, m J/ψ π + π 0 (m µ + µ − K + ), is required to be between 5100 MeV/c 2 and 5700 MeV/c 2 (5180 MeV/c 2 and 5400 MeV/c 2 ). For the B + → J/ψ K + channel, this selection results in a good signal purity, and no further selection criteria are needed.
Two vetoes are applied for the B + → J/ψ ρ + channel to reject B + → J/ψ K + and B + → J/ψ π + decays combined with a random π 0 : the mass of the combination of the J/ψ and the charged pion candidate, evaluated under both the pion and kaon mass hypotheses, must be outside a 50 MeV/c 2 window around the known B + mass.
In addition to the preselection discussed above, a multivariate selection is performed on the B + → J/ψ ρ + channel using an artificial neural network from the TMVA package [22,23], mainly to reduce combinatorial background. The 12 input variables are as follows: the flight distance and direction angle of the B + candidate, defined as the angle between its momentum and the vector connecting its primary and decay vertices; the transverse momentum of the B + and ρ + candidates; the maximum transverse momentum of the two muons, the χ 2 IP of the charged pion; the minimum χ 2 IP of the two muons; the quality of the B + candidate vertex; the change in the vertex quality when adding the closest track that is not part of the signal candidate; the quality of a kinematic fit of the full decay chain; the quality of the π 0 identification; and the p T asymmetry in a cone around the flight direction of the B + meson, defined as where i runs over all final state tracks of B + → J/ψ ρ + and j over all other tracks in a cone around the B + meson flight direction. The classifier is trained using background events from both a lower (4800-5000 MeV/c 2 ) and upper (5700-6000 MeV/c 2 ) sideband of the B + candidate mass and simulated signal decays, where the p T distribution of the B + meson and the number of tracks per event are weighted to match the corresponding distributions in a B 0 → J/ψ K * 0 data sample, where K * 0 → K + π − . The B 0 → J/ψ K * 0 candidates are obtained using the same preselection on the charged particles as for the B + → J/ψ ρ + decay, except that a looser particle-identification criterion is used for the pion. Moderate criteria are placed on the kaon p, p T and particle identification to obtain a good signal purity. For quantities depending on the kinematics of the π 0 meson, the channel B + → J/ψ K * + , with K * + → K + π 0 , is used to compare the distributions of the input variables between simulation and data. For all variables good agreement is found.
The neural network selection criterion is optimized by maximizing the figure of merit N sig / N sig + N bg , where N sig is the expected number of signal decays and N bg is the estimated background yield, both between 5200 MeV/c 2 and 5450 MeV/c 2 . The value of N sig is calculated using the ratio of the previously measured B + → J/ψ ρ + and B + → J/ψ K + branching fractions [14], the observed number of B + → J/ψ K + decays, and the efficiencies of the B + → J/ψ ρ + and B + → J/ψ K + channels. The value of N bg is obtained by extrapolating the shape of the background into the signal region. The optimized cut rejects 99.4% of background events in both, the 2011 and 2012 data samples, while retaining 49% (45%) of signal events in the 2011 (2012) data samples.
After applying the full selection for B + → J/ψ ρ + decays, about 9% of events have more than one candidate. Most of these events contain a genuine B + → J/ψ ρ + decay, along with a candidate comprised of the charged particles from the signal decay combined with a prompt π 0 meson. The latter constitutes a peaking background that is difficult to model, and therefore, events with multiple candidates are removed from the analysis. The efficiency of this rejection is evaluated on simulated samples. The sample of B + → J/ψ K + decays contains only 0.2% events with multiple candidates and no rejection is required.

Efficiencies
The branching ratio of the decay B + → J/ψ ρ + is calculated using The efficiencies ε B + →J/ψ K + and ε B + →J/ψ ρ + are composed of the geometrical acceptance, trigger, reconstruction, particle identification and selection efficiencies. In addition, there is an efficiency due to the removal of multiple candidates in the B + → J/ψ ρ + sample. The efficiency for the decay products of a B + → J/ψ ρ + or B + → J/ψ K + decay to be within the acceptance of the LHCb detector is taken from simulation. The efficiency to trigger on one or both muons from the J/ψ decay is also taken from simulation. When forming the ratio of branching fractions between the B + → J/ψ ρ + and B + → J/ψ K + decays, the ratio of trigger efficiencies is close to unity.
The charged-particle reconstruction efficiency is taken from simulation with kinematicsdependent correction factors applied that are determined using a tag-and-probe technique applied on a detached J/ψ → µ + µ − data sample [24]. The π 0 reconstruction efficiency is also obtained from simulation, with a p T -dependent correction factor obtained from the deviation between the known and observed ratio of branching fractions of B + → J/ψ K + and B + → J/ψ K * + decays, with K * + → K + π 0 followed by π 0 → γγ [25]. Data samples of B + → J/ψ K + and B + → J/ψ K * + decays, collected in 2011 and 2012 by the LHCb experiment, are fitted to obtain the number of observed decays. After correcting for the charged-particle selection efficiencies, the double ratio between the ratio of observed decays and the ratio of known branching fractions is the efficiency to reconstruct a π 0 meson. This result is compared to the π 0 efficiency obtained from simulated samples to determine the correction factors in intervals of the p T of the π 0 meson. Using this procedure, the uncertainty on the branching fraction of the B + → J/ψ K + decay cancels in the measurement of the B + → J/ψ ρ + branching fraction, and only the uncertainty on the branching fraction of the B + → J/ψ K * + decay contributes.
The charged-particle identification efficiency is evaluated using a tag-and-probe technique on dedicated calibration data samples with clean signatures [26]. Given the similar kinematics for the muons from B + → J/ψ ρ + and B + → J/ψ K + decays, the muon identification efficiency fully cancels when forming the ratio of branching fractions. The charged-pion identification efficiency is 72% for the 2011 data and 74% for the 2012 data. The efficiencies for identifying positively and negatively charged pions are compatible within their statistical uncertainties. The efficiency for identifying the kaon is above 95% for both data taking periods.
The remaining offline selection efficiencies are taken from simulation, where all kinematic distributions for the signal decay are found to be compatible with the corresponding distributions observed in data for B 0 → J/ψ K * 0 and B + → J/ψ K * + decays.

Invariant mass fits
The yield of B + → J/ψ K + decays is determined using an extended unbinned maximumlikelihood fit to the m µ + µ − K + distribution. The signal shape is described using the sum of two Crystal Ball functions [27] and a Gaussian function, where all three functions share the peak position value. The tail parameters of the Crystal Ball functions are fixed to the values obtained from simulation. An exponential function is used to model the combinatorial background. The fit yields a total of 362 739 ± 992 B + → J/ψ K + decays in 2011 and 816 197 ± 1545 in 2012. The fit is performed separately for the two magnet polarities; Figure 2 shows the fit to the data taken with down polarity.
The yield of B + → J/ψ ρ + decays is determined using a simultaneous two-dimensional extended unbinned maximum-likelihood fit to the m J/ψ π + π 0 and m π + γγ distributions in the 2011 and 2012 data sets. The signal shape in m J/ψ π + π 0 is modelled with the sum of two Crystal Ball functions with a shared peak position value. The values of the tail parameters are taken from simulation. The signal shape in m π + γγ is modelled with a relativistic Breit-Wigner function with Here, P J/ψ is the momentum of the J/ψ meson in the B + rest frame, q(m π + γγ ) is the pion momentum in the dipion rest frame, L eff is the relative angular momentum of the ρ + meson with respect to the J/ψ , q 0 = q(m ρ + ), Γ 0 is the nominal width, and R is a barrier factor radius. The tail parameters of the Crystal Ball function, along with m ρ + and Γ 0 of the Breit-Wigner, are fixed to values obtained from a fit to simulated B + → J/ψ ρ + decays, generated using the same functional form with input values m ρ + = 768.5 MeV/c 2 and Γ 0 = 151 MeV/c 2 . This strategy therefore folds in the effect of the m π + γγ resolution into the Breit-Wigner model, which results in a value of Γ 0 = 175 MeV/c 2 . The relative angular momentum L eff is fixed to 0. Non-resonant J/ψ π + π 0 decays are described with the same shape as the signal in m J/ψ π + π 0 , while in m π + γγ a three-body phase-space distribution multiplied by P 2 J/ψ is used, which is motivated by angular momentum conservation. Given the slowly varying shape in m π + γγ , no description of the detector resolution is needed. The combinatorial background is described by an exponential function in m J/ψ π + π 0 and a first-order polynomial in m π + γγ .
Two partially reconstructed backgrounds are included in the B + → J/ψ ρ + fit. The first is the decay B + → J/ψ K * + , with K * + → K 0 S π + and K 0 S → π 0 π 0 , where one π 0 is not reconstructed. The second is the decay B 0 s → J/ψ φ, with φ → π + π − π 0 , where one charged pion is not reconstructed. As the shapes of both contributions exhibit significant correlations between m J/ψ π + π 0 and m π + γγ , they are described by two-dimensional kernel density estimators [28] with adaptive kernels determined using simulation.
The final component is the background from B + → J/ψ K * + decays, with K * + → K + π 0 , where the kaon is misidentified as a pion. Despite the stringent particle-identification criterion applied, a small amount of these decays is present in the final sample. The product of two one-dimensional kernel density estimators is used to describe the shape in the two invariant mass distributions. This background yield is fixed relative to that of the B + → J/ψ K * + decay, with K * + → K 0 S π + , using the known ratio of the branching fractions [14].
The value of A CP is calculated using where A CP raw is the raw asymmetry, determined from a fit to the 2011 and 2012 data sets, split by the charge of the B meson, where all fit components are modelled as in the branching fraction measurement; and A prod is the production asymmetry of B + mesons in LHCb. For 2011 and 2012 they were measured to be (−0.41 ± 0.49 ± 0.11)% and (−0.53 ± 0.31 ± 0.10)%, respectively [29], with the first uncertainty being statistical and the second systematic.
For the background contribution from the B + → J/ψ K * + decay with K * + → K 0 S π + , the CP asymmetry is fixed to the known value of (−4.8 ± 3.3)% [14] after correcting for the production asymmetry. For the partially reconstructed background from B 0 s → J/ψ φ decays, no charge asymmetry is expected, as the final state is identical for both B 0 s and B 0 s mesons.

.1 Uncertainties on the branching fraction
The systematic uncertainties for the branching fraction measurement are summarized in Table 1. The trigger efficiencies are derived from simulation. The ratio of efficiencies for B + → J/ψ K + and B + → J/ψ ρ + decays has a small deviation from unity due to the slightly different kinematical distributions of the J/ψ mesons. To account for potential mismodelling of the impact of the K + and ρ + on the J/ψ trigger efficiency in simulation, the p T distributions of the B + and J/ψ mesons in B + → J/ψ K + decays are weighted to match those in B + → J/ψ ρ + decays, and the trigger efficiency is reevaluated. The resulting difference between this ratio of trigger efficiencies and unity is taken as a systematic uncertainty. For the charged-particle reconstruction efficiency, the correction factors between simulation and data are varied within their uncertainties, and the effect on the ratio of B + → J/ψ ρ + and B + → J/ψ K + decays is evaluated. The uncertainties consist of both statistical and systematic components, where the latter are due to the limited precision on the knowledge of the LHCb material budget. Only the contribution of the material budget uncertainty and the different interaction cross-section of pions and kaons with the material significantly contribute to the uncertainty.
The uncertainty on the π 0 reconstruction efficiency is dominated by the uncertainty on the branching fraction of the B + → J/ψ K * + decay, which is used in its calculation, and the number of B + → J/ψ K * + candidates in each kinematic bin. These values are added in quadrature to the systematic uncertainty of the method.
The systematic uncertainty for the charged-pion identification efficiency is evaluated by calculating the difference in efficiency between the nominal method, which is performed in intervals of p T , pseudorapidity, and the number of tracks in the event, and an alternative, unbinned method. In addition, a different quantity for the event multiplicity and a different scheme of interval boundaries are used. All three uncertainties are added in quadrature. A similar procedure is used for the kaon and muon identification efficiencies, resulting in smaller systematic uncertainties.
For B + → J/ψ K + decays, the uncertainty on the selection efficiency is dominated by the limited size of the simulated data set. For the B + → J/ψ ρ + selection, the largest component of this uncertainty arises from potential discrepancies in the π 0 identification variable between simulation and data. The decay B + → J/ψ K * + is used to weight the distribution of this variable, and the change in the multivariate classifier efficiency with respect to the baseline value is taken as a systematic uncertainty. A smaller contribution comes from the limited size of the simulated data set.
The uncertainty on the procedure to remove multiple-candidate events in the B + → J/ψ ρ + channel is determined on simulation by calculating the efficiency to select events that contain only one candidate, where the simulated events are weighted so that the event occupancy matches that observed in data. In addition, the invariant mass distributions in data are fitted for all events and events that only contain one candidate, where the signal shape is unchanged for both fits, i.e. the peaking background arising when the charged particles from a genuine signal decay are combined with a prompt π 0 meson is ignored. The difference in efficiency between data and simulation is taken as a systematic uncertainty.
The two-dimensional kernel density estimators use adaptive kernel widths. The effect of the kernel width is tested by setting it to a constant value, either higher or lower than the default value, and taking the resulting difference in the branching fraction as a systematic uncertainty. The tail parameters of the Crystal Ball function in B + → J/ψ ρ + decays are varied by ±20% with respect to the baseline values to account for the uncertainties in the fit to simulation. The value of m ρ + of the Breit-Wigner function is left free to vary, instead of being fixed to its baseline value. Furthermore, to take a possible difference in the experimental resolution of the ρ + mass between data and simulation into account, the value of Γ 0 of the Breit-Wigner shape is altered by the difference in the width of the B + mass peak between data and simulation in the decay B + → J/ψ K * + . All differences with respect to the branching fraction of the baseline fit are added in quadrature and added to the overall systematic uncertainty. To take possible correlations in the signal shape between m J/ψ π + π 0 and m π + γγ into account, the baseline model is replaced by a two-dimensional kernel density estimator whose shape is obtained from simulation. The resulting difference in the branching fraction with respect to the baseline model is taken as a systematic uncertainty.
The polarization of the decay products of the B + meson in the B + → J/ψ ρ + decay is unknown. Using simulated decays where L eff is set to 1 or 2, the branching fraction is recalculated and the difference in the observed branching fraction with respect to the baseline result is taken as a systematic uncertainty. As a consistency test, simulated decays with L eff = 0 are studied. A small bias with respect to the baseline branching fraction is observed, which is corrected for and taken as a systematic uncertainty.
To test the effect of different polarization amplitudes of B + → J/ψ ρ + in simulation and data on the efficiency to reconstruct the decay, the values of the helicity amplitudes from B 0 → J/ψ ρ 0 [30] are used in the simulation of B + → J/ψ ρ + decays and varied within their uncertainties. The largest deviation with respect to the baseline value is taken as a systematic uncertainty.
In order to estimate the systematic uncertainty from the chosen fit range, the fit to determine the branching fraction is repeated 1000 times with random intervals in m J/ψ π + π 0 and m π + γγ larger or smaller than those of the baseline fit. The width of the resulting distribution of the measured branching fractions of B + → J/ψ ρ + decays for each interval is taken as a systematic uncertainty. A small bias with respect to the baseline result was observed. The bias is corrected for and also taken as a systematic uncertainty.
As an alternative modelling for the nonresonant contribution in m π + γγ , the shape is modelled with the form (P J/ψ /m B ) 2 [31], with m B the known mass of the B + meson [14]. The difference of the branching fraction with respect to the baseline fit is taken as a systematic uncertainty.
A contribution from interference between the nonresonant B + → J/ψ π + π 0 and signal B + → J/ψ ρ + decays could arise due to an asymmetric efficiency of the angular distribution of the π + π 0 system. To assess the effect of not taking the interference term in the branching fraction fit into account, several samples of one million simulated decays are generated including the interference term, where for each sample a different fixed phase for the nonresonant contribution is chosen. The shape of the angular acceptance is taken from the full simulation of B + → J/ψ ρ + decays. Each sample is fitted with the baseline description, which does not include the interference term. The largest relative difference in the signal yield with respect to the generated value is taken as a systematic uncertainty. To investigate possible exotic resonance contributions, the invariant masses of J/ψ π + and J/ψ π 0 are inspected and no excess compared to the expectation is found. Given the nature of the systematic uncertainties on the π 0 and charged-particle reconstruction efficiencies, the selection efficiency and the removal of multiple candidates, their correlation in 2011 and 2012 is set to 1. All other uncertainties either result from a common fit to the combined data sets of 2011 and 2012 or are treated as uncorrelated.

Uncertainties on the CP asymmetry
Most systematic uncertainties cancel when calculating the A CP ratio. The remaining contributions are listed in Table 2. The largest contributions come from the uncertainty on the knowledge of the direct CP asymmetry of the B + → J/ψ K * + decay for the partially reconstructed background, and the limited knowledge of the production asymmetry for B + mesons in the 2011 and 2012 data sets. The signal model is again replaced by a two-dimensional kernel density estimator to take possible correlations into account, taking the difference in the CP -asymmetry results as a systematic uncertainty.
To evaluate a possible bias on the asymmetry result, a permutation test is performed where the data set is split 1000 times randomly, instead of by the charge of the B + meson, and the asymmetry is evaluated. As an additional check, simulated decays are added to each of these randomly split samples, corresponding to ±5% and ±10% asymmetry, to assess the robustness of the asymmetry fit for a non-zero A CP value. No bias in the resulting distributions is observed. The total systematic uncertainty for A CP is formed by
Both results are the most precise to date and are consistent with previous measurements. Furthermore, the measured value of A CP is consistent with the corresponding measurement using B 0 → J/ψ ρ 0 decays, as expected from isospin symmetry [2]. LHCb collaboration