BSMCharacterisation: SM_MB.fr

File SM_MB.fr, 24.7 KB (added by Benjamin Fuks, 9 years ago)

SM FeynRules model file slightly modified for the MB model requirements

Line 
1(***************************************************************************************************************)
2(****** This is the FeynRules mod-file for the Standard model ******)
3(****** ******)
4(****** Authors: N. Christensen, C. Duhr, B. Fuks ******)
5(****** ******)
6(****** Choose whether Feynman gauge is desired. ******)
7(****** If set to False, unitary gauge is assumed. ****)
8(****** Feynman gauge is especially useful for CalcHEP/CompHEP where the calculation is 10-100 times faster. ***)
9(****** Feynman gauge is not supported in MadGraph and Sherpa. ****)
10(***************************************************************************************************************)
11
12(* ************************** *)
13(* ***** Information ***** *)
14(* ************************** *)
15M$ModelName = "Standard Model";
16
17M$Information = {
18 Authors -> {"N. Christensen", "C. Duhr", "B. Fuks"},
19 Version -> "1.4.6",
20 Date -> "15. 04. 2014",
21 Institutions -> {"Michigan State University", "Universite catholique de Louvain (CP3)", "IPHC Strasbourg / University of Strasbourg"},
22 Emails -> {"neil@pa.msu.edu", "claude.duhr@uclouvain.be", "benjamin.fuks@cnrs.in2p3.fr"},
23 URLs -> "http://feynrules.phys.ucl.ac.be/view/Main/StandardModel"
24};
25
26FeynmanGauge = True;
27
28(* ************************** *)
29(* ***** NLO Variables ****** *)
30(******************************)
31
32FR$LoopSwitches = {{Gf, MW}};
33FR$RmDblExt = { ymb -> MB, ymc -> MC, ymdo -> MD, yme -> Me,
34 ymm -> MMU, yms -> MS, ymt -> MT, ymtau -> MTA, ymup -> MU};
35
36(* ************************** *)
37(* ***** Change log ***** *)
38(* ************************** *)
39
40(* v1.4.6: NLO variable added. *)
41(* v1.4.5: Added widths for ghosts. *)
42(* v1.4.4: Changed widths of goldstone bosons to be the same as for the W and Z bosons *)
43(* v1.4.3: Updated conventions for the symmetric structure constants of SU3. *)
44(* v1.4.2: Set FeynmanGauge=True as default again. *)
45(* v1.4: Added SU(2) representation. *)
46(* -> Modification in the field declarations (doublets are added) *)
47(* -> Modification in the Lagrangian (much simpler). *)
48(* v1.3: Added yukawa couplings for all fermions for gauge invariance. *)
49(* Added yukawa couplings for 1st generation fermions to Massless.rst. *)
50(* Updated parameters to PDG 2010. *)
51(* v1.2: Set FeynmanGauge=True as default. *)
52(* Set Gluonic ghosts to be included in both gauges. *)
53(* v1.1: Fixed yukawa couplings in Feynman gauge. *)
54(* Changed yd[n] CKM[n,m] to yd[m] CKM[n,m]. *)
55(* Changed yu[n] Conjugate[CKM[m,n]] to yu[m] Conjugate[CKM[m,n]]. *)
56
57(* ************************** *)
58(* ***** vevs ***** *)
59(* ************************** *)
60M$vevs = { {Phi[2],vev} };
61
62(* ************************** *)
63(* ***** Gauge groups ***** *)
64(* ************************** *)
65M$GaugeGroups = {
66 U1Y == {
67 Abelian -> True,
68 CouplingConstant -> g1,
69 GaugeBoson -> B,
70 Charge -> Y
71 },
72 SU2L == {
73 Abelian -> False,
74 CouplingConstant -> gw,
75 GaugeBoson -> Wi,
76 StructureConstant -> Eps,
77 Representations -> {Ta,SU2D},
78 Definitions -> {Ta[a_,b_,c_]->PauliSigma[a,b,c]/2, FSU2L[i_,j_,k_]:> I Eps[i,j,k]}
79 },
80 SU3C == {
81 Abelian -> False,
82 CouplingConstant -> gs,
83 GaugeBoson -> G,
84 StructureConstant -> f,
85 Representations -> {T,Colour},
86 SymmetricTensor -> dSUN
87 }
88};
89
90
91(* ************************** *)
92(* ***** Indices ***** *)
93(* ************************** *)
94
95IndexRange[Index[SU2W ]] = Unfold[Range[3]];
96IndexRange[Index[SU2D ]] = Unfold[Range[2]];
97IndexRange[Index[Gluon ]] = NoUnfold[Range[8]];
98IndexRange[Index[Colour ]] = NoUnfold[Range[3]];
99IndexRange[Index[Generation]] = Range[3];
100
101IndexStyle[SU2W, j];
102IndexStyle[SU2D, k];
103IndexStyle[Gluon, a];
104IndexStyle[Colour, m];
105IndexStyle[Generation, f];
106
107
108(* ************************** *)
109(* *** Interaction orders *** *)
110(* *** (as used by mg5) *** *)
111(* ************************** *)
112
113M$InteractionOrderHierarchy = {
114 {QCD, 1},
115 {QED, 2}
116};
117
118
119(* ************************** *)
120(* **** Particle classes **** *)
121(* ************************** *)
122M$ClassesDescription = {
123
124(* Gauge bosons: physical vector fields *)
125 V[1] == {
126 ClassName -> A,
127 SelfConjugate -> True,
128 Mass -> 0,
129 Width -> 0,
130 ParticleName -> "a",
131 PDG -> 22,
132 PropagatorLabel -> "a",
133 PropagatorType -> W,
134 PropagatorArrow -> None,
135 FullName -> "Photon"
136 },
137 V[2] == {
138 ClassName -> Z,
139 SelfConjugate -> True,
140 Mass -> {MZ, 91.1876},
141 Width -> {WZ, 2.4952},
142 ParticleName -> "Z",
143 PDG -> 23,
144 PropagatorLabel -> "Z",
145 PropagatorType -> Sine,
146 PropagatorArrow -> None,
147 FullName -> "Z"
148 },
149 V[3] == {
150 ClassName -> W,
151 SelfConjugate -> False,
152 Mass -> {MW, 79.82436},
153 Width -> {WW, 2.085},
154 ParticleName -> "W+",
155 AntiParticleName -> "W-",
156 QuantumNumbers -> {Q -> 1},
157 PDG -> 24,
158 PropagatorLabel -> "W",
159 PropagatorType -> Sine,
160 PropagatorArrow -> Forward,
161 FullName -> "W"
162 },
163 V[4] == {
164 ClassName -> G,
165 SelfConjugate -> True,
166 Indices -> {Index[Gluon]},
167 Mass -> 0,
168 Width -> 0,
169 ParticleName -> "g",
170 PDG -> 21,
171 PropagatorLabel -> "G",
172 PropagatorType -> C,
173 PropagatorArrow -> None,
174 FullName -> "G"
175 },
176
177(* Ghosts: related to physical gauge bosons *)
178 U[1] == {
179 ClassName -> ghA,
180 SelfConjugate -> False,
181 Ghost -> A,
182 QuantumNumbers -> {GhostNumber -> 1},
183 Mass -> 0,
184 Width -> 0,
185 PropagatorLabel -> "uA",
186 PropagatorType -> GhostDash,
187 PropagatorArrow -> Forward
188 },
189 U[2] == {
190 ClassName -> ghZ,
191 SelfConjugate -> False,
192 Ghost -> Z,
193 QuantumNumbers -> {GhostNumber -> 1},
194 Mass -> {MZ,91.1876},
195 Width -> {WZ, 2.4952},
196 PropagatorLabel -> "uZ",
197 PropagatorType -> GhostDash,
198 PropagatorArrow -> Forward
199 },
200 U[31] == {
201 ClassName -> ghWp,
202 SelfConjugate -> False,
203 Ghost -> W,
204 QuantumNumbers -> {GhostNumber -> 1, Q -> 1},
205 Mass -> {MW, 79.82436},
206 Width -> {WW, 2.085},
207 PropagatorLabel -> "uWp",
208 PropagatorType -> GhostDash,
209 PropagatorArrow -> Forward
210 },
211 U[32] == {
212 ClassName -> ghWm,
213 SelfConjugate -> False,
214 Ghost -> Wbar,
215 QuantumNumbers -> {GhostNumber -> 1, Q -> -1},
216 Mass -> {MW, 79.82436},
217 Width -> {WW, 2.085},
218 PropagatorLabel -> "uWm",
219 PropagatorType -> GhostDash,
220 PropagatorArrow -> Forward
221 },
222 U[4] == {
223 ClassName -> ghG,
224 SelfConjugate -> False,
225 Indices -> {Index[Gluon]},
226 Ghost -> G,
227 PDG -> 82,
228 QuantumNumbers ->{GhostNumber -> 1},
229 Mass -> 0,
230 Width -> 0,
231 PropagatorLabel -> "uG",
232 PropagatorType -> GhostDash,
233 PropagatorArrow -> Forward
234 },
235
236(* Gauge bosons: unphysical vector fields *)
237 V[11] == {
238 ClassName -> B,
239 Unphysical -> True,
240 SelfConjugate -> True,
241 Definitions -> { B[mu_] -> -sw Z[mu]+cw A[mu]}
242 },
243 V[12] == {
244 ClassName -> Wi,
245 Unphysical -> True,
246 SelfConjugate -> True,
247 Indices -> {Index[SU2W]},
248 FlavorIndex -> SU2W,
249 Definitions -> { Wi[mu_,1] -> (Wbar[mu]+W[mu])/Sqrt[2], Wi[mu_,2] -> (Wbar[mu]-W[mu])/(I*Sqrt[2]), Wi[mu_,3] -> cw Z[mu] + sw A[mu]}
250 },
251
252(* Ghosts: related to unphysical gauge bosons *)
253 U[11] == {
254 ClassName -> ghB,
255 Unphysical -> True,
256 SelfConjugate -> False,
257 Ghost -> B,
258 Definitions -> { ghB -> -sw ghZ + cw ghA}
259 },
260 U[12] == {
261 ClassName -> ghWi,
262 Unphysical -> True,
263 SelfConjugate -> False,
264 Ghost -> Wi,
265 Indices -> {Index[SU2W]},
266 FlavorIndex -> SU2W,
267 Definitions -> { ghWi[1] -> (ghWp+ghWm)/Sqrt[2], ghWi[2] -> (ghWm-ghWp)/(I*Sqrt[2]), ghWi[3] -> cw ghZ+sw ghA}
268 } ,
269
270(* Fermions: physical fields *)
271 F[1] == {
272 ClassName -> vl,
273 ClassMembers -> {ve,vm,vt},
274 Indices -> {Index[Generation]},
275 FlavorIndex -> Generation,
276 SelfConjugate -> False,
277 Mass -> 0,
278 Width -> 0,
279 QuantumNumbers -> {LeptonNumber -> 1},
280 PropagatorLabel -> {"v", "ve", "vm", "vt"} ,
281 PropagatorType -> S,
282 PropagatorArrow -> Forward,
283 PDG -> {12,14,16},
284 ParticleName -> {"ve","vm","vt"},
285 AntiParticleName -> {"ve~","vm~","vt~"},
286 FullName -> {"Electron-neutrino", "Mu-neutrino", "Tau-neutrino"}
287 },
288 F[2] == {
289 ClassName -> l,
290 ClassMembers -> {e, mu, ta},
291 Indices -> {Index[Generation]},
292 FlavorIndex -> Generation,
293 SelfConjugate -> False,
294 Mass -> {Ml, {Me,5.11*^-4}, {MMU,0.10566}, {MTA,1.777}},
295 Width -> 0,
296 QuantumNumbers -> {Q -> -1, LeptonNumber -> 1},
297 PropagatorLabel -> {"l", "e", "mu", "ta"},
298 PropagatorType -> Straight,
299 PropagatorArrow -> Forward,
300 PDG -> {11, 13, 15},
301 ParticleName -> {"e-", "mu-", "ta-"},
302 AntiParticleName -> {"e+", "mu+", "ta+"},
303 FullName -> {"Electron", "Muon", "Tau"}
304 },
305 F[3] == {
306 ClassName -> uq,
307 ClassMembers -> {u, c, t},
308 Indices -> {Index[Generation], Index[Colour]},
309 FlavorIndex -> Generation,
310 SelfConjugate -> False,
311 Mass -> {Mu, {MU, 2.55*^-3}, {MC,1.27}, {MT,172}},
312 Width -> {0, 0, {WT,1.50833649}},
313 QuantumNumbers -> {Q -> 2/3},
314 PropagatorLabel -> {"uq", "u", "c", "t"},
315 PropagatorType -> Straight,
316 PropagatorArrow -> Forward,
317 PDG -> {2, 4, 6},
318 ParticleName -> {"u", "c", "t" },
319 AntiParticleName -> {"u~", "c~", "t~"},
320 FullName -> {"u-quark", "c-quark", "t-quark"}
321 },
322 F[4] == {
323 ClassName -> dq,
324 ClassMembers -> {d, s, b},
325 Indices -> {Index[Generation], Index[Colour]},
326 FlavorIndex -> Generation,
327 SelfConjugate -> False,
328 Mass -> {Md, {MD,5.04*^-3}, {MS,0.101}, {MB,4.7}},
329 Width -> 0,
330 QuantumNumbers -> {Q -> -1/3},
331 PropagatorLabel -> {"dq", "d", "s", "b"},
332 PropagatorType -> Straight,
333 PropagatorArrow -> Forward,
334 PDG -> {1,3,5},
335 ParticleName -> {"d", "s", "b" },
336 AntiParticleName -> {"d~", "s~", "b~"},
337 FullName -> {"d-quark", "s-quark", "b-quark"}
338 },
339
340(* Fermions: unphysical fields *)
341 F[11] == {
342 ClassName -> LL,
343 Unphysical -> True,
344 Indices -> {Index[SU2D], Index[Generation]},
345 FlavorIndex -> SU2D,
346 SelfConjugate -> False,
347 QuantumNumbers -> {Y -> -1/2},
348 Definitions -> { LL[sp1_,1,ff_] :> Module[{sp2}, ProjM[sp1,sp2] vl[sp2,ff]], LL[sp1_,2,ff_] :> Module[{sp2}, ProjM[sp1,sp2] l[sp2,ff]] }
349 },
350 F[12] == {
351 ClassName -> lR,
352 Unphysical -> True,
353 Indices -> {Index[Generation]},
354 FlavorIndex -> Generation,
355 SelfConjugate -> False,
356 QuantumNumbers -> {Y -> -1},
357 Definitions -> { lR[sp1_,ff_] :> Module[{sp2}, ProjP[sp1,sp2] l[sp2,ff]] }
358 },
359 F[13] == {
360 ClassName -> QL,
361 Unphysical -> True,
362 Indices -> {Index[SU2D], Index[Generation], Index[Colour]},
363 FlavorIndex -> SU2D,
364 SelfConjugate -> False,
365 QuantumNumbers -> {Y -> 1/6},
366 Definitions -> {
367 QL[sp1_,1,ff_,cc_] :> Module[{sp2}, ProjM[sp1,sp2] uq[sp2,ff,cc]],
368 QL[sp1_,2,ff_,cc_] :> Module[{sp2,ff2}, CKM[ff,ff2] ProjM[sp1,sp2] dq[sp2,ff2,cc]] }
369 },
370 F[14] == {
371 ClassName -> uR,
372 Unphysical -> True,
373 Indices -> {Index[Generation], Index[Colour]},
374 FlavorIndex -> Generation,
375 SelfConjugate -> False,
376 QuantumNumbers -> {Y -> 2/3},
377 Definitions -> { uR[sp1_,ff_,cc_] :> Module[{sp2}, ProjP[sp1,sp2] uq[sp2,ff,cc]] }
378 },
379 F[15] == {
380 ClassName -> dR,
381 Unphysical -> True,
382 Indices -> {Index[Generation], Index[Colour]},
383 FlavorIndex -> Generation,
384 SelfConjugate -> False,
385 QuantumNumbers -> {Y -> -1/3},
386 Definitions -> { dR[sp1_,ff_,cc_] :> Module[{sp2}, ProjP[sp1,sp2] dq[sp2,ff,cc]] }
387 },
388
389(* Higgs: physical scalars *)
390 S[1] == {
391 ClassName -> H,
392 SelfConjugate -> True,
393 Mass -> {MH,125},
394 Width -> {WH,0.00407},
395 PropagatorLabel -> "H",
396 PropagatorType -> D,
397 PropagatorArrow -> None,
398 PDG -> 25,
399 ParticleName -> "H",
400 FullName -> "H"
401 },
402
403(* Higgs: physical scalars *)
404 S[2] == {
405 ClassName -> G0,
406 SelfConjugate -> True,
407 Goldstone -> Z,
408 Mass -> {MZ, 91.1876},
409 Width -> {WZ, 2.4952},
410 PropagatorLabel -> "Go",
411 PropagatorType -> D,
412 PropagatorArrow -> None,
413 PDG -> 250,
414 ParticleName -> "G0",
415 FullName -> "G0"
416 },
417 S[3] == {
418 ClassName -> GP,
419 SelfConjugate -> False,
420 Goldstone -> W,
421 Mass -> {MW, 79.82436},
422 QuantumNumbers -> {Q -> 1},
423 Width -> {WW, 2.085},
424 PropagatorLabel -> "GP",
425 PropagatorType -> D,
426 PropagatorArrow -> None,
427 PDG -> 251,
428 ParticleName -> "G+",
429 AntiParticleName -> "G-",
430 FullName -> "GP"
431 },
432
433(* Higgs: unphysical scalars *)
434 S[11] == {
435 ClassName -> Phi,
436 Unphysical -> True,
437 Indices -> {Index[SU2D]},
438 FlavorIndex -> SU2D,
439 SelfConjugate -> False,
440 QuantumNumbers -> {Y -> 1/2},
441 Definitions -> { Phi[1] -> -I GP, Phi[2] -> (vev + H + I G0)/Sqrt[2] }
442 }
443};
444
445
446(* ************************** *)
447(* ***** Gauge ***** *)
448(* ***** Parameters ***** *)
449(* ***** (FeynArts) ***** *)
450(* ************************** *)
451
452GaugeXi[ V[1] ] = GaugeXi[A];
453GaugeXi[ V[2] ] = GaugeXi[Z];
454GaugeXi[ V[3] ] = GaugeXi[W];
455GaugeXi[ V[4] ] = GaugeXi[G];
456GaugeXi[ S[1] ] = 1;
457GaugeXi[ S[2] ] = GaugeXi[Z];
458GaugeXi[ S[3] ] = GaugeXi[W];
459GaugeXi[ U[1] ] = GaugeXi[A];
460GaugeXi[ U[2] ] = GaugeXi[Z];
461GaugeXi[ U[31] ] = GaugeXi[W];
462GaugeXi[ U[32] ] = GaugeXi[W];
463GaugeXi[ U[4] ] = GaugeXi[G];
464
465
466(* ************************** *)
467(* ***** Parameters ***** *)
468(* ************************** *)
469M$Parameters = {
470
471 (* External parameters *)
472 aEWM1 == {
473 ParameterType -> External,
474 BlockName -> SMINPUTS,
475 OrderBlock -> 1,
476 Value -> 127.9,
477 InteractionOrder -> {QED,-2},
478 Description -> "Inverse of the EW coupling constant at the Z pole"
479 },
480 Gf == {
481 ParameterType -> External,
482 BlockName -> SMINPUTS,
483 OrderBlock -> 2,
484 Value -> 1.16637*^-5,
485 InteractionOrder -> {QED,2},
486 TeX -> Subscript[G,f],
487 Description -> "Fermi constant"
488 },
489 aS == {
490 ParameterType -> External,
491 BlockName -> SMINPUTS,
492 OrderBlock -> 3,
493 Value -> 0.1184,
494 InteractionOrder -> {QCD,2},
495 TeX -> Subscript[\[Alpha],s],
496 Description -> "Strong coupling constant at the Z pole"
497 },
498 ymdo == {
499 ParameterType -> External,
500 BlockName -> YUKAWA,
501 OrderBlock -> 1,
502 Value -> 5.04*^-3,
503 Description -> "Down Yukawa mass"
504 },
505 ymup == {
506 ParameterType -> External,
507 BlockName -> YUKAWA,
508 OrderBlock -> 2,
509 Value -> 2.55*^-3,
510 Description -> "Up Yukawa mass"
511 },
512 yms == {
513 ParameterType -> External,
514 BlockName -> YUKAWA,
515 OrderBlock -> 3,
516 Value -> 0.101,
517 Description -> "Strange Yukawa mass"
518 },
519 ymc == {
520 ParameterType -> External,
521 BlockName -> YUKAWA,
522 OrderBlock -> 4,
523 Value -> 1.27,
524 Description -> "Charm Yukawa mass"
525 },
526 ymb == {
527 ParameterType -> External,
528 BlockName -> YUKAWA,
529 OrderBlock -> 5,
530 Value -> 4.7,
531 Description -> "Bottom Yukawa mass"
532 },
533 ymt == {
534 ParameterType -> External,
535 BlockName -> YUKAWA,
536 OrderBlock -> 6,
537 Value -> 172,
538 Description -> "Top Yukawa mass"
539 },
540 yme == {
541 ParameterType -> External,
542 BlockName -> YUKAWA,
543 OrderBlock -> 11,
544 Value -> 5.11*^-4,
545 Description -> "Electron Yukawa mass"
546 },
547 ymm == {
548 ParameterType -> External,
549 BlockName -> YUKAWA,
550 OrderBlock -> 13,
551 Value -> 0.10566,
552 Description -> "Muon Yukawa mass"
553 },
554 ymtau == {
555 ParameterType -> External,
556 BlockName -> YUKAWA,
557 OrderBlock -> 15,
558 Value -> 1.777,
559 Description -> "Tau Yukawa mass"
560 },
561 RCKM== { ParameterType->External, Indices->{Index[Generation],Index[Generation]}, BlockName->VCKM,
562 Value -> {RCKM[i_?NumericQ,i_?NumericQ]->1, RCKM[i_?NumericQ, j_?NumericQ]:>0 /; i=!=j}, Description->"CKM mixing matrix (real part)"},
563 ICKM== { ParameterType->External, Indices->{Index[Generation],Index[Generation]}, BlockName->IMVCKM, Value -> {ICKM[i_,j_] -> 0}, Description->"CKM mixing matrix (imaginary part)"},
564 CKM == { TeX->Superscript[V,ckm], ParameterType->Internal, ComplexParameter->True, Indices->{Index[Generation],Index[Generation]}, Unitary->True,
565 Value->{CKM[i_,j_]:>RCKM[i,j]+I*ICKM[i,j]}, Description-> "CKM mixing matrix"},
566
567 (* Internal Parameters *)
568 aEW == {
569 ParameterType -> Internal,
570 Value -> 1/aEWM1,
571 InteractionOrder -> {QED,2},
572 TeX -> Subscript[\[Alpha], EW],
573 Description -> "Electroweak coupling contant"
574 },
575(* MW == {
576 ParameterType -> Internal,
577 Value -> Sqrt[MZ^2/2+Sqrt[MZ^4/4-Pi/Sqrt[2]*aEW/Gf*MZ^2]],
578 TeX -> Subscript[M,W],
579 Description -> "W mass"
580 },*)
581 sw2 == {
582 ParameterType -> Internal,
583 Value -> 1-(MW/MZ)^2,
584 Description -> "Squared Sin of the Weinberg angle"
585 },
586 ee == {
587 ParameterType -> Internal,
588 Value -> Sqrt[4 Pi aEW],
589 InteractionOrder -> {QED,1},
590 TeX -> e,
591 Description -> "Electric coupling constant"
592 },
593 cw == {
594 ParameterType -> Internal,
595 Value -> Sqrt[1-sw2],
596 TeX -> Subscript[c,w],
597 Description -> "Cosine of the Weinberg angle"
598 },
599 sw == {
600 ParameterType -> Internal,
601 Value -> Sqrt[sw2],
602 TeX -> Subscript[s,w],
603 Description -> "Sine of the Weinberg angle"
604 },
605 gw == {
606 ParameterType -> Internal,
607 Definitions -> {gw->ee/sw},
608 InteractionOrder -> {QED,1},
609 TeX -> Subscript[g,w],
610 Description -> "Weak coupling constant at the Z pole"
611 },
612 g1 == {
613 ParameterType -> Internal,
614 Definitions -> {g1->ee/cw},
615 InteractionOrder -> {QED,1},
616 TeX -> Subscript[g,1],
617 Description -> "U(1)Y coupling constant at the Z pole"
618 },
619 gs == {
620 ParameterType -> Internal,
621 Value -> Sqrt[4 Pi aS],
622 InteractionOrder -> {QCD,1},
623 TeX -> Subscript[g,s],
624 ParameterName -> G,
625 Description -> "Strong coupling constant at the Z pole"
626 },
627 vev == {
628 ParameterType -> Internal,
629 Value -> 2*MW*sw/ee,
630 InteractionOrder -> {QED,-1},
631 Description -> "Higgs vacuum expectation value"
632 },
633 lam == {
634 ParameterType -> Internal,
635 Value -> MH^2/(2*vev^2),
636 InteractionOrder -> {QED, 2},
637 Description -> "Higgs quartic coupling"
638 },
639 muH == {
640 ParameterType -> Internal,
641 Value -> Sqrt[vev^2 lam],
642 TeX -> \[Mu],
643 Description -> "Coefficient of the quadratic piece of the Higgs potential"
644 },
645 yl == {
646 ParameterType -> Internal,
647 Indices -> {Index[Generation], Index[Generation]},
648 Definitions -> {yl[i_?NumericQ, j_?NumericQ] :> 0 /; (i =!= j)},
649 Value -> {yl[1,1] -> Sqrt[2] yme / vev, yl[2,2] -> Sqrt[2] ymm / vev, yl[3,3] -> Sqrt[2] ymtau / vev},
650 InteractionOrder -> {QED, 1},
651 ParameterName -> {yl[1,1] -> ye, yl[2,2] -> ym, yl[3,3] -> ytau},
652 TeX -> Superscript[y, l],
653 Description -> "Lepton Yukawa couplings"
654 },
655 yu == {
656 ParameterType -> Internal,
657 Indices -> {Index[Generation], Index[Generation]},
658 Definitions -> {yu[i_?NumericQ, j_?NumericQ] :> 0 /; (i =!= j)},
659 Value -> {yu[1,1] -> Sqrt[2] ymup/vev, yu[2,2] -> Sqrt[2] ymc/vev, yu[3,3] -> Sqrt[2] ymt/vev},
660 InteractionOrder -> {QED, 1},
661 ParameterName -> {yu[1,1] -> yup, yu[2,2] -> yc, yu[3,3] -> yt},
662 TeX -> Superscript[y, u],
663 Description -> "Up-type Yukawa couplings"
664 },
665 yd == {
666 ParameterType -> Internal,
667 Indices -> {Index[Generation], Index[Generation]},
668 Definitions -> {yd[i_?NumericQ, j_?NumericQ] :> 0 /; (i =!= j)},
669 Value -> {yd[1,1] -> Sqrt[2] ymdo/vev, yd[2,2] -> Sqrt[2] yms/vev, yd[3,3] -> Sqrt[2] ymb/vev},
670 InteractionOrder -> {QED, 1},
671 ParameterName -> {yd[1,1] -> ydo, yd[2,2] -> ys, yd[3,3] -> yb},
672 TeX -> Superscript[y, d],
673 Description -> "Down-type Yukawa couplings"
674 }
675};
676
677(* ************************** *)
678(* ***** Lagrangian ***** *)
679(* ************************** *)
680
681LGauge := Block[{mu,nu,ii,aa},
682 ExpandIndices[-1/4 FS[B,mu,nu] FS[B,mu,nu] - 1/4 FS[Wi,mu,nu,ii] FS[Wi,mu,nu,ii] - 1/4 FS[G,mu,nu,aa] FS[G,mu,nu,aa], FlavorExpand->SU2W]];
683
684LFermions := Block[{mu},
685 ExpandIndices[I*(
686 QLbar.Ga[mu].DC[QL, mu] + LLbar.Ga[mu].DC[LL, mu] + uRbar.Ga[mu].DC[uR, mu] + dRbar.Ga[mu].DC[dR, mu] + lRbar.Ga[mu].DC[lR, mu]),
687 FlavorExpand->{SU2W,SU2D}]/.{CKM[a_,b_] Conjugate[CKM[a_,c_]]->IndexDelta[b,c], CKM[b_,a_] Conjugate[CKM[c_,a_]]->IndexDelta[b,c]}];
688
689LHiggs := Block[{ii,mu, feynmangaugerules},
690 feynmangaugerules = If[Not[FeynmanGauge], {G0|GP|GPbar ->0}, {}];
691
692 ExpandIndices[DC[Phibar[ii],mu] DC[Phi[ii],mu] + muH^2 Phibar[ii] Phi[ii] - lam Phibar[ii] Phi[ii] Phibar[jj] Phi[jj], FlavorExpand->{SU2D,SU2W}]/.feynmangaugerules
693 ];
694
695LYukawa := Block[{sp,ii,jj,cc,ff1,ff2,ff3,yuk,feynmangaugerules},
696 feynmangaugerules = If[Not[FeynmanGauge], {G0|GP|GPbar ->0}, {}];
697
698 yuk = ExpandIndices[
699 -yd[ff2, ff3] CKM[ff1, ff2] QLbar[sp, ii, ff1, cc].dR [sp, ff3, cc] Phi[ii] -
700 yl[ff1, ff3] LLbar[sp, ii, ff1].lR [sp, ff3] Phi[ii] -
701 yu[ff1, ff2] QLbar[sp, ii, ff1, cc].uR [sp, ff2, cc] Phibar[jj] Eps[ii, jj], FlavorExpand -> SU2D];
702 yuk = yuk /. { CKM[a_, b_] Conjugate[CKM[a_, c_]] -> IndexDelta[b, c], CKM[b_, a_] Conjugate[CKM[c_, a_]] -> IndexDelta[b, c]};
703 yuk+HC[yuk]/.feynmangaugerules
704 ];
705
706LGhost := Block[{LGh1,LGhw,LGhs,LGhphi,mu, generators,gh,ghbar,Vectorize,phi1,phi2,togoldstones,doublet,doublet0},
707 (* Pure gauge piece *)
708 LGh1 = -ghBbar.del[DC[ghB,mu],mu];
709 LGhw = -ghWibar.del[DC[ghWi,mu],mu];
710 LGhs = -ghGbar.del[DC[ghG,mu],mu];
711
712 (* Scalar pieces: see Peskin pages 739-742 *)
713 (* phi1 and phi2 are the real degrees of freedom of GP *)
714 (* Vectorize transforms a doublet in a vector in the phi-basis, i.e. the basis of real degrees of freedom *)
715 gh = {ghB, ghWi[1], ghWi[2], ghWi[3]};
716 ghbar = {ghBbar, ghWibar[1], ghWibar[2], ghWibar[3]};
717 generators = {-I/2 g1 IdentityMatrix[2], -I/2 gw PauliSigma[1], -I/2 gw PauliSigma[2], -I/2 gw PauliSigma[3]};
718 doublet = Expand[{(-I phi1 - phi2)/Sqrt[2], Phi[2]} /. MR$Definitions /. vev -> 0];
719 doublet0 = {0, vev/Sqrt[2]};
720 Vectorize[{a_, b_}]:= Simplify[{Sqrt[2] Re[Expand[a]], Sqrt[2] Im[Expand[a]], Sqrt[2] Re[Expand[b]], Sqrt[2] Im[Expand[b]]}/.{Im[_]->0, Re[num_]->num}];
721 togoldstones := {phi1 -> (GP + GPbar)/Sqrt[2], phi2 -> (-GP + GPbar)/(I Sqrt[2])};
722 LGhphi=Plus@@Flatten[Table[-ghbar[[kkk]].gh[[lll]] Vectorize[generators[[kkk]].doublet0].Vectorize[generators[[lll]].(doublet+doublet0)],{kkk,4},{lll,4}]] /.togoldstones;
723
724ExpandIndices[ LGhs + If[FeynmanGauge, LGh1 + LGhw + LGhphi,0], FlavorExpand->SU2W]];
725
726LSM:= LGauge + LFermions + LHiggs + LYukawa + LGhost;