3D Histo Fitting Study

Preliminary fitting results Kenji Hamano

Last modified : Mar 27, 2006 Back to Home Back to privious page

R16b skims Run3 was used.
Data luminosities:

	luminosity	# events	# files
OnPeak	30.630 fb-1	17527801/438622047	877
OffPeak	2.394	659867/31451100	33

MC luminosity

MC mode	cross section	# events	Luminosity	# files
BpBm (1235)	0.525e6 fb	23854196/86492000	164.75	1193
B0B0bar (1237)	0.525e6	17081965/61158000	116.49	855
ccbar (1005)	1.30e6	6127861/57700000	44.385	307
uds (998)	2.09e6	5639210/124040000	59.349	282
tautau (3429)	0.94e6	544465/69852000	74.311	28

Bhabha veto is applied.
KMicroNotPionGTL is applied.
FitDprob > 0.01 and FitBprob > 0.01 are applied.
Thrust cut |cosDlnonDl| < 0.88 is applied

Fitted result (Run3)

Branching fractions	SP5/SP6	Fitted raw number	New BR
B->Dlnu	0.0210	0.02237 +- 0.00036 (1.6 %)	0.0231
B->D*lnu	0.0560	0.05984 +- 0.00053 (0.89 %)	0.0618
B->D1lnu	0.0056	0.00293 +- 0.00038 (13 %)	0.0030
B->D0*lnu	0.0020	(0.00105)	0.0011
B->D1'lnu	0.0037	(0.00194)	0.0020
B->D2*lnu	0.0037	(0.00194)	0.0020
B->DPilnu	0.0060	0.00302 +- 0.00038 (12 %)	0.0031
B->D*Pilnu	0.0020	(0.00375)	0.0039
Sum	0.1000	0.09684	0.1000
Chi-square	999.499	471.991	na

	SP5/SP6	Fitted raw number	D0* only	D1 only	D1' only	D2* only
B->Dlnu	0.0210	0.02237	0.02163 (-3.3 %)	0.02276 (+1.7 %)	0.02288 (+2.3 %)	0.02146 (-4.1 %)
B->D*lnu	0.0560	0.05984	0.06078 (+1.6 %)	0.05997 (+0.2 %)	0.05982 (-0.03 %)	0.06018 (+0.6 %)
B->D0*lnu	0.0056	0.00293	0.00641	0.0	0.0	0.0
B->D1lnu	0.0020	(0.00105)	0.0	0.0634	0.0	0.0
B->D1'lnu	0.0037	(0.00194)	0.0	0.0	0.00572	0.0
B->D2*lnu	0.0037	(0.00194)/0.00786	0.0	0.0	0.0	0.00897
B->DPilnu	0.0060	0.00302	0.0401	0.0315	0.00344	0.00319
B->D*Pilnu	0.0020	(0.00375)	(0.00497)	(0.00391)	(0.00427)	(0.00396)
Sum	0.1000	0.09684	0.09789	0.09613	0.09613	0.09776

This fitting was performed with constraints:
Br(B->D*Pilnu) = Br(B->DPilnu)*67/54
This ratio 67/54 is based on the most recent BELLE measurement.

Covariance matrix and details are here

Float two components of D**lnu

1. D0* and D1+D1'+D2*
Details are here

2. D2* and D0*+D1+D1'
Details are here

3. D0*+D2* and D1+D1'
Details are here

Background dependence

Vary backgrond and fit

	SP5/SP6	Fitted raw number	New BR
B->Dlnu	0.0210	0.02287 +- 0.00039 (1.7 %)	0.0231
B->D*lnu	0.0560	0.06065 +- 0.00056 (0.92 %)	0.0613
B->D1lnu	0.0056	0.00284 +- 0.00041 (15 %)	0.0029
B->D0*lnu	0.0020	(0.00101)	0.0010
B->D1'lnu	0.0037	(0.00188)	0.0019
B->D2*lnu	0.0037	(0.00188)	0.0019
B->DPilnu	0.0060	0.00349 +- 0.00050 (14 %)	0.0035
B->D*Pilnu	0.0020	(0.00433)	0.0044
Sum	0.1000	0.09895	0.1000
Diff B	1	0.863 +- 0.044 (5.1 %)	na
Casc L	1	1.131 +- 0.050 (4.4 %)	na
L MisID	1	0.0 +- 0.144	na
Real D	1	1.008 +- 0.050 (5.0 %)	na

