DO-TH 93/29 OHSTPY-HEP-T-93-011 September 1993 4 9 9 1 A PARTON MODEL FOR INCLUSIVE n SEMILEPTONIC B MESON DECAYS a J 4 1 1 v C. H. Jina1, W. F. Palmerb2 and E. A. Paschosa1 5 5 2 aInstitut fu¨r Physik 1 Universit¨at Dortmund 0 D–44221 Dortmund, Germany 4 9 / bDepartment of Physics, The Ohio State University h p Columbus, Ohio 43210, USA - p e h : v i X r Abstract a The parton model for semileptonic B meson decays is studied with special attention to the decay distributions. We find that the spectra show dramatic variations when we introduce cuts on the hadronic energyorinvariantmassofhadrons. Resultsforbothb → uandb → c decays are presented. The detailed spectra may help to separate the two types of decays. 1Supported by Bundesministerium fu¨r Forschung und Technologie,05-6DO93P,Bonn, FRG and by the CEC Science Project no SC1-CT91-0729. 2Supported in part by the US Department of Energy under contractDOE/ER/01545- 605. The decays ofB-mesonshave beenstudiedextensively andhavebeenvery useful in extracting properties of the weak interactions. In particular, the b → u transition is still an active field of investigation. On the experimental side there is a large effort hindered by the difficulty of identifying the decays of B-mesons to light quarks. This difficulty of identifying b → u events forced investigators to concentrate on the end-point spectrum of the electron energy or study exclusive decays[1]. In this letter we present an extended study of the semileptonic B-decays. We are interested in the decay − B → X +e ν¯ (1) u for which we wish to give explicit distributions in several kinematic variables. With the decay spectrum completely determined by a single parameter in the distribution function, it may be compared to the competing process in − B → X + e ν¯. With the help of these distributions it should be possible c to test characteristics of the decay products with various kinematical cuts, and, in particular, collect events typical for this process. The hope is that with the results presented here the experimental analyses could incorporate many more events than those confined to the endpoint region of the electron spectrum. There are three models for the inclusive decays of reaction (1). The early modelofAltarellietal.[2], tobedenotedasACM,uses aspectatormodeland adistribution ofquarks within theB-mesondescribed by theFermimotionof the spectator quark. This approach treats the phase space effects correctly, butisrathercrude,astheauthorsstate[2],becauseitdependsonanunknown distribution function for the spectator quark. The model pictures the overall decay as the disintegration of the B-meson into the spectator quark plus the decay products of the heavy quark. Thesecond approach[3,4]usesthepartonmodelinaninfinitemomentum frame (IMF). The probability of finding a b-quark in a B-meson carrying a fraction x of the mesons momentum in the IMF is given by the distribution function f(x). The distribution function f(x) has the functional form sug- gested by theoretical arguments[5, 6, 7] and it peaks at large values of x. The kinematics for the decay of a heavy to a light quark involve the corre- lation of two currents at short distances, where the incoherence of the decay products is justified. This model has recently been improved and applied 2 to double differential inclusive decay distributions in reference [8] where a detailed treatment of the kinematics and the allowed physical regions can be found. A special feature of the parton model is the distribution of b quark in the B-meson, which provides a continuum spectrum for the mass of the recoiling hadrons. Here we only summarize the main results and point out new features before turning to applications. A third model, suggested very recently[9], improves the free quark model by including 1/m corrections from the heavy quark theory. These correc- Q tions turn out to be rather small. Obviously, the three models incorporate dynamical corrections in a different manner. In the parton model dynamic corrections are included through the experimentally determined distribution function and contributions from two scaling variables x+ and x−, to be de- scribed below. In the parton model[3, 4] the decay kinematics are different from those in deep inelastic scattering. As a consequence, the energy-momentum con- serving δ-function has now two roots, that is, it gives two scaling variables q0 ±|~q| q± x± = = (2) M M B B which are the light-cone variables of the current-momentum. Thus we have again scaling of the distribution function but now in terms of two variables. We define the kinematic variables: P = momentum of B-meson; P ,E = B e e momentum, energy of electron; P = momentum of neutrino; P = momen- ν X tum of hadrons; q = P +P = momentum of current; and M = invariant e ν X mass of the final hadronic system. In terms of these the decay is in general defined as: G2|V |2 d3P d3P dΓ = ub LµνW e ν. (3) (2π)5M µν 2E 2E B e ν The leptonic tensor has the simple form Lµν = 2(PµPν +PµPν −gµνP ·P +iεµν PαPβ). (4) e ν ν e e ν αβ e ν The hadronic tensor is defined in analogy to deep inelastic scattering, P P Pαqβ W (P ,q) = −g W (q2,q·P )+ Bµ BνW (q2,q·P )−iε B W (q2,q·P )+··· µν B µν 1 B M2 2 B µναβ M2 3 B B B (5) 3 wheretheintegrationoverthehadronicvariableshasalreadybeenperformed. Aftercarryingouttheintegrationsoverxwearriveatthestructurefunctions3 2 W1(q ,q·PB) = 2[f(x+)−f(x−)], (6) 4 2 2 W2(q ,q ·PB)/MB = M |~q|[x+f(x+)+x−f(x−)], (7) B 2 2 2 W3(q ,q ·PB)/MB = −M |~q|[f(x+)+f(x−)]. (8) B and hence the b → u triple differential decay rate dΓ G2|V |2q −E ub 0 e dE dq2ds = 8π3M |~q| {x+f(x+)(2Ee −MBx−)+(x+ ↔ x−)}. (9) e B This formula is simple and shows the dependence on the two light-cone vari- ables x± . The distribution function f(x) can be taken from references [5, 6, 7]. Its functionalformisvery similar tothefragmentationfunctionof abquark into a B-meson. It has also been argued, on physical grounds, that the distribu- tion and fragmentation functions for heavy mesons are the same (reciprocity relation)[3, 10, 11]. The fragmentation function was measured in several experiments and can be represented[11, 12] as x(1−x)2 f(x) = N (10) [(1−x)2 +ε x]2 p with ε a parameter andN thenormalization factor. We will use three values p ε = 0.003,0.006. p Finally, we must comment on the kinematic regions where the model is valid. Two criteria must be fulfilled. (i) The decay involves the correlation of two currents. We are more confident in applying the model when this distance between the two currents is close to the light-cone or at short distances. 3 Note that the structure functions in the present article are different from those in [4]. 4 (ii) The model should apply in the region where many particles are pro- duced. Thus, we must avoid the edge of the phase space where one or two pions are produced. To sum up, we feel rather confident in using the parton model for the − decay B → X +e ν¯ away from the edges of the phase space. For the decay u − B → X + e ν¯, we are less confident, because the distances involved are c − larger. Furthermore, for the decays B → X +e ν¯ we must keep the c-quark c mass, which modifies equation (9), and it makes a difference whether we choose the minimum value of the final hadronic mass as m or M , as we c D will show in figures 1b and 2a. As we mentioned, one of the purposes of this work is to present detailed spectra which can eventually be compared with experiments. The simplest parameter to measure is the electron energy. In addition one may be able to measure the total hadronic energy E = M −q . In Fig. 1a we show the X B 0 b → udoubledifferentialdecay rateinE andq . Thespectra showastriking e 0 dependence of q . Most of the events occur for 0.3M ≤ E ≤ 0.5M . 0 B X B For smaller hadronic energies the spectra shift to higher electron energies, E ≥ 2.0 GeV. This correlation of events may help to isolate b → u events. e − For comparison we also calculated the distribution for the B → X +e ν¯ c decay. In Fig. 1b we show the double differential decay rate in E and q . e 0 Comparing with Fig. 1a we note that this channel runs out of events at E = 2.25 GeV and the events above 2 GeV are very few. By making a cut e in q > M −M = 0.65M there are no events left for b → c. This can be 0 B D B used as another criterion for isolating b → u events. TheintegratedspectrumdΓ/dE isshowninFig.2afortwofragmentation e parameters ε . In the same figure we show the corresponding spectra for p b → c. For these curves we have chosen |V | = |V |. It is clear that the cb ub b → c spectrum is softer and a very small fraction lies at E > 2.0 GeV. We e show two spectra for the b → c decays, by varying the minimum value of the final hadronic mass. We have chosen two values (M ) = m = 1.5 GeV X min c and (M ) = M = 1.86 GeV. We notice that the total decay rate varies X min D considerablywhenwevary(M ) . Wefeelthatthefinalquarkmassshould X min 5 be replaced by the running charm quark mass4. As a result, the choice of (M ) influences the determinations of |V |. X min cb For comparison we calculated the b → u spectrum dΓ/dE in the ACM e model[2]. We show in Fig. 2b the ACM and parton spectra together. (In these spectra we include QCD radiative corrections which appear as a multi- plicative factor, asdescribed in[8]after thework of[2,13].) Wenote thatthe parton spectrum is lower and falls off less steeply than in the ACM model. Consequently, the values for V extracted from the end-point energy will be ub larger in the parton model with the experimental fragmentation function of eqn (10). Of course we can reproduce the free quark model by modifying the fragmentation function so that it has the limit δ(1−x) when the width parameter tends to zero, but this is not what is being measured for the frag- mentation function [11, 12, 7]. This means that the semileptonic rate of a free b quark is substantially greater than that of a b quark confined in a B meson. This is analogous to deep inelastic scattering, where moments of the distribution functions were measured to be smaller than one[14]. Thus the spectra in Fig. 2b differ not only in shape but also in overall scale. The double differential b → u decay rate dΓ/dE dM depends strongly e X on the invariant mass of the hadrons M . We show in Fig. 3 dΓ/dE dM for X e X various values of M . The decay spectrum is larger in the range 0.2M ≤ X B M ≤ 0.4M . The characteristic feature again is that the spectra shift to X B larger values of E as M decreases. After integration over E the sum of e X e the curves gives the mass distribution, which agrees with mass distributions published before[3] and is similar but not identical with the curves in ref. [15]. The parton model provides an interesting alternative for analysing inclu- sive semileptonic decays in important kinematical regions where the decay is dominated by short distance physics. We have produced decay spectra for a variety of interesting and physically accessible observables which probe the decay dynamics in a much more complete way than the simple electron spectrum. These spectra are obtained from a simple and compact formula (eq. (9)) based on a one parameter parton model. As high statistics data 4This topic is still under investigation. 6 become available this model may play an increasingly important role in sep- arating b → u from b → c decays and lead to a better determination of the V /V . In fact, the distribution function f(x) can be measured, in principle, ub cb in the b → u decays and then be used as an input to calculate the b → c spectra. CHJ wishes to thank Deutscher Akademischer Austauschdienst (DAAD) for financial support. WFP wishes to thank the North Atlantic Treaty Or- ganization for a Travel Grant. EAP wishes to acknowledge the hospitality of The Ohio State University, where this work was initiated. References [1] R. Fulton et al. (CLEO), Phys. Rev. Lett. 64 (1990) 16; H. Albrecht et al. (ARGUS), Phys. Lett. B234 (1990) 409; D. Besson (CLEO Collaboration), Proceedings of the XVI International Symposium on Lepton-Photon Interactions, Cornell University, Ithaca, NY, August, 1993 [2] G. Altarelli, N. Cabibbo, G. Corbo, L. Maiani and G. Martinelli, Nucl. Phys. B208 (1982) 365 [3] A. Bareiss and E.A. Paschos, Nucl. Phys. B327 (1989) 353 [4] A. Bareiss, Z. Phys. C53 (1992) 311 [5] M. Suzuki, Phys. Lett. B71 (1977) 139 [6] J.D. Bjorken, Phys. Rev. D17 (1978) 171 [7] G.L.Balayan,A.G.OganesyanandA.Yu.Khodzhamiryan, Sov. J. Nucl. Phys. 49 (1989) 697 [8] C. H. Jin, W. F. Palmer and E. A. Paschos, report DO-TH 93/21 and OHSTPY-HEP-T-93-011, September 1993 [9] I.I. Bigi, M. Shifman, N.G. Uraltsev and A. Vainshtein, Phys. Rev. Lett. 71 (1993) 496 7 [10] S.J. Brodsky, C. Peterson and N. Sakai, Phys. Rev. D23 (1981) 2745 [11] C. Peterson, D. Schlatter, J. Schmitt and P.M. Zerwas, Phys. Rev. D27 (1983) 105 [12] J. Chrin, Z. Phys. C36 (1987) 163 [13] G. Corbo, Nucl. Phys. B212 (1983) 99 [14] Particle Data Group, Phys. Rev. D45 (1992) Part II [15] V. Barger, C.S. Kim and R.J.N. Phillips, Phys. Lett. B251 (1990) 629 8 Figure Captions 1a. The b → u double differential decay rate dΓ/dE dq vs E for various e 0 e values of q in the rest frame of the B meson. We set m = 0,M = 0,M = 0 u π B 5.3 GeV,α = 0 and ε = 0.006. s p 1b. The b → c double differential decay rate dΓ/dE dq vs E for various e 0 e values of q in the rest frame of the B meson. For the solid lines we use 0 m = 1.5 GeV,M = 5.3 GeV,(M ) = M = 1.86 GeV,α = 0 and c B X min D s ε = 0.006. Thedashedlineshavethesameparametersexceptfor(M ) = p X min m = 1.5 GeV. For q = 0.4M the two curves coincide. c 0 B 2a. dΓ/dE for semileptonic B meson decays in the rest frame of the B e meson. The values for the parameters are m = 0,m = 1.5 GeV,M = u c B 5.3 GeV,M = 0,M = 1.86 GeV,α = 0 and various values of ε as shown. π D s p For the b → c decays the low and high lines correspond to (M ) = M X min D and (M ) = m , respectively. X min c 2b. The b → u spectrum in the ACM model and the parton model. Both curves include QCD radiative corrections with α = 0.24. s 3. The b → u double differential decay rate dΓ/dE dM vs E for various e X e values of M inthe rest frameof the B meson. The values for the parameters X are the same as in Fig. 1a. 9 This figure "fig1-1.png" is available in "png"(cid:10) format from: http://arXiv.org/ps/hep-ph/9401255v1