Charged Multiplicities in Z Decays into u, d, and s Quarks

About 4.4 million hadronic decays of Z bosons, recorded by the OPAL detector at LEP at a centre-of-mass energy of around sqrt(s) = 91.2 GeV, are used to determine the mean charged particle multiplicities for the three light quark flavours. Events from primary u, d, and s quarks are tagged by selecting characteristic particles which carry a large fraction of the beam energy. The charged particle multiplicities are measured in the hemispheres opposite to these particles. An unfolding procedure is applied to obtain these multiplicities for each primary light quark flavour. This yields= 17.77 +- 0.51 +0.86 -1.20,= 21.44 +- 0.63 +1.46 -1.17,= 20.02 +- 0.13 +0.39 -0.37, where statistical and systematic errors are given. The results forandare almost fully statistically anti-correlated. Within the errors the result is consistent with the flavour independence of the strong interaction for the particle multiplicities in events from the light up, down, and strange quarks.


Introduction
The flavour independence of the strong coupling is a fundamental property of quantum chromodynamics (QCD). A breaking of the flavour symmetry should only occur due to calculable mass effects. These mass effects have been observed for bottom quarks [1,2] using event shapes and jet rates in the final state of electron-positron annihilation into bb. Another observable which can be employed to test flavour independence is the multiplicity of charged hadrons in jets originating from quarks of a specific flavour. According to the local parton hadron duality hypothesis (LPHD) [3], the particle multiplicity is related to the gluon multiplicity inside a jet which depends on the value of the strong coupling constant [4]. A dependence of the multiplicity of charged hadrons on the quark flavour has been found for heavy quarks [5]. Up to now only few measurements exist on the flavour independence of the strong interaction in the light quark sector (up, down, strange quarks) [2]. Within the large uncertainties due to the limited statistics, flavour independence of the strong interaction is supported.
This paper presents a new, high statistics investigation of the flavour dependence of the strong interaction. It is based on the mean charged multiplicity determined separately for events of primary up, down, and strange quarks in e + e − annihilation at centre-of-mass energies close to the mass of the Z boson. To identify the flavour of the primary quark in the Z decay, the leading particle effect is exploited [6]. This assumption of a correlation between the flavour of the primary quark and the type of the hadron carrying the largest momentum has recently received further support by measurements of the SLD [7] and OPAL collaborations [8].
The study presented in this paper uses three different selections of events with leading K 0 S and K ± mesons, and highly energetic stable charged particles, denoted "high x E " in the following. Due to the different fractions of primary up, down, and strange quarks in these samples a statistical decomposition of the contributions from each of the three quark flavours is possible. The leading K 0 S and K ± meson selections, even though dominantly stemming from primary strange quarks, yield the separation of up and down quark events since a leading K 0 S (K ± ) meson is rarely formed from a primary up (down) quark while both K 0 S and K ± are equally likely to be produced in a strange quark event. Together with the high x E selection, containing up, down, and strange quark events in approximately equal proportions, properties like the mean charged particle multiplicity in events of each of these three light quark flavours can be determined by statistical unfolding. The mean charged multiplicity is determined from all long-lived charged particles in the hemisphere opposite to the leading particle where the two hemispheres of an event are defined by the plane perpendicular to the thrust axis which contains the main interaction point. Choosing the opposite hemisphere minimizes the bias of the measured multiplicity due to the high energy of the leading particle.
Before the three selections of leading particles are presented in Section 3, the OPAL detector and the data samples used for this study are briefly introduced in Section 2. The determination of the charged particle multiplicities is discussed in Section 4. The investigation of various sources of systematic uncertainty is detailed in Section 5. Section 6 presents a cross-check of the analysis using flavour fractions obtained from the data. The results of this analysis are presented in Section 7.

