Particle multiplicity of unbiased gluon jets from e+e- three-jet events

The charged particle multiplicities of two- and three-jet events from the reaction e+e- ->Z0 ->hadrons are measured for Z0 decays to light quark (uds) flavors. Using recent theoretical expressions to account for biases from event selection, results corresponding to unbiased gluon jets are extracted over a range of jet energies from about 11 to 30 GeV. We find consistency between these results and direct measurements of unbiased gluon jet multiplicity from upsilon and Z0 decays. The unbiased gluon jet data including the direct measurements are compared to corresponding results for quark jets. We perform fits based on analytic expressions for particle multiplicity in jets to determine the ratio r = Ng/Nq of multiplicities between gluon and quark jets as a function of energy. We also determine the ratio of slopes, r(1) = (dNg/dy)/(dNq/dy), and of curvatures, r(2) = (d2Ng/dy2)/(d2Nq/dy2), where y specifies the energy scale. At 30 GeV, we find r = 1.422 +/- 0.051, r(1) = 1.761 +/- 0.071 and r(2) = 1.98 +/- 0.13, where the uncertainties are the statistical and systematic terms added in quadrature. These results are in general agreement with theoretical predictions. In addition, we use the measurements of the energy dependence of Ng and Nq to determine an effective value of the ratio of QCD color factors, CA/CF. Our result, CA/CF = 2.23 +/- 0.14 (total), is consistent with the QCD value of 2.25.


Introduction
The mean charged particle multiplicity of a gluon jet has often been measured in the annihilation of an electron and positron to hadrons, e + e − → hadrons.The usual method (see for example [1]- [4]) is to select three-jet quark-antiquark-gluon qqg final states for which the events and individual jets are defined using a jet algorithm such as the Durham [5] or Luclus [6] jet finder.The particle multiplicity of a jet determined with this technique is found to depend on which algorithm is employed.Therefore, these jets and the associated qqg events are called "biased." In contrast, theoretical calculations usually define gluon jet multiplicity inclusively, by the particles in hemispheres of gluon-gluon (gg) systems in an overall color singlet.Quark jets are defined analogously as hemispheres of quark-antiquark (qq) systems.The hemisphere definition of jets yields results which are independent of a jet finder.Therefore, these jets are called "unbiased."Unbiased gluon jet multiplicity has so far been measured only in Υ [7,8] and Z 0 [9]- [11] decays, corresponding to jet energies of about 5 and 40 GeV, respectively.In contrast, unbiased quark jet multiplicity is easy to measure and has been determined at many scales (for a recent review, see [12]).
It is of interest to measure unbiased gluon jet multiplicity at other scales.Such measurements would allow a test of recent predictions [13,14] from Quantum Chromodynamics (QCD) for the scale dependence of the multiplicity ratio r between gluon and quark jets and for the related ratios of slopes r (1) and of curvatures r (2) : r (1) ≡ dN g /dy dN q /dy , where N g and N q are the mean particle multiplicities of gluon and quark jets, y = ln (Q/Λ), Q is the jet energy and Λ is the QCD scale parameter.Recently, a method to extract unbiased gluon jet multiplicity from biased e + e − → qqg events was proposed [15,16], extending earlier formalism [17].This method provides an indirect means to determine the particle multiplicity of gluon jets, in a manner which corresponds to the hemisphere definition of jets.By combining measurements of N g found from this method with the unbiased measurements for N q , the ratios r, r (1) and r (2) can be determined at a variety of scales and used to test the corresponding QCD results in a quantitative manner.
The purpose of the present study is to test this method [15,16] to obtain unbiased measurements of N g and to apply the results to determine the ratios, eqs.( 1)- (3).The data were collected with the OPAL detector operating at the LEP accelerator at CERN.We determine the ratios (1)-(3) as a function of scale and compare the results to recent QCD predictions [13,14].Previously published measurements of the scale dependence of r [4,18,19] and r (1)  [18,19] were based on biased jet samples, whereas experimental results for the ratio of curvatures r (2) have not previously been presented.We also use the results on the scale dependence of unbiased gluon and quark jet multiplicities to derive a measurement of the ratio of QCD color factors C A /C F .

Theoretical Framework
Analytic expressions for the mean particle multiplicity of e + e − three-jet events, valid to the nextto-leading order of perturbation theory (NLO, also called MLLA [20]), were recently presented in [16]: These results, along with the others presented in this section, are based on massless quarks.
The reason for the two different expressions, eqs.( 4) and (5), is that there is an ambiguity in the definition of the gluon jet transverse momentum when the gluon radiation is hard.k ⊥,Lu and k ⊥,Le are the transverse momenta of the gluon with respect to the quark-antiquark system using the definition of either the Lund (k ⊥,Lu ) [15] or Leningrad (k ⊥,Le ) [17] groups, see [16] and references therein for further discussion.Note that due to QCD coherence, the multiplicity of a gluon jet produced in e + e − annihilations depends on its transverse momentum and not its energy (see for example [17]).N qqg is the particle multiplicity of a three-jet event selected using a jet algorithm (see below).N qq and N gg are the multiplicities of two-jet qq and gg systems, given by twice N q and N g , respectively.The scales L, L qq , k ⊥,Lu and k ⊥,Le are defined as follows: p ⊥,Le = s qg s qg s qq , with s = E 2 c.m. , s ig = p i •p g (i = q, q), and p q , p q and p g the 4-momenta of the q, q and g.L specifies the e + e − c.m. energy (E c.m. ) and L qq the energy of the qq system in the qq rest frame.
The multiplicity of the gg system in this formalism, N gg , depends only on a single scale: k ⊥,Lu in eq. ( 4) or k ⊥,Le in eq. ( 5).This dependence on a single scale is a statement that N gg is unbiased, i.e.N gg (k ⊥ ) in eq. ( 4) or ( 5) is equivalent to the inclusive multiplicity of a gg event from a color singlet source produced at the same scale k ⊥ , to NLO accuracy.In contrast, the multiplicity of the qq system, N qq , depends on two scales: L and k ⊥,Lu in eq. ( 4) or L qq and k ⊥,Lu in eq. ( 5).This dependence on two scales implies that the multiplicity of the qq system is biased, depending on a hard scale L or L qq for the production of the quark and antiquark jets and on a cutoff k ⊥ ≈ k ⊥,Lu below which the gluon jet is not resolved.This accounting for the bias in quark jet multiplicity due to the jet finder criteria used to select the qqg events is the difference between the theoretical results used here [15] and those presented previously [17].
The above expressions are valid for jet algorithms which employ a transverse momentum cutoff k ⊥ to resolve "two-jet" qq from "three-jet" qqg events.Examples of such algorithms are the Durham [5], Cambridge [21] and Luclus [6] jet finders.The formalism is not expected to be valid for the Jade [22] or cone jet finders [23], which use invariant masses or angles rather than transverse momenta to separate the two-and three-jet event classes.
According to the prescription in [16], the resolution scale k ⊥ should be adjusted separately for each event so exactly three jets are reconstructed.In contrast, a fixed value of the resolution scale results in truncation of higher order radiation, biasing the gluon jet's properties.Note that any radiation ("sub-jet") emitted within the quark or gluon jet must necessarily have a smaller transverse momentum than the gluon jet itself: otherwise the role of the "gluon jet" and "sub-jet" would be reversed.Thus the transverse momentum of the gluon jet, k ⊥,gluon , defines an effective cutoff for sub-jet radiation, i.e.
In the formalism of [15], the biased quark jet multiplicities N qq (J, k ⊥,Lu ) in eqs.( 4) and ( 5) (with J representing either L or L qq ) can be obtained from the corresponding unbiased multiplicities N qq (J ′ ), which depend only on a single scale J ′ , through the relation: where J ′ = k ⊥,Lu + c q with c q = 3/2.The unbiased term N qq (K) is equivalent to the mean multiplicity of inclusive e + e − → hadrons events as a function of K = E c.m. .
Finally, a relationship between the scale evolution of unbiased gluon and quark jet multiplicities is given [15]: where

