TopHeavyRes: 4topEFT.fr

File 4topEFT.fr, 6.7 KB (added by LucDarme, 3 years ago)
Line 
1(***************************************************************************************************************)
2(****** FeynRules for pure four top EFT production ******)
3(****** ******)
4(****** Authors: By L. Darme (for EW-breaking part) ******)
5(****** EW-perserving part derived from C. Degrande (https://feynrules.irmp.ucl.ac.be/wiki/TopEffTh ******)
6(****** ******)
7(***************************************************************************************************************)
8
9M$ModelName = "4topEFT";
10
11
12M$Information = {Authors -> {" L. Darme"},
13 Version -> "1.0",
14 Date -> "26.02. 2021",
15 Institutions -> {"INFN"},
16 Emails -> {"luc.darme@lnf.infn.it"}};
17
18
19 M$InteractionOrderHierarchy = { {QCD, 1}, {QED, 2}, {NP,1} };
20 M$InteractionOrderLimit = { {QCD, 99}, {QED, 99}, {NP,99} };
21
22
23
24(**************** Parameters *************)
25
26M$Parameters = {
27
28(* External parameters *)
29
30LamNP == {
31 ParameterType -> External,
32 BlockName -> DIM6,
33 Value -> 1,
34 TeX -> \[CapitalLambda],
35 Description -> "cHigh scale EFT "},
36
37 Cuu == {
38 ParameterType -> External,
39 BlockName -> DIM6,
40 InteractionOrder -> { NP, 2},
41 Value -> 1,
42 TeX -> Subscript[C,uu],
43 Description -> "coefficient of Ouu "},
44
45 Cqq== {
46 ParameterType -> External,
47 BlockName -> DIM6,
48 InteractionOrder -> { NP, 2},
49 Value -> 1,
50 TeX -> Subscript[C,qq],
51 Description -> "coefficient of Oqq "},
52
53
54 C1qu== {
55 ParameterType -> External,
56 BlockName -> DIM6,
57 InteractionOrder -> { NP, 2},
58 Value -> 1,
59 TeX -> Subsuperscript[C,qu,1],
60 Description -> "coefficient of O1qu "},
61
62 C8qu== {
63 ParameterType -> External,
64 BlockName -> DIM6,
65 InteractionOrder -> { NP, 2},
66 Value -> 1,
67 TeX -> Subsuperscript[C,qu,8],
68 Description -> "coefficient of O8qu "},
69
70
71 C1s== {
72 ParameterType -> External,
73 BlockName -> DIM6,
74 InteractionOrder -> { NP, 2},
75 Value -> 1,
76 TeX -> Subsuperscript[C,s,1],
77 Description -> "coefficient of O1ss "},
78
79 C8s== {
80 ParameterType -> External,
81 BlockName -> DIM6,
82 InteractionOrder -> { NP, 2},
83 Value -> 1,
84 TeX -> Subsuperscript[C,s,8],
85 Description -> "coefficient of O8ss "},
86
87 C1ps== {
88 ParameterType -> External,
89 BlockName -> DIM6,
90 InteractionOrder -> { NP, 2},
91 Value -> 1,
92 TeX -> Subsuperscript[C,ps,1],
93 Description -> "coefficient of O1psps "},
94
95 C8ps== {
96 ParameterType -> External,
97 BlockName -> DIM6,
98 InteractionOrder -> { NP, 2},
99 Value -> 1,
100 TeX -> Subsuperscript[C,ps,8],
101 Description -> "coefficient of O8psps "},
102
103
104
105 CuuLam == {
106 ParameterType -> Internal,
107 InteractionOrder -> { NP, 2},
108 Value -> Cuu/LamNP/LamNP,
109 TeX -> Subscript[C,uu]/\[CapitalLambda]^2,
110 Description -> "coefficient of Ouu in TeV-2"},
111
112 CqqLam== {
113 ParameterType -> Internal,
114InteractionOrder -> { NP, 2},
115 Value -> Cqq/LamNP/LamNP,
116 TeX -> Subscript[C,qq]/\[CapitalLambda]^2,
117 Description -> "coefficient of Oqq in TeV-2"},
118
119
120 C1quLam== {
121 ParameterType -> Internal,
122InteractionOrder -> { NP, 2},
123 Value -> C1qu/LamNP/LamNP,
124 TeX -> Subsuperscript[C,qu,1]/\[CapitalLambda]^2,
125 Description -> "coefficient of O1qu in TeV-2"},
126
127 C8quLam== {
128 ParameterType -> Internal,
129InteractionOrder -> { NP, 2},
130 Value -> C8qu/LamNP/LamNP,
131 TeX -> Subsuperscript[C,qu,8]/\[CapitalLambda]^2,
132 Description -> "coefficient of O8qu in TeV-2"},
133
134
135 C1sLam== {
136 ParameterType -> Internal,
137 InteractionOrder -> { NP, 2},
138 Value -> C1s/LamNP/LamNP,
139 TeX -> Subsuperscript[C,s,1]/\[CapitalLambda]^2,
140 Description -> "coefficient of O1ss in TeV-2"},
141
142 C8sLam== {
143 ParameterType -> Internal,
144 InteractionOrder -> { NP, 2},
145 Value -> C8s/LamNP/LamNP,
146 TeX -> Subsuperscript[C,s,8]/\[CapitalLambda]^2,
147 Description -> "coefficient of O8ss in TeV-2"},
148
149 C1psLam== {
150 ParameterType -> Internal,
151 InteractionOrder -> { NP, 2},
152 Value -> C1ps/LamNP/LamNP,
153 TeX -> Subsuperscript[C,ps,1]/\[CapitalLambda]^2,
154 Description -> "coefficient of O1psps in TeV-2"},
155
156 C8psLam== {
157 ParameterType -> Internal,
158 InteractionOrder -> { NP, 2},
159 Value -> C8ps/LamNP/LamNP,
160 TeX -> Subsuperscript[C,ps,8]/\[CapitalLambda]^2,
161 Description -> "coefficient of O8psps in TeV-2"}
162
163
164};
165
166
167 L4topVector := Expand[ExpandIndices[
168 (C8quLam* Module[{ii,jj,kk,kk1,jj1,mu,aa,cc,sp1,sp2,sp3,sp4,cc1,cc2,cc3,cc4}, ( T[aa,cc,cc2]T[aa,cc3,cc4] uqbar[sp1,3,cc].uq[sp2,3,cc2] Ga[mu,sp1,spa] ProjM[spa,sp2]) uqbar[sp3,3,cc3].uq[sp4,3,cc4] Ga[mu,sp3,spb] ProjP[spb,sp4]]+
169 C1quLam* Module[{ii,jj,kk,kk1,jj1,mu,aa,cc,sp1,sp2,sp3,sp4,cc1,cc2,cc3}, (uqbar[sp1,3,cc].uq[sp2,3,cc] Ga[mu,sp1,spa] ProjM[spa,sp2]) uqbar[sp3,3,cc3].uq[sp4,3,cc3] Ga[mu,sp3,spb] ProjP[spb,sp4] ]+
170 CqqLam* Module[{ii,jj,kk,kk1,jj1,mu,aa,cc,sp1,sp2,sp3,sp4,cc1,cc2,cc3}, (uqbar[sp1,3,cc].uq[sp2,3,cc] Ga[mu,sp1,spa] ProjM[spa,sp2]) uqbar[sp3,3,cc3].uq[sp4,3,cc3] Ga[mu,sp3,spb] ProjM[spb,sp4] ]+
171 CuuLam* Module[{ii,jj,kk,kk1,jj1,mu,aa,cc,sp1,sp2,sp3,sp4,cc1,cc2,cc3}, (uqbar[sp1,3,cc].uq[sp2,3,cc] Ga[mu,sp1,spa] ProjP[spa,sp2]) uqbar[sp3,3,cc3].uq[sp4,3,cc3] Ga[mu,sp3,spb] ProjP[spb,sp4] ])/1000^2 ]];
172
173
174 L4topScalar := Expand[ExpandIndices[(C1sLam* Module[{mm,mm2,sp,sp1,sp2,sp3,sp4}, uqbar[sp,3,mm].uq[sp,3,mm]* uqbar[sp3,3,mm2].uq[sp3,3,mm2]] +
175 C8sLam* Module[{aa,mm,mm2,nn,nn2,sp,sp1,sp2,sp3,sp4,ff,gg}, T[aa,mm,nn] T[aa,mm2,nn2]( ProjM[sp,sp2] + ProjP[sp,sp2]) uqbar[sp,3,mm].uq[sp2,3,nn]*( ProjM[sp3,sp4] + ProjP[sp3,sp4]) uqbar[sp3,3,mm2].uq[sp4,3,nn2]] +
176 C1psLam* Module[{mm,mm2,sp,sp1,sp2,sp3,sp4}, ( -ProjM[sp,sp2] + ProjP[sp,sp2]) uqbar[sp,3,mm].uq[sp2,3,mm]*( -ProjM[sp3,sp4] + ProjP[sp3,sp4]) uqbar[sp3,3,mm2].uq[sp4,3,mm2]] +
177 C8psLam* Module[{aa,mm,mm2,nn,nn2,sp,sp1,sp2,sp3,sp4,ff,gg}, T[aa,mm,nn] T[aa,mm2,nn2]( -ProjM[sp,sp2] + ProjP[sp,sp2]) uqbar[sp,3,mm].uq[sp2,3,nn]*( -ProjM[sp3,sp4] + ProjP[sp3,sp4]) uqbar[sp3,3,mm2].uq[sp4,3,nn2]])/1000^2]];
178
179
180LagNP := L4topVector + L4topScalar