The OPAL detector, data, and Monte Carlo simulation
The OPAL detector has been described in detail elsewhere [9]. The analysis presented here relies mainly on the reconstruction of charged particles in the large volume jet chamber whose performance has been presented in [10]. A solenoidal coil surrounds all tracking detectors. It provides a field of 0.435 T along the beam axis 1 . Tracks of charged particles are reconstructed with up to 159 space points. Their momenta in the plane transverse to the beam-axis, p xy , can be determined to a precision of σ pxy /p xy = 0.02 2 + (p xy · 0.0015/GeV) 2 . This resolution degrades towards the acceptance boundary of | cos θ | ≈ 0.98. From the specific energy loss dE/dx, which is measured from up to 159 samples with a resolution of σ(dE/dx)/(dE/dx) ≈ 0.035, the type of the charged particle can be identified over a wide momentum range.
The analysis is based on data recorded with the OPAL detector between 1990 and √ s around 91.2 GeV (LEP I). The events considered for this study were preselected by a standard selection for high multiplicity events [11] which relies on a minimum number of measured tracks of charged particles in the tracking detectors and clusters of energy deposited in the electromagnetic calorimeter. The remaining background of two-photon processes and τ -pair events are estimated to be 0.07 % and 0.11 %, respectively.
Tracks to be used for the reconstruction or identification of the leading particles and for the determination of the charged multiplicity were required to pass the quality selection cuts as detailed in [12]. These tracks had to have at least 20 hits in the jet chamber and a closest approach to the interaction point in the plane perpendicular to the beam axis of |d 0 | < 5 cm. The track momenta had to be between 0.1 and 65 GeV. Events with τ -pairs are further suppressed by requiring at least seven tracks.
The Monte Carlo simulation comprised about 8.5 million hadronic events simulated by version 7.4 of the JETSET program [13] which has been tuned to describe the OPAL data [14]. The generated events were passed through a detailed simulation of the OPAL detector [15] and processed using the same reconstruction and selection algorithm as the measured data.

Selection of flavour enriched data samples
Every event is divided into two hemispheres by the plane perpendicular to the thrust axis containing the interaction point. In each hemisphere we searched for K 0 S and K ± mesons, and highly energetic stable charged particles, whose energy fraction x E = 2E/ √ s is above a certain threshold. Due to the cuts, in particular on x E , these three samples are enriched in events from up, down, and strange quarks in different proportions (see Table 1). The x E selection cuts were chosen to have samples of comparable size.
3.1 K 0 S selection K 0 S mesons were selected via their decay into two charged pions using the procedure described in [16]. It was adapted for large momenta of the K 0 S mesons by dropping the specific d 0 cut imposed on the two pions in the default procedure. Any two oppositely charged particles were combined. The invariant mass m ππ was calculated assuming the pion mass for the particles. K 0 S candidates were required to have |m K 0 S − m ππ | < 60 MeV using the mass of a K 0 S meson, m K 0 S , given in [17]. Photon conversions were rejected by demanding m ee > 100 MeV if the electron mass is assigned to both tracks. The impact of detector effects close to the acceptance boundary is reduced by restricting the polar angle of K 0 S candidates to | cos θ | < 0.9. To enrich the sample in primary strange and down quarks and to suppress up, charm, and bottom events, a requirement on the scaled energy of the K 0 S candidate of x E > 0.4 was applied. Candidates with x E > 1.07 were rejected. This takes into account the 7 % momentum resolution for particles with the beam energy.   Table 1 gives the number of tagged hemispheres in the data as well as the expected number and the flavour composition of the sample taken from the Monte Carlo simulation plus their statistical errors. The sizeable charm quark contribution is due to D mesons decaying into K 0 S .

K ± selection
Charged kaons were identified using the energy loss measurements in the jet chamber [10].
To effectively reject the contribution of charged pions and protons, only tracks in the central region of the jet chamber, | cos θ | < 0.72, were considered which have at least 130 hits for the momentum measurement and a minimum of 100 hits for the determination of the energy loss. Tracks azimuthally closer than 0.5 • to an anode wire plane of the jet chamber were not accepted to avoid biases from the significantly degraded resolution of both momentum and energy loss measurement close to these planes.
Each track is assigned five dE/dx weights, w h , calculated from the probabilities of the measured energy loss of the track to be consistent with that expected for a particle of type h, where h can be p(p), K ± , π ± , µ ± or e ± . The weight is positive signed if dE/dx(measured)>dE/dx(expected) according to a particular hypothesis, and negative otherwise. K ± were selected by requiring: |w K | > max (|w p |, |w π |, |w µ |, |w e |) |w K | > 0.1 and |w π | < 0.1 This yielded, according to the simulation, a composition of the sample of tagged particles of about 77 % kaons, 13 % pions, and 9.5 % protons. The remainder is due to muons and charged hyperons. It was not attempted to correct the K ± selection for the pion and proton contributions directly.  shows the energy spectrum of the K ± candidates. To enrich the sample in events of primary strange and up quarks, a cut on the minimum scaled energy of x E > 0.5 was applied as shown in Figure 1 (e). Badly reconstructed tracks were rejected by requiring x E < 1.07. Table 1 summarizes the number of tagged hemispheres, the expected number and the primary flavour composition taken from the Monte Carlo simulation.

