next up previous

Measurement of the Moments of the Hadronic Invariant Mass Distribution in Semileptonic B Decays
Brief Description of the Analysis -- CDF note 6973

The CDF Collaboration

May 18, 2004

Abstract:

Using 180 pbtex2html_wrap_inline768 of data collected with the CDFII detector at the Tevatron, we present a measurement of the first two moments of the hadronic invariant mass squared distribution in semileptonic B decays. By combining a direct measurement of the tex2html_wrap_inline772 piece with the D and tex2html_wrap_inline776 pieces taken from PDG, we find
displaymath756

displaymath757
where ``BR'' refers to the uncertainty coming from the branching ratios needed for the combination of the D, tex2html_wrap_inline776 and tex2html_wrap_inline772 pieces. tex2html_wrap_inline784 and tex2html_wrap_inline786 are used to determine the dominant parameters in the HQET expansion in the pole-mass scheme,
displaymath758

displaymath759
and in the 1S scheme
displaymath760

displaymath761
The results are compatible with previous determinations. The uncertainties are comparable or slightly better, although we have tried to produce a measurement with as little model dependence as possible.

Description

In order to constrain the length of the side opposite the angle tex2html_wrap_inline788 in the Cabibbo-Kobayashi-Maskawa (CKM) unitarity triangle, a precise measurement of tex2html_wrap_inline790 is needed. tex2html_wrap_inline790 is generally extracted from semileptonic B decays. Currently, the most precise method is based on the measurement of the inclusive semi-leptonic partial width into charm, tex2html_wrap_inline796. The Operator Product Expansion (OPE) applied to Heavy Quark Effective Theory (HQET) relates the experimental determination of tex2html_wrap_inline798 with tex2html_wrap_inline790. The relationship takes the form of an expansion in inverse powers of the B mass, tex2html_wrap_inline804. At each order in the expansion, new free non-perturbative parameters (vacuum expectation values of products of HQET operators with the right dimensionality) enter. There is one such free parameter (tex2html_wrap_inline806) at order tex2html_wrap_inline808, two (tex2html_wrap_inline810) at order tex2html_wrap_inline812, six at order tex2html_wrap_inline814, etc. In order to extract tex2html_wrap_inline790 from tex2html_wrap_inline798 some external information on these parameters is needed.

The same theoretical framework that predicts tex2html_wrap_inline798 predicts any weighted integral of tex2html_wrap_inline822, provided the weight is a smooth function of tex2html_wrap_inline824. By using as weight functions tex2html_wrap_inline826 and tex2html_wrap_inline828, with tex2html_wrap_inline830, the spin-averaged D mass, one defines the so-called first two moments of the hadronic mass distribution:
 eqnarray97
which are nothing but the (shifted) mean and variance of the tex2html_wrap_inline834 distribution in semileptonic charmed decays of B mesons. The moments are not sensitive to tex2html_wrap_inline790 but they are more sensitive to the non-perturbative parameters of HQET than tex2html_wrap_inline798 is. Therefore, measuring the moments provides a useful constraint on tex2html_wrap_inline840 that improves the overall precision on tex2html_wrap_inline790 as determined from tex2html_wrap_inline798. This is the goal of this analysis.

The tex2html_wrap_inline846 distribution in tex2html_wrap_inline848 decays can be split into three contributions corresponding to tex2html_wrap_inline850, where tex2html_wrap_inline852 stands for any neutral charmed state, resonant or not, other than tex2html_wrap_inline854, tex2html_wrap_inline856:
 equation126
where tex2html_wrap_inline798 is the inclusive tex2html_wrap_inline860 semileptonic width, tex2html_wrap_inline862 and tex2html_wrap_inline864 are the exclusive tex2html_wrap_inline860 partial widths to tex2html_wrap_inline868 and tex2html_wrap_inline870 respectively, and tex2html_wrap_inline872 is the (normalized) hadronic invariant mass distribution in the tex2html_wrap_inline852 channel. We will take tex2html_wrap_inline798, tex2html_wrap_inline862, tex2html_wrap_inline864, tex2html_wrap_inline882 and tex2html_wrap_inline884 from the Particle Data Group and concentrate on measuring tex2html_wrap_inline872. In this way, we only have to measure the invariant mass distribution for the tex2html_wrap_inline852 component without having to determine the tex2html_wrap_inline854, tex2html_wrap_inline856 components or the relative normalization between those and the tex2html_wrap_inline852 channel. The tex2html_wrap_inline852 spectrum is not well known, and includes, at least, two narrow and two wide states, together with a possible non-resonant tex2html_wrap_inline898 contribution. Its measurement is the main task of this analysis.

Only tex2html_wrap_inline900 decays (and charge conjugated) are reconstructed. Channels with neutrals are included using isospin relationships. Feed-down from one channel to another because of missing neutrals is subtracted statistically using the data themselves and isospin. It is assumed that the tex2html_wrap_inline898 decays of tex2html_wrap_inline860 saturate the difference between the inclusive and exclusive (to D and tex2html_wrap_inline776) semileptonic branching ratios.

tex2html_wrap_inline910 and tex2html_wrap_inline912 events are reconstructed using the following decay channels: tex2html_wrap_inline914, tex2html_wrap_inline916. Events are selected based on the recontructed tex2html_wrap_inline854 and tex2html_wrap_inline920 mass and the tex2html_wrap_inline922 mass difference. The tex2html_wrap_inline924 channel is reconstructed from the satellite peak in the tex2html_wrap_inline926 mass distribution. The tex2html_wrap_inline928 vertex (the B vertex) is reconstructed in 3D, and an additional pion, tex2html_wrap_inline932, is searched for, compatible with coming from the B vertex and incompatible with coming from the beam line.

Background subtraction is performed using mostly the data: side-bands are used to assess combinatorial background under the D mass and mass-difference peaks, and wrong-sing pion-lepton combinations (tex2html_wrap_inline938) parameterize prompt background to the tex2html_wrap_inline940 candidates. A small possible physics background, coming mostly from tex2html_wrap_inline942 decays with the tex2html_wrap_inline944 decaying semileptonically, is subtracted using Monte Carlo. Since only the shape of the mass distribution for the tex2html_wrap_inline772 component matters, the only efficiency corrections needed are those that can bias the mass distribution, together with the relative efficiency for the tex2html_wrap_inline948 and tex2html_wrap_inline920 components of the tex2html_wrap_inline772 piece. They are both obtained from a detailed, realistic Monte Carlo simulation after a small correction obtained using non-prompt tracks from data. In order to compare with theoretical predictions, the moments have to be measured with a well defined cut on tex2html_wrap_inline954, the lepton momentum in the B rest frame. While we cannot access this variable directly in data, our measurement has been corrected to a value tex2html_wrap_inline958 GeV/c using Monte Carlo.

The corrected tex2html_wrap_inline928 mass distribution is used to determine tex2html_wrap_inline962 and tex2html_wrap_inline964, the first and second moments of the tex2html_wrap_inline772 component of the mass-squared distribution, by simply computing the mean and variance of the corrected mass-squared distribution, without any assumption about the shape or rate of its several components.

The full moments of the hadronic mass-squared distribution, tex2html_wrap_inline784 and tex2html_wrap_inline786, are obtained by combining tex2html_wrap_inline962 and tex2html_wrap_inline964 with the D and tex2html_wrap_inline776 pieces, computed using values from PDG03. For the exclusive branching ratios to D and tex2html_wrap_inline776, all available information coming from charged and neutral B decays has been properly combined, using isospin invariance: the assumption is made that the isospin-related partial widths (not the branching ratios) of tex2html_wrap_inline986 and tex2html_wrap_inline988 are identical. Finally, the HQET parameters tex2html_wrap_inline806 and tex2html_wrap_inline992 are determined in the pole and 1S schemes, after applying constraints on the other HQET parameters coming from the known B and D hyperfine mass splittings.

Statistical uncertainties dominate the measurements of tex2html_wrap_inline962 and tex2html_wrap_inline964. Uncertainties in the inclusive and exclusive semileptonic branching ratios become important when combining tex2html_wrap_inline962 and tex2html_wrap_inline964 with the tex2html_wrap_inline854 and tex2html_wrap_inline856 pieces to obtain tex2html_wrap_inline784 and tex2html_wrap_inline786, the moments of the whole charm mass distribution. Finally, theoretical uncertanities become dominant in the extraction of the HQET parameters tex2html_wrap_inline806 and tex2html_wrap_inline992. Experimental systematics are subdominant in all cases.

