Wprime: weff.fr

File weff.fr, 4.4 KB (added by Benjamin Fuks, 7 years ago)

the main model file.

Line 
1(* ********************************************************* *)
2(* ***** ***** *)
3(* ***** FeynRules model file: SM + Wprime ***** *)
4(* ***** Author: B. Fuks ***** *)
5(* ***** ***** *)
6(* ********************************************************* *)
7
8(* ************************** *)
9(* ***** Information ***** *)
10(* ************************** *)
11M$ModelName = "WEff";
12M$Information = { Authors->{"Benjamin Fuks"},
13 Emails->{"benjamin.fuks@iphc.cnrs.fr"},
14 Institutions->{"IPHC Strasbourg / University of Strasbourg"},
15 Date->"03.07.12", Version->"1.0.0",
16 References->{"Z. Sullivan, Phys. Rev. D 66 (2002) 075011."},
17 URLs->{""} };
18FeynmanGauge = True;
19
20
21
22
23
24(* ************************** *)
25(* ***** Fields ***** *)
26(* ************************** *)
27M$ClassesDescription = {
28(* Additional W' boson *)
29 V[7] == {
30 ClassName->Wp, SelfConjugate->False,
31 Mass->{MWp,1000.},
32 Width->{WWp,10.},
33 ParticleName->"Wp+",
34 AntiParticleName->"Wp-",
35 QuantumNumbers -> {Q->1},
36 PDG->34,
37 PropagatorLabel->"Wp",
38 PropagatorType->Sine,
39 PropagatorArrow->Forward}
40};
41
42(* ************************** *)
43(* ***** Parameters ***** *)
44(* ************************** *)
45M$Parameters = {
46 CRq == { TeX->Subsuperscript[C,q,R], ParameterType->External, ComplexParameter->False, Indices->{Index[Generation],Index[Generation]}, BlockName->CRq,
47 Value-> {
48 CRq[1,1]->1.0, CRq[1,2]->.225773, CRq[1,3]->0.,
49 CRq[2,1]->0, CRq[2,2]->.97418, CRq[2,3]->0.,
50 CRq[3,1]->0, CRq[3,2]->0., CRq[3,3]->1. },
51 Description->"Right-handed W' couplings to quarks"},
52 CRl == { TeX->Subsuperscript[C,l,R], ParameterType->External, ComplexParameter->False, Indices->{Index[Generation],Index[Generation]}, BlockName->CRl,
53 Value-> {
54 CRl[1,1]->1., CRl[1,2]->0., CRl[1,3]->0.,
55 CRl[2,1]->0., CRl[2,2]->1., CRl[2,3]->0.,
56 CRl[3,1]->0., CRl[3,2]->0., CRl[3,3]->1. },
57 Description->"Right-handed W' couplings to leptons"},
58 CLq == { TeX->Subsuperscript[C,q,L], ParameterType->External, ComplexParameter->False, Indices->{Index[Generation],Index[Generation]}, BlockName->CLq,
59 Value-> {
60 CLq[1,1]->.97418, CLq[1,2]->.225773, CLq[1,3]->0.,
61 CLq[2,1]->-.225773, CLq[2,2]->.97418, CLq[2,3]->0.,
62 CLq[3,1]->0., CLq[3,2]->0., CLq[3,3]->1. },
63 Description->"Left-handed W' couplings to quarks"},
64 CLl == { TeX->Subsuperscript[C,l,L], ParameterType->External, ComplexParameter->False, Indices->{Index[Generation],Index[Generation]}, BlockName->CLl,
65 Value-> {
66 CLl[1,1]->1., CLl[1,2]->0., CLl[1,3]->0.,
67 CLl[2,1]->0., CLl[2,2]->1., CLl[2,3]->0.,
68 CLl[3,1]->0., CLl[3,2]->0., CLl[3,3]->1. },
69 Description->"Left-handed W' couplings to leptons"
70 },
71
72
73 kL == { ParameterType -> External,
74 Value -> 0.0,
75 InteractionOrder -> {NP,1},
76 BlockName->WPCOUP,
77 OrderBlock->1,
78 TeX -> Subscript[kappa,L],
79 Description -> "Left-handed W' coupling constant scale factor"},
80
81 kR == { ParameterType -> External,
82 Value -> 1.0,
83 InteractionOrder -> {NP,1},
84 TeX -> Subscript[kappa,R],
85 BlockName->WPCOUP,
86 OrderBlock->2,
87 Description -> "Right-handed W' coupling constant scale factor"
88 },
89
90
91 gL == { ParameterType -> Internal,
92 Definitions -> {gL->kL*ee/sw},
93 InteractionOrder -> {NP,1},
94 TeX -> Subscript[g,L],
95 Description -> "Left-handed W' coupling constant"
96 },
97
98 gR == { ParameterType -> Internal,
99 Definitions -> {gR->kR*ee/sw},
100 InteractionOrder -> {NP,1},
101 TeX -> Subscript[g,R],
102 Description -> "Right-handed W' coupling constant"}
103};
104
105(* ************************** *)
106(* *** Interaction orders *** *)
107(* ************************** *)
108M$InteractionOrderHierarchy = { {QCD, 1}, {QED, 2} };
109
110
111(* New Lagrangian interaction terms *)
112
113LWprimeTmp := ExpandIndices[
114 gL/Sqrt[2]*(vlbar.Ga[mu].ProjM.CLl.l + uqbar.Ga[mu].ProjM.CLq.dq)*Wp[mu] +
115 gR/Sqrt[2]*(vlbar.Ga[mu].ProjP.CRl.l + uqbar.Ga[mu].ProjP.CRq.dq)*Wp[mu],
116 FlavorExpand -> {SU2D}];
117
118LWprime := LWprimeTmp + HC[LWprimeTmp];
119