High x E particle selection
Even though the fragmentation of bottom and charm quarks is harder than that of the light quarks, the cascade decays of the b-and c-hadrons lead to many stable particles which consequently have a significantly softer energy spectrum. Thus the selection of stable charged particles carrying a high energy fraction x E , calculated using the assumption that the particle is a pion, depletes an event sample of primary charm and bottom quarks. Besides the cuts on the scaled energy of 0.7 < x E < 1.07, additional cuts concerned the selection of well-reconstructed tracks which were required to have at least 130 hits used for the reconstruction, a polar angle of | cos θ | < 0.72 and an azimuthal angular distance from the closest anode wire plane of at least 0.5 • . Moreover, the χ 2 rϕ per degree of freedom of the fit of the track to the hits in the plane perpendicular to the beam axis had to be less than 2. This tight cut rejected badly reconstructed tracks, including those containing kinks due to an unrecognized decay in flight.
The energy spectrum of the high x E candidates and the flavour composition depending on the cut on the minimum energy of the high x E particles are shown in Figures 1 (c) and (f). The slightly rising fraction of heavy quark events at very large x E is due to the remaining contribution of badly reconstructed tracks in this energy region. Such tracks occur with the same low probability for each primary flavour.
In Table 1 the number of tagged hemispheres found in the data are given and compared with the expected number from the simulation. It should be pointed out that about 10 % of the high x E hemispheres are also selected with the K ± selection cuts. No rejection was applied, but it was checked that eliminating the 10 % had no significant effect on the result of this analysis. The small statistical correlation was taken into account.

