Search for the decay B0→a1±ρ

We present a search for the rare B -meson decay B 0 ! a (cid:1) 1 (cid:1) (cid:2) with a (cid:1) 1 ! (cid:2) (cid:3) (cid:2) (cid:4) (cid:2) (cid:1) . We use (cid:5) 110 (cid:1) 1 : 2 (cid:6) (cid:7) 10 6 (cid:3) (cid:5) 4 S (cid:6) ! B (cid:1) B decays collected with the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. We obtain an upper limit of 30 (cid:7) 10 (cid:4) 6 (cid:5) 90%C : L : (cid:6) for the branching fraction product B (cid:5) B 0 ! a (cid:1) 1 (cid:1) (cid:2) (cid:6) B (cid:5) a (cid:1) 1 ! (cid:2) (cid:3) (cid:2) (cid:4) (cid:2) (cid:1) (cid:6) , where we assume that the a (cid:1) 1 decays exclusively to (cid:1) 0 (cid:2) (cid:1) .

We present a search for the rare B-meson decay B 0 ! a 1 with a 1 ! ÿ . We use 110 1:2 10 6 4S ! B B decays collected with the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. We obtain an upper limit of 30 10 ÿ6 90%C:L: for the branching fraction product BB 0 ! a 1  In the standard model, CP-violating effects in the B-meson system arise from a single phase in the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix [1]. The decay B 0 ! a 1 proceeds via a b ! u ud transition [2], and interference between direct decay and decay after B 0 B 0 mixing results in a time-dependent decay-rate asymmetry that is sensitive to the angle argÿV td V tb =V ud V ub [3] in the unitarity triangle of the CKM matrix. An additional motivation for studying B 0 ! a 1 is that this is a significant background to B ! decays, e.g. [4 -8], which currently provide the most accurate measurement of . The ARGUS experiment previously searched for the decay B 0 ! a 1 , which resulted in an upper limit of BB 0 ! a 1 < 3:4 10 ÿ3 (90% C.L.) [9]. This paper presents the result of a search for B 0 ! a 1 with a 1 ! ÿ , where we assume that the a 1 decays exclusively to 0 . A theoretical prediction of the branching fraction BB 0 ! a 1 Ba 1 ! 3 has been made by Bauer, Stech and Wirbel [10] within the framework of the factorization model. They predict a value of 43 10 ÿ6 , assuming jV ub =V cb j 0:08.
The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC during the years 2003-2004. This represents a total integrated luminosity of 100 fb ÿ1 taken at the 4S resonance (on-peak), corresponding to a sample of 110 1:2 million B B pairs. An additional 21:6 fb ÿ1 of data, collected at approximately 40 MeV below the 4S resonance (off-peak), were used to study background from e e ÿ ! q q q u; d; s; c continuum events. The BABAR detector is described in detail elsewhere [11]. Surrounding the interaction point is a silicon vertex tracker (SVT) with 5 double-sided layers which measures the impact parameters of charged particle tracks in both the plane transverse to, and along the beam direction. A 40 layer drift chamber (DCH) surrounds the SVT and provides measurements of the transverse momenta for charged particles. Both the SVT and the DCH operate in the magnetic field of a 1.5 T solenoid. Charged hadron identification is achieved through measurements of particle energy-loss in the tracking system and the Cherenkov angle obtained from a detector of internally reflected Cherenkov light. A CsI(Tl) electromagnetic calorimeter (EMC) provides photon detection, electron identification, and 0 reconstruction. Finally, the instrumented flux return of the magnet allows discrimination of muons from pions.
We reconstruct B 0 ! a 1 ÿ candidates from combinations of a 1 ! ÿ and ÿ ! 0 ÿ candidates. The a 1 1260 ! 3 decay proceeds mainly through the intermediate states and [12]. We do not distinguish between the dominant P-wave and S-wave in the channel ÿ . The Monte Carlo (MC) signal events are simulated as B 0 decays to a 1 1260 ÿ with a 1 ! 0 using the GEANT4-based [13] BABAR MC simulation. Possible contributions from B 0 decays to a 2 1320 ÿ and 1300 ÿ are investigated. We only consider events that have a minimum of one 0 and four charged tracks, where the charged tracks are required to be inconsistent with lepton, proton and kaon hypotheses.
We form 0 ! candidates from pairs of photon candidates that have been identified as localized energy deposits in the EMC that have the lateral energy distribution expected for a photon. Each photon is required to have an energy E > 50 MeV, and the 0 is required to have an invariant mass of 0:10 < m < 0:16 GeV=c 2 .
The ÿ mesons are formed from one track that is consistent with a ÿ and the aforementioned 0 candidate. The candidate ÿ is required to have an invariant mass of 0:5 < m ÿ < 1:1 GeV=c 2 . We also constrain the cosine of the angle between the 0 momentum and the direction opposite to the B 0 in the ÿ rest frame ( cos ÿ ) to be between ÿ0:9 and 0.98. This removes backgrounds which peak at the extremes of the distribution where the signal reconstruction efficiency also falls off.
We form the a 1 candidate from combinations of three charged pions. We first form a 0 ! ÿ candidate from two oppositely charged tracks. This combination is required to have an invariant mass of 0:4 < m 0 < 1:1 GeV=c 2 . The a 1 candidate is then formed by adding another charged track to the 0 , and requiring that the mass of the a 1 satisfies 0:6 < m a 1 < 1:5 GeV=c 2 . The vertex of the B-candidate is constrained to originate from the beam spot. In order to reduce background from continuum events we require that j cos T j < 0:7, where T is the angle between the B thrust axis and that of the rest of the event (ROE).
We use two kinematic variables, m ES and E, in order to isolate any signal. We define the beam-energy substituted is the e e ÿ center-of-mass (CM) energy. The second kinematic vari-able, E, is the difference between the B-candidate energy and the beam energy in the CM frame. We require m ES > 5:25 GeV=c 2 and ÿ0:15 < E < 0:1 GeV.
Additional separation between signal and continuum is obtained by combining several kinematic and topological variables into a Fisher discriminant (F ) [14]. The variables L 0 , L 2 , and j cos TR j, and the output of a multivariate tagging algorithm [15] are used as inputs to F . L 0 and L 2 are defined as where the sum is over the ROE, p i is the particle momentum in the CM frame. i is the angle of the particle direction relative to the thrust axis of the B-candidate, and cos TR is the cosine of the angle between the B thrust axis and the beam axis. The multivariate tagging algorithm identifies the flavor of the other B in the event to be either a B 0 or B 0 . The output of this algorithm is ranked into categories of different signal purity.
We expect the polarization of the a 1 ÿ final state to be predominantly longitudinal, as was found in the similar decay B ! [4 -8]. We have used both longitudinal and transverse polarized signal MC simulated data in this analysis. After applying the selection cuts above, we have 2.8 (2.3) longitudinal (transverse) polarized signal MC simulated data candidates per event.
We define as self-cross-feed (SCF) the set of candidates that were incorrectly reconstructed from particles in events that contain a true signal candidate. We select one B candidate per event in which the mass of the reconstructed 0 is closest to that of the true 0 mass [16]. Choosing the candidate using the 0 mass reduces the SCF fraction by 18% relative to a random selection. To avoid potentially biasing our final result, we do not use information from the 0 meson in the remainder of the analysis. After all selection cuts have been applied, the longitudinal and transverse signal SCF fractions are 0.58 and 0.42, respectively. The selection efficiency of longitudinal (transverse) signal is 9.44% (10.15%).
Besides the continuum background we also have background from B decays. We divide the B-background into the following four categories according to B-meson charge and the charm content of the final states: (i) B 0 ! charm, (ii) B 0 ! charmless, (iii) B ! charm and (iv) B ! charmless. From large samples of inclusive MC simulated data we expect 2394, 424, 3281 and 215 events of these background types, respectively. In addition, a number of exclusive B-background modes that have a similar final state to the signal were studied. This includes those that have an intermediate a 1 meson in the decay. None of these modes were seen to have a significant efficiency after the selection cuts had been applied.
We perform an extended unbinned maximum likelihood fit to the data. The likelihood model has the following types: (i)-(iv) the four aforementioned inclusive B-background categories, (v) true signal, (vi) SCF signal and the (vii) e e ÿ ! q q q u; d; s; c continuum background. The probability density function (PDF) for each event i has the form P i;c P i;c m ES ; E; F ; m a 1 ; m ÿ ; cos ÿ . From these individual PDFs the total likelihood is constructed. The parameters n and n 0 are the numbers of selected on-peak events, and the sum of the yields N c , where c is one of the seven types in the likelihood model. Correlations between the m ÿ and cos ÿ variables are taken into account for the real (T) and fake combinatorial (F) ÿ candidates. All other correlations between fit variables are seen to be small. However, the effect of ignoring them results in a bias on the fitted signal yield which is discussed below. Each of the signal distributions has a signal yield and a polarization fraction that are denoted by N sig (N 0 sig ) and f L (f 0 L ), for the true (SCF) signal, such that the sum of the true and SCF signal is described by where f 0 L is an effective parameter. The continuum yield,  [18], which is allowed to vary in the fit, Pn Polynomial of order n, BW Breit-Wigner, helicity cos 2 ÿ or sin 2 ÿ depending on partial wave which is modified by a quadratic acceptance function, BG m-hel Background cos ÿ and m ÿ PDF of Eq. (4), off-peak PDF taken from off-peak data, and 1D smoothed 1D histogram. N sig , N 0 sig , and the parameters of the continuum m ES and E PDFs are allowed to vary in the fit. Under the assumption that no significant signal is observed, the value of f L is fixed to 1.0 in the fit. The value of f 0 L is also fixed to 1.0 in the fit since it is highly correlated with N sig , N 0 sig and f L . Only the fitted value of N sig is used to derive the final result. We also fix the B-background yields to the aforementioned values.
The PDFs used for each component are given in Table I. The signal and B-backgrounds are parameterized using MC. We use a nonparametric smoothing algorithm [19] when defining some of the background PDFs (as indicated by the abbreviation NP). We account for the difference between F and T ÿ ! ÿ 0 distributions in the background using the product of 1D PDFs, denoted in Table I as 'BG m-hel', such that P i m ÿ ; cos ÿ 1 ÿ a T P F;i m ÿ P F;i cos ÿ a T P T;i m ÿ P T;i cos ÿ ; (4) where a T is the fraction of T events. The continuum shape for cos ÿ (m ÿ ) is derived from off-peak (on-peak) data. The true ÿ resonance Breit Wigner shape uses m ÿ 0:77 GeV=c 2 , and ÿ 0:150 GeV [16]. The parametrizations used for this PDF are summarized in Table II. The results from the fit are N sig 90 38stat, N 0 sig 42 98 and a continuum yield of 25798 182 events. The bias on the fitted signal yield is evaluated from ensembles of mock experiments obtained from samples of signal MC events combined with background events generated from the PDFs. The bias is found to be 22 events (24%), resulting in a corrected signal yield of 68 38stat. We use the event-weighting technique described in Ref. [20] to compare the true signal and continuum PDF shapes to the data. For each of the plots in Fig. 1, we perform a fit excluding the variable being plotted, and the covariance matrix is used to determine the weight that each event is true signal or q q background. The resulting distributions can be directly compared with the PDFs used in the fit. The distributions shown in Fig. 1 are not corrected for fit bias, and the uncertainty on each of the data points is statistical. No change in signal yield is seen when a 2 1320 ÿ and 1300 ÿ components are included in the fit. Table III summarizes the systematic uncertainties on the signal yield. Each entry in the table indicates one systematic effect, and the contributions are added in quadrature to give the total presented. The uncertainty due to PDF parametrization is evaluated by variation of both the signal and the background PDF parameters within their uncertainties about their nominal values. The uncertainty from the continuum m ES and E PDFs which are allowed to vary in the fit, are only included in the quoted statistical  1 (color online). The true signal (top) and continuum (bottom) distributions for (left to right) m ES , E, m a 1 using the weighting technique described in Ref. [20]. The points represent the weighted data, and solid curves represent the corresponding PDFs.  Table I.
uncertainty. We assign a systematic uncertainty due to fit bias, evaluated as half of the fit bias correction on the signal yield. To validate the expected B-background yields, and to assign a systematic uncertainty we perform a number of cross-checks in which we allow the background yields to vary in turn when fitting the data. We use a control sample of B ! D events to determine the systematic uncertainty in the fraction of SCF signal events. The effect of exclusive B meson decays to final states including a 1 -mesons were evaluated using ensembles of mock experiments. In particular, the systematic uncertainty on the signal yield from neglecting B ! a 1 a 1 modes in the fit is 6 events. We assign a systematic uncertainty from using a relativistic Breit-Wigner with a Blatt-Weisskopf form factor with a range parameter of 3:0 GeV ÿ1 for the a 1 meson line shape. In the fit we assume that the a 1 meson width, ÿ a 1 , is 400 MeV. We evaluate a systematic uncertainty due to this assumption by varying ÿ a 1 over the experimentally allowed range: 250 -600 MeV [16]. The difference in the distribution of F between data and MC is evaluated with a large sample of B ! D ? decays. The systematic uncertainties that contribute to the branching fraction only through the efficiency come from charged particle identification (6.0%), 0 meson reconstruction (3.0%), tracking efficiency (3.2%), and the number of B meson pairs (1.1%). The systematic error contribution from MC statistics is negligible.
When the fit bias correction of ÿ22 events is applied to the signal yield, and one accounts for systematic uncertainties, the significance of the result is 0.95 standard deviations. Figure 2 shows the distribution of ÿ lnL=L max for the fit, with and without these systematic errors. L max is the value of the likelihood corresponding to the nominal fit result. The branching fraction value for the fit-bias-corrected signal yield of 68 38stat 45 ÿ56 syst is BB 0 ! a 1 ÿ Ba 1 ! ÿ 15:7 8:7stat 10:3 ÿ12:8 syst 10 ÿ6 . This assumes that f L 1:0 and that the branching fraction of a 1 ! ÿ 0:5. As the signal yield obtained is not significant, we calculate the upper limit x UL , by integrating the likelihood function (including systematic uncertainties) from 0 to x UL , for different physically allowed values of f L , such that the C.L. of the upper limit is 90%. As the signal efficiency is a function of f L , we report the most conservative upper limit obtained, which corresponds to f L 1:0. On doing this, an upper limit of 30 10 ÿ6 90%C:L: is obtained. We have performed a search for the decay B 0 ! a 1 in a data sample of 100 fb ÿ1 . After correcting for fit bias and accounting for systematic uncertainties, the signal yield is 68 38stat 45 ÿ56 syst events, with a significance of 0:95. As there is no significant evidence for a signal, we place an upper limit of 30 10 ÿ6 90%C:L: on BB 0 ! a 1 ÿ Ba 1 ! ÿ , where we assume that the a 1 decays exclusively to 0 . Assuming Ba 1 ! ÿ is equal to Ba 1 ! 0 0 , and Ba 1 ! 3 is equal to 100%, we obtain BB 0 ! a 1 ÿ < 61 10 ÿ6 90%C:L:. This upper limit corresponds to a significant improvement over the previous bound and is compatible with theoretical expectations [10]. This result is a significant improvement in constraining an important B background contribution in B ! decays.
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.  2 (color online). The ÿ lnL=L max distribution from the fit to data with f L 1:0 (fixed). This distribution has been corrected for fit bias. The solid curve is for statistical errors only, and the dashed curve includes systematic errors.