WZPrimeAtNLO: vPrimeNLO.fr

File vPrimeNLO.fr, 7.3 KB (added by Richard Ruiz, 8 years ago)

W' Z' FeynRules model file in Feynman gauge

Line 
1(* ********************************************************* *)
2(* ***** ***** *)
3(* ***** FeynRules model file: SM + W_SSM + Z_SSM ***** *)
4(* ***** Authors: B. Fuks, R. Ruiz ***** *)
5(* ***** ***** *)
6(* ********************************************************* *)
7
8(* ************************** *)
9(* ***** Information ***** *)
10(* ************************** *)
11M$ModelName = "vPrimeNLO";
12M$Information = { Authors->{"B. Fuks, R. Ruiz"},
13 Emails->{"benjamin.fuks@iphc.cnrs.fr,rruiz@durham.ac.uk"},
14 Institutions->{"IPHC Strasbourg / University of Strasbourg, IPPP/University of Durham"},
15 Date->"2017 January 20",
16 Version->"1.1",
17 References->{"B. Fuks, R. Ruiz [arXiv:1701.YYYYY], and references therein."},
18 URLs->{"feynrules.irmp.ucl.ac.be/"} };
19FeynmanGauge = True;
20
21
22
23
24
25(* ************************** *)
26(* ***** Fields ***** *)
27(* ************************** *)
28M$ClassesDescription = {
29(* Sequential SM Z prime boson *)
30 V[32] == {
31 ClassName -> Zp,
32 SelfConjugate -> True,
33 Mass -> {MZp, 3000.00},
34 Width -> {WZp, 89.59},
35 ParticleName -> "Zp",
36 PDG -> 32,
37 PropagatorLabel -> "Zp",
38 PropagatorType -> Sine,
39 PropagatorArrow -> None,
40 FullName -> "Zp"
41 },
42(* Sequential SM W prime boson *)
43 V[34] == {
44 ClassName ->Wp,
45 SelfConjugate ->False,
46 Mass ->{MWp, 3000.00},
47 Width ->{WWp, 101.27},
48 ParticleName ->"Wp+",
49 AntiParticleName->"Wp-",
50 QuantumNumbers ->{Q->1},
51 PDG ->34,
52 PropagatorLabel ->"Wp",
53 PropagatorType ->Sine,
54 PropagatorArrow ->Forward
55 FullName ->"Wp"
56 }
57};
58
59(* ************************** *)
60(* ***** Parameters ***** *)
61(* ************************** *)
62M$Parameters = {
63 CRq == { TeX->Subsuperscript[C,q,R], ParameterType->External, ComplexParameter->False, Indices->{Index[Generation],Index[Generation]}, BlockName->CRq,
64 Value-> {
65 CRq[1,1]->1.0, CRq[1,2]->.225773, CRq[1,3]->0.,
66 CRq[2,1]->0, CRq[2,2]->.97418, CRq[2,3]->0.,
67 CRq[3,1]->0, CRq[3,2]->0., CRq[3,3]->1. },
68 Description->"Right-handed W' couplings to quarks"},
69 CRl == { TeX->Subsuperscript[C,l,R], ParameterType->External, ComplexParameter->False, Indices->{Index[Generation],Index[Generation]}, BlockName->CRl,
70 Value-> {
71 CRl[1,1]->1., CRl[1,2]->0., CRl[1,3]->0.,
72 CRl[2,1]->0., CRl[2,2]->1., CRl[2,3]->0.,
73 CRl[3,1]->0., CRl[3,2]->0., CRl[3,3]->1. },
74 Description->"Right-handed W' couplings to leptons"},
75 CLq == { TeX->Subsuperscript[C,q,L], ParameterType->External, ComplexParameter->False, Indices->{Index[Generation],Index[Generation]}, BlockName->CLq,
76 Value-> {
77 CLq[1,1]->.97418, CLq[1,2]->.225773, CLq[1,3]->0.,
78 CLq[2,1]->-.225773, CLq[2,2]->.97418, CLq[2,3]->0.,
79 CLq[3,1]->0., CLq[3,2]->0., CLq[3,3]->1. },
80 Description->"Left-handed W' couplings to quarks"},
81 CLl == { TeX->Subsuperscript[C,l,L], ParameterType->External, ComplexParameter->False, Indices->{Index[Generation],Index[Generation]}, BlockName->CLl,
82 Value-> {
83 CLl[1,1]->1., CLl[1,2]->0., CLl[1,3]->0.,
84 CLl[2,1]->0., CLl[2,2]->1., CLl[2,3]->0.,
85 CLl[3,1]->0., CLl[3,2]->0., CLl[3,3]->1. },
86 Description->"Left-handed W' couplings to leptons"
87 },
88
89
90 kL == { ParameterType -> External,
91 Value -> 1.0,
92 BlockName->SSMCOUP,
93 OrderBlock->1,
94 TeX -> Subscript[kappa,L],
95 Description -> "Left-handed W' coupling constant scale factor"
96 },
97
98 kR == { ParameterType -> External,
99 Value -> 0.0,
100 TeX -> Subscript[kappa,R],
101 BlockName->SSMCOUP,
102 OrderBlock->2,
103 Description -> "Right-handed W' coupling constant scale factor"
104 },
105
106(* Z_SSM vector coupling to fermions: gZp,f_V = (1/2)TL3 - Qf*sW2 *)
107 gZpuV == {
108 ParameterType -> Internal,
109 Value -> (1/2)*( 1/2) - ( 2/3)*sw2,
110 Description -> "Vector Zp-u-uBar Coupling"
111 },
112 gZpdV == {
113 ParameterType -> Internal,
114 Value -> (1/2)*(-1/2) - (-1/3)*sw2,
115 Description -> "Vector Zp-d-dBar Coupling"
116 },
117 gZpvV == {
118 ParameterType -> Internal,
119 Value -> (1/2)*( 1/2) - (0)*sw2,
120 Description -> "Vector Zp-v-vBar Coupling"
121 },
122 gZpeV == {
123 ParameterType -> Internal,
124 Value -> (1/2)*(-1/2) - (-1)*sw2,
125 Description -> "Vector Zp-e-eBar Coupling"
126 },
127
128(* Z_SSM axial-vector coupling to fermions: gZp,f_A = (-1/2)TL3 *)
129 gZpuA == {
130 ParameterType -> Internal,
131 Value -> (-1/2)*(1/2),
132 Description -> "Axial-vector Zp-u-uBar Coupling"
133 },
134 gZpdA == {
135 ParameterType -> Internal,
136 Value -> (-1/2)*(-1/2),
137 Description -> "Axial-vector Zp-d-dBar Coupling"
138 },
139 gZpvA == {
140 ParameterType -> Internal,
141 Value -> (-1/2)*(1/2),
142 Description -> "Axial-vector Zp-v-vBar Coupling"
143 },
144 gZpeA == {
145 ParameterType -> Internal,
146 Value -> (-1/2)*(-1/2),
147 Description -> "Axial-vector Zp-e-eBar Coupling"
148 },
149
150(* Z_SSM LH coupling to fermions: gZp,f_L = gZpfV - gZpfA *)
151 gZpuL == {
152 ParameterType -> Internal,
153 Value -> gZpuV - gZpuA,
154 Description -> "LH Zp-u-uBar Coupling"
155 },
156 gZpdL == {
157 ParameterType -> Internal,
158 Value -> gZpdV - gZpdA,
159 Description -> "LH Zp-d-dBar Coupling"
160 },
161 gZpvL == {
162 ParameterType -> Internal,
163 Value -> gZpvV - gZpvA,
164 Description -> "LH Zp-v-vBar Coupling"
165 },
166 gZpeL == {
167 ParameterType -> Internal,
168 Value -> gZpeV - gZpeA,
169 Description -> "LH Zp-e-eBar Coupling"
170 },
171
172(* Z_SSM RH coupling to fermions: gZp,f_R = gZpfV + gZpfA *)
173 gZpuR == {
174 ParameterType -> Internal,
175 Value -> gZpuV + gZpuA,
176 Description -> "LH Zp-u-uBar Coupling"
177 },
178 gZpdR == {
179 ParameterType -> Internal,
180 Value -> gZpdV + gZpdA,
181 Description -> "LH Zp-d-dBar Coupling"
182 },
183 gZpvR == {
184 ParameterType -> Internal,
185 Value -> gZpvV + gZpvA,
186 Description -> "LH Zp-v-vBar Coupling"
187 },
188 gZpeR == {
189 ParameterType -> Internal,
190 Value -> gZpeV + gZpeA,
191 Description -> "LH Zp-e-eBar Coupling"
192 }
193
194};
195
196(* ************************** *)
197(* *** Interaction orders *** *)
198(* ************************** *)
199
200(* ZSSM Currents *)
201LZSSMTmpU := -kL*ee/sw/cw *Zp[mu]*(gZpuR* uqbar.Ga[mu].ProjP.uq \
202 + gZpuL* uqbar.Ga[mu].ProjM.uq );
203LZSSMTmpD := -kL*ee/sw/cw *Zp[mu]*(gZpdR* dqbar.Ga[mu].ProjP.dq \
204 + gZpdL* dqbar.Ga[mu].ProjM.dq );
205LZSSMTmpV := -kL*ee/sw/cw *Zp[mu]*(gZpvR* vlbar.Ga[mu].ProjP.vl \
206 + gZpvL* vlbar.Ga[mu].ProjM.vl );
207LZSSMTmpE := -kL*ee/sw/cw *Zp[mu]*(gZpeR* lbar.Ga[mu].ProjP.l \
208 + gZpeL* lbar.Ga[mu].ProjM.l );
209LZSSM := LZSSMTmpU + LZSSMTmpD + LZSSMTmpV + LZSSMTmpE;
210
211(* WSSM Currents *)
212LWSSMTmp := -kL*ee/sw/Sqrt[2]*(vlbar.Ga[mu].ProjM.CLl.l + uqbar.Ga[mu].ProjM.CLq.dq)*Wp[mu] \
213 -kR*ee/sw/Sqrt[2]*(vlbar.Ga[mu].ProjP.CRl.l + uqbar.Ga[mu].ProjP.CRq.dq)*Wp[mu];
214
215LWSSM := LWSSMTmp + HC[LWSSMTmp];
216
217(* Combine Everything *)
218LSSM := LWSSM + LZSSM;
219LFull := LSM + LSSM;