Checks of the flavour composition
A crucial ingredient in this analysis is the knowledge of the flavour composition of the three tagged event samples. The JETSET Monte Carlo simulation using the tuned parameters of [14] predicts for each of the tagged particles in the region of high values of x E a larger differential cross-section (1/N) · (dN/dx E ) than measured in the data. This discrepancy was accounted for by applying a flavour independent reweighting to each of the three tagged event samples in the simulation. The reweighting factors were derived from the ratio of the x E spectra of the tagged particles in the data and the simulation. They were parametrized by a function f (x E ) = c · (1 − x E ) δ whose form is motivated by the Lund symmetric fragmentation function [18]. Each tag type had its own reweighting function, Table 2: Number of doubly tagged events when: (a) both tags belong to the lower x E region, (b) at least one tag belongs to a higher x E region. The upper rows are the measured data. The lower rows present the reweighted simulation data (MC). Also the total number of tags are given for data and the simulation (third column).
but δ was found to be consistent for all three tags. The impact of this reweighting on the flavour fractions is negligible and was considered as an additional systematic uncertainty in Section 5.
The flavour composition and its sensitivity to changes in the fragmentation function can be checked with measured data. To this end the number of events with a single and a double tag have been counted and compared in measured and simulated data following the proposal of [19]. Here a double-tagged event has exactly one tagged particle in each hemisphere, whereas single-tagged events have at least one tag. These counts depend on the sums and products of the flavour fractions. Due to the small number of tagging samples it is not possible to find a unique solution in this analysis. However, the singletagged and double-tagged events can be used to cross-check the flavour fractions obtained from the Monte Carlo simulation using measured data only.
Double-tagged events were studied in different x E ranges of the tagged particles. Each of the three samples was divided into two x E regions. For each sample the lower x E region was chosen so that the impact from the choice of the parameters of the fragmentation function is small. The remaining part corresponds to the higher x E region which is particularly sensitive to the hardness of the fragmentation function. In detail, the lower x E ranges were chosen for the K 0 S and K ± samples to be 0.4 < x E < 0.6, and for the high x E sample to be 0.6 < x E < 0.7. The higher x E regions are 0.6 < x E < 1.07 for both kaon selections and 0.7 < x E < 1.07 for the high x E sample. For both the lower and the higher x E regions the number of double-tagged events were determined. The events were classified according to whether (a) both tags belong to the lower x E region, or (b) at least one tag belongs to a higher x E region. Table 2 lists the double tag counts for data and simulation, normalized to the integrated luminosity of the data after the reweighting. Good agreement between the double tag counts in data and simulation is found for both x E regions, as can be seen from Table 2. This can be quantified using the χ 2 per degree of freedom of the double tag counts which is 3.6/8 (5.3/8) in the higher (lower) x E regions. Since the double tag counts depend on the flavour fractions, the good agreement for the absolute number of tags and for the number of double tags leads to the conclusion that the estimate of the flavour fractions by the JETSET generator is reliable. Different multiplicities for events from primary u, d, or s quarks could imply different fragmentation functions for these quarks. To test the impact of different fragmentation functions another cross-check was conducted which considers a different reweighting factor f (x E ) = c · (1 − x E ) δ for one of the three light quarks. This corresponds to the case of a different fragmentation function for one quark flavour with respect to the other two. The flavour dependence was modelled by allowing in turn for a different value of δ for one of the three light quarks while maintaining as a constraint in the fit of the δ parameters the agreement of both the number of single tag counts and the x E spectra in data and simulation as before.
When this reweighting was applied to the simulation, the light flavour fractions of the two kaon selections changed by less than 0.08, while for the high x E selection the fractions changed by up to a factor of two. The comparison of the double tag counts after the flavour dependent reweighting revealed, however, a worse description of the data by the simulation than in the case of the flavour independent reweighting. This can be inferred from the χ 2 per degree of freedom for the number of double tag counts in the whole x E regions considered which increased to 27.7/7, 13.7/7, or 36.5/7, respectively for each sample, when allowing for a different value of δ and, hence, a different fragmentation function for up, down, or strange quarks. These large values of χ 2 per degree of freedom strongly disfavour a flavour dependent reweighting of the fragmentation functions and, therefore, the substantial changes of the flavour fractions associated with this reweighting. This gives further evidence that the flavour fractions derived from the simulation can be considered to be an accurate estimate of the flavour fractions in the data. We, therefore, use for our main result the flavour fractions obtained from JETSET without application of the reweighting procedure. The small effect on the flavour fractions due to the reweighting will yield a systematic uncertainty which is discussed in Section 5. 4 Determination of charged particle multiplicities 4

.1 Determination of flavour dependent multiplicities
The flavour dependent multiplicities are defined as the average number of charged particles with life-times τ > 300 ps emerging from a decay of a Z into one of the three light quark flavours, u, d, or s, as for the inclusive multiplicities in [20]. To obtain these multiplicities from the charged particle tracks measured for each of the three selections, several corrections must be applied. In particular, the finite detector resolution and acceptance, and biases due to the tag selections must be accounted for.
The requirement of a high energy of the tagging particles in all three selections affects the multiplicity in its hemisphere due to a reduced phase space for the production of other particles. In the opposite hemisphere the effect is much smaller. This can be seen in Figure 2 where the mean charged particle multiplicities per hemisphere in the three selections are plotted in bins of the tagging particle's energy. To keep the bias of the multiplicity due to the tagging small we determine the multiplicities only in the hemispheres opposite to the tagged particles.
The remaining small dependence of the multiplicity in the opposite hemisphere on the energy of the tagging particle is due to the impact of the high energy of this particle. The higher the energy required of the leading particle the more the phase space is reduced, which is available for particle production, and, therefore, the mean charged multiplicity in the opposite hemisphere. It should also be noted from Figure 2 that the Monte Carlo simulation satisfactorily describes the effects in the opposite hemisphere due to the tag, although it tends to slightly overestimate the multiplicity in particular if the energy of the tagging particle is very high. This is due to a slight underestimation of the rate of badly reconstructed tracks in the simulation which particularly affects the hemisphere of the tagging particle if the measured energy fraction, x E , of such a track is close to or above the physical limit of 1. Since the differential cross-section is tiny for leading particles with x E close to 1, the effect of events selected due to badly reconstructed tracks on the mean charged multiplicity in the opposite hemisphere is negligible when considering the full range of x E used for the tagging. Any remaining discrepancies between data and Monte Carlo simulation in the opposite hemisphere are taken as a source of systematic error, which is treated in Section 5.   Table 3 lists the averages of the distributions for the hemisphere opposite to the tag. In the K 0 S sample good agreement between data and the Monte Carlo is found, but in the events selected by a highly energetic K ± or high x E tag the simulation predicts a slightly larger multiplicity.

