LeptoQuark: VLferm.fr

File VLferm.fr, 13.2 KB (added by JavierFuentes, 4 years ago)

Changed PDG code for vector-like fermions

Line 
1(* ************************************************************************* *)
2(* ***** ***** *)
3(* ***** FeynRules model file: VLferm ***** *)
4(* ***** Authors: J. Fuentes-Martin and D. Faroughy ***** *)
5(* ***** ***** *)
6(* ************************************************************************* *)
7
8(* ************************** *)
9(* ***** Information ***** *)
10(* ************************** *)
11
12M$ModelName = "VLferm";
13
14M$Information = {Authors -> {"J. Fuentes-Martin", "D. Faroughy"},
15 Version -> "1.0",
16 Date -> "21.06.2022",
17 Institutions -> {"University of Granada", "University of Zurich"},
18 Emails -> {"javier.fuentes@ugr.es", "faroughy@physik.uzh.ch"}
19 };
20
21(* ************************** *)
22(* *** Interaction orders *** *)
23(* ************************** *)
24
25M$InteractionOrderHierarchy = { {QCD,1}, {NP,1}, {QED,2} };
26
27
28(* ************************** *)
29(* **** Particle classes **** *)
30(* ************************** *)
31
32M$ClassesDescription = {
33(* Vector-like fermions: physical fields *)
34 F[5] == {
35 ClassName -> Ln,
36 SelfConjugate -> False,
37 Mass -> {MLn, 1000},
38 Width -> {WLn, 10},
39 QuantumNumbers -> {LeptonNumber -> 1},
40 PropagatorLabel -> "Ln",
41 PropagatorType -> Straight,
42 PropagatorArrow -> Forward,
43 PDG -> 72,
44 ParticleName -> "Ln",
45 AntiParticleName -> "Ln~",
46 FullName -> "Neutral vector-like lepton"
47 },
48 F[6] == {
49 ClassName -> Le,
50 SelfConjugate -> False,
51 Mass -> {MLe, 1000},
52 Width -> {WLe, 10},
53 QuantumNumbers -> {Q -> -1, LeptonNumber -> 1},
54 PropagatorLabel -> "Le",
55 PropagatorType -> Straight,
56 PropagatorArrow -> Forward,
57 PDG -> 71,
58 ParticleName -> "Le-",
59 AntiParticleName -> "Le+",
60 FullName -> "Charged vector-like lepton"
61 },
62 F[7] == {
63 ClassName -> Qu,
64 SelfConjugate -> False,
65 Indices -> {Index[Colour]},
66 Mass -> {MQu, 2000},
67 Width -> {WQu, 10},
68 QuantumNumbers -> {Q -> 2/3},
69 PropagatorLabel -> "Qu",
70 PropagatorType -> Straight,
71 PropagatorArrow -> Forward,
72 PDG -> 62,
73 ParticleName -> "Qu",
74 AntiParticleName -> "Qu~",
75 FullName -> "Up-type vector-like quark"
76 },
77 F[8] == {
78 ClassName -> Qd,
79 SelfConjugate -> False,
80 Indices -> {Index[Colour]},
81 Mass -> {MQd, 2000},
82 Width -> {WQd, 10},
83 QuantumNumbers -> {Q -> -1/3},
84 PropagatorLabel -> "Qd",
85 PropagatorType -> Straight,
86 PropagatorArrow -> Forward,
87 PDG -> 61,
88 ParticleName -> "Qd",
89 AntiParticleName -> "Qd~",
90 FullName -> "Down-type vector-like quark"
91 },
92
93(* Vector-like fermions: unphysical fields *)
94
95 F[16] == {
96 ClassName -> chiL,
97 Unphysical -> True,
98 Indices -> {Index[SU2D]},
99 FlavorIndex -> SU2D,
100 SelfConjugate -> False,
101 QuantumNumbers -> {Y -> -1/2},
102 Definitions -> {
103 chiL[sp1_,1] :> Ln[sp1],
104 chiL[sp1_,2] :> Le[sp1]
105 }
106 },
107 F[17] == {
108 ClassName -> chiQ,
109 Unphysical -> True,
110 Indices -> {Index[SU2D], Index[Colour]},
111 FlavorIndex -> SU2D,
112 SelfConjugate -> False,
113 QuantumNumbers -> {Y -> 1/6},
114 Definitions -> {
115 chiQ[sp1_,1,cc_] :> Qu[sp1,cc],
116 chiQ[sp1_,2,cc_] :> Qd[sp1,cc]
117 }
118 }
119};
120
121(* ************************** *)
122(* ***** Parameters ***** *)
123(* ************************** *)
124
125M$Parameters = {
126(* External Parameters *)
127 betaLQL == {
128 ParameterType -> External,
129 Value -> 1.0,
130 InteractionOrder -> {NP,0},
131 TeX -> Subsuperscript[\[Beta],L,QL],
132 BlockName -> NPLQCOUP,
133 OrderBlock -> 8,
134 Description -> "LQ Q L left-handed coupling"
135 },
136 betaRQL == {
137 ParameterType -> External,
138 Value -> 1.0,
139 InteractionOrder -> {NP,0},
140 TeX -> Subsuperscript[\[Beta],R,QL],
141 BlockName -> NPLQCOUP,
142 OrderBlock -> 9,
143 Description -> "LQ Q L right-handed coupling"
144 },
145 betaL3L == {
146 ParameterType -> External,
147 Value -> 1.0,
148 InteractionOrder -> {NP,0},
149 TeX -> Subsuperscript[\[Beta],L,3L],
150 BlockName -> NPLQCOUP,
151 OrderBlock -> 10,
152 Description -> "LQ q3 L coupling"
153 },
154 betaLQ3 == {
155 ParameterType -> External,
156 Value -> 1.0,
157 InteractionOrder -> {NP,0},
158 TeX -> Subsuperscript[\[Beta],L,Q3],
159 BlockName -> NPLQCOUP,
160 OrderBlock -> 11,
161 Description -> "LQ Q l3 coupling"
162 },
163 betaL2L == {
164 ParameterType -> External,
165 Value -> 0.2,
166 InteractionOrder -> {NP,0},
167 TeX -> Subsuperscript[\[Beta],L,2L],
168 BlockName -> NPLQCOUP,
169 OrderBlock -> 12,
170 Description -> "LQ q2 L coupling"
171 },
172 betaLQ2 == {
173 ParameterType -> External,
174 Value -> 0.2,
175 InteractionOrder -> {NP,0},
176 TeX -> Subsuperscript[\[Beta],L,Q2],
177 BlockName -> NPLQCOUP,
178 OrderBlock -> 13,
179 Description -> "LQ Q l2 coupling"
180 },
181
182 zetaLQQ == {
183 ParameterType -> External,
184 Value -> 1.0,
185 InteractionOrder -> {NP,0},
186 TeX -> Subsuperscript[\[Zeta],L,QQ],
187 BlockName -> NPZPCOUP,
188 OrderBlock -> 13,
189 Description -> "Zp Q Q left-handed coupling"
190 },
191 zetaRQQ == {
192 ParameterType -> External,
193 Value -> 1.0,
194 InteractionOrder -> {NP,0},
195 TeX -> Subsuperscript[\[Zeta],R,QQ],
196 BlockName -> NPZPCOUP,
197 OrderBlock -> 14,
198 Description -> "Zp Q Q right-handed coupling"
199 },
200 zetaLQ2 == {
201 ParameterType -> External,
202 Value -> 0.2,
203 InteractionOrder -> {NP,0},
204 TeX -> Subsuperscript[\[Zeta],L,Q2],
205 BlockName -> NPZPCOUP,
206 OrderBlock -> 15,
207 Description -> "Zp Q q2 left-handed coupling"
208 },
209 zetaLLL == {
210 ParameterType -> External,
211 Value -> 1.0,
212 InteractionOrder -> {NP,0},
213 TeX -> Subsuperscript[\[Zeta],L,LL],
214 BlockName -> NPZPCOUP,
215 OrderBlock -> 16,
216 Description -> "Zp L L left-handed coupling"
217 },
218 zetaRLL == {
219 ParameterType -> External,
220 Value -> 1.0,
221 InteractionOrder -> {NP,0},
222 TeX -> Subsuperscript[\[Zeta],R,LL],
223 BlockName -> NPZPCOUP,
224 OrderBlock -> 17,
225 Description -> "Zp L L right-handed coupling"
226 },
227 zetaLL2 == {
228 ParameterType -> External,
229 Value -> 0.2,
230 InteractionOrder -> {NP,0},
231 TeX -> Subsuperscript[\[Zeta],L,L2],
232 BlockName -> NPZPCOUP,
233 OrderBlock -> 18,
234 Description -> "Zp L l2 left-handed coupling"
235 },
236
237 kappaLQQ == {
238 ParameterType -> External,
239 Value -> 1.0,
240 InteractionOrder -> {NP,0},
241 TeX -> Subsuperscript[\[Kappa],L,QQ],
242 BlockName -> NPGPCOUP,
243 OrderBlock -> 10,
244 Description -> "Gp Q Q left-handed coupling"
245 },
246 kappaRQQ == {
247 ParameterType -> External,
248 Value -> 1.0,
249 InteractionOrder -> {NP,0},
250 TeX -> Subsuperscript[\[Kappa],R,QQ],
251 BlockName -> NPGPCOUP,
252 OrderBlock -> 11,
253 Description -> "Gp Q Q left-handed coupling"
254 },
255 kappaLQ2 == {
256 ParameterType -> External,
257 Value -> 0.2,
258 InteractionOrder -> {NP,0},
259 TeX -> Subsuperscript[\[Kappa],L,Q2],
260 BlockName -> NPGPCOUP,
261 OrderBlock -> 12,
262 Description -> "Gp Q q2 left-handed coupling"
263 },
264
265 (* Internal Parameters *)
266 betaLchiQ == {
267 ParameterType -> Internal,
268 Indices -> {Index[Generation]},
269 Definitions -> {betaLchiQ[1] -> 0},
270 Value -> {betaLchiQ[2] -> betaLQ2, betaLchiQ[3] -> betaLQ3},
271 InteractionOrder -> {NP,0},
272 TeX -> Subsuperscript[\[Beta],L,Ql],
273 Description -> "Q-SM-LQ coupling vector"
274 },
275 betaLchiL == {
276 ParameterType -> Internal,
277 Indices -> {Index[Generation]},
278 Definitions -> {betaLchiL[1] -> 0},
279 Value -> {betaLchiL[2] -> betaL2L, betaLchiL[3] -> betaL3L},
280 InteractionOrder -> {NP,0},
281 TeX -> Subsuperscript[\[Beta],L,qL],
282 Description -> "SM-L-Zp coupling vector"
283 },
284 zetaLchiQ == {
285 ParameterType -> Internal,
286 Indices -> {Index[Generation]},
287 Definitions -> {zetaLchiQ[1] -> 0, zetaLchiQ[3] -> 0},
288 Value -> {zetaLchiQ[2] -> zetaLQ2},
289 InteractionOrder -> {NP,0},
290 TeX -> Subsuperscript[\[Zeta],L,Qq],
291 Description -> "SM-Q-Zp coupling vector"
292 },
293 zetaLchiL == {
294 ParameterType -> Internal,
295 Indices -> {Index[Generation]},
296 Definitions -> {zetaLchiL[1] -> 0, zetaLchiL[3] -> 0},
297 Value -> {zetaLchiL[2] -> zetaLL2},
298 InteractionOrder -> {NP,0},
299 TeX -> Subsuperscript[\[Zeta],L,Ll],
300 Description -> "SM-L-Zp coupling vector"
301 },
302 kappaLchiQ == {
303 ParameterType -> Internal,
304 Indices -> {Index[Generation]},
305 Definitions -> {kappaLchiQ[1] -> 0, kappaLchiQ[3] -> 0},
306 Value -> {kappaLchiQ[2] -> kappaLQ2},
307 InteractionOrder -> {NP,0},
308 TeX -> Subsuperscript[\[Kappa],L,Qq],
309 Description -> "SM-Q-Gp coupling vector"
310 }
311};
312
313
314(* ************************** *)
315(* ***** Lagrangian ***** *)
316(* ************************** *)
317
318LVLfermkin := Block[{mu},
319 ExpandIndices[
320 I*(chiQbar.Ga[mu].DC[chiQ, mu] + chiLbar.Ga[mu].DC[chiL, mu]),
321 FlavorExpand->{SU2W,SU2D}]
322 ];
323
324LU1VL := Block[{ff1,ff2,s1,s2,s3,cc1,ii1,mu},
325 ExpandIndices[
326 gU/Sqrt[2] * VLQ[mu,cc1] *
327 (
328 betaLQL * chiQbar[s1,ii1,cc1] * Ga[mu,s1,s2] * ProjM[s2,s3] * chiL[s3,ii1]
329 + betaRQL * chiQbar[s1,ii1,cc1] * Ga[mu,s1,s2] * ProjP[s2,s3] * chiL[s3,ii1]
330 + betaLchiL[ff1] * newCKM[ff2,ff1] * uqbar[s1,ff2,cc1] * Ga[mu,s1,s2] * ProjM[s2,s3] * Ln[s3]
331 + betaLchiL[ff1] * dqbar[s1,ff1,cc1] * Ga[mu,s1,s2] * ProjM[s2,s3] * Le[s3]
332 + betaLchiQ[ff1] * chiQbar[s1,ii1,cc1] * Ga[mu,s1,s2] * ProjM[s2,s3] * LL[s3,ii1,ff1]
333 )
334 ,FlavorExpand-> SU2D]
335 ];
336
337LZpVL := Block[{s1,s2,s3,mu,cc1,ii1,ff1},
338 ExpandIndices[
339 gZp/(2 * Sqrt[6]) * Zp[mu] *
340 (
341 zetaLQQ * chiQbar[s1,ii1,cc1] * Ga[mu,s1,s2] * ProjM[s2,s3] * chiQ[s3,ii1,cc1]
342 + zetaRQQ * chiQbar[s1,ii1,cc1] * Ga[mu,s1,s2] * ProjP[s2,s3] * chiQ[s3,ii1,cc1]
343 -3 * zetaLLL * chiLbar[s1,ii1] * Ga[mu,s1,s2] * ProjM[s2,s3] * chiL[s3,ii1]
344 -3 * zetaRLL * chiLbar[s1,ii1] * Ga[mu,s1,s2] * ProjP[s2,s3] * chiL[s3,ii1]
345 )
346 ,FlavorExpand-> SU2D]
347 ];
348
349LZpVLMix := Block[{s1,s2,s3,mu,cc1,ii1,ff1,ff2},
350 ExpandIndices[
351 gZp/(2 * Sqrt[6]) * Zp[mu] *
352 (
353 zetaLchiQ[ff1] * HC[newCKM[ff2,ff1]] * Qubar[s1,cc1] * Ga[mu,s1,s2] * ProjM[s2,s3] * uq[s3,ff2,cc1]
354 + zetaLchiQ[ff1] * Qdbar[s1,cc1] * Ga[mu,s1,s2] * ProjM[s2,s3] * dq[s3,ff1,cc1]
355 -3 * zetaLchiL[ff1] * chiLbar[s1,ii1] * Ga[mu,s1,s2] * ProjM[s2,s3] * LL[s3,ii1,ff1]
356 )
357 ,FlavorExpand-> SU2D]
358 ];
359
360LGpVL := Block[{s1,s2,s3,aa1,cc1,cc2,mu,ii1,ff1},
361 ExpandIndices[
362 gGp * Gp[mu,aa1] *
363 (
364 kappaLQQ * chiQbar[s1,ii1,cc1] * T[aa1,cc1,cc2] * Ga[mu,s1,s2] * ProjM[s2,s3] * chiQ[s3,ii1,cc2]
365 + kappaRQQ * chiQbar[s1,ii1,cc1] * T[aa1,cc1,cc2] * Ga[mu,s1,s2] * ProjP[s2,s3] * chiQ[s3,ii1,cc2]
366 )
367 ,FlavorExpand-> SU2D]
368 ];
369
370
371LGpVLMix := Block[{s1,s2,s3,aa1,cc1,cc2,mu,ff1,ff2},
372 ExpandIndices[
373 gGp * Gp[mu,aa1] *
374 (
375 kappaLchiQ[ff1] * HC[newCKM[ff2,ff1]] * Qubar[s1,cc1] * T[aa1,cc1,cc2] * Ga[mu,s1,s2] * ProjM[s2,s3] * uq[s3,ff2,cc2]
376 + kappaLchiQ[ff1] * Qdbar[s1,cc1] * T[aa1,cc1,cc2] * Ga[mu,s1,s2] * ProjM[s2,s3] * dq[s3,ff1,cc2]
377 )
378 ,FlavorExpand-> SU2D]
379 ];
380
381LVLferm := LVLfermkin + LU1VL + HC[LU1VL] + LZpVL + LZpVLMix + HC[LZpVLMix] + LGpVL + LGpVLMix + HC[LGpVLMix];