Measurement of Vub using b hadron semileptonic decays

The magnitude of the CKM matrix element Vub is determined by measuring the inclusive charmless semileptonic branching fraction of beauty hadrons at OPAL based on b ->Xu l nu event topology and kinematics. This analysis uses OPAL data collected between 1991 and 1995, which correspond to about four million hadronic Z decays. We measure Br(b ->Xu l) to be (1.63 +/- 0.53 +0.55/-0.62) x 10^(-3). The first uncertainty is the statistical error and the second is the systematic error. From this analysis, Vub is determined to be: |Vub| = (4.00 +/- 0.65(stat) +0.67/-0.76(sys) +/- 0.19(HQE)) x 10^(-3). The last error represents the theoretical uncertainties related to the extraction of |Vub| from Br(b ->Xu l) using the Heavy Quark Expansion.


Introduction
The CKM matrix [1] describes the relation between quark weak and mass eigenstates, with the element V ub describing decays of the b to u quark.Its magnitude, |V ub |, can be calculated by measuring the inclusive b → u semileptonic decay rate.Given that the branching fraction of inclusive b → u semileptonic decay is of order 10 −3 , a large number of b hadrons are required to measure |V ub |.The dominant background to b → X u ℓν comes from b → X c ℓν decays because the branching ratio of b → X c ℓν is more than 50 times greater than that of b → X u ℓν.Here the lepton ℓ refers to either an electron or a muon, and b denotes all weakly decaying b hadrons (1) .X u and X c represent hadronic states resulting from a b quark semileptonic decay to a u or c quark respectively.The determination of |V ub | depends on the b to u and b to c semileptonic decay models.
The inclusive method developed by ARGUS [2] and CLEO [3] is to extract |V ub |/|V cb | from the excess of events in the 2.3 to 2.6 GeV/c region of the lepton momentum spectrum in the B meson rest frame, where the b → X c ℓν contributions vanish.This technique uses only a small fraction of the lepton phase space and so has considerable model dependence in extrapolating to the entire lepton spectrum in the B rest frame.In addition, since the LEP experiments can not precisely determine the B meson rest frame, this method is not appropriate for the LEP experiments.Instead, at LEP, |V ub | or |V ub |/|V cb | is extracted using a larger portion of the lepton spectrum as well as other kinematic variables.The inclusive measurement of the branching fraction of the b → X u ℓν decay has been performed at LEP by ALEPH [4], DELPHI [5] and L3 [6].
The theoretical uncertainty for the value of |V ub | extracted from a measurement of inclusive b → X u ℓν branching fraction differs from that extracted from measurements of exclusive b → u semileptonic decay rates.A recent theoretical study concludes that there is a 5% theoretical uncertainty on |V ub | values derived from b → X u ℓν inclusive measurements [7], using the Heavy Quark Expansion.There is a 15% theoretical uncertainty associated with |V ub | values extracted from measurements of the exclusive branching fractions B → πℓν or B → ρℓν [8], interpreted within the framework of the Heavy Quark Effective Theory (HQET).
In this paper, we describe the determination of |V ub | using the inclusive b → X u ℓν decay rate from the OPAL data taken at center of mass energies near the Z resonance.The event preselection, the b → X u ℓν decay models and the neural network used to separate b → X u ℓν from the background will be discussed in detail in the following sections.

The OPAL detector, data and Monte Carlo samples
The OPAL detector is a multi-purpose 4π spectrometer incorporating excellent charged and neutral particle detection capabilities.The OPAL detector is described in detail elsewhere [9].A brief description is given here.The central tracking system consists of a silicon microvertex detector, a vertex chamber, a jet chamber and z chambers.The momentum of tracks and the primary and secondary vertex position are reconstructed by the central tracking system, which is located inside a solenoid.The solenoid provides a magnetic field of 0.435T.Outside the solenoid is the electromagnetic calorimeter, which is composed of lead glass blocks and is used to measure the energies and positions of electrons and photons.The hadron calorimeter lies outside the electromagnetic calorimeter and is used to measure the energy of hadrons emerging from the electromagnetic calorimeter and assists in the identification of muons.The outermost OPAL detector is the muon detector which consists of a system of barrel and endcap muon chambers.A large fraction of muons with momenta less than 2 GeV/c are absorbed by the other detectors or the iron shielding before reaching the muon chambers.
The current analysis uses OPAL 1991 to 1995 data, collected near the Z resonance, comprising about four million hadronic Z decays.Monte Carlo simulated events were generated using the JETSET 7.4 [10] generator, with parameters described in [11].Approximately five million hadronic Z → b b decays were generated to study the b → X c ℓν decay and the b → c → ℓ cascade decay.Six million hadronic Z → qq (where q can be u, d, s, c and b) decays were generated to study the leptons from primary charm quarks and light quarks.Two hundred thousand events from a b → X u ℓν hybrid model [12] were produced to simulate the b → u semileptonic decay.The hybrid model will be described in detail in Section 3.1.