Unfolding procedure
The three tag samples are composed of events with all five primary quark flavours accessible at LEP I. A statistical unfolding procedure was applied to obtain the mean charged particle multiplicities for each light quark flavour. The procedure simultaneously corrects for effects due to the selection cuts and for detector effects.   Table 4: Mean generated charged multiplicities n q, MC had. as obtained from the Monte Carlo generator JETSET 7.4. Only charged particles with life-times τ > 300 ps were considered.
The measured mean charged multiplicities, n tag , in the hemisphere opposite to the tagged particle are expressed as Here f tag q is the fraction of events with primary quark flavour q in a sample denoted by "tag". It is listed in Table 1. The term in square brackets was also taken from the Monte Carlo simulation. It is the ratio of the mean charged multiplicity of the opposite hemisphere for events with primary quark flavour q in the "tag" sample, n tag q, MC , listed in Table 3, and the number for both hemispheres obtained without detector simulation, n q, MC had. , which is given in Table 4. This ratio corrects (i) for the biases due to the tagging cuts, (ii) for the restriction to the hemisphere opposite to the tag particle when determining the charged multiplicity, and (iii) for effects from the selection of multihadronic events, the finite resolution and acceptance of the detector. On average the correction factor in square brackets is 1/(2 × 1.007) for all flavours varying at most by ±3.7 %, and where the factor of 2 is due to considering a single hemisphere only for the determination of the multiplicities denoted by "tag".
The flavour dependent mean charged multiplicities, n q , can be determined by solving the equation system (1). The values of n c and n b were taken from our previous determination [21]: where the errors are statistical and systematic, respectively. The solution of the equation system (1) using the measured multiplicities in the hemisphere opposite to the tagging particles from Table 3  and using the branching ratios of Z → uu, dd, ss from [22] one obtains n uds = 19.90 ± 0.09 (stat.). This is in agreement within systematic errors with the results of [21,23] which were obtained using a different method. The mean charged multiplicities n u and n d are almost fully anti-correlated, since the statistical separation power for u and d quarks is due to the small but different fractions of these quarks in the K 0 S and the K ± samples.

Test of the unfolding procedure
It was checked that the unfolding procedure described in Section 4.2 is able to determine the flavour dependent multiplicities. For this test the Monte Carlo data sample was split into two disjoint parts. The smaller, comprising about the same statistics as the measured data, was used in place of the measured data which were to be corrected using the larger part of the Monte Carlo events. The multiplicities n u = 20.05 ± 0.63 (stat.), n d = 19.98 ± 0.81 (stat.), and n s = 19.89 ± 0.16 (stat.) are in excellent agreement with the mean multiplicities obtained from the Monte Carlo generator (cp. Table 4).

Systematic study of the particle multiplicities
In addition to the consistency check based on the double tag rates in Section 3.4 and the test of the unfolding procedure in Section 4.3, many sources of systematic uncertainties were investigated which may be subdivided into five groups, namely (i) fluctuations due to limited Monte Carlo statistics, (ii) the experimental precision of the multiplicity in charm and bottom events, (iii) impacts of detector effects and event selection and the dependence of the mean charged multiplicity on the choice of the hadronization model, (iv) variations of tagging cuts, and finally, (v) impacts on the flavour fractions, f tag q , due to variations of hadronization parameters of the Monte Carlo generator, the matching of the Monte Carlo x E spectra and hemisphere correlations to the data, and changes of the hadronization model. Apart from the items (i) and (ii), whose error contributions were estimated by error propagation, the general procedure to derive the uncertainty associated with the source was to repeat the complete analysis with one of the cuts or parameters varied. Any deviation from the result found for the standard set of cuts and parameters was interpreted as a systematic uncertainty.
The respective error contributions of the uncertainties derived for the groups (i) to (v) were quadratically added to estimate the total systematic uncertainties. The individual error contributions are listed in Table 5. In the following sections the sources (iii) to (v) of systematic uncertainties will be discussed in more detail.