Detector and data sample
The OPAL detector is described in detail elsewhere [24,25].The tracking system consists of a silicon microvertex detector, an inner vertex chamber, a large volume jet chamber, and specialized chambers at the outer radius of the jet chamber which improve the measurements in the z-direction. 1The tracking system covers the region | cos θ| < 0.98 and is enclosed by a solenoidal magnet coil with an axial field of 0.435 T. Electromagnetic energy is measured by a lead-glass calorimeter located outside the magnet coil, which also covers | cos θ| < 0.98.
The present analysis is based on a sample of about 2 283 000 hadronic annihilation events, corresponding to the OPAL sample collected within 0.3 GeV of the Z 0 peak that includes readout of both the r-φ and z coordinates of the silicon strip microvertex detector [25].The procedures for identifying hadronic annihilation events are described in [26].
We employ charged tracks measured in the tracking chambers and clusters measured in the electromagnetic calorimeter.Charged tracks are required to have at least 20 measured points (of 159 possible) in the jet chamber, or at least 50% of the points expected based on the track's polar angle, whichever is smaller.In addition, the tracks are required to have a momentum in the direction perpendicular to the beam axis greater than 0.05 GeV/c, to lie in the region | cos θ| < 0.96, to point to the origin to within 5 cm in the r-φ plane and 30 cm in the z direction, and to yield a reasonable χ2 per degree-of-freedom (d.o.f.) for the track fit in the r-φ plane.Clusters are required to have an energy greater than 0.10 GeV if they are in the barrel section of the detector (| cos θ| < 0.82) or 0.25 GeV if they are in the endcap section (0.82 < | cos θ| < 0.98).A matching algorithm is employed to reduce double counting of energy in cases where charged tracks point towards electromagnetic clusters.Specifically, the expected calorimeter energy of the associated tracks is subtracted from the cluster energy.If the energy of a cluster is smaller than that expected for the associated tracks, the cluster is not used.Each accepted track and cluster is considered to be a particle.Tracks are assigned the pion mass.Clusters are assigned zero mass since they originate mostly from photons.
To eliminate residual background and events in which a significant number of particles is lost near the beam direction, the number of accepted charged tracks in an event is required to be at least five and the thrust axis [27] of the event, calculated using the particles, is required to satisfy | cos(θ thrust )| < 0.90, where θ thrust is the angle between the thrust and beam axes.The residual background to the sample of hadronic events from all sources is estimated to be less than 1% [26].

Selection of uds events
The theoretical formalism we employ is based on massless quarks, as stated in Sect. 2. Therefore, beginning with the data sample described in Sect.3, we select light quark (u, d and s) events for our study.
The uds event tagging is based on the signed impact parameter values of charged tracks with respect to the primary event vertex, d sign .The magnitude of d sign is given by the distance of closest approach of the track to the event vertex in the r-φ plane.The sign of d sign is determined as follows.Jets are reconstructed in each event using the cone jet finder [23]. 2 For each track assigned to a jet by the jet finder, hemispheres are defined by the plane perpendicular to the corresponding jet axis.The sign of d sign is determined using the crossing point in the r-φ plane between the track and jet axis to which it is assigned.d sign is positive if the crossing point is in the forward hemisphere and negative otherwise.The forward hemisphere is the hemisphere which contains the jet.Because of the relatively long lifetime of hadrons containing c and b quarks, the distribution of d sign is strongly skewed toward positive values for c and b events but not for uds events.
Note that the jet finding procedure described in the previous paragraph is employed solely for the purpose of identifying uds events and not for the analysis of gluon jet multiplicity presented below in Sects.5-6.
Charged tracks are selected for the uds tagging procedure if they have r-φ coordinate information from at least one of the two silicon detector layers, a momentum of 1 GeV/c or larger, and a maximum distance of closest approach to the primary event vertex in the r-φ plane of 0.3 cm with a maximum uncertainty on this quantity of 0.1 cm.In addition they must be assigned to a jet using the cone jet finding procedure mentioned above.If no track in an event satisfies these requirements (1.5% of the events), the event is eliminated.The number of tracks which meet these requirements and have d sign /σ d sign > 3.0 is determined, where σ d sign is the uncertainty associated with d sign .An event is tagged as a uds event if this number is zero.In total 1 109 017 events are tagged.The number of events in our final event sample and the corresponding uds purity are presented in Sect. 5.

