Measurement of Ab_FB using inclusive b-hadron decays

Based on a sample of four million events collected by ALEPH from 1991 to 1995, a measurement of the forward-backward asymmetry in Z ->b-barb decays using inclusive final states is presented. High-performance tagging of b-barb events in a wide angular range is achieved using neural network techniques. An optimal hemisphere charge estimator is built by merging primary and secondary vertex information, leading kaon identification and jet charge in a neural network. The average charge asymmetry, the flavour tagging efficiencies and mean b-hemisphere charges are measured from data and used to extract the pole b asymmetry in the Standard Model A^{0,b}_{FB} = 0.1009 +- 0.0027 (stat) +- 0.0015 (syst), corresponding to a value of the effective weak mixing angle of sin^2(theta_W^{eff}) = 0.23193 +- 0.00056.


Introduction
The measurements of the forward-backward asymmetry of b-quarks from e + e − → Z → b b production offer the most precise determination of the weak mixing angle at LEP.This gives a sensitive test of Standard Model [1] predictions of electroweak radiative corrections and hence constrain the allowed range of the Higgs boson mass.
The polar angle distribution of the b-quark in the e + e − → Z → b b process is forwardbackward asymmetric: where θ b is the angle between the incoming electron and the outgoing b-quark.In this paper an improved measurement of A b FB using inclusive b-hadron decays is presented, based on the data sample recorded by ALEPH on the Z from 1991 to 1995.The analysis takes advantage of a reprocessing of LEP I data with improved charged particle tracking and of neural network techniques making maximum use of the available event information.
The method of selecting events containing b b quark pairs is upgraded with respect to the previous analysis [2].The b-tagging algorithm based on the impact parameters of charged particle tracks is complemented with information from displaced secondary vertices, event shape variables and lepton identification.This leads to a 15% reduction of the statistical uncertainty.The estimate of the b-quark direction is based upon the hemisphere charge method [3], but also incorporates information from fast kaon tagging and separated primary and secondary vertex charge estimators.The fraction of incorrect charge tags is reduced by 10% with respect to the previous measurement, implying a 20% further reduction of the statistical uncertainty.Another new feature is the control of systematic uncertainties by use of double tag methods for both flavour and charge tags, which yields a reduction of the systematic uncertainty by about a factor of two.
Although lepton identification is used to select events in the b-tag algorithm, no use is made in the hemisphere charge method of the lepton beyond that of an ordinary charged track in the detector.This reduces the correlation with the forthcoming ALEPH A b FB analysis based on semileptonic final states.

The method
The measurement of A b FB requires knowledge of the direction of the b-quark from Z → b b decay.The quark-antiquark axis is estimated using the reconstructed thrust axis, the direction of which is given throughout this paper by a positive value for the cosine of the thrust polar angle, cos θ.Each event is then divided into two hemispheres, F and B, by a plane perpendicular to the thrust axis.The forward hemisphere, F, is defined as the one into which the incoming electron points.The F-B orientation of the b-quark is determined on a statistical basis by estimating the hemisphere charges.
Hemisphere charges, Q F and Q B , are formed using a neural network designed to optimise the separation power between b-and b-quarks.The neural network also performs well for lighter quark flavours, albeit not in an optimal way.A forward-backward asymmetry for flavour f at a given value of cos θ is then proportional to the mean charge flow, Q f FB , 1 between forward and backward hemispheres in pure f f events : where n f F (n f B ) is the number of events with the primary quark from Z decay emitted into the forward (backward) hemisphere, and δ f = Q f − Q f is the average difference between the charges measured in the hemispheres of the quark and anti-quark, called the charge separation for flavour f .As shown in [4], the same sample of events used to measure Q f FB can also be used to extract δ f .This can be understood by considering a single hemisphere charge measurement, Q f , which can be written as : where R represents the measurement fluctuation due to fragmentation and detector effects.
The product of the two hemisphere charges then averages to : The measurement fluctuation correlation, R f R f , arises from effects of charge conservation, sharing a common event axis and crossover of particles close to the hemisphere boundary.It is then useful to define : where Q f tot is the total charge measured in an f f event.Thus, the observable δf is equal to δ f to within a correction term, k f , depending on the polar angle and taking an average value of 9% for heavy flavours in this analysis.Compared with the charge separation itself, the k f correction term is less sensitive to the details of quark fragmentation and detector resolution [5].Relying on the Monte Carlo prediction of k f , it is possible to extract δ f by fitting the quantities δf to a range of measurements of the flavour combined δ using the relation: where P f are flavour purities of the event sample under study.As described in Section 8, these purities are also measured in the data.
For each value of cos θ the mean charge flow, the purities and the charge separations are measured.The b-quark asymmetry is then determined according to Equation (1), averaged over quark flavours: where (5) Rearranging Equation (4), The b-quark forward-backward asymmetry, A b FB , is extracted from a simultaneous fit to the measurements in nine angular bins in the range 0 < cos θ < 0.95.A tiny, but non-zero, correction is applied to take into account the flavour dependence of the acceptance in each angular bin.