Uncertainties from detector simulation and event selection
The detector simulation was needed to correct for the fraction of multiplicity that was not recorded by the detector either owing to limited acceptance, biases of the event selection, or interactions in the detector such as δ-electrons or hadronic interactions, all of which have been studied in great detail in [20]. In total an error of about 1 % due to this source of systematic uncertainty was estimated in [20]. This was also adopted for this measurement.

Uncertainties from tagging
The measured multiplicity depends on the energy of the tagged particle. Although this effect is reduced by measuring the multiplicity in the hemisphere opposite to the tag residual effects have to be taken into account.
• Variations of the x E cuts in steps of 0.05 by as much as ±0.1 around the default cut value were done. The largest up-and downward excursions yielded the error estimates for the multiplicities.
• Also a more stringent cut on the maximum accepted momentum, x E < 0.9, for all tags did not change the measured multiplicities by more than 0.4 %.
For each of the three tag types the impact of the choice of the most important selection cuts was studied.
• For the K 0 S the allowed mass window around the world average mass value [17] was altered from 120 MeV within the range of 80 to 200 MeV.
Since the K ± selection relies very much on the capability to precisely measure the energy loss in the jet chamber, the relevant cuts were investigated.
• As the precision of the measurement of dE/dx degrades towards the very forward and backward regions, the | cos θ K ± | < 0.72 requirement was varied in small steps between 0.6 and 0.8, taking the r.m.s. to estimate the contribution to the systematic uncertainty.
• The minimum number of both jet chamber hits and dE/dx hits was relaxed from 130 to 80.
• The required number of dE/dx hits was increased from 100 to 110.
• From the various modifications to the energy loss weight cuts, namely choosing a tighter cut on the kaon dE/dx weight of 0.5 instead of 0.1, adding an additional rejection cut on the proton dE/dx weight of 0.4, or considering a better pion rejection by requiring its dE/dx weight to be larger than −0.06, only the first one contributes a significant systematic error.
The high x E selection critically depends on the ability of the jet chamber to reliably measure high momentum particles.
• As for K ± , the allowed range of | cos θ high x E | was varied from less than 0.72 to either less than 0.6 or 0.8.
• A variation of the minimum number of hits required for a high x E particle in the jet chamber by ±10 around the standard value of 130 had only a negligible effect on all flavours.
• Tightening the cut on the χ 2 rϕ of the track fit from 2 to 1.5 affected the multiplicity by less than 0.5 %.