Analysis method
Three-jet events are defined using a jet algorithm.For our standard analysis we employ the Durham jet finder [5].As a systematic check, we use the Cambridge [21] and Luclus [6] jet finders (see Sect. 7).We choose these algorithms because they employ a cutoff in transverse momentum to separate two-and three-jet events, in correspondence to the theoretical formalism of [15,16].In particular, this means that eqs.(4), ( 5), ( 13) and ( 14) are expected to be applicable for these jet algorithms.
The resolution scale of the Durham jet finder, y cut , is adjusted separately for each tagged uds event so exactly three jets are reconstructed.For the Cambridge jet finder, the resolution scale is again y cut .For Luclus, the corresponding parameter is d join .We apply this procedure rather than use a fixed resolution scale to avoid introducing a bias in the gluon jet properties as discussed in Sect. 2. The jets are assigned energies based on the angles between them, assuming massless kinematics (see for example [28]).The jets are ordered from 1 to 3 such that jet 1 has the highest energy.Events are retained if the angles between the highest energy jet and the other two are the same to within 3 • , the so-called "Y events." 3 A schematic diagram of a Y event is shown in Fig. 1.For Y events, the scales L qq , k ⊥,Lu and k ⊥,Le (eqs.( 7)-( 9)) depend only on E c.m. and one inter-jet angle, conveniently chosen to be θ 1 (see Fig. 1).For 35 • ≤ θ 1 ≤ 120 • , the range of θ 1 we employ for our gluon jet analysis (Sect.6.2), 22 365 events are selected: this is our final event sample.We require θ 1 ≥ 35 • in order that the two lower energy jets be clearly separated from each other.Note that since θ 2 and θ 3 (see Fig. 1) are not exactly equal, the energies of jets 2 and 3 are not exactly the same.
The uds event purity of this sample is estimated using the Jetset Monte Carlo multihadronic 3 Y events were first studied in [28,29]; the overall charged multiplicity of Y events was first presented in [18].
Figure 1: Schematic representation of a three-jet quark-antiquark-gluon qqg event with a Y event topology [28,29], in which the angle between the highest energy jet and each of the two lower energy jets is about the same.The jets are ordered in energy such that jet 1 has the highest energy.The angle θ 1 opposite the highest energy jet is used to specify the event topology.In our analysis, the gluon jet is assumed to be jet 3.
event generator [30] including initial state photon radiation, detector simulation [31], and the same analysis procedures as are applied to the data.We use a sample of about 6 000 000 Monte Carlo events generated with version 7.4 of the program and the parameters given in [32].The estimated uds purity is found to be (78.5 ± 0.2 (stat.))%.The Monte Carlo predicts that 70% of the background events are c events and 30% are b events.The ratio of the number of events in the final sample to the number in the initial multihadronic Z 0 decay sample is (0.980 ± 0.007 (stat.))× 10 −2 for the data and (1.036 ± 0.004 (stat.))× 10 −2 for the Monte Carlo.
For simplicity, we identify the gluon jet by assuming it is the lowest energy jet in an event, i.e. jet 3.For the majority of Y events, θ 1 << θ 2 , θ 3 .For these events, jet 1 has an overwhelming probability (typically about 97%) to be a quark (or antiquark) jet, while jets 2 and 3 have about an equal probability to arise from a quark or gluon.Since θ 2 ≈ θ 3 (see Fig. 1), the scales L qq , k ⊥,Lu and k ⊥,Le are almost unchanged if the gluon is actually jet 2 rather than jet 3. Similarly, as the energy of jet 3 increases so that θ 1 ≈ θ 2 ≈ θ 3 , all three jets have about the same probability to be the gluon jet and the determination of scales is again insensitive to whether the gluon jet assignment is correct or not.This relative insensitivity to gluon jet misidentification is our motivation for employing Y events.
The scales L qq , k ⊥,Lu and k ⊥,Le and the three-jet event charged particle multiplicity N ch.qqg are measured as a function of θ 1 .L qq , k ⊥,Lu and k ⊥,Le are determined using all accepted charged and neutral particles.For the QCD scale parameter Λ, used in the definition of these last three variables (also for the definition of L, see eq. ( 6)), we choose a value of 0.20 GeV.We verified our results are not sensitive to this choice.
The measured distributions of N ch.qqg , L qq , k ⊥,Lu and k ⊥,Le are corrected for experimental acceptance and resolution, the effects of initial state photon radiation, and misidentification of uds events, using bin-by-bin multiplicative factors derived from Jetset.The corrections for L qq , k ⊥,Lu and k ⊥,Le also account for gluon jet misidentification, i.e. for the contributions of events in which the gluon jet is not jet 3. The data are corrected to the hadron level.The hadron level does not include detector response or initial state radiation and treats all charged and neutral particles with lifetimes greater than 3 × 10 −10 s as stable: hence charged particles from the decays of K 0 S and weakly decaying hyperons are included in the corrected distributions.The correction factors are determined as described in [33].The corrections are typically 1.08  1: Measurements of the mean charged particle multiplicity N ch.qqg of three-jet uds flavor Y events from Z 0 decays, as a function of the angle θ 1 between the two lowest energy jets.The jets are defined using the Durham jet-finder.The results for the quark jet scale √ s qq and the gluon jet scales p ⊥,Lu and p ⊥,Le are also given.The uncertainties are statistical only.These data are shown in Fig. 2.
for N ch.qqg , 0.98 for L qq , 1.03 for k ⊥,Lu , and 1.05 for k ⊥,Le , and do not vary much with θ 1 .There is good agreement between the distributions of the uncorrected data and the Monte Carlo sample when the latter includes initial-state radiation and detector simulation, with the same analysis procedures applied to the simulated events as to the data.
Corrected distributions of N ch.qqg , √ s qq , p ⊥,Lu and p ⊥,Le are shown in Fig. 2. The last three variables are obtained from the corrected distributions of L qq , k ⊥,Lu and k ⊥,Le using eqs.( 7)- (9).We prefer to present √ s qq , p ⊥,Lu and p ⊥,Le rather than L qq , k ⊥,Lu and k ⊥,Le because they provide a clearer illustration of the physical scales in our analysis.These data are listed in Table 1.Besides the results from our standard analysis, the results found using the Cambridge and Luclus jet finders are shown in Fig. 2. For about 2% of the events, the resolution parameter y cut of the Cambridge algorithm could not be adjusted to yield three reconstructed jets, see [34] for a discussion of this feature of the algorithm.These events were rejected for the Cambridge  The mean charged particle multiplicity of three-jet uds flavor Y events from Z 0 decays, selected using the Durham [5], Cambridge [21] and Luclus [6] jet finders, as a function of the opening angle θ 1 ; (b) the corresponding scales √ s qq , p ⊥,Lu and p ⊥,Le , see eqs. ( 10)- (12).The data have been corrected for initial state photon radiation, detector response, misidentification of uds events, and (for part (b)) misidentification of gluon jets.The uncertainties (too small to be visible) are statistical only.jet finder based study.
We invert eqs.( 4) and ( 5) to obtain expressions for the unbiased gluon jet charged particle multiplicities N ch.
gg (k ⊥ ): where N ch.qqg , L qq , k ⊥,Lu and k ⊥,Le are determined from the Y events as described above.
To find the N ch.qq (L, k ⊥,Lu ) terms in eq. ( 15), we employ two methods.First, for the standard analysis, we perform a direct measurement.Specifically we determine the particle multiplicity of two-jet uds events from Z 0 decays (i.e. at a fixed value of L corresponding to E c.m. = 91.2GeV) as a function of the jet resolution scale k ⊥,Lu .This result is presented in Sect.6.1.Second, as a systematic check (see Sect. 7), we evaluate the analytic expression, eq. ( 13), as explained in the following paragraph.To find the N ch.qq (L qq , k ⊥,Lu ) terms in eq. ( 16), we utilize only the second of these methods, i.e. the one based on eq.(13).A direct measurement is not straightforward in this case because L qq , unlike L, is a variable quantity.A direct measurement of N ch.qq (L qq , k ⊥,Lu ) at the scales relevant for our analysis would in principle require a determination of the particle multiplicity in two-jet events as a function of p ⊥,Lu for c.m. energies below the Z 0 .So far, no such measurements are available.
Our procedure to determine N ch.qq (J, k ⊥,Lu ) in eqs.( 15) and ( 16) using the analytic expression, eq. ( 13), is as follows.We begin with the compilation of the mean charged particle multiplicity of e + e − hadronic annihilations, N ch.
e + e − versus E c.m. , given in [12].These data are shown as a function of E c.m. by the open symbols in Fig. 3. (For simplicity, not all the LEP-averaged results compiled in [12] are displayed in Fig. 3; however, these data are included in the procedure described below.)We correct these results for the contributions of c and b events so they correspond to our treatment of N ch.qqg as explained in Sect.4, i.e. to uds events only.This correction is a bin-by-bin multiplicative factor determined from Jetset, given by the ratio of N ch.
e + e − between uds and udscb events as a function of E c.m. .Jetset has been found to provide an accurate description of the multiplicity of both the uds flavor and flavor inclusive samples in e + e − hadronic annihilations (see for example [9,10,35]), justifying its use for this correction.The predictions of Jetset for N ch.
e + e − in udscb and uds flavor events are shown by the dashed and dash-dotted curves in Fig. 3.The correction factors have values between 0.92 and 0.96.[12] for the original references) of the inclusive charged particle multiplicity in e + e − annihilations (open symbols) as a function of E c.m. .The LEP results are averages over the LEP experiments (see [12]).The solid points show the corresponding results after corrections to remove the contributions of c and b events.For simplicity the uncertainties of the corrected data are not shown: these are essentially the same as those of the corresponding uncorrected data points.The solid curve is a fit of the next-to-next-to-next-to-leading order (3NLO) expression for quark jet multiplicity [14,36] to the corrected data.The predictions of the Jetset event generator for udscb and uds flavor events are also shown.
The corrected results for N ch.
e + e − , corresponding to the terms N qq (J ′ ) in eq. ( 13) and which we denote N ch.qq (L) in the following, are shown by the solid points in Fig. 3.We fit these corrected data using the next-to-next-to-next-to-leading order (3NLO) analytic expression for the scale dependence of quark jet multiplicity N q [14,36] assuming n f = 5.Two parameters are fitted: an effective QCD scale parameter, Λ ef f. , and an overall normalization constant, K.The fit results are given in the first row of Table 2.The result of the fit is shown by the solid curve in Fig. 3. Using this fitted curve and its derivative, we numerically evaluate eq. ( 13) as a function of θ 1 to determine the biased quark terms N ch.qq (J, k ⊥,Lu ) in eqs.( 15) and ( 16).The derivative is evaluated by differentiating the analytic expression for N q with respect to y = L/2 and using the fitted parameter values for N ch.qq (L) from Table 2. √ y cut with E vis. the event energy.The data have been corrected for initial state photon radiation, detector response, and misidentification of uds events.The uncertainties are the statistical and systematic terms added in quadrature.The statistical uncertainties are smaller than the symbol sizes.The solid curve shows the prediction [15] of the next-to-leading order (NLO) expression, eq. ( 13).The prediction is seem to agree with the data up to scales of about 40 GeV.The dashed curve is the NLO result found by altering the choice of data points used in the fit of N ch.
e + e − versus scale.

