| 1 | (* ********************************************************* *)
|
|---|
| 2 | (* ***** ***** *)
|
|---|
| 3 | (* ***** FeynRules model file: SM + singlet vector boson ***** *)
|
|---|
| 4 | (* ***** Author: B. Fuks, L. Darme ***** *)
|
|---|
| 5 | (* ***** ***** *)
|
|---|
| 6 | (* ********************************************************* *)
|
|---|
| 7 |
|
|---|
| 8 | (* ************************** *)
|
|---|
| 9 | (* ***** Information ***** *)
|
|---|
| 10 | (* ************************** *)
|
|---|
| 11 | M$ModelName = "SingletVector";
|
|---|
| 12 | M$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 | M$InteractionOrderHierarchy = { {QCD, 1}, {QED, 2}, {NP,1} };
|
|---|
| 25 | M$InteractionOrderLimit = { {QCD, 99}, {QED, 99}, {NP,99} };
|
|---|
| 26 |
|
|---|
| 27 | (* ************************** *)
|
|---|
| 28 | (* ***** Fields ***** *)
|
|---|
| 29 | (* ************************** *)
|
|---|
| 30 |
|
|---|
| 31 | M$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 | (* ************************** *)
|
|---|
| 49 | M$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 | (* ************************** *)
|
|---|
| 75 | L1 := ExpandIndices[-1/4 FS[Zp,mu,nu] FS[Zp,mu,nu] + 1/2 MZp^2 Zp[mu] Zp[mu]];
|
|---|
| 76 |
|
|---|
| 77 | LupL := 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] ]];
|
|---|
| 78 | LupR := 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 |
|
|---|
| 80 | LagNP := L1 + LupL + LupR ;
|
|---|