Search for the decays B0→e+e-γ and B0→μ+μ-γ

We present results of a search for the decays B 0 ! ‘ (cid:1) ‘ (cid:2) (cid:1) ( ‘ (cid:3) e or (cid:2) ). The search is performed using 320 (cid:4) 10 6 B (cid:1) B pairs collected at the (cid:2) (cid:5) 4 S (cid:6) resonance with the BABAR detector at the PEP-II B Factory at SLAC. We ﬁnd no signiﬁcant signal and set the following branching fraction upper limits at the 90% conﬁdence level: B (cid:5) B 0 ! e (cid:1) e (cid:2) (cid:1) (cid:6) 1 : 2 (cid:4) 10 (cid:2) 7 and B (cid:5) B 0 ! (cid:2) (cid:1) (cid:2) (cid:2) (cid:1) (cid:6) < : 6 (cid:4) 10 .

We present results of a search for the decays B 0 ! ' ' ÿ (' e or ). The search is performed using 320 10 6 B B pairs collected at the 4S resonance with the BABAR detector at the PEP-II B Factory at SLAC. We find no significant signal and set the following branching fraction upper limits at the 90% confidence level: BB 0 ! e e ÿ < 1:2 10 ÿ7 and BB 0 ! ÿ < 1: 6 10 ÿ7 . DOI: 10.1103/PhysRevD.77.011104 PACS numbers: 13.20.He, 14.40.Nd Radiative leptonic decays of neutral B mesons, B 0 ! ' ' ÿ with ' e or [1], are flavor-changing neutralcurrent transitions that are forbidden at the tree level in the standard model (SM). In the SM, such processes are described by penguin and box diagrams to leading order, as shown in Fig. 1. The largest contributions arise from processes in which a photon is emitted from one of the initial quarks, thus avoiding the helicity suppression of the purely leptonic decays B 0 ! ' ' ÿ . A search for the processes B 0 ! ' ' ÿ has been performed by the BABAR collaboration and others [2], but there is no previous search for the B 0 ! ' ' ÿ decays. The expected branching fraction for B 0 ! ' ' ÿ processes is of order 10 ÿ10 [3]. A measured branching fraction bigger than this would be evidence for new physics.
The analysis described in this paper uses a sample of 320 10 6 B B pairs recorded with the BABAR detector at the PEP-II asymmetric energy e e ÿ storage rings. This corresponds to an integrated luminosity of 292 fb ÿ1 collected at the 4S resonance.
A detailed description of the BABAR detector can be found elsewhere [4]. Charged-particle trajectories are measured by a five-layer silicon vertex tracker and a 40layer drift chamber operating in a 1.5 T magnetic field. A detector of internally reflected Cherenkov light is used for charged hadron identification. Surrounding this is an electromagnetic calorimeter (EMC) consisting of 6580 CsI(Tl) crystals, and the instrumented flux return for the solenoid, which consists of layers of steel interspersed with resistive plate chambers or limited streamer tubes.
A full BABAR Monte Carlo (MC) simulation using GEANT4 [5] is used to evaluate signal efficiencies and to identify and study background sources. The signal MC sample is based on a calculation where the B 0 ! ' ' ÿ transitions depend on three Wilson coefficients C 7 , C 9 , and C 10 at leading order [6].
We reconstruct the B 0 signal candidates by combining two oppositely charged leptons and a photon. The B 0 vertex is fitted using a Kalman Filter method [7]. The leptons are required to originate from a common vertex, and the B 0 candidate is required to be consistent with coming from the beam interaction point.
To minimize the number of misidentified particles, the leptons are required to satisfy stringent particle identification criteria [8]. For the electron candidates, bremsstrahlung photons are detected as energy deposits (clusters) in the EMC close to the intersection of their tracks with the EMC and the electron momentum is corrected for this energy loss. The corrected electron momentum is used in the calculation of the invariant mass of the reconstructed B mesons. For photon clusters, the transverse shower shape is required to be consistent with an electromagnetic shower. Leptons and photons are required to reside fully in the geometric acceptance of the detector.
Since the signal event contains two neutral B mesons and no additional particles, the total energy of each B meson in the center-of-mass (CM) frame must be equal to half of the total beam energy in the CM where E beam is the beam energy in the CM frame, p i and m i are the momenta in the CM frame, and the masses of the daughter particles i (i ' ; ' ÿ ; ), respectively. E beam is used instead of the measured B meson energy in the CM frame because E beam is more precisely known. For correctly reconstructed B 0 mesons, the m ES distribution has a maximum at the B 0 mass with a standard deviation of about 3 MeV=c 2 and the E distribution has a maximum near zero with a standard deviation of about 30 MeV.
The B 0 ! ' ' ÿ candidates are selected by requiring ÿ0:5 E 0:5 GeV and 5:0 m ES 5:3 GeV=c 2 . These ranges include both background-and signaldominated regions. As shown in Fig. 2, five backgrounddominated regions (sideband areas) are used for the background estimation. To avoid experimenter's bias, the events in the signal-dominated region (signal box) and in the shaded area covering the signal box are not included in the analysis until the final selection criteria have been optimized and the background estimation has been finalized. The size of the signal box is chosen to contain more than 80% of the signal: ÿ0:146ÿ0:112 E 0:082 GeV for the e e ÿ ( ÿ ) mode, and 5:270 m ES 5:289 GeV=c 2 for both modes. The asymmetry in E is due to the unrecovered energy losses of the charged leptons and the photon in the detector.
The dominant backgrounds are: (1) unmodeled higherorder QED and hadron production by two-photon processes for the e e ÿ mode; (2) B decays where the photon comes from a 0 decay, or one or both leptons are from a J= or 2S decay; and (3) continuum background from e e ÿ ! f f (where f u; d; s; c, or ) processes at the parton level.
To suppress higher-order QED processes and hadron production by two-photon annihilation, we tighten the fiducial requirements for electrons and photons and introduce additional selection criteria for the e e ÿ candidates: we require the cosine of the polar angle of e ÿ () to be between ÿ0:74 ( ÿ 0:62) and 0.81 (0.80) in the laboratory frame, the energy of the photon to be 0:3 GeV also in the laboratory frame, the number of charged tracks (EMC clusters) in the event to be 510, and the ratio of the second-to-zeroth order Fox-Wolfram moments (R 2 ) [9], which is calculated with the charged tracks and neutral clusters in the rest of the event (ROE), to be 0:7.
To reduce the number of events where the photon is from a 0 decay, we veto photon candidates that can be combined with any other photon in the event to form a 0 candidate with a mass within 3 standard deviations ( 20 MeV=c 2 ) of the nominal 0 mass.
We veto lepton candidates that form a suitable J= or 2S using the invariant mass of the lepton pair and E, as described in Ref. [10].
To suppress the continuum background, we require R 2 , calculated from all charged tracks and neutral clusters, to be less than 0.35, and the absolute value of the cosine of the angle between the thrust axis of the B 0 candidate and that of the ROE to be less than 0.8. These variables are used in a neural network combined with the following variables: (1) the absolute value of the cosine of the angle between the B 0 direction and the beam axis, (2) the absolute value of the cosine of the angle between the thrust axis of the B 0 candidate's decay products and the beam axis, (3) the ratio of second order to zeroth order Legendre moments of all charged tracks and neutral clusters, and (4) the invariant mass of the dileptons. The neural network rejects 20(36)% of the background while keeping 95(89)% of the signal, for the e e ÿ ( ÿ ) mode. All the selection criteria are optimized with MC samples to discriminate signal from background.
After all requirements are applied, there are on average 1.01(1.07) candidates per event for the e e ÿ ( ÿ ) mode. In events with multiple candidates, the one with the highest probability for the vertex fit is retained. The signal efficiency is 7.4(5.2)% for the e e ÿ ( ÿ ) mode. The e e ÿ mode has higher efficiency because electrons have higher detection efficiency than muons.
To assess possible background contributions that peak in the signal box, we examined 32 exclusive hadronic and semileptonic B decays using MC, including events where both B mesons decay semileptonically, and found no significant contribution.
A variety of methods to estimate the background in the signal box have been studied, including fitting and counting methods in various m ES and E sideband areas with different conditions. All methods yield results that are compatible within uncertainties.
The chosen method is model independent, is based on data only, and has a small systematic uncertainty. To estimate the background level in the signal box, five different sideband areas are used, as indicated in Fig. 2. The ratio R M est is the estimated ratio of the yield in the signal box to the yield in the M1 box. The expected background in the signal box (n exp bg ) is calculated by multiplying R M est by the yield in the M1 box. We estimate R M est as the mean of two ratios R U and R L , where R UL N U2L2 =N U1L1 , and where N X is the yield in box X. This assumes that the changes in the ratio R L , R M est , and R U are linear in E. To test our assumption of this linearity, we use MC samples and calculate the ratio R M by dividing the yield in the signal box by the yield in the M1 box. The relative difference between R M and R M est in MC samples is assigned as a systematic uncertainty. The estimated background is 1:75 1:38 0:36 (2:66 1:40 1:58) events for the e e ÿ ( ÿ ) mode, where the stated errors represent the statistical and systematic uncertainties, respectively.
The dominant source of systematic uncertainty in the signal yield is the calculation used for the signal MC [6]. The three theoretical input parameters, the Wilson coefficients C 7 , C 9 , and C 10 , used in the calculation are varied by 10%, as recommended by the authors of Ref. [6]. This variation changes the kinematics of the signal events and can thereby impact the detection efficiency. The largest relative change in signal efficiency by this variation is assigned as a systematic uncertainty.
The uncertainty in the branching fraction B4S ! B 0 B 0 [11] is taken as systematic uncertainty. We have studied e e ÿ ! ÿ decays in data to assess the systematic uncertainty in photon reconstruction. The systematic uncertainty from the lepton identification has been determined using an independent control sample of J= decays. The uncertainty in the number of B B events is 1.1% [12].
The systematic uncertainty related to an imperfect simulation of the ROE is studied using a control sample of B 0 ! J= K 0 S events. The same continuum background suppression requirements are applied on this sample and the signal efficiency is calculated. The relative difference in the signal efficiencies between data and MC samples is assigned as a systematic uncertainty.
The systematic uncertainty related to the tracking efficiency is determined from e e ÿ ! ÿ interactions, with one decaying leptonically and the other to three charged hadrons. All the contributions to the systematic uncertainties are added in quadrature and summarized in Table I.
After applying the selection criteria we find one event in the signal box for each mode, as shown in Fig. 3 and Table II. These numbers are compatible with the expected background for both modes.
An upper limit on the branching fraction is computed from where N UL is the 90% confidence level (C.L.) upper limit for the signal yield, determined by taking into account the one observed event in the signal box and the estimated background, using the frequentist method described in Ref. [13] including both statistical and systematic uncertainties, N B 0 is the number of neutral B mesons and sig is the signal reconstruction efficiency. The systematic uncertainties are included in sig . We used B4S ! B 0 B 0 49:4 0:8% [11]. The 90% C.L. branching fraction upper limits obtained are BB 0 ! e e ÿ < 1:2 10 ÿ7 and BB 0 ! ÿ < 1:6 10 ÿ7 .
In summary, a search for B 0 ! ' ' ÿ (' e or ) decays has been performed based on 320 10 6 B B events. We obtain 90% C.L. upper limits for the branching fractions of BB 0 ! e e ÿ < 1:2 10 ÿ7 and BB 0 ! ÿ < 1:6 10 ÿ7 , which represent the first limits placed on these decay channels. These are well above the SM expectations.
We are grateful for the excellent luminosity and machine conditions provided by our PEP-II colleagues, and for the substantial dedicated effort from the computing organizations that support BABAR. The collaborating institutions wish to thank SLAC for its support and kind hospitality. This work is supported by DOE and NSF (USA), NSERC  TABLE II. Summary of the results where n obs (n exp bg ) is the number of observed (expected background) events in the signal box, sig is the efficiency, N UL is the 90% C.L. upper limit for the signal yield, and B UL is the upper limit on the branching fraction at the 90% C.L. The stated uncertainties in n exp bg are statistical and systematic, and the uncertainty in sig is systematic.