6.1
The biased quark jet multiplicities N ch.qq (J, k ⊥,Lu ) In this section we present our direct measurement of the biased quark jet terms N ch.qq (L, k ⊥,Lu ) in eq.(15).We then use these data to test the analytic expression, eq. ( 13).
We apply the Durham jet finder to the sample of uds flavor events described in Sect. 4 using a fixed value of the resolution parameter y cut (see below).Events are retained if exactly two jets are reconstructed.The solid points in Fig. 4 show the charged particle multiplicity of the selected two-jet events for different choices of y cut .We relate y cut to the transverse momentum p ⊥,Lu through p ⊥,Lu = E vis.
√ y cut [16], with E vis. the sum of the particle energy in the event.
The data are corrected for initial state radiation, the effects of the detector, and uds event p ⊥,Lu (GeV)  3: Measurements of the mean charged particle multiplicity of biased two-jet uds flavor events from Z 0 decays as a function of the transverse momentum cutoff p ⊥,Lu used to separate two-and three-jet events.The jets are defined using the Durham jet-finder.The first uncertainty is statistical and the second systematic.These data are shown in Fig. 4.
misidentification, and are plotted as a function of p ⊥,Lu .The solid points in Fig. 4 therefore represent the direct measurements of the biased quark jet terms N ch.qq (L, k ⊥,Lu ) mentioned in Sect. 5.The vertical lines show the total uncertainties, with statistical and systematic terms added in quadrature.The evaluation of systematic uncertainties is discussed in Sect.7. The statistical uncertainties are negligible compared to the systematic terms.These data are listed in Table 3.As part of our evaluation of systematic uncertainties, results analogous to those shown in Fig. 4 are determined using the Cambridge jet finder.For the Cambridge algorithm, the resolution parameter y cut is related to p ⊥,Lu by the expression given above.We do not determine analogous results using the Luclus jet finder because it is not obvious how to relate the resolution parameter d join to p ⊥,Lu in this case.
The solid curve in Fig. 4 shows the prediction for N ch.qq (L, k ⊥,Lu ) found using eq.( 13), evaluated as explained in Sect. 5.The curve is seen to describe the direct measurements fairly well.In particular, the shape of the analytic curve is generally similar to the data, at least for values of scale below about 40 GeV.This suggests that the NLO result, eq. ( 13), is adequate for the purposes of our study.We return to this question in Sect 6.2.
The difference between the results found using the Durham and Cambridge jet finders is the dominant source of systematic uncertainty for the data in Fig. 4 (Sect.7).The results from the Durham algorithm (standard analysis) lie closer to the NLO curve than those from the Cambridge algorithm, especially for small values of θ 1 .
As a systematic check of the NLO prediction for N ch.qq (L, k ⊥,Lu ), we re-evaluated the prediction after excluding the ARGUS and 91 GeV LEP measurements from the fit of N ch.
e + e − versus E c.m. (see Fig. 3).More details are given in Sect.7.This alternate result is shown by the dashed curve in Fig. 4. The difference between the solid and dashed curves is seen to be similar to the difference between the solid curve and symbols.This suggests that the agree-ment of the analytic expression eq. ( 13) with the direct measurements is quite reasonable, once sources of systematic uncertainty are considered.

6.2
Unbiased gluon jet multiplicity versus scale Fig. 5 shows our results for the unbiased charged particle multiplicities of gg events, N ch.gg .The solid points in Fig. 5a are obtained from eq. ( 15) using the direct measurements of N ch.qq (L, k ⊥,Lu ), i.e. the solid points in Fig. 4. 4 The small horizontal bars indicate the statistical uncertainties.The vertical lines show the total uncertainties, with statistical and systematic terms added in quadrature.The evaluation of systematic uncertainties is discussed in Sect.7. The asterisks and open symbols in Fig. 5b show the corresponding results from eqs. ( 15) and ( 16) using the calculated expressions for N ch.qq (J, k ⊥,Lu ) from eq. ( 13).Thus the asterisks in Fig. 5b are derived using the solid curve in Fig. 4, for example.The two sets of results from eq. ( 15) (solid points in Fig. 5a and asterisks in Fig. 5b) are seen to be very similar.This is consistent with the observation in Sect.6.1 that eq. ( 13) is generally adequate for the purposes of our study.
The results from eq. ( 15) are plotted as a function of Q = p ⊥,Lu and those from eq. ( 16) as a function of Q = p ⊥,Le .The variable p ⊥,Lu extends over a smaller range of scale than p ⊥,Le because the factor of s in the denominator of the Lund definition of transverse momentum is always larger in Y events than the factor of s qq in the Leningrad definition, see eqs.(11) and (12).This effect is also visible in Fig. 2b and Table 1, i.e. p ⊥,Lu has a maximum of 30.5 GeV compared to 53.1 GeV for p ⊥,Le .Note that with θ 1 ≥ 35 • as in our standard analysis, the lower bound on the gluon jet scale is 11.1 GeV for p ⊥,Lu and 16.8 GeV for p ⊥,Le , see Fig. 2b.Included in Fig. 5 are direct measurements of the inclusive charged particle multiplicity of unbiased gluon jets from the CLEO and OPAL Collaborations, currently the only results which utilize a hemisphere definition of gluon jets.These results are shown by the triangular symbols.The three data points at scales of about 4-7 GeV are derived from the hadronic component of Υ(1S) → γgg events [8].The scale Q for these data is given by the invariant mass of the hadronic system.Similarly, Υ(3S) → γχ b2 (2P) → γgg events provide the result at Q = 10.3GeV [7], with the scale given by the χ b2 (2P) mass.The measurement at Q = 80.2 GeV [11] is determined using hadronic Z 0 decays: Z 0 → qqg incl., where g incl.refers to a gluon jet hemisphere recoiling against two almost collinear identified quark jets q and q in the opposite hemisphere [37].For our purposes the g incl.measurement from [11] has been multiplied by a factor of two both for the multiplicity and energy scale so it corresponds to gg "two-jet events" analogous to the other data in Fig. 5. Fig. 5 also shows the predictions of the Herwig [38,39] and Jetset Monte Carlo event generators for the inclusive charged particle multiplicity of gg events versus Q = E c.m. .For these results we use version 6.2 [39] of Herwig with the following changes to the default parameter set: the parameter specifying the angular smearing of udsc quark clusters, CLSMR (1), is set to 0.40, the parameter controlling the mass spectrum of b quark clusters, PSPLT(2), is set to 0.33, and the parameter controlling the production of decuplet baryons, DECWT, is set  (a) Results from eq. ( 15) using the measured biased quark jet terms N ch.qq (L, k ⊥,Lu ).The small horizontal bars indicate the statistical uncertainties.The vertical lines show the total uncertainties, with statistical and systematic terms added in quadrature.The solid curve is a fit of the 3NLO expression for gluon jet multiplicity [14,36] to the data.(b) The corresponding results from eqs. ( 15) and ( 16) using the NLO expression for the biased quark jet terms N ch.qq (J, k ⊥,Lu ), eq. ( 13), where J = L for the N ch.gg (k ⊥,Lu ) results and J = L qq for the N ch.gg (k ⊥,Le ) results.For both (a) and (b), the triangles show measurements of the inclusive charged particle multiplicity of unbiased gluon jets from the CLEO [7,8] and OPAL [11] Collaborations.The predictions of the Herwig and Jetset Monte Carlo event generators are also shown.to 0.70.Herwig 6.2 with these changes is found to provide a good description of the global features of the OPAL data including those of flavor selected samples.
From Fig. 5 it is seen that the Monte Carlo predictions describe the direct measurements of unbiased gluon jet multiplicity (triangle symbols) well.The results from the present analysis based on eq. ( 15) (solid points in Fig. 5a and asterisks in Fig. 5b) are also well described by the Monte Carlo curves.In contrast, the results from eq. ( 16) (open symbols in Fig. 5b) are generally well above the Monte Carlo predictions, and -if extrapolated to lower and higher energies -appear inconsistent with the direct measurements from CLEO and OPAL as well.We note that both Jetset and Herwig have been found to provide a good description of gluon jet properties, including multiplicity, in other studies (see for example [2]- [4], [11], [40]).We therefore conclude that the equation based on the Lund definition of the gluon jet scale, eq. ( 4), yields results which are more consistent with other studies than the equation based on the Leningrad definition, eq. ( 5).Our data thus allow a fairly clear discrimination between the two possibilities eqs.( 4) and ( 5).We henceforth restrict our analysis of gluon jets to the former set of results.Numerical values for the gluon jet measurements based on eq. ( 4) (and ( 15)) are presented in Table 4.The analogous results based on eq. ( 5) (and ( 16)) are presented in Table 5.For information, we include results corresponding to all θ 1 bins in these tables, not just the bins above 35 • used in our standard analysis.Note that the theoretical formalism is not expected to be valid for small values of θ 1 .This explains the negative values of multiplicity determined for the first bin, 0 ≤ θ 1 ≤ 5 • .

