811 likes | 954 Views
Hadronic B Decays To Double-Charm Final States. SERGIO GRANCAGNOLO L.Lanceri – J.P.Lees BINP Novosibirsk Particle Physics Seminar. Outline. Introduction The BaBar Detector at PEP-II The D sJ observations Theoretical Interpretations of D sJ Analysis of B D (*) D sJ decays
E N D
Hadronic B DecaysTo Double-Charm Final States SERGIO GRANCAGNOLO L.Lanceri – J.P.Lees BINP Novosibirsk Particle Physics Seminar
Outline • Introduction • The BaBar Detector at PEP-II • The DsJ observations • Theoretical Interpretations of DsJ • Analysis of BD(*)DsJ decays • Results: branching fractions and angular distributions • Comparison with models and conclusions Sergio Grancagnolo
The Standard Model • Fundamental particles: • 6 quark , 6 leptons • 4 interactions • The model works well but there are several issues to be understood, for instance: • Higgs boson • Supersymmetry • Strong interactions W,Zbosons Sergio Grancagnolo
Quantum Numbers Of The Quarks Quark Property Sergio Grancagnolo
CKM Matrix and Unitary Triangle qi=u,c,t Unitary relationship CKM W+ Vij qj=d,s,b A complex phase in the V matrix can be a source of CP violation in B decays VV†=I VudVub*+VcdVcb*+VtdVtb*=0 a VtdVtb* VudVub* Unitary triangle g b VcdVcb* Sergio Grancagnolo
_ _ _ _ _ Mesons in the Quark Model • Quarks exist only in baryons and mesons • Mesons are made of a quark-antiquark pair • As an example: • Mesons are not stable • Mass, charge and lifetime are main characteristics • Meson width~ 1/lifetime depends on the allowed decay modes Sergio Grancagnolo
sQ sq q _ Q Heavy Quark Approximation ℓ In the heavy quark approximation mq<<mQ,, mQ sQ, j conserved However J, P good quantum numbers Sergio Grancagnolo
Charmed Mesons Spectroscopy _ _ • States with ℓ=1 can decay strongly with emission of a pseudoscalar meson • j=1/2 emission in s-wave • j=3/2 emission in d-wave • D*0,D´1observed by CLEO, Focus and Belle • Broad resonances as expected ℓ=0 ℓ=1 broad ~100 MeV narrow ~10 MeV Sergio Grancagnolo
_ The expected cs Meson Spectra M.Di Pierro, E.Eichten Phys. Rev. D64, 114004 (2001) 2.51 GeV 2.36 GeV States expected but not observed • Masses over threshold DK(*) • Broad states (large widths) * Sergio Grancagnolo
B Meson Decay • Spectator quark model the other u,d quark enters the final state without participating to the interaction • In hadronic decays, could be tested the factorization hypothesis: the final hadrons are produced independently Since mb >>mu,d the B meson decay dominantly through the disintegration of the b quark. The main transition W* ℓn semileptonic is the weak decay bcW* where _ hadronic W* qiqj W* virtual boson Sergio Grancagnolo
_ - _ _ _ _ Exclusive Hadronic B decays • In exclusive decays all particles in final state are reconstructed • Double charm decays contains two mesons with charm quarks • Examples: Ds- BDsD _ D(*)0,D(*)+ B-,B0 D(*)0 K(*)- B DDK B-,B0 D(*)0,D(*)+ Sergio Grancagnolo
The PEP-II B-factory at SLAC PEP-II is a high luminosity, asymmetric, e+e- collider Integrated luminosity Lint=254 fb-1 113fb-1 Ldesign = 3 x 1033 cm-2s-1 Lpeak = 9.21 x 1033 cm-2s-1 Sergio Grancagnolo year
B-factory Cross Sections E(e+) = 3.1 GeV E(e-) = 9.0 GeV The boost allows a separation of the two B vertices. boost: bg=0.56 Ecm=10.58 GeV _ _ U(4S)BB _ _ _ s[e+e- hadrons](nb) _ √s(GeV) _ _ e+e- bb on-resonance BB _ “coontinuum”e+e- cchigh momentum charmed particles Sergio Grancagnolo
BABAR Detector 1.5 T solenoid Electromagnetic Calorimeter e+ (3.1 GeV) Cerenkov Detector (DIRC) e-(9 GeV) Drift Chamber Silicon Vertex Tracker Instrumented Flux Return Sergio Grancagnolo
DsJ(2317) Discovery * + • BaBar discovered a new particle decaying into Dsp0 • c and s quarks • Mass < DK threshold • Width < 10 MeV • Seen by Belle and CLEO • Is this the expected Ds0? BaBar collaboration Phys.Rev.Lett. 90, 242001 (2003) + _ Dsp0Invariant mass + *+ m=2.317GeV GeV Inclusive selection of high momentum charmed meson from coontinuume+e- cc _ Sergio Grancagnolo
DsJ(2460) Discovery + • CLEO observed another state decaying to Ds p0! • c and s quarks • Mass < (DK)* threshold • Width < 10 MeV • Observed also decay modes: • Dsg, Dsp+p- • Is this the expected Ds1? CLEO collaboration Phys. Rev. D68, 032002 (2003) *+ _ Ds p0Invariant mass *+ 80 60 m=2.460 GeV Events/7 MeV/c2 40 20 0 + + 2.25 2.5 2.75 GeV + Seen by BaBar and Belle Sergio Grancagnolo
_ The Observed cs Meson Spectra 2.51 GeV 2.36 GeV New states observed • Masses below threshold DK(*) • Narrow states * Sergio Grancagnolo
Isospin Violation in These Decays • Isospin symmetry is not exact • Violation already observed in Ds* Dsp0 decay _ _ _ _ _ _ _ DsJ Dsp0 Invoked hp oscillation P.L.Cho, M.B.Wise Phys.Rev.D49: 6228-6231,1994 Sergio Grancagnolo
Theoretical Interpretations of DsJ Standard interpretations Exotic interpretations
Standard interpretations Entia non sunt multiplicanda praeter necessitatem (G.Occam) • Quark models • Potential: coulombian • (0-,1-),(0+,1+) chiral partners • doublets mass splitting via chiral symmetry breaking • transitions via scalar meson + linear Cahn, Jackson + spherical not linear Lucha, Schoberl need to adjust a posteriori input parameters, predict mass higher than observed or not reproduce non-strange charmed mesons spectra Bardeen, Eichten, Hill hyperfine splitting for charmed mesons (D, D*, etc.) marginally compatible with experiments Sergio Grancagnolo
Standard interpretations • Unitarized chiral models • generalization replacing a light quark with an heavy quark • Non-perturbative methods • lattice QCD • QCD sum rules Beveren, Rupp several new mesons predicted not observed Bali initial difficulties to reproduces masses, reproduces mass splitting Dai, Huang, Liu, Zhu low accuracy Sergio Grancagnolo
D K Ds p Exotic Interpretations Barnes, Close, Lipkin Dsp molecule cs DK 4-q mixing DK molecule Szczepaniak Browder, Pakvasa, Petrov D _ _ _ _ qq qqqq K di-quark pairs _ qq Maiani, Piccinini, Polosa, Riquer _ qq Sergio Grancagnolo
Analysis of BD(*)DsJ decays Branching ratios: Method Event selection Signal and Backgrounds Efficiency and “cross-feed”
BD(*)DsJ Decays • ExclusiveDsJ production: expected to be dominant • Allow to measure DsJquantum numbers • In principle, allow to discriminate between conventional and multi-quark scenarios compared with other B decays such as BD(*)Dsand BD(*)D(*)K • If the DsJ is the conventional cs state should be produced in the following graph: _ _ Weak external W emission DsJ- - _ _ _ _ _ D(*)0,D(*)+ B-,B0 Same graph as BD(*)Ds similar branching ratios could be expected Sergio Grancagnolo
BD(*)DsJ Decays (II) • We search for DsJ particles looking at the 12 combinations: • With DsJ decays: • We measure branching ratios, quantum numbers JP Sergio Grancagnolo
Subdecay Modes Intermediate particles are reconstructed in the following modes: Green::clean modes Total: 60 different submodes combined to give the 12 combinations Sergio Grancagnolo
Analysis Goal and Method • We aim to measure branching ratiosBri (i=1…12) of the exclusive double charm two body production of DsJ(2317)+and DsJ(2460)+in B0and B+ • nisignumber of signal candidates for mode i • after combinatorial background subtraction • nixfdnumber of crossfeed events for mode i • contains background from other signal modes • eireconstruction efficiency from simulation • NBB= [122.0 ± 0.6(stat) ± 1.3(syst)] 106 (113 fb-1) * _ Sergio Grancagnolo
_ D0 Ds + DsJ(2460)+ * A specific example: B0D*-DsJ(2460)+ * • Reconstruct the chain: • Reconstruct tracks (K,p) and photons (g) • Select D0, Ds , f, p0computing invariant masses • Use beam energykinematicconstraint • Fit nisig in Dsg invariant mass distribution K+ K+ p- f D*- p- K- p+ B0 g Sergio Grancagnolo
Event Selection: Invariant Masses Invariant mass: D0 Kp Dsfp D* D0p fKK 40000 20000 0 0.99 1.02 1.04 m(f)(GeV/c2) Particles masses are set to their nominal values (mass constraint) Sergio Grancagnolo
Event Selection: B candidates • Compute p*B and E*Bfrom selected D*, Ds, g • Use the B-factory constraint E*beam to compute: 5.272<mES<5.288 GeV mES “Signal box”: |DE|<32MeV Use of beam kinematic variables better resolution ΔE uncorrelation Sidebandsto estimate background outside signal box Sergio Grancagnolo
DE resolution • Same resolution for all the submodes • A systematic error will take in account differences between data and simulation Simulation of signal events Data candidates in mES signal region s(DE)=16.1 s(DE)=18.9 Cross-hatched background from sidebands Missing energy effect Sergio Grancagnolo
DE resolution (II) Final values used in selection (MeV) Better resolution for modes with a p0 (mass constraint) Sergio Grancagnolo
Background Rejection • Reduction of the combinatorial background • Simulated signal events selected in signal region • Background from data events selected in DsJ mass sideband region • Curves represent fraction of events cut by m(D0g)> mcut(D0g) • Optimal cut set at the maximum separation between two samples m(D*g) cut Events rejected: 25% signal 75% backgrd m(D*g)>2.4GeV/c2 Gev/c2 Sergio Grancagnolo
Optimization • Maximized the significance ratio: S = simulated signal events in signal region B = background from data in m(DsJ) sidebands Tried different cut levels for D and Dsusing PID, vertexing and helicity cut f cos(qhel) f mass 5000 40000 2500 20000 Tried different numbers of s cut for variables: DE, m(Ds), m(D) -1 1 1.94 2.0 cos(qhel) m(f) Cleaner modes require less stringent cuts Sergio Grancagnolo
Fit nisig in DsJ(2460)+Ds+g m(Dsg) • Finally, in selected candidates: m(Dsg) • Fit the background shape with a polynomial • Fit the signal peak with a Gaussianof fixed width • s=12 MeV • estimated in data • Events in the signal peak: Entries/10 Mev/c2 GeV/c2 significance=11.7 nisig = 53.0±7.7 Sergio Grancagnolo
Efficiency and Cross-feed • From gi=60k simulated signal events for each mode i • Efficiency: nisim= number of B0D*-DsJ(2460)+ events reconstructed in the corresponding simulated sample • Total cross-feed: nijsim= number of B0D*-DsJ(2460)+ events reconstructed in the simulated sample (mode j) fij= cross-feed from the mode j to the mode i Typical efficiency range: 1-10% depending on the presence of photons, soft tracks, stringent cuts, etc. ; Sergio Grancagnolo
Narrow Cross-feed Efficiency Reconstructed mode: B0D*0Ds1- [Ds-g] m(DsJ) Generated mode: B0D*0Ds1- [Ds-g] nisim= 2778 gi=60000 ei=(4.63±0.08)% Cross-feed Generated mode: B0D*+Ds1- [Ds+g] nijsim= 24 gj=60000 fij=(0.82±0.04)% Narrow: sxfdssig Sergio Grancagnolo GeV/c2
Wide Cross-feed Efficiency Reconstructed mode: B0D0Ds1- [Ds-g] m(DsJ) Generated mode: B0D0Ds1- [Ds*-p0] nisim= 1350 gi=60000 Cross-feed ei=(2.25±0.07)% Generated mode: B0D*+Ds0- [Ds-p0] nisim= 144 gi=60000 fij=(0.24±0.02)% Cross-feed nisim= 162 gi=60000 Generated mode: B0D*0Ds0- [Ds-p0] fij=(0.27±0.02)% Wide: sxfd 2.5 ssig Sergio Grancagnolo GeV/c2
Branching Ratios and Cross-feed An iterative procedure is needed: • Compute for each mode i without considering cross-feed • Estimate nixfdusing Brj and the cross-feed fij from all the modes • Subtract the number of cross-feed events • Compute the corrected branching ratio • Recompute the cross-feed iterating point 2-4 until convergence. __ Sergio Grancagnolo
Fit Results And Significance s=5.5 s=4.2 s=5.0 s=5.2 s=7.4 s=11.7 s=3.1 s=5.1 s=4.3 s=6.0 s=7.7 s=2.5 Sergio Grancagnolo
Main Systematic Errors Depends on the tracks or photons number • Tracking efficiency 9% • g/p0 efficiency 5% • Background fitting model5% • Tried exponential instead of polynomial to fit background • DE width 5% • Changed the width of the DE signal region by ±3 MeV • DsJ width 3% • Varied by ±1 MeV the s of the Gaussian (12 MeV) that fit the signal Modes with D*0 more affected Sergio Grancagnolo
Branching Ratios Results Phys.Rev.Lett.93:181801,2004 NEW! NEW! NEW! NEW! NEW! NEW! Sergio Grancagnolo Measurements with significance>5
DsJ(2460)+ Angular Analysis (I) _ • Use B0DsJ+D-and B+DsJ+D0 with DsJ+Dsg • B DDsJ+is a transition 0- 0- JP so DsJis polarized • Compute the helicity angleqh of DsJ+Dsgand compare with the predictions for JP=1+ and JP=2+(0+forbidden) Sergio Grancagnolo
DsJ(2460)+ Angular Analysis (II) DsJ events are fitted separately in 5cos(qh) bins not used cut m(Dg)>2.3 Sergio Grancagnolo Simulation is used to correct for detector acceptance
DsJ(2460)+ Angular Analysis (III) • Expected distribution for JP=1+ is: 1-cos2(qh) • Distribution compatible with this case • c2/d.o.f.=3.9/4 • Supporting the Ds1+ hypothesis for this state • Comparison with JP=2+ hypothesis is also provided • c2/d.o.f.=34.5/4 Sergio Grancagnolo
Some Comparisons With Models • Branching ratios smaller than the corresponding BD(*)Ds(*) • Factorization effects could be important and could not cancel in the ratios RD0,1 • support a multiquark hypothesis • Observation of electromagnetic DsJ(2460)+ decay • supports a conventional cs picture • In agreement with prediction from chiral multiplets we measure: Colangelo, De Fazio, Ferrandes: Mod.Phys.Lett. A19:2083,2004 Godfrey Phys.Lett. B568:254,2003 _ Bardeen, Eichten, Hill: Phys.Rev. D68:054024,2003 Sergio Grancagnolo