Hadron multiplicity in e+e-events induced by top quark pairs at the ILC energy
© Kisselev and Petrov 2008
Received: 05 December 2007
Accepted: 04 April 2008
Published: 04 April 2008
The average charged hadron multiplicity in the e+e- events with the primary -pair at the collision energy 500 GeV, as well as the average multiplicity of charged hadrons from the top quark are calculated in QCD to be 86.7 ± 1.11 and 41.0 ± 0.54, respectively.
PACS Codes: 14.65.Ha, 13.66.Bc, 12.38.Bx
Experiments at LEP and SLAC revealed, besides other important results, quite interesting feature of the hadron multiple production dependent on the mass of the "primary" (anti)quarks which launch the process of the QCD evolution. It appeared that differences between the light and heavy quark-induced multiplicities become energy-independent. QCD calculations describe the phenomenon quite well.
Certainly, LEP could not give the information on the events induced by the top quarks. Recent discussions of the ILC project give us occasion to provide QCD predictions concerning the hadron multiple production in the events with primary t-quarks.
We also theoretically calculated the average hadronic multiplicity from the top quark:
n t = 41.03 ± 0.27.
Both values correspond to the average value of the top mass m t = 170.9. Everywhere below, it is assumed that we deal with average multiplicities of charged hadrons.
The paper is organized as follows. In order to make our calculations of the hadron multiplicity in top quark events more easy for understanding, we consider first the multiple hadron production in () events. The hadron multiplicities in e+e- events associated with the -pair production are calculated in Section 3 in the framework of perturbative QCD. In Section 4 the numerical estimations and our main results are presented.
2 Hadron multiplicity in e+e- annihilation associated with or -pair production
where q means a type of quarks produced in the process of e+e- annihilation into hadrons at the collision energy W. In what follows, the notation q = Q (heavy quark) will mean charm or beauty quark, while the notation q = l (light quark) will correspond to a massless case (when a pair of u, d or s-quarks is produced, whose masses are assumed to be equal to zero). The top quark production (q = t) will be studied in Sections 3 and 4.
The first term in the r.h.s. of Eq. (3), 2n q , is the multiplicity of primary (anti)quark of the type q (i.e. the multiplicity from the leading hadron which contains this (anti)quark). It is taken from an analysis of the data (2n c = 5.2, 2n b = 11.1 , and 2n l = 2.4 ).
The quantity n g (k2) in (3) is the mean multiplicity of the gluon jet with a virtuality k2, for which we will take a QCD-based parametric form, with parameters fit to data, while E q (k2/W2) is the inclusive spectrum of the gluon jet emitted by primary quarks. It was explained in detail in Ref.  that one should not consider this mechanism of hadron production via gluon jets as due to "a single cascading gluon". That quantity E(k2/W2) is an inclusive spectrum of the gluon jets is seen, e.g., from the fact that the average number of jets ∫dk2/k2 E q (k2/W2) ≠ 1.
Q0 is a phenomenological parameter denoting the scale at which "precon-finement" of the off-shell partons occurs (as explained in Ref. ).
In particular, (Y) means the multiplicity of hadrons in light quark events, while (Y) denotes the multiplicity of hadrons in a process when a pair of the heavy quarks is produced.
is the following. It describes the average number of hadrons produced in virtual gluon jets emitted by the primary quark and antiquark of the type q. In other words, it is the multiplicity in event except for multiplicity of the decay products of the primary quarks at the final stage of hadronization (the terms 2n q in (7)).
For the massless case, the function E ≡ E l was calculated in our paper . In terms of variable
σ = exp(-η),
At W ≫ m Q , one can neglect small power-like corrections O(m2/W2). In such a case, the quantity δ Ql is defined by 
δ Ql = 2(n Q - n l ) - ΔN Q (Y Q ),
are introduced. The lower limit of integration in Eq. (15), -ln(W2/), is taken -∞ because of the fast convergence of the integral at negative y.
Let us use another dimensionless variable
ρ = exp(-y).
with the multiplicity difference ΔN Q defined by Eq. (15).
Our calculations  of the multiplicity differences δ Ql = - (Q = b, c) with the use of formula (23) appeared to be in a good agreement with the data. Recently we have reconsidered the QCD upper limit on quantity δ bl  which appeared to be very close to all present experimental data on δ bl . The data on N ll as well as on δ bl at different energies corrected for detector effects as well as for initial state radiation were recently cited in .
3 Hadron multiplicity in e+e- annihilation associated with -pair production
The goal of this paper is to calculate , the average multiplicity of hadrons produced in e+e- events with the primary pair. We consider the case when the top (antitop) decay mode is pure hadronic. As a byproduct, we will calculate n t , the hadron multiplicity of the on-shell top decay products.
where denotes the virtual top quark (antiquark).
The factorization of the matrix element (24) means that there is no significant space-time overlap in the decay products of the on-shell t and -quarks. Note that the off-shell t and -quarks fragment into hadrons through the emission of the gluon jets in a coherent way (the first term in the r.h.s of Eq. (24)). The QCD non-singlet evolution of the primary virtual t-quark is very slow because the difference of virtualities in logarithmic scale is very small down to the top quark mass. In other words, the virtual t-quark becomes "real" after just a few gluon radiations.
The effect of possible color reconnection was investigated by comparing hadronic multiplicities in e+e- → W+W- → and e+e-→ W+W- → events. No evidence for final state interactions was found by measuring the difference [8, 9]. From the space-time point of view W bosons and t-quarks behave in a similar way, ie. the latter manage to cover the distance Δ l ~ 1/Γ t , where Γ t is the full width of the top. Since Γ t ≃ Γ w , we expect no interference effects in the decays of the on-shell t and -quarks.
This formula has been derived by calculating QCD diagrams in the first order in the strong coupling constant (see our comments after Eq. (3)). In the massless case (m = 0), we immediately come to the function E(k2/q2) (10), while by neglecting small corrections O(m2/q2), one can derive (after variables are properly changed) the explicit form of the function ΔE Q (k2/m2) = E Q (q2, k2, m2) - E(q2, k2) (19). In our case (q2 = W2, m = m t ) we will estimate the integrals in Eqs. (26), (27) numerically (for details, see Section 4).
Thus, the multiplicity n t is a sum of three terms:
n t = n W + n tb + n b .
The quantity n b is experimentally measurable one . The first term in Eq. (29), n W , is the hadron multiplicity of the W boson decay products. The second term, n tb , is the hadron multiplicity in the gluon jets emitted by the on-shell top quark before its weak decay as well as by the bottom quark after the top decay.
3.1 Multiplicity of W boson decay products
The W+ boson can decay either into two light quarks ( and pairs) or into () pair. The former case is treated analogously to the light quark event in e+e- annihilation taken at the collision energy W = m W . Here we will study the latter case.
The function (Y) in (32) is the hadronic multiplicity in light quark events. Thus, we need to find an expression for ΔN cl . Note that the formulae (15), (19) from Section 2 correspond to the case when a pair of heavy or pair of light quarks is produced. Now we have to study the case when hadrons are produced in association with a single heavy quark (namely, c-quark) and one light quark.
3.2 Multiplicity of top and bottom decay products
As was already said above, the on-shell top quark can emit jets before it weakly decays into W+b. After the weak decay of the top, the off-shell b-quark "throws off" its virtuality by emitting massive gluon jets. The fragmentation of these massive gluon jets into hadrons results in the average hadron multiplicity n tb .
with k2 being the gluon jet invariant mass, (m t - m W - m b )2 its upper bound.
3.3 Associated multiplicity of hadrons in events
where n tb is defined above (41). In what follows, we will use the value
δ cl = 1.03 ± 0.34
from Ref. .
4 Numerical estimates of hadron multiplicities
where c = 1.63, and k2 is the invariant mass of the jet. We have got the following values of the parameters:
a = 3.89, b = 0.01, Q0 = 0.87 GeV.
m t = 170.9 ± 1.8 GeV.
As for the bottom quark, its pole mass is quoted in  to be m b = 4.7 – 5.0 GeV. We will take the average value
m b = 4.85 ± 0.15 GeV.
By using our fit, we obtain (m W ) = 19.09. Then we get (see Eqs. (40), (49), (41)):
n W = 19.34 ± 0.10,
n tb = 16.14 ± 0.24.
The error in Eq. (54) is defined by that of the multiplicity difference δ cl (49), while that in Eq. (55) comes from uncertainties of the quark masses m t (52) and m b (53).
Our result (54) is in a nice agreement with the experimental values from Ref. ,
n W = 19.3 ± 0.3 ± 0.3,
and Ref. ,
n W = 19.44 ± 0.13 ± 0.12.
Now let us calculate the associated hadron multiplicity in event (48) at fixed energy W = 500 GeV. To do this, we need to estimate the multiplicity N t (W, m t ) by using formulae (26) and (27):
N t (W = 500 GeV) = 4.61 ± 0.11.
The errors in (58) come from top quark mass errors. It follows From Eqs. (61), (58) that = 86.67 ± 0.55.
It appeared that the data on the average multiplicity in light quark events can be fitted well by using this expression (with A = 4.21, B = 0.012 and Q0 = 0.93 GeV). In particular, we have obtained the following average values for the hadronic multiplicities: n W = 19.52, n tb = 16.43, N t = 4.59. Thus, theoretical uncertainties can be estimated to be 0.47 and 0.96 for n t and , respectively.
Taking into account the phenomenological value of n b ,
n b = 5.55 ± 0.09,
we obtain from (47):
n t (t → hadrons) = 41.03 ± 0.54.
All estimations (63)–(65) correspond to the collision energy W = 500 GeV. We can mention the estimation of the hadron multiplicity from Ref. , , which was obtained for W = 390 GeV and m t = 175 GeV. For the same values of W and m t , our formulae give .
The formulae (61)–(65) is our main result. We hope that the hadron multiplicities of the top decay products (Eqs. (61) and (62)) will be measured at the LHC.
Here we present some formulae for the case when the W boson decays into hadrons via production of (or ) pair. Since the total width of the W boson, Γ W , is much less than its mass, and the hadron multiplicity is a smooth function of energy, we will use zero width approximation and take the multiplicity at W = m W . It can be shown that the account of the W boson width results in corrections which are numerically small (less than 1.6%, see Appendix B).
where m t and m b are masses of the top and beauty quark, respectively. In what follows, we will neglect power corrections of the type O(m c /m t ) and O(m b /m t ).
Note that the parts of tensors Π μν and D μν antisymmetric in indices (as the last term in Eq. (A.1)) give no contribution after integration in lk.
The analytic expressions for the functions D1(k2, qk), D2(k2, qk) in (A.2) are rather complicated to be shown here. That is why we present only the final results of our QCD calculations based on the formulae of this Appendix (see Eqs. (34), (35) in the main text).
In this Appendix we will demonstrate that the account of the W boson width results in only small corrections to the hadronic multiplicities.
In both <n h > (B.2) and N (B.3) the factor m W Γ W is introduced, while common constants are omitted.
we do obtain that the mean multiplicity is equal to n h ().
The numerical calculations with the use of formulae (B.2), (B.3) result in the following values:
<n w > = 19.04,
<n tb > = 16.37.
Thus, the account of non-zero width of the W boson slightly changes the average value of the multiplicities. Namely, n W (54) has gone down by 0.3, while n tb (55) has gone up by 0.23, but their sum remains almost unchanged.
We are thankful to the referee for his comments and critical remarks that helped us to improve the presentation of some our results.
- Petrov VA, Kisselev AV: Z Phys C. 1995, 66: 453-10.1007/BF01556371.View ArticleADSGoogle Scholar
- Petrov VA, Kisselev AV: Eur Phys B (Proc Suppl). 1995, 39 (2-3): 364-10.1016/0920-5632(95)00102-F.View ArticleADSGoogle Scholar
- Kisselev AV, Petrov VA: Eur Phys J C. 2007, 50: 21-10.1140/epjc/s10052-006-0197-z.View ArticleADSGoogle Scholar
- Schumm BA, Dokshitzer YuL, Khoze VA, Koetke DS: Phys Rev Lett. 1992, 69: 3025-10.1103/PhysRevLett.69.3025.View ArticleADSGoogle Scholar
- Chrin J, DELPHI Collaboration, et al: Proceedings of the 27th International Conference on High Energy Physics: 20–27. Edited by: Bussey PJ, Knowles IG. 1994, 893-July ; Glasgow, UKGoogle Scholar
- Bassetto A, Ciafaloni M, Marchesini G: Nucl Phys B. 1980, 163: 477-10.1016/0550-3213(80)90413-7.View ArticleADSGoogle Scholar
- Dokshitzer YuL, Fabbri F, Khoze VA, Ochs W: Eur Phys J C. 2006, 45: 387-10.1140/epjc/s2005-02424-5.View ArticleADSGoogle Scholar
- Abbiendi G, OPAL Collaboration, et al: Phys Lett B. 1999, 453: 153-10.1016/S0370-2693(99)00309-3.View ArticleADSGoogle Scholar
- Abreu P, DELPHI Collaboration, et al: Eur Phys J C. 2000, 18: 203-10.1007/s100520000528.View ArticleADSGoogle Scholar
- Pleier M-A: Talk presented at the "XXVII Physics in Collision": 26–29. 2007, June ; Annecy, FranceGoogle Scholar
- Review of Particle Physics (Particle Data Group): Nucl & Part Phys G. 2006, 33: 1-
- Ballestrero A, Khoze VA, Maina E, Moretti S, Stirling WJ: Z Phys C. 1996, 72: 71-10.1007/s002880050224.ADSGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.