Signal and background simulation
The b to u semileptonic decay and background simulation are described below.The b to u semileptonic decay and background simulated events are passed through the full OPAL detector simulation [13] to produce the corresponding response.For this paper, the production fractions of B + , B 0 , B 0 s and Λ b in Z decays were adjusted to reproduce those given by the Particle Data Group [14].

The b → X u ℓν hybrid model
Several theoretical models have been proposed for the b → X u ℓν decay.Exclusive boundstate models [15][16][17][18] approximate the inclusive b → X u ℓν lepton spectra by summing contributions from all the exclusive final states.These exclusive models do not include all the possible final states nor any non-resonant states and therefore yield an incomplete prediction of the inclusive lepton momentum distribution, especially in the region of high hadronic invariant mass.The inclusive free quark models [19][20][21][22][23] treat the heavy quark as a free quark and the final state as a quark plus gluons.Free quark models are known to give poor agreement with experiments at low u quark recoil momentum.Therefore, a hybrid model [12] has been proposed to model the b → X u ℓν decay by using the exclusive model in the lower hadronic invariant mass region and using the inclusive model in the higher hadronic invariant mass region.The ISGW2 model [18] is used as the exclusive part of the hybrid model.The ACCMM model [19], combined with the W decay model [24] plus JETSET fragmentation, is used as the inclusive part of the hybrid model.Since the ISGW2 exclusive model includes the exclusive resonant final states 1S, 2S and 1P up to 1.5 GeV/c 2 in the hadronic mass, the boundary between the inclusive and exclusive parts of the hybrid model is placed at the hadronic invariant mass of 1.5 GeV/c 2 .The relative normalization of the inclusive and exclusive parts of the hybrid model is determined by the inclusive model.This hybrid model is only applied to decays of B mesons.There are no theoretical predictions for b to u semileptonic transitions of b baryons.The exclusive transitions of the b baryons in the OPAL tune of JETSET [10,11] are used.
In order to estimate systematic uncertainties due to modeling of the inclusive spectrum, alternative models are also studied.Signal events were generated with the QCD universal function [20][21][22] and parton [23] models.The invariant mass distributions of the hadronic recoil uq system are shown in Figure 1a for the QCD universal function, ACCMM and parton models.The invariant mass distribution of the hadronic recoil uq system for the hybrid model is shown in Figure 1b.

Background simulation
The ACCMM model [19] is used to describe the lepton spectrum of b → X c ℓν and b → c → ℓ decays.The fragmentation function of Peterson et al. [25] is used to describe the b quark and c quark fragmentation.The branching fractions of Xℓν were modified to reproduce those given by the Particle Data Group [14].The Λ b lepton momentum spectrum corresponding to -56% polarization [26] was used.
A hadronic event selection [27] and detector performance requirements are applied to the data.The thrust polar angle (2) | cos θ| is required to be less than 0.9 to ensure that the events are well contained within the acceptance of the detector.The selected events must pass the b identification, the lepton selection and the b semileptonic decay selection.All these selections are described in detail in the following sections.After all these preselections, the b → X u ℓν decay purity is 1.3% and the main background is from b → X c ℓν decays.

b identification
A neural network algorithm [28] based on charged particle vertex information is used to separate the b flavour events from the other flavour events in each hemisphere.If either hemisphere passes this neural network selection, the event is selected.After this neural network selection, the b purity is more than 91% and the b identification selection efficiency is approximately 30% per hemisphere from the Monte Carlo simulation in which a branching fraction of 1.0 × 10 −3 for the b → X u ℓν transition is incorporated.Both hemispheres are searched for electron and muon candidates after the b identification.