The ALEPH detector
The ALEPH detector is described in detail in [6], and its performance in [7].The tracking system consists of two layers of double-sided silicon vertex-detector (VDET), an inner tracking chamber (ITC) and a time projection chamber (TPC).The VDET single hit resolution is 12 µm at normal incidence for both the rφ and rz projections and 22 µm at maximum polar angle.The polar angle coverage of the inner and outer layers are | cos θ| < 0.84 and | cos θ| < 0.69 respectively.The ITC provides up to 8 rφ hits at radii 16 to 26 cm relative to the beam with an average resolution of 150 µm and has an angular coverage of | cos θ| < 0.97.The TPC measures up to 21 space points per track at radii between 38 and 171 cm, with an rφ resolution of 170 µm and a z resolution of 740 µm and with an angular coverage of | cos θ| < 0.96.In addition, the TPC wire planes provide up to 338 samples of ionisation energy loss (dE/dx).
Tracks are reconstructed using the TPC, ITC and VDET which are immersed in a 1.5T axial magnetic field.This provides a transverse momentum resolution of σ(1/p T ) = 0.0006 (GeV/c) −1 for 45 GeV muons.Multiple scattering dominates at low momentum and adds a constant term of 0.005 to σ(p T )/p T .
Outside of the TPC, the electromagnetic calorimeter (ECAL) consists of 45 layers of lead interleaved with proportional wire chambers.The ECAL is used to identify photons and electrons and gives an energy resolution σ(E)/E = 0.18/ E/GeV + 0.009 for isolated particles.The hadron calorimeter (HCAL) is formed by the iron of the magnet return yoke interleaved with 23 layers of streamer tubes.It is used to measure hadronic energy and, together with two surrounding layers of muon chambers, to identify muons.
The information from the subdetectors is combined in an energy flow algorithm [7] which gives a list of charged and neutral track momenta.
Recently the LEP I data have been reprocessed using improved reconstruction algorithms.In particular, the VDET hits are distributed among tracks according to a global χ 2 minimisation procedure which improves the hit association efficiency by more than 2%.Information from TPC wires is now used in addition to the pad information to improve the coordinate resolution by a factor of two in z, and by 30% in rφ for low momentum tracks.Similarly the pad information is added to the dE/dx information from TPC wires, increasing the fraction of tracks with a useful dE/dx measurement to almost 100% and providing on average a two sigma separation between pions and kaons with momenta above 2.5 GeV/c.

Monte Carlo simulation
The analysis makes use of a Monte Carlo sample (MC) of 8.1 million simulated hadronic Z decays as well as two dedicated heavy flavour samples of 4.9 million Z → b b decays and 2.4 million Z → cc decays.The simulation is based on JETSET [8] with parameters tuned to reproduce inclusive particle spectra and event shape distributions measured in hadronic Z decays [9].
In the MC simulation, the most relevant physics input parameters have been adjusted, using the re-weighting technique, to the recently measured values [10,11] shown in Table 1.The fragmentation of heavy quarks into hadrons is assumed to follow the model of Peterson et al. [12] with parameters tuned to match the values of the mean heavy hadron fractions of the beam energy reported in the  [10,11] to the Monte Carlo simulation which are used for re-weighting.

