TopHeavyRes: ScalarOctet.fr

File ScalarOctet.fr, 2.3 KB (added by LucDarme, 3 years ago)
Line 
1(* ********************************************************* *)
2(* ***** ***** *)
3(* ***** FeynRules model file: SM + colored scalars ***** *)
4(* ***** Author: B. Fuks, L. Darme ***** *)
5(* ***** ***** *)
6(* ********************************************************* *)
7
8(* ************************** *)
9(* ***** Information ***** *)
10(* ************************** *)
11M$ModelName = "ColoredScalarsAndPScalars";
12M$Information = {
13 Authors -> {"Luc Darme","Benjamin Fuks"},
14 Emails -> {"luc.darme@lnf.infn.it","fuks@lpthe.jussieu.fr"},
15 Institutions -> {"LPTHE / UPMC","LNF, INFN"},
16 Date -> "26.02.2021",
17 Version -> "1",
18 References -> "xxx"};
19
20(* ************************** *)
21(* ***** Orders ***** *)
22(* ************************** *)
23M$InteractionOrderHierarchy = { {QCD, 1}, {QED, 2}, {NP,1} };
24M$InteractionOrderLimit = { {QCD, 99}, {QED, 99}, {NP,99} };
25
26(* ************************** *)
27(* ***** Fields ***** *)
28(* ************************** *)
29M$ClassesDescription = {
30 S[4] == {
31 ClassName -> sig8,
32 SelfConjugate -> True,
33 Indices -> {Index[Gluon]},
34 Mass -> {m8,400},
35 Width -> {w8,10},
36 ParticleName -> "sig8",
37 PDG -> 9000001,
38 PropagatorLabel -> "sig8",
39 PropagatorType -> ScalarDash,
40 PropagatorArrow -> Forward
41 }
42};
43
44(* ************************** *)
45(* ***** Parameters ***** *)
46(* ************************** *)
47M$Parameters = {
48(**** slguons top interaction ****)
49 yS8 == {
50 TeX -> Subsuperscript[y,"8S"],
51 ComplexParameter -> False,
52 ParameterType -> External,
53 BlockName -> NPyS8,
54 Value -> 0.3,
55 InteractionOrder -> {NP,1},
56 Description -> "Sgluon coupling to quarks"
57 }
58};
59
60(* ************************** *)
61(* ***** Lagrangian ***** *)
62(* ************************** *)
63L1 := 1/2 DC[sig8[aa],mu] DC[sig8[aa],mu] - 1/2 m8^2 sig8[aa] sig8[aa];
64
65L2a := Block[{sp,sp2,sp3,ff,gg,mm,nn,aa},ExpandIndices[ sig8[aa] T[aa,mm,nn] yS8*(ProjM[sp,sp2] + ProjP[sp,sp2]) uqbar[sp,3,mm].uq[sp2,3,nn]]];
66
67(*L2b = g8g*dSUN[a,b,c] sig8[a] FS[G,mu,nu,b] FS[G,mu,nu,c];*)
68
69LagNP := L1 + L2a + HC[Expand[L2a]] ;