350 likes | 710 Views
Dalitz plot analysis in D meson decays. Hai-Bo Li IHEP Workshop on Partial Wave Analysis and Dalitz Plot Analysis Beijing, Jan. 25-26. D 0 ! K B 0 ! K. Rich structures in D three body decays. Interference effects More important in D than in B decay
E N D
Dalitz plot analysis in D meson decays Hai-Bo Li IHEP Workshop on Partial Wave Analysis and Dalitz Plot Analysis Beijing, Jan. 25-26
D0! K B0! K Rich structures in D three body decays Interference effects More important in D than in B decay (usually!), since there are many resonance below D mass region (scalar)! Or rich final state interaction.
PDG2006 by D. Asner History of D Dalitz decays • Charm three body analyses from several experiments summarized in PDG2006 review by D. Asner. • Variety of Physics probed • Doubly-Cabibbo suppressed decays • CP Violation • charm mixing • Properties of light mesons • Properties pp & Kp S-wave • Usual analysis technique for P 3P decays is Dalitz-plot analysis technique
Dalitz-plot Analysis Technique(1) • Want to describe the internal dynamics of D0 ABC (P 3P) decay • Daughter 4-momenta: 12 parameters • Conservation of 4-momentum: 4 constraints • Masses of Decay products: 3 constraints • D is spin-0 : 3 orientations uninteresting • Decay described by 2 degrees of freedom • 3 Lorentz invariants (MAB)2, (MAC)2, (MBC)2 • Related by • (MD)2 + (MA)2 + (MB)2 + (MC)2 = (MAB)2 + (MAC)2 + (MBC)2 • Dalitz plot is (MAB)2 vs (MBC)2 • Phase space is “flat” in these variables
Dalitz-plot Analysis Technique(2) • Structure on Dalitz plot is due to • Non-uniform reconstruction efficiency • Background • Internal dynamics of the decay D ABC • Analyze structure on the Dalitz plot to elucidate a broad range of physics topics Phase space “triangular” for massless daughters becomes rounder as Q decreases BKSK+K-
BESIII Opportunity in Dalitz plot analysis • e+e- (3770) DD, 18 M D0D0 , 14M D+D- /year • 1. D mesons are produced almost at rest; clean events; • 2. Flavor tag • 3. CP tag • 4. 20 modes D P1P2P3 ( P= p, K) good for DP analysis • e+e- s=4170 DSD*S, 2 M DSD*S /year • DS Dalitz plot can also be done at BESIII. • Huge C0, C from (2S) and J/ Decays.
Technical Issues on Dalitz Plot analyses • Parameterization of the angular distributions for 3-body Dalitz plots. • Parameterization of the angular distributions for 4-body Dalitz plots. • Parameterization of the resonant line-shapes. • Parameterization of Non-Resonant contributions. • Normalization of the amplitudes (analytic? toyMC technique? ) • Efficiency. • Backgrounds. • Mass resolution. • Visualization of fit results (mass projections and angular distributions). • Estimators of the fit quality (chi^2, likelihood, adaptive binning). • How to tabulate Results (fractions, amplitudes and phases). • Common systematic errors evaluation. • Validations (use of toyMC, of simulated data).
Dalitz Model ---Parameterization of Dynamics • But: very often not just an isolated resonance • Mass dependence • Behaviour near decay-channel thresholds • Strongly overlapping resonances • and non-resonant dynamics • (e.g. crossed channel) By S. Spanier K matrix Isobar Model starting from 2-Body Scattering LASS parameterizations for S wave Relativistic Breit-Wigner form
a H b k H b Fk FH a c c Methodology? • We made no recommendation on the type of fit to make, but we do expect that the isobar model will still play a large role in many analyses. – Let’s start there • Recommend reading D. Asner’s review • hep-ex/0410014 v1, 5 Oct 2004 Isobar Model 2 d2G + ak e ifk / dm2dm2 k bc ab Isobars NR From D. Asner
Isobar Model • Each isobar term can be written as Mk = (-2pq)L YL0(cosqbc) TL(m2) x FH(q,rH) x Fk(m2) (for resonance J=L decay to 0+0) • TL(m2) = BWL(m2) = sin dk eidk • tan dk = M0G(m)/(M02-m2) • G(m) = G0 (p/p0)2L+1 (M0/m) Fk(p)2 • Form factor Fk usually taken as Blatt-Weisskopf damping factor. • A suggestion – define it to be unity at M0 for each resonance. • ! Easier comparison of complex coefficients ak .
k j • Best (only?) method for coupled channel analyses • Single pole in single channel -> BW • Distant, narrow poles -> sum of BW’s • One pole in >1 channel -> Flatte BW • Production vector allocates poles produced in three body decays • Recommend: D. Asner, hep-ex/0410014 v1, 5 Oct 2004 (again) • S U. Chung et al., Annalen Phys.4:404-430,1995 K Matrix Fits Klaus Peters talk in this meeting
K matrix parameterization For review see: Partial wave analysis in K matrix formalism, S.U. Chung et al., Annalen Phys.4:404-430,1995
Resonances in K-Matrix Formalism • The Breit-Wigner form is an approximation for isolated resonances far from • the opening of thresholds. • The K-matrix parameterization of multi-body interaction in terms of well • constrained 2-body interactions (isobar model) is a way to build on existing • amplitude measurements. • The K-matrix can be parameterized as an extension of the relativistic Breit • -Wigner form for s-channel resonance dominance. • The K-matrix is a good starting point to construct Lorentz invariant, 2-body • unitary, analytic (crossing) amplitudes. • Rem: • In some models even K-matrix poles are interpreted • as particles !?
Interest in K-+ s-wave—E791 hep-ex/0506040 B. Meadows for E791 ~138 % c2/d.o.f. = 2.7 Flat “NR” term does not give good description of data. Traditional iso-bar model
E791 D+!K-p+p+ PRL 89, 121801 hep-ex/0506040 B. Meadows for E791 ~89 % c2/d.o.f. = 0.73 (95 %) Probability Mk = 797 § 19 § 42 MeV/c2 Gk = 410 § 43 § 85 MeV/c2 The precise line-shape of K* is very important for this analysis. J/ KKp0 can definitely do that at BESIII.
E791 Model Independent PWA Compare MIPWA & LASS Data Compare MIPWA & BW Isobar Amplitude Phase Phase S-wave P-wave D-wave
E791 D+!p-p+p+ PRL, 86, 770 (2001) E791 collaboration
D+ p+p-p+ analysis from CLEOc hep-ex/0607069 CLEO collaboration 281 pb-1/1.8M D pairs KS removed Leading fractions shown on the projection plots.
81496 events with 91.5 fb-1(EPS05 paper) Much better fit quality than BW fit with only known resonances Similar quality to BW fit that included extra pp S-wave resonances Compare pp S-wave content CLEO BW Model - 15.1% BW fit with s1, s2– 30.4% A/S K-matrix fit - 16.2% AS param: K-matrix Analysis D0KSp+p- BaBar BaBar Data: K-matrix Model
Efficiency over the DP Branching fractions: B = N_signal/N_DD* Is the efficiency is constant across the DP? In general, no, there can be significant efficiency variation across the DP, from PID/tracking efficiency. Solution: Signal events and efficiency found bin by bin. However, we need to find the number of bins such that: 1. Not to few: efficiency variation are covered; 2. Not too many bins: Fewer statistics. Look at the variation of <> versus bin area. Efficiency variations from PID/tracking.
The momentum resolution effect Events does not respect dalitz kinematics boundary! Due to the momentumresolution, some signal events fall outside the Kinematics boundary. Since some signals are being cut away, the amplitude/Branching Fraction is being altered!
Rectangular Dalitz Plot Homogenization B0 p+p-p0 BaBar Easy to do efficiency correction, and bin by bin fit.
Efficiency and Resolution on B0p+p-p0 Dalitz Plot Fraction of mis-reconstructed events in Dalitz plane Large variations in the Dalitz plot corners Selection efficiency on Dalitz plane 00 Efficiency roughly flat Fraction of SCF events highly dependent on the Dalitz plot position +– 00 +– –+ –+ Truth-matched mass res.: double Gaussian core: 7÷8 MeV/c2 Truth-matched events Misreconstructed events SCF events: large variations in resolution over the DP 4D-convolution functions introduced G(m’rec, q’rec, m’true, q’true)
D0-D0 mixing from D0K+p-p0 Decay (1) • Event-level tagged • Prominent K* peak in DCS Mode Charm2006
D0-D0 mixing from D0K-p+p0 Decay (2) • Event-level tagged • Prominent r peak in CF Mode Charm2006
Resonant structures on the DP sensitive to D mixing • The resonance amplitudes are different for DCS and CF- there is more sensitivity to mixing • In D0 K-p+p0 the main resonance is K-r+ • In D0 K+p-p0 the main resonance is K*+ p- The ratio will depend on the point On the DP, full DP analysis will be very sensitive to mixing parameters, especially RM.
bc transition bu transition i δB -i γ e e Vub Vcb γ can be measured from the interference between decays with bcus and bucs transitions D0K0+- decay & CKM angle γ Interference occurs when some final state is accessible by both D0 and D0Giri-Grossman-Soffer-Zupan: PRD68, 054018 (2003):Final state = Ks0π+π- Dalitz Plot Analysis
BESIII Impact on g • Simultaneous CP tagged & flavor tagged analysis of D0Ks+- [correlated D’s,(3770)DD, (4140)DD(n)g(m)p0] • One can write • We will extract as well as in a model independent way. • This is exactly what the g analyses need. • CLEO-c has enough data in hand to reduce systematic error to 2o with 20 fb-1at BESIII.
CP-tagged D KSp+p-decay at BESIII CP=-1 CP=+1 Toy study from A. Bondar BESIII sensitivity A model independent way is proposed by A. Bondar and A. Poluektov hep-ph/0510246 50 ab-1 at super-B will allow a model independent g/f3 measurement at accuracy 2o. However, 10 fb-1 at y(3770/4170) data needed to accompany the analysis @BES-III: 10, 000, KSp+p- ,7,500 p+ p- p0,1,900 KSK+K- The d(cos(dD)) 2% 10 –20
Dalitz Fitter--Tools RooFit BaBar ……..many RooFit based fitters Charmfitter Miniut based, BaBar …. Charmfitter
Summary • Learn more experiences from BaBar, Belle, CLEO, FOCUS, E791 ….. Towards Uniformity Discussions with CLEO and BELLE • Methodology