Uncertainties due to impacts on flavour fractions
The Monte Carlo event generator requires several parameters to be adjusted for a proper description of the measured data. Changing these hadronization parameters affects the flavour composition, f tag q , and, therefore, tests the sensitivity of the result of this analysis on the particular choice of these parameters. In the following we consider those parameters which are expected to have an impact on the flavour fractions. These were varied about their tuned values within the intervals quoted in [14] apart from from the b and γ s parameters where larger variations were investigated, in particular the range for γ s given in [24]. Most of the parameter variations did not significantly affect the flavour dependent multiplicities, n q, MC had. , or the flavour fractions, f tag q , of Eq. (1). We therefore mention only the significant changes of the flavour fractions due to the parameter variations.
• The amount of gluon radiation was modified by varying the Λ LLA parameter and the cut-off of the parton shower, Q 0 .
• During the hadronization step a hadron receives extra transverse momentum whose size is controlled by the parameter σ q which was changed from its tuned value.
• The fraction of energy and momentum transferred from a heavy c or b quark to the related hadron also has an impact on the background contributions from these heavy quarks. The relevant ε parameters of the Peterson et al. [25] fragmentation function were varied. Only the change of ε c , except for the charm quark fraction itself, had a noticeable effect on the strange quark fraction of about −0.01 to +0.005 for the two kaon selections.
• To account for a remaining uncertainty from the hardness of the light quark fragmentation functions despite the impact of the reweighting procedure, the parameter b of the Lund symmetric fragmentation function in JETSET was varied in several steps in the range of 0.26 to 0.92. This changed in particular the strange quark fraction in the range of −0.06 to +0.04 for the K 0 S and in the range of −0.04 to +0.03 for the K ± selection. The u and d fractions varied owing to this variation of b by ±0.012 at most in all three selections.
• Since the statistical separation between u and d quarks is based on the K 0 S and K ± production, the relative amount of ss quark pairs created from the vacuum, γ s , was changed in steps in the range of 0.21 to 0.41 [24]. This led to changes of the light quark fractions in the K ± (K 0 S ) selection ranging from −0.05, −0.03, +0.06 (−0.02, −0.03, +0.04) to +0.06, +0.02, −0.05 (+0.01, +0.02, −0.04), respectively, for u, d, s quarks. For the high x E selection only changes of the u and s quark fractions were observed at a level of ±0.01.
• For baryons with strangeness an extra suppression factor (us/ud)/(s/d) exists in JETSET whose variation resulted in only negligible changes of both the flavour composition and the flavour dependent multiplicities.
• The production of vector mesons was considered. Although the down and strange quark fractions of the high x E selection changed by about ±0.02, the relative portion V u,d of vector mesons from u and d quarks was of minor importance for the flavour dependent multiplicities.
• The fraction of vector mesons from s quarks, V s , substantially affects the yields of leading K 0 S and K ± mesons. The variation of this parameter mostly changed the fraction from up and down quarks in the K ± selection by about ±0.01, and at the same level also that from down and strange quarks of the high x E selection. Since the statistical u-d quark separation strongly relies on the u and d quark fractions in the two kaon selections, a change of this parameter led to a large error contribution as can be seen from Table 5.
• Similarly, the decay products of tensor mesons have a softer x E spectrum. However, the yields of these are not well known. Thus, the impact was estimated by analyzing about 3.5 million fully simulated events from the JETSET generator tuned to describe the data, but without the production of tensor mesons. This yielded a substantial reduction of the strange quark fraction by −0.025 in the K 0 S , −0.035 in the K ± , and −0.065 in the high x E selections. The up and down quark fractions in the K 0 S and K ± selections changed by −0.01 and −0.005, whereas the d quark fraction in the high x E selection increased by about +0.03. Thus the d and s quark fractions in the high x E selection are about equal if tensor mesons do not contribute. The difference in the results with simulated tensor mesons is assigned as a systematic error.
• The uncertainty of the inclusive branching ratios of charm and bottom to K 0 S and K ± was considered. From the data in [17] the relative errors on the B(c → K 0 S X, K ± X) branching ratios have been determined to be ±11.5 %, ±13.9 %, and on B(b → K 0 S X, K ± X) to be ±21.6 %, ±10.0 %, respectively. Their impact on the flavour dependent multiplicities was assessed by varying the charm and bottom contributions in the K 0 S and K ± tagged event samples according to these percentages.
• To account for the poor simulation of the x E spectra of the tagged particles in the tuned Monte Carlo the reweighting functions of Section 3.4, f (x E ) = c · (1 − x E ) δ , were employed. Applying these specific reweighting functions to all three selections resulted in a small reduction of the strange quark fraction by 1 to 2 % while the heavy quark contribution is increased accordingly. Moreover, it leads to an anticorrelated change of the order of 1 to 3 % of the u and d quark fractions but in opposite directions for the K 0 S and the K ± event samples. The final multiplicities of u and d quarks change by up to 2 %.
• Finally, the small discrepancies between simulation and data noticed in the correlations of the opposite hemispheres shown in Figure 2 were considered. Depending on the x E of the tagged particle the Monte Carlo events were reweighted to match the data. This had only a negligible effect on both the flavour fractions and the flavour dependent multiplicities changing them by less than 0.1 %.
In addition to the JETSET Monte Carlo event generator also the HERWIG program [26] was examined. While HERWIG yields a good description of the shape of the x E spectra of the tagged particles, it does not reproduce the absolute rates. For instance, HERWIG predicts a much larger contribution to the high x E sample from d quarks and less from u quarks due to a larger number of protons produced in HERWIG. Since the HERWIG expectation of the total proton yield is not in agreement with the measurements [17,27], one might suspect that the flavour fractions obtained from this Monte Carlo generator are not reliable. Thus, no systematic error contribution from using HER-WIG is quoted.
6 Cross-check with flavour fractions derived from double tag counts The flavour fractions, f tag q , can be derived from data only using the number of single and double tag counts since they are related to the sum and products of flavour fractions [19]. The non-linearity of the equation system requires for its solution more than the six relations which are given by the numbers of single and double tags considered in this analysis even when the heavy flavour contributions are kept fixed.
Such an investigation has been done in [8] using the single and double tag counts of charged pions, kaons, protons, K 0 S and Λ for x E > 0.4, 0.5, and 0.6. This analysis yielded the flavour dependent production rates of these hadrons where only the SU(2) isospin flavour symmetry and the flavour independence of the strong interaction were assumed for the production of charged pions, and charged and neutral kaons.
To compare with the results from [8], the flavour fractions, f tag q , were calculated from the quoted production rates neglecting correlations for this cross-check. The calculation took into account that these rates had been corrected for mis-identification of the tagging particle. The high x E sample was approximated by adding the production rates of the charged particles, i.e. pions, kaons, and protons. The comparison of the flavour fractions at an x E -cut of 0.4 for K 0 S , 0.5 for K ± and 0.6 in the case of the high x E sample revealed some substantial differences with respect to the JETSET prediction. The K 0 S fraction from up quarks is 0.4 times and that from down quarks 1.7 times the value obtained from the Monte Carlo simulation used for this analysis. The up quark contribution to the K ± sample is 1.23 times higher, while that of down quarks is about 0.75 times that obtained from the simulation. The fraction of high x E particles originating from up quarks is about 10 % less than expected from the simulation. The strange quark fractions agree within ±0.02 for each tagged event sample. These changes in the up and down quark flavour fractions are considerably larger than what was found from the variation of the JETSET parameters in Section 5.3.
Although the analysis in [8] explicitly assumed flavour independence of the strong coupling and SU(2) isospin symmetry, the flavour fractions obtained from [8] were used for the unfolding procedure described in Section 4.2 as a cross-check of the results. This yielded changes of the flavour dependent multiplicities of ∆ n u = +0.60, ∆ n d = −1.41, and ∆ n s = +0.32. Since these changes are essentially covered by the systematic errors discussed in Section 5 they are not considered as an additional systematic uncertainty.

