1 | (***************************************************************************************************************)
|
---|
2 | (****** FeynRules mod-file for Model Independent searches of top partners ******)
|
---|
3 | (****** X(5/3), T(2/3), B(-1/3) & Y(-4/3) with arbitrary couplings ******)
|
---|
4 | (****** ******)
|
---|
5 | (****** Authors: M. Buchkremer, G. Cacciapaglia, A. Deandrea,L. Panizzi ******)
|
---|
6 | (****** ******)
|
---|
7 | (***************************************************************************************************************)
|
---|
8 |
|
---|
9 | M$ModelName = "VLQ";
|
---|
10 |
|
---|
11 |
|
---|
12 | M$Information = {Authors -> {"M. Buchkremer","G. Cacciapaglia","A. Deandrea","L. Panizzi"},
|
---|
13 | Version -> "1.2.5",
|
---|
14 | Date -> "10. 04. 2013",
|
---|
15 | Institutions -> {"Universite catholique de Louvain (CP3)","Universite de Lyon (CNRS/IN2P3)","University of Southampton (School of Physics and Astronomy)"},
|
---|
16 | Emails -> {"mathieu.buchkremer@uclouvain.be", "g.cacciapaglia@ipnl.in2p3.fr","deandrea@ipnl.in2p3.fr", "l.panizzi@soton.ac.uk"}};
|
---|
17 |
|
---|
18 |
|
---|
19 | (******* Index definitions ********)
|
---|
20 |
|
---|
21 | IndexRange[ Index[Generation] ] = Range[3]
|
---|
22 | IndexRange[ Index[Colour] ] = NoUnfold[Range[3]]
|
---|
23 | IndexRange[ Index[Gluon] ] = NoUnfold[Range[8]]
|
---|
24 | IndexRange[ Index[SU2W] ] = Unfold[Range[3]]
|
---|
25 | IndexStyle[Colour, i]
|
---|
26 | IndexStyle[Generation, f]
|
---|
27 | IndexStyle[Gluon ,a]
|
---|
28 | IndexStyle[SU2W ,k]
|
---|
29 |
|
---|
30 | (******* Gauge parameters (for FeynArts) ********)
|
---|
31 |
|
---|
32 | GaugeXi[ V[1] ] = GaugeXi[A];
|
---|
33 | GaugeXi[ V[2] ] = GaugeXi[Z];
|
---|
34 | GaugeXi[ V[3] ] = GaugeXi[W];
|
---|
35 | GaugeXi[ V[4] ] = GaugeXi[G];
|
---|
36 | GaugeXi[ S[1] ] = 1;
|
---|
37 | GaugeXi[ S[2] ] = GaugeXi[Z];
|
---|
38 | GaugeXi[ S[3] ] = GaugeXi[W];
|
---|
39 | GaugeXi[ U[1] ] = GaugeXi[A];
|
---|
40 | GaugeXi[ U[2] ] = GaugeXi[Z];
|
---|
41 | GaugeXi[ U[31] ] = GaugeXi[W];
|
---|
42 | GaugeXi[ U[32] ] = GaugeXi[W];
|
---|
43 | GaugeXi[ U[4] ] = GaugeXi[G];
|
---|
44 |
|
---|
45 | (**************** Parameters *************)
|
---|
46 |
|
---|
47 | M$Parameters = {
|
---|
48 |
|
---|
49 | (* External parameters, SM *)
|
---|
50 |
|
---|
51 | \[Alpha]EWM1== {
|
---|
52 | ParameterType -> External,
|
---|
53 | BlockName -> SMINPUTS,
|
---|
54 | ParameterName -> aEWM1,
|
---|
55 | InteractionOrder -> {QED, -2},
|
---|
56 | Value -> 127.9,
|
---|
57 | Description -> "Inverse of the electroweak coupling constant"},
|
---|
58 |
|
---|
59 | Gf == {
|
---|
60 | ParameterType -> External,
|
---|
61 | BlockName -> SMINPUTS,
|
---|
62 | TeX -> Subscript[G, f],
|
---|
63 | InteractionOrder -> {QED, 2},
|
---|
64 | Value -> 1.16600 * 10^(-5),
|
---|
65 | Description -> "Fermi constant"},
|
---|
66 |
|
---|
67 | \[Alpha]S == {
|
---|
68 | ParameterType -> External,
|
---|
69 | BlockName -> SMINPUTS,
|
---|
70 | TeX -> Subscript[\[Alpha], s],
|
---|
71 | ParameterName -> aS,
|
---|
72 | InteractionOrder -> {QCD, 2},
|
---|
73 | Value -> 0.118,
|
---|
74 | Description -> "Strong coupling constant at the Z pole."},
|
---|
75 |
|
---|
76 | ymdo == {
|
---|
77 | ParameterType -> External,
|
---|
78 | BlockName -> YUKAWA,
|
---|
79 | Value -> 5.04*10^(-3),
|
---|
80 | OrderBlock -> {1},
|
---|
81 | Description -> "Down Yukawa mass"},
|
---|
82 |
|
---|
83 | ymup == {
|
---|
84 | ParameterType -> External,
|
---|
85 | BlockName -> YUKAWA,
|
---|
86 | Value -> 2.55*10^(-3),
|
---|
87 | OrderBlock -> {2},
|
---|
88 | Description -> "Up Yukawa mass"},
|
---|
89 |
|
---|
90 | yms == {
|
---|
91 | ParameterType -> External,
|
---|
92 | BlockName -> YUKAWA,
|
---|
93 | Value -> 0.101,
|
---|
94 | OrderBlock -> {3},
|
---|
95 | Description -> "Strange Yukawa mass"},
|
---|
96 |
|
---|
97 | ymc == {
|
---|
98 | ParameterType -> External,
|
---|
99 | BlockName -> YUKAWA,
|
---|
100 | Value -> 1.25,
|
---|
101 | OrderBlock -> {4},
|
---|
102 | Description -> "Charm Yukawa mass"},
|
---|
103 |
|
---|
104 | ymb == {
|
---|
105 | ParameterType -> External,
|
---|
106 | BlockName -> YUKAWA,
|
---|
107 | Value -> 4.2,
|
---|
108 | OrderBlock -> {5},
|
---|
109 | Description -> "Bottom Yukawa mass"},
|
---|
110 |
|
---|
111 | ymt == {
|
---|
112 | ParameterType -> External,
|
---|
113 | BlockName -> YUKAWA,
|
---|
114 | Value -> 174.3,
|
---|
115 | OrderBlock -> {6},
|
---|
116 | Description -> "Top Yukawa mass"},
|
---|
117 |
|
---|
118 | yme == {
|
---|
119 | ParameterType -> External,
|
---|
120 | BlockName -> YUKAWA,
|
---|
121 | Value -> 5.11*10^(-4),
|
---|
122 | OrderBlock -> {11},
|
---|
123 | Description -> "Electron Yukawa mass"},
|
---|
124 |
|
---|
125 | ymm == {
|
---|
126 | ParameterType -> External,
|
---|
127 | BlockName -> YUKAWA,
|
---|
128 | Value -> 0.10566,
|
---|
129 | OrderBlock -> {13},
|
---|
130 | Description -> "Muon Yukawa mass"},
|
---|
131 |
|
---|
132 | ymtau == {
|
---|
133 | ParameterType -> External,
|
---|
134 | BlockName -> YUKAWA,
|
---|
135 | Value -> 1.777,
|
---|
136 | OrderBlock -> {15},
|
---|
137 | Description -> "Tau Yukawa mass"},
|
---|
138 |
|
---|
139 | CKM == {
|
---|
140 | ParameterType -> External,
|
---|
141 | BlockName -> CKMBlock,
|
---|
142 | ComplexParameter -> False,
|
---|
143 | Indices -> {Index[Generation], Index[Generation]},
|
---|
144 | TensorClass -> CKM,
|
---|
145 | Unitary -> True,
|
---|
146 | Value -> {CKM[1,1] -> 0.97428,
|
---|
147 | CKM[1,2] -> 0.2253,
|
---|
148 | CKM[1,3] -> 0.00347,
|
---|
149 | CKM[2,1] -> 0.2252,
|
---|
150 | CKM[2,2] -> 0.97345,
|
---|
151 | CKM[2,3] -> 0.0410,
|
---|
152 | CKM[3,1] -> 0.00862,
|
---|
153 | CKM[3,2] -> 0.0403,
|
---|
154 | CKM[3,3] -> 0.999152},
|
---|
155 | Description -> "SM CKM Matrix"},
|
---|
156 |
|
---|
157 | (* External parameters, VLQ *)
|
---|
158 |
|
---|
159 | KX == {
|
---|
160 | ParameterType -> External,
|
---|
161 | BlockName -> Kappa,
|
---|
162 | ComplexParameter -> False,
|
---|
163 | Value -> 1,
|
---|
164 | Description -> "Kappa_X parameter"},
|
---|
165 |
|
---|
166 | KT == {
|
---|
167 | ParameterType -> External,
|
---|
168 | BlockName -> Kappa,
|
---|
169 | ComplexParameter -> False,
|
---|
170 | Value -> 1,
|
---|
171 | Description -> "Kappa_T parameter"},
|
---|
172 |
|
---|
173 | KB == {
|
---|
174 | ParameterType -> External,
|
---|
175 | BlockName -> Kappa,
|
---|
176 | ComplexParameter -> False,
|
---|
177 | Value -> 1,
|
---|
178 | Description -> "Kappa_B parameter"},
|
---|
179 |
|
---|
180 | KY == {
|
---|
181 | ParameterType -> External,
|
---|
182 | BlockName -> Kappa,
|
---|
183 | ComplexParameter -> False,
|
---|
184 | Value -> 1,
|
---|
185 | Description -> "Kappa_Y parameter"},
|
---|
186 |
|
---|
187 | xitpw == {
|
---|
188 | ParameterType -> External,
|
---|
189 | BlockName -> Xi,
|
---|
190 | ComplexParameter -> False,
|
---|
191 | Value -> 0.4,
|
---|
192 | Description -> "Branching ratio of T in W"},
|
---|
193 |
|
---|
194 | xitpz == {
|
---|
195 | ParameterType -> External,
|
---|
196 | BlockName -> Xi,
|
---|
197 | ComplexParameter -> False,
|
---|
198 | Value -> 0.3,
|
---|
199 | Description -> "Branching ratio of T in Z"},
|
---|
200 |
|
---|
201 | xitph == {
|
---|
202 | ParameterType -> External,
|
---|
203 | BlockName -> Xi,
|
---|
204 | ComplexParameter -> False,
|
---|
205 | Value -> 0.3,
|
---|
206 | Description -> "Branching ratio of T in H"},
|
---|
207 |
|
---|
208 | xibpw == {
|
---|
209 | ParameterType -> External,
|
---|
210 | BlockName -> Xi,
|
---|
211 | ComplexParameter -> False,
|
---|
212 | Value -> 0.4,
|
---|
213 | Description -> "Branching ratio of B in W"},
|
---|
214 |
|
---|
215 | xibpz == {
|
---|
216 | ParameterType -> External,
|
---|
217 | BlockName -> Xi,
|
---|
218 | ComplexParameter -> False,
|
---|
219 | Value -> 0.3,
|
---|
220 | Description -> "Branching ratio of B in Z"},
|
---|
221 |
|
---|
222 | xibph == {
|
---|
223 | ParameterType -> External,
|
---|
224 | BlockName -> Xi,
|
---|
225 | ComplexParameter -> False,
|
---|
226 | Value -> 0.3,
|
---|
227 | Description -> "Branching ratio of B in H"},
|
---|
228 |
|
---|
229 | zetaXuL == {
|
---|
230 | ParameterType -> External,
|
---|
231 | BlockName -> Zeta,
|
---|
232 | ComplexParameter -> False,
|
---|
233 | Value -> 0.3,
|
---|
234 | Description -> "X-u mixing (left-handed)"},
|
---|
235 |
|
---|
236 | zetaXcL == {
|
---|
237 | ParameterType -> External,
|
---|
238 | BlockName -> Zeta,
|
---|
239 | ComplexParameter -> False,
|
---|
240 | Value -> 0.3,
|
---|
241 | Description -> "X-c mixing (left-handed)"},
|
---|
242 |
|
---|
243 | zetaXtL == {
|
---|
244 | ParameterType -> External,
|
---|
245 | BlockName -> Zeta,
|
---|
246 | ComplexParameter -> False,
|
---|
247 | Value -> 0.4,
|
---|
248 | Description -> "X-t mixing (left-handed)"},
|
---|
249 |
|
---|
250 | zetaTuL == {
|
---|
251 | ParameterType -> External,
|
---|
252 | BlockName -> Zeta,
|
---|
253 | ComplexParameter -> False,
|
---|
254 | Value -> 0.3,
|
---|
255 | Description -> "T-u mixing (left-handed)"},
|
---|
256 |
|
---|
257 | zetaTcL == {
|
---|
258 | ParameterType -> External,
|
---|
259 | BlockName -> Zeta,
|
---|
260 | ComplexParameter -> False,
|
---|
261 | Value -> 0.3,
|
---|
262 | Description -> "T-c mixing (left-handed)"},
|
---|
263 |
|
---|
264 | zetaTtL == {
|
---|
265 | ParameterType -> External,
|
---|
266 | BlockName -> Zeta,
|
---|
267 | ComplexParameter -> False,
|
---|
268 | Value -> 0.4,
|
---|
269 | Description -> "T-t mixing (left-handed)"},
|
---|
270 |
|
---|
271 | zetaBdL == {
|
---|
272 | ParameterType -> External,
|
---|
273 | BlockName -> Zeta,
|
---|
274 | ComplexParameter -> False,
|
---|
275 | Value -> 0.3,
|
---|
276 | Description -> "B-d mixing (left-handed)"},
|
---|
277 |
|
---|
278 | zetaBsL == {
|
---|
279 | ParameterType -> External,
|
---|
280 | BlockName -> Zeta,
|
---|
281 | ComplexParameter -> False,
|
---|
282 | Value -> 0.3,
|
---|
283 | Description -> "B-s mixing (left-handed)"},
|
---|
284 |
|
---|
285 | zetaBbL == {
|
---|
286 | ParameterType -> External,
|
---|
287 | BlockName -> Zeta,
|
---|
288 | ComplexParameter -> False,
|
---|
289 | Value -> 0.4,
|
---|
290 | Description -> "B-b mixing (left-handed)"},
|
---|
291 |
|
---|
292 | zetaYdL == {
|
---|
293 | ParameterType -> External,
|
---|
294 | BlockName -> Zeta,
|
---|
295 | ComplexParameter -> False,
|
---|
296 | Value -> 0.3,
|
---|
297 | Description -> "Y-d mixing (left-handed)"},
|
---|
298 |
|
---|
299 | zetaYsL == {
|
---|
300 | ParameterType -> External,
|
---|
301 | BlockName -> Zeta,
|
---|
302 | ComplexParameter -> False,
|
---|
303 | Value -> 0.3,
|
---|
304 | Description -> "Y-s mixing (left-handed)"},
|
---|
305 |
|
---|
306 | zetaYbL == {
|
---|
307 | ParameterType -> External,
|
---|
308 | BlockName -> Zeta,
|
---|
309 | ComplexParameter -> False,
|
---|
310 | Value -> 0.4,
|
---|
311 | Description -> "Y-b mixing (left-handed)"},
|
---|
312 |
|
---|
313 |
|
---|
314 | zetaXuR == {
|
---|
315 | ParameterType -> External,
|
---|
316 | BlockName -> Zeta,
|
---|
317 | ComplexParameter -> False,
|
---|
318 | Value -> 0,
|
---|
319 | Description -> "X-u mixing (right-handed)"},
|
---|
320 |
|
---|
321 | zetaXcR == {
|
---|
322 | ParameterType -> External,
|
---|
323 | BlockName -> Zeta,
|
---|
324 | ComplexParameter -> False,
|
---|
325 | Value -> 0,
|
---|
326 | Description -> "X-c mixing (right-handed)"},
|
---|
327 |
|
---|
328 | zetaXtR == {
|
---|
329 | ParameterType -> External,
|
---|
330 | BlockName -> Zeta,
|
---|
331 | ComplexParameter -> False,
|
---|
332 | Value -> 0,
|
---|
333 | Description -> "X-t mixing (right-handed)"},
|
---|
334 |
|
---|
335 | zetaTuR == {
|
---|
336 | ParameterType -> External,
|
---|
337 | BlockName -> Zeta,
|
---|
338 | ComplexParameter -> False,
|
---|
339 | Value -> 0,
|
---|
340 | Description -> "T-u mixing (right-handed)"},
|
---|
341 |
|
---|
342 | zetaTcR == {
|
---|
343 | ParameterType -> External,
|
---|
344 | BlockName -> Zeta,
|
---|
345 | ComplexParameter -> False,
|
---|
346 | Value -> 0,
|
---|
347 | Description -> "T-c mixing (right-handed)"},
|
---|
348 |
|
---|
349 | zetaTtR == {
|
---|
350 | ParameterType -> External,
|
---|
351 | BlockName -> Zeta,
|
---|
352 | ComplexParameter -> False,
|
---|
353 | Value -> 0,
|
---|
354 | Description -> "T-t mixing (right-handed)"},
|
---|
355 |
|
---|
356 | zetaBdR == {
|
---|
357 | ParameterType -> External,
|
---|
358 | BlockName -> Zeta,
|
---|
359 | ComplexParameter -> False,
|
---|
360 | Value -> 0,
|
---|
361 | Description -> "B-d mixing (right-handed)"},
|
---|
362 |
|
---|
363 | zetaBsR == {
|
---|
364 | ParameterType -> External,
|
---|
365 | BlockName -> Zeta,
|
---|
366 | ComplexParameter -> False,
|
---|
367 | Value -> 0,
|
---|
368 | Description -> "B-s mixing (right-handed)"},
|
---|
369 |
|
---|
370 | zetaBbR == {
|
---|
371 | ParameterType -> External,
|
---|
372 | BlockName -> Zeta,
|
---|
373 | ComplexParameter -> False,
|
---|
374 | Value -> 0,
|
---|
375 | Description -> "B-b mixing (right-handed)"},
|
---|
376 |
|
---|
377 | zetaYdR == {
|
---|
378 | ParameterType -> External,
|
---|
379 | BlockName -> Zeta,
|
---|
380 | ComplexParameter -> False,
|
---|
381 | Value -> 0,
|
---|
382 | Description -> "Y-d mixing (right-handed)"},
|
---|
383 |
|
---|
384 | zetaYsR == {
|
---|
385 | ParameterType -> External,
|
---|
386 | BlockName -> Zeta,
|
---|
387 | ComplexParameter -> False,
|
---|
388 | Value -> 0,
|
---|
389 | Description -> "Y-s mixing (right-handed)"},
|
---|
390 |
|
---|
391 | zetaYbR == {
|
---|
392 | ParameterType -> External,
|
---|
393 | BlockName -> Zeta,
|
---|
394 | ComplexParameter -> False,
|
---|
395 | Value -> 0,
|
---|
396 | Description -> "Y-b mixing (right-handed)"},
|
---|
397 |
|
---|
398 |
|
---|
399 | (* Internal Parameters, SM *)
|
---|
400 |
|
---|
401 | \[Alpha]EW == {
|
---|
402 | ParameterType -> Internal,
|
---|
403 | Value -> 1/\[Alpha]EWM1,
|
---|
404 | TeX -> Subscript[\[Alpha], EW],
|
---|
405 | ParameterName -> aEW,
|
---|
406 | InteractionOrder -> {QED, 2},
|
---|
407 | Description -> "Electroweak coupling contant"},
|
---|
408 |
|
---|
409 |
|
---|
410 | MW == {
|
---|
411 | ParameterType -> Internal,
|
---|
412 | Value -> Sqrt[MZ^2/2+Sqrt[MZ^4/4-Pi/Sqrt[2]*\[Alpha]EW/Gf*MZ^2]],
|
---|
413 | TeX -> Subscript[M, W],
|
---|
414 | Description -> "W mass"},
|
---|
415 |
|
---|
416 | sw2 == {
|
---|
417 | ParameterType -> Internal,
|
---|
418 | Value -> 1-(MW/MZ)^2,
|
---|
419 | Description -> "Squared Sin of the Weinberg angle"},
|
---|
420 |
|
---|
421 | ee == {
|
---|
422 | TeX -> e,
|
---|
423 | ParameterType -> Internal,
|
---|
424 | Value -> Sqrt[4 Pi \[Alpha]EW],
|
---|
425 | InteractionOrder -> {QED, 1},
|
---|
426 | Description -> "Electric coupling constant"},
|
---|
427 |
|
---|
428 | cw == {
|
---|
429 | TeX -> Subscript[c, w],
|
---|
430 | ParameterType -> Internal,
|
---|
431 | Value -> Sqrt[1 - sw2],
|
---|
432 | Description -> "Cos of the Weinberg angle"},
|
---|
433 |
|
---|
434 | sw == {
|
---|
435 | TeX -> Subscript[s, w],
|
---|
436 | ParameterType -> Internal,
|
---|
437 | Value -> Sqrt[sw2],
|
---|
438 | Description -> "Sin of the Weinberg angle"},
|
---|
439 |
|
---|
440 | gw == {
|
---|
441 | TeX -> Subscript[g, w],
|
---|
442 | ParameterType -> Internal,
|
---|
443 | Value -> ee / sw,
|
---|
444 | InteractionOrder -> {QED, 1},
|
---|
445 | Description -> "Weak coupling constant"},
|
---|
446 |
|
---|
447 | g1 == {
|
---|
448 | TeX -> Subscript[g, 1],
|
---|
449 | ParameterType -> Internal,
|
---|
450 | Value -> ee / cw,
|
---|
451 | InteractionOrder -> {QED, 1},
|
---|
452 | Description -> "U(1)Y coupling constant"},
|
---|
453 |
|
---|
454 | gs == {
|
---|
455 | TeX -> Subscript[g, s],
|
---|
456 | ParameterType -> Internal,
|
---|
457 | Value -> Sqrt[4 Pi \[Alpha]S],
|
---|
458 | InteractionOrder -> {QCD, 1},
|
---|
459 | ParameterName -> G,
|
---|
460 | Description -> "Strong coupling constant"},
|
---|
461 |
|
---|
462 | v == {
|
---|
463 | ParameterType -> Internal,
|
---|
464 | Value -> 2*MW*sw/ee,
|
---|
465 | InteractionOrder -> {QED, -1},
|
---|
466 | Description -> "Higgs VEV"},
|
---|
467 |
|
---|
468 | \[Lambda] == {
|
---|
469 | ParameterType -> Internal,
|
---|
470 | Value -> MH^2/(2*v^2),
|
---|
471 | InteractionOrder -> {QED, 2},
|
---|
472 | ParameterName -> lam,
|
---|
473 | Description -> "Higgs quartic coupling"},
|
---|
474 |
|
---|
475 | muH == {
|
---|
476 | ParameterType -> Internal,
|
---|
477 | Value -> Sqrt[v^2 \[Lambda]],
|
---|
478 | TeX -> \[Mu],
|
---|
479 | Description -> "Coefficient of the quadratic piece of the Higgs potential"},
|
---|
480 |
|
---|
481 | yl == {
|
---|
482 | TeX -> Superscript[y, l],
|
---|
483 | Indices -> {Index[Generation]},
|
---|
484 | AllowSummation -> True,
|
---|
485 | ParameterType -> Internal,
|
---|
486 | Value -> {yl[1] -> Sqrt[2] yme / v, yl[2] -> Sqrt[2] ymm / v, yl[3] -> Sqrt[2] ymtau / v},
|
---|
487 | ParameterName -> {yl[1] -> ye, yl[2] -> ym, yl[3] -> ytau},
|
---|
488 | InteractionOrder -> {QED, 1},
|
---|
489 | ComplexParameter -> False,
|
---|
490 | Description -> "Lepton Yukawa coupling"},
|
---|
491 |
|
---|
492 | yu == {
|
---|
493 | TeX -> Superscript[y, u],
|
---|
494 | Indices -> {Index[Generation]},
|
---|
495 | AllowSummation -> True,
|
---|
496 | ParameterType -> Internal,
|
---|
497 | Value -> {yu[1] -> Sqrt[2] ymup / v, yu[2] -> Sqrt[2] ymc / v, yu[3] -> Sqrt[2] ymt / v},
|
---|
498 | ParameterName -> {yu[1] -> yup, yu[2] -> yc, yu[3] -> yt},
|
---|
499 | InteractionOrder -> {QED, 1},
|
---|
500 | ComplexParameter -> False,
|
---|
501 | Description -> "U-quark Yukawa coupling"},
|
---|
502 |
|
---|
503 | yd == {
|
---|
504 | TeX -> Superscript[y, d],
|
---|
505 | Indices -> {Index[Generation]},
|
---|
506 | AllowSummation -> True,
|
---|
507 | ParameterType -> Internal,
|
---|
508 | Value -> {yd[1] -> Sqrt[2] ymdo / v, yd[2] -> Sqrt[2] yms / v, yd[3] -> Sqrt[2] ymb / v},
|
---|
509 | ParameterName -> {yd[1] -> ydo, yd[2] -> ys, yd[3] -> yb},
|
---|
510 | InteractionOrder -> {QED, 1},
|
---|
511 | ComplexParameter -> False,
|
---|
512 | Description -> "D-quark Yukawa coupling"},
|
---|
513 |
|
---|
514 |
|
---|
515 | (************** Internal Parameters, VLQ **************)
|
---|
516 | (* X couplings *)
|
---|
517 |
|
---|
518 | KXuL == {
|
---|
519 | ParameterType -> Internal,
|
---|
520 | BlockName -> WIDTH,
|
---|
521 | ComplexParameter -> False,
|
---|
522 | Value -> (ee/sw*Sqrt[zetaXuL/gamma0xw])/Sqrt[2],
|
---|
523 | InteractionOrder -> {QED, 1},
|
---|
524 | Description -> "XuW coupling (left-handed)"},
|
---|
525 |
|
---|
526 | KXcL == {
|
---|
527 | ParameterType -> Internal,
|
---|
528 | BlockName -> WIDTH,
|
---|
529 | ComplexParameter -> False,
|
---|
530 | Value -> (ee/sw*Sqrt[zetaXcL/gamma0xw])/Sqrt[2],
|
---|
531 | InteractionOrder -> {QED, 1},
|
---|
532 | Description -> "XcW coupling (left-handed)"},
|
---|
533 |
|
---|
534 | KXtL == {
|
---|
535 | ParameterType -> Internal,
|
---|
536 | BlockName -> WIDTH,
|
---|
537 | ComplexParameter -> False,
|
---|
538 | Value -> (ee/sw*Sqrt[zetaXtL/gamma0xw])/Sqrt[2],
|
---|
539 | InteractionOrder -> {QED, 1},
|
---|
540 | Description -> "XtW coupling (left-handed)"},
|
---|
541 |
|
---|
542 | KXuR == {
|
---|
543 | ParameterType -> Internal,
|
---|
544 | BlockName -> WIDTH,
|
---|
545 | ComplexParameter -> False,
|
---|
546 | Value -> (ee/sw*Sqrt[zetaXuR/gamma0xw])/Sqrt[2],
|
---|
547 | InteractionOrder -> {QED, 1},
|
---|
548 | Description -> "XuW coupling (right-handed)"},
|
---|
549 |
|
---|
550 | KXcR == {
|
---|
551 | ParameterType -> Internal,
|
---|
552 | BlockName -> WIDTH,
|
---|
553 | ComplexParameter -> False,
|
---|
554 | Value -> (ee/sw*Sqrt[zetaXcR/gamma0xw])/Sqrt[2],
|
---|
555 | InteractionOrder -> {QED, 1},
|
---|
556 | Description -> "XcW coupling (right-handed)"},
|
---|
557 |
|
---|
558 | KXtR == {
|
---|
559 | ParameterType -> Internal,
|
---|
560 | BlockName -> Kappa,
|
---|
561 | ComplexParameter -> False,
|
---|
562 | Value -> (ee/sw*Sqrt[zetaXtR/gamma0xw])/Sqrt[2],
|
---|
563 | InteractionOrder -> {QED, 1},
|
---|
564 | Description -> "XtW coupling (right-handed)"},
|
---|
565 |
|
---|
566 | (* Y couplings *)
|
---|
567 |
|
---|
568 | KYdL == {
|
---|
569 | ParameterType -> Internal,
|
---|
570 | BlockName -> Kappa,
|
---|
571 | ComplexParameter -> False,
|
---|
572 | Value -> (ee/sw*Sqrt[zetaYdL/gamma0yw])/Sqrt[2],
|
---|
573 | InteractionOrder -> {QED, 1},
|
---|
574 | Description -> "YdW coupling (left-handed)"},
|
---|
575 |
|
---|
576 | KYsL == {
|
---|
577 | ParameterType -> Internal,
|
---|
578 | BlockName -> Kappa,
|
---|
579 | ComplexParameter -> False,
|
---|
580 | Value -> (ee/sw*Sqrt[zetaYsL/gamma0yw])/Sqrt[2],
|
---|
581 | InteractionOrder -> {QED, 1},
|
---|
582 | Description -> "YsW coupling (left-handed)"},
|
---|
583 |
|
---|
584 | KYbL == {
|
---|
585 | ParameterType -> Internal,
|
---|
586 | BlockName -> Kappa,
|
---|
587 | ComplexParameter -> False,
|
---|
588 | Value -> (ee/sw*Sqrt[zetaYbL/gamma0yw])/Sqrt[2],
|
---|
589 | InteractionOrder -> {QED, 1},
|
---|
590 | Description -> "YbW coupling (left-handed)"},
|
---|
591 |
|
---|
592 | KYdR == {
|
---|
593 | ParameterType -> Internal,
|
---|
594 | BlockName -> Kappa,
|
---|
595 | ComplexParameter -> False,
|
---|
596 | Value -> (ee/sw*Sqrt[zetaYdR/gamma0yw])/Sqrt[2],
|
---|
597 | InteractionOrder -> {QED, 1},
|
---|
598 | Description -> "YdW coupling (right-handed)"},
|
---|
599 |
|
---|
600 | KYsR == {
|
---|
601 | ParameterType -> Internal,
|
---|
602 | BlockName -> Kappa,
|
---|
603 | ComplexParameter -> False,
|
---|
604 | Value -> (ee/sw*Sqrt[zetaYsR/gamma0yw])/Sqrt[2],
|
---|
605 | InteractionOrder -> {QED, 1},
|
---|
606 | Description -> "YsW coupling (right-handed)"},
|
---|
607 |
|
---|
608 | KYbR == {
|
---|
609 | ParameterType -> Internal,
|
---|
610 | BlockName -> Kappa,
|
---|
611 | ComplexParameter -> False,
|
---|
612 | Value -> (ee/sw*Sqrt[zetaYbR/gamma0yw])/Sqrt[2],
|
---|
613 | InteractionOrder -> {QED, 1},
|
---|
614 | Description -> "YbW coupling (right-handed)"},
|
---|
615 |
|
---|
616 | (* T couplings *)
|
---|
617 |
|
---|
618 | KTuLw == {
|
---|
619 | ParameterType -> Internal,
|
---|
620 | BlockName -> Kappa,
|
---|
621 | ComplexParameter -> False,
|
---|
622 | Value -> (ee/sw*Sqrt[zetaTuL*xitpw/gamma0tpw])/Sqrt[2],
|
---|
623 | InteractionOrder -> {QED, 1},
|
---|
624 | Description -> "TuW coupling (left-handed)"},
|
---|
625 |
|
---|
626 | KTcLw == {
|
---|
627 | ParameterType -> Internal,
|
---|
628 | BlockName -> Kappa,
|
---|
629 | ComplexParameter -> False,
|
---|
630 | Value -> (ee/sw*Sqrt[zetaTcL*xitpw/gamma0tpw])/Sqrt[2],
|
---|
631 | InteractionOrder -> {QED, 1},
|
---|
632 | Description -> "TcW coupling (left-handed)"},
|
---|
633 |
|
---|
634 | KTtLw == {
|
---|
635 | ParameterType -> Internal,
|
---|
636 | BlockName -> Kappa,
|
---|
637 | ComplexParameter -> False,
|
---|
638 | Value -> (ee/sw*Sqrt[zetaTtL*xitpw/gamma0tpw])/Sqrt[2],
|
---|
639 | InteractionOrder -> {QED, 1},
|
---|
640 | Description -> "TtW coupling (left-handed)"},
|
---|
641 |
|
---|
642 | KTuRw == {
|
---|
643 | ParameterType -> Internal,
|
---|
644 | BlockName -> Kappa,
|
---|
645 | ComplexParameter -> False,
|
---|
646 | Value -> (ee/sw*Sqrt[zetaTuR*xitpw/gamma0tpw])/Sqrt[2],
|
---|
647 | InteractionOrder -> {QED, 1},
|
---|
648 | Description -> "TuW coupling (right-handed)"},
|
---|
649 |
|
---|
650 | KTcRw == {
|
---|
651 | ParameterType -> Internal,
|
---|
652 | BlockName -> Kappa,
|
---|
653 | ComplexParameter -> False,
|
---|
654 | Value -> (ee/sw*Sqrt[zetaTcR*xitpw/gamma0tpw])/Sqrt[2],
|
---|
655 | InteractionOrder -> {QED, 1},
|
---|
656 | Description -> "TcW coupling (right-handed)"},
|
---|
657 |
|
---|
658 | KTtRw == {
|
---|
659 | ParameterType -> Internal,
|
---|
660 | BlockName -> Kappa,
|
---|
661 | ComplexParameter -> False,
|
---|
662 | Value -> (ee/sw*Sqrt[zetaTtR*xitpw/gamma0tpw])/Sqrt[2],
|
---|
663 | InteractionOrder -> {QED, 1},
|
---|
664 | Description -> "TtW coupling (right-handed)"},
|
---|
665 |
|
---|
666 | KTuLz == {
|
---|
667 | ParameterType -> Internal,
|
---|
668 | BlockName -> Kappa,
|
---|
669 | ComplexParameter -> False,
|
---|
670 | Value -> (ee/sw*Sqrt[zetaTuL*xitpz/gamma0tpz])/2/cw,
|
---|
671 | InteractionOrder -> {QED, 1},
|
---|
672 | Description -> "TuZ coupling (left-handed)"},
|
---|
673 |
|
---|
674 | KTcLz == {
|
---|
675 | ParameterType -> Internal,
|
---|
676 | BlockName -> Kappa,
|
---|
677 | ComplexParameter -> False,
|
---|
678 | Value -> (ee/sw*Sqrt[zetaTcL*xitpz/gamma0tpz])/2/cw,
|
---|
679 | InteractionOrder -> {QED, 1},
|
---|
680 | Description -> "TcZ coupling (left-handed)"},
|
---|
681 |
|
---|
682 | KTtLz == {
|
---|
683 | ParameterType -> Internal,
|
---|
684 | BlockName -> Kappa,
|
---|
685 | ComplexParameter -> False,
|
---|
686 | Value -> (ee/sw*Sqrt[zetaTtL*xitpz/gamma0tpz])/2/cw,
|
---|
687 | InteractionOrder -> {QED, 1},
|
---|
688 | Description -> "TtZ coupling (left-handed)"},
|
---|
689 |
|
---|
690 | KTuRz == {
|
---|
691 | ParameterType -> Internal,
|
---|
692 | BlockName -> Kappa,
|
---|
693 | ComplexParameter -> False,
|
---|
694 | Value -> (ee/sw*Sqrt[zetaTuR*xitpz/gamma0tpz])/2/cw,
|
---|
695 | InteractionOrder -> {QED, 1},
|
---|
696 | Description -> "TuZ coupling (right-handed)"},
|
---|
697 |
|
---|
698 | KTcRz == {
|
---|
699 | ParameterType -> Internal,
|
---|
700 | BlockName -> Kappa,
|
---|
701 | ComplexParameter -> False,
|
---|
702 | Value -> (ee/sw*Sqrt[zetaTcR*xitpz/gamma0tpz])/2/cw,
|
---|
703 | InteractionOrder -> {QED, 1},
|
---|
704 | Description -> "TcZ coupling (right-handed)"},
|
---|
705 |
|
---|
706 | KTtRz == {
|
---|
707 | ParameterType -> Internal,
|
---|
708 | BlockName -> Kappa,
|
---|
709 | ComplexParameter -> False,
|
---|
710 | Value -> (ee/sw*Sqrt[zetaTtR*xitpz/gamma0tpz])/2/cw,
|
---|
711 | InteractionOrder -> {QED, 1},
|
---|
712 | Description -> "TtZ coupling (right-handed)"},
|
---|
713 |
|
---|
714 | KTuLh == {
|
---|
715 | ParameterType -> Internal,
|
---|
716 | BlockName -> Kappa,
|
---|
717 | ComplexParameter -> False,
|
---|
718 | Value -> (Sqrt[zetaTuL*xitph/gamma0tph]),
|
---|
719 | InteractionOrder -> {QED, 0},
|
---|
720 | Description -> "TuH coupling (left-handed)"},
|
---|
721 |
|
---|
722 | KTcLh == {
|
---|
723 | ParameterType -> Internal,
|
---|
724 | BlockName -> Kappa,
|
---|
725 | ComplexParameter -> False,
|
---|
726 | Value -> (Sqrt[zetaTcL*xitph/gamma0tph]),
|
---|
727 | InteractionOrder -> {QED, 0},
|
---|
728 | Description -> "TcH coupling (left-handed)"},
|
---|
729 |
|
---|
730 | KTtLh == {
|
---|
731 | ParameterType -> Internal,
|
---|
732 | BlockName -> Kappa,
|
---|
733 | ComplexParameter -> False,
|
---|
734 | Value -> (Sqrt[zetaTtL*xitph/gamma0tph]),
|
---|
735 | InteractionOrder -> {QED, 0},
|
---|
736 | Description -> "TtH coupling (left-handed)"},
|
---|
737 |
|
---|
738 | KTuRh == {
|
---|
739 | ParameterType -> Internal,
|
---|
740 | BlockName -> Kappa,
|
---|
741 | ComplexParameter -> False,
|
---|
742 | Value -> (Sqrt[zetaTuR*xitph/gamma0tph]),
|
---|
743 | InteractionOrder -> {QED, 0},
|
---|
744 | Description -> "TuH coupling (right-handed)"},
|
---|
745 |
|
---|
746 | KTcRh == {
|
---|
747 | ParameterType -> Internal,
|
---|
748 | BlockName -> Kappa,
|
---|
749 | ComplexParameter -> False,
|
---|
750 | Value -> (Sqrt[zetaTcR*xitph/gamma0tph]),
|
---|
751 | InteractionOrder -> {QED, 0},
|
---|
752 | Description -> "TcH coupling (right-handed)"},
|
---|
753 |
|
---|
754 | KTtRh == {
|
---|
755 | ParameterType -> Internal,
|
---|
756 | BlockName -> Kappa,
|
---|
757 | ComplexParameter -> False,
|
---|
758 | Value -> (Sqrt[zetaTtR*xitph/gamma0tph]),
|
---|
759 | InteractionOrder -> {QED, 0},
|
---|
760 | Description -> "TtH coupling (right-handed)"},
|
---|
761 |
|
---|
762 | (* B couplings *)
|
---|
763 |
|
---|
764 | KBdLw == {
|
---|
765 | ParameterType -> Internal,
|
---|
766 | BlockName -> Kappa,
|
---|
767 | ComplexParameter -> False,
|
---|
768 | Value -> (ee/sw*Sqrt[zetaBdL*xibpw/gamma0bpw])/Sqrt[2],
|
---|
769 | InteractionOrder -> {QED, 1},
|
---|
770 | Description -> "BdW coupling (left-handed)"},
|
---|
771 |
|
---|
772 | KBsLw == {
|
---|
773 | ParameterType -> Internal,
|
---|
774 | BlockName -> Kappa,
|
---|
775 | ComplexParameter -> False,
|
---|
776 | Value -> (ee/sw*Sqrt[zetaBsL*xibpw/gamma0bpw])/Sqrt[2],
|
---|
777 | InteractionOrder -> {QED, 1},
|
---|
778 | Description -> "BsW coupling (left-handed)"},
|
---|
779 |
|
---|
780 | KBbLw == {
|
---|
781 | ParameterType -> Internal,
|
---|
782 | BlockName -> Kappa,
|
---|
783 | ComplexParameter -> False,
|
---|
784 | Value -> (gw*Sqrt[zetaBbL*xibpw/gamma0bpw])/Sqrt[2],
|
---|
785 | InteractionOrder -> {QED, 1},
|
---|
786 | Description -> "BbW coupling (left-handed)"},
|
---|
787 |
|
---|
788 | KBdRw == {
|
---|
789 | ParameterType -> Internal,
|
---|
790 | BlockName -> Kappa,
|
---|
791 | ComplexParameter -> False,
|
---|
792 | Value -> (ee/sw*Sqrt[zetaBdR*xibpw/gamma0bpw])/Sqrt[2],
|
---|
793 | InteractionOrder -> {QED, 1},
|
---|
794 | Description -> "BdW coupling (right-handed)"},
|
---|
795 |
|
---|
796 | KBsRw == {
|
---|
797 | ParameterType -> Internal,
|
---|
798 | BlockName -> Kappa,
|
---|
799 | ComplexParameter -> False,
|
---|
800 | Value -> (gw*Sqrt[zetaBsR*xibpw/gamma0bpw])/Sqrt[2],
|
---|
801 | InteractionOrder -> {QED, 1},
|
---|
802 | Description -> "BsW coupling (right-handed)"},
|
---|
803 |
|
---|
804 | KBbRw == {
|
---|
805 | ParameterType -> Internal,
|
---|
806 | BlockName -> Kappa,
|
---|
807 | ComplexParameter -> False,
|
---|
808 | Value -> (gw*Sqrt[zetaBbR*xibpw/gamma0bpw])/Sqrt[2],
|
---|
809 | InteractionOrder -> {QED, 1},
|
---|
810 | Description -> "BbW coupling (right-handed)"},
|
---|
811 |
|
---|
812 | KBdLz == {
|
---|
813 | ParameterType -> Internal,
|
---|
814 | BlockName -> Kappa,
|
---|
815 | ComplexParameter -> False,
|
---|
816 | Value -> (gw*Sqrt[zetaBdL*xibpz/gamma0bpz])/2/cw,
|
---|
817 | InteractionOrder -> {QED, 1},
|
---|
818 | Description -> "BdZ coupling (left-handed)"},
|
---|
819 |
|
---|
820 | KBsLz == {
|
---|
821 | ParameterType -> Internal,
|
---|
822 | BlockName -> Kappa,
|
---|
823 | ComplexParameter -> False,
|
---|
824 | Value -> (gw*Sqrt[zetaBsL*xibpz/gamma0bpz])/2/cw,
|
---|
825 | InteractionOrder -> {QED, 1},
|
---|
826 | Description -> "BsZ coupling (left-handed)"},
|
---|
827 |
|
---|
828 | KBbLz == {
|
---|
829 | ParameterType -> Internal,
|
---|
830 | BlockName -> Kappa,
|
---|
831 | ComplexParameter -> False,
|
---|
832 | Value -> (gw*Sqrt[zetaBbL*xibpz/gamma0bpz])/2/cw,
|
---|
833 | InteractionOrder -> {QED, 1},
|
---|
834 | Description -> "BbZ coupling (left-handed)"},
|
---|
835 |
|
---|
836 | KBdRz == {
|
---|
837 | ParameterType -> Internal,
|
---|
838 | BlockName -> Kappa,
|
---|
839 | ComplexParameter -> False,
|
---|
840 | Value -> (gw*Sqrt[zetaBdR*xibpz/gamma0bpz])/2/cw,
|
---|
841 | InteractionOrder -> {QED, 1},
|
---|
842 | Description -> "BdZ coupling (right-handed)"},
|
---|
843 |
|
---|
844 | KBsRz == {
|
---|
845 | ParameterType -> Internal,
|
---|
846 | BlockName -> Kappa,
|
---|
847 | ComplexParameter -> False,
|
---|
848 | Value -> (gw*Sqrt[zetaBsR*xibpz/gamma0bpz])/2/cw,
|
---|
849 | InteractionOrder -> {QED, 1},
|
---|
850 | Description -> "BsZ coupling (right-handed)"},
|
---|
851 |
|
---|
852 | KBbRz == {
|
---|
853 | ParameterType -> Internal,
|
---|
854 | BlockName -> Kappa,
|
---|
855 | ComplexParameter -> False,
|
---|
856 | Value -> (gw*Sqrt[zetaBbR*xibpz/gamma0bpz])/2/cw,
|
---|
857 | InteractionOrder -> {QED, 1},
|
---|
858 | Description -> "BbZ coupling (right-handed)"},
|
---|
859 |
|
---|
860 | KBdLh == {
|
---|
861 | ParameterType -> Internal,
|
---|
862 | BlockName -> Kappa,
|
---|
863 | ComplexParameter -> False,
|
---|
864 | Value -> (Sqrt[zetaBdL*xibph/gamma0bph]),
|
---|
865 | InteractionOrder -> {QED, 0},
|
---|
866 | Description -> "BdH coupling (left-handed)"},
|
---|
867 |
|
---|
868 | KBsLh == {
|
---|
869 | ParameterType -> Internal,
|
---|
870 | BlockName -> Kappa,
|
---|
871 | ComplexParameter -> False,
|
---|
872 | Value -> (Sqrt[zetaBsL*xibph/gamma0bph]),
|
---|
873 | InteractionOrder -> {QED, 0},
|
---|
874 | Description -> "BsH coupling (left-handed)"},
|
---|
875 |
|
---|
876 | KBbLh == {
|
---|
877 | ParameterType -> Internal,
|
---|
878 | BlockName -> Kappa,
|
---|
879 | ComplexParameter -> False,
|
---|
880 | Value -> (Sqrt[zetaBbL*xibph/gamma0bph]),
|
---|
881 | InteractionOrder -> {QED, 0},
|
---|
882 | Description -> "BbH coupling (left-handed)"},
|
---|
883 |
|
---|
884 | KBdRh == {
|
---|
885 | ParameterType -> Internal,
|
---|
886 | BlockName -> Kappa,
|
---|
887 | ComplexParameter -> False,
|
---|
888 | Value -> (Sqrt[zetaBdR*xibph/gamma0bph]),
|
---|
889 | InteractionOrder -> {QED, 0},
|
---|
890 | Description -> "BdH coupling (right-handed)"},
|
---|
891 |
|
---|
892 | KBsRh == {
|
---|
893 | ParameterType -> Internal,
|
---|
894 | BlockName -> Kappa,
|
---|
895 | ComplexParameter -> False,
|
---|
896 | Value -> (Sqrt[zetaBsR*xibph/gamma0bph]),
|
---|
897 | InteractionOrder -> {QED, 0},
|
---|
898 | Description -> "BsH coupling (right-handed)"},
|
---|
899 |
|
---|
900 | KBbRh == {
|
---|
901 | ParameterType -> Internal,
|
---|
902 | BlockName -> Kappa,
|
---|
903 | ComplexParameter -> False,
|
---|
904 | Value -> (Sqrt[zetaBbR*xibph/gamma0bph]),
|
---|
905 | InteractionOrder -> {QED, 0},
|
---|
906 | Description -> "BbH coupling (right-handed)"},
|
---|
907 |
|
---|
908 | (* Internal Width functions *)
|
---|
909 |
|
---|
910 | gamma0tpw == {
|
---|
911 | ParameterType -> Internal,
|
---|
912 | BlockName -> WIDTH,
|
---|
913 | ComplexParameter -> False,
|
---|
914 | Value -> (1-MW^2/MTP^2)*(1+MW^2/MTP^2-2*MW^4/MTP^4),
|
---|
915 | Description -> "T partial width for T>Wq (massless q)"},
|
---|
916 |
|
---|
917 | gamma0tpz == {
|
---|
918 | ParameterType -> Internal,
|
---|
919 | BlockName -> WIDTH,
|
---|
920 | ComplexParameter -> False,
|
---|
921 | Value -> 1/2*(1-MZ^2/MTP^2)*(1+MZ^2/MTP^2-2*MZ^4/MTP^4),
|
---|
922 | Description -> "T partial width for T>Zq (massless q)"},
|
---|
923 |
|
---|
924 | gamma0tph == {
|
---|
925 | ParameterType -> Internal,
|
---|
926 | BlockName -> WIDTH,
|
---|
927 | ComplexParameter -> False,
|
---|
928 | Value -> 1/2*(1-MH^2/MTP^2)^2,
|
---|
929 | Description -> "T partial width for T>Hq (massless q)"},
|
---|
930 |
|
---|
931 | gamma0bpw == {
|
---|
932 | ParameterType -> Internal,
|
---|
933 | BlockName -> WIDTH,
|
---|
934 | ComplexParameter -> False,
|
---|
935 | Value -> (1-MW^2/MBP^2)*(1+MW^2/MBP^2-2*MW^4/MBP^4),
|
---|
936 | Description -> "B partial width for B>Wq (massless q)"},
|
---|
937 |
|
---|
938 | gamma0bpz == {
|
---|
939 | ParameterType -> Internal,
|
---|
940 | BlockName -> WIDTH,
|
---|
941 | ComplexParameter -> False,
|
---|
942 | Value -> 1/2*(1-MZ^2/MBP^2)*(1+MZ^2/MBP^2-2*MZ^4/MBP^4),
|
---|
943 | Description -> "B partial width for B>Zq (massless q)"},
|
---|
944 |
|
---|
945 | gamma0bph == {
|
---|
946 | ParameterType -> Internal,
|
---|
947 | BlockName -> WIDTH,
|
---|
948 | ComplexParameter -> False,
|
---|
949 | Value -> 1/2*(1-MH^2/MBP^2)^2,
|
---|
950 | Description -> "B partial width for B>Hq (massless q)"},
|
---|
951 |
|
---|
952 | gamma0xw == {
|
---|
953 | ParameterType -> Internal,
|
---|
954 | BlockName -> WIDTH,
|
---|
955 | ComplexParameter -> False,
|
---|
956 | Value -> (1-MW^2/MX^2)*(1+MW^2/MX^2-2*MW^4/MX^4),
|
---|
957 | Description -> "X partial width for X>Wq (massless q)"},
|
---|
958 |
|
---|
959 | gamma0yw == {
|
---|
960 | ParameterType -> Internal,
|
---|
961 | BlockName -> WIDTH,
|
---|
962 | ComplexParameter -> False,
|
---|
963 | Value -> (1-MW^2/MY^2)*(1+MW^2/MY^2-2*MW^4/MY^4),
|
---|
964 | Description -> "Y partial width for Y>Wq (massless q)"}}
|
---|
965 |
|
---|
966 | (************** Gauge Groups ******************)
|
---|
967 |
|
---|
968 | M$GaugeGroups = {
|
---|
969 |
|
---|
970 | U1Y == {
|
---|
971 | Abelian -> True,
|
---|
972 | GaugeBoson -> B,
|
---|
973 | Charge -> Y,
|
---|
974 | CouplingConstant -> g1},
|
---|
975 |
|
---|
976 | SU2L == {
|
---|
977 | Abelian -> False,
|
---|
978 | GaugeBoson -> Wi,
|
---|
979 | StructureConstant -> Eps,
|
---|
980 | CouplingConstant -> gw},
|
---|
981 |
|
---|
982 | SU3C == {
|
---|
983 | Abelian -> False,
|
---|
984 | GaugeBoson -> G,
|
---|
985 | StructureConstant -> f,
|
---|
986 | SymmetricTensor -> dSUN,
|
---|
987 | Representations -> {T, Colour},
|
---|
988 | CouplingConstant -> gs}
|
---|
989 | }
|
---|
990 |
|
---|
991 | (********* Particle Classes **********)
|
---|
992 |
|
---|
993 | M$ClassesDescription = {
|
---|
994 |
|
---|
995 | (********** Fermions ************)
|
---|
996 | (* Leptons (neutrino): I_3 = +1/2, Q = 0 *)
|
---|
997 | F[1] == {
|
---|
998 | ClassName -> vl,
|
---|
999 | ClassMembers -> {ve,vm,vt},
|
---|
1000 | FlavorIndex -> Generation,
|
---|
1001 | SelfConjugate -> False,
|
---|
1002 | Indices -> {Index[Generation]},
|
---|
1003 | Mass -> 0,
|
---|
1004 | Width -> 0,
|
---|
1005 | QuantumNumbers -> {LeptonNumber -> 1},
|
---|
1006 | PropagatorLabel -> {"v", "ve", "vm", "vt"} ,
|
---|
1007 | PropagatorType -> S,
|
---|
1008 | PropagatorArrow -> Forward,
|
---|
1009 | PDG -> {12,14,16},
|
---|
1010 | FullName -> {"Electron-neutrino", "Mu-neutrino", "Tau-neutrino"} },
|
---|
1011 |
|
---|
1012 | (* Leptons (electron): I_3 = -1/2, Q = -1 *)
|
---|
1013 | F[2] == {
|
---|
1014 | ClassName -> l,
|
---|
1015 | ClassMembers -> {e, m, tt},
|
---|
1016 | FlavorIndex -> Generation,
|
---|
1017 | SelfConjugate -> False,
|
---|
1018 | Indices -> {Index[Generation]},
|
---|
1019 | Mass -> {Ml, {Me, 5.11 * 10^(-4)}, {MM, 0.10566}, {MTA, 1.777}},
|
---|
1020 | Width -> 0,
|
---|
1021 | QuantumNumbers -> {Q -> -1, LeptonNumber -> 1},
|
---|
1022 | PropagatorLabel -> {"l", "e", "m", "tt"},
|
---|
1023 | PropagatorType -> Straight,
|
---|
1024 | ParticleName -> {"e-", "m-", "tt-"},
|
---|
1025 | AntiParticleName -> {"e+", "m+", "tt+"},
|
---|
1026 | PropagatorArrow -> Forward,
|
---|
1027 | PDG -> {11, 13, 15},
|
---|
1028 | FullName -> {"Electron", "Muon", "Tau"} },
|
---|
1029 |
|
---|
1030 | (* Quarks (u): I_3 = +1/2, Q = +2/3 *)
|
---|
1031 | F[3] == {
|
---|
1032 | ClassMembers -> {u, c, t},
|
---|
1033 | ClassName -> uq,
|
---|
1034 | FlavorIndex -> Generation,
|
---|
1035 | SelfConjugate -> False,
|
---|
1036 | Indices -> {Index[Generation], Index[Colour]},
|
---|
1037 | Mass -> {Mu, {MU, 2.55*10^(-3)}, {MC, 1.40}, {MT, 174.3}},
|
---|
1038 | Width -> {0, 0, {WT, 1.51013490}},
|
---|
1039 | QuantumNumbers -> {Q -> 2/3},
|
---|
1040 | PropagatorLabel -> {"uq", "u", "c", "t"},
|
---|
1041 | PropagatorType -> Straight,
|
---|
1042 | PropagatorArrow -> Forward,
|
---|
1043 | PDG -> {2, 4, 6},
|
---|
1044 | FullName -> {"u-quark", "c-quark", "t-quark"}},
|
---|
1045 |
|
---|
1046 | (* Quarks (d): I_3 = -1/2, Q = -1/3 *)
|
---|
1047 | F[4] == {
|
---|
1048 | ClassMembers -> {d, s, b},
|
---|
1049 | ClassName -> dq,
|
---|
1050 | FlavorIndex -> Generation,
|
---|
1051 | SelfConjugate -> False,
|
---|
1052 | Indices -> {Index[Generation], Index[Colour]},
|
---|
1053 | Mass -> {Md, {MD, 5.04*10^(-3)}, {MS, 0.101}, {MB, 4.2}},
|
---|
1054 | Width -> 0,
|
---|
1055 | QuantumNumbers -> {Q -> -1/3},
|
---|
1056 | PropagatorLabel -> {"dq", "d", "s", "b"},
|
---|
1057 | PropagatorType -> Straight,
|
---|
1058 | PropagatorArrow -> Forward,
|
---|
1059 | PDG -> {1,3,5},
|
---|
1060 | FullName -> {"d-quark", "s-quark", "b-quark"} },
|
---|
1061 |
|
---|
1062 | (* VLQ Quarks X, Q=5/3*)
|
---|
1063 | F[5] == {
|
---|
1064 | ClassMembers -> {x},
|
---|
1065 | ClassName -> xq,
|
---|
1066 | SelfConjugate -> False,
|
---|
1067 | Indices -> {Index[Colour]},
|
---|
1068 | Mass -> {{MX,600}},
|
---|
1069 | Width -> {{WX, 1}},
|
---|
1070 | QuantumNumbers -> {Q -> 5/3},
|
---|
1071 | PropagatorLabel -> {"x"},
|
---|
1072 | PropagatorType -> Straight,
|
---|
1073 | PropagatorArrow -> Forward,
|
---|
1074 | PDG -> {6000005},
|
---|
1075 | FullName -> {"X-quark"}},
|
---|
1076 |
|
---|
1077 | (* VLQ Quarks T, Q=2/3 *)
|
---|
1078 | F[6] == {
|
---|
1079 | ClassName -> tpq,
|
---|
1080 | ClassMembers -> {tp},
|
---|
1081 | SelfConjugate -> False,
|
---|
1082 | Indices -> {Index[Colour]},
|
---|
1083 | Mass -> {{MTP,600}},
|
---|
1084 | Width -> {{WTP,1}},
|
---|
1085 | QuantumNumbers -> {Q -> 2/3},
|
---|
1086 | PropagatorLabel -> {"tp"},
|
---|
1087 | PropagatorType -> Straight,
|
---|
1088 | PropagatorArrow -> Forward,
|
---|
1089 | PDG -> {6000006},
|
---|
1090 | FullName -> {"T-quark"}},
|
---|
1091 |
|
---|
1092 | (* VLQ Quarks B, Q=-1/3 *)
|
---|
1093 | F[7] == {
|
---|
1094 | ClassName -> bpq,
|
---|
1095 | ClassMembers -> {bp},
|
---|
1096 | SelfConjugate -> False,
|
---|
1097 | Indices -> {Index[Colour]},
|
---|
1098 | Mass -> {{MBP,600}},
|
---|
1099 | Width -> {{WBP, 1}},
|
---|
1100 | QuantumNumbers -> {Q -> -1/3},
|
---|
1101 | PropagatorLabel -> {"bp"},
|
---|
1102 | PropagatorType -> Straight,
|
---|
1103 | PropagatorArrow -> Forward,
|
---|
1104 | PDG -> {6000007},
|
---|
1105 | FullName -> {"B-quark"}},
|
---|
1106 |
|
---|
1107 | (* VLQ Quarks Y, Q=-4/3 *)
|
---|
1108 | F[8] == {
|
---|
1109 | ClassMembers -> {y},
|
---|
1110 | ClassName -> yq,
|
---|
1111 | SelfConjugate -> False,
|
---|
1112 | Indices -> {Index[Colour]},
|
---|
1113 | Mass -> {{MY,600}},
|
---|
1114 | Width -> {{WY, 1}},
|
---|
1115 | QuantumNumbers -> {Q -> -4/3},
|
---|
1116 | PropagatorLabel -> {"y"},
|
---|
1117 | PropagatorType -> Straight,
|
---|
1118 | PropagatorArrow -> Forward,
|
---|
1119 | PDG -> {6000008},
|
---|
1120 | FullName -> {"Y-quark"}},
|
---|
1121 |
|
---|
1122 | (********** Ghosts **********)
|
---|
1123 | U[1] == {
|
---|
1124 | ClassName -> ghA,
|
---|
1125 | SelfConjugate -> False,
|
---|
1126 | Indices -> {},
|
---|
1127 | Ghost -> A,
|
---|
1128 | Mass -> 0,
|
---|
1129 | QuantumNumbers -> {GhostNumber -> 1},
|
---|
1130 | PropagatorLabel -> uA,
|
---|
1131 | PropagatorType -> GhostDash,
|
---|
1132 | PropagatorArrow -> Forward},
|
---|
1133 |
|
---|
1134 | U[2] == {
|
---|
1135 | ClassName -> ghZ,
|
---|
1136 | SelfConjugate -> False,
|
---|
1137 | Indices -> {},
|
---|
1138 | Mass -> {MZ, 91.1876},
|
---|
1139 | Ghost -> Z,
|
---|
1140 | QuantumNumbers -> {GhostNumber -> 1},
|
---|
1141 | PropagatorLabel -> uZ,
|
---|
1142 | PropagatorType -> GhostDash,
|
---|
1143 | PropagatorArrow -> Forward},
|
---|
1144 |
|
---|
1145 | U[31] == {
|
---|
1146 | ClassName -> ghWp,
|
---|
1147 | SelfConjugate -> False,
|
---|
1148 | Indices -> {},
|
---|
1149 | Mass -> {MW, Internal},
|
---|
1150 | Ghost -> W,
|
---|
1151 | QuantumNumbers -> {Q-> 1, GhostNumber -> 1},
|
---|
1152 | PropagatorLabel -> uWp,
|
---|
1153 | PropagatorType -> GhostDash,
|
---|
1154 | PropagatorArrow -> Forward},
|
---|
1155 |
|
---|
1156 | U[32] == {
|
---|
1157 | ClassName -> ghWm,
|
---|
1158 | SelfConjugate -> False,
|
---|
1159 | Indices -> {},
|
---|
1160 | Mass -> {MW, Internal},
|
---|
1161 | Ghost -> Wbar,
|
---|
1162 | QuantumNumbers -> {Q-> -1, GhostNumber -> 1},
|
---|
1163 | PropagatorLabel -> uWm,
|
---|
1164 | PropagatorType -> GhostDash,
|
---|
1165 | PropagatorArrow -> Forward},
|
---|
1166 |
|
---|
1167 | U[4] == {
|
---|
1168 | ClassName -> ghG,
|
---|
1169 | SelfConjugate -> False,
|
---|
1170 | Indices -> {Index[Gluon]},
|
---|
1171 | Ghost -> G,
|
---|
1172 | Mass -> 0,
|
---|
1173 | QuantumNumbers -> {GhostNumber -> 1},
|
---|
1174 | PropagatorLabel -> uG,
|
---|
1175 | PropagatorType -> GhostDash,
|
---|
1176 | PropagatorArrow -> Forward},
|
---|
1177 |
|
---|
1178 | U[5] == {
|
---|
1179 | ClassName -> ghWi,
|
---|
1180 | Unphysical -> True,
|
---|
1181 | Definitions -> {ghWi[1] -> (ghWp + ghWm)/Sqrt[2],
|
---|
1182 | ghWi[2] -> (ghWm - ghWp)/Sqrt[2]/I,
|
---|
1183 | ghWi[3] -> cw ghZ + sw ghA},
|
---|
1184 | SelfConjugate -> False,
|
---|
1185 | Ghost -> Wi,
|
---|
1186 | Indices -> {Index[SU2W]},
|
---|
1187 | FlavorIndex -> SU2W},
|
---|
1188 |
|
---|
1189 | U[6] == {
|
---|
1190 | ClassName -> ghB,
|
---|
1191 | SelfConjugate -> False,
|
---|
1192 | Definitions -> {ghB -> -sw ghZ + cw ghA},
|
---|
1193 | Indices -> {},
|
---|
1194 | Ghost -> B,
|
---|
1195 | Unphysical -> True},
|
---|
1196 |
|
---|
1197 | (************ Gauge Bosons ***************)
|
---|
1198 | (* Gauge bosons: Q = 0 *)
|
---|
1199 | V[1] == {
|
---|
1200 | ClassName -> A,
|
---|
1201 | SelfConjugate -> True,
|
---|
1202 | Indices -> {},
|
---|
1203 | Mass -> 0,
|
---|
1204 | Width -> 0,
|
---|
1205 | PropagatorLabel -> "a",
|
---|
1206 | PropagatorType -> W,
|
---|
1207 | PropagatorArrow -> None,
|
---|
1208 | PDG -> 22,
|
---|
1209 | FullName -> "Photon" },
|
---|
1210 |
|
---|
1211 | V[2] == {
|
---|
1212 | ClassName -> Z,
|
---|
1213 | SelfConjugate -> True,
|
---|
1214 | Indices -> {},
|
---|
1215 | Mass -> {MZ, 91.1876},
|
---|
1216 | Width -> {WZ, 2.44639985},
|
---|
1217 | PropagatorLabel -> "Z",
|
---|
1218 | PropagatorType -> Sine,
|
---|
1219 | PropagatorArrow -> None,
|
---|
1220 | PDG -> 23,
|
---|
1221 | FullName -> "Z" },
|
---|
1222 |
|
---|
1223 | (* Gauge bosons: Q = -1 *)
|
---|
1224 | V[3] == {
|
---|
1225 | ClassName -> W,
|
---|
1226 | SelfConjugate -> False,
|
---|
1227 | Indices -> {},
|
---|
1228 | Mass -> {MW, Internal},
|
---|
1229 | Width -> {WW, 2.03535570},
|
---|
1230 | QuantumNumbers -> {Q -> 1},
|
---|
1231 | PropagatorLabel -> "W",
|
---|
1232 | PropagatorType -> Sine,
|
---|
1233 | PropagatorArrow -> Forward,
|
---|
1234 | ParticleName ->"W+",
|
---|
1235 | AntiParticleName ->"W-",
|
---|
1236 | PDG -> 24,
|
---|
1237 | FullName -> "W" },
|
---|
1238 |
|
---|
1239 | V[4] == {
|
---|
1240 | ClassName -> G,
|
---|
1241 | SelfConjugate -> True,
|
---|
1242 | Indices -> {Index[Gluon]},
|
---|
1243 | Mass -> 0,
|
---|
1244 | Width -> 0,
|
---|
1245 | PropagatorLabel -> G,
|
---|
1246 | PropagatorType -> C,
|
---|
1247 | PropagatorArrow -> None,
|
---|
1248 | PDG -> 21,
|
---|
1249 | FullName -> "G" },
|
---|
1250 |
|
---|
1251 | V[5] == {
|
---|
1252 | ClassName -> Wi,
|
---|
1253 | Unphysical -> True,
|
---|
1254 | Definitions -> {Wi[mu_, 1] -> (W[mu] + Wbar[mu])/Sqrt[2],
|
---|
1255 | Wi[mu_, 2] -> (Wbar[mu] - W[mu])/Sqrt[2]/I,
|
---|
1256 | Wi[mu_, 3] -> cw Z[mu] + sw A[mu]},
|
---|
1257 | SelfConjugate -> True,
|
---|
1258 | Indices -> {Index[SU2W]},
|
---|
1259 | FlavorIndex -> SU2W,
|
---|
1260 | Mass -> 0,
|
---|
1261 | PDG -> {1,2,3}},
|
---|
1262 |
|
---|
1263 | V[6] == {
|
---|
1264 | ClassName -> B,
|
---|
1265 | SelfConjugate -> True,
|
---|
1266 | Definitions -> {B[mu_] -> -sw Z[mu] + cw A[mu]},
|
---|
1267 | Indices -> {},
|
---|
1268 | Mass -> 0,
|
---|
1269 | Unphysical -> True},
|
---|
1270 |
|
---|
1271 |
|
---|
1272 | (************ Scalar Fields **********)
|
---|
1273 | (* physical Higgs: Q = 0 *)
|
---|
1274 | S[1] == {
|
---|
1275 | ClassName -> H,
|
---|
1276 | SelfConjugate -> True,
|
---|
1277 | Mass -> {MH, 125},
|
---|
1278 | Width -> {WH, 0.00679485838},
|
---|
1279 | PropagatorLabel -> "H",
|
---|
1280 | PropagatorType -> D,
|
---|
1281 | PropagatorArrow -> None,
|
---|
1282 | PDG -> 25,
|
---|
1283 | TeXParticleName -> "\\phi",
|
---|
1284 | TeXClassName -> "\\phi",
|
---|
1285 | FullName -> "H" },
|
---|
1286 |
|
---|
1287 | S[2] == {
|
---|
1288 | ClassName -> phi,
|
---|
1289 | SelfConjugate -> True,
|
---|
1290 | Mass -> {MZ, 91.5445065},
|
---|
1291 | Width -> Wphi,
|
---|
1292 | PropagatorLabel -> "Phi",
|
---|
1293 | PropagatorType -> D,
|
---|
1294 | PropagatorArrow -> None,
|
---|
1295 | ParticleName ->"phi0",
|
---|
1296 | PDG -> 250,
|
---|
1297 | FullName -> "Phi",
|
---|
1298 | Goldstone -> Z },
|
---|
1299 |
|
---|
1300 | S[3] == {
|
---|
1301 | ClassName -> phi2,
|
---|
1302 | SelfConjugate -> False,
|
---|
1303 | Mass -> {MW, Internal},
|
---|
1304 | Width -> Wphi2,
|
---|
1305 | PropagatorLabel -> "Phi2",
|
---|
1306 | PropagatorType -> D,
|
---|
1307 | PropagatorArrow -> None,
|
---|
1308 | ParticleName ->"phi+",
|
---|
1309 | AntiParticleName ->"phi-",
|
---|
1310 | PDG -> 251,
|
---|
1311 | FullName -> "Phi2",
|
---|
1312 | TeXClassName -> "\\phi^+",
|
---|
1313 | TeXParticleName -> "\\phi^+",
|
---|
1314 | TeXAntiParticleName -> "\\phi^-",
|
---|
1315 | Goldstone -> W,
|
---|
1316 | QuantumNumbers -> {Q -> 1}}
|
---|
1317 | }
|
---|
1318 |
|
---|
1319 |
|
---|
1320 | (*****************************************************************************************)
|
---|
1321 |
|
---|
1322 | (* SM Lagrangian *)
|
---|
1323 |
|
---|
1324 | (******************** Gauge F^2 Lagrangian terms*************************)
|
---|
1325 | (*Sign convention from Lagrangian in between Eq. (A.9) and Eq. (A.10) of Peskin & Schroeder.*)
|
---|
1326 | LGauge = -1/4 (del[Wi[nu, i1], mu] - del[Wi[mu, i1], nu] + gw Eps[i1, i2, i3] Wi[mu, i2] Wi[nu, i3])*
|
---|
1327 | (del[Wi[nu, i1], mu] - del[Wi[mu, i1], nu] + gw Eps[i1, i4, i5] Wi[mu, i4] Wi[nu, i5]) -
|
---|
1328 |
|
---|
1329 | 1/4 (del[B[nu], mu] - del[B[mu], nu])^2 -
|
---|
1330 |
|
---|
1331 | 1/4 (del[G[nu, a1], mu] - del[G[mu, a1], nu] + gs f[a1, a2, a3] G[mu, a2] G[nu, a3])*
|
---|
1332 | (del[G[nu, a1], mu] - del[G[mu, a1], nu] + gs f[a1, a4, a5] G[mu, a4] G[nu, a5]);
|
---|
1333 |
|
---|
1334 |
|
---|
1335 | (********************* Fermion Lagrangian terms*************************)
|
---|
1336 | (*Sign convention from Lagrangian in between Eq. (A.9) and Eq. (A.10) of Peskin & Schroeder.*)
|
---|
1337 | LFermions = Module[{Lkin, LQCD, LEWleft, LEWright},
|
---|
1338 |
|
---|
1339 | Lkin = I uqbar.Ga[mu].del[uq, mu] +
|
---|
1340 | I dqbar.Ga[mu].del[dq, mu] +
|
---|
1341 | I lbar.Ga[mu].del[l, mu] +
|
---|
1342 | I vlbar.Ga[mu].del[vl, mu];
|
---|
1343 |
|
---|
1344 | LQCD = gs (uqbar.Ga[mu].T[a].uq +
|
---|
1345 | dqbar.Ga[mu].T[a].dq)G[mu, a];
|
---|
1346 |
|
---|
1347 | LBright =
|
---|
1348 | -2ee/cw B[mu]/2 lbar.Ga[mu].ProjP.l + (*Y_lR=-2*)
|
---|
1349 | 4ee/3/cw B[mu]/2 uqbar.Ga[mu].ProjP.uq - (*Y_uR=4/3*)
|
---|
1350 | 2ee/3/cw B[mu]/2 dqbar.Ga[mu].ProjP.dq; (*Y_dR=-2/3*)
|
---|
1351 |
|
---|
1352 | LBleft =
|
---|
1353 | -ee/cw B[mu]/2 vlbar.Ga[mu].ProjM.vl - (*Y_LL=-1*)
|
---|
1354 | ee/cw B[mu]/2 lbar.Ga[mu].ProjM.l + (*Y_LL=-1*)
|
---|
1355 | ee/3/cw B[mu]/2 uqbar.Ga[mu].ProjM.uq + (*Y_QL=1/3*)
|
---|
1356 | ee/3/cw B[mu]/2 dqbar.Ga[mu].ProjM.dq ; (*Y_QL=1/3*)
|
---|
1357 |
|
---|
1358 | LWleft = ee/sw/2(
|
---|
1359 | vlbar.Ga[mu].ProjM.vl Wi[mu, 3] - (*sigma3 = ( 1 0 )*)
|
---|
1360 | lbar.Ga[mu].ProjM.l Wi[mu, 3] + (* ( 0 -1 )*)
|
---|
1361 |
|
---|
1362 | Sqrt[2] vlbar.Ga[mu].ProjM.l W[mu] +
|
---|
1363 | Sqrt[2] lbar.Ga[mu].ProjM.vl Wbar[mu]+
|
---|
1364 |
|
---|
1365 | uqbar.Ga[mu].ProjM.uq Wi[mu, 3] - (*sigma3 = ( 1 0 )*)
|
---|
1366 | dqbar.Ga[mu].ProjM.dq Wi[mu, 3] + (* ( 0 -1 )*)
|
---|
1367 |
|
---|
1368 | Sqrt[2] uqbar.Ga[mu].ProjM.CKM.dq W[mu] +
|
---|
1369 | Sqrt[2] dqbar.Ga[mu].ProjM.HC[CKM].uq Wbar[mu]
|
---|
1370 | );
|
---|
1371 |
|
---|
1372 | Lkin + LQCD + LBright + LBleft + LWleft];
|
---|
1373 |
|
---|
1374 |
|
---|
1375 | (** Note : Modifications to the SM W and Z currents should be considered here above **)
|
---|
1376 |
|
---|
1377 | (******************** Higgs Lagrangian terms****************************)
|
---|
1378 | Phi := If[FeynmanGauge, {-I phi2, (v + H + I phi)/Sqrt[2]}, {0, (v + H)/Sqrt[2]}];
|
---|
1379 | Phibar := If[FeynmanGauge, {I phi2bar, (v + H - I phi)/Sqrt[2]} ,{0, (v + H)/Sqrt[2]}];
|
---|
1380 |
|
---|
1381 |
|
---|
1382 |
|
---|
1383 | LHiggs := Block[{PMVec, WVec, Dc, Dcbar, Vphi},
|
---|
1384 |
|
---|
1385 | PMVec = Table[PauliSigma[i], {i, 3}];
|
---|
1386 | Wvec[mu_] := {Wi[mu, 1], Wi[mu, 2], Wi[mu, 3]};
|
---|
1387 |
|
---|
1388 | (*Y_phi=1*)
|
---|
1389 | Dc[f_, mu_] := I del[f, mu] + ee/cw B[mu]/2 f + ee/sw/2 (Wvec[mu].PMVec).f;
|
---|
1390 | Dcbar[f_, mu_] := -I del[f, mu] + ee/cw B[mu]/2 f + ee/sw/2 f.(Wvec[mu].PMVec);
|
---|
1391 | Vphi[Phi_, Phibar_] := -muH^2 Phibar.Phi + \[Lambda] (Phibar.Phi)^2;
|
---|
1392 |
|
---|
1393 | (Dcbar[Phibar, mu]).Dc[Phi, mu] - Vphi[Phi, Phibar]];
|
---|
1394 |
|
---|
1395 |
|
---|
1396 | (*************** Yukawa Lagrangian***********************)
|
---|
1397 | LYuk := If[FeynmanGauge,
|
---|
1398 |
|
---|
1399 | Module[{s,r,n,m,i}, -
|
---|
1400 | yd[m] CKM[n,m] uqbar[s,n,i].ProjP[s,r].dq[r,m,i] (-I phi2) -
|
---|
1401 | yd[n] dqbar[s,n,i].ProjP[s,r].dq[r,n,i] (v+H +I phi)/Sqrt[2] -
|
---|
1402 |
|
---|
1403 | yu[n] uqbar[s,n,i].ProjP[s,r].uq[r,n,i] (v+H -I phi)/Sqrt[2] + (*This sign from eps matrix*)
|
---|
1404 | yu[m] Conjugate[CKM[m,n]] dqbar[s,n,i].ProjP[s,r].uq[r,m,i] ( I phi2bar) -
|
---|
1405 |
|
---|
1406 | yl[n] vlbar[s,n].ProjP[s,r].l[r,n] (-I phi2) -
|
---|
1407 | yl[n] lbar[s,n].ProjP[s,r].l[r,n] (v+H +I phi)/Sqrt[2]
|
---|
1408 | ],
|
---|
1409 |
|
---|
1410 | Module[{s,r,n,m,i}, -
|
---|
1411 | yd[n] dqbar[s,n,i].ProjP[s,r].dq[r,n,i] (v+H)/Sqrt[2] -
|
---|
1412 | yu[n] uqbar[s,n,i].ProjP[s,r].uq[r,n,i] (v+H)/Sqrt[2] -
|
---|
1413 | yl[n] lbar[s,n].ProjP[s,r].l[r,n] (v+H)/Sqrt[2]
|
---|
1414 | ]
|
---|
1415 | ];
|
---|
1416 |
|
---|
1417 | LYukawa := LYuk + HC[LYuk];
|
---|
1418 |
|
---|
1419 | (** Note : Modifications to the SM H currents should be considered here above **)
|
---|
1420 |
|
---|
1421 | (**************Ghost terms**************************)
|
---|
1422 | (* Now we need the ghost terms which are of the form: *)
|
---|
1423 | (* - g * antighost * d_BRST G *)
|
---|
1424 | (* where d_BRST G is BRST transform of the gauge fixing function. *)
|
---|
1425 |
|
---|
1426 | LGhost := If[FeynmanGauge,
|
---|
1427 | Block[{dBRSTG,LGhostG,dBRSTWi,LGhostWi,dBRSTB,LGhostB},
|
---|
1428 |
|
---|
1429 | (***********First the pure gauge piece.**********************)
|
---|
1430 | dBRSTG[mu_,a_] := 1/gs Module[{a2, a3}, del[ghG[a], mu] + gs f[a,a2,a3] G[mu,a2] ghG[a3]];
|
---|
1431 | LGhostG := - gs ghGbar[a].del[dBRSTG[mu,a],mu];
|
---|
1432 |
|
---|
1433 | dBRSTWi[mu_,i_] := sw/ee Module[{i2, i3}, del[ghWi[i], mu] + ee/sw Eps[i,i2,i3] Wi[mu,i2] ghWi[i3] ];
|
---|
1434 | LGhostWi := - ee/sw ghWibar[a].del[dBRSTWi[mu,a],mu];
|
---|
1435 |
|
---|
1436 | dBRSTB[mu_] := cw/ee del[ghB, mu];
|
---|
1437 | LGhostB := - ee/cw ghBbar.del[dBRSTB[mu],mu];
|
---|
1438 |
|
---|
1439 | (***********Next the piece from the scalar field.************)
|
---|
1440 | LGhostphi := - ee/(2*sw*cw) MW ( - I phi2 ( (cw^2-sw^2)ghWpbar.ghZ + 2sw*cw ghWpbar.ghA ) +
|
---|
1441 | I phi2bar ( (cw^2-sw^2)ghWmbar.ghZ + 2sw*cw ghWmbar.ghA ) ) -
|
---|
1442 | ee/(2*sw) MW ( ( (v+H) + I phi) ghWpbar.ghWp + ( (v+H) - I phi) ghWmbar.ghWm ) -
|
---|
1443 | I ee/(2*sw) MZ ( - phi2bar ghZbar.ghWp + phi2 ghZbar.ghWm ) -
|
---|
1444 | ee/(2*sw*cw) MZ (v+H) ghZbar.ghZ ;
|
---|
1445 |
|
---|
1446 |
|
---|
1447 | (***********Now add the pieces together.********************)
|
---|
1448 | LGhostG + LGhostWi + LGhostB + LGhostphi]
|
---|
1449 |
|
---|
1450 | ,
|
---|
1451 |
|
---|
1452 | (*If unitary gauge, only include the gluonic ghost.*)
|
---|
1453 | Block[{dBRSTG,LGhostG},
|
---|
1454 |
|
---|
1455 | (***********First the pure gauge piece.**********************)
|
---|
1456 | dBRSTG[mu_,a_] := 1/gs Module[{a2, a3}, del[ghG[a], mu] + gs f[a,a2,a3] G[mu,a2] ghG[a3]];
|
---|
1457 | LGhostG := - gs ghGbar[a].del[dBRSTG[mu,a],mu];
|
---|
1458 |
|
---|
1459 | (***********Now add the pieces together.********************)
|
---|
1460 | LGhostG]
|
---|
1461 |
|
---|
1462 | ];
|
---|
1463 |
|
---|
1464 | (*********SM Lagrangian*******)
|
---|
1465 | LSM := LGauge + LHiggs + LFermions + LYukawa + LGhost;
|
---|
1466 |
|
---|
1467 |
|
---|
1468 | (*********VLQ Lagrangians*******)
|
---|
1469 | (** We assume that the physical and mass eigenstates match for vector-like quarks **)
|
---|
1470 |
|
---|
1471 | (*********LT, EW interactions*******)
|
---|
1472 |
|
---|
1473 | LTW :=
|
---|
1474 | +KT*KTuLw*(tpbar.W[mu].Ga[mu].ProjM.d)+KT*KTuRw*(tpbar.W[mu].Ga[mu].ProjP.d)+KT*KTuLw*(dbar.Wbar[mu].Ga[mu].ProjM.tp)+KT*KTuRw*(dbar.Wbar[mu].Ga[mu].ProjP.tp)+KT*KTcLw*(tpbar.W[mu].Ga[mu].ProjM.s)+KT*KTcRw*(tpbar.W[mu].Ga[mu].ProjP.s)+KT*KTcLw*(sbar.Wbar[mu].Ga[mu].ProjM.tp)+KT*KTcRw*(sbar.Wbar[mu].Ga[mu].ProjP.tp)+KT*KTtLw*(tpbar.W[mu].Ga[mu].ProjM.b)+KT*KTtRw*(tpbar.W[mu].Ga[mu].ProjP.b)+KT*KTtLw*(bbar.Wbar[mu].Ga[mu].ProjM.tp)+KT*KTtRw*(bbar.Wbar[mu].Ga[mu].ProjP.tp);
|
---|
1475 |
|
---|
1476 | LTZ :=+KT*KTuLz*(tpbar.Z[mu].Ga[mu].ProjM.u)+KT*KTuRz*(tpbar.Z[mu].Ga[mu].ProjP.u)+KT*KTuLz*(ubar.Z[mu].Ga[mu].ProjM.tp)+KT*KTuRz*(ubar.Z[mu].Ga[mu].ProjP.tp)+KT*KTcLz*(tpbar.Z[mu].Ga[mu].ProjM.c)+KT*KTcRz*(tpbar.Z[mu].Ga[mu].ProjP.c)+KT*KTcLz*(cbar.Z[mu].Ga[mu].ProjM.tp)+KT*KTcRz*(cbar.Z[mu].Ga[mu].ProjP.tp)+KT*KTtLz*(tpbar.Z[mu].Ga[mu].ProjM.t)+KT*KTtRz*(tpbar.Z[mu].Ga[mu].ProjP.t)+KT*KTtLz*(tbar.Z[mu].Ga[mu].ProjM.tp)+KT*KTtRz*(tbar.Z[mu].Ga[mu].ProjP.tp);
|
---|
1477 |
|
---|
1478 | LTH:=-KT*MTP*KTuLh*(tpbar.H.ProjP.u)/v-KT*MTP*KTuLh*(ubar.H.ProjM.tp)/v-KT*MTP*KTuRh*(tpbar.H.ProjM.u)/v-KT*MTP*KTuRh*(ubar.H.ProjP.tp)/v-KT*MTP*KTcLh*(tpbar.H.ProjP.c)/v-KT*MTP*KTcLh*(cbar.H.ProjM.tp)/v-KT*MTP*KTcRh*(tpbar.H.ProjM.c)/v-KT*MTP*KTcRh*(cbar.H.ProjP.tp)/v-KT*MTP*KTtLh*(tpbar.H.ProjP.t)/v-KT*MTP*KTtLh*(tbar.H.ProjM.tp)/v-KT*MTP*KTtRh*(tpbar.H.ProjM.t)/v-KT*MTP*KTtRh*(tbar.H.ProjP.tp)/v;
|
---|
1479 |
|
---|
1480 |
|
---|
1481 |
|
---|
1482 | (*********LB, EW interactions*******)
|
---|
1483 |
|
---|
1484 | LBW :=+KB*KBdLw*(bpbar.Wbar[mu].Ga[mu].ProjM.u)+KB*KBdRw*(bpbar.Wbar[mu].Ga[mu].ProjP.u)+KB*KBdLw*(ubar.W[mu].Ga[mu].ProjM.bp)+KB*KBdRw*(ubar.W[mu].Ga[mu].ProjP.bp)+KB*KBsLw*(bpbar.Wbar[mu].Ga[mu].ProjM.c)+KB*KBsRw*(bpbar.Wbar[mu].Ga[mu].ProjP.c)+KB*KBsLw*(cbar.W[mu].Ga[mu].ProjM.bp)+KB*KBsRw*(cbar.W[mu].Ga[mu].ProjP.bp)+KB*KBbLw*(bpbar.Wbar[mu].Ga[mu].ProjM.t)+KB*KBbRw*(bpbar.Wbar[mu].Ga[mu].ProjP.t)+KB*KBbLw*(tbar.W[mu].Ga[mu].ProjM.bp)+KB*KBbRw*(tbar.W[mu].Ga[mu].ProjP.bp);
|
---|
1485 |
|
---|
1486 | LBZ := +KB*KBdLz*(bpbar.Z[mu].Ga[mu].ProjM.d)+KB*KBdRz*(bpbar.Z[mu].Ga[mu].ProjP.d)+KB*KBdLz*(dbar.Z[mu].Ga[mu].ProjM.bp)+KB*KBdRz*(dbar.Z[mu].Ga[mu].ProjP.bp)+KB*KBsLz*(bpbar.Z[mu].Ga[mu].ProjM.s)+KB*KBsRz*(bpbar.Z[mu].Ga[mu].ProjP.s)+KB*KBsLz*(sbar.Z[mu].Ga[mu].ProjM.bp)+KB*KBsRz*(sbar.Z[mu].Ga[mu].ProjP.bp)+KB*KBbLz*(bpbar.Z[mu].Ga[mu].ProjM.b)+KB*KBbRz*(bpbar.Z[mu].Ga[mu].ProjP.b)+KB*KBbLz*(bbar.Z[mu].Ga[mu].ProjM.bp)+KB*KBbRz*(bbar.Z[mu].Ga[mu].ProjP.bp);
|
---|
1487 |
|
---|
1488 | LBH:=-KB*MBP*KBdLh*(bpbar.H.ProjP.d)/v-KB*MBP*KBdLh*(dbar.H.ProjM.bp)/v-KB*MBP*KBdRh*(bpbar.H.ProjM.d)/v-KB*MBP*KBdRh*(dbar.H.ProjP.bp)/v-KB*MBP*KBsLh*(bpbar.H.ProjP.s)/v-KB*MBP*KBsLh*(sbar.H.ProjM.bp)/v-KB*MBP*KBsRh*(bpbar.H.ProjM.s)/v-KB*MBP*KBsRh*(sbar.H.ProjP.bp)/v-KB*MBP*KBbLh*(bpbar.H.ProjP.b)/v-KB*MBP*KBbLh*(bbar.H.ProjM.bp)/v-KB*MBP*KBbRh*(bpbar.H.ProjM.b)/v-KB*MBP*KBbRh*(bbar.H.ProjP.bp)/v;
|
---|
1489 |
|
---|
1490 | (*********LX, EW interactions*******)
|
---|
1491 |
|
---|
1492 |
|
---|
1493 | LXW :=
|
---|
1494 | KX*KXuL*(xbar.W[mu].Ga[mu].ProjM.u)+KX*KXuR*(xbar.W[mu].Ga[mu].ProjP.u)+KX*KXuL*(ubar.Wbar[mu].Ga[mu].ProjM.x)+KX*KXuR*(ubar.Wbar[mu].Ga[mu].ProjP.x)+KX*KXcL*(xbar.W[mu].Ga[mu].ProjM.c)+KX*KXcR*(xbar.W[mu].Ga[mu].ProjP.c)+KX*KXcL*(cbar.Wbar[mu].Ga[mu].ProjM.x)+KX*KXcR*(cbar.Wbar[mu].Ga[mu].ProjP.x)+KX*KXtL*(xbar.W[mu].Ga[mu].ProjM.t)+KX*KXtR*(xbar.W[mu].Ga[mu].ProjP.t)+KX*KXtL*(tbar.Wbar[mu].Ga[mu].ProjM.x)+KX*KXtR*(tbar.Wbar[mu].Ga[mu].ProjP.x);
|
---|
1495 |
|
---|
1496 |
|
---|
1497 |
|
---|
1498 | (*********LY, EW interactions*******)
|
---|
1499 |
|
---|
1500 | LYW :=
|
---|
1501 | +KY*KYdL*(ybar.Wbar[mu].Ga[mu].ProjM.d)+KY*KYdR*(ybar.Wbar[mu].Ga[mu].ProjP.d)+KY*KYdL*(dbar.W[mu].Ga[mu].ProjM.y)+KY*KYdR*(dbar.W[mu].Ga[mu].ProjP.y)+KY*KYsL*(ybar.Wbar[mu].Ga[mu].ProjM.s)+KY*KYsR*(ybar.Wbar[mu].Ga[mu].ProjP.s)+KY*KYsL*(sbar.W[mu].Ga[mu].ProjM.y)+KY*KYsR*(sbar.W[mu].Ga[mu].ProjP.y)+KY*KYbL*(ybar.Wbar[mu].Ga[mu].ProjM.b)+KY*KYbR*(ybar.Wbar[mu].Ga[mu].ProjP.b)+KY*KYbL*(bbar.W[mu].Ga[mu].ProjM.y)+KY*KYbR*(bbar.W[mu].Ga[mu].ProjP.y);
|
---|
1502 |
|
---|
1503 |
|
---|
1504 | (*********Kinetic, mass & QCD lagrangians for VLQ*******)
|
---|
1505 |
|
---|
1506 | LTK := I tpbar.Ga[mu].del[tp, mu];
|
---|
1507 | LBK := I bpbar.Ga[mu].del[bp, mu];
|
---|
1508 | LXK := I xbar.Ga[mu].del[x, mu];
|
---|
1509 | LYK := I ybar.Ga[mu].del[y, mu];
|
---|
1510 |
|
---|
1511 | LTM := -MTP.tpbar.tp;
|
---|
1512 | LBM := -MBP.bpbar.bp;
|
---|
1513 | LXM := -MX.xbar.x;
|
---|
1514 | LYM := -MY.ybar.y;
|
---|
1515 |
|
---|
1516 |
|
---|
1517 | LTG := gs (tpbar.Ga[mu].T[a].tp)G[mu, a];
|
---|
1518 | LBG := gs (bpbar.Ga[mu].T[a].bp)G[mu, a];
|
---|
1519 | LXG := gs (xbar.Ga[mu].T[a].x)G[mu, a];
|
---|
1520 | LYG := gs (ybar.Ga[mu].T[a].y)G[mu, a];
|
---|
1521 |
|
---|
1522 |
|
---|
1523 | LTA := 2*ee/3 (tpbar.Ga[mu].tp)A[mu];
|
---|
1524 | LBA := -1*ee/3 (bpbar.Ga[mu].bp)A[mu];
|
---|
1525 | LXA := 5*ee/3 (xbar.Ga[mu].x)A[mu];
|
---|
1526 | LYA := -4*ee/3 (ybar.Ga[mu].y)A[mu];
|
---|
1527 |
|
---|
1528 |
|
---|
1529 | LT := LTW + LTZ + LTH + LTK + LTM + LTG +LTA ;
|
---|
1530 | LB := LBW + LBZ + LBH + LBK + LBM + LBG +LBA ;
|
---|
1531 | LX := LXW + LXK + LXM + LXG + LXA ;
|
---|
1532 | LY := LYW + LYK + LYM + LYG + LYA ;
|
---|
1533 |
|
---|
1534 | LVLQ := LT + LB + LX + LY;
|
---|
1535 |
|
---|
1536 |
|
---|
1537 |
|
---|
1538 | (*********Total Lagrangian*******)
|
---|
1539 |
|
---|
1540 | L := LSM + LVLQ;
|
---|
1541 |
|
---|
1542 |
|
---|