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 |
|
---|
9 | M$ModelName = "4topEFT";
|
---|
10 |
|
---|
11 |
|
---|
12 | M$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 |
|
---|
26 | M$Parameters = {
|
---|
27 |
|
---|
28 | (* External parameters *)
|
---|
29 |
|
---|
30 | LamNP == {
|
---|
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,
|
---|
114 | InteractionOrder -> { 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,
|
---|
122 | InteractionOrder -> { 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,
|
---|
129 | InteractionOrder -> { 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 |
|
---|
180 | LagNP := L4topVector + L4topScalar
|
---|