Results
Adding the individual contributions to the systematic uncertainty from Monte Carlo statistics, the mean charged multiplicity in c and b quark events, the correction of detector acceptance and resolution, the selection cuts, and the parameters of the hadronization model in quadrature, the total systematic error of the flavour dependent mean charged multiplicities is about ±7 % for up and down, and ±2 % for strange quark events. The final results are All errors apart from selection cuts, detector effects, and hadronization models were obtained by error propagation and, therefore, obey the statistical correlation as quoted in Section 4.2. The correlation coefficients for the remaining two groups of systematic uncertainties were determined separately for each individual systematic variation.
These flavour dependent mean charged multiplicities agree within the total errors, even though deviations may be expected due to decays of hadrons. For instance, the decay of the φ(1020) to charged kaons is enhanced over the decay to neutral kaons. Thus, the results from the JETSET generator displayed in Table 4 show a 1% spread between the charged multiplicities from up and strange quark events.
Some of the flavour dependent mean charged multiplicities are strongly correlated which must be considered when comparing the results. Taking these statistical and systematic correlations into account the ratios of the multiplicities are which are all close to unity.

Summary and conclusions
A measurement of the charged multiplicity in events of u, d, and s quarks from e + e − annihilation at the Z pole is presented. Exploiting the leading particle effect to identify the primary quark flavour, samples of events differently enriched in u, d, and s were selected by tagging highly energetic K 0 S , K ± and charged particles. The multiplicities per light quark flavour were obtained from a statistical unfolding of the multiplicities measured in the hemispheres opposite to the tagged particles. Hemisphere correlations were observed, but were found to be small and well-described by the Monte Carlo simulation. The final corrected multiplicities were determined to be The ratios of pairs of these multiplicities which take the correlations into account are all close to unity. This agrees with the expectation of QCD that the charged particle multiplicities are flavour independent for the light up, down, and strange quarks apart from small effects due to particle decays.