|
|
|
---
|
|
|
|
title: Already defined Operators in FlavorKit
|
|
|
|
permalink: /Already_defined_Operators_in_FlavorKit/
|
|
|
|
---
|
|
|
|
|
|
|
|
[Category:FlavorKit](/Category:FlavorKit "wikilink")
|
|
|
|
|
|
|
|
Lagrangian
|
|
|
|
==========
|
|
|
|
|
|
|
|
In this section we present our notation and conventions for the operators (and their corresponding Wilson coefficients) implemented in PreSARAH. Although a more complete list of flavor violating operators can be built, we will concentrate on those implemented in PreSARAH. If necessary, the user can extend it by adding his/her own operators.
|
|
|
|
|
|
|
|
The interaction Lagrangian relevant for flavor violating processes can be written as
|
|
|
|
|
|
|
|
$\\label{eq:L-total}
|
|
|
|
{\\cal L}_{\\text{FV}} = {\\cal L}_{\\text{LFV}} + {\\cal L}_{\\text{QFV}} \\, .$
|
|
|
|
|
|
|
|
The first piece contains the operators that can trigger lepton flavor violation whereas the second piece contains the operators responsible for quark flavor violation.
|
|
|
|
|
|
|
|
The general Lagrangian relevant for lepton flavor violation can be written as
|
|
|
|
|
|
|
|
$\\label{eq:L-LFV}
|
|
|
|
{\\cal L}_{\\text{LFV}} = {\\cal L}_{\\ell \\ell \\gamma} + {\\cal L}_{\\ell \\ell Z}
|
|
|
|
+ {\\cal L}_{\\ell \\ell h} + {\\cal L}_{4 \\ell} + {\\cal L}_{2 \\ell 2q} \\, .$
|
|
|
|
|
|
|
|
The first term contains theℓ − ℓ − *γ* interaction, given by
|
|
|
|
|
|
|
|
$\\label{eq:L-llg}
|
|
|
|
{\\cal L}_{\\ell \\ell \\gamma} = e \\, \\bar \\ell_\\beta \\left\[ \\gamma^\\mu \\left(K_1^L P_L + K_1^R P_R \\right) + i m_{\\ell_\\alpha} \\sigma^{\\mu \\nu} q_\\nu \\left(K_2^L P_L + K_2^R P_R \\right) \\right\] \\ell_\\alpha A_\\mu + h.c.$
|
|
|
|
|
|
|
|
Here*e* is the electric charge,*q* the photon momentum,$P_{L,R} =
|
|
|
|
\\frac{1}{2} (1 \\mp \\gamma_5)$ are the usual chirality projectors andℓ<sub>*α*, *β*</sub> denote the lepton flavors. For practical reasons, we will always consider the photonic contributions independently, and we will not include them in other vector operators. On the contrary, the*Z*- and Higgs boson contributions will be included whenever possible. Therefore, theℓ − ℓ − *Z* andℓ − ℓ − *h* interaction Lagrangians will only be used for observables involving real*Z*- and Higgs bosons. These two Lagrangians can be written as
|
|
|
|
|
|
|
|
$\\label{eq:L-llZ}
|
|
|
|
{\\cal L}_{\\ell \\ell Z} = \\bar \\ell_\\beta \\left\[ \\gamma^\\mu \\left(R_1^L P_L + R_1^R P_R \\right) + p^\\mu \\left(R_2^L P_L + R_2^R P_R \\right) \\right\] \\ell_\\alpha Z_\\mu \\, ,$
|
|
|
|
|
|
|
|
where*p* is theℓ<sub>*β*</sub> 4-momentum, and
|
|
|
|
|
|
|
|
$\\label{eq:L-llh}
|
|
|
|
{\\cal L}_{\\ell \\ell h} = \\bar \\ell_\\beta \\left(S_L P_L + S_R P_R \\right) \\ell_\\alpha h \\, .$
|
|
|
|
|
|
|
|
The general4ℓ 4-fermion interaction Lagrangian can be written as
|
|
|
|
|
|
|
|
${\\cal L}_{4 \\ell} = \\sum_{\\substack{I=S,V,T\\\\X,Y=L,R}} A_{XY}^I \\bar \\ell_\\beta \\Gamma_I P_X \\ell_\\alpha \\bar \\ell_\\delta \\Gamma_I P_Y \\ell_\\gamma + h.c. \\, , \\label{eq:L-4L}$
|
|
|
|
|
|
|
|
whereℓ<sub>*α*, *β*, *γ*, *δ*</sub> denote the lepton flavors and*Γ*<sub>*S*</sub> = 1,*Γ*<sub>*V*</sub> = *γ*<sub>*μ*</sub> and*Γ*<sub>*T*</sub> = *σ*<sub>*μ**ν*</sub>. We omit flavor indices in the Wilson coefficients for the sake of clarity. This Lagrangian contains the most general form compatible with Lorentz invariance. The Wilson coefficients*A*<sub>*L**R*</sub><sup>*S*</sup> and*A*<sub>*R**L*</sub><sup>*S*</sup> were included in , but absent in . As previously stated, the coefficients in Eq. do not include photonic contributions, but they include Z-boson and scalar ones. Finally, the general2ℓ2*q* four fermion interaction Lagrangian at the quark level is given by
|
|
|
|
|
|
|
|
${\\cal L}_{2 \\ell 2q} = {\\cal L}_{2 \\ell 2d} + {\\cal L}_{2 \\ell 2u} \\label{eq:L-2L2Q}$
|
|
|
|
|
|
|
|
where
|
|
|
|
|
|
|
|
$\\begin{aligned}
|
|
|
|
{\\cal L}_{2 \\ell 2d} =
|
|
|
|
&& \\sum_{\\substack{I=S,V,T\\\\X,Y=L,R}} B_{XY}^I \\bar \\ell_\\beta \\Gamma_I P_X \\ell_\\alpha \\bar d_\\gamma \\Gamma_I P_Y d_\\gamma + h.c. \\label{eq:L-2L2D} \\\\
|
|
|
|
{\\cal L}_{2 \\ell 2u} = && \\left. {\\cal L}_{2 \\ell 2d} \\right|_{d \\to u, \\, B \\to C} \\label{eq:L-2L2U} \\, .\\end{aligned}$
|
|
|
|
|
|
|
|
Here*d*<sub>*γ*</sub> denotes the d-quark flavor.
|
|
|
|
|
|
|
|
Let us now consider the Lagrangian relevant for quark flavor violation. This can be written as
|
|
|
|
|
|
|
|
$\\label{eq:L-QFV}
|
|
|
|
{\\cal L}_{\\text{QFV}} = {\\cal L}_{q q \\gamma} + {\\cal L}_{q q g} + {\\cal L}_{4 d} +
|
|
|
|
{\\cal L}_{2d2l} + {\\cal L}_{2d2\\nu} + {\\cal L}_{du\\ell\\nu} + {\\cal L}_{d d H} \\, .$
|
|
|
|
|
|
|
|
The first two terms correspond to operators that couple quark bilinears to massless gauge bosons. These are
|
|
|
|
|
|
|
|
$\\begin{aligned}
|
|
|
|
{\\cal L}_{q q \\gamma} = && e \\left\[ \\bar d_\\beta \\sigma_{\\mu \\nu} \\left( m_{d_\\beta} Q_1^L P_L + m_{d_\\alpha} Q_1^R P_R \\right) d_\\alpha \\right\] F^{\\mu \\nu} \\label{eq:L-qqgamma} \\\\
|
|
|
|
{\\cal L}_{q q g} = && g_s \\left\[ \\bar d_\\beta \\sigma_{\\mu \\nu} \\left( m_{d_\\beta} Q_2^L P_L + m_{d_\\alpha} Q_2^R P_R \\right) T^a d_\\alpha \\right\] G_a^{\\mu \\nu} \\label{eq:L-qqgluon} \\, .\\end{aligned}$
|
|
|
|
|
|
|
|
Here*T*<sup>*a*</sup> are*S**U*(3) matrices. The Wilson coefficients*Q*<sub>1, 2</sub><sup>*L*, *R*</sup> can be easily related to the usual*C*<sub>7, 8</sub><sup>(′)</sup> coefficients, sometimes normalized with an additional$\\frac{1}{16 \\pi^2}$ factor. The4*d* four fermion interaction Lagrangian can be written as
|
|
|
|
|
|
|
|
${\\cal L}_{4 d} = \\sum_{\\substack{I=S,V,T\\\\X,Y=L,R}} D_{XY}^I \\bar d_\\beta \\Gamma_I P_X d_\\alpha \\bar d_\\delta \\Gamma_I P_Y d_\\gamma + h.c. \\, , \\label{eq:L-4D}$
|
|
|
|
|
|
|
|
where*d*<sub>*α*, *β*, *γ*, *δ*</sub> denote the lepton flavors. Again, we omit flavor indices in the Wilson coefficients for the sake of clarity. The2*d*2ℓ four fermion interaction Lagrangian is given by
|
|
|
|
|
|
|
|
${\\cal L}_{2d 2 \\ell} = \\sum_{\\substack{I=S,V,T\\\\X,Y=L,R}} E_{XY}^I \\bar d_\\beta \\Gamma_I P_X d_\\alpha \\bar \\ell_\\gamma \\, \\Gamma_I P_Y \\ell_\\gamma + hc \\label{eq:L-2D2L} \\, .$
|
|
|
|
|
|
|
|
Hereℓ<sub>*γ*</sub> denotes the lepton flavor.${\\cal L}_{2d 2 \\ell}$ should not be confused with${\\cal L}_{2 \\ell 2d}$. In the former case one has QFV operators, whereas in the latter one has LFV operators. This distinction has been made for practical reasons. The2*d*2*ν* and*d**u*ℓ*ν* terms of the QFV Lagrangian are
|
|
|
|
|
|
|
|
$\\begin{aligned}
|
|
|
|
{\\cal L}_{2d 2 \\nu} &=& \\sum_{X,Y=L,R} F_{XY}^V \\bar d_\\beta \\gamma_\\mu P_X d_\\alpha \\bar \\nu_\\gamma \\gamma^\\mu P_Y \\nu_\\gamma + hc \\label{eq:L-2D2V} \\\\
|
|
|
|
{\\cal L}_{du\\ell\\nu} &=& \\sum_{\\substack{I=S,V\\\\X,Y=L,R}} G_{XY}^I \\bar d_\\beta \\Gamma_I P_X u_\\alpha \\bar \\ell_\\gamma \\, \\Gamma_I P_Y \\nu_\\gamma + hc \\label{eq:L-DULV} \\, .\\end{aligned}$
|
|
|
|
|
|
|
|
Note that we have not introduced scalar or tensor2*d*2*ν* operators, nor tensor*d**u*ℓ*ν* ones, and that lepton flavor (denoted by the index*γ*) is conserved in these operators. Finally, we have also included a term in the Lagrangian accounting for operators of the type(*d̄**Γ**d*)*S* and(*d̄**Γ**d*)*P*, where*S* (</math>P</math>) is a virtual [1] scalar (pseudoscalar) state. This piece can be written as
|
|
|
|
|
|
|
|
$\\label{eq:L-ddSP}
|
|
|
|
{\\cal L}_{d d H} = \\bar d_\\beta \\left(H_L^S P_L + H_R^S P_R \\right) d_\\alpha S
|
|
|
|
+ \\bar d_\\beta \\left(H_L^P P_L + H_R^P P_R \\right) d_\\alpha P \\, .$
|
|
|
|
|
|
|
|
Operators available by default in the SPheno output of SARAH
|
|
|
|
============================================================
|
|
|
|
|
|
|
|
The operators presented abvoe have been implemented by using the results of PreSARAH in SARAH. Those are exported to SPheno. In the following a list of all internal names for these operators, which can be used in the calculation of new flavor observables is given.
|
|
|
|
|
|
|
|
2-Fermion-1-Boson operators
|
|
|
|
---------------------------
|
|
|
|
|
|
|
|
These operators are arrays with either two or three elements. While operators involving vector bosons have always dimension3 × 3, those with scalars have dimension3 × 3 × *n*<sub>*g*</sub>.*n*<sub>*g*</sub> is the number of generations of the considered scalar and for*n*<sub>*g*</sub> = 1 the last index is dropped.
|
|
|
|
|
|
|
|
(*d̄*<sub>*β*</sub>*σ*<sub>*μ**ν*</sub>*Γ**d*<sub>*α*</sub>)*F*<sup>*μ**ν*</sup> and(*d̄*<sub>*β*</sub>*σ*<sub>*μ**ν*</sub>*Γ**d*<sub>*α*</sub>)*G*<sup>*μ**ν*</sup>
|
|
|
|
|
|
|
|
| Variable | Operator | Name |
|
|
|
|
|----------|-----------------------------------------------------------------------------------------------------------------------------------------------|-------------------------------|
|
|
|
|
| CC7 | *e**m*<sub>*d*<sub>*β*</sub></sub>(*d̄*<sub>*β*</sub>*σ*<sub>*μ**ν*</sub>*P*<sub>*L*</sub>*d*<sub>*α*</sub>)*F*<sup>*μ**ν*</sup> | *Q*<sub>1</sub><sup>*L*</sup> |
|
|
|
|
| CC7p | *e**m*<sub>*d*<sub>*α*</sub></sub>(*d̄*<sub>*β*</sub>*σ*<sub>*μ**ν*</sub>*P*<sub>*R*</sub>*d*<sub>*α*</sub>)*F*<sup>*μ**ν*</sup> | *Q*<sub>1</sub><sup>*R*</sup> |
|
|
|
|
| CC8 | *g*<sub>*s*</sub>*m*<sub>*d*<sub>*β*</sub></sub>(*d̄*<sub>*β*</sub>*σ*<sub>*μ**ν*</sub>*P*<sub>*L*</sub>*d*<sub>*α*</sub>)*G*<sup>*μ**ν*</sup> | *Q*<sub>2</sub><sup>*L*</sup> |
|
|
|
|
| CC8p | *g*<sub>*s*</sub>*m*<sub>*d*<sub>*α*</sub></sub>(*d̄*<sub>*β*</sub>*σ*<sub>*μ**ν*</sub>*P*<sub>*R*</sub>*d*<sub>*α*</sub>)*G*<sup>*μ**ν*</sup> | *Q*<sub>2</sub><sup>*R*</sup> |
|
|
|
|
|
|
|
|
These operators are derived by PreSARAH with the following input files
|
|
|
|
|
|
|
|
NameProcess="Gamma2Q";
|
|
|
|
|
|
|
|
ConsideredProcess = "2Fermion1Vector";
|
|
|
|
FermionOrderExternal={1,2};
|
|
|
|
NeglectMasses={3};
|
|
|
|
|
|
|
|
|
|
|
|
ExternalFields= {bar[BottomQuark], BottomQuark,Photon};
|
|
|
|
CombinationGenerations = {{3,2}};
|
|
|
|
|
|
|
|
|
|
|
|
AllOperators={
|
|
|
|
{OA2qSL,Op[7] Pair[ec[3],k[1]]},
|
|
|
|
{OA2qSR,Op[6] Pair[ec[3],k[1]]},
|
|
|
|
{OA2qVL,Op[7,ec[3]]},
|
|
|
|
{OA2qVR,Op[6,ec[3]]}
|
|
|
|
};
|
|
|
|
|
|
|
|
OutputFile = "Gamma2Q.m";
|
|
|
|
|
|
|
|
Filters = {};
|
|
|
|
|
|
|
|
NameProcess="Gluon2Q";
|
|
|
|
|
|
|
|
ConsideredProcess = "2Fermion1Vector";
|
|
|
|
FermionOrderExternal={1,2};
|
|
|
|
NeglectMasses={3};
|
|
|
|
|
|
|
|
|
|
|
|
ExternalFields= {bar[BottomQuark], BottomQuark,Gluon};
|
|
|
|
CombinationGenerations = {{3,2}};
|
|
|
|
|
|
|
|
|
|
|
|
AllOperators={
|
|
|
|
{OG2qSL,Op[7] Pair[ec[3],k[1]]},
|
|
|
|
{OG2qSR,Op[6] Pair[ec[3],k[1]]}
|
|
|
|
};
|
|
|
|
|
|
|
|
OutputFile = "Gluon2Q.m";
|
|
|
|
|
|
|
|
Filters = {};
|
|
|
|
|
|
|
|
The normalization is changed to match the standard definitions by
|
|
|
|
|
|
|
|
ProcessWrapper = True;
|
|
|
|
NameProcess = "Gamma2Q"
|
|
|
|
ExternalFields = {bar[BottomQuark], BottomQuark, Photon};
|
|
|
|
|
|
|
|
SumContributionsOperators["Gamma2Q"] = {
|
|
|
|
{CC7, OA2qSL},
|
|
|
|
{CC7p, OA2qSR}
|
|
|
|
};
|
|
|
|
|
|
|
|
NormalizationOperators["Gamma2Q"] ={
|
|
|
|
"CC7(2,:) = 0.25_dp*CC7(2,:)/sqrt(Alpha_160*4*Pi)/MFd(2)",
|
|
|
|
"CC7(3,:) = 0.25_dp*CC7(3,:)/sqrt(Alpha_160*4*Pi)/MFd(3)",
|
|
|
|
"CC7p(2,:) = 0.25_dp*CC7p(2,:)/sqrt(Alpha_160*4*Pi)/MFd(2)",
|
|
|
|
"CC7p(3,:) = 0.25_dp*CC7p(3,:)/sqrt(Alpha_160*4*Pi)/MFd(3)",
|
|
|
|
|
|
|
|
"CC7SM(2,:) = 0.25_dp*CC7SM(2,:)/sqrt(Alpha_160*4*Pi)/MFd(2)",
|
|
|
|
"CC7SM(3,:) = 0.25_dp*CC7SM(3,:)/sqrt(Alpha_160*4*Pi)/MFd(3)",
|
|
|
|
"CC7pSM(2,:) = 0.25_dp*CC7pSM(2,:)/sqrt(Alpha_160*4*Pi)/MFd(2)",
|
|
|
|
"CC7pSM(3,:) = 0.25_dp*CC7pSM(3,:)/sqrt(Alpha_160*4*Pi)/MFd(3)"
|
|
|
|
};
|
|
|
|
|
|
|
|
ProcessWrapper = True;
|
|
|
|
NameProcess = "Gluon2Q"
|
|
|
|
ExternalFields = {bar[BottomQuark], BottomQuark, Gluon};
|
|
|
|
|
|
|
|
SumContributionsOperators["Gluon2Q"] = {
|
|
|
|
{CC8, OG2qSL},
|
|
|
|
{CC8p, OG2qSR}};
|
|
|
|
|
|
|
|
NormalizationOperators["Gluon2Q"] ={
|
|
|
|
"CC8(2,:) = 0.25_dp*CC8(2,:)/sqrt(AlphaS_160*4*Pi)/MFd(2)",
|
|
|
|
"CC8(3,:) = 0.25_dp*CC8(3,:)/sqrt(AlphaS_160*4*Pi)/MFd(3)",
|
|
|
|
"CC8p(2,:) = 0.25_dp*CC8p(2,:)/sqrt(AlphaS_160*4*Pi)/MFd(2)",
|
|
|
|
"CC8p(3,:) = 0.25_dp*CC8p(3,:)/sqrt(AlphaS_160*4*Pi)/MFd(3)",
|
|
|
|
|
|
|
|
"CC8SM(2,:) = 0.25_dp*CC8SM(2,:)/sqrt(AlphaS_160*4*Pi)/MFd(2)",
|
|
|
|
"CC8SM(3,:) = 0.25_dp*CC8SM(3,:)/sqrt(AlphaS_160*4*Pi)/MFd(3)",
|
|
|
|
"CC8pSM(2,:) = 0.25_dp*CC8pSM(2,:)/sqrt(AlphaS_160*4*Pi)/MFd(2)",
|
|
|
|
"CC8pSM(3,:) = 0.25_dp*CC8pSM(3,:)/sqrt(AlphaS_160*4*Pi)/MFd(3)"
|
|
|
|
|
|
|
|
};
|
|
|
|
|
|
|
|
$\\bar \\ell_\\beta \\left( q^2 \\gamma^\\mu + i m_{\\ell_\\alpha} \\sigma^{\\mu \\nu} q_\\nu \\right) \\ell_\\alpha A_\\mu$
|
|
|
|
|
|
|
|
| Variable | Operator | Name |
|
|
|
|
|----------|-------------------------------------------------------------------------------------------------------|-------------------------------|
|
|
|
|
| K2L | $e m_{\\ell_\\alpha} (\\bar \\ell_\\beta \\sigma_{\\mu\\nu} P_L \\ell_\\alpha) q^{\\nu} A^\\mu$ | *K*<sub>2</sub><sup>*L*</sup> |
|
|
|
|
| K1L | $q^2 (\\bar \\ell_\\beta \\gamma_\\mu P_L \\ell_\\alpha) A^\\mu$ | *K*<sub>1</sub><sup>*L*</sup> |
|
|
|
|
| K2R | $e m_{\\ell_\\alpha} (\\bar \\ell_\\beta \\sigma_{\\mu\\nu} P_R \\ell_\\alpha) q^{\\nu} A^\\mu$ | *K*<sub>2</sub><sup>*L*</sup> |
|
|
|
|
| K1R | $q^2 (\\bar \\ell_\\beta \\gamma_\\nu P_R \\ell_\\alpha) A^\\mu$ | *K*<sub>1</sub><sup>*R*</sup> |
|
|
|
|
|
|
|
|
These operators are derived by PreSARAH with the following input files
|
|
|
|
|
|
|
|
NameProcess="Gamma2l";
|
|
|
|
|
|
|
|
ConsideredProcess = "2Fermion1Vector";
|
|
|
|
FermionOrderExternal={1,2};
|
|
|
|
NeglectMasses={3};
|
|
|
|
|
|
|
|
|
|
|
|
ExternalFields= {bar[ChargedLepton], ChargedLepton,Photon};
|
|
|
|
CombinationGenerations = {{2,1},{3,1},{3,2}};
|
|
|
|
|
|
|
|
|
|
|
|
AllOperators={
|
|
|
|
{OA2lSL,Op[6] Pair[ec[3],k[1]]},
|
|
|
|
{OA2lSR,Op[7] Pair[ec[3],k[1]]},
|
|
|
|
{OA1L,Op[6,ec[3]] Pair[k[3],k[3]]},
|
|
|
|
{OA1R,Op[7,ec[3]] Pair[k[3],k[3]]}
|
|
|
|
};
|
|
|
|
|
|
|
|
OutputFile = "Gamma2l.m";
|
|
|
|
|
|
|
|
Filters = {};
|
|
|
|
|
|
|
|
The normalization is changed to match the standard definitions by
|
|
|
|
|
|
|
|
ProcessWrapper = True;
|
|
|
|
NameProcess = "Gamma2l"
|
|
|
|
ExternalFields = {bar[ChargedLepton], ChargedLepton, Photon};
|
|
|
|
|
|
|
|
SumContributionsOperators["Gamma2l"] = {
|
|
|
|
{K1L, OA1L},
|
|
|
|
{K1R, OA1R},
|
|
|
|
{K2L, OA2lSL},
|
|
|
|
{K2R, OA2lSR}};
|
|
|
|
|
|
|
|
NormalizationOperators["Gamma2l"] ={
|
|
|
|
"K1L = K1L/sqrt(Alpha_MZ*4*Pi)",
|
|
|
|
"K1R = K1R/sqrt(Alpha_MZ*4*Pi)",
|
|
|
|
"K2L(2,:) = -0.5_dp*K2L(2,:)/sqrt(Alpha_MZ*4*Pi)/MFe(2)",
|
|
|
|
"K2L(3,:) = -0.5_dp*K2L(3,:)/sqrt(Alpha_MZ*4*Pi)/MFe(3)",
|
|
|
|
"K2R(2,:) = -0.5_dp*K2R(2,:)/sqrt(Alpha_MZ*4*Pi)/MFe(2)",
|
|
|
|
"K2R(3,:) = -0.5_dp*K2R(3,:)/sqrt(Alpha_MZ*4*Pi)/MFe(3)"
|
|
|
|
};
|
|
|
|
|
|
|
|
$(\\bar \\ell \\Gamma \\ell) Z$
|
|
|
|
|
|
|
|
| Variable | Operator | Name |
|
|
|
|
|----------|------------------------------------------------------|-------------------------------|
|
|
|
|
| OZ2lVL | $(\\bar{\\ell} \\, \\gamma^\\mu P_L \\ell) Z_\\mu$ | *R*<sub>1</sub><sup>*L*</sup> |
|
|
|
|
| OZ2lSL | $(\\bar{\\ell} p^\\mu P_L \\ell) Z_\\mu$ | *R*<sub>2</sub><sup>*L*</sup> |
|
|
|
|
| OZ2lVR | $(\\bar{\\ell} \\, \\gamma^\\mu P_R \\ell) Z_\\mu$ | *R*<sub>1</sub><sup>*R*</sup> |
|
|
|
|
| OZ2lSR | $(\\bar{\\ell} p^\\mu P_R \\ell) Z_\\mu$ | *R*<sub>2</sub><sup>*R*</sup> |
|
|
|
|
|
|
|
|
In the following we omit flavor indices for the sake of simplicity. These operators are derived by PreSARAH with the following input files
|
|
|
|
|
|
|
|
NameProcess="Z2l";
|
|
|
|
|
|
|
|
ConsideredProcess = "2Fermion1Vector";
|
|
|
|
FermionOrderExternal={1,2};
|
|
|
|
NeglectMasses={1,2};
|
|
|
|
|
|
|
|
|
|
|
|
ExternalFields= {ChargedLepton,bar[ChargedLepton],Zboson};
|
|
|
|
CombinationGenerations = {{1,2},{1,3},{2,3}};
|
|
|
|
|
|
|
|
|
|
|
|
AllOperators={
|
|
|
|
{OZ2lSL,Op[7] Pair[ec[3],k[1]]}, {OZ2lSR,Op[6] Pair[ec[3],k[1]]},
|
|
|
|
{OZ2lVL,Op[7,ec[3]]}, {OZ2lVR,Op[6,ec[3]]}
|
|
|
|
};
|
|
|
|
|
|
|
|
OutputFile = "Z2l.m";
|
|
|
|
|
|
|
|
Filters = {};
|
|
|
|
|
|
|
|
$(\\bar{\\ell} \\Gamma \\ell) h$
|
|
|
|
|
|
|
|
| Variable | Operator | Name |
|
|
|
|
|----------|---------------------------------|-------------------|
|
|
|
|
| OH2lSL | $\\bar{\\ell} P_L \\ell \\, h$ | *S*<sub>*L*</sub> |
|
|
|
|
| OH2lSR | $\\bar{\\ell} P_R \\ell \\, h$ | *S*<sub>*R*</sub> |
|
|
|
|
|
|
|
|
These operators are derived by PreSARAH with the following input files
|
|
|
|
|
|
|
|
NameProcess="H2l";
|
|
|
|
|
|
|
|
ConsideredProcess = "2Fermion1Scalar";
|
|
|
|
FermionOrderExternal={1,2};
|
|
|
|
NeglectMasses={1,2};
|
|
|
|
|
|
|
|
|
|
|
|
ExternalFields= {ChargedLepton,bar[ChargedLepton],HiggsBoson};
|
|
|
|
CombinationGenerations = {{1,2,ALL},{1,3,ALL},{2,3,ALL}};
|
|
|
|
|
|
|
|
|
|
|
|
AllOperators={{OH2lSL,Op[7]},
|
|
|
|
{OH2lSR,Op[6]}
|
|
|
|
};
|
|
|
|
|
|
|
|
OutputFile = "H2l.m";
|
|
|
|
|
|
|
|
Filters = {};
|
|
|
|
|
|
|
|
(*d̄**Γ**d*)*S* and(*d̄**Γ**d*)*P*
|
|
|
|
|
|
|
|
| Variable | Operator | Name |
|
|
|
|
|----------|-----------------------------|---------------------------------|
|
|
|
|
| OH2qSL | *d̄**P*<sub>*L*</sub>*d* *S* | </math>H_L^S</math> |
|
|
|
|
| OH2qSR | *d̄**P*<sub>*R*</sub>*d* *S* | *H*<sub>*R*</sub><sup>*S*</sup> |
|
|
|
|
| OAh2qSL | *d̄**P*<sub>*L*</sub>*d* *P* | </math>H_L^P</math> |
|
|
|
|
| OAh2qSR | *d̄**P*<sub>*R*</sub>*d* *P* | *H*<sub>*R*</sub><sup>*P*</sup> |
|
|
|
|
|
|
|
|
These auxiliary [2] operators are derived by PreSARAH with the following input files
|
|
|
|
|
|
|
|
NameProcess="H2q";
|
|
|
|
|
|
|
|
(* operators needed for double penguins with internal scalars *)
|
|
|
|
(* we neglect therefore the mass of the scalar in the loop functions *)
|
|
|
|
(* and treat it as massless *)
|
|
|
|
|
|
|
|
ConsideredProcess = "2Fermion1Scalar";
|
|
|
|
FermionOrderExternal={2,1};
|
|
|
|
NeglectMasses={3};
|
|
|
|
|
|
|
|
|
|
|
|
ExternalFields= {DownQuark,bar[DownQuark],HiggsBoson};
|
|
|
|
CombinationGenerations = {{2,1,ALL},{3,1,ALL},{3,2,ALL}};
|
|
|
|
|
|
|
|
|
|
|
|
AllOperators={{OH2qSL,Op[7]},
|
|
|
|
{OH2qSR,Op[6]}
|
|
|
|
};
|
|
|
|
|
|
|
|
OutputFile = "H2q.m";
|
|
|
|
|
|
|
|
Filters = {};
|
|
|
|
|
|
|
|
NameProcess="A2q";
|
|
|
|
|
|
|
|
(* operators needed for double penguins with internal scalars *)
|
|
|
|
(* we neglect therefore the mass of the scalar in the loop functions *)
|
|
|
|
(* and treat it as massless *)
|
|
|
|
|
|
|
|
ConsideredProcess = "2Fermion1Scalar";
|
|
|
|
FermionOrderExternal={2,1};
|
|
|
|
NeglectMasses={3};
|
|
|
|
|
|
|
|
|
|
|
|
ExternalFields= {DownQuark,bar[DownQuark],PseudoScalar};
|
|
|
|
CombinationGenerations = {{2,1,ALL},{3,1,ALL},{3,2,ALL}};
|
|
|
|
|
|
|
|
|
|
|
|
AllOperators={{OAh2qSL,Op[7]},
|
|
|
|
{OAh2qSR,Op[6]}
|
|
|
|
};
|
|
|
|
|
|
|
|
OutputFile = "A2q.m";
|
|
|
|
|
|
|
|
Filters = {};
|
|
|
|
|
|
|
|
4-Fermion operators
|
|
|
|
-------------------
|
|
|
|
|
|
|
|
All operators listed below carry four indices and have dimension3 × 3 × 3 × 3. In addition, the user can access the different contributions of all operators from tree-level diagrams, as well as penguin and box diagrams. The name conventions are as follows: for each operator `op` the additional parameter exist
|
|
|
|
|
|
|
|
- `TSop` : tree-level contributions with scalar propagator
|
|
|
|
- `TVop` : tree-level contributions with scalar propagator
|
|
|
|
- `PSop` : sum of penguin and self-energy contributions with scalar propagator
|
|
|
|
- `PVop` : sum of penguin and self-energy contributions with scalar propagator
|
|
|
|
- `Bop` : box contributions.
|
|
|
|
|
|
|
|
We will denote the 4-fermion operators involving two leptons and two down-type quarks depending on whether they lead to LFV or to QFV processes
|
|
|
|
ℓℓ*d**d*
|
|
|
|
for LFV and*d**d*ℓℓ for QFV.
|
|
|
|
|
|
|
|
$(\\bar{d} \\Gamma d) (\\bar{\\ell} \\Gamma^\\prime \\ell)$ and(*d̄**Γ**d*)(*ν̄**Γ*<sup>′</sup>*ν*)
|
|
|
|
|
|
|
|
| Variable | Operator | Name |
|
|
|
|
|----------|--------------------------------------------------------------------------------------|------------------------------------|
|
|
|
|
| OddllSLL | $(\\bar{d} P_L d) (\\bar{\\ell} P_L \\ell)$ | *E*<sub>*L**L*</sub><sup>*S*</sup> |
|
|
|
|
| OddllSRR | $(\\bar{d} P_R d) (\\bar{\\ell} P_R \\ell)$ | *E*<sub>*R**R*</sub><sup>*S*</sup> |
|
|
|
|
| OddllSLR | $(\\bar{d} P_L d) (\\bar{\\ell} P_R \\ell)$ | *E*<sub>*L**R*</sub><sup>*S*</sup> |
|
|
|
|
| OddllSRL | $(\\bar{d} P_R d) (\\bar{\\ell} P_L \\ell)$ | *E*<sub>*R**L*</sub><sup>*S*</sup> |
|
|
|
|
| OddllVLL | $(\\bar{d} \\gamma_\\mu P_L d) (\\bar{\\ell} \\gamma^\\mu P_L \\ell)$ | *E*<sub>*L**L*</sub><sup>*V*</sup> |
|
|
|
|
| OddvvVLL | (*d̄**γ*<sub>*μ*</sub>*P*<sub>*L*</sub>*d*)(*ν̄**γ*<sup>*μ*</sup>*P*<sub>*R*</sub>*ν*) | *F*<sub>*L**L*</sub><sup>*V*</sup> |
|
|
|
|
| OddllVRR | $(\\bar{d} \\gamma_\\mu P_R d) (\\bar{\\ell} \\gamma^\\mu P_R \\ell)$ | *E*<sub>*R**R*</sub><sup>*V*</sup> |
|
|
|
|
| OddvvVRR | (*d̄**γ*<sub>*μ*</sub>*P*<sub>*R*</sub>*d*)(*ν̄**γ*<sup>*μ*</sup>*P*<sub>*R*</sub>*ν*) | *F*<sub>*R**R*</sub><sup>*V*</sup> |
|
|
|
|
| OddllVLR | $(\\bar{d} \\gamma_\\mu P_L d) (\\bar{\\ell} \\gamma^\\mu P_R \\ell)$ | *E*<sub>*L**R*</sub><sup>*V*</sup> |
|
|
|
|
| OddvvVLR | (*d̄**γ*<sub>*μ*</sub>*P*<sub>*L*</sub>*d*)(*ν̄**γ*<sup>*μ*</sup>*P*<sub>*R*</sub>*ν*) | *F*<sub>*L**R*</sub><sup>*V*</sup> |
|
|
|
|
| OddllVRL | $(\\bar{d} \\gamma_\\mu P_R d) (\\bar{\\ell} \\gamma^\\mu P_L \\ell)$ | *E*<sub>*R**L*</sub><sup>*V*</sup> |
|
|
|
|
| OddvvVRL | (*d̄**γ*<sub>*μ*</sub>*P*<sub>*R*</sub>*d*)(*ν̄**γ*<sup>*μ*</sup>*P*<sub>*L*</sub>*ν*) | *F*<sub>*R**L*</sub><sup>*V*</sup> |
|
|
|
|
| OddllTLL | $(\\bar{d} \\sigma_{\\mu\\nu} P_L d) (\\bar{\\ell} \\sigma^{\\mu\\nu} P_L \\ell)$ | *E*<sub>*L**L*</sub><sup>*T*</sup> |
|
|
|
|
| OddllTRR | $(\\bar{d} \\sigma_{\\mu\\nu} P_R d) (\\bar{\\ell} \\sigma^{\\mu\\nu} P_R \\ell)$ | *E*<sub>*R**R*</sub><sup>*T*</sup> |
|
|
|
|
| OddllTLR | $(\\bar{d} \\sigma_{\\mu\\nu} P_L d) (\\bar{\\ell} \\sigma^{\\mu\\nu} P_R \\ell)$ | *E*<sub>*L**R*</sub><sup>*T*</sup> |
|
|
|
|
| OddllTRL | $(\\bar{d} \\sigma_{\\mu\\nu} P_R d) (\\bar{\\ell} \\sigma^{\\mu\\nu} P_L \\ell)$ | *E*<sub>*R**L*</sub><sup>*T*</sup> |
|
|
|
|
|
|
|
|
These operators are derived by PreSARAH with the following input files
|
|
|
|
|
|
|
|
NameProcess="2d2L";
|
|
|
|
|
|
|
|
ConsideredProcess = "4Fermion";
|
|
|
|
FermionOrderExternal={2,1,4,3};
|
|
|
|
NeglectMasses={1,2,3,4};
|
|
|
|
|
|
|
|
|
|
|
|
ExternalFields= {DownQuark,bar[DownQuark],ChargedLepton,bar[ChargedLepton]};
|
|
|
|
|
|
|
|
CombinationGenerations = {{3,1,1,1}, {3,1,2,2}, {3,1,3,3},
|
|
|
|
{3,2,1,1}, {3,2,2,2}, {3,2,3,3}};
|
|
|
|
|
|
|
|
|
|
|
|
AllOperators={{OddllSLL,Op[7].Op[7]},
|
|
|
|
{OddllSRR,Op[6].Op[6]},
|
|
|
|
{OddllSRL,Op[6].Op[7]},
|
|
|
|
{OddllSLR,Op[7].Op[6]},
|
|
|
|
|
|
|
|
{OddllVRR,Op[7,Lor[1]].Op[7,Lor[1]]},
|
|
|
|
{OddllVLL,Op[6,Lor[1]].Op[6,Lor[1]]},
|
|
|
|
{OddllVRL,Op[7,Lor[1]].Op[6,Lor[1]]},
|
|
|
|
{OddllVLR,Op[6,Lor[1]].Op[7,Lor[1]]},
|
|
|
|
|
|
|
|
{OddllTLL,Op[-7,Lor[1],Lor[2]].Op[-7,Lor[1],Lor[2]]},
|
|
|
|
{OddllTLR,Op[-7,Lor[1],Lor[2]].Op[-6,Lor[1],Lor[2]]},
|
|
|
|
{OddllTRL,Op[-6,Lor[1],Lor[2]].Op[-7,Lor[1],Lor[2]]},
|
|
|
|
{OddllTRR,Op[-6,Lor[1],Lor[2]].Op[-6,Lor[1],Lor[2]]}
|
|
|
|
};
|
|
|
|
|
|
|
|
NameProcess="2d2nu";
|
|
|
|
|
|
|
|
ConsideredProcess = "4Fermion";
|
|
|
|
FermionOrderExternal={2,1,4,3};
|
|
|
|
NeglectMasses={1,2,3,4};
|
|
|
|
|
|
|
|
ExternalFields= {DownQuark,bar[DownQuark],Neutrino,bar[Neutrino]};
|
|
|
|
|
|
|
|
CombinationGenerations = Flatten[Table[{{2,1, neutrino1, neutrino2},
|
|
|
|
{3,1, neutrino1, neutrino2},{3,2, neutrino1, neutrino2}},
|
|
|
|
{neutrino1,1,3},{neutrino2,1,3}],2];
|
|
|
|
|
|
|
|
|
|
|
|
AllOperators={{OddvvVRR,Op[7,Lor[1]].Op[7,Lor[1]]},
|
|
|
|
{OddvvVLL,Op[6,Lor[1]].Op[6,Lor[1]]},
|
|
|
|
{OddvvVRL,Op[7,Lor[1]].Op[6,Lor[1]]},
|
|
|
|
{OddvvVLR,Op[6,Lor[1]].Op[7,Lor[1]]}
|
|
|
|
};
|
|
|
|
|
|
|
|
$(\\bar{\\ell} \\Gamma \\ell) (\\bar{d} \\Gamma^\\prime d)$ and$(\\bar{\\ell} \\Gamma \\ell) (\\bar{u} \\Gamma^\\prime u)$
|
|
|
|
|
|
|
|
| Variable | Operator | Name |
|
|
|
|
|----------|--------------------------------------------------------------------------------------|------------------------------------|
|
|
|
|
| OllddSLL | $(\\bar{\\ell} P_L \\ell) (\\bar{d} P_L d)$ | *B*<sub>*L**L*</sub><sup>*S*</sup> |
|
|
|
|
| OlluuSLL | $(\\bar{\\ell} P_L \\ell) (\\bar{u} P_L u)$ | *C*<sub>*L**L*</sub><sup>*S*</sup> |
|
|
|
|
| OllddSRR | $(\\bar{\\ell} P_R \\ell) (\\bar{d} P_R d)$ | *B*<sub>*R**R*</sub><sup>*S*</sup> |
|
|
|
|
| OlluuSRR | $(\\bar{\\ell} P_R \\ell) (\\bar{u} P_R u)$ | *C*<sub>*R**R*</sub><sup>*S*</sup> |
|
|
|
|
| OllddSRL | $(\\bar{\\ell} P_R \\ell) (\\bar{d} P_L d)$ | *B*<sub>*R**L*</sub><sup>*S*</sup> |
|
|
|
|
| OlluuSRL | $(\\bar{\\ell} P_R \\ell) (\\bar{u} P_L u)$ | *C*<sub>*R**L*</sub><sup>*S*</sup> |
|
|
|
|
| OllddSLR | $(\\bar{\\ell} P_L \\ell) (\\bar{d} P_R d)$ | *B*<sub>*L**R*</sub><sup>*S*</sup> |
|
|
|
|
| OlluuSLR | $(\\bar{\\ell} P_L \\ell) (\\bar{u} P_R u)$ | *C*<sub>*L**R*</sub><sup>*S*</sup> |
|
|
|
|
| OllddVLL | $(\\bar{\\ell} \\gamma_\\mu P_L \\ell) (\\bar{d} \\gamma^\\mu P_L d)$ | *B*<sub>*L**L*</sub><sup>*V*</sup> |
|
|
|
|
| OlluuVLL | $(\\bar{\\ell} \\gamma_\\mu P_L \\ell) (\\bar{u} \\gamma^\\mu P_L u)$ | *C*<sub>*L**L*</sub><sup>*V*</sup> |
|
|
|
|
| OllddVRR | $(\\bar{\\ell} \\gamma_\\mu P_R \\ell) (\\bar{d} \\gamma^\\mu P_R d)$ | *B*<sub>*R**R*</sub><sup>*V*</sup> |
|
|
|
|
| OlluuVRR | $(\\bar{\\ell} \\gamma_\\mu P_R \\ell) (\\bar{u} \\gamma^\\mu P_R u)$ | *C*<sub>*R**R*</sub><sup>*V*</sup> |
|
|
|
|
| OllddVLR | $(\\bar{\\ell} \\gamma_\\mu P_L \\ell) (\\bar{d} \\gamma^\\mu P_R d)$ | *B*<sub>*L**R*</sub><sup>*V*</sup> |
|
|
|
|
| OlluuVLR | $(\\bar{\\ell} \\gamma_\\mu P_L \\ell) (\\bar{u} \\gamma^\\mu P_R u)$ | *C*<sub>*L**R*</sub><sup>*V*</sup> |
|
|
|
|
| OllddVRL | $(\\bar{\\ell} \\gamma_\\mu P_R \\ell) (\\bar{d} \\gamma^\\mu P_L d)$ | *B*<sub>*R**L*</sub><sup>*V*</sup> |
|
|
|
|
| OlluuVRL | $(\\bar{\\ell} \\gamma_\\mu P_R \\ell) (\\bar{u} \\gamma^\\mu P_L u)$ | *C*<sub>*R**L*</sub><sup>*V*</sup> |
|
|
|
|
| OllddTLL | $(\\bar{\\ell} \\sigma_{\\mu\\nu} P_L \\ell) (\\bar{d} \\sigma^{\\mu\\nu} P_L d)$ | *B*<sub>*L**L*</sub><sup>*T*</sup> |
|
|
|
|
| OlluuTLL | $(\\bar{\\ell} \\sigma_{\\mu\\nu} P_L \\ell) (\\bar{u} \\sigma^{\\mu\\nu} P_L u)$ | *C*<sub>*L**L*</sub><sup>*T*</sup> |
|
|
|
|
| OllddTRR | $(\\bar{\\ell} \\sigma_{\\mu\\nu} P_R \\ell) (\\bar{d} \\sigma^{\\mu\\nu} P_R d)$ | *B*<sub>*R**R*</sub><sup>*T*</sup> |
|
|
|
|
| OlluuTRR | $(\\bar{\\ell} \\sigma_{\\mu\\nu} P_R \\ell) (\\bar{u} \\sigma^{\\mu\\nu} P_R u)$ | *C*<sub>*R**R*</sub><sup>*T*</sup> |
|
|
|
|
| OllddTLR | $(\\bar{\\ell} \\sigma_{\\mu\\nu} P_L \\ell) (\\bar{d} \\sigma^{\\mu\\nu} P_R d)$ | *B*<sub>*L**R*</sub><sup>*T*</sup> |
|
|
|
|
| OlluuTLR | $(\\bar{\\ell} \\sigma_{\\mu\\nu} P_L \\ell) (\\bar{u} \\sigma^{\\mu\\nu} P_R u)$ | *C*<sub>*L**R*</sub><sup>*T*</sup> |
|
|
|
|
| OllddTRL | $(\\bar{\\ell} \\sigma_{\\mu\\nu} P_R \\ell) (\\bar{d} \\sigma^{\\mu\\nu} P_L d)$ | *B*<sub>*R**L*</sub><sup>*T*</sup> |
|
|
|
|
| OlluuTRL | $(\\bar{\\ell} \\sigma_{\\mu\\nu} P_R \\ell) (\\bar{u} \\sigma^{\\mu\\nu} P_L u)$ | *C*<sub>*R**L*</sub><sup>*T*</sup> |
|
|
|
|
|
|
|
|
NameProcess="2L2d";
|
|
|
|
|
|
|
|
ConsideredProcess = "4Fermion";
|
|
|
|
FermionOrderExternal={2,1,4,3};
|
|
|
|
NeglectMasses={1,2,3,4};
|
|
|
|
|
|
|
|
|
|
|
|
ExternalFields= {ChargedLepton,bar[ChargedLepton],DownQuark,bar[DownQuark]};
|
|
|
|
CombinationGenerations = {{2,1,1,1}, {3,1,1,1}, {3,2,1,1},
|
|
|
|
{2,1,2,2}, {3,1,2,2}, {3,2,2,2}};
|
|
|
|
|
|
|
|
|
|
|
|
AllOperators={{OllddSLL,Op[7].Op[7]},
|
|
|
|
{OllddSRR,Op[6].Op[6]},
|
|
|
|
{OllddSRL,Op[6].Op[7]},
|
|
|
|
{OllddSLR,Op[7].Op[6]},
|
|
|
|
|
|
|
|
{OllddVRR,Op[7,Lor[1]].Op[7,Lor[1]]},
|
|
|
|
{OllddVLL,Op[6,Lor[1]].Op[6,Lor[1]]},
|
|
|
|
{OllddVRL,Op[7,Lor[1]].Op[6,Lor[1]]},
|
|
|
|
{OllddVLR,Op[6,Lor[1]].Op[7,Lor[1]]},
|
|
|
|
|
|
|
|
{OllddTLL,Op[-7,Lor[1],Lor[2]].Op[-7,Lor[1],Lor[2]]},
|
|
|
|
{OllddTLR,Op[-7,Lor[1],Lor[2]].Op[-6,Lor[1],Lor[2]]},
|
|
|
|
{OllddTRL,Op[-6,Lor[1],Lor[2]].Op[-7,Lor[1],Lor[2]]},
|
|
|
|
{OllddTRR,Op[-6,Lor[1],Lor[2]].Op[-6,Lor[1],Lor[2]]}
|
|
|
|
};
|
|
|
|
|
|
|
|
NameProcess="2L2u";
|
|
|
|
|
|
|
|
ConsideredProcess = "4Fermion";
|
|
|
|
FermionOrderExternal={2,1,4,3};
|
|
|
|
NeglectMasses={1,2,3,4};
|
|
|
|
|
|
|
|
|
|
|
|
ExternalFields= {ChargedLepton,bar[ChargedLepton],UpQuark,bar[UpQuark]};
|
|
|
|
CombinationGenerations = {{2,1,1,1},{3,1,1,1},{3,2,1,1}};
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
AllOperators={{OlluuSLL,Op[7].Op[7]},
|
|
|
|
{OlluuSRR,Op[6].Op[6]},
|
|
|
|
{OlluuSRL,Op[6].Op[7]},
|
|
|
|
{OlluuSLR,Op[7].Op[6]},
|
|
|
|
|
|
|
|
{OlluuVRR,Op[7,Lor[1]].Op[7,Lor[1]]},
|
|
|
|
{OlluuVLL,Op[6,Lor[1]].Op[6,Lor[1]]},
|
|
|
|
{OlluuVRL,Op[7,Lor[1]].Op[6,Lor[1]]},
|
|
|
|
{OlluuVLR,Op[6,Lor[1]].Op[7,Lor[1]]},
|
|
|
|
|
|
|
|
{OlluuTLL,Op[-7,Lor[1],Lor[2]].Op[-7,Lor[1],Lor[2]]},
|
|
|
|
{OlluuTLR,Op[-7,Lor[1],Lor[2]].Op[-6,Lor[1],Lor[2]]},
|
|
|
|
{OlluuTRL,Op[-6,Lor[1],Lor[2]].Op[-7,Lor[1],Lor[2]]},
|
|
|
|
{OlluuTRR,Op[-6,Lor[1],Lor[2]].Op[-6,Lor[1],Lor[2]]}
|
|
|
|
};
|
|
|
|
|
|
|
|
(*d̄**Γ**d*)(*d̄**Γ*<sup>′</sup>*d*) and$(\\bar{\\ell} \\Gamma \\ell) (\\bar{\\ell} \\Gamma^\\prime \\ell)$
|
|
|
|
|
|
|
|
| Variable | Operator | Name |
|
|
|
|
|----------|----------------------------------------------------------------------------------------------|------------------------------------|
|
|
|
|
| O4dSLL | (*d̄**P*<sub>*L*</sub>*d*)(*d̄**P*<sub>*L*</sub>*d*) | *D*<sub>*L**L*</sub><sup>*S*</sup> |
|
|
|
|
| O4lSLL | $(\\bar{\\ell} P_L \\ell) (\\bar{\\ell} P_L \\ell)$ | *A*<sub>*L**L*</sub><sup>*S*</sup> |
|
|
|
|
| O4dSRR | (*d̄**P*<sub>*R*</sub>*d*)(*d̄**P*<sub>*R*</sub>*d*) | *D*<sub>*R**R*</sub><sup>*S*</sup> |
|
|
|
|
| O4lSRR | $(\\bar{\\ell} P_R \\ell) (\\bar{\\ell} P_R \\ell)$ | *A*<sub>*R**R*</sub><sup>*S*</sup> |
|
|
|
|
| O4dSLR | (*d̄**P*<sub>*L*</sub>*d*)(*d̄**P*<sub>*R*</sub>*d*) | *D*<sub>*L**R*</sub><sup>*S*</sup> |
|
|
|
|
| O4lSLR | $(\\bar{\\ell} P_L \\ell) (\\bar{\\ell} P_R \\ell)$ | *A*<sub>*L**R*</sub><sup>*S*</sup> |
|
|
|
|
| O4dSRL | (*d̄**P*<sub>*R*</sub>*d*)(*d̄**P*<sub>*L*</sub>*d*) | *D*<sub>*R**L*</sub><sup>*S*</sup> |
|
|
|
|
| O4lSRL | $(\\bar{\\ell} P_R \\ell) (\\bar{\\ell} P_L \\ell)$ | *A*<sub>*R**L*</sub><sup>*S*</sup> |
|
|
|
|
| O4dVLL | (*d̄**γ*<sub>*μ*</sub>*P*<sub>*L*</sub>*d*)(*d̄**γ*<sup>*μ*</sup>*P*<sub>*L*</sub>*d*) | *D*<sub>*L**L*</sub><sup>*V*</sup> |
|
|
|
|
| O4lVLL | $(\\bar{\\ell} \\gamma_\\mu P_L \\ell) (\\bar{\\ell} \\gamma^\\mu P_L \\ell)$ | *A*<sub>*L**L*</sub><sup>*V*</sup> |
|
|
|
|
| O4dVRR | (*d̄**γ*<sub>*μ*</sub>*P*<sub>*R*</sub>*d*)(*d̄**γ*<sup>*μ*</sup>*P*<sub>*R*</sub>*d*) | *D*<sub>*R**R*</sub><sup>*V*</sup> |
|
|
|
|
| O4lVRR | $(\\bar{\\ell} \\gamma_\\mu P_R \\ell) (\\bar{\\ell} \\gamma^\\mu P_R \\ell)$ | *A*<sub>*R**R*</sub><sup>*V*</sup> |
|
|
|
|
| O4dVLR | (*d̄**γ*<sub>*μ*</sub>*P*<sub>*L*</sub>*d*)(*d̄**γ*<sup>*μ*</sup>*P*<sub>*R*</sub>*d*) | *D*<sub>*L**R*</sub><sup>*V*</sup> |
|
|
|
|
| O4lVLR | $(\\bar{\\ell} \\gamma_\\mu P_L \\ell) (\\bar{\\ell} \\gamma^\\mu P_R \\ell)$ | *A*<sub>*L**R*</sub><sup>*V*</sup> |
|
|
|
|
| O4dVRL | (*d̄**γ*<sub>*μ*</sub>*P*<sub>*R*</sub>*d*)(*d̄**γ*<sup>*μ*</sup>*P*<sub>*L*</sub>*d*) | *D*<sub>*R**L*</sub><sup>*V*</sup> |
|
|
|
|
| O4lVRL | $(\\bar{\\ell} \\gamma_\\mu P_R \\ell) (\\bar{\\ell} \\gamma^\\mu P_L \\ell)$ | *A*<sub>*R**L*</sub><sup>*V*</sup> |
|
|
|
|
| O4dTLL | (*d̄**σ*<sub>*μ**ν*</sub>*P*<sub>*L*</sub>*d*)(*d̄**σ*<sup>*μ**ν*</sup>*P*<sub>*L*</sub>*d*) | *D*<sub>*L**L*</sub><sup>*T*</sup> |
|
|
|
|
| O4lTLL | $(\\bar{\\ell} \\sigma_{\\mu\\nu} P_L \\ell) (\\bar{\\ell} \\sigma^{\\mu\\nu} P_L \\ell)$ | *A*<sub>*L**L*</sub><sup>*T*</sup> |
|
|
|
|
| O4dTRR | (*d̄**σ*<sub>*μ**ν*</sub>*P*<sub>*R*</sub>*d*)(*d̄**σ*<sup>*μ**ν*</sup>*P*<sub>*R*</sub>*d*) | *D*<sub>*R**R*</sub><sup>*T*</sup> |
|
|
|
|
| O4lTRR | $(\\bar{\\ell} \\sigma_{\\mu\\nu} P_R \\ell) (\\bar{\\ell} \\sigma^{\\mu\\nu} P_R \\ell)$ | *A*<sub>*R**R*</sub><sup>*T*</sup> |
|
|
|
|
| O4dTLR | (*d̄**σ*<sub>*μ**ν*</sub>*P*<sub>*L*</sub>*d*)(*d̄**σ*<sup>*μ**ν*</sup>*P*<sub>*R*</sub>*d*) | *D*<sub>*L**R*</sub><sup>*T*</sup> |
|
|
|
|
| O4lTLR | $(\\bar{\\ell} \\sigma_{\\mu\\nu} P_L \\ell) (\\bar{\\ell} \\sigma^{\\mu\\nu} P_R \\ell)$ | *A*<sub>*L**R*</sub><sup>*T*</sup> |
|
|
|
|
| O4dTRL | (*d̄**σ*<sub>*μ**ν*</sub>*P*<sub>*R*</sub>*d*)(*d̄**σ*<sup>*μ**ν*</sup>*P*<sub>*L*</sub>*d*) | *D*<sub>*R**L*</sub><sup>*T*</sup> |
|
|
|
|
| O4lTRL | $(\\bar{\\ell} \\sigma_{\\mu\\nu} P_R \\ell) (\\bar{\\ell} \\sigma^{\\mu\\nu} P_L \\ell)$ | *A*<sub>*R**L*</sub><sup>*T*</sup> |
|
|
|
|
|
|
|
|
NameProcess="4d";
|
|
|
|
|
|
|
|
ConsideredProcess = "4Fermion";
|
|
|
|
FermionOrderExternal={2,1,4,3};
|
|
|
|
NeglectMasses={1,2,3,4};
|
|
|
|
|
|
|
|
|
|
|
|
ExternalFields= {DownQuark,bar[DownQuark],DownQuark,bar[DownQuark]};
|
|
|
|
|
|
|
|
ColorFlow = ColorDelta[1,2] ColorDelta[3,4];
|
|
|
|
|
|
|
|
CombinationGenerations = {{3,1,3,1},{3,2,3,2},{2,1,2,1}};
|
|
|
|
|
|
|
|
|
|
|
|
AllOperators={{O4dSLL,Op[7].Op[7]},
|
|
|
|
{O4dSRR,Op[6].Op[6]},
|
|
|
|
{O4dSRL,Op[6].Op[7]},
|
|
|
|
{O4dSLR,Op[7].Op[6]},
|
|
|
|
|
|
|
|
{O4dVRR,Op[7,Lor[1]].Op[7,Lor[1]]},
|
|
|
|
{O4dVLL,Op[6,Lor[1]].Op[6,Lor[1]]},
|
|
|
|
{O4dVRL,Op[7,Lor[1]].Op[6,Lor[1]]},
|
|
|
|
{O4dVLR,Op[6,Lor[1]].Op[7,Lor[1]]},
|
|
|
|
|
|
|
|
{O4dTLL,Op[-7,Lor[1],Lor[2]].Op[-7,Lor[1],Lor[2]]},
|
|
|
|
{O4dTLR,Op[-7,Lor[1],Lor[2]].Op[-6,Lor[1],Lor[2]]},
|
|
|
|
{O4dTRL,Op[-6,Lor[1],Lor[2]].Op[-7,Lor[1],Lor[2]]},
|
|
|
|
{O4dTRR,Op[-6,Lor[1],Lor[2]].Op[-6,Lor[1],Lor[2]]}
|
|
|
|
};
|
|
|
|
|
|
|
|
Filters = {NoPenguins};
|
|
|
|
|
|
|
|
NameProcess="4L";
|
|
|
|
|
|
|
|
ConsideredProcess = "4Fermion";
|
|
|
|
FermionOrderExternal={2,1,4,3};
|
|
|
|
NeglectMasses={1,2,3,4};
|
|
|
|
|
|
|
|
ExternalFields= {ChargedLepton,bar[ChargedLepton],ChargedLepton,bar[ChargedLepton]};
|
|
|
|
CombinationGenerations = {{2,1,1,1},{3,1,1,1},{3,2,2,2}};
|
|
|
|
|
|
|
|
|
|
|
|
AllOperators={{O4lSLL,Op[7].Op[7]},
|
|
|
|
{O4lSRR,Op[6].Op[6]},
|
|
|
|
{O4lSRL,Op[6].Op[7]},
|
|
|
|
{O4lSLR,Op[7].Op[6]},
|
|
|
|
|
|
|
|
{O4lVRR,Op[7,Lor[1]].Op[7,Lor[1]]},
|
|
|
|
{O4lVLL,Op[6,Lor[1]].Op[6,Lor[1]]},
|
|
|
|
{O4lVRL,Op[7,Lor[1]].Op[6,Lor[1]]},
|
|
|
|
{O4lVLR,Op[6,Lor[1]].Op[7,Lor[1]]},
|
|
|
|
|
|
|
|
{O4lTLL,Op[-7,Lor[1],Lor[2]].Op[-7,Lor[1],Lor[2]]},
|
|
|
|
{O4lTLR,Op[-7,Lor[1],Lor[2]].Op[-6,Lor[1],Lor[2]]},
|
|
|
|
{O4lTRL,Op[-6,Lor[1],Lor[2]].Op[-7,Lor[1],Lor[2]]},
|
|
|
|
{O4lTRR,Op[-6,Lor[1],Lor[2]].Op[-6,Lor[1],Lor[2]]}
|
|
|
|
};
|
|
|
|
|
|
|
|
Filters = {NoCrossedDiagrams};
|
|
|
|
|
|
|
|
$(\\bar{d} \\Gamma u) (\\bar{\\ell} \\Gamma^\\prime \\nu)$
|
|
|
|
|
|
|
|
| Variable | Operator | Name |
|
|
|
|
|----------|-------------------------------------------------------------------------|------------------------------------|
|
|
|
|
| OdulvVLL | $(\\bar{d} \\gamma_\\mu P_L u) (\\bar{\\ell} \\gamma^\\mu P_L \\nu)$ | *G*<sub>*L**L*</sub><sup>*V*</sup> |
|
|
|
|
| OdulvSLL | $(\\bar{d} P_L u) (\\bar{\\ell} P_L \\nu)$ | *G*<sub>*L**L*</sub><sup>*S*</sup> |
|
|
|
|
| OdulvVRR | $(\\bar{d} \\gamma_\\mu P_R u) (\\bar{\\ell} \\gamma^\\mu P_R \\nu)$ | *G*<sub>*R**R*</sub><sup>*V*</sup> |
|
|
|
|
| OdulvSRR | $(\\bar{d} P_R u) (\\bar{\\ell} P_R \\nu)$ | *G*<sub>*R**R*</sub><sup>*S*</sup> |
|
|
|
|
| OdulvVLR | $(\\bar{d} \\gamma_\\mu P_L u) (\\bar{\\ell} \\gamma^\\mu P_R \\nu)$ | *G*<sub>*L**R*</sub><sup>*V*</sup> |
|
|
|
|
| OdulvSLR | $(\\bar{d} P_L u) (\\bar{\\ell} P_R \\nu)$ | *G*<sub>*L**R*</sub><sup>*S*</sup> |
|
|
|
|
| OdulvVRL | $(\\bar{d} \\gamma_\\mu P_R u) (\\bar{\\ell} \\gamma^\\mu P_L \\nu)$ | *G*<sub>*R**L*</sub><sup>*V*</sup> |
|
|
|
|
| OdulvSRL | $(\\bar{d} P_R u) (\\bar{\\ell} P_L \\nu)$ | *G*<sub>*R**L*</sub><sup>*S*</sup> |
|
|
|
|
|
|
|
|
NameProcess="dulv";
|
|
|
|
|
|
|
|
ConsideredProcess = "4Fermion";
|
|
|
|
FermionOrderExternal={2,1,3,4};
|
|
|
|
NeglectMasses={1,2,3,4};
|
|
|
|
|
|
|
|
|
|
|
|
ExternalFields= {DownQuark,bar[UpQuark], Neutrino, bar[ChargedLepton]};
|
|
|
|
|
|
|
|
CombinationGenerations =
|
|
|
|
Flatten[Table[{{3,1,i,j},{3,2,i,j},{2,2,i,j},{2,1,i,j}},{i,1,3},{j,1,3}],2];
|
|
|
|
|
|
|
|
Clear[i,j];
|
|
|
|
|
|
|
|
|
|
|
|
AllOperators={{OdulvSLL,Op[7].Op[7]},
|
|
|
|
{OdulvSRR,Op[6].Op[6]},
|
|
|
|
{OdulvSRL,Op[6].Op[7]},
|
|
|
|
{OdulvSLR,Op[7].Op[6]},
|
|
|
|
|
|
|
|
{OdulvVRR,Op[7,Lor[1]].Op[7,Lor[1]]},
|
|
|
|
{OdulvVLL,Op[6,Lor[1]].Op[6,Lor[1]]},
|
|
|
|
{OdulvVRL,Op[7,Lor[1]].Op[6,Lor[1]]},
|
|
|
|
{OdulvVLR,Op[6,Lor[1]].Op[7,Lor[1]]}
|
|
|
|
};
|
|
|
|
|
|
|
|
Filters = {NoBoxes, NoPenguins};
|
|
|
|
|
|
|
|
See also
|
|
|
|
========
|
|
|
|
|
|
|
|
[1] We would like to emphasize that our implementation of these operators is only valid for virtual scalars and pseudoscalars. They have been introduced in order to provide the 1-loop vertices necessary for the computation of the double penguin contributions to*Δ**M*<sub>*B*<sub>*q*</sub></sub>. Therefore, they are not valid for observables in which the scalar or pseudoscalar states are real particles.
|
|
|
|
|
|
|
|
[2] The(*d̄**Γ**d*)*S* and(*d̄**Γ**d*)*P* operators have been introduced to compute double penguin corrections to*Δ**M*<sub>*B*<sub>*q*</sub></sub>, where*S* and*P* appear as intermediate (virtual) particles. They should not be used in processes where the scalar or pseudoscalar states are real particles because the loop functions are calculated with vanishing external momenta. |
|
|
|
\ No newline at end of file |