The neural net b-tag
Due to the long lifetime and high mass of b-hadrons, b-jets have several characteristic properties.Six discriminating variables are combined using a neural network to tag b-quark jets.A similar scheme is used by ALEPH to identify b-quark jets in searches for neutral Higgs bosons conducted at LEP II [13].
The jets are clustered with the JADE algorithm using a y cut value of 0.02, and for each jet six variables are defined: two of them are lifetime-based; a third one is based on the transverse momentum of identified leptons; the last three are based on jet-shape properties.These quantities are 1.P jet : probability of the jet being a light quark (uds) jet based upon impact parameters of tracks in the jet [14].
2. ∆χ 2 svx : the vertex χ 2 difference between assigning tracks in the jet both to the secondary and primary vertices compared to assigning all tracks to the primary vertex.This is based upon a secondary vertex pattern recognition algorithm which searches for displaced vertices via a three-dimensional grid point search [15]; 3. p T of identified leptons [16] with respect to the jet axis; 4. S b : the boosted sphericity of the jet defined to be the sphericity of final state charged and neutral particles in the rest frame of the jet; 5. Multiplicity/ln(E jet ): the charged and neutral particle multiplicity of the jet divided by the natural logarithm of the jet energy in GeV; 6. Σp 2 T : the sum of squared transverse momenta, p 2 T , of each charged or neutral particle with respect to the jet axis.
Among the six input variables, the lifetime tags have the largest weight.However, the inclusion of the other variables increases the b-quark discriminating power, especially at the most forward angles that are not covered by the vertex detector.A hemisphere b-tag is defined as the maximum neural net output value among the jets in the hemisphere with energies exceeding 10 GeV.
Prior to the application of the lifetime reconstruction algorithms in the MC, the track parameters have been given additional smearing and the VDET additional hit inefficiency according to a procedure described in [17].This is necessary in order to render the MC distributions of the lifetime tags in good agreement with data.Small discrepancies still remain and combine in the neural network to give the differences between the distributions of the b-tag shown in Figure 1.Deviations up to about 5% between data and simulation are seen in some b-tag bins and for this reason the b-quark purity as a function of b-tag is determined from the data themselves.

The neural net charge tag
A neural network is also used to determine hemisphere charges.This technique has been employed earlier for CP violation studies in B 0 / B0 decays as described in [18].In the present analysis information regarding lepton identification is left out from the tag.
In order to achieve high tagging performance independently of the b-hadron species, momentum and polar angle, several charge estimators are combined.These charge tags are complemented with other variables, which do not carry information about the b-hadron charge in themselves, but are correlated with the relative tagging power of the charge estimators, thereby allowing an optimal combination to be achieved in the whole phase space.The two types of input variables are described separately in Sections 6.1 and 6.2.