Determination of the ratios
Having established that the formalism of [15] yields consistent results with the direct measurements of unbiased gluon jet multiplicity if the equation for three-jet event multiplicity based on the Lund definition of transverse momentum is employed (Sect.6.2), we proceed to a comparison of the unbiased gluon and quark jet terms and to a test of the theoretical expressions for the ratios r, r (1) and r (2) , see eqs. ( 1)- (3).For this purpose we use the gluon jet data shown by the solid points in Fig. 5a since these results are derived using the direct measurements of N ch.qq (L, k ⊥,Lu ) rather than the NLO expression, eq. ( 13).They therefore rely on fewer theoretical assumptions.
We employ the following procedure to determine r, r (1) and r (2) from experiment.The corrected unbiased quark and gluon jet multiplicities, shown by the solid points in Figs. 3  and 5a, are separately fitted using the 3NLO expressions for N q and N g [14,36], respectively.To better constrain the results, the measurements of inclusive gluon jet multiplicity from CLEO at 10.3 GeV and OPAL at 80.2 GeV are included in the gluon jet fit.The fitted parameters are Λ ef f. and the normalization constant K. (For the quark jet data, this is the same fit discussed in Sect.5.) The results of the fits are shown by the solid curves in Figs. 3 and 5a.The fitted parameters and corresponding χ 2 per degree-of-freedom are listed in Table 2.The χ 2 per d.o.f. and overall description of the data are seen to be reasonable in both cases.The ratio of the fitted expressions for N g and N q defines r, eq. ( 1).We calculate the first and second derivatives of the analytic equations for N g and N q with respect to y and evaluate the resulting expressions gg as a function of energy scale Q = p ⊥,Lu , obtained using eq.( 15).The corresponding bins of θ 1 in Y events are also indicated.The results for θ 1 ≥ 35 • , listed in the bottom part of the table, are used in our standard analysis.These data are shown by the solid points in Fig. 5a.For information, we also list the results obtained for smaller values of θ 1 .Note that the theoretical formalism used to determine these results is not expected to be valid for small θ 1 .This explains the negative value of multiplicity derived for θ 1 < 5 • .The corresponding results for the multiplicity ratio r between gluon and quark jets are also given.The results for r for θ 1 ≥ 35 • are shown by the solid points in Figs. 6  and 8a.gg as a function of energy scale Q = p ⊥,Le , obtained using eq.( 16).Note that the theoretical formalism used to determine these results is not expected to be valid for small θ 1 .This explains the negative value of multiplicity derived for θ 1 < 5  The triangle and square symbols show results based on direct measurements of multiplicity in unbiased jets from the CLEO [7,8], ARGUS [41] and OPAL [11] Collaborations.The dashdotted and long-dashed curves show theoretical predictions based on analytic [14,15] and numerical [13] techniques.
Our results for r are shown in Fig. 6, and those for r (1) and r (2) in Fig. 7.The central results are indicated by solid curves.The shaded bands show the statistical uncertainties.The overall uncertainties, with statistical and systematic terms added in quadrature, are shown by the open bands.The bands cover the energy range of data in the gluon jet fit, from 10.3 to 80.2 GeV.The short-dashed curves show extrapolations of the results to lower and higher scales.At 30 GeV, a typical scale for the gluon jet data in these figures, the ratios are determined to be: r (30 GeV) = 1.422 ± 0.006 ± 0.051 , (17) r (1) (30 GeV) = 1.761 ± 0.013 ± 0.070 , (18) r (2) where the first uncertainty is statistical and the second systematic.At 80 GeV, corresponding to the highest energy gluon jet data in the fits, the results are: The ratios of slopes r (1) and of curvatures r (2) between unbiased gluon and uds quark jets as a function of scale.The curves and bands have the same meaning as in Fig. 6.
Analogous results for r versus scale were previously published in [4,18,19] and for r (1) in [18,19].These earlier studies were based on biased jet samples, however.Therefore the quantitative results from those studies cannot be compared to theoretical predictions without ambiguity.
Our result for r (2) is the first experimental determination of that quantity.
Included in Fig. 6 are the results for r derived from the individual gluon jet measurements in Fig. 5a.The solid points in Fig. 6 are obtained by dividing the results shown by the solid points in Fig. 5a by the results of the 3NLO parametrization of quark jet multiplicity at the corresponding scales.Numerical values of these results for r are compiled in Table 4.The triangle and square symbols in Fig. 6 are obtained as follows: (i) for the OPAL data point at 80.2 GeV, the published result [11] is used, based on uds quark events analogous to the present study; (ii) for the three CLEO points at 4-7 GeV, the published results [8] are used after applying a correction obtained from Jetset to remove the contributions of c events from the quark jet samples; (iii) for the point labelled "CLEO/ARGUS," the CLEO measurement of gluon jet multiplicity at 10.3 GeV [7] is divided by the ARGUS measurement of N ch.
e + e − at 10.5 GeV [41] (see Fig. 3), after applying a correction to account for c and b events in the quark jet sample and for the small energy difference between the gluon and quark jet results.
As a consistency check, we also determined the ratios r and r (1) using the technique employed in [19].The corrected quark jet results (solid points in Fig. 3) were fitted using a polynomial expression: with L given in eq. ( 6), and where a 0 , a 1 and a 2 are fitted parameters.We find a 0 = 13.20±0.67,a 1 = −2.91±0.12 and a 2 = 0.2818±0.0057,where the uncertainties are statistical.The χ 2 (d.o.f.) of the fit is 3.74 (15).We then fitted the 19 gluon jet measurements in our analysis (solid points in Fig. 5a and the 10.3 and 80.2 GeV results from CLEO and OPAL) using the expression: with the results b 0 = −4.07± 0.27, b 1 = 1.759 ± 0.030 and χ 2 (d.o.f.) = 31.5(17).The result for r, given by the ratio of eqs.( 24) and ( 23), is indistinguishable from the solid curve in Fig. 6.The result for r (1) , obtained by differentiation of eq. ( 24), is, by construction, a constant: ).This result is consistent with that shown in Fig. 7a, see also eq. ( 18).Not only is r (1) constant with this method, but r (2) ≡ b 1 = r (1) , i.e. no distinction is made between the ratios of slopes and curvatures.Therefore this alternate technique does not provide as much information as our standard procedure described above.
To obtain a consistency check of our result for r (2) , we applied the following procedure.The 19 gluon jet measurements were clustered into 16 groups of four contiguous measurements.Thus, labelling the data points from 1 to 19 in order of increasing Q, the first group is composed of data points 1-4, the second of 2-5, etc.A straight line fit was performed for each group: with a g and b g the fitted parameters.By fitting to four points at a time, rather than e.g. to two points only, the influence of bin-to-bin fluctuations is reduced.We then evaluated the fitted polynomial expression eq. ( 23) at the scale values of the gluon jet measurements to obtain 19 corresponding results for quark jets, which were then fitted in the same manner: A scale value was associated with each slope point (a g ) i or (a q ) i using the arithmetic mean of the four Q values in the fits.As a systematic check, we also defined these scales using the weighted means of the Q values of the four associated gluon jet measurements.This alternate method yielded consistent results with those presented below.The slopes (a g ) i and (a q ) i were in turn fitted using straight lines: with c i and d i (i = q, g) the fitted terms.The ratio r (2) found using this technique is a constant given by r (2) = c g /c q .We find r (2) = 2.25 ± 0.16 (stat.),consistent with our result in Fig. 7b, see also eq. ( 19).
The procedure described in the previous paragraph can also be used to obtain a second consistency check for r (1) .Taking the weighted mean of the 16 results for the ratio of slopes r i = (a g /a q ) i yields r (1) = 1.68 ± 0.02 (stat.),consistent with the result in eq. ( 18) to within about 1.0 standard deviation of the total uncertainties.We note that the method to determine r (1) and r (2) based on eqs.( 25)-( 28) is more direct than our standard procedure because no functional form is assumed for how the gluon jet multiplicity varies with scale.The resulting uncertainties are considerably larger, however, and no information on the scale dependence of r (1) and r (2) is obtained.Therefore we retain the analysis based on the 3NLO parametrizations of multiplicities for our standard results.
Analytic expressions for r, r (1) and r (2) , valid to 3NLO, were recently presented by Capella et al. [14].All three ratios are predicted to converge to the ratio of QCD color factors, C A /C F = 2.25, in the asymptotic limit of large jet energies.This convergence is predicted to be more rapid for r (1) than r [14,18] and yet more rapid for r (2)  [14].Our data are in agreement with this prediction, i.e. r < r (1) < r (2) < 2.25 for the scales accessible in our study, see eqs. ( 17)- (22).Note that we also observe the hierarchy r < r (1) < r (2) using the alternate analysis strategies summarized by eqs.( 23)- (28).
To make the comparison of our results for r, r (1) and r (2) more clear, the data from Figs. 6 and 7 are displayed together in Fig. 8a.The dominant uncertainties in our analysis are systematic (see Sect. 7) and are correlated between the three ratios.Therefore, in Fig. 8b, we present results for the differences r (2) − r, r (1) − r and r (2) − r (1) , for which correlated systematic uncertainties partially cancel.At 30 GeV, we find: r (1) − r (30 GeV) = 0.341 ± 0.014 ± 0.075 , r (2) − r (1) (30 GeV) = 0.221 ± 0.024 ± 0.070 , where the first uncertainty is statistical and the second systematic.These results demonstrate that r (2) exceeds r (1) by 3.0 standard deviations of the total uncertainties, with yet more significant deviations from zero observed for the other two differences of ratios.
Figure 8: (a) The experimental results for r, r (1) and r (2) from Figs. 6 and 7, collected together; (b) the differences between ratios r (2) − r, r (1) − r and r (2) − r (1) as a function of scale.The central results are shown by solid curves; the shaded and open bands show the corresponding statistical and total uncertainties, respectively.

