JOLANTA BRODZICKA Institute of Nuclear Physics, Krakow November 05, 2003 Doubly charmed B decays B D(*)D(*)K ( for ~140 fb -1 )
Jan 23, 2016
JOLANTA BRODZICKA
Institute of Nuclear Physics, Krakow
November 05, 2003
Doubly charmed B decays
B D(*)D(*)K
( for ~140 fb-1 )
JOLANTA BRODZICKA November 05, 2003
B D(*) D(*) K
Plan
b c c s transition “wrong-sign” D production physics motivations analysis details preliminary results ( for ~140 fb-1 ) Dalitz-plot analysis summary
JOLANTA BRODZICKA November 05, 2003
B D(*) D(*) K
Mechanisms of BD(*)D(*)K decays
( I ) through external W-emission amplitudes
( II ) internal W-emission amplitudes (color-suppressed)
( III ) external +internal W-emission amplitudes
b cW - c c s + dd (uu) pair creation
B+ D(*)0 D(*)+ K0
B0 D(*)- D(*)0 K+
B+ D(*)+ D(*)- K+
B0 D(*)0 D(*)0 K0
B0 D(*)- D(*)+ K0
B+ D(*)0 D(*)0 K+
22 decay modes + c.c
JOLANTA BRODZICKA November 05, 2003
B D(*) D(*) K
Physics motivations
B D(*)D(*)K : good place to explore spectroscopy:
X DD cc-bar states above DD threshold scarcely known ((3770) (4040) …) molecular charmonia ( X(3872) ? ) , ccqq states, ccg hybrid states… Y D(*)K from W-vertex
c s : L= 0 0- Ds(1970) 1- D*s
(2112) well known
L= 1 jP = 1/2+ 0+ DsJ
± (2317) 1+ DsJ± (2457) seen, do not decay
to DK
( chiral doublet to Ds± Ds*± )
jP = 3/2+ 1+ Ds1± (2536) 2+ DsJ± (2573) not seen in B
decays ( do chiral partners exist?) charm deficit B D(*) D(*) K : significant (~1/2) contribution to b ccs transitions precise BF’s measurements needed
B0 D(*)+ D(*)- K0S to probe both sin2 and cos2
JOLANTA BRODZICKA November 05, 2003
B D(*) D(*) K
Analysis details at least 5 good tracks with abs(IP_dz)< 5cm abs(IP_dr)< 0.4cm R2< 0.3 K± : P(K/) > 0.4 ± : P(/K) > 0.1 electron veto: el_id < 0.95 K0
S abs( M(+ -) - MKs ) < 15 MeV only good_Ks
0 E > 50 MeV abs( M( ) -M0 ) < 15MeV
D±, D0 and D*±, D*0 reconstruction ( N(D)>1 : p(D) < 2. GeV in (4S) system )
B±, B0 reconstruction : all possible (22 + c.c) physical combinations D(*)D(*)K B vertex fit: with IP and B lifetime constraints Mbc > 5.2 GeV -0.40 < E < 0.35 GeV
JOLANTA BRODZICKA November 05, 2003
B D(*) D(*) K
D0
reconstruction
D0 K D0 K3 D0 K0
D0 Ks D0 KKD0 KsKK
D(*) plots for ~ 11fb-1 (exp17)
• vertex fit (refit 0 )
• M (DREC) cut: ± 20MeV
• (D0 K0 : - 50MeV)
• constraint_mass fit
BF ~ 28%
for skimmed sample (with DDK candidate )
(doesn`t applied for these plots)
JOLANTA BRODZICKA November 05, 2003
B D(*) D(*) K
D± reconstruction
D± K D± Ks
D± Ks0 D± KsK D± KK
BF ~ 15%
• reconstruction procedure like in D0 case
D(*) plots for ~ 11fb-1 (exp17)
D±
JOLANTA BRODZICKA November 05, 2003
B D(*) D(*) K
D*± reconstructionD(*) plots for ~ 11fb-1 (exp17)
D*± D0 ± BF ~ 68%
• vertex fit
(constraints :IP, B)
• abs(M(D*)-M(D)-mPDG)< 2.5MeV
( D*± D± 0 hopeless )
D0 KD0 K3 D0 K0
D0 Ks
D0 KK D0 KsKK
JOLANTA BRODZICKA November 05, 2003
B D(*) D(*) K
D*0 reconstructionD(*) plots for ~ 11fb-1 (exp17)
• vertex fit
(constraints :IP, B) refit 0 (vtx_constraint)• abs(M(D*)-M(D)-mPDG)< 5MeV
D*0 D0 0
BF ~ 62%
D0 K D0 Ks D0 K0
JOLANTA BRODZICKA November 05, 2003
B D(*) D(*) K
Multi-candidates events treatment
more than 1 candidates in the same B sub-mode• multi D(*) candidates per event • D(*)D(*) combinations with different K`s
~3. B cand per event ( all sub-modes )
D, D* probabilities (LR):
B(B(MMDD**)) LR_DLR_D** ((MMDD**) =) =
S(S(MMDD**))
S(S(MMDD**))++
S(MD), B(MD)parameterization from MD fitB(MB(MDD))
LR_D (MLR_D (MDD) ) =S(MS(MDD))
S(MS(MDD))++
S(MD*), B( MD*) parameterization
from MD* fit
JOLANTA BRODZICKA November 05, 2003
B D(*) D(*) K
Choice of the best B candidate
( for each B decay sub-mode separately )
LR_B = LR_DLR_B = LR_D(*) (*) ×× LR_D LR_D(*)(*)
• best B candidate : with max LR_B • equal LR_B case ( B`s differ only in K ) : larger K±_ID or better K0
S mass candidate chosen
B probability (LR_B)
S/(S+B) choice method “combines” both criteria:(M-MPDG)/ and S/B ratio
LR_B used also for background discrimination ( very efficient )
JOLANTA BRODZICKA November 05, 2003B D(*)D(*)K
B+ D0D0K+
S=94.4 ± 13.0
eff = ( 4.80 ± 0.14 ) *10-4
BF = ( 1 .30 ± 0.18 ± 0.21 ) * 10-3
E
N/7
.5M
eV
N/2
MeV
Mbc
LR > 0.04
Plot for Mbc>5.27 GeV
plot for abs(E)<25MeV
Fully reconstructed signal (from 2dim likelihood fit)
XNSIG_1 94.4 ± 12.9
XE0_1 -0.32E-02 ± 0.11E-02
XEW_1 0.59E-02 ± 0.07E-02
XM0_1 5.2815 ± 0.28E-03
XMW_1 0.25E-02 ± 0.04E-02
(external+internal W-emission)
JOLANTA BRODZICKA November 05, 2003B D(*)D(*)K
B0 D-D0K+
S= 127.6 ± 15.3
eff = ( 5.70 ± 0.14 ) *10-4
BF = ( 1 .75 ± 0.21 ± 0.29 ) * 10-3
E
N/7
.5M
eV
N/2
MeV
Mbc
LR > 0.01
Plot for Mbc>5.27 GeV
plot for abs(E)<25MeV
Fully reconstructed signal (from 2dim likelihood fit)
XNSIG_1 127.6 ± 15.3
XE0_1 -0.19E-02 ± 0.09E-02
XEW_1 0.55E-02 ± 0.07E-02
XM0_1 5.2812 ± 0.39E-03
XMW_1 0.26E-02 ± 0.03E-02
(external W-emission)
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
B0 D*-D0 K+ LR > 0.05
Fully reconstructed signal(from 2-dim Mbc vs.E likelihood fit)
XNSIG_1 86.9 ± 10.6
XE0_1 -0.19E-02 ± 0.09E-02
XEW_1 0.69E-02 ± 0.09E-02
XM0_1 5.2807 ± 0.41E-03
XMW_1 0.30E-02 ± 0.03E-02
eff = ( 2.84 ± 0.07 ) *10-4
BF = ( 2.02 ± 0.25 ± 0.34 ) * 10-3
Plot for Mbc>5.27GeV
Plot for abs(E)<25MeV
E
Mbc
N/7
.5M
eV
N/2
MeV
(external W-emission)
S = 86.9 ± 10.6
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
B+ D0D*0K+ ( right D*K flavor comb.) LR > 0.01 (ext+int)
Plot for Mbc>5.265GeV
Plot for abs(E)<45MeV
E
Mbc
Fully reconstructed signal (from 2-dim likelihood fit)
XNSIG_1 49.4 ± 11.6
XE0_1 0.104E-02 ± 0.43E-02
XEW_1 0.18E-01 ± 0.04E-01
XM0_1 5.2819 ± 0.98E-03
XMW_1 0.40E-02 ± 0.08E-02
S = 49.4 ± 11.6
Signif.= 7.0
eff = ( 0.64 ± 0.03) *10-4
BF = ( 5 .30 ± 0.93 ± 0.87) * 10-
3
N/7
.5M
eV
N/2
MeV
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
B+ D0D*+K0
s LR > 0.005
Fit results for directly reconstr. signal (from 2-dim Mbc vs.E likelihood fit)
XNSIG_1 40.2 ± 10.1
XE0_1 -0.13E-02 ± 0.25E-02
XEW_1 0.11E-01 ± 0.03E-01
XM0_1 5.2808 ± 0.80E-03
XMW_1 0.28E-02 ± 0.07E-02
S = 40.2 ± 10.1
Signif. = 7.5
Plot for Mbc>5.27GeV
Plot for abs(E)<25MeV
N/2
MeV
N/7
.5M
eV
E
Mbc
(ext)
BF = ( 1 .79 ± 0.45 ± 0.25) * 10-3
eff = ( 1.48 ± 0.06) *10-4
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
B0 D*-D*0K+
XNSIG_1 43.4 ± 10.1
XE0_1 -0.34E-02 ± 0.74E-02
XEW_1 0.31E-01 ± 0.07E-02
XM0_1 5.2814 ± 0.84E-03
XMW_1 0.32E-02 ± 0.07E-02E
Mbc
LR > 0.0 (ext)
Fitting results (2dim fit)
Plot for Mbc>5.27GeV
Plot for abs(E)<45MeV
S = 49.4 ± 11.6
Signif.= 7.1
eff = ( 0.49± 0.03) *10-4
BF = ( 5 .82 ± 1.35 ± 0.96) * 10-
3
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
B0 D*-D*+K0s
LR > 0.0 (ext+int)
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
Dalitz plot for
for signal-box events: abs(E)<25MeV Mbc>5.27GeV
M2 ( D0K+ ) right
M2 (
D0
B D
0 )
M ( D0K+ ) right
M ( D0B D0 ) M ( D0B K+ ) wrong
N/1
0M
eV
LR > 0.04
(3770)
B+ D0D0K+
Hatched histogram: Mbc sideband events abs(E+0.16)<50MeV Mbc<5.26GeVnormalized to backgr.in pion_lost signal box
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
M ( D0B D0 )
M ( D0K+ ) right
M ( D0B K+ ) wrong
mass distr. for B+ D0D0K+ + missing for lost signal-box events: abs(E+0.16)<50MeV Mbc>5.265GeV
N/1
0M
eV
LR > 0.1
Hatched histogram: Mbc sideband events abs(E+0.16)<50MeV Mbc<5.26GeVnormalized to backgr.in pion_lost signal box
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
B+ D0D0K+
N /
25M
eV
N /
50M
eV
N /
50M
eV
M ( D0K+ ) M ( D0BK+ )
M ( D0B D0 )
• fitting B signal in DD, DK invariant
mass bins inv. mass distributions without background
fitted B yield PhaseSpace distr. (3-body MC)
(3770)
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
Dalitz plot for
for signal-box events: abs(E)<25MeV Mbc>5.27GeV
M2 ( D0K+ )
M2 (
D0
D-)
M ( D0K+ )
M ( D0D- ) M ( D-K+ )
N/1
0M
eV
LR > 0.04B0 D-D0K+
Hatched histogram: Mbc sideband events abs(E+0.16)<50MeV Mbc<5.26GeVnormalized to backgr.in pion_lost signal box
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
B0 D-D0K+
SIG_1 15.3 ± 5.6
M0_1 2.577 ± 0.014
MW_1 0.025 ± 0.007
SIG_2 32.5 ± 7.1
M0_2 2.717 ± 0.008
MW_2 0.029 ± 0.006
N /
50M
eV
N /
50M
eV
M ( D0K+ )
M ( D0K+ )
DsJ± (2573) ?
DsJ chiral partner???
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
M2 (
D0
K+
)
M2 ( D*-D0 )
M ( D*-D0 )
M ( D0K+ )
M ( D*-K+ )
Dalitz plot for B0 D*-D0 K+ for signal-box events: abs(E)<25MeV Mbc>5.27GeV
N/1
0M
eV
LR > 0.005
?
Hatched histogram: Mbc sideband events abs(E)<80MeV Mbc<5.26GeVnormalized to backgr.in signal box
M ( D*-K+ )
also for LR > 0.01
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
M ( D*-D0 )
M ( D0K+ )
M ( D*-K+ )
mass distr. for + missing for lost signal-box events: abs(E+0.16)<50MeV Mbc>5.265GeV
N/1
0M
eV
LR > 0.05
Ds1 ? Hatched histogram: Mbc sideband events abs(E+0.16)<50MeV Mbc<5.26GeVnormalized to backgr.in pion_lost signal box
B0 D*-D0 K+
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
B0 D*-D0K+
N
/ 5
0M
eV
N /
50M
eV
N /
50M
eV
M ( D0K+ ) M ( D*-K+ )
M ( D*-D0 )
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
Dalitz plot for for signal-box events: abs(E)<25MeV Mbc>5.27GeV
M2
( D
*+K
0s )
M2 ( D*+D0B )
M ( D0K0s )M ( D*+D0B )
M ( D*+K0s )
Hatched histogram: Mbc sideband events abs(E)<80MeV Mbc<5.26GeVnormalized to backgroundin signal box
N/1
0M
eV
LR > 0.002B+ D0D*+K0s
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
M ( D*+D0B )
M ( D*+K0s )
M ( D0B K0s )
mass distr. for + missing for lost signal-box events: abs(E+0.16)<50MeV Mbc>5.265GeV
N/1
0M
eV
LR > 0.05
Ds1 ? Hatched histogram: Mbc sideband events abs(E+0.16)<50MeV Mbc<5.26GeVnormalized to backgr.in pion_lost signal box
B+ D0D*+K0s
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
B+ D0D*+K0s
N
/ 5
0M
eV
N /
50M
eV
N /
150M
eV
N /
100
MeV
M ( D0Ks )
M ( D*+K0s )
M ( D*+K0s )
M ( D*+D0B )
for signal lost box
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
Dalitz plot for B+ D0D*0K+ ( both D*0K flavor comb.)
for signal-box events: abs(E)<45MeV Mbc>5.265 GeV
M2 ( D*0 D0B )
M2 (
D*0
K+
)
M ( D*0K+ ) right
M ( D*0 D0B ) M ( D0K+ ) right
3.87GeVN/1
0M
eV
LR > 0.01
Hatched histogram: Mbc sideband events abs(E)<80MeV Mbc<5.26GeVnormalized to backgr.in signal box
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
M ( D*0D0B )
M ( D*0K+ ) right
M ( D0K+ ) right
mass distr. for B+ D0D*0K+ + missing for lost signal-box events: abs(E+0.16)<50MeV Mbc>5.265GeV
N/1
0M
eV
LR > 0.01
Hatched histogram: Mbc sideband events abs(E+0.16)<50MeV Mbc<5.26GeVnormalized to backgr.in pion_lost signal box
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
B+ D0D*0K+
N
/ 5
0M
eV
N /
50M
eV
N /
50M
eV
N /
50M
eV
M ( D0K+ )
M ( D0B D*0 )M ( D*0 K+ )
M ( D*0 K+ )
for signal lost box
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
Dalitz plot for LR > 0.0 for signal box: Mbc>5.27GeV abs(E)<45MeV
M2 ( D*-D*0 )
M2 (
D*0
K+
)
M ( D*-D*0 )
M ( D*0K+ )
M ( D*-K+ )
~10 events at the threshold at 2.54 GeV come from Ds1± (2536)-> D*0K± ?
B0 D*-D*0K+
JOLANTA BRODZICKA November 05, 2003
B D(*)D(*)K
B0 D*-D*0K+
N
/ 5
0M
eV
N /
50M
eV
N /
100
MeV
M ( D*0K+ )M ( D*-D*0 )
M ( D*-K+ )
JOLANTA BRODZICKA November 05, 2003
B D(*) D(*) K
LR_D*± vs. M plots
D0 K D0 K3
D0 K0D0 Ks
D0 KK D0 KsKK
LR_D0 vs. D_mass plots
JOLANTA BRODZICKA March 12, 2003B DDK
best B candidate choice method (see my December ICPV talk for more details)
B(MB(MDD))LR_D (MLR_D (MDD) ) =
S(MS(MDD))
S(MS(MDD))++D probability :
S(MD), B(MD)parameterization from MD fit
B probability for B± D0D0BK± : LR_B = LR_DLR_B = LR_D11 ××
LR_DLR_D22• best B candidate : with max LR_B
LR_DLR_D11 vs. vs. LR_DLR_D2 2 plots plots
(kpi)(kpi) (kpi)(no kpi)
choice method efficiency (for (kpi)(kpi) )
best = 82 %(estimated from rejected MC events DE, Mbc distributions and correlations between accepted and rejected B candidates )
for B± D0D0BK± (kpi)(kpi)reconstr. eff is 15% efficiency lost due to B_best choice method is 0.9%
JOLANTA BRODZICKA March 12, 2003B DDK
LR_B cut(good for bckg reduction and s/b improvement)
LR_D1 * LR_D2 cut
S /
sqrt
(S +
B )
LR_D1 * LR_D2 cut
S /
B
Signal MC: B± D0 D0B K±
Background: B± D0 D0B K± data sideband
( for BF = 1.5 * 10-3 )
JOLANTA BRODZICKA December 12, 2002
B D(*) D(*) K
LR_B cut efficiency LR_B cut: good for background reduction and s /b ratio
improvement
LR_D1 * LR_D2 cut LR_D1 * LR_D2 cut
Back
gro
und
reduct
ion (
DA
TA
)
Effi
cien
cy
reduct
ion (
sig
nal
MC
)
Signal MC: B± D0 D0B K± for all D0D0B combinations
Background: B± D0 D0B K± sideband
allall
JOLANTA BRODZICKA October 10, 2003B D(*)D(*)K
Possibility of Ds1+(2536) detection
Generated 20K MC decays:
B+ Ds1+(2536) D0B (Ds1)=5MeV
Ds1+(2536) D*0 K+ D0,D0B K
For directly reconstr. signal: B+ D*0 K+ D0B
M ( D*0K±)
N 140M0 2536 7.1
For partially reconstr. signal: B+ D0 K+ D0B + missing pion
M ( D0K±)
N 1420M0 2393 5.0
-143MeV shift
N/1
MeV
N/1
MeV
LR > 0.01
JOLANTA BRODZICKA October 10, 2003B D(*)D(*)K
Possibility of X(3872) detection (if any ) Generated 20K MC decays:
B+ X(3872) K+ (X)=5MeV Mmin=M(D*0+D0)X(3872) D*0D0B D0,D0B K
N/1
MeV
N/1
MeV
M ( D*0D0B) M ( D*0D0B)
N 1100 N 3400
For directly reconstr. signal: B+ D*0 D0B K+
For partially reconstr. signal: B+ D0 D0BK+ + missing pion
LR > 0.01