Lepton selection
Electrons are identified by a neural network [28] using the track and calorimeter information.The electron momentum is required to be greater than 2 GeV/c.Electrons from photon conversions, γ → e + e − , contribute a significant background to the prompt electron samples.Another neural network is used to reject this background [28].The photon conversion background is reduced by 94% after the photon conversion neural network selection, whilst retaining 98% of the selected prompt electrons.After all these requirements, the resulting electron efficiency is approximately 74% with a purity of 94% within the geometrical acceptance.
Muons are identified using reconstructed track segments in the muon chambers [28].The muon momentum is required to be greater than 3 GeV/c.The reconstructed tracks in the central detector are extrapolated to the muon chambers to see if they match the track segments reconstructed in the external muon chambers.The measured energy loss dE/dx is also required to be consistent with the expected value for a muon.After all these requirements, the muon selection efficiency is approximately 90% and the muon purity approximately 93% within the geometrical acceptance.
Electron and muon momenta transverse to the direction of the jet containing the lepton are required to be greater than 0.5 GeV/c in order to reject leptons from light quark decay.The lepton is included in the calculation of the jet direction.The jet finding is based on the cone algorithm [29].

b semileptonic decay selection
A neural network [30] based on lepton information is used to separate the b hadron semileptonic decays, b → X c ℓν and b → X u ℓν, from non b semileptonic decays.The distributions of the neural network output variable are shown in Figure 2.After this neural network b semileptonic decay selection, the b hadron semileptonic decay purity is 97% and the efficiency is 65% for this neural network; the c → ℓ events, where c is a primary quark, and b → c → ℓ events are suppressed.

b → X u ℓν neural network
Because of the dominant b → X c ℓν background, it is difficult to enrich the sample in b → X u ℓν decays using only one kinematic variable.A multi-layered feed-forward artificial neural network based on the JETNET 3.0 program [31] is used to enrich the sample in b → X u ℓν decays.There are four layers in this neural network.The neural network structure is 7-10-10-1.In the first layer, seven variables are used as inputs to the neural network.The last layer is the neural network output variable.A figure of merit [32] is used to determine the discrimination power of these seven variables in separating two classes of events, i.e. signal and background.The higher the figure of merit, the better the separation between the two classes.Over twenty kinematic variables were initially considered as inputs to the b → X u ℓν neural network.Only seven variables are selected as inputs to the b → X u ℓν neural network based on good separation between b → X u ℓν and background and good agreement between data and Monte Carlo simulated events.These seven input variables, in order of decreasing figure of merit, are: 1. the invariant mass of the most energetic final state particle combined with the lepton, 2. the lepton energy in the b hadron rest frame, where the b hadron energy and momentum are estimated using the techniques described in [33], 3. the lepton momentum transverse to the jet axis (the jet axis calculation includes the lepton), 4. the transverse momentum of the most energetic final state hadron with respect to the lepton direction (assuming all hadrons are pions), 5. the rapidity of the most energetic final state hadron calculated with respect to the lepton direction (assuming all hadrons are pions), 6. the fraction of the reconstructed b hadron energy carried by the lepton, 7. the reconstructed hadronic invariant mass, M x , which is calculated by: where i denotes all hadronic tracks and clusters.W i is the probability that the i th hadronic track or unassociated cluster comes from b decay and is calculated using the techniques described in [34].E i and p i are the energy and momentum of the i th hadronic track or neutral cluster.
Only the tracks and clusters from the same jet as the lepton are included in the calculation of these seven input variables.The seven input variable distributions for the b → X u ℓν and the background in the Monte Carlo simulation are shown in Figure 3.The agreement between the Monte Carlo simulated events and the OPAL data for these seven variables is shown in Figure 4.
Twelve thousand b → X u ℓν events, which were simulated with the hybrid model and have passed the event preselection, and the same number of background events from the multi-hadron Z → qq Monte Carlo simulation after the preselection are used to train the b → X u ℓν neural network.Two other samples of signal and background events of the same size are used to test the neural network performance.The neural network output distributions from b → X u ℓν and background are shown in Figure 5a.
The background composition from the b → X u ℓν neural network is shown in Figure 5b.Ninety percent of the background in this analysis comes from the b → X c ℓν decay, 6.8% from the b → c decay with the c subsequently decaying to a lepton.Another 0.6% comes from the c → ℓ decay in which the c quark is the primary quark.Other background processes make up the remaining 2.6%, of which 36% is from the b → τ decay with the τ subsequently decaying to an electron or a muon, and most of the rest of the background is from a pion or a kaon misidentified as an electron or a muon.

