Module flavio.physics.ddecays.formfactors.bcl

BCL parametrization of $D\to \pi,K$ form factors.

Taken from flavio.physics.bdecays.formfactors.b_p.bcl.py

Functions

def ff(process, q2, par, n=3, t0=None)

Expand source code

def ff(process, q2, par, n=3, t0=None):
    r"""Central value of $D\to P$ form factors in the standard convention
    and BCL parametrization (arXiv:0807.2722).

    The standard convention defines the form factors $f_+$, $f_0$, and $f_T$.
    """
    flavio.citations.register("Bourrely:2008za")
    pd = process_dict[process]
    mpl = par[process + ' BCL m+']
    m0 = par[process + ' BCL m0']
    mB = par['m_'+pd['D']]
    mP = par['m_'+pd['P']]
    ff = {}
    a = {}
    for i in ['f+', 'fT']:
        a[i] = [par[process + ' BCL' + ' a' + str(j) + '_' + i] for j in range(n)]
    # only the first n-1 parameters for f0 are taken from par
    # the nth one is chosen to fulfill the kinematic constraint f+(0)=f0(0)
    a['f0'] = [par[process + ' BCL' + ' a' + str(j) + '_f0'] for j in range(n-1)]
    fplus_q20 = pole('f+', mpl, 0) * param_fplusT(mB, mP, a['f+'], 0, t0)
    f0_q20 = pole('f0', m0, 0) * param_f0(mB, mP, a['f0'], 0, t0)
    an_f0 = (f0_q20-fplus_q20)/z(mB, mP, 0, t0)**(n-1)
    a['f0'].append(an_f0)
    # evaluate FFs
    ff['f+'] = pole('f+', mpl, q2) * param_fplusT(mB, mP, a['f+'], q2, t0)
    ff['fT'] = pole('fT', mpl, q2) * param_fplusT(mB, mP, a['fT'], q2, t0)
    ff['f0'] = pole('f0', m0, q2) * param_f0(mB, mP, a['f0'], q2, t0)
    return ff

Central value of $D\to P$ form factors in the standard convention and BCL parametrization (arXiv:0807.2722).

The standard convention defines the form factors $f_+$, $f_0$, and $f_T$.