Comparison to theoretical expressions
The 3NLO analytic predictions of Capella et al. [14] for r, r (1) and r (2) are shown by the longdashed curves in Figs. 6 and 7.For Q = 30 GeV, the ratios are predicted to be 1.73, 1.87 and 2.00, respectively, assuming n f = 5 and Λ = 0.20 GeV.Thus, the analytic predictions for r and r (1) exceed the corresponding experimental results in eqs.( 17)-( 18) by about 22% and 6%, while the theory agrees with the data, eq. ( 19), for r (2) .Therefore the analytic predictions agree better with the data the higher the level of differentiation of the multiplicity curves.This suggests that higher order corrections are smaller for r (2) than r (1) , and for r (1) than r.
We note there are ambiguities in the analytic calculation related to the effective number of active quark flavors and the scale Q at which to evaluate the strong coupling strength.The 3NLO predictions shown here are based on n f = 5 (as stated above) and Q = E c.m. .Other plausible choices such as n f = 4 and Q = E c.m. /4 change the analytic prediction for r by as much as about 10% and can bring the theoretical curve into better agreement with the data, see Fig. 2 of [9] for example.
A second QCD prediction for r versus scale was recently presented by Lupia and Ochs [13].This prediction is based on numerical, rather than analytic, techniques, incorporating more accurate phase space limits for the emission of soft gluons and a more complete treatment of energy conservation than the analytic expression.Unlike the analytic result, a prediction is presented only for r, however, not for r (1) or r (2) .We derive a prediction for r (1) from the numerical calculation in the following manner.We begin with numerical predictions for the individual gluon and quark jet multiplicities, 5 given in steps of 0.02 units in y with y = ln (Q/Q c ) where Q c = 0.507 GeV [13].We make linear interpolations between the points and use the slopes of these lines to estimate dN g /dy and dN q /dy, which are then used to determine r (1) as a function of scale Q.We attempted to extend this technique to obtain a prediction for r (2)  as well.The derived results for r (2) were very scattered from point to point as a function of Q because of limited numerical precision.We therefore do not quote a prediction for r (2) based on the calculation of [13].The predictions of the numerical calculation for r and r (1) are shown by the short-dasheddotted curves in Figs. 6 and 7a.At 30 GeV, the numerical calculation predicts r to be 1.45, in agreement with our measurement at that scale, eq. ( 17).The corresponding result for r (1)  is 1.64, about 7% below the data, eq. ( 18).For r, the agreement with data is therefore better for the numerical method than for the analytic calculation.The level of agreement of the two predictions with data is similar for r (1) .The difference between the analytic and numerical results for r (long-dashed and short-dash-dotted curves in Fig. 6) is believed to arise from the differences in phase space limits and energy conservation noted above.For further discussion, see [12].
Finally, theoretical predictions for r, r (1) and r (2) can be derived from the formalism of Edén and Gustafson [15], using eq.( 14).To obtain a prediction for r (1) , we begin with the experimental results for the slope of quark jets, dN qq (L)/dL, determined in the manner described in Sect.6.3.1 for dN qq (y)/dy, with y = L/2.We insert this result in eq. ( 14) to obtain a prediction for dN gg (L)/dL.Dividing this result by the measured dN qq (L)/dL yields the prediction for r (1) .To obtain a prediction for r (2) , we differentiate eq. ( 14) with respect to L and follow an analogous procedure to that described for r (1) .To obtain a prediction for r, we fix the gluon jet multiplicity at 80.2 GeV to the g incl.measurement from OPAL [11].We then use the prediction for dN gg (L)/dL determined as described above to evaluate the gluon jet multiplicity at other scales and divide the result by the 3NLO expression for the quark jet multiplicity fitted to data (solid curve in Fig. 3).
The predictions obtained in this manner are shown by the long-dash-dotted curves in Figs. 6  and 7.The results are seen to be in good overall agreement with the data.We note, however, that these predictions are based on the experimental measurements of quark jet multiplicities.Furthermore, the prediction for r utilizes the measured gluon jet multiplicity at Q = 80.2 GeV as stated above.Therefore, the predictions we derive based on [15] are not independent of the data.

