return function
# ... and the same for the interpolated version (see qcdf_interpolate.py)
def ha_qcdf_interpolate_function(B, V, contribution='all'):
scale = config['renormalization scale']['bvll']
def function(wc_obj, par_dict, q2, cp_conjugate):
return flavio.physics.bdecays.bvll.qcdf_interpolate.helicity_amps_qcdf(q2, par_dict, B, V, cp_conjugate, contribution)
return function
# loop over hadronic transitions and lepton flavours
# BTW, it is not necessary to loop over tau: for tautau final states, the minimum
# q2=4*mtau**2 is so high that QCDF is not valid anymore anyway!
for had in [('B0','K*0'), ('B+','K*+'), ('B0','rho0'), ('B+','rho+'), ('Bs','phi'), ]:
process = had[0] + '->' + had[1] + 'll' # e.g. B0->K*0mumu
quantity = process + ' spectator scattering'
a = AuxiliaryQuantity(name=quantity, arguments=['q2', 'cp_conjugate'])
a.description = ('Contribution to ' + process + ' helicity amplitudes from'
' non-factorizable spectator scattering.')
# Implementation: QCD factorization
iname = process + ' QCDF'
i = Implementation(name=iname, quantity=quantity,
function=ha_qcdf_function(B=had[0], V=had[1]))
i.set_description("QCD factorization")
# Implementation: interpolated QCD factorization
iname = process + ' QCDF interpolated'
i = Implementation(name=iname, quantity=quantity,
function=ha_qcdf_interpolate_function(B=had[0], V=had[1]))
i.set_description("Interpolated version of QCD factorization")
def fct_deltaC9_constant(B, V):
def fct(wc_obj, par_dict, q2, cp_conjugate):
par = par_dict.copy()
return HelicityAmpsDeltaC_9_shift(B, V, par, q2)()
return fct
# AuxiliaryQuantity & Implementation: subleading effects at LOW q^2
for had in [('B0','K*0'), ('B+','K*+'), ('Bs','phi'), ]:
process = had[0] + '->' + had[1] + 'll' # e.g. B0->K*0mumu
quantity = process + ' subleading effects at low q2'
a = AuxiliaryQuantity(name=quantity, arguments=['q2', 'cp_conjugate'])
a.description = ('Contribution to ' + process + ' helicity amplitudes from'
' subleading hadronic effects (i.e. all effects not included'
r'elsewhere) at $q^2$ below the charmonium resonances')
# Implementation: C7-C7'-polynomial
iname = process + ' deltaC7, 7p polynomial'
i = Implementation(name=iname, quantity=quantity,
function=fct_deltaC7C7p_polynomial(B=had[0], V=had[1]))
i.set_description(r"Effective shift in the Wilson coefficient $C_7(\mu_b)$"
r" (in the $0$ and $-$ helicity amplitudes) and"
r" $C_7'(\mu_b)$ (in the $+$ helicity amplitude)"
r" as a first-order polynomial in $q^2$.")
# AuxiliaryQuantity & Implementation: subleading effects at HIGH q^2
return transversity_amps_deltaC7_polynomial(q2, par)
def fct_deltaC9_constant(wc_obj, par_dict, q2, cp_conjugate):
par = par_dict.copy()
if cp_conjugate:
par = conjugate_par(par)
return transversity_amps_deltaC9_constant(q2, par_dict)
# AuxiliaryQuantity & Implementatation: subleading effects at LOW q^2
quantity = 'Lambdab->Lambdall subleading effects at low q2'
a = AuxiliaryQuantity(name=quantity, arguments=['q2', 'cp_conjugate'])
a.description = (
r'Contribution to $\Lambda_b\to \Lambda \ell^+\ell^-$ transversity amplitudes from'
r' subleading hadronic effects (i.e. all effects not included'
r' elsewhere) at $q^2$ below the charmonium resonances')
# Implementation: C7-polynomial
iname = 'Lambdab->Lambdall deltaC7 polynomial'
i = Implementation(name=iname,
quantity=quantity,
function=fct_deltaC7_polynomial)
i.set_description(r"Effective shift in the Wilson coefficient $C_7(\mu_b)$"
r" as a first-order polynomial in $q^2$.")
# AuxiliaryQuantity & Implementatation: subleading effects at HIGH q^2
quantity = 'Lambdab->Lambdall subleading effects at high q2'
a = AuxiliaryQuantity(name=quantity, arguments=['q2', 'cp_conjugate'])
a.description = (
# ... and the same for the interpolated version (see qcdf_interpolate.py)
def ha_qcdf_interpolate_function(B, V, contribution="all"):
scale = config["renormalization scale"]["bvll"]
def function(wc_obj, par_dict, q2, cp_conjugate):
return flavio.physics.bdecays.bvll.qcdf_interpolate.helicity_amps_qcdf(
q2, par_dict, B, V, cp_conjugate, contribution
)
return function
# loop over hadronic transitions and lepton flavours
# BTW, it is not necessary to loop over tau: for tautau final states, the minimum
# q2=4*mtau**2 is so high that QCDF is not valid anymore anyway!
for had in [("B0", "K*0"), ("B+", "K*+"), ("B0", "rho0"), ("B+", "rho+"), ("Bs", "phi")]:
process = had[0] + "->" + had[1] + "ll" # e.g. B0->K*0mumu
quantity = process + " spectator scattering"
a = AuxiliaryQuantity(name=quantity, arguments=["q2", "cp_conjugate"])
a.description = "Contribution to " + process + " helicity amplitudes from" " non-factorizable spectator scattering."
# Implementation: QCD factorization
iname = process + " QCDF"
i = Implementation(name=iname, quantity=quantity, function=ha_qcdf_function(B=had[0], V=had[1]))
i.set_description("QCD factorization")
# Implementation: interpolated QCD factorization
iname = process + " QCDF interpolated"
i = Implementation(name=iname, quantity=quantity, function=ha_qcdf_interpolate_function(B=had[0], V=had[1]))
i.set_description("Interpolated version of QCD factorization")
def fct_deltaC9_constant(wc_obj, par_dict, q2, cp_conjugate):
par = par_dict.copy()
if cp_conjugate:
par = conjugate_par(par)
return transversity_amps_deltaC9_constant(q2, par_dict)
# AuxiliaryQuantity & Implementatation: subleading effects at LOW q^2
quantity = 'Lambdab->Lambdall subleading effects at low q2'
a = AuxiliaryQuantity(name=quantity,
arguments=['q2', 'cp_conjugate'])
a.description = (r'Contribution to $\Lambda_b\to \Lambda \ell^+\ell^-$ transversity amplitudes from'
r' subleading hadronic effects (i.e. all effects not included'
r' elsewhere) at $q^2$ below the charmonium resonances')
# Implementation: C7-polynomial
iname = 'Lambdab->Lambdall deltaC7 polynomial'
i = Implementation(name=iname, quantity=quantity,
function=fct_deltaC7_polynomial)
i.set_description(r"Effective shift in the Wilson coefficient $C_7(\mu_b)$"
r" as a first-order polynomial in $q^2$.")
# AuxiliaryQuantity & Implementatation: subleading effects at HIGH q^2
quantity = 'Lambdab->Lambdall subleading effects at high q2'
a = AuxiliaryQuantity(name=quantity, arguments=['q2', 'cp_conjugate'])
a.description = ('Contribution to $\Lambda_b\to \Lambda \ell^+\ell^-$ transversity amplitudes from'
' subleading hadronic effects (i.e. all effects not included'
