TopHeavyRes: VectorOctet.fr

File VectorOctet.fr, 2.7 KB (added by LucDarme, 4 years ago)
Line 
1(* ********************************************************* *)
2(* ***** ***** *)
3(* ***** FeynRules model file: SM + colored vector boson ***** *)
4(* ***** Author: B. Fuks, L. Darme ***** *)
5(* ***** ***** *)
6(* ********************************************************* *)
7
8(* ************************** *)
9(* ***** Information ***** *)
10(* ************************** *)
11M$ModelName = "ColoredVector";
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(* ************************** *)
22(* ***** Orders ***** *)
23(* ************************** *)
24
25M$InteractionOrderHierarchy = { {QCD, 1}, {QED, 2}, {NP,1} };
26M$InteractionOrderLimit = { {QCD, 99}, {QED, 99}, {NP,99} };
27
28(* ************************** *)
29(* ***** Fields ***** *)
30(* ************************** *)
31
32M$ClassesDescription = {
33 V[6] == {
34 ClassName -> Vo,
35 SelfConjugate -> True,
36 Indices -> {Index[Gluon]},
37 Mass -> {MVo, 3000},
38 Width -> {WVo, 500},
39 ParticleName -> "Vo",
40 PropagatorLabel -> "Vo",
41 PropagatorType -> Sine,
42 PropagatorArrow -> None,
43 FullName -> "Voctet"
44 }
45};
46
47
48(* ************************** *)
49(* ***** Parameters ***** *)
50(* ************************** *)
51M$Parameters = {
52
53
54 g8uL == {
55 TeX -> Subscript[g,"8L"],
56 ParameterType -> External,
57 ComplexParameter -> False,
58 BlockName -> NPG8L,
59 InteractionOrder -> {NP,1},
60 Value -> 0.3,
61 InteractionOrder -> {NP,1},
62 Description -> "Vector octet coupling to Left-handed up quark"},
63 g8uR == {
64 TeX -> Subscript[g,"8R"],
65 ParameterType -> External,
66 ComplexParameter -> False,
67 BlockName -> NPG8R,
68 InteractionOrder -> {NP,1},
69 Value -> 0.3,
70 Description -> "Vector octet coupling to Right-handed up quark"}
71};
72
73(* ************************** *)
74(* ***** Lagrangian ***** *)
75(* ************************** *)
76L1 := ExpandIndices[-1/4 (DC[Vo[mu,cc],nu]-DC[Vo[nu,cc],mu])^2 + 1/2 MVo^2 Vo[mu,aa] Vo[mu,aa]];
77
78LupL := Block[{sp,sp2,sp3,ff,gg,mm,nn},ExpandIndices[ Vo[mu,a] T[a,mm,nn] *g8uL*ProjM[sp2,sp3]*Ga[mu,sp,sp2] uqbar[sp,3,mm].uq[sp3,3,nn] ]];
79LupR := Block[{sp,sp2,sp3,ff,gg,mm,nn},ExpandIndices[ Vo[mu,a] T[a,mm,nn] *g8uR*ProjP[sp2,sp3]*Ga[mu,sp,sp2] uqbar[sp,3,mm].uq[sp3,3,nn] ]];
80
81LagNP := L1 + LupL + LupR ;