Measurement of C
Measurements of the energy dependence of unbiased particle multiplicities in gluon and quark jets can also be used to determine an effective value of the ratio of QCD color factors C A /C F .Such a possibility was first discussed in [18].To perform this measurement, we integrate eq. ( 14) to obtain: with N gg (L 0 ) an integration constant which serves as the reference point for N gg (L), and where We evaluate eq. ( 33) numerically using the experimental results for dN qq (Q)/dQ discussed in Sect.6.3.1.As the reference point N gg (L 0 ) in eq. ( 32), we use the direct measurement of gluon jet multiplicity from OPAL [11].Thus L 0 = 2 ln (Q 0 /Λ), with Q 0 = 80.2 GeV and N gg (L 0 ) = 28.56 [11].For Λ we used the fitted result for Λ ef f. for quark jets from Table 2, Λ = 0.19 GeV.Note that the results for C A /C F are independent of this choice.
A one parameter fit of eq. ( 32) is made to the gluon jet measurements in Fig. 5a, excluding the CLEO data below 7 GeV, to obtain where the first uncertainty is statistical.The second uncertainty is a systematic term evaluated as described in Sect.7. The third uncertainty is an additional systematic term defined by the difference observed if the CLEO measurement of N gg at 10.3 GeV is used as the reference point, rather than the OPAL measurement at 80.2 GeV.The χ 2 (d.o.f.) of the fit is 29.5 (17).The result of our fit, eq. ( 34), is consistent with the QCD value C A /C F = 2.25.
The fitted expression of eq. ( 32) yields a curve which is almost identical to the 3NLO curve in Fig. 5a: therefore we do not display the fitted curve of eq. ( 32) in addition.In Fig. 9 we present the result for the event multiplicity in Y events.The data in this figure are the The mean charged particle multiplicity of three-jet uds flavor Y events from Z 0 decays compared to the analytic expression, eq. ( 4), for which the gluon jet multiplicity term N gg (k ⊥,Lu ) is determined from a one parameter fit of eq. ( 32) to unbiased gluon jet data and the quark jet term N qq (L, k ⊥,Lu ) is measured.The open band shows the evaluated uncertainty of the analytic result.Note the theory is not expected to be applicable for small values of θ 1 because the events in this case do not have a well defined three jet event structure.
same as in Fig. 2a.The curve is obtained from eq. ( 4), with N gg (k ⊥,Lu ) the fitted expression of eq. ( 32) and N qq (L, k ⊥,Lu ) the measured result from Sect.6.1.The open band shows the uncertainty of the analytic result, evaluated using the uncertainty from eq. ( 34) in addition to those shown for the N qq (L, k ⊥,Lu ) terms in Fig. 4. The curve is seen to provide an accurate description of the data to within its uncertainty for θ 1 values larger than approximately 10 • .This represents a considerable improvement compared to the results of [42] based on earlier theoretical formalism [17], in which the corresponding fitted analytic expression was found to describe the data for θ 1 values larger than about 60 • only (see Fig. 4 of [42]).
To test our technique of determining C A /C F , we apply our analysis to parton level Monte Carlo events generated with a c.m. energy of 10 TeV.The parton level is defined using finalstate partons, i.e. those which are present after the termination of the parton shower.We use Herwig version 6.2 with the parameters discussed in Sect.6.2.Applying our analysis to parton level events with a very large energy tests our method to determine C A /C F since the parton level asymptotic result must be consistent with 2.25 for the analysis to be considered valid, see [42] for further discussion.The parton level Herwig events at 10 TeV yield C A /C F = 2.285 ± 0.010 (stat.)± 0.047 (ref. ), where the central value is defined using the Herwig parton level prediction for N gg at 80 GeV as the reference point.The first uncertainty is statistical while the second uncertainty is defined by the difference found if the Herwig prediction for N gg at 10 GeV is used as the reference point instead (i.e. this is analogous to how the third uncertainty of eq. ( 34) is defined).Thus our one parameter fit method yields an asymptotic result of 2.25 to within the uncertainties associated with the analysis procedure, implying that the similarity of our measurement eq. ( 34) with the QCD value of C A /C F is not merely coincidental.For information, the parton level Herwig result at E c.m. = 91 GeV is 2.03 ± 0.01 (stat.)if the prediction for N gg at 80 GeV is used as the reference point and 2.35±0.02(stat.)if the prediction at 10 GeV is used.The analogous results for hadron level Herwig events are 2.34 ± 0.01 and 2.07 ± 0.02, respectively.We therefore observe a much stronger dependence on the choice of the normalization point for the 91 GeV Monte Carlo results than for the 10 TeV results or for the data.The Herwig results at 91 GeV are nonetheless consistent with 2.25, once the systematic uncertainty related to choice of the normalization point is considered.