Monte Carlo
The branching fraction of b → X u ℓν decay can be obtained from the best fit of the Monte Carlo simulated events to OPAL data based on the b → X u ℓν neural network output distributions.Br(b → X u ℓν) is extracted from the b → X u ℓν neural network output distributions by minimizing: where N data k is the number of events from the data in the k th bin of the neural network output.N data is the total number of events in the data after preselection.The free parameter x is the fraction of signal events in the data after preselection, which can be converted to Br(b → X u ℓν) based on the number of signal events and the number of background events in the Monte Carlo simulation after preselection.f MC bu k is the fraction of simulated signal events in the k th bin of the b → X u ℓν neural network output with respect to the total number of simulated signal events after preselection.f is the fraction of simulated background events in the k th bin of the b → X u ℓν neural network output with respect to the total number of simulated background events after preselection.Here the background includes b → X c ℓν, b → c → ℓ, c → ℓ and other contributions.The sum over the index k is performed from the neural network cut to the last bin in the neural network output distribution.The Br(b → X u ℓν) from the fit result x, as well as its statistical and systematic errors, depends on the b → X u ℓν neural network cut.The resulting Br(b → X u ℓν) is stable, with variations less than 0.1 × 10 −3 , as the neural network cut varies in value from 0.3 to 0.7.A neural network cut of 0.7 is chosen to minimize the total relative errors and yields where the uncertainty is the statistical error only.
In Figure 6a, the neural network output distribution from data and the Monte Carlo simulation events with no b → X u ℓν semileptonic decay is shown and the excess of events in the data can be seen in the last bin.Here the distribution of Monte Carlo simulated events with no b → X u ℓν transitions is normalized to the same number of entries as the data.The χ 2 /ndf is 14.6/9, which corresponds to a 10% confidence level, assuming no contributions from b → X u ℓν transition.Here the χ 2 is calculated by summing over all bins in the neural network output distribution.When the b → X u ℓν transition is incorporated in the Monte Carlo simulation with a branching fraction of 1.63 × 10 −3 , the Monte Carlo simulation agrees much better with the data, as can be seen in Figure 6b.The χ 2 /ndf is then 8.3/8, corresponding to a 41% confidence level.
The data after subtracting the background from the Monte Carlo simulated events agree well with the simulated b → X u ℓν signal within statistical errors, which is shown in Figure 7.A series of cross checks, dividing the lepton samples into electron and muon samples and dividing the data into two samples for the years 1991 to 1993 and the years 1994 to 1995, are performed.The Br(b → X u ℓν) results are consistent within statistical errors for these cross checks.

