Already_defined_Operators_in_FlavorKit

Already defined Operators in FlavorKit

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

{\mathcal L}_{\text{FV}} = {\mathcal L}_{\text{LFV}} + {\mathcal 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

{\mathcal L}_{\text{LFV}} = {\mathcal L}_{\ell \ell \gamma} + {\mathcal L}_{\ell \ell Z} + {\mathcal L}_{\ell \ell h} + {\mathcal L}_{4 \ell} + {\mathcal L}_{2 \ell 2q} \, .

The first term contains theℓ − ℓ − γ interaction, given by

{\mathcal 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.

Heree is the electric charge,q the photon momentum,$P_{L,R} = \frac{1}{2} (1 \mp \gamma_5)$ are the usual chirality projectors andℓ_α, β 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, theZ- and Higgs boson contributions will be included whenever possible. Therefore, theℓ − ℓ − Z andℓ − ℓ − h interaction Lagrangians will only be used for observables involving realZ- and Higgs bosons. These two Lagrangians can be written as

{\mathcal 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 \, ,

wherep is theℓ_β 4-momentum, and

{\mathcal 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

{\mathcal 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. \, ,

whereℓ_{α, β, γ, δ} denote the lepton flavors andΓ_S = 1,Γ_V = γ_μ andΓ_T = σ_μ**ν. 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 coefficientsA_L**R^S andA_R**L^S 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ℓ2q four fermion interaction Lagrangian at the quark level is given by

{\mathcal L}_{2 \ell 2q} = {\mathcal L}_{2 \ell 2d} + {\mathcal L}_{2 \ell 2u}

where

{\mathcal 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. \\ {\mathcal L}_{2 \ell 2u} = \left. {\mathcal L}_{2 \ell 2d} \right|_{d \to u, \, B \to C} \, .

Hered_γ denotes the d-quark flavor.

Let us now consider the Lagrangian relevant for quark flavor violation. This can be written as

{\mathcal L}_{\text{QFV}} = {\mathcal L}_{q q \gamma} + {\mathcal L}_{q q g} + {\mathcal L}_{4 d} + {\mathcal L}_{2d2l} + {\mathcal L}_{2d2\nu} + {\mathcal L}_{du\ell\nu} + {\mathcal L}_{d d H} \, .

The first two terms correspond to operators that couple quark bilinears to massless gauge bosons. These are

{\mathcal 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} \\ {\mathcal 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} \, .

HereT^a areS**U(3) matrices. The Wilson coefficientsQ_1, 2^L, R can be easily related to the usualC_7, 8^(′) coefficients, sometimes normalized with an additional\frac{1}{16 \pi^2} factor. The4d four fermion interaction Lagrangian can be written as

{\mathcal 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. \, ,

whered_{α, β, γ, δ} denote the lepton flavors. Again, we omit flavor indices in the Wilson coefficients for the sake of clarity. The2d2ℓ four fermion interaction Lagrangian is given by

{\mathcal 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 \, .

Hereℓ_γ denotes the lepton flavor.{\mathcal L}_{2d 2 \ell} should not be confused with{\mathcal 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. The2d2ν andd**uℓν terms of the QFV Lagrangian are

{\mathcal 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 \\ {\mathcal 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 \, .

Note that we have not introduced scalar or tensor2d2ν operators, nor tensord**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, whereS (P) is a virtual [1] scalar (pseudoscalar) state. This piece can be written as

{\mathcal 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_g.n_g is the number of generations of the considered scalar and forn_g = 1 the last index is dropped.

(d̄_βσ_μ**νΓ**d_α)F^μ**ν and(d̄_βσ_μ**νΓ**d_α)G^μ**ν

Variable	Operator	Name
CC7	em_{d_β}(d̄_βσ_μνP_Ld_α)F^μν	Q₁^L
CC7p	em_{d_α}(d̄_βσ_μνP_Rd_α)F^μν	Q₁^R
CC8	g_sm_{d_β}(d̄_βσ_μνP_Ld_α)G^μν	Q₂^L
CC8p	g_sm_{d_α}(d̄_βσ_μνP_Rd_α)G^μν	Q₂^R

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₂^L
K1L	`q^2 (\bar \ell_\beta \gamma_\mu P_L \ell_\alpha) A^\mu`	K₁^L
K2R	`e m_{\ell_\alpha} (\bar \ell_\beta \sigma_{\mu\nu} P_R \ell_\alpha) q^{\nu} A^\mu`	K₂^L
K1R	`q^2 (\bar \ell_\beta \gamma_\nu P_R \ell_\alpha) A^\mu`	K₁^R

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₁^L
OZ2lSL	`(\bar{\ell} p^\mu P_L \ell) Z_\mu`	R₂^L
OZ2lVR	`(\bar{\ell} \, \gamma^\mu P_R \ell) Z_\mu`	R₁^R
OZ2lSR	`(\bar{\ell} p^\mu P_R \ell) Z_\mu`	R₂^R

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_L
OH2lSR	`\bar{\ell} P_R \ell \, h`	S_R

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_Ld S	H_L^S
OH2qSR	d̄P_Rd S	H_R^S
OAh2qSL	d̄P_Ld P	H_L^P
OAh2qSR	d̄P_Rd P	H_R^P

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 andd**dℓℓ for QFV.

(\bar{d} \Gamma d) (\bar{\ell} \Gamma^\prime \ell) and(d̄Γd)(ν̄**Γ^′ν)

Variable	Operator	Name
OddllSLL	`(\bar{d} P_L d) (\bar{\ell} P_L \ell)`	E_L**L^S
OddllSRR	`(\bar{d} P_R d) (\bar{\ell} P_R \ell)`	E_R**R^S
OddllSLR	`(\bar{d} P_L d) (\bar{\ell} P_R \ell)`	E_L**R^S
OddllSRL	`(\bar{d} P_R d) (\bar{\ell} P_L \ell)`	E_R**L^S
OddllVLL	`(\bar{d} \gamma_\mu P_L d) (\bar{\ell} \gamma^\mu P_L \ell)`	E_L**L^V
OddvvVLL	(d̄γ_μP_Ld)(ν̄γ^μP_Rν)	F_L**L^V
OddllVRR	`(\bar{d} \gamma_\mu P_R d) (\bar{\ell} \gamma^\mu P_R \ell)`	E_R**R^V
OddvvVRR	(d̄γ_μP_Rd)(ν̄γ^μP_Rν)	F_R**R^V
OddllVLR	`(\bar{d} \gamma_\mu P_L d) (\bar{\ell} \gamma^\mu P_R \ell)`	E_L**R^V
OddvvVLR	(d̄γ_μP_Ld)(ν̄γ^μP_Rν)	F_L**R^V
OddllVRL	`(\bar{d} \gamma_\mu P_R d) (\bar{\ell} \gamma^\mu P_L \ell)`	E_R**L^V
OddvvVRL	(d̄γ_μP_Rd)(ν̄γ^μP_Lν)	F_R**L^V
OddllTLL	`(\bar{d} \sigma_{\mu\nu} P_L d) (\bar{\ell} \sigma^{\mu\nu} P_L \ell)`	E_L**L^T
OddllTRR	`(\bar{d} \sigma_{\mu\nu} P_R d) (\bar{\ell} \sigma^{\mu\nu} P_R \ell)`	E_R**R^T
OddllTLR	`(\bar{d} \sigma_{\mu\nu} P_L d) (\bar{\ell} \sigma^{\mu\nu} P_R \ell)`	E_L**R^T
OddllTRL	`(\bar{d} \sigma_{\mu\nu} P_R d) (\bar{\ell} \sigma^{\mu\nu} P_L \ell)`	E_R**L^T

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_L**L^S
OlluuSLL	`(\bar{\ell} P_L \ell) (\bar{u} P_L u)`	C_L**L^S
OllddSRR	`(\bar{\ell} P_R \ell) (\bar{d} P_R d)`	B_R**R^S
OlluuSRR	`(\bar{\ell} P_R \ell) (\bar{u} P_R u)`	C_R**R^S
OllddSRL	`(\bar{\ell} P_R \ell) (\bar{d} P_L d)`	B_R**L^S
OlluuSRL	`(\bar{\ell} P_R \ell) (\bar{u} P_L u)`	C_R**L^S
OllddSLR	`(\bar{\ell} P_L \ell) (\bar{d} P_R d)`	B_L**R^S
OlluuSLR	`(\bar{\ell} P_L \ell) (\bar{u} P_R u)`	C_L**R^S
OllddVLL	`(\bar{\ell} \gamma_\mu P_L \ell) (\bar{d} \gamma^\mu P_L d)`	B_L**L^V
OlluuVLL	`(\bar{\ell} \gamma_\mu P_L \ell) (\bar{u} \gamma^\mu P_L u)`	C_L**L^V
OllddVRR	`(\bar{\ell} \gamma_\mu P_R \ell) (\bar{d} \gamma^\mu P_R d)`	B_R**R^V
OlluuVRR	`(\bar{\ell} \gamma_\mu P_R \ell) (\bar{u} \gamma^\mu P_R u)`	C_R**R^V
OllddVLR	`(\bar{\ell} \gamma_\mu P_L \ell) (\bar{d} \gamma^\mu P_R d)`	B_L**R^V
OlluuVLR	`(\bar{\ell} \gamma_\mu P_L \ell) (\bar{u} \gamma^\mu P_R u)`	C_L**R^V
OllddVRL	`(\bar{\ell} \gamma_\mu P_R \ell) (\bar{d} \gamma^\mu P_L d)`	B_R**L^V
OlluuVRL	`(\bar{\ell} \gamma_\mu P_R \ell) (\bar{u} \gamma^\mu P_L u)`	C_R**L^V
OllddTLL	`(\bar{\ell} \sigma_{\mu\nu} P_L \ell) (\bar{d} \sigma^{\mu\nu} P_L d)`	B_L**L^T
OlluuTLL	`(\bar{\ell} \sigma_{\mu\nu} P_L \ell) (\bar{u} \sigma^{\mu\nu} P_L u)`	C_L**L^T
OllddTRR	`(\bar{\ell} \sigma_{\mu\nu} P_R \ell) (\bar{d} \sigma^{\mu\nu} P_R d)`	B_R**R^T
OlluuTRR	`(\bar{\ell} \sigma_{\mu\nu} P_R \ell) (\bar{u} \sigma^{\mu\nu} P_R u)`	C_R**R^T
OllddTLR	`(\bar{\ell} \sigma_{\mu\nu} P_L \ell) (\bar{d} \sigma^{\mu\nu} P_R d)`	B_L**R^T
OlluuTLR	`(\bar{\ell} \sigma_{\mu\nu} P_L \ell) (\bar{u} \sigma^{\mu\nu} P_R u)`	C_L**R^T
OllddTRL	`(\bar{\ell} \sigma_{\mu\nu} P_R \ell) (\bar{d} \sigma^{\mu\nu} P_L d)`	B_R**L^T
OlluuTRL	`(\bar{\ell} \sigma_{\mu\nu} P_R \ell) (\bar{u} \sigma^{\mu\nu} P_L u)`	C_R**L^T

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̄**Γ^′d) and(\bar{\ell} \Gamma \ell) (\bar{\ell} \Gamma^\prime \ell)

Variable	Operator	Name
O4dSLL	(d̄P_Ld)(d̄P_Ld)	D_L**L^S
O4lSLL	`(\bar{\ell} P_L \ell) (\bar{\ell} P_L \ell)`	A_L**L^S
O4dSRR	(d̄P_Rd)(d̄P_Rd)	D_R**R^S
O4lSRR	`(\bar{\ell} P_R \ell) (\bar{\ell} P_R \ell)`	A_R**R^S
O4dSLR	(d̄P_Ld)(d̄P_Rd)	D_L**R^S
O4lSLR	`(\bar{\ell} P_L \ell) (\bar{\ell} P_R \ell)`	A_L**R^S
O4dSRL	(d̄P_Rd)(d̄P_Ld)	D_R**L^S
O4lSRL	`(\bar{\ell} P_R \ell) (\bar{\ell} P_L \ell)`	A_R**L^S
O4dVLL	(d̄γ_μP_Ld)(d̄γ^μP_Ld)	D_L**L^V
O4lVLL	`(\bar{\ell} \gamma_\mu P_L \ell) (\bar{\ell} \gamma^\mu P_L \ell)`	A_L**L^V
O4dVRR	(d̄γ_μP_Rd)(d̄γ^μP_Rd)	D_R**R^V
O4lVRR	`(\bar{\ell} \gamma_\mu P_R \ell) (\bar{\ell} \gamma^\mu P_R \ell)`	A_R**R^V
O4dVLR	(d̄γ_μP_Ld)(d̄γ^μP_Rd)	D_L**R^V
O4lVLR	`(\bar{\ell} \gamma_\mu P_L \ell) (\bar{\ell} \gamma^\mu P_R \ell)`	A_L**R^V
O4dVRL	(d̄γ_μP_Rd)(d̄γ^μP_Ld)	D_R**L^V
O4lVRL	`(\bar{\ell} \gamma_\mu P_R \ell) (\bar{\ell} \gamma^\mu P_L \ell)`	A_R**L^V
O4dTLL	(d̄σ_μ*νP_Ld)(d̄*σ^μ**νP_Ld)	D_L**L^T
O4lTLL	`(\bar{\ell} \sigma_{\mu\nu} P_L \ell) (\bar{\ell} \sigma^{\mu\nu} P_L \ell)`	A_L**L^T
O4dTRR	(d̄σ_μ*νP_Rd)(d̄*σ^μ**νP_Rd)	D_R**R^T
O4lTRR	`(\bar{\ell} \sigma_{\mu\nu} P_R \ell) (\bar{\ell} \sigma^{\mu\nu} P_R \ell)`	A_R**R^T
O4dTLR	(d̄σ_μ*νP_Ld)(d̄*σ^μ**νP_Rd)	D_L**R^T
O4lTLR	`(\bar{\ell} \sigma_{\mu\nu} P_L \ell) (\bar{\ell} \sigma^{\mu\nu} P_R \ell)`	A_L**R^T
O4dTRL	(d̄σ_μ*νP_Rd)(d̄*σ^μ**νP_Ld)	D_R**L^T
O4lTRL	`(\bar{\ell} \sigma_{\mu\nu} P_R \ell) (\bar{\ell} \sigma^{\mu\nu} P_L \ell)`	A_R**L^T

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_L**L^V
OdulvSLL	`(\bar{d} P_L u) (\bar{\ell} P_L \nu)`	G_L**L^S
OdulvVRR	`(\bar{d} \gamma_\mu P_R u) (\bar{\ell} \gamma^\mu P_R \nu)`	G_R**R^V
OdulvSRR	`(\bar{d} P_R u) (\bar{\ell} P_R \nu)`	G_R**R^S
OdulvVLR	`(\bar{d} \gamma_\mu P_L u) (\bar{\ell} \gamma^\mu P_R \nu)`	G_L**R^V
OdulvSLR	`(\bar{d} P_L u) (\bar{\ell} P_R \nu)`	G_L**R^S
OdulvVRL	`(\bar{d} \gamma_\mu P_R u) (\bar{\ell} \gamma^\mu P_L \nu)`	G_R**L^V
OdulvSRL	`(\bar{d} P_R u) (\bar{\ell} P_L \nu)`	G_R**L^S

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};

GitLab