Contributor info: Dr Matteo Fael

Matteo Fael, Ilija S. Milutin, K. Keri Vos

SciPost Phys. Codebases 55 (2025) · published 25 March 2025 | Toggle abstract · pdf

We introduce the code kolya, an open-source tool for phenomenological analyses of inclusive semileptonic $B\to X_c l\bar \nu_l$ decays. It contains a library to compute predictions for the total rate and various kinematic moments within the framework of the heavy quark expansion, utilizing the so-called kinetic scheme. The library currently includes power corrections up to $1/m_b^5$. All available QCD perturbative corrections are implemented via interpolation grids for fast numerical evaluation. We also include effects from new physics parameterized as Wilson coefficients of dimension-six operators in the weak effective theory below the electroweak scale. The library is interfaced to CRunDec for easy evaluation of the quark masses and strong coupling constant at different renormalization scales. The library is developed in Python and does not require compilation. The code is flexible and can be extended to other kinds of inclusive decays such as $B$ meson lifetimes, $B \to X_u l\bar \nu_l$ and $D$ meson decays.

Matteo Fael, Ilija S. Milutin, K. Keri Vos

SciPost Phys. Codebases 55-r1.0 (2025) · published 25 March 2025 | Toggle abstract · src

We introduce the code kolya, an open-source tool for phenomenological analyses of inclusive semileptonic $B\to X_c l\bar \nu_l$ decays. It contains a library to compute predictions for the total rate and various kinematic moments within the framework of the heavy quark expansion, utilizing the so-called kinetic scheme. The library currently includes power corrections up to $1/m_b^5$. All available QCD perturbative corrections are implemented via interpolation grids for fast numerical evaluation. We also include effects from new physics parameterized as Wilson coefficients of dimension-six operators in the weak effective theory below the electroweak scale. The library is interfaced to CRunDec for easy evaluation of the quark masses and strong coupling constant at different renormalization scales. The library is developed in Python and does not require compilation. The code is flexible and can be extended to other kinds of inclusive decays such as $B$ meson lifetimes, $B \to X_u l\bar \nu_l$ and $D$ meson decays.

Matteo Fael, Fabian Lange, Kay Schönwald, Matthias Steinhauser

SciPost Phys. Proc. 7, 041 (2022) · published 22 June 2022 | Toggle abstract · pdf

We discuss a practical approach to compute master integrals entering physical quantities which depend on one parameter. As an example we consider four-loop QCD corrections to the relation between a heavy quark mass defined in the $\overline{\rm MS}$ and on-shell scheme in the presence of a second heavy quark.

Matteo Fael

SciPost Phys. Proc. 7, 026 (2022) · published 21 June 2022 | Toggle abstract · pdf

In these proceedings we discuss the relation between the kinetic and the on-shell schemes for the bottom and the charm quarks and present the methods for the calculation of the mass relation to higher orders in perturbative QCD. The bottom mass in the kinetic scheme is a pivotal input parameter in the inclusive determination of $|V_{cb}|$ from $B\to X_c \ell \nu_\ell$ decays. By combining the relation between the kinetic and the on-shell mass with well-know results for the $\overline{\mathrm{MS}}$-on-shell conversion, we obtain a prediction for $m_b^\mathrm{kin}$ based on precise determinations of $\overline{m}_b(\overline{m}_b)$.

Matteo Fael

SciPost Phys. Proc. 1, 009 (2019) · published 18 February 2019 | Toggle abstract · pdf

These proceedings present the differential decay rates and the branching ratios of the tau and muon decays $\tau \to \ell \ell' \ell' \nu \bar\nu$ (with $\ell,\ell'=\mu,e$) and $\mu \to e e e \nu \bar \nu$ in the Standard Model at NLO. These five-body leptonic decays are a tool to study the Lorentz structure of weak interactions and to test lepton flavour universality. They are also a source of SM background to searches for the lepton-flavour-violating decays $\mu \to e e e$ and $\tau \to \ell \ell' \ell'$. Even if the shift in the branching ratios induced by radiative corrections turns out to be small and of order 1% - mainly due to a running effect of the fine structure constant - locally in the phase space these corrections can reach the 5 - 10% level, depending on the applied cuts. We found for instance that in the phase space region where the neutrino energies are small, and the momenta of the three charged leptons have a similar signature as in $\mu \to eee$ and $\tau \to \ell \ell'\ell'$, the NLO corrections decrease the leading-order prediction by about 10 - 20%.