Systematic errors
The list of systematic errors is given in Table 1.Unless otherwise specified, the systematic errors are estimated by varying each parameter described by one standard deviation and taking the corresponding largest errors.The resulting systematic errors in Table 1  (see text) ±0.19 b → X u ℓν modeling error (inclusive) Parton [23], QCD [20] ±0.14 b → X u ℓν modeling error (exclusive) ISGW2 [18]   within the experimental range 0.702 ± 0.008 recommended by the LEP Electroweak Working Group [35].This value is consistent with a recent determination of x E b = 0.714 ± 0.009 from SLD [37].The systematic error is also estimated from the Collins and Spiller fragmentation function [38] and Kartvelishvili fragmentation function [39].The uncertainties in c quark fragmentation can be neglected because the background from c → ℓ, where c is a primary quark, is very small.b → X c ℓν lepton momentum spectrum modeling: Different decay models are used to predict the lepton spectrum in the b hadron rest frame for the b → X c ℓν decay.Although all models are derived for B 0 and B + semileptonic decay only, they are extrapolated to the B s and Λ b semileptonic decay.This will be correct in the simple spectator model and is a reasonable approximation for this analysis.The lepton spectrum from the ACCMM model [19] is used as a base model for the b → X c ℓν decay.The systematic errors due to b → X c ℓν lepton momentum spectrum modeling are estimated from the ISGW [16] and ISGW** [36] models as prescribed by the LEP Electroweak Working Group [35].
The lepton spectrum from the b → c → ℓ decay in the ACCMM model is different from the lepton spectrum in the ISGW model.The systematic error due to the shape of the b → c → ℓ lepton spectrum is calculated and is found to be negligible.The lepton spectrum from the c → ℓ decay is varied from the ACCMM model to the ISGW model, where the c quark is a primary quark from Z decay.The systematic error is calculated and found to be negligible.2.
The Br( B → D * * ℓν) in Table 2 is obtained by averaging the Br( B → D * * ℓν) from ARGUS [40], ALEPH [41], DELPHI [42] and the total B semileptonic decay branching fraction subtracting the contribution from B to D and D * semileptonic decay, described by the LEP, CDF and SLD Heavy Flavour Working Group [43].For the decay of B → D * * ℓν, in which D * * refers to D 1 , D * 2 , D 2 and D * 1 , the branching ratio for each specific D * * final state is not well measured.For this analysis, the narrow final states of D * * in B → D * * ℓν are replaced by the broad states and then vice-versa to check the sensitivity of the Br(b → X u ℓν) to the relative ratio of the narrow and broad states of D * * in B → D * * ℓν.The effect on the Br(b → X u ℓν) is found to be negligible.
c hadron decay multiplicity: The systematic error of the Br(b → X u ℓν) associated with the c hadron decay charge multiplicity is estimated using the average charged track multiplicity of D + , D 0 , D + s decays as measured by MARK III [44].The systematic uncertainty of the Br(b → X u ℓν) is ±0.07×10 −3 from the uncertainty of c hadron decay multiplicity.
Λ b production rate: The PDG [14] gives the production fraction of B + , B 0 , B 0 s and Λ b in Z decay as (38.9 ± 1.3)%, (38.9 ± 1.3)%, (10.7 ± 1.4)% and (11.6 ± 2.0)% respectively.The neural network output variable distributions among B + , B 0 and B 0 s are similar and the systematic effects caused by the uncertainties of the production fraction of B + , B 0 and B 0 s are negligible.Due to the difference of the neural network output variable distributions between Λ b and B mesons, the fraction of Λ b is varied by one standard deviation to determine the corresponding systematic error.Λ b polarization: The lepton momentum spectrum from Λ b semileptonic decays depends on the degree of Λ b polarization.The systematic uncertainties are estimated by using the Λ b polarization ranging from -13% to -87%, as 95% [26] confidence level limits, which are converted to one standard deviation in Table 1.
Lepton identification efficiency: The number of selected events in the signal and background depends on the electron identification efficiency and the muon identification efficiency.The electron identification efficiency has been studied using control samples of electrons from e + e − → e + e − events and photon conversions, and is modeled to a precision of 4% [28].The muon identification efficiency has been studied by using the muon pairs produced in two-photon collisions and Z → µ + µ − events, giving an uncertainty of 2% [30].
Lepton fake rate: Fake electrons in the electron sample are primarily from charged hadrons (mainly charged pions) misidentified as electrons and from untagged photon conversions.The uncertainty associated with electron misidentification is ±21% [28].The muon fake rate is studied from K 0 s → π + π − and τ → 3π decay.The uncertainty of the fake muon rate is estimated to be ±8%.b → Xτ ντ branching ratio: One important composition in the "other" background in Figure 5b results from a b quark semileptonic decay to a τ lepton, with the τ lepton subsequently decaying to an electron or a muon.The branching ratio of b → Xτ ντ is (2.6 ± 0.4)% [14].The systematic error is estimated using the uncertainties of the b → Xτ ντ branching ratio.
Uncertainty of the b lifetime The average b hadron lifetime is measured to be (1.564 ± 0.014) ps [14].The uncertainty in b lifetime results in a negligible uncertainty in Br(b → X u ℓν).
Uncertainty of R b : The fraction of Z → b b events in hadronic Z decays, R b , is measured to be 0.2178 ± 0.0017 [14].The uncertainty in R b results in a negligible uncertainty in the background composition.
The first error 0.53 is the statistical error from the data only.The errors associated with the limited statistics of the Monte Carlo sample are included in the systematic error.This result is consistent with similar measurements from ALEPH, DELPHI and L3, the other three LEP experiments, shown in Table 3.In Table 3, the first error in Br(b → X u ℓν) combines the statistical error from the data and limited Monte Carlo statistics as well as the uncorrelated systematic uncertainties due to experimental systematic errors, such as detector resolution and lepton identification efficiency.The second error contains the systematic uncertainties from the b → X c ℓν Monte Carlo simulation models.The third error contains the systematic uncertainties from the b → X u ℓν models.The second and third errors are correlated between the various experiments.|V ub | can be obtained from Br(b → X u ℓν) [45,46] with inputs slightly revised as described by the LEP Heavy Flavour Working Group [43] in the context of the Heavy Quark Expansion [7]:

2 )Figure 1 :
Figure 1: a,b.The uq invariant mass distributions, a using the QCD universal function, ACCMM and parton inclusive models, b using the hybrid model.Only the portion of the uq invariant mass above 1.5 GeV/c 2 from the inclusive model in a is used in the hybrid model.The boundary between the exclusive model (left arrow) and the inclusive model (right arrow) in the hybrid model is indicated by the dashed line in b.

