TopHeavyRes: VectorSinglet.fr

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