Charge estimators
Eight charge estimators are used to provide information on the initial charge state of b-quark hemispheres as inclusively as possible.These are: 1.The jet charge.This is the charge estimator used in the previous analysis [2]: the weighted sum of particle charges in a hemisphere, where the weights are the particle momenta along the thrust axis, p L , raised to the power κ.Four different values of κ are used for inputs to the neural net: 0.0, 0.5, 1.1 and 2.2.These focus on correspondingly higher values of track momenta but are of course highly correlated.Such charge estimators work well for all species of b-hadrons.
2. The secondary vertex charge.Using a topological vertexing algorithm combining information from all charged tracks in the hemisphere [18], a best estimate for a secondary vertex position is obtained in each hemisphere.An estimate of the charge of this vertex is calculated as: where w i , calculated by a dedicated neural net [18], is the probability that track i comes from the secondary vertex and q i is the track charge.This estimator provides a high quality charge tag for B ± hemispheres, and also helps to indicate which hemispheres are more likely to contain a neutral B meson.
to the B 0 have some correlation with the primary b-quark charge.Therefore a charge estimator, Q Pvtx , is calculated according to : 4. The weighted secondary vertex charge, Q Svtx , is calculated in a similar manner, replacing the weight (1 3 with the aim of improving the tagging for both charged and neutral B-decays via leading particle effects from the secondary decay.

5.
Fast kaon identification is formed by another dedicated sub-net [18] trained to identify charged kaons from b-hadron decays.The output for the most kaon-like particle in the hemisphere is signed by the charge of the particle and used as a charge estimator.

Control variables
The following variables are used as inputs to the neural net in order to provide it with some topological and b-hadron specific information: 1. | cos θ| is included since it is correlated with the quality of the secondary vertex reconstruction and hence with the tagging power of the estimators relying on the separation between tracks from the primary and secondary vertex.4. The spread of track separation weights i.e. the width of the w i factors distribution in a given hemisphere.This allows the net to de-weight those charge estimators which suffer from an ambiguous allocation of tracks to the primary and secondary vertices in cases of high charged multiplicities and/or poor vertexing.
For the asymmetry calculation, the difference between the neural net outputs of the two hemispheres is used.This is shown in Figure 2 together with the MC prediction.The charge separation is seen to be reasonably well simulated, but with a slightly larger width in data than in MC (there is also a small difference in the average value of the two distributions, but this is too small to be visible on the plot).

Event selection
The data set used for this analysis consists of four million hadronic Z decays recorded by ALEPH during the period 1991 to 1995 in a centre-of-mass energy range of M Z ± 3 GeV.Events are selected according to the standard ALEPH hadronic event selection based on charged track information [19].This selection has an efficiency of 97.5 % and the backgrounds from Z decays to τ + τ − and γγ interactions are estimated to be ∼ 0.3% each.These backgrounds are reduced to the 10 −5 level after the application of the subsequent cuts and can be safely neglected.
The average beam-spot position is determined every 75 events and used to constrain the event-by-event interaction point.In order to remove lifetime information from the tracks, these are projected onto the plane perpendicular to their parent jets (unless they point behind the interaction point found in a first iteration).Combining these projections with the beam-spot position fixes the interaction point to a precision of 50 × 10 × 60 µm 3 in horizontal, vertical and beam directions respectively.The thrust axis is determined in each event from all charged tracks and calorimeter clusters using the ALEPH energy flow algorithm [7].To ensure that a good fraction of the event is contained within the detector acceptance, the cosine of the thrust axis polar angle is required to be less than 0.95.In order to be able to evaluate the hemisphere b-tag and charge tag it is required that each hemisphere has at least one charged track and at least one jet with energy greater than 10 GeV.With these selection cuts a sample of 3,734,425 hadronic Z decays is retained for the remainder of the analysis.
In order to make use of double-tag methods for determining flavour purities, the b-tag range between zero and one is divided into eight sub-intervals as shown in Figure 3.The 672,111 events which are located in the region indicated by the black area in Figure 3, labelled E, are the ones accepted for the A FB analysis.The acceptance is less strict in the forward region which is not covered by the vertex detector and therefore has a broader distribution of b-tags.

Determination of flavour purities
The fractions of primary quark flavours in the selected event sample are measured in the data using the double tag method, that has been used previously in e.g. the ALEPH R b measurement [17,20].
In each of nine angular intervals, eight b-tag sub-intervals are defined as indicated in Figure 3, out of which only seven are mutually independent and the eighth is considered as "not tagged".These give rise to 28 doubly tagged event classes and 7 singly tagged classes: n ij and n i8 , which are related to the tagging efficiencies by: where ǫ f i is the probability for a primary quark of flavour f to be located in interval i, c f ij describes correlations between the b-tag values in the two hemispheres, R f is the Z → f f fraction of the Z hadronic decays and n had is the total number of events in the sample under consideration.The flavour index, f , runs over three flavour classes only (uds, c and b), hence efficiencies and correlation factors are averaged over the three light flavours.
Various assumptions are made to reduce the number of unknown parameters in Equation (9).The hemisphere correlations are taken from MC simulation and the values of R f are taken from previous measurements.Furthermore, it is necessary to fix the ratio of the two smallest efficiencies, where ǫ uds 1 ≈ 0.0008 and ǫ uds 2 ≈ 0.0013 according to the simulation.The ratio between these numbers is taken from simulation.
The remaining 20 unknown efficiencies ǫ f i are determined from a fit using Equations ( 9) leaving 15 remaining degrees of freedom in the fit.
The fitted efficiencies are combined to give event tag efficiencies by summing over the fiducial region labelled E in Figure 3: and from these efficiencies the three purities are constructed as: where the used R f fractions are shown in Table 6.The result of the fit is shown in Table 2.
The covariance matrix of the fitted efficiencies is propagated to the event-tag purities with diagonal elements also shown in Table 2 3 as obtained from the double tag fit.The values are given in units of 10 −4 .The second error is from MC statistics contributing through the assumed hemisphere tag correlations.
+45% and −80% for the uds-c, the uds-b and the c-b purity correlations, respectively.The χ 2 of the nine independent fits at different cos θ averages 0.99 per degree of freedom.
Already from Figure 1 it is seen that the efficiencies in the highest b-tag bins are higher in the data than in the MC.The fit determines that this enhancement is slightly larger for light and charm quarks than for b-quarks, resulting in a b purity which is about 0.5% lower than predicted by the MC in the selected sample.

Determination of quark charge separations
The mean values of the hemisphere charge separations, δ f , are determined from data in each bin of cos θ using the following procedure.
The data sample is subdivided into 14 bins with flavour compositions ranging from almost pure uds to almost pure b flavour.The sum of squares of the two hemisphere b-tags is used to define the bins and the flavour compositions of these new bins are derived from the fit to the data described in the previous section.
In each bin of this event b-tag variable, δ is measured according to Equation (3).A fit is performed using the 14 δ measurements and the right side of Equation ( 3) with the flavour specific values δb , δc and δuds as free parameters.Since the b-tags in the 14 bins bias the charge estimator differently, additional assumptions taken from MC are needed in order to compare the 14 δ measurements.For each flavour, the differences between the δ in a given bin and the δ in the E-region, β = δE − δbin , are held fixed at their MC values.Thus, for each flavour, the fit shifts δ by a fixed amount in all bins, minimising the following χ 2 : where i labels the b-tag bins and f the flavours (b, c and uds).The δf appearing in the χ 2 refers to MC events in the E-region and the correction factors r f are the three free parameters of the fit.
The inputs to the fit, averaged over all cos θ, are shown in Table 3 for b-and c-quarks together with the flavour combined output and the MC prediction.The fit reproduces accurately the b-tag dependence of δ.The angular dependence of the fitted parameters is shown in Figure 4 and, since the ratios between fitted and predicted values are consistent with being constant over angles, their average values are quoted below: r b = 0.995 ± 0.004 , r c = 1.031 ± 0.011 , r uds = 1.009 ± 0.009 .The average correlation coefficients of the error matrix are −64% between uds and c, 25% between uds and b and −49% between c and b.The average χ 2 probability is 0.72 over the independent fits in nine bins of polar angle.
The fitted values of δf are translated into the charge separations, δ f by Equation ( 2), where the charge correlation correction factors, (1 + k f ), are taken from MC.These factors vary with polar angle, averaging 1.086 for c-quarks and 1.092 for b-quarks.
In the fit of A b FB from Equation ( 6), the values for the charge separations, δ f (θ), are thus the MC values corrected by the r f (θ) factors.The ratios between individual light quark charge separations are taken from MC.In summary, the following values are used (averaging over polar angle and showing only the diagonal elements of the statistical error matrix):

Systematic errors
The determinations of the sample purity, the charge separation and the forward-backward asymmetry itself all rely to some extent on samples of Monte Carlo events.The features of the simulation most relevant for the analysis have been identified and varied within their errors in order to estimate the systematic uncertainty of the measurement as described in detail in the following sections.

The uds purity at large b-tags
In the fit yielding the flavour purities of the selected event sample it was assumed that ratio of the tiny efficiencies for selecting a light flavour with b-tag > 0.98 and for selecting a light flavour with 0.94 < b-tag < 0.98 was correctly predicted by the MC to be 0.62, averaged over polar angle.A systematic uncertainty on this ratio of ±30% is estimated by varying the track parameter smearing and the gluon splitting rates to heavy flavours in the MC.Propagating this uncertainty to the asymmetry yields a variation of ±0.00037.

Correlations between hemisphere b-tags
The correlation factors (c f ij in Equation ( 9)) are taken from Monte Carlo simulation.The influence of the correlations on the purities, via the fitted flavour efficiencies, is limited by Equation ( 11 In order to estimate the uncertainty in the correlation coefficients, the input parameters to the MC given in Table 1 are changed, one by one by the re-weighting technique, to the one σ higher values.Two additional re-weighted MC samples are studied in order to cover effects of QCD gluon radiation and of fragmentation.In one sample the angles separating the two selected b-jets are forced to agree with data and in another sample the charged multiplicity of the primary vertex is forced to agree with data.The largest effects on the b purities arise from changing the x b beam energy fraction, the charge particle multiplicities and the inter-jet angles, ranging from typically ±10 −5 in the central region to ±10 −4 at the most forward angles.The individual changes are added in quadrature to produce the total variation in the asymmetry of 0.00001.
There could be effects other than b-physics parameters causing an inaccurate simulation of the correlation factors.The possibilities include inaccuracies in the simulation of vertexing and of QCD effects.Therefore, a study of additional systematic errors is carried out along the lines of the ALEPH R b measurement [17,20].Here, the contribution to the correlation from a given variable v is calculated as: where f i (v) is the fraction of hemispheres, among all hemispheres with the variable situated between v and v + dv, having the tag i. Similarly g i (v) is the corresponding fraction with the opposite hemisphere tagged by tag i.The chosen variables, v, are the thrust polar and azimuthal angles, an estimator of the b-hadron momentum and the six input variables to the b-tagging neural network.Since the most important correlations are those involving b b events, the ten most significant of those are studied.In order to avoid having to subtract a very large non-b background in some bins, a soft b-tag is placed before evaluating the integrals.This soft b-tag leaves a 64% pure b-sample.
For the azimuthal angle, where no significant correlation is expected, the result of Equation ( 13) is indeed unity within a precision of 10 −3 .The rest of the variables do produce correlations, but the differences between data and MC are all between 10 −3 and 10 −2 and hence their contribution to the error on the measured asymmetry (via the flavour purities) is again very small, in total 0.00001.
Until now only correlations in b b events have been studied.Although the individual correlations for the lighter quarks have very little impact, the possibility of a combined effect of a large number of correlations has been investigated by simply ignoring all of them in the fit (except at the very most forward angle where geometrical effects are important).This results in a change in A b FB of +0.00002.A summary of all errors connected with the estimation of flavour purities is given in Table 4.

Statistical error in the purity determination ±0.00039
Light quark b-tag efficiencies ±0.00037 Hemisphere b-tag correlations ±0.00003Total systematic error due to purities ±0.00054 Table 4: Summary of systematic errors on the A b FB measurement from the determination of flavour purities.

Systematic uncertainties on charge separations
Various assumptions contribute to systematic uncertainties on the determination of the charge separations, δ f .In a first step, δb , δc and δuds are extracted from a fit to data, where the assumed flavour purities and charge bias corrections in the different b-tag bins contribute to the systematic error together with the statistical errors of the fit.In a second step, the constant of proportionality, (1+k f ), between δf and δ f contributes to the uncertainty.These systematic uncertainties are summarised in Table 5 and discussed in detail in the following.
• Flavour purities -This is evaluated by using the purities predicted by the MC instead of the fitted purities.The full changes in the fitted values of δf are propagated to the error on the asymmetry.
• Corrections for b-tag bias -Two overlapping procedures are used to estimate this uncertainty: 1.Each of the b-physics parameters listed in Table 1 are increased one at the time by one standard deviation, using re-weighted versions of the MC, and the differences in the fitted δ b and δ c are propagated to the asymmetry and entered as systematic errors in Table 5 .These estimates include the effect of the parameters on the charge correlation corrections discussed below.As in the case of hemisphere tag correlations, the effect of re-weighting the MC to agree with the observed primary vertex multiplicity and the observed distribution of inter-jet angles are also included.Finally, the V ector/(V ector + P seudoscalar) ratio in D-meson production is varied within its experimental error of 0.05 [21].
2. The values of δc and δuds are fixed at their MC values instead of being left free in the fit.In this case Equation (3) can be solved directly for δb in the E-region without any assumptions concerning bias corrections.The full difference between the result of this and the default result is taken as an independent contribution to the systematic error.This corresponds to an excursion in the fitted value of δc by almost two times the other uncertainties on this quantity.
• Charge correlation corrections -One contribution to the systematic error from the (1 + k f ) factors is found by varying the input parameters of the MC as described above.
Another contribution is estimated, like in the case of hemisphere b-tag correlations, from the difference between data and MC when building on the E-sample a "projection" of the hemisphere charge correlations along a series of observables, as follows: where Q ± (v) is the average positive or negative hemisphere charge, given the variable v in either the same or the opposite hemisphere, and P (v) is the probability density of v. Eight variables are chosen that are weakly correlated while covering the variables relevant for the charge determination.These are the thrust, the thrust axis direction, some inputs to the charge tagging neural net, the b-tag and the b-hadron momentum.
The differences in (1 + k v ) between data and MC are propagated to Table 5.
This procedure relates to the simulation of charge correlations for events in the selected sample consisting mostly of b b events.Previous studies using the jet-charge technique [5,3] have assessed the uncertainty in k c separately.Using a weighting power of κ = 0.5 the value k c = 0.085 ± 0.022 was found, where the error includes systematic uncertainties.This value is equal to the charm charge correlation in the present study and its error is propagated as an additional contribution to Table 5.
• Light quark charge separations -The light quark charge separations, δ u , δ d and δ s , are taken from Monte Carlo, but corrected by the overall correction factor from the fit.Systematic uncertainties are ascribed according to the results of [5] which represent relative uncertainties of 2.1%, 4.0% and 1.5% for u, d and s respectively.These errors are much larger than the error on the combined light quark δuds from the fit.However, when propagated through to A b FB , the systematic uncertainties arising from light quark charges are small, in total 0.00003.
• Track parameter smearing -The extra smearing and hit inefficiency mentioned in Section 5 also affects the charge tag.The effect of dropping these corrections to the MC altogether results in a change in the asymmetry of +0.00016 which is taken as an additional contribution to the systematic error.

Experimental systematic errors on Q FB
For the measuring apparatus itself to produce a fake Q FB , a forward-backward asymmetric bias in the charge measurement is needed.Such an effect could come from an asymmetry in effect.However, the effect on A b FB from dropping the correction is very small (0.00002) and this is added to the error.

Electroweak observables
Assumptions regarding electroweak observables other than A b FB , e.g.R b and A c FB , belong to a special class, since they are important for the Standard Model interpretation of the result.Their values used in this analysis are taken from the Standard Model and are listed in Table 6 together with their impacts on the measured value of A b FB .In case of the branching fractions R f , the variation in one flavour is compensated by changes in the light quark R f 's (f = uds) according to their relative magnitudes.

QCD correction
The radiation of gluons from the original quark pair can in principle lower the asymmetry.The effect is expected to be small with the present method, where the thrust-defined hemisphere charges are calibrated back to the charges of the primary quarks in these hemispheres using partly data and partly MC.However, this built-in correction may not be sufficient, since the hemispheres defined by the q q axis are not identical to those defined by the thrust.Therefore the asymmetry measured by repeating the analysis on the high statistics b b MC is compared with the generator level asymmetry in the same selected sample of MC events.The generator level asymmetry is found to be higher by a factor 1.0027 ± 0.0009.The most accurate calculation available of the full QCD correction in an unbiased sample of b b events results in a correction factor 1 + C QCD = 1.0354 ± 0.0063 [22,23].The ratio of C QCD to the same quantity calculated with the ALEPH MC is used to scale the residual QCD correction.Hereby the residual correction factor becomes 1.0034 ± 0.0019, where the error takes into account the statistical uncertainty and the fraction of the theoretical error which affects this analysis.
A summary of all the systematic errors is given in The selection efficiency fit which yields the flavour purities of the selected event sample has been checked by performing the fit on MC and verifying that the input efficiencies are indeed returned exactly.It has furthermore been checked that the distribution of b-tags in the MC agrees with that measured in the data after re-weighting the MC by the ratio of fitted to simulated flavour efficiencies.Finally, the overall b b fraction, R b , is left free in the fit.In order to get a stable fit, it is necessary to tie down some of the previously free parameters, and therefore the small bhemisphere efficiencies for b-tags less than 0.6 are taken from MC in this case.The result (for the three million preselected hadronic Z decays at peak energy) is R b = 0.2174±0.0008(stat).This measured value is consistent with the world average (0.21643±0.00073[23]).The change in the asymmetry from letting R b float in the fit is −0.0005, in agreement with the sensitivity shown in Table 6 and consistent with the uncertainty of the fit.

Lepton tagged events
For the δ fit, it has also been checked that it returns the input δf 's when applied to the MC.
The sample of events with an identified high-p T lepton offers an additional check on the assumptions entering the b-charge tag.This check is, however, limited by uncertainties in the lepton tagging due to c → l, b → c → l and B 0 → B0 transitions as well as lepton misidentification.These uncertainties are reduced to below the 10 −2 level by placing a high-p T cut of 1.5 GeV/c on the lepton sample.The distribution of the signed difference between the two hemisphere charges in such events is shown in Figure 5.The sign of the lepton provides the b b orientation in the event and allows a measurement of δ b = −0.3395± 0.0016, after subtracting the 3.8% non-b background in this sample.This is a factor 0.990 ± 0.006 smaller than the same quantity measured in the q q MC.The corresponding ratio between the δb 's was found to be 0.995 ± 0.004 in Section 9.This then provides a check on the MC prediction of k b : δb δMC b = 0.995 ± 0.007 (stat) .

Stability checks
In order to study the stability of the analysis the accepted region in the plane of the two b-tags is varied over a wide range.The results for peak data shown in Table 8 are consistent   8: Results on peak data for various event selections characterised by the sizes of the samples and their (predicted) b purities.The event selection shown in bold-face is the chosen one, the E-sample.
The measurements performed for each LEP I year in the data recorded at Z peak energies are shown in Table 9.They are seen to be consistent with being constant.

Figure 1 :
Figure 1: Distribution of the neural net hemisphere b-tags and the ratio between data and MC.The shaded region is the distribution for b-hemispheres in the MC.

2 . 3 .
The reconstructed b-hadron momentum is included because the relative accuracy of the various charge estimators depends on the b-hadron momentum.An estimator of this momentum is constructed from the jet closest to the line-of-flight of the b-hadron.The estimator is the sum of the charged track momenta in the jet, weighted by the probability that they come from the secondary vertex, and the projections of the neutral momenta onto the line-of-flight of the b-hadron.The missing energy in the hemisphere is also added.The reconstructed proper time of the b-hadron is used based on the reconstructed b-hadron momentum and the measured decay length.The intention here is to incorporate the increased probability of B 0 d mixing at long proper times.

Figure 2 :
Figure 2: The charge difference between forward and backward hemispheres measured in the selected event sample labelled E in Figure 3.The MC distribution corresponds to an asymmetry of A b FB = 0.0967.

1 EFigure 3 :
Figure 3: Upper plots: Subdivision into b-tag intervals with numbers above the plots enumerating the intervals.The shaded histogram is the expected b-tag distribution for b-hemispheres, shown separately for two regions of polar angle.Lower plots: The event acceptance, shown shaded and labelled E, in the plane of the two hemisphere b-tags.

Figure 4 :
Figure 4: The fitted values of δb and δc as a function of polar angle together with the values predicted by the MC.
) because the b-quark purity is high, 88% on the average.The correlation coefficient with the largest impact on the b purity is the negative correlation for observing both hemispheres in the highest b-tag bin: c b 11 − 1 = −0.02 in the central region of the detector.It is due to the pull on the shared reconstructed primary vertex exerted by a secondary vertex from b-hadron decay.If that coefficient is changed by an amount ±0.01, the b purity changes by the amount ±2.3 • 10 −4 in the central region.The other coefficients carry impacts on the flavour purities of size 10 −4 or less for a one percent change in the coefficient.

Figure 5 :
Figure 5: The difference between the two hemisphere charges, signed by the lepton charge, for events in the E sample containing a high-p T lepton.

8 Figure 6 :Figure 7 :
Figure 6: The forward-backward b b asymmetry as a function of thrust polar angle at three center-of-mass energies.The errors shown are statistical only.

Table 1 :
table.List of physics input parameters

Table 2 :
Flavour purities of the event sample labelled E in Figure and correlation coefficients being typically −90% ,

Table 3 :
From left to right: bin-by-bin flavour purities (P in percent), biases (β = δE − δbin ), measured, MC expected and fitted values of δbin .The numbers in this table are integrated over all polar angles.

Table 6 :
Summary of Standard Model inputs used in the measurement and their influence on the measured A b FB .

Table 7 :
Summary of systematic errors on the A b FB measurement.