Figure 2 :
Figure 2: a,b.The b hadron semileptonic decay neural network output distributions.a for the b semileptonic decays and the scaled background from the Monte Carlo simulated events, here the background indicates all events excluding b → X c ℓν and b → X u ℓν , b comparison between the OPAL data and the Monte Carlo simulated events.The selected region is shown by the arrow in b.

Figure 3 : 1 - 7 .
Figure 3: 1-7.Comparison between the signal b → X u ℓν and the background in the Monte Carlo simulation for the seven b → X u ℓν neural network input variables.The b → X u ℓν signal and background are normalized to unity.The input variables in plots 1 to 7 are in the same order as the input variables defined in the text of Section 5.

Figure 4 : 1 - 7 .Figure 5 :
Figure 4: 1-7.Comparison between the OPAL data and the Monte Carlo simulated events for the seven neural network input variables.A branching fraction of 1.63 × 10 −3 for the b → X u ℓν transition is incorporated in the Monte Carlo simulation for comparison with the OPAL data.The input variables in plots 1 to 7 are in the same order as the input variables defined in the text of Section 5.

Figure 6 :
Figure 6: a,b.The neural network output distributions for data and Monte Carlo simulated events, a with no b → X u ℓν transition in the Monte Carlo simulated events, b with a branching fraction of 1.63 × 10 −3 b → X u ℓν decay incorporated.The distribution of Monte Carlo simulated events is normalized to the data for both plots.

Figure 7 :
Figure 7: a,b.The neural network output distributions.a Data after subtracting the background from the Monte Carlo simulated events (points) show agreement with the simulated b → X u ℓν signal (solid histogram).b Data after subtracting the Monte Carlo simulated events with a branching fraction of 1.63 × 10 −3 b → X u ℓν decay incorporated.Here the error bars include the statistical error from data and Monte Carlo simulated events.

b
and c hadron semileptonic decay branching fractions: The systematic error is estimated from the uncertainties of the branching fractions of B → Dℓν, B → D * ℓν, B → D * * ℓν and Λ b → Λ c Xℓν.There is a 6.8% background contribution from the b → c → ℓ decays and a 0.6% background contribution from the c → ℓ decays.The systematic error is also estimated from the uncertainties of the branching fractions of the b → c → ℓ decays.A summary of these systematic errors from the uncertainties of b hadron and c hadron semileptonic decay branching ratios is shown in Table

( 3 . 19 Table 2 :
04 ± 0.44)%[43] ±0.16 Br(b → c → ℓ) (8.4 +0.42 −0.39 )%[30] ∓0.02 Br(Λ b → Λ c Xℓν) (7.9 ±1.9)%[14] ∓0.06Total ±0.Systematic errors for Br(b → X u ℓν) from uncertainties of the b hadron and c hadron semileptonic decay branching ratios.MonteCarlo modeling errors: The systematic error for the Monte Carlo modeling errors is estimated by re-weighting each input variable distribution in the Monte Carlo simulation to agree with the corresponding data distributions.A branching fraction of 1.63 × 10 −3 for the b → X u ℓν transition is incorporated in the Monte Carlo simulation as shown in Figure 4.This gives a conservative estimation of the systematic uncertainty due to the modeling of the input variables.b → X u ℓν modeling error from the hybrid model: The boundary between the inclusive and exclusive regions in the hybrid model is varied from 1.5 GeV/c 2 to 0.9 GeV/c 2 .This conservatively estimates the systematic error arising from the placement of the boundary between the inclusive and exclusive models.This produces a uncertainty of ±0.19×10 −3 for Br(b → X u ℓν).b → X u ℓν inclusive model: The ACCMM model is the base inclusive model.The QCD universal function model and the parton model are used to estimate the systematic errors in the inclusive part of the b → X u ℓν hybrid model.This gives a change of -0.14 × 10 −3 for the QCD universal function model and +0.02 × 10 −3 for the parton model for the branching ratio of b → X u ℓν.The largest variation of Br(b → X u ℓν) from these models is taken as the systematic uncertainty of Br(b → X u ℓν) from the inclusive model.b → X u ℓν exclusive model: The ISGW2 model is the base exclusive model.The model implemented in the OPAL tune of JETSET [11] Monte Carlo simulation, which has the u quark and the spectator quark forming one single hadron in the final state, is used to estimate the systematic error in the exclusive part of the b → X u ℓν hybrid model.