B0 → D(*)π and B0 → Dρ and determination of sin(2β+γ)



sin(2β+γ) can be extracted from time dependent asymmetries in B decays to D(*) π final states, looking at the interference effect between the decay amplitudes implying b→c and b→u transitions. The time dependent rates can be written as:

R(B0→ D-π+) = Ne-ΓΔt {1+C cos(ΔmdΔt) +S sin(ΔmdΔt)}
R(B0→ D-π+) = Ne-ΓΔt {1-C cos(ΔmdΔt) -S sin(ΔmdΔt)}
R(B0→ D+π-) = Ne-ΓΔt {1+C cos(ΔmdΔt) -S sin(ΔmdΔt)}
R(B0→ D+π-) = Ne-ΓΔt {1-C cos(ΔmdΔt) +S sin(ΔmdΔt)}

where S and C parameters are defined as

S = 2r/(1+r2) sin(2β+γ-δ)
S = 2r/(1+r2) sin(2β+γ+δ)
C = (1-r2)/(1+r2)

and r and δ are the absolute value and the phase of the amplitude ratio A(B0→D-π+) / A(B0→D-π+). It is clear that the ratio r is quite small being of the order of λ |Vub/Vcb| ~0.02. Analogous expressions can be written for D(*) π decays.


Available analyses and measurements

BaBar and Belle provided three different measurements of this channel, with total (both) or partial (BaBar only) reconstruction of the final state.
In hep-ex/0303030 O. Long et al. showed the existence of a correlation between tag side and reconstruction side in time dependent CP measurements at the B-Factories. This is related to the fact that interference between b→c and b→u transitions in B→D X decays can occur also in the tag side. In the case of sin(2β+γ) this is particularly important, since the effect studied on reconstruction side is exactly the same. In order to include it in the fit, S and S entering the time dependent rates are replaced by

a = 2r(sin(2β+γ)cos(δ)
b = 2r'(sin(2β+γ)cos(δ')
c = 2cos(2β+γ)(r sin(δ)-r'sin(δ'))

where r' and δ' are the analogous of r and δ for the tag side. It is important to stress the fact that this interference effect on the tag side can not occur in case the B meson are tagged using semileptonic decays. In other words, r'= 0 when only semileptonic decays are considered. In the following we will consider the fitted quantities a, c(lepton),a(*) and c(*)(lepton), which are functions of r,r(*), δ, δ(*) and 2β+γ.

With the present experimental data, a determination of 2β+γ cannot be obtained from D(*)π(ρ) modes alone. The number of free parameters exceeds the available constraints as shown by equations above. Without further input, one can only find correlations among 2β+γ and the hadronic parameters. The only information one can extract in a model independent way comes from the r vs 2β+γ correlation plot for the various channels.

Values taken from HFAG Summer 2004 (ICHEP 2004, Beijing and FPCP 2004, Daegu)
Observable a(Dπ) c(Dπ) a(D*π) c(D*π) a(Dρ) c(Dρ)
average -0.026 ± 0.021 -0.034 ± 0.030 -0.028 ± 0.012 -0.001 ± 0.017 -0.024 ± 0.033 -0.098 ± 0.058


Distribution of r(Dπ) vs 2β+γ
(EPS), [JPG]

Distribution of r(Dπ) vs 2β+γ

Distribution of r(D*π) vs 2β+γ
(EPS), [JPG]

Distribution of r(D*π) vs 2β+γ

Distribution of r(Drho) vs 2β+γ
(EPS), [JPG]

Distribution of r(Dρ) vs 2β+γ

An independent information on the hadronic parameters would allow a determination of 2β+γ. For instance, assuming SU(3) flavour symmetry and neglecting annihilation contributions, one can estimate BR(B0 → D*-π+) from BR(B0 → Ds*-π+), up to a theoretical uncertainty related to SU(3) breaking effect and to the size of annihilation contributions. We take into account the combination of these two effects assigning a ± 100% flat error to be convolved with the experimental one. In this way, we obtain as input values r(D*π) = 0.015 ± 0.006 ± 0.015 (where the first error is statistical and the second is the flat theoretical error). The same approach for the other two channels give r(Dπ) = 0.020 ± 0.003 ± 0.020 and r(Dρ) = 0.003 ± 0.006 ± 0.003. Under these assumptions, we get a constraint on 2β+γ as shown in the figure below.


sin(2β+γ) distribution from Dπ, D*π and Dρ
(EPS), [JPG]

constraint on 2β+γ combining
Dπ, D*π and Dρ
2β+γ = (± 90 ± 46)o

bound on ρη plane from Dπ, D*π and Dρ
(EPS), [JPG]


constraint on the (ρ ; η) plane
combining Dπ D*π and Dρ.

prior distribution for r(Dπ)
(EPS), [JPG]

prior distribution for r(Dπ)
r(Dπ) = 0.021 ± 0.013 ([0.002,0.041] @ 95% Prob.)

prior distribution for r(D*π)
(EPS), [JPG]

prior distribution for r(D*π)
r(D*π) = 0.017 ± 0.010 ([0.001,0.037] @ 95% Prob.)

prior distribution for r(Dρ)
(EPS), [JPG]

prior distribution for r(Dρ)
r(Dρ) = 0.006 ± 0.004 ([0.000,0.016] @ 95% Prob.)

We would like to thank Riccardo Faccini, Marie Legendre, Cecilia Voena and Marco Zito from BaBar Collaboration for the useful discussions on the interpretation of experimental results.


Valid HTML 4.01!