Systematic uncertainties
To evaluate systematic uncertainties, the analysis was repeated with the following changes to the standard analysis.The full differences between the standard results and those found using each of these changes were used to define symmetric systematic uncertainties.For items 1, 5 and 7, listed below, the largest of the described differences with respect to the standard result was assigned as the systematic uncertainty.The uncertainties were added in quadrature to define the total systematic uncertainties.
For the measurement of the gluon jet multiplicity N ch.gg (Fig. 5), the following changes were made.Items 1-3 and 5 in this list were also applied to the study of two-jet events, see Sect.6.1.The items are listed in roughly decreasing order of the size of their contribution to the total uncertainties of the N ch.
1.The jet finding was performed using the Cambridge and Luclus jet finders, rather than the Durham jet finder.For the Luclus jet finder, we did not perform a direct measurement of the biased quark jet terms N ch.qq (J, k ⊥,Lu ) since it is not obvious how to relate the resolution scale d join to p ⊥,Lu as mentioned in Sect.6.1.Instead we used the analytic expression, eq. ( 13), for the determination of the gluon jet multiplicities N ch.
gg (L).The difference between the Durham and Cambridge results was used to define the systematic uncertainty due to jet finding for the two jet event measurements shown in Fig. 4.
2. Tracks selected for the uds tagging procedure were required to have a signed impact parameter which satisfied d sign /σ d sign > 1.5, rather than d sign /σ d sign > 3.0; at the same time we required | cos(θ thrust )| to be less than 0.8, rather than 0.9.The purpose of these two changes was to increase the purity of tagged uds events; the estimated uds event purity increased to 85.1% from 78.5% while the number of selected events decreased to 8610 from 22 365.The more restrictive condition on | cos(θ thrust )| yields a higher uds purity because a larger fraction of tracks in an event traverses the silicon microvertex detector.
3. The Herwig Monte Carlo event generator was used to determine the corrections for initialstate photon radiation, detector response, and uds event and gluon jet misidentification, rather than Jetset.For this purpose, we used version 5.9 of the program [38] with the parameter values in [32], except that the cluster mass cutoff CLMAX was increased from 3.40 to 3.75 GeV/c 2 to improve the model's description of the mean charged particle multiplicity in inclusive hadronic Z 0 decays.
4. For our measurements of the gluon jet multiplicities N ch.gg based on eq. ( 15), we used the analytic expression, eq. ( 13), rather than the direct measurement of N ch.qq (L, k ⊥,Lu ).
5. Charged tracks alone were used for the data and Monte Carlo samples with detector simulation, rather than charged tracks plus electromagnetic clusters.As an additional check on the track selection, the standard analysis (tracks and clusters) was varied by increasing the minimum transverse momentum of charged tracks with respect to the beam axis from from 0.05 GeV/c to 0.15 GeV/c.
For the two-jet events (Fig. 4), the dominant systematic term was from item 1 at low p ⊥,Lu and from item 3 at intermediate and high p ⊥,Lu .
In addition to the above, the following changes were made to evaluate systematic uncertainties for the ratios r, r (1) and r (2) .6.The ARGUS and 91 GeV LEP measurements of N ch.
e + e − were excluded from the fit of quark jet multiplicity (solid curve in Fig. 3).These two data points have the smallest uncertainties of the results in the fit and were excluded to potentially maximize a systematic variation.The results of the alternate fit are Λ ef f. = 0.261 ± 0.061, K = 0.154 ± 0.014 and χ 2 (d.o.f.) = 1.5 (14).Besides the evaluation of a systematic uncertainty for r, r (1) and r (2) , this alternate result for N ch.qq (L) was used to derive the NLO prediction shown by the dashed curve in Fig. 4.
7. Y events with 50 • ≤ θ 1 ≤ 120 • were used in the fit of gluon jet multiplicity versus scale, rather than 35 • ≤ θ 1 ≤ 120 • ; as an alternate check, the ratios were determined after increasing the inclusive gluon jet measurement from CLEO at 10.3 GeV by its one standard deviation total uncertainty, then decreasing it by one standard deviation, and repeating the same operation for the inclusive gluon jet measurement from OPAL at 80.2 GeV.
For the ratios, the largest systematic terms were from items 1, 6 and 7.The contributions of items 2 and 5 were by comparison negligible, with the other items intermediate.For our measurement of C A /C F (Sect.6.4), increasing and decreasing the inclusive OPAL measurement at 80.2 GeV by one standard deviation (item 7) provided the largest contribution to the overall uncertainty.Item 1 provided the second largest contribution while the contributions of the other items were much smaller.
For the data listed in Tables 4 and 5 and the corresponding results in Figs. 5, 6 and 8a, the systematic uncertainty evaluated for each bin was averaged with the results from its two neighbors to reduce the effect of bin-to-bin fluctuations.The single neighbor was used for bins on the endpoints of the distributions.

10 π 2
/27 − 3/2, and n f equals the number of active quark flavors, taken to be n f = 5.C A and C F are the QCD color factors, 3 and 4/3, respectively.

Figure 3 :
Figure3: Measurements (see[12] for the original references) of the inclusive charged particle multiplicity in e + e − annihilations (open symbols) as a function of E c.m. .The LEP results are averages over the LEP experiments (see[12]).The solid points show the corresponding results after corrections to remove the contributions of c and b events.For simplicity the uncertainties of the corrected data are not shown: these are essentially the same as those of the corresponding uncorrected data points.The solid curve is a fit of the next-to-next-to-next-to-leading order (3NLO) expression for quark jet multiplicity[14,36] to the corrected data.The predictions of the Jetset event generator for udscb and uds flavor events are also shown.

Figure 4 :
Figure 4: The mean charged particle multiplicity of biased two-jet uds flavor events from Z 0 decays selected using the Durham jet finder.The results are presented as a function of p ⊥,Lu , related to the resolution parameter y cut of the jet finder by p ⊥,Lu = E vis.√ y cut with E vis. the

Figure 5 :
Figure 5: The mean charged particle multiplicity of unbiased gg events as a function of scale.(a) Results from eq. (15) using the measured biased quark jet terms N ch.qq (L, k ⊥,Lu ).The small horizontal bars indicate the statistical uncertainties.The vertical lines show the total uncertainties, with statistical and systematic terms added in quadrature.The solid curve is a fit of the 3NLO expression for gluon jet multiplicity[14,36] to the data.(b) The corresponding results from eqs. (15) and (16) using the NLO expression for the biased quark jet terms N ch.qq (J, k ⊥,Lu ), eq.(13), where J = L for the N ch.gg (k ⊥,Lu ) results and J = L qq for the N ch.gg (k ⊥,Le ) results.For both (a) and (b), the triangles show measurements of the inclusive charged particle multiplicity of unbiased gluon jets from the CLEO[7,8] and OPAL[11] Collaborations.The predictions of the Herwig and Jetset Monte Carlo event generators are also shown. •

Figure 6 :
Figure6: Results for the ratio r of the mean charged particle multiplicities between unbiased gluon and uds quark jets as a function of scale.The central result is shown by the solid curve; the shaded and open bands show the corresponding statistical and total uncertainties, respectively (the statistical uncertainties are barely visible).The short-dashed curve is an extrapolation of the fit.The solid points show measurements based on gluon jet results from the present analysis.The triangle and square symbols show results based on direct measurements of multiplicity in unbiased jets from the CLEO[7,8], ARGUS[41] and OPAL[11] Collaborations.The dashdotted and long-dashed curves show theoretical predictions based on analytic[14,15] and numerical[13] techniques.

Figure 7 :
Figure 7: (a) and (b):The ratios of slopes r(1) and of curvatures r(2) between unbiased gluon and uds quark jets as a function of scale.The curves and bands have the same meaning as in Fig.6.

Figure 9 :
Figure9: The mean charged particle multiplicity of three-jet uds flavor Y events from Z 0 decays compared to the analytic expression, eq.(4), for which the gluon jet multiplicity term N gg (k ⊥,Lu ) is determined from a one parameter fit of eq.(32) to unbiased gluon jet data and the quark jet term N qq (L, k ⊥,Lu ) is measured.The open band shows the evaluated uncertainty of the analytic result.Note the theory is not expected to be applicable for small values of θ 1 because the events in this case do not have a well defined three jet event structure.

Table 2 :
Results of fits of the 3NLO expressions for quark and gluon jet multiplicities to measurements of charged particle multiplicity in unbiased qq and gg events, see the solid curves in Figs.3 and 5a.The fits are performed assuming five active quark flavors.The uncertainties are statistical only.The χ 2 per degree-of-freedom (d.o.f.) of the fits are also given.The χ 2 values are determined using the statistical uncertainties of the data and not the systematic terms.

Table 4 :
Measurements of unbiased gluon jet multiplicity N ch.

Table 5 :
Measurements of unbiased gluon jet multiplicity N ch.