This fitting was performed with constraints:
      Br(B->D0*lnu) = Br(B->D1lnu)*20/56
      Br(B->D1'lnu) = Br(B->D1lnu)*37/56
      Br(B->D2*lnu) = Br(B->D1lnu)*37/56
      Br(B->D*Pilnu) = Br(B->DPilnu)*67/54
This ratio 67/54 is based on the most recent BELLE measurement.

Covariance matrix and details are here

Difference is +2.1 % for Dlnu and +1.3 % for D*lnu

Plots
D0KPi mode
Lepton momentum (D0KPi)
Lepton momentum : bin by bin (D0KPi)
D momentum (D0KPi)
D momentum : bin by bin (D0KPi)
cosBY (D0KPi)
cosBY : bin by bin (D0KPi)

DchKPiPi mode
Lepton momentum (DchKPiPi)
Lepton momentum : bin by bin (DchKPiPi)
D momentum (DchKPiPi)
D momentum : bin by bin (DchKPiPi)
cosBY (DchKPiPi)
cosBY : bin by bin (DchKPiPi)

f+-/f00, t+-/t00, D+/D0 dependence

Float (vary) the ratios f+-/f00, t+-/t00 and D+/D0
Result is here

D0KPi mode
Lepton momentum (D0KPi)
Lepton momentum : bin by bin (D0KPi)
D momentum (D0KPi)
D momentum : bin by bin (D0KPi)
cosBY (D0KPi)
cosBY : bin by bin (D0KPi)

DchKPiPi mode
Lepton momentum (DchKPiPi)
Lepton momentum : bin by bin (DchKPiPi)
D momentum (DchKPiPi)
D momentum : bin by bin (DchKPiPi)
cosBY (DchKPiPi)
cosBY : bin by bin (DchKPiPi)

Run differences

Varying 4 branching fractions
Chi squre is proportional to number of events in the first approximation, I normalised to Run3 using luminosity.
Run1 result
chisq = 293.148 -> 293.148*30.630/17.801 = 504.4
Dlnu = 0.02214 +- 2.3 % (-1.0 %)
D*lnu = 0.05919 +- 1.3 % (-1.1 %)
Run2 result
chisq = 668.177 -> 668.177*30.630/59.418 = 344.4
Dlnu = 0.02354 +- 1.2 % (+5.2 %)
D*lnu = 0.06164 +- 0.7 % (+3.0 %)
Run3 result
chisq = 471.991
Dlnu = 0.02237 +- 1.6 % (0 %)
D*lnu = 0.05984 +- 0.9 % (0 %)
Run4 result
chisq = 818.33 -> 818.33*30.630/99.357 = 252.3
Dlnu = 0.02265 +- 0.9 % (+1.25 %)
D*lnu = 0.05863 +- 0.5 % (-2.0 %)
RunAll result
chisq = 1710.97 -> 1710.97*30.630/207.206 = 252.9
Dlnu = 0.02251 +- 0.6 % (+0.6 %)
D*lnu = 0.05945 +- 0.4 % (-0.65 %)

Floating background etc...
Run1 result
chisq = 252.486
Dlnu = 0.02672 +- 6.4 % (+5.0 %)
D*lnu = 0.06628 +- 3.4 % (+4.6 %)
Run2 result
chisq = 538.283
Dlnu = 0.02622 +- 4.4 % (+3.0 %)
D*lnu = 0.06445 +- 2.6 % (+1.7 %)
Run3 result
chisq = 393.738
Dlnu = 0.02545 +- 5.7 % (0 %)
D*lnu = 0.06336 +- 3.4 % (0 %)
Run4 result
chisq = 684.572
Dlnu = 0.02533 +- 2.9 % (-0.5 %)
D*lnu = 0.06261 +- 2.1 % (-1.2 %)
RunAll result
chisq = 1345.18
Dlnu = 0.02534 +- 2.2 % (-0.4 %)
D*lnu = 0.06270 +- 1.2 % (-1.0 %)

Fitting method comparison (Run3)

Compare following settings
1. Variables : B_B1, B_B2, B_B32, B_B41
    Constraints: B_B31 = B_B32*20/56, B_B33 = B_B34 = B_B32*37/56, B_B42 = B_B41*67/54
2. Variables : B_B1, B_B2, B_B31, B_B32, B_B4
    Constraints: B_B33 = B_B34 = B_B32*37/56, B_B42 = B_B41*67/54
3. Variables : B_B1, B_B2, B_B31, B_B32, B_B4, Uncorr, CascL, FakeL, ContD
    Constraints: B_B33 = B_B34 = B_B32*37/56, B_B42 = B_B41*67/54
4. Variables : B_B1, B_B2, B_B31, B_B32, B_B4, Uncorr, CascL, FakeL, ContD, f_+0, c_D+0
    Constraints: B_B33 = B_B34 = B_B32*37/56, B_B42 = B_B41*67/54
        Gaussian constraint on f_+0 and c_D+0
4. Variables : B_B1, B_B2, B_B32, B_B4, Uncorr, CascL, FakeL, ContD, f_+0, c_D+0
    Constraints: B_B33 = B_B34 = B_B32*37/56, B_B42 = B_B41*67/54
        Gaussian constraint on f_+0 and c_D+0
t_+-/t_00 = 1.086 is used

	SP5/SP6	method 1	method 2	method 3	method 4	method 5	new BR
Chi-sq	936.735	445.493	443.242	404.455	401.336	406.595
B->Dlnu	0.0210	0.02372 +- 0.00038 (1.6 %)	0.02343 +- 0.00042 (1.8 %)	0.02404 +- 0.00045 (1.9 %)	0.02419 +-0.00048 (2.0 %)	0.02449 +- 0.00040 (1.6 %)	0.0236
B->D*lnu	0.0560	0.06227 +- 0.00056 (0.9 %)	0.06237 +- 0.00057 (0.9 %)	0.06310 +- 0.00060 (1.0 %)	0.06273 +- 0.00067 (1.1 %)	0.06271 +- 0.00067 (1.1 %)	0.0604
B->D0*lnu	0.0020	(0.00112)	0.00313 +- 0.00135 (43 %)	0.00364 +- 0.00137 (38 %)	0.00441 +- 0.00143 (32 %)	(0.00112)	0.0011
B->D1lnu	0.0056	0.00313 +- 0.00039 (12 %)	0.00260 +- 0.00053 (20 %)	0.00245 +- 0.00056 (23 %)	0.00228 +- 0.00056 (25 %)	0.00315 +- 0.00041 (13 %)	0.0030
B->D1'lnu	0.0037	(0.00207)	(0.00172)	(0.00162)	(0.00151)	(0.00208)	0.0020
B->D2*lnu	0.0037	(0.00207)	(0.00172)	(0.00162)	(0.00151)	(0.00208)	0.0020
B->DPilnu	0.0060	0.00321 +- 0.00039 (12 %)	0.00318 +- 0.00039 (12 %)	0.00373 +- 0.00052 (14 %)	0.00375 +- 0.00052 (14 %)	0.00367 +- 0.00052 (14 %)	0.0035
B->D*Pilnu	0.0020	(0.00398)	(0.00394)	(0.00463)	(0.00465)	(0.00455)	0.0044
Sum	0.1000	0.10157	0.10209	0.10483	0.10503	0.10385	0.100
Uncorr	1	1	1	0.878 +- 0.045 (5.1 %)	0.883 +- 0.044 (5.0 %)	0.877 +- 0.044 (5.0 %)
CascL	1	1	1	1.104 +- 0.052 (4.7 %)	1.094 +- 0.052 (4.8 %)	1.119 +- 0.051 (4.6 %)
FakeL	1	1	1	0.00 +- 0.12	0.00 +- 0.111	0.00 +- 0.123
ContD	1	1	1	0.973 +- 0.050 (5.1 %)	0.965 +- 0.049 (5.1 %)	0.968 +- 0.048 (5.0 %)
f_+0	1	1	1	1	1.0071 +- 0.0093 (0.9 %)	1.0061 +- 0.0075 (0.8 %)
c_D0+	1	1	1	1	0.9870 +- 0.0088 (0.9 %)	0.9931 +- 0.0074 (0.7 %)
detail	na	here	here	here	here	here

Plots of method 5

D0KPi mode
Lepton momentum (D0KPi)
Lepton momentum : bin by bin (D0KPi)
D momentum (D0KPi)
D momentum : bin by bin (D0KPi)
cosBY (D0KPi)
cosBY : bin by bin (D0KPi)

DchKPiPi mode
Lepton momentum (DchKPiPi)
Lepton momentum : bin by bin (DchKPiPi)
D momentum (DchKPiPi)
D momentum : bin by bin (DchKPiPi)
cosBY (DchKPiPi)
cosBY : bin by bin (DchKPiPi)

Dlnu FF slope

Result is here
This fitting was performed with constraints:
      Br(B->D0*lnu) = Br(B->D1lnu)*20/56
      Br(B->D1'lnu) = Br(B->D1lnu)*37/56
      Br(B->D2*lnu) = Br(B->D1lnu)*37/56
      Br(B->D*Pilnu) = Br(B->DPilnu)*67/54
This ratio 67/54 is based on the most recent BELLE measurement.
Gaussian constraint on
   f_+0
   c_D+0
t_0+ = 1/1.086;

Branching fractions	SP5/SP6	Fitted raw number	New BR
B->Dlnu	0.0210	0.08973 +- 0.00488 ( %)	0.0
B->D*lnu	0.0560	0.06066 +- 0.00100 ( %)	0.0
B->D1lnu	0.0056	0.00229 +- 0.00054 ( %)	0.00
B->D0*lnu	0.0020	0.00460 +- 0.00015 ( %)	0.00
B->D1'lnu	0.0037	(0.00)	0.00
B->D2*lnu	0.0037	(0.00)	0.00
B->DPilnu	0.0060	0.00406 +- 0.00048 ( %)	0.00
B->D*Pilnu	0.0020	(0.00)	0.00
Sum	0.1000	0.	0.1000
FF slope	0	1.145 +- 0.024 ( %)	na
Uncorr	1	0.914 +- 0.044 ( %)	na
Casc L	1	0.981 +- 0.059 ( %)	na
Fake L	1	0.198 +- 0.230 ( %)	na
Cont D	1	1.074 +- 0.050 ( %)	na
f_+0	1	1.016 +- 0.026 ( %)	na
c_D+0	1	0.952 +- 0.018 ( %)	na
Chi-square	1940.3	427.2	na

D0KPi mode
Lepton momentum (D0KPi)
Lepton momentum : bin by bin (D0KPi)
D momentum (D0KPi)
D momentum : bin by bin (D0KPi)
cosBY (D0KPi)
cosBY : bin by bin (D0KPi)

DchKPiPi mode
Lepton momentum (DchKPiPi)
Lepton momentum : bin by bin (DchKPiPi)
D momentum (DchKPiPi)
D momentum : bin by bin (DchKPiPi)
cosBY (DchKPiPi)
cosBY : bin by bin (DchKPiPi)

New parametrization of slope

Result is here
This fitting was performed with constraints:
      Br(B->D1'lnu) = Br(B->D1lnu)*37/56
      Br(B->D2*lnu) = Br(B->D1lnu)*37/56
      Br(B->D*Pilnu) = Br(B->DPilnu)*67/54
This ratio 67/54 is based on the most recent BELLE measurement.
Gaussian constraint on
   f_+0
   c_D+0
t_0+ = 1/1.086;

Branching fractions	SP5/SP6	Fitted raw number	New BR
B->Dlnu	0.0210	0.08973 +- 0.00488 ( %)	0.0
B->D*lnu	0.0560	0.06066 +- 0.00100 ( %)	0.0
B->D1lnu	0.0056	0.00229 +- 0.00054 ( %)	0.00
B->D0*lnu	0.0020	0.00460 +- 0.00015 ( %)	0.00
B->D1'lnu	0.0037	(0.00)	0.00
B->D2*lnu	0.0037	(0.00)	0.00
B->DPilnu	0.0060	0.00406 +- 0.00048 ( %)	0.00
B->D*Pilnu	0.0020	(0.00)	0.00
Sum	0.1000	0.	0.1000
FF slope	0	1.145 +- 0.024 ( %)	na
Uncorr	1	0.914 +- 0.044 ( %)	na
Casc L	1	0.981 +- 0.059 ( %)	na
Fake L	1	0.198 +- 0.230 ( %)	na
Cont D	1	1.074 +- 0.050 ( %)	na
f_+0	1	1.016 +- 0.026 ( %)	na
c_D+0	1	0.952 +- 0.018 ( %)	na
Chi-square	1940.3	427.2	na