Blessed Numbers

This section summarizes the numbers which have been blessed for this analysis.

Yields

Signal yields after the selection of tex2html_wrap_inline1020 and tex2html_wrap_inline912 combinations are given in Table 1.

  table178
Table: Yields after the selection of tex2html_wrap_inline928 events. The number of events in each case corresponds to the sum of both charge combinations, tex2html_wrap_inline1026. For the tex2html_wrap_inline1028, tex2html_wrap_inline1030 combined yield, a small (tex2html_wrap_inline1032) background from tex2html_wrap_inline1034 events has been subtracted.

Results

Systematic Uncertainties

Systematic uncertainties are tabulated in Table 2 and Table 3.

 table249
Table 2: Systematic uncertainties in the measurements of the tex2html_wrap_inline772 and full moments and in the extraction of tex2html_wrap_inline806 and tex2html_wrap_inline992 in the pole scheme. 

 table258
Table 3: Theoretical systematic uncertainties in the extraction of tex2html_wrap_inline1064 and tex2html_wrap_inline992 in the 1S scheme. All other systematics are as for the pole mass case. 

Blessed Figures

The following figures have been blessed for this analysis.

  figure267
Figure 1: Distributions of tex2html_wrap_inline1138 in tex2html_wrap_inline1140 combinations. Charge conjugated combinations are also included.

  figure274
Figure 4: Distributions of tex2html_wrap_inline1142 and tex2html_wrap_inline1144 in tex2html_wrap_inline1146 combinations. Charge conjugated combinations are also included.

  figure284
Figure: Distributions of tex2html_wrap_inline1142 and tex2html_wrap_inline1144 in tex2html_wrap_inline1152 combinations. Charge conjugated combinations are also included.

  figure294
Figure: Distributions of tex2html_wrap_inline1142 and tex2html_wrap_inline1144 in tex2html_wrap_inline1158 combinations. Charge conjugated combinations are also included.

 figure304
Figure 5: Raw invariant mass distribution for the tex2html_wrap_inline1160 channels. Charge conjugated combinations are also included. The mass region is limited at 3.5 GeVtex2html_wrap_inline1162 for illustration. 

 figure310
Figure 6: Raw invariant mass distribution for the tex2html_wrap_inline1164 channel. Charge conjugated combinations are also included. The mass region is limited at 3.5 GeVtex2html_wrap_inline1162 for illustration. 

  figure315
Figure 7: Fully corrected invariant mass distribution tex2html_wrap_inline1168. Charge conjugated combinations are also included. The number of events in each bin has been background subtracted and corrected for mass-dependent and tex2html_wrap_inline1170 relative efficiency corrections. The plotted errors take into account all corrections and subtractions

 figure323
Figure 8:   68% CL error ellipses for tex2html_wrap_inline806-tex2html_wrap_inline992 in the pole mass scheme including only statistical errors, statistical plus experimental systematics and total errors (including theoretical errors).

 figure327
Figure 10:   Comparison between the CDF measurement of tex2html_wrap_inline1176 and previous determinations. The HQET prediction has been constrained to the CDF measurements of tex2html_wrap_inline784 and tex2html_wrap_inline786. Note that points from the same experiment for different lepton energy cuts are highly correlated.

 figure333
Figure:   Comparison between the CDF measurement of tex2html_wrap_inline1182 and previous determinations. The HQET prediction has been constrained to the CDF measurements of tex2html_wrap_inline784 and tex2html_wrap_inline786. Note that points from the same experiment for different lepton energy cuts are highly correlated.

About this document ...

Measurement of the Moments of the Hadronic Invariant Mass Distribution in Semileptonic B Decays
Brief Description of the Analysis -- CDF note 6973

This document was generated using the LaTeX2HTML translator Version 96.1-h (September 30, 1996) Copyright © 1993, 1994, 1995, 1996, Nikos Drakos, Computer Based Learning Unit, University of Leeds.

The command line arguments were:
latex2html blessed_bhadr-moments.tex.

The translation was initiated by Marg Shapiro on Tue May 18 11:06:31 CDT 2004


next up previous

Marg Shapiro
Tue May 18 11:06:31 CDT 2004