ALRM: ALRM_LO.fr

File ALRM_LO.fr, 71.3 KB (added by Mustafa Ashry, 8 months ago)

Model File

Line 
1 (***************************************************************************************************************)
2 (****** This is the FEYNRULES model file for an Alternative Left-Right Symmetric Model @ LO ******)
3 (****** ******)
4 (****** Author: Mustafa Ashry ******)
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 = "Alternative Left-Right Symmetric Model";
16M$Information = {
17 Authors -> {"Mustafa Ashry"},
18 Version -> "1.05",
19 Date -> "14. 03. 2024",
20 Institutions -> {"Department of Mathematics, Faculty of Science, Cairo University, 12613 Giza, Egypt","Center for Fundamental Physics, Zewail City of Science and Technology, 6th of October City, Giza 12578, Egypt"},
21 Emails -> {"mustafa@sci.cu.edu.eg","mashry@zewailcity.edu.eg"},
22 References -> {"M. Ashry and S. Khalil, Phenomenological aspects of a TeV-scale alternative left-right model, Physical Review D 91, 015009 (2015) (arXiv:1310.3315)", "Mustafa Ashry, MSc Thesis, TeV-Scale Left-Right Symmetric Model with Minimal Higgs Sector, Cairo University (2015)"},
23 URLs -> {"https://feynrules.irmp.ucl.ac.be/wiki/ALRM","https://github.com/Mustafa-Ashry/ALRM-FeynRules"}};
24
25 (* ************************** *)
26 (* ***** Change log ***** *)
27 (* ************************** *)
28(* v1.01 14.12.16 : Implemented with broken SU(2)_R symmetry and model-dependent Higgs diagonolization matrix *)
29(* v1.02 14.12.17 : Restored the SU(2)_R symmetry into left-right symmetry *)
30(* v1.03 14.12.18 : Added the the model-dependent hYY and hgg effective Lagrangian *)
31(* v1.04 14.12.19 : Added the Higgs diagonolization and the hYY and hgg effective Lagrangian generically and model-independent *)
32(* v1.05 14.03.24 : Added the generic gauge bosons diagonolization and taking the left and right SU(2) gauge couplings (gL and gR), left and right CKM matrices (CKML and CKMR) and left and right MNS matrices (MNSL and MNSR) NOT equal in general, i.e., generically, nonmanifest left-right symmetry *)
33(* v1.06 14.12.24 : To Add the seesaw mechanism for neutrino masses(?) *)
34
35FeynmanGauge = True;
36
37 (* ************************** *)
38 (* ***** vevs ***** *)
39 (* ************************** *)
40M$vevs = {{Phi[2,2],k},{ChiL[2],vL},{ChiR[2],vR}};
41
42 (* ************************** *)
43 (* ***** Gauge groups ***** *)
44 (* ************************** *)
45M$GaugeGroups = {
46 U1BL == {
47 Abelian -> True,
48 CouplingConstant -> gBL,
49 GaugeBoson -> B,
50 Charge -> BL },
51 SU2L == {
52 Abelian -> False,
53 CouplingConstant -> gL,
54 GaugeBoson -> Wi,
55 StructureConstant -> Eps,
56 Representations -> {LTa,LSU2D},
57 Definitions -> {LTa[a_,b_,c_]->PauliSigma[a,b,c]/2, FSU2L[i_,j_,k_]:> I Eps[i,j,k]} },
58 SU2R == {
59 Abelian -> False,
60 CouplingConstant -> gR,
61 GaugeBoson -> Wpi,
62 StructureConstant -> Eps,
63 Representations -> {RTa,RSU2D},
64 Definitions -> {RTa[a_,b_,c_]->PauliSigma[a,b,c]/2, FSU2R[i_,j_,k_]:> I Eps[i,j,k]} },
65 SU3C == {
66 Abelian -> False,
67 CouplingConstant -> gs,
68 GaugeBoson -> G,
69 StructureConstant -> f,
70 Representations -> {T,Colour},
71 SymmetricTensor -> dSUN } };
72
73 (* ************************** *)
74 (* ***** Indices ***** *)
75 (* ************************** *)
76
77IndexRange[Index[SU2L ]] = Unfold[Range[3]];
78IndexRange[Index[SU2R ]] = Unfold[Range[3]];
79IndexRange[Index[LSU2D ]] = Unfold[Range[2]];
80IndexRange[Index[RSU2D ]] = Unfold[Range[2]];
81IndexRange[Index[Gluon ]] = NoUnfold[Range[8]];
82IndexRange[Index[Colour ]] = NoUnfold[Range[3]];
83IndexRange[Index[Generation]] = Range[3];
84IndexRange[Index[higgs ]] = Range[3];
85IndexRange[Index[ngauge ]] = Range[3];
86IndexRange[Index[bjorkeneff]] = Range[5];
87
88IndexStyle[SU2L, j];
89IndexStyle[SU2R, o];
90IndexStyle[LSU2D, k];
91IndexStyle[RSU2D, p];
92IndexStyle[Gluon, a];
93IndexStyle[Colour, m];
94IndexStyle[Generation, f];
95IndexStyle[higgs, h];
96IndexStyle[ngauge g];
97IndexStyle[bjorkeneff, b];
98
99 (* ************************** *)
100 (* *** Interaction orders *** *)
101 (* *** (as used by mg5) *** *)
102 (* ************************** *)
103
104M$InteractionOrderHierarchy = { {QCD, 1}, {QED, 2}, {HIG, 2}, {HIW, 2} };
105M$InteractionOrderLimit = { {QCD, 99}, {QED, 99}, {HIG, 1}, {HIW, 1} };
106
107 (* ************************** *)
108 (* **** Particle classes **** *)
109 (* ************************** *)
110M$ClassesDescription = {
111
112 (* Gauge bosons: physical vector fields *)
113 V[1] == {
114 ClassName -> A,
115 SelfConjugate -> True,
116 Mass -> 0,
117 Width -> 0,
118 ParticleName -> "A",
119 PDG -> 22,
120 PropagatorLabel -> "A",
121 PropagatorType -> W,
122 PropagatorArrow -> None,
123 FullName -> "Photon" },
124 V[2] == {
125 ClassName -> Z,
126 SelfConjugate -> True,
127 Mass -> {MZ,91.1876},
128 Width -> {WZ,2.52},
129 ParticleName -> "Z",
130 PDG -> 23,
131 PropagatorLabel -> "Z",
132 PropagatorType -> Sine,
133 PropagatorArrow -> None,
134 FullName -> "Z" },
135 V[3] == {
136 ClassName -> W,
137 SelfConjugate -> False,
138 Mass -> {MW,80.379},
139 Width -> {WW,2.11},
140 ParticleName -> "W+",
141 AntiParticleName -> "W-",
142 QuantumNumbers -> {Q -> 1},
143 PDG -> 24,
144 PropagatorLabel -> "W",
145 PropagatorType -> Sine,
146 PropagatorArrow -> Forward,
147 FullName -> "W" },
148 V[4] == {
149 ClassName -> G,
150 SelfConjugate -> True,
151 Indices -> {Index[Gluon]},
152 Mass -> 0,
153 Width -> 0,
154 ParticleName -> "G",
155 PDG -> 21,
156 PropagatorLabel -> "G",
157 PropagatorType -> C,
158 PropagatorArrow -> None,
159 FullName -> "G" },
160 V[5] == {
161 ClassName -> Zp,
162 SelfConjugate -> True,
163 Mass -> {MZp,Internal},
164 Width -> {WZp,13.2},
165 ParticleName -> "Zp",
166 PDG -> 100,
167 PropagatorLabel -> "Zp",
168 PropagatorType -> Sine,
169 PropagatorArrow -> None,
170 FullName -> "Zp"
171 },
172 V[6] == {
173 ClassName -> Wp,
174 SelfConjugate -> False,
175 Mass -> {MWp,Internal},
176 Width -> {WWp,3.04},
177 ParticleName -> "Wp+",
178 AntiParticleName -> "Wp-",
179 QuantumNumbers -> {Q -> 1},
180 PDG -> 101,
181 PropagatorLabel -> "Wp",
182 PropagatorType -> Sine,
183 PropagatorArrow -> Forward,
184 FullName -> "Wp" },
185
186 (* Ghosts: related to physical gauge bosons *)
187 U[1] == {
188 ClassName -> ghA,
189 SelfConjugate -> False,
190 Ghost -> A,
191 QuantumNumbers -> {GhostNumber -> 1},
192 Mass -> 0,
193 Width -> 0,
194 PropagatorLabel -> "uA",
195 PropagatorType -> GhostDash,
196 PropagatorArrow -> Forward },
197 U[2] == {
198 ClassName -> ghZ,
199 SelfConjugate -> False,
200 Ghost -> Z,
201 QuantumNumbers -> {GhostNumber -> 1},
202 Mass -> {MZ,91.1876},
203 Width -> {WZ,2.52},
204 PropagatorLabel -> "uZ",
205 PropagatorType -> GhostDash,
206 PropagatorArrow -> Forward },
207 U[31] == {
208 ClassName -> ghWp,
209 SelfConjugate -> False,
210 Ghost -> W,
211 QuantumNumbers -> {GhostNumber -> 1, Q -> 1},
212 Mass -> {MW,80.385},
213 Width -> {WW,2.11},
214 PropagatorLabel -> "uWp",
215 PropagatorType -> GhostDash,
216 PropagatorArrow -> Forward },
217 U[32] == {
218 ClassName -> ghWm,
219 SelfConjugate -> False,
220 Ghost -> Wbar,
221 QuantumNumbers -> {GhostNumber -> 1, Q -> -1},
222 Mass -> {MW,80.385},
223 Width -> {WW,2.11},
224 PropagatorLabel -> "uWm",
225 PropagatorType -> GhostDash,
226 PropagatorArrow -> Forward },
227 U[4] == {
228 ClassName -> ghG,
229 SelfConjugate -> False,
230 Indices -> {Index[Gluon]},
231 Ghost -> G,
232 QuantumNumbers ->{GhostNumber -> 1},
233 Mass -> 0,
234 Width -> 0,
235 PropagatorLabel -> "uG",
236 PropagatorType -> GhostDash,
237 PropagatorArrow -> Forward },
238 U[5] == {
239 ClassName -> ghZp,
240 SelfConjugate -> False,
241 Ghost -> Zp,
242 QuantumNumbers -> {GhostNumber -> 1},
243 Mass -> {MZp,Internal},
244 Width -> {WZp,13.2},
245 PropagatorLabel -> "uZp",
246 PropagatorType -> GhostDash,
247 PropagatorArrow -> Forward },
248 U[61] == {
249 ClassName -> ghWpp,
250 SelfConjugate -> False,
251 Ghost -> Wp,
252 QuantumNumbers -> {GhostNumber -> 1, Q -> 1},
253 Mass -> {MWp,Internal},
254 Width -> {WWp,3.04},
255 PropagatorLabel -> "uWpp",
256 PropagatorType -> GhostDash,
257 PropagatorArrow -> Forward },
258 U[62] == {
259 ClassName -> ghWpm,
260 SelfConjugate -> False,
261 Ghost -> Wpbar,
262 QuantumNumbers -> {GhostNumber -> 1, Q -> -1},
263 Mass -> {MWp,Internal},
264 Width -> {WWp,3.04},
265 PropagatorLabel -> "uWpm",
266 PropagatorType -> GhostDash,
267 PropagatorArrow -> Forward },
268
269 (* Gauge bosons: unphysical vector fields *)
270 V[11] == {
271 ClassName -> B,
272 Unphysical -> True,
273 SelfConjugate -> True,
274 Definitions -> { B[mu_] -> TAZZp[3,1] (-sth Z[mu] + cth Zp[mu]) + TAZZp[3,2] A[mu] + TAZZp[3,3] (cth Z[mu] + sth Zp[mu])} },
275 V[12] == {
276 ClassName -> Wi,
277 Unphysical -> True,
278 SelfConjugate -> True,
279 Indices -> {Index[SU2L]},
280 FlavorIndex -> SU2L,
281 Definitions -> { Wi[mu_,1] -> (Wbar[mu]+W[mu])/Sqrt[2], Wi[mu_,2] -> (Wbar[mu]-W[mu])/(I*Sqrt[2]),
282 Wi[mu_,3] -> TAZZp[1,1] (-sth Z[mu] + cth Zp[mu]) + TAZZp[1,2] A[mu]} },
283 V[13] == {
284 ClassName -> Wpi,
285 Unphysical -> True,
286 SelfConjugate -> True,
287 Indices -> {Index[SU2R]},
288 FlavorIndex -> SU2R,
289 Definitions -> { Wpi[mu_,1] -> (Wpbar[mu]+Wp[mu])/Sqrt[2], Wpi[mu_,2] -> (Wpbar[mu]-Wp[mu])/(I*Sqrt[2]),
290 Wpi[mu_,3] -> TAZZp[2,1] (-sth Z[mu] + cth Zp[mu]) + TAZZp[2,2] A[mu] + TAZZp[2,3] (cth Z[mu] + sth Zp[mu]) } },
291 (* Ghosts: related to unphysical gauge bosons *)
292 U[11] == {
293 ClassName -> ghB,
294 Unphysical -> True,
295 SelfConjugate -> False,
296 Ghost -> B,
297 Definitions -> {ghB -> TAZZp[3,1] (-sth ghZ + cth ghZp) + TAZZp[3,2] ghA + TAZZp[3,3] (cth ghZ + sth ghZp)}},
298 U[12] == {
299 ClassName -> ghWi,
300 Unphysical -> True,
301 SelfConjugate -> False,
302 Ghost -> Wi,
303 Indices -> {Index[SU2L]},
304 FlavorIndex -> SU2L,
305 Definitions -> {ghWi[1] -> (ghWp+ghWm)/Sqrt[2], ghWi[2] -> (ghWm-ghWp)/(I*Sqrt[2]), ghWi[3] -> TAZZp[1,1] (-sth ghZ + cth ghZp) + TAZZp[1,2] ghA } },
306 U[13] == {
307 ClassName -> ghWpi,
308 Unphysical -> True,
309 SelfConjugate -> False,
310 Ghost -> {Wpi},
311 Indices -> {Index[SU2R]},
312 FlavorIndex -> SU2R,
313 Definitions -> { ghWpi[1] -> (ghWpp+ghWpm)/Sqrt[2], ghWpi[2] -> (ghWpm-ghWpp)/(I*Sqrt[2]),
314 ghWpi[3] -> TAZZp[2,1] (-sth ghZ + cth ghZp) + TAZZp[2,2] ghA + TAZZp[2,3] (cth ghZ + sth ghZp) } } ,
315 (* Fermions: physical fields *)
316 F[1] == {
317 ClassName -> v,
318 ClassMembers -> {ve,vm,vt},
319 Indices -> {Index[Generation]},
320 FlavorIndex -> Generation,
321 SelfConjugate -> False,
322 Mass -> {Mv,{Mve,1.0*^-12},{Mvm,0.0089*^-9},{Mvt,0.0504*^-9}},
323 Width -> 0,
324 QuantumNumbers -> {LeptonNumber -> 1},
325 PropagatorLabel -> {"v", "ve", "vm", "vt"} ,
326 PropagatorType -> S,
327 PropagatorArrow -> Forward,
328 PDG -> {12,14,16},
329 ParticleName -> {"ve","vm","vt"},
330 AntiParticleName -> {"ve~","vm~","vt~"},
331 FullName -> {"Electron-neutrino", "Mu-neutrino", "Tau-neutrino"} },
332 F[2] == {
333 ClassName -> l,
334 ClassMembers -> {e,mu,ta},
335 Indices -> {Index[Generation]},
336 FlavorIndex -> Generation,
337 SelfConjugate -> False,
338 Mass -> {Ml, {Me,5.11*^-4}, {MMU,0.10566}, {MTA,1.777}},
339 Width -> 0,
340 QuantumNumbers -> {Q -> -1, LeptonNumber -> 1},
341 PropagatorLabel -> {"l", "e", "mu", "ta"},
342 PropagatorType -> Straight,
343 PropagatorArrow -> Forward,
344 PDG -> {11,13,15},
345 ParticleName -> {"e-", "mu-", "ta-"},
346 AntiParticleName -> {"e+", "mu+", "ta+"},
347 FullName -> {"Electron", "Muon", "Tau"} },
348 F[3] == {
349 ClassName -> uq,
350 ClassMembers -> {u,c,t},
351 Indices -> {Index[Generation], Index[Colour]},
352 FlavorIndex -> Generation,
353 SelfConjugate -> False,
354 Mass -> {Mu, {MU, 2.55*^-3}, {MC,1.27}, {MT,172}},
355 Width -> {0,0,{WT,1.50833649}},
356 QuantumNumbers -> {Q -> 2/3},
357 PropagatorLabel -> {"uq", "u", "c", "t"},
358 PropagatorType -> Straight,
359 PropagatorArrow -> Forward,
360 PDG -> {2,4,6},
361 ParticleName -> {"u", "c", "t" },
362 AntiParticleName -> {"u~", "c~", "t~"},
363 FullName -> {"u-quark", "c-quark", "t-quark"} },
364 F[4] == {
365 ClassName -> dq,
366 ClassMembers -> {d, s, b},
367 Indices -> {Index[Generation], Index[Colour]},
368 FlavorIndex -> Generation,
369 SelfConjugate -> False,
370 Mass -> {Md, {MD,5.04*^-3}, {MS,0.101}, {MB,4.7}},
371 Width -> 0,
372 QuantumNumbers -> {Q -> -1/3},
373 PropagatorLabel -> {"dq", "d", "s", "b"},
374 PropagatorType -> Straight,
375 PropagatorArrow -> Forward,
376 PDG -> {1,3,5},
377 ParticleName -> {"d", "s", "b" },
378 AntiParticleName -> {"d~", "s~", "b~"},
379 FullName -> {"d-quark", "s-quark", "b-quark"} },
380 F[5] == {
381 ClassName -> n,
382 ClassMembers -> {ne,nm,nt},
383 Indices -> {Index[Generation]},
384 FlavorIndex -> Generation,
385 SelfConjugate -> False,
386 Mass -> {Mn,{Mne,300},{Mnm,500},{Mnt,700}},
387 Width -> {0,{Wnm,7.09*^-4},{Wnt,1.13}},
388 QuantumNumbers -> {Q->0,LeptonNumber -> 1},
389 PropagatorLabel -> {"n", "ne", "nm", "nt"} ,
390 PropagatorType -> S,
391 PropagatorArrow -> Forward,
392 PDG -> {103,104,105},
393 ParticleName -> {"ne","nm","nt"},
394 AntiParticleName -> {"ne~","nm~","nt~"},
395 FullName -> {"Electron-scotino", "Mu-scotino", "Tau-scotino"} },
396 F[6] == {
397 ClassName -> ddq,
398 ClassMembers -> {dd, ds, db},
399 Indices -> {Index[Generation], Index[Colour]},
400 FlavorIndex -> Generation,
401 SelfConjugate -> False,
402 Mass -> {Mdd, {MDD,300}, {MDS,600}, {MDB,1000}},
403 Width -> {{Wdd,0.}, {Wds,0.262}, {Wdb,9.85}},
404 QuantumNumbers -> {Q -> -1/3},
405 PropagatorLabel -> {"ddq", "dd", "ds", "db"},
406 PropagatorType -> Straight,
407 PropagatorArrow -> Forward,
408 PDG -> {106,107,108},
409 ParticleName -> {"dd", "ds", "db" },
410 AntiParticleName -> {"dd~", "ds~", "db~"},
411 FullName -> {"dd-exquark", "ds-exquark", "db-exquark"} },
412 (* Fermions: unphysical fields *)
413 F[11] == {
414 ClassName -> LL,
415 Unphysical -> True,
416 Indices -> {Index[LSU2D], Index[Generation]},
417 FlavorIndex -> LSU2D,
418 SelfConjugate -> False,
419 QuantumNumbers -> {BL -> -1/2},
420 Definitions -> { LL[sp1_,1,ff_] :> Module[{sp2,ff2}, MNSL[ff, ff2] ProjM[sp1,sp2] v[sp2,ff2]],
421 LL[sp1_,2,ff_] :> Module[{sp2}, ProjM[sp1,sp2] l[sp2,ff]] } },
422 F[12] == {
423 ClassName -> LR,
424 Unphysical -> True,
425 Indices -> {Index[RSU2D], Index[Generation]},
426 FlavorIndex -> RSU2D,
427 SelfConjugate -> False,
428 QuantumNumbers -> {BL -> -1/2},
429 Definitions -> { LR[sp1_,1,ff_] :> Module[{sp2,ff2}, MNSR[ff, ff2] ProjP[sp1,sp2] n[sp2,ff2]],
430 LR[sp1_,2,ff_] :> Module[{sp2}, ProjP[sp1,sp2] l[sp2,ff]] } },
431 F[13] == {
432 ClassName -> nuR,
433 Unphysical -> True,
434 Indices -> {Index[Generation]},
435 FlavorIndex -> Generation,
436 SelfConjugate -> False,
437 QuantumNumbers -> {BL -> 0},
438 Definitions -> { nuR[sp1_,ff_] :> Module[{sp2}, ProjP[sp1,sp2] v[sp2,ff]] } },
439 F[14] == {
440 ClassName -> nL,
441 Unphysical -> True,
442 Indices -> {Index[Generation]},
443 FlavorIndex -> Generation,
444 SelfConjugate -> False,
445 QuantumNumbers -> {BL -> 0},
446 Definitions -> { nL[sp1_,ff_] :> Module[{sp2}, ProjM[sp1,sp2] n[sp2,ff]] } },
447 F[15] == {
448 ClassName -> QL,
449 Unphysical -> True,
450 Indices -> {Index[LSU2D], Index[Generation], Index[Colour]},
451 FlavorIndex -> LSU2D,
452 SelfConjugate -> False,
453 QuantumNumbers -> {BL -> 1/6},
454 Definitions -> { QL[sp1_,1,ff_,cc_] :> Module[{sp2}, ProjM[sp1,sp2] uq[sp2,ff,cc]],
455 QL[sp1_,2,ff_,cc_] :> Module[{sp2,ff2}, CKML[ff,ff2] ProjM[sp1,sp2] dq[sp2,ff2,cc]] } },
456 F[16] == {
457 ClassName -> QR,
458 Unphysical -> True,
459 Indices -> {Index[RSU2D], Index[Generation], Index[Colour]},
460 FlavorIndex -> RSU2D,
461 SelfConjugate -> False,
462 QuantumNumbers -> {BL -> 1/6},
463 Definitions -> { QR[sp1_,1,ff_,cc_] :> Module[{sp2}, ProjP[sp1,sp2] uq[sp2,ff,cc]],
464 QR[sp1_,2,ff_,cc_] :> Module[{sp2,ff2}, CKMR[ff,ff2] ProjP[sp1,sp2] ddq[sp2,ff2,cc]] } },
465 F[17] == {
466 ClassName -> dR,
467 Unphysical -> True,
468 Indices -> {Index[Generation], Index[Colour]},
469 FlavorIndex -> Generation,
470 SelfConjugate -> False,
471 QuantumNumbers -> {BL -> -1/3},
472 Definitions -> { dR[sp1_,ff_,cc_] :> Module[{sp2}, ProjP[sp1,sp2] dq[sp2,ff,cc]] } },
473 F[18] == {
474 ClassName -> ddL,
475 Unphysical -> True,
476 Indices -> {Index[Generation], Index[Colour]},
477 FlavorIndex -> Generation,
478 SelfConjugate -> False,
479 QuantumNumbers -> {BL -> -1/3},
480 Definitions -> { ddL[sp1_,ff_,cc_] :> Module[{sp2}, ProjM[sp1,sp2] ddq[sp2,ff,cc]] } },
481
482 (* Higgs: physical scalars *)
483 S[1] == {
484 ClassName -> H0,
485 SelfConjugate -> True,
486 Mass -> {mh0,125.},
487 Width -> {WH0,0.00407},
488 PropagatorLabel -> "H",
489 PropagatorType -> D,
490 PropagatorArrow -> None,
491 PDG -> 25,
492 ParticleName -> "H0",
493 FullName -> "H0" },
494 S[2] == {
495 ClassName -> H1,
496 SelfConjugate -> True,
497 Mass -> {mh1,Internal},
498 Width -> {WH1,0.95},
499 PropagatorLabel -> "H1",
500 PropagatorType -> D,
501 PropagatorArrow -> None,
502 PDG -> 26,
503 ParticleName -> "H1",
504 FullName -> "H1" },
505 S[3] == {
506 ClassName -> H2,
507 SelfConjugate -> True,
508 Mass -> {mh2,Internal},
509 Width -> {WH2,1.46},
510 PropagatorLabel -> "H2",
511 PropagatorType -> D,
512 PropagatorArrow -> None,
513 PDG -> 27,
514 ParticleName -> "H2",
515 FullName -> "H2" },
516 S[4] == {
517 ClassName -> H3,
518 SelfConjugate -> True,
519 Mass -> {mh3,Internal},
520 Width -> {WH3,321.0},
521 PropagatorLabel -> "H3",
522 PropagatorType -> D,
523 PropagatorArrow -> None,
524 PDG -> 28,
525 ParticleName -> "H3",
526 FullName -> "H3" },
527 S[5] == {
528 ClassName -> A1,
529 SelfConjugate -> True,
530 Mass -> {ma1,Internal},
531 Width -> {WA1,0.95},
532 PropagatorLabel -> "A1",
533 PropagatorType -> D,
534 PropagatorArrow -> None,
535 PDG -> 29,
536 ParticleName -> "A1",
537 FullName -> "A1" },
538 S[6] == {
539 ClassName -> A2,
540 SelfConjugate -> True,
541 Mass -> {ma2,Internal},
542 Width -> {WA2,1.45},
543 PropagatorLabel -> "A2",
544 PropagatorType -> D,
545 PropagatorArrow -> None,
546 PDG -> 30,
547 ParticleName -> "A2",
548 FullName -> "A2" },
549 S[7] == {
550 ClassName -> HP1,
551 SelfConjugate -> False,
552 Mass -> {mhp1,Internal},
553 QuantumNumbers -> {Q -> 1},
554 Width -> {WHP1,1.00},
555 PropagatorLabel -> "HP1",
556 PropagatorType -> D,
557 PropagatorArrow -> None,
558 PDG -> 31,
559 ParticleName -> "H1+",
560 AntiParticleName -> "H1-",
561 FullName -> "HP1" },
562 S[8] == {
563 ClassName -> HM2,
564 SelfConjugate -> False,
565 Mass -> {mhm2,Internal},
566 QuantumNumbers -> {Q -> -1},
567 Width -> {WHM2,0.325},
568 PropagatorLabel -> "HM2",
569 PropagatorType -> D,
570 PropagatorArrow -> None,
571 PDG -> 32,
572 ParticleName -> "H2-",
573 AntiParticleName -> "H2+",
574 FullName -> "HM2" },
575 (* Higgs: Goldstone *)
576 S[9] == {
577 ClassName -> G10,
578 SelfConjugate -> True,
579 Goldstone -> Z,
580 Mass -> {MZ,91.1876},
581 Width -> {WZ,2.52},
582 PropagatorLabel -> "G1o",
583 PropagatorType -> D,
584 PropagatorArrow -> None,
585 PDG -> 33,
586 ParticleName -> "G10",
587 FullName -> "G10" },
588 S[10] == {
589 ClassName -> G20,
590 SelfConjugate -> True,
591 Goldstone -> Zp,
592 Mass -> {MZp,Internal},
593 Width -> {WZp,16.7},
594 PropagatorLabel -> "G2o",
595 PropagatorType -> D,
596 PropagatorArrow -> None,
597 PDG -> 34,
598 ParticleName -> "G20",
599 FullName -> "G20" },
600 S[11] == {
601 ClassName -> GP1,
602 SelfConjugate -> False,
603 Goldstone -> W,
604 Mass -> {MW,80.385},
605 Width -> {WW,2.11},
606 QuantumNumbers -> {Q -> 1},
607 PropagatorLabel -> "GP1",
608 PropagatorType -> D,
609 PropagatorArrow -> None,
610 PDG -> 35,
611 ParticleName -> "G1+",
612 AntiParticleName -> "G1-",
613 FullName -> "GP1" },
614 S[12] == {
615 ClassName -> GM2,
616 SelfConjugate -> False,
617 Goldstone -> Wp,
618 Mass -> {MWp,Internal},
619 Width -> {WWp,6.33},
620 QuantumNumbers -> {Q -> -1},
621 PropagatorLabel -> "GM2",
622 PropagatorType -> D,
623 PropagatorArrow -> None,
624 PDG -> 36,
625 ParticleName -> "G2-",
626 AntiParticleName -> "G2+",
627 FullName -> "GM2" },
628
629 (* Higgs: unphysical scalars *)
630 S[131] == {
631 ClassName -> Phi,
632 Unphysical -> True,
633 Indices -> {Index[LSU2D],Index[RSU2D]},
634 SelfConjugate -> False,
635 QuantumNumbers -> {BL -> 0},
636 Definitions -> { Phi[1,1] -> (H1 +I A1)/Sqrt[2], Phi[1,2] -> (cb HP1 +sb GP1), Phi[2,1] -> cz HM2 + sz GM2,
637 Phi[2,2] -> (k+THP0 H0+THP2 H2+THP3 H3+I((1/Sqrt[1+tb^2+tz^2])A2-(tz/Sqrt[1+tz^2])G10+(tb/Sqrt[(1+tz^2)(1+tb^2+tz^2)])G20))/Sqrt[2] } },
638 S[141] == {
639 ClassName -> ChiL,
640 Unphysical -> True,
641 Indices -> {Index[LSU2D]},
642 SelfConjugate -> False,
643 QuantumNumbers -> {BL -> 1/2},
644 Definitions -> { ChiL[1] -> (-sb HP1 + cb GP1) ,
645 ChiL[2] -> (vL+THL0 H0+THL2 H2+THL3 H3 + I ((-tb/Sqrt[1+tb^2+tz^2]) A2 +(Sqrt[(1+tz^2)/(1+tb^2+tz^2)]) G20))/Sqrt[2]} },
646 S[151] == {
647 ClassName -> ChiR,
648 Unphysical -> True,
649 Indices -> {Index[RSU2D]},
650 SelfConjugate -> False,
651 QuantumNumbers -> {BL -> 1/2},
652 Definitions -> { ChiR[1] -> -sz HM2bar + cz GM2bar,
653 ChiR[2] -> (vR+THR0 H0+THR2 H2+THR3 H3+I ((tz/Sqrt[1+tb^2+tz^2])A2 +(1/Sqrt[1+tz^2])G10+(tb tz/Sqrt[(1+tz^2)(1+tb^2+tz^2)])G20))/Sqrt[2]} } };
654
655 (* ************************** *)
656 (* ***** Gauge ***** *)
657 (* ***** Parameters ***** *)
658 (* ***** (FeynArts) ***** *)
659 (* ************************** *)
660
661GaugeXi[ V[1] ] = GaugeXi[A];
662GaugeXi[ V[2] ] = GaugeXi[Z];
663GaugeXi[ V[3] ] = GaugeXi[W];
664GaugeXi[ V[4] ] = GaugeXi[G];
665GaugeXi[ V[5] ] = GaugeXi[Zp];
666GaugeXi[ V[6] ] = GaugeXi[Wp];
667GaugeXi[ S[1] ] = 1;
668GaugeXi[ S[2] ] = 1;
669GaugeXi[ S[3] ] = 1;
670GaugeXi[ S[4] ] = 1;
671GaugeXi[ S[5] ] = 1;
672GaugeXi[ S[6] ] = 1;
673GaugeXi[ S[7] ] = 1;
674GaugeXi[ S[8] ] = 1;
675GaugeXi[ S[9] ] = GaugeXi[Z];
676GaugeXi[ S[11] ] = GaugeXi[W];
677GaugeXi[ S[10] ] = GaugeXi[Zp];
678GaugeXi[ S[12] ] = GaugeXi[Wp];
679GaugeXi[ U[1] ] = GaugeXi[A];
680GaugeXi[ U[2] ] = GaugeXi[Z];
681GaugeXi[ U[31] ] = GaugeXi[W];
682GaugeXi[ U[32] ] = GaugeXi[W];
683GaugeXi[ U[4] ] = GaugeXi[G];
684GaugeXi[ U[5] ] = GaugeXi[Zp];
685GaugeXi[ U[61] ] = GaugeXi[Wp];
686GaugeXi[ U[62] ] = GaugeXi[Wp];
687
688 (* ************************** *)
689 (* ***** Lagrangian ***** *)
690 (* ************************** *)
691
692 (* ***** Gauge Lagrangian ***** *)
693
694LGauge := Block[{mu,nu,ii,jj,aa},
695 ExpandIndices[-1/4 FS[B,mu,nu] FS[B,mu,nu]-1/4 FS[Wi,mu,nu,ii] FS[Wi,mu,nu,ii]
696 -1/4 FS[Wpi,mu,nu,jj] FS[Wpi,mu,nu,jj]-1/4 FS[G,mu,nu,aa] FS[G,mu,nu,aa], FlavorExpand->{SU2L,SU2R}]];
697
698 (* ***** Fermion Lagrangian ***** *)
699
700LFermions := Block[{mu},
701 ExpandIndices[I*(QLbar.Ga[mu].DC[QL,mu]+QRbar.Ga[mu].DC[QR,mu]+LLbar.Ga[mu].DC[LL,mu]+LRbar.Ga[mu].DC[LR,mu]
702 +ddLbar.Ga[mu].DC[ddL,mu]+dRbar.Ga[mu].DC[dR, mu]+nLbar.Ga[mu].DC[nL,mu]+nuRbar.Ga[mu].DC[nuR,mu]),
703 FlavorExpand->{SU2L,SU2R,LSU2D,RSU2D}]/.{CKML[a_,b_] Conjugate[CKML[a_,c_]]->IndexDelta[b,c],CKML[b_,a_] Conjugate[CKML[c_,a_]]-> IndexDelta[b,c], CKMR[a_,b_] Conjugate[CKMR[a_,c_]]->IndexDelta[b,c],CKMR[b_,a_] Conjugate[CKMR[c_,a_]]-> IndexDelta[b,c], MNSL[a_,b_] Conjugate[MNSL[a_,c_]]->IndexDelta[b,c], MNSL[b_,a_] Conjugate[MNSL[c_,a_]]->IndexDelta[b,c], MNSR[a_,b_] Conjugate[MNSR[a_,c_]]->IndexDelta[b,c], MNSR[b_,a_] Conjugate[MNSR[c_,a_]]->IndexDelta[b,c]}];
704
705 (* ***** Higgs Lagrangian ***** *)
706
707(* Higgs Kinetic *)
708
709LHiggskin := Block[{ii,jj,kk,ll,aa,bb,mu, feynmangaugerules},
710 feynmangaugerules=If[Not[FeynmanGauge],{G10|G20|GP1|GM2|GP1bar|GM2bar->0},{}];
711 ExpandIndices[
712Sum[(del[Phibar[ii,jj],mu]+(I/2)(gL Wi[mu,bb] HC[PauliSigma[bb,ii,ll]] Phibar[ll,jj] - gR Phibar[ii,ll] Wpi[mu,bb] HC[PauliSigma[bb,ll,jj]])) (del[Phi[ii,jj], mu]-(I/2) (gL Wi[mu,aa] PauliSigma[aa,ii,kk] Phi[kk,jj] - gR Phi[ii,kk] Wpi[mu,aa] PauliSigma[aa,kk,jj])),{ii,2},{jj,2}] +DC[ChiLbar[ii],mu] DC[ChiL[ii],mu] + DC[ChiRbar[jj],mu] DC[ChiR[jj],mu],FlavorExpand->{LSU2D,RSU2D,SU2L,SU2R}]/.feynmangaugerules];
713
714(* Higgs Potential *)
715
716LHiggspot := Block[{ii,jj,kk,kk1,kk2,kk3,kk4,ll,feynmangaugerules},
717 feynmangaugerules = If[Not[FeynmanGauge], {G10|G20|GP1|GM2|GP1bar|GM2bar ->0}, {}];
718-ExpandIndices[
719-mu1 (Phibar[ii,jj]Phi[ii,jj]) + lm1 (Phibar[ii,jj]Phi[ii,jj])(Phibar[kk,ll]Phi[kk,ll]) + lm2 (Phibar[ii,jj] (-Eps[ii,kk1]HC[Phi[kk1,kk2]]Eps[kk2,jj])) ((-Eps[kk,kk1]Phi[kk1,kk2]Eps[kk2,ll]) Phi[kk,ll]) - mu2 (ChiLbar[ii] ChiL[ii] + ChiRbar[ii] ChiR[ii]) + lm3 ((ChiLbar[ii] ChiL[ii])(ChiLbar[jj] ChiL[jj]) + (ChiRbar[kk] ChiR[kk])(ChiRbar[ll] ChiR[ll])) + 2lm4 (ChiLbar[ii] ChiL[ii]) (ChiRbar[jj] ChiR[jj]) + 2al1 (Phibar[ii,jj]Phi[ii,jj]) (ChiLbar[kk] ChiL[kk] + ChiRbar[kk] ChiR[kk]) + 2al2 (ChiLbar[ii] Phi[ii,jj] Phibar[kk,jj] ChiL[kk] + ChiRbar[ii] Phibar[jj,ii] Phi[jj,kk] ChiR[kk]) + 2al3 Sum[(ChiLbar[ii] (-Eps[ii,kk1]HC[Phi[kk1,kk2]]Eps[kk2,jj]) (HC[PauliMatrix[2].Table[HC[Phi[kk3,kk4]],{kk3,2},{kk4,2}].PauliMatrix[2]][[jj,kk]]) ChiL[kk] + ChiRbar[ii] (HC[PauliMatrix[2].Table[HC[Phi[kk3,kk4]],{kk3,2},{kk4,2}].PauliMatrix[2]][[ii,jj]]) (-Eps[jj,kk1]HC[Phi[kk1,kk2]]Eps[kk2,kk]) ChiR[kk]),{ii,2},{jj,2},{kk,2},{kk1,2},{kk2,2}] + mu3 (ChiLbar[ii] Phi[ii,jj] ChiR[jj] + ChiRbar[ii] Phibar[jj,ii] ChiL[jj]),FlavorExpand->{LSU2D,RSU2D}]/.feynmangaugerules ];
720
721LHiggs := LHiggskin + LHiggspot;
722
723 (* ***** Yukawa Lagrangian ***** *)
724
725LYukawa := Block[{sp,ii,jj,cc,ff1,ff2,ff3,yuk,feynmangaugerules},
726 feynmangaugerules = If[Not[FeynmanGauge], {G10|G20|GP1|GM2|GP1bar|GM2bar ->0}, {}];
727 yuk = ExpandIndices[
728 (* Quarks *) -yu[ff1,ff2] QLbar[sp,ii,ff1,cc].QR[sp,jj,ff2,cc] (-Eps[ii,kk1]HC[Phi[kk1,kk2]]Eps[kk2,jj])
729 -yd[ff2,ff3] CKML[ff1,ff2] QLbar[sp,ii,ff1,cc].dR[sp,ff3,cc] ChiL[ii]
730 -ydd[ff2,ff3] CKMR[ff1,ff2] QRbar[sp,ii,ff1,cc].ddL[sp,ff3,cc] ChiR[ii]
731 (* Leptons *)-yl[ff1,ff2] LLbar[sp,ii,ff1].LR[sp,jj,ff2] Phi[ii,jj]
732 -ynu[ff2,ff3] MNSL[ff1,ff2] LLbar[sp,ii,ff1].nuR[sp,ff3] ChiLbar[jj] Eps[ii, jj]
733 -yn[ff2,ff3] MNSR[ff1,ff2] LRbar[sp,ii,ff1].nL[sp,ff3] ChiRbar[jj] Eps[ii, jj]
734 (* seesaw *)(*+ MN nuRbar nuR*), FlavorExpand->{LSU2D,RSU2D}];
735 yuk=yuk/.{CKML[a_,b_] Conjugate[CKML[a_,c_]]->IndexDelta[b,c], CKML[b_,a_] Conjugate[CKML[c_,a_]]->IndexDelta[b,c], CKMR[a_,b_] Conjugate[CKMR[a_,c_]]->IndexDelta[b,c], CKMR[b_,a_] Conjugate[CKMR[c_,a_]]->IndexDelta[b,c],
736 MNSL[a_,b_]Conjugate[MNSL[a_,c_]]->IndexDelta[b,c],MNSL[b_,a_]Conjugate[MNSL[c_,a_]]->IndexDelta[b,c],
737 MNSR[a_,b_]Conjugate[MNSR[a_,c_]]->IndexDelta[b,c],MNSR[b_,a_]Conjugate[MNSR[c_,a_]]->IndexDelta[b,c]};yuk+HC[yuk]/.feynmangaugerules];
738
739 (* ***** Ghost Lagrangian ***** *)
740
741LGhost := Block[{LGhs,LGhB,LGhw,LGhwp,LGhphi,mu,generators,gh,ghbar,Vectorize,Vectorize1,phi11,phi12,phi21,phi22,phi31,phi32,phi41,phi42,LGhChiL,LGhChiR,LGhPhiLL,LGhPhiRR,LGhPhiLR,LGhPhiRL,togoldstones,bidoublet,bidoublet0,doubletchiL,doubletchiL0,doubletchiR,doubletchiR0},
742 (* Pure gauge piece *)
743 LGhs = -ghGbar.del[DC[ghG,mu],mu];
744 LGhB = -ghBbar.del[DC[ghB,mu],mu];
745 LGhw = -ghWibar[x1].del[del[ghWi[x1],mu],mu]- gL ghWibar[jjj].(del[Wi[mu,i1],mu] Eps[jjj,i1,i2].ghWi[i2]+Eps[jjj,i1,i2].del[ghWi[i2],mu] Wi[mu,i1]);
746 LGhwp = -ghWpibar[x2].del[del[ghWpi[x2],mu],mu]- gL ghWpibar[jjj].(del[Wpi[mu,i1],mu] Eps[jjj,i1,i2].ghWpi[i2]+ Eps[jjj,i1,i2].del[ghWpi[i2],mu] Wpi[mu,i1]);
747
748 (* Scalar pieces: see Peskin pages 739-742 *)
749 (* phi11, phi12 are the real degrees of freedom of (Phi) HP1, GP1 and phi21, phi22 are the real degrees of freedom of (Phi) HM2, GM2 *)
750 (* phi31, phi32 are the real degrees of freedom of (ChiL) HP1, GP1 and phi41, phi42 are the real degrees of freedom of (ChiR) HM2, GM2 *)
751 (* Vectorize transforms a doublet in a vector in the phi-basis, i.e. the basis of real degrees of freedom *)
752 (* Vectorize1 transforms the bidoublet in a vector in the phi-basis, i.e. the basis of real degrees of freedom *)
753
754 (* Higgs-ghost piece *)
755 gh = {ghB, ghWi[1], ghWi[2], ghWi[3], ghWpi[1], ghWpi[2], ghWpi[3]};
756 ghbar = {ghBbar, ghWibar[1], ghWibar[2], ghWibar[3], ghWpibar[1], ghWpibar[2], ghWpibar[3]};
757 generators = {-I/2 gBL IdentityMatrix[2], -I/2 gL PauliSigma[1], -I/2 gL PauliSigma[2], -I/2 gL PauliSigma[3], -I/2 gR PauliSigma[1], -I/2 gR PauliSigma[2], -I/2 gR PauliSigma[3]};
758 bidoublet = Expand[{{ Phi[1,1] , (phi11+I phi12)/Sqrt[2] },{ (phi21+I phi22)/Sqrt[2] , Phi[2,2]}} /. MR$Definitions /. k -> 0];
759 bidoublet0 = {{0, 0},{0, k/Sqrt[2]}};
760 doubletchiL = Expand[{ (phi31+I phi32)/Sqrt[2] , ChiL[2]} /. MR$Definitions /. vL -> 0];
761 doubletchiL0 = {0, vL/Sqrt[2]};
762 doubletchiR = Expand[{ (phi41+I phi42)/Sqrt[2] , ChiR[2]} /. MR$Definitions /. vR -> 0];
763 doubletchiR0 = {0, vR/Sqrt[2]};
764 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}];
765 Vectorize1[{{a_, b_},{c_, d_}}]:= Simplify[{Sqrt[2] Re[Expand[a]], Sqrt[2] Im[Expand[a]], Sqrt[2] Re[Expand[b]], Sqrt[2] Im[Expand[b]],Sqrt[2] Re[Expand[c]],Sqrt[2] Im[Expand[c]], Sqrt[2] Re[Expand[d]], Sqrt[2] Im[Expand[d]]}/.{Im[_]->0, Re[num_]->num}];
766 togoldstones := {phi11 -> cb (HP1 + HP1bar)/Sqrt[2] + sb (GP1 + GP1bar)/Sqrt[2], phi12 -> cb (HP1 - HP1bar)/(I Sqrt[2]) + sb (GP1 - GP1bar)/(I Sqrt[2]), phi21 -> cz (HM2 + HM2bar)/Sqrt[2] + sz (GM2 + GM2bar)/Sqrt[2], phi22 -> cz (HM2 - HM2bar)/(I Sqrt[2]) + sz (GM2 - GM2bar)/(I Sqrt[2]),phi31 -> -sb (HP1 + HP1bar)/Sqrt[2] + cb (GP1 + GP1bar)/Sqrt[2], phi32 -> -sb (HP1 - HP1bar)/(I Sqrt[2]) + cb (GP1 - GP1bar)/(I Sqrt[2]), phi41 -> -sz (HM2 + HM2bar)/Sqrt[2] + cz (GM2 + GM2bar)/Sqrt[2], phi42 -> -sz (-HM2 + HM2bar)/(I Sqrt[2]) + cz (-GM2 + GM2bar)/(I Sqrt[2])};
767 LGhChiL =Plus@@Flatten[Table[-ghbar[[kkk]].gh[[lll]] Vectorize[generators[[kkk]].doubletchiL0].Vectorize[generators[[lll]].(doubletchiL0+doubletchiL)],{kkk,4},{lll,4}]] /.togoldstones;
768 LGhChiR =Plus@@Flatten[Table[-ghbar[[kkk]].gh[[lll]] Vectorize[generators[[kkk]].doubletchiR0].Vectorize[generators[[lll]].(doubletchiR0+doubletchiR)],{kkk,{1,5,6,7}},{lll,{1,5,6,7}}]] /.togoldstones;
769 LGhPhiLL =Plus@@Flatten[Table[-ghbar[[kkk]].gh[[lll]] Vectorize1[generators[[kkk]].bidoublet0].Vectorize1[generators[[lll]].(bidoublet0+bidoublet)],{kkk,2,4},{lll,2,4}]] /.togoldstones;
770 LGhPhiRR =Plus@@Flatten[Table[-ghbar[[kkk]].gh[[lll]] Vectorize1[bidoublet0.generators[[kkk]]].Vectorize1[(bidoublet0+bidoublet).generators[[lll]]],{kkk,5,7},{lll,5,7}]] /.togoldstones;
771 LGhPhiLR =Plus@@Flatten[Table[+ghbar[[kkk]].gh[[lll]] Vectorize1[generators[[kkk]].bidoublet0].Vectorize1[(bidoublet0+bidoublet).generators[[lll]]],{kkk,2,4},{lll,5,7}]] /.togoldstones;
772 LGhPhiRL =Plus@@Flatten[Table[+ghbar[[kkk]].gh[[lll]] Vectorize1[bidoublet0.generators[[kkk]]].Vectorize1[generators[[lll]].(bidoublet0+bidoublet)],{kkk,5,7},{lll,2,4}]] /.togoldstones;
773 LGhPhi =LGhPhiLL+LGhPhiRR+LGhPhiLR+LGhPhiRL;
774 LGhphi =LGhChiL+LGhChiR+LGhPhi;
775ExpandIndices[ LGhs + If[FeynmanGauge, LGhB + LGhw + LGhwp + LGhphi,0], FlavorExpand->{SU2L,SU2R}]];
776
777 (* ***** Effective Higgs-gamma-gamma and Higgs-gluon-gluon Lagrangian ***** *)
778
779LEff := Block[{mu,nu,aa}, ExpandIndices[ H0 ( Ghaa FS[A,mu,nu] FS[A,mu,nu] + Ghgg FS[G,mu,nu,aa] FS[G,mu,nu,aa] ) ]];
780
781 (* ***** Full Lagrangian ***** *)
782
783LALRM:= LGauge + LFermions + LHiggs + LYukawa + LGhost + LEff;
784
785(* ************************** *)
786(* ***** Parameters ***** *)
787(* ************************** *)
788
789M$Parameters = {
790
791 (* External Parameters *)
792
793 aEWM1 == {
794 ParameterType -> External,
795 BlockName -> SMINPUTS,
796 OrderBlock -> 1,
797 Value -> 132.50698,
798 InteractionOrder -> {QED,-2},
799 TeX -> Superscript[Subscript[\[Alpha],EW],-1],
800 Description -> "Inverse of the EW coupling constant at the Z pole" },
801 Gf == {
802 ParameterType -> External,
803 BlockName -> SMINPUTS,
804 OrderBlock -> 2,
805 Value -> 1.16637*^-5,
806 InteractionOrder -> {QED,2},
807 TeX -> Subscript[G,F],
808 Description -> "Fermi constant" },
809 aS == {
810 ParameterType -> External,
811 BlockName -> SMINPUTS,
812 OrderBlock -> 3,
813 Value -> 0.1184,
814 InteractionOrder -> {QCD,2},
815 TeX -> Subscript[\[Alpha],s],
816 Description -> "Strong coupling constant at the Z pole" },
817 sw2 == {
818 ParameterType -> External,
819 BlockName -> SMINPUTS,
820 OrderBlock -> 4,
821 Value -> 0.22343,
822 InteractionOrder -> {QED,0},
823 TeX -> Superscript[Subscript[s,w],2],
824 Description -> "Squared Sine of the Weinberg angle" },
825 lm2 == {
826 ParameterType -> External,
827 BlockName -> ALRMINPUTS,
828 OrderBlock -> 1,
829 Value -> 0.05,
830 InteractionOrder -> {QED,2},
831 TeX -> Subscript[\[Lambda],2],
832 Description -> "Higgs quartic coupling" },
833 lm3 == {
834 ParameterType -> External,
835 BlockName -> ALRMINPUTS,
836 OrderBlock -> 2,
837 Value -> 0.006,
838 InteractionOrder -> {QED,2},
839 TeX -> Subscript[\[Lambda],3],
840 Description -> "Higgs quartic coupling" },
841 al1 == {
842 ParameterType -> External,
843 BlockName -> ALRMINPUTS,
844 OrderBlock -> 3,
845 Value -> 0.003,
846 InteractionOrder -> {QED,2},
847 TeX -> Subscript[\[Alpha],1],
848 Description -> "Higgs quartic coupling" },
849 al2 == {
850 ParameterType -> External,
851 BlockName -> ALRMINPUTS,
852 OrderBlock -> 4,
853 Value -> -0.01,
854 InteractionOrder -> {QED,2},
855 TeX -> Subscript[\[Alpha],2],
856 Description -> "Higgs quartic coupling" },
857 al3 == {
858 ParameterType -> External,
859 BlockName -> ALRMINPUTS,
860 OrderBlock -> 5,
861 Value -> 0.005,
862 InteractionOrder -> {QED,2},
863 TeX -> Subscript[\[Alpha],3],
864 Description -> "Higgs quartic coupling" },
865 tb == {
866 ParameterType -> External,
867 BlockName -> ALRMINPUTS,
868 OrderBlock -> 6,
869 Value -> 10.,
870 InteractionOrder -> {QED,0},
871 TeX -> Subscript[t,\[Beta]],
872 Description -> "vevs ratio: Tan beta" },
873 mu3 == {
874 ParameterType -> External,
875 BlockName -> ALRMINPUTS,
876 OrderBlock -> 7,
877 Value -> -0.22, (*Most relevant to mh2*)
878 InteractionOrder -> {QED,1},
879 TeX -> Subscript[\[Mu],3],
880 Description -> "Higgs trilinear coupling" },
881 gR == {
882 ParameterType -> External,
883 BlockName -> ALRMINPUTS,
884 OrderBlock -> 8,
885 Value -> 0.642,
886 InteractionOrder -> {QED,1},
887 TeX -> Subscript[g,R],
888 Description -> "Right Weak coupling constant at the Z pole" },
889 ymdo == {
890 ParameterType -> External,
891 BlockName -> YUKAWA,
892 OrderBlock -> 1,
893 Value -> 5.04*^-3,
894 TeX -> Subscript[m,d],
895 Description -> "Down Yukawa mass" },
896 ymup == {
897 ParameterType -> External,
898 BlockName -> YUKAWA,
899 OrderBlock -> 2,
900 Value -> 2.55*^-3,
901 TeX -> Subscript[m,u],
902 Description -> "Up Yukawa mass" },
903 yms == {
904 ParameterType -> External,
905 BlockName -> YUKAWA,
906 OrderBlock -> 3,
907 Value -> 0.101,
908 TeX -> Subscript[m,s],
909 Description -> "Strange Yukawa mass" },
910 ymc == {
911 ParameterType -> External,
912 BlockName -> YUKAWA,
913 OrderBlock -> 4,
914 Value -> 1.27,
915 TeX -> Subscript[m,c],
916 Description -> "Charm Yukawa mass" },
917 ymb == {
918 ParameterType -> External,
919 BlockName -> YUKAWA,
920 OrderBlock -> 5,
921 Value -> 4.7,
922 TeX -> Subscript[m,b],
923 Description -> "Bottom Yukawa mass" },
924 ymt == {
925 ParameterType -> External,
926 BlockName -> YUKAWA,
927 OrderBlock -> 6,
928 Value -> 172,
929 TeX -> Subscript[m,t],
930 Description -> "Top Yukawa mass" },
931 ymddo == {
932 ParameterType -> External,
933 BlockName -> YUKAWA,
934 OrderBlock -> 7,
935 Value -> 3*^2,
936 TeX -> Subscript[m,d'],
937 Description -> "Exotic-Down Yukawa mass" },
938 ymds == {
939 ParameterType -> External,
940 BlockName -> YUKAWA,
941 OrderBlock -> 8,
942 Value -> 6*^2,
943 TeX -> Subscript[m,s'],
944 Description -> "Exotic-Strange Yukawa mass" },
945 ymdb == {
946 ParameterType -> External,
947 BlockName -> YUKAWA,
948 OrderBlock -> 9,
949 Value -> 1*^3,
950 TeX -> Subscript[m,b'],
951 Description -> "Exotic-Bottom Yukawa mass" },
952 yme == {
953 ParameterType -> External,
954 BlockName -> YUKAWA,
955 OrderBlock -> 11,
956 Value -> 5.11*^-4,
957 TeX -> Subscript[m,e],
958 Description -> "Electron Yukawa mass" },
959 ymnue == {
960 ParameterType -> External,
961 BlockName -> YUKAWA,
962 OrderBlock -> 12,
963 Value -> 1.*^-12,
964 TeX -> Subscript[m,Subscript[\[Nu],e]],
965 Description -> "Electron-Neutrino Yukawa mass" },
966 ymm == {
967 ParameterType -> External,
968 BlockName -> YUKAWA,
969 OrderBlock -> 13,
970 Value -> 0.10566,
971 TeX -> Subscript[m,\[Mu]],
972 Description -> "Muon Yukawa mass" },
973 ymnum == {
974 ParameterType -> External,
975 BlockName -> YUKAWA,
976 OrderBlock -> 14,
977 Value -> 0.0089*^-9,
978 TeX -> Subscript[m,Subscript[\[Nu],\[Mu]]],
979 Description -> "Muon-Neutrino Yukawa mass" },
980 ymtau == {
981 ParameterType -> External,
982 BlockName -> YUKAWA,
983 OrderBlock -> 15,
984 Value -> 1.777,
985 TeX -> Subscript[m,\[Tau]],
986 Description -> "Tau Yukawa mass" },
987 ymnutau == {
988 ParameterType -> External,
989 BlockName -> YUKAWA,
990 OrderBlock -> 16,
991 Value -> 0.0504*^-9,
992 TeX -> Subscript[m,Subscript[\[Nu],\[Tau]]],
993 Description -> "Tau-Neutrino Yukawa mass" },
994 ymne == {
995 ParameterType -> External,
996 BlockName -> YUKAWA,
997 OrderBlock -> 17,
998 Value -> 300,
999 TeX -> Subscript[m,Subscript[n,e]],
1000 Description -> "Electron-Scotino Yukawa mass" },
1001 ymnm == {
1002 ParameterType -> External,
1003 BlockName -> YUKAWA,
1004 OrderBlock -> 18,
1005 Value -> 500,
1006 TeX -> Subscript[m,Subscript[n,\[Mu]]],
1007 Description -> "Muon-Scotino Yukawa mass" },
1008 ymntau == {
1009 ParameterType -> External,
1010 BlockName -> YUKAWA,
1011 OrderBlock -> 19,
1012 Value -> 700,
1013 TeX -> Subscript[m,Subscript[n,\[Tau]]],
1014 Description -> "Tau-Scotino Yukawa mass" },
1015 sL12 == {
1016 ParameterType -> External,
1017 BlockName -> CKMBLOCK,
1018 OrderBlock -> 1,
1019 Value -> 0.221,
1020 TeX -> Subscript[sL,12],
1021 Description -> "Left Quark mixing: SinTheta12,PDG-94"},
1022 sL23 == {
1023 ParameterType -> External,
1024 BlockName -> CKMBLOCK,
1025 OrderBlock -> 2,
1026 Value -> 0.040,
1027 TeX -> Subscript[sL,23],
1028 Description -> "Left Quark mixing: SinTheta23,PDG-94"},
1029 sL13 == {
1030 ParameterType -> External,
1031 BlockName -> CKMBLOCK,
1032 OrderBlock -> 3,
1033 Value -> 0.0035,
1034 TeX -> Subscript[sL,13],
1035 Description -> "Left Quark mixing: SinTheta13,PDG-94"},
1036
1037 delL13 == {
1038 ParameterType -> External,
1039 BlockName -> CKMBLOCK,
1040 OrderBlock -> 4,
1041 Value -> 1.20,
1042 TeX -> Subscript[\[Delta]L,13],
1043 Description -> "Left CKM CP-Violation Phase"},
1044
1045 sR12 == {
1046 ParameterType -> External,
1047 BlockName -> CKMBLOCK,
1048 OrderBlock -> 5,
1049 Value -> 0.225,
1050 TeX -> Subscript[sR,12],
1051 Description -> "Right Quark mixing: SinThetaR12,PDG-194"},
1052 sR23 == {
1053 ParameterType -> External,
1054 BlockName -> CKMBLOCK,
1055 OrderBlock -> 6,
1056 Value -> 0.040,
1057 TeX -> Subscript[sR,23],
1058 Description -> "Right Quark mixing: SinThetaR23,PDG-194"},
1059 sR13 == {
1060 ParameterType -> External,
1061 BlockName -> CKMBLOCK,
1062 OrderBlock -> 7,
1063 Value -> 0.0035,
1064 TeX -> Subscript[sR,13],
1065 Description -> "Right Quark mixing: SinThetaR13,PDG-194"},
1066 delR13 == {
1067 ParameterType -> External,
1068 BlockName -> CKMBLOCK,
1069 OrderBlock -> 8,
1070 Value -> 1.20,
1071 TeX -> Subscript[\[Delta]R,13],
1072 Description -> "Right CKM CP-Violation Phase"},
1073
1074 ssL12 == {
1075 ParameterType -> External,
1076 BlockName -> MNSBLOCK,
1077 OrderBlock -> 1,
1078 Value -> 0.5506,
1079 TeX -> Subscript[ssL,12],
1080 Description -> "Left Lepton mixing: SinTheta12"},
1081 ssL23 == {
1082 ParameterType -> External,
1083 BlockName -> MNSBLOCK,
1084 OrderBlock -> 2,
1085 Value -> 0.7558,
1086 TeX -> Subscript[ssL,23],
1087 Description -> "Left Lepton mixing: SinTheta23"},
1088 ssL13 == {
1089 ParameterType -> External,
1090 BlockName -> MNSBLOCK,
1091 OrderBlock -> 3,
1092 Value -> 0.1485,
1093 TeX -> Subscript[ssL,13],
1094 Description -> "Left Lepton mixing: SinTheta13"},
1095
1096 delLCP == {
1097 ParameterType -> External,
1098 BlockName -> MNSBLOCK,
1099 OrderBlock -> 4,
1100 Value -> 3.4383,
1101 TeX -> Subscript[\[Delta]L,CP],
1102 Description -> "Left MNS CP-Violation Phase"},
1103
1104 ssR12 == {
1105 ParameterType -> External,
1106 BlockName -> MNSBLOCK,
1107 OrderBlock -> 5,
1108 Value -> 0.551,
1109 TeX -> Subscript[ssR,12],
1110 Description -> "Right Lepton mixing: SinTheta12"},
1111 ssR23 == {
1112 ParameterType -> External,
1113 BlockName -> MNSBLOCK,
1114 OrderBlock -> 6,
1115 Value -> 0.756,
1116 TeX -> Subscript[ssR,23],
1117 Description -> "Right Lepton mixing: SinTheta23"},
1118 ssR13 == {
1119 ParameterType -> External,
1120 BlockName -> MNSBLOCK,
1121 OrderBlock -> 7,
1122 Value -> 0.149,
1123 TeX -> Subscript[ssR,13],
1124 Description -> "Right Lepton mixing: SinTheta13"},
1125
1126 delRCP == {
1127 ParameterType -> External,
1128 BlockName -> MNSBLOCK,
1129 OrderBlock -> 8,
1130 Value -> 3.438,
1131 TeX -> Subscript[\[Delta]R,CP],
1132 Description -> "Right MNS CP-Violation Phase"},
1133
1134
1135 (* Internal Parameters *)
1136
1137 cL12 == {
1138 ParameterType -> Internal,
1139 BlockName -> CKMBLOCK,
1140 Value -> Sqrt[1-sL12^2],
1141 TeX -> Subscript[cL,12],
1142 Description -> "Quark mixing: CosTheta12,PDG-94"},
1143 cL23 == {
1144 ParameterType -> Internal,
1145 BlockName -> CKMBLOCK,
1146 Value -> Sqrt[1-sL23^2],
1147 TeX -> Subscript[cL,23],
1148 Description -> "Quark mixing: CosTheta23,PDG-94"},
1149 cL13 == {
1150 ParameterType -> Internal,
1151 BlockName -> CKMBLOCK,
1152 Value -> Sqrt[1-sL13^2],
1153 TeX -> Subscript[cL,13],
1154 Description -> "Quark mixing: CosTheta13,PDG-94"},
1155
1156 cR12 == {
1157 ParameterType -> Internal,
1158 BlockName -> CKMBLOCK,
1159 Value -> Sqrt[1-sR12^2],
1160 TeX -> Subscript[cR,12],
1161 Description -> "Right Quark mixing: CosThetaR12,PDG-194"},
1162 cR23 == {
1163 ParameterType -> Internal,
1164 BlockName -> CKMBLOCK,
1165 Value -> Sqrt[1-sR23^2],
1166 TeX -> Subscript[cR,23],
1167 Description -> "Right Quark mixing: CosThetaR23,PDG-194"},
1168 cR13 == {
1169 ParameterType -> Internal,
1170 BlockName -> CKMBLOCK,
1171 Value -> Sqrt[1-sR13^2],
1172 TeX -> Subscript[cR,13],
1173 Description -> "Right Quark mixing: CosThetaR13,PDG-194"},
1174
1175 ccL12 == {
1176 ParameterType -> Internal,
1177 BlockName -> MNSBLOCK,
1178 Value -> Sqrt[1-ssL12^2],
1179 TeX -> Subscript[ccL,12],
1180 Description -> "Left Lepton mixing: CosTheta12,PDG-94"},
1181 ccL23 == {
1182 ParameterType -> Internal,
1183 BlockName -> MNSBLOCK,
1184 Value -> Sqrt[1-ssL23^2],
1185 TeX -> Subscript[ccL,23],
1186 Description -> "Left Lepton mixing: CosTheta23,PDG-94"},
1187 ccL13 == {
1188 ParameterType -> Internal,
1189 BlockName -> MNSBLOCK,
1190 Value -> Sqrt[1-ssL13^2],
1191 TeX -> Subscript[ccL,13],
1192 Description -> "Left Lepton mixing: CosTheta13,PDG-94"},
1193
1194 ccR12 == {
1195 ParameterType -> Internal,
1196 BlockName -> MNSBLOCK,
1197 Value -> Sqrt[1-ssR12^2],
1198 TeX -> Subscript[ccR,12],
1199 Description -> "Right Lepton mixing: CosThetaR12,PDG-194"},
1200 ccR23 == {
1201 ParameterType -> Internal,
1202 BlockName -> MNSBLOCK,
1203 Value -> Sqrt[1-ssR23^2],
1204 TeX -> Subscript[ccR,23],
1205 Description -> "Right Lepton mixing: CosThetaR23,PDG-194"},
1206 ccR13 == {
1207 ParameterType -> Internal,
1208 BlockName -> MNSBLOCK,
1209 Value -> Sqrt[1-ssR13^2],
1210 TeX -> Subscript[ccR,13],
1211 Description -> "Right Lepton mixing: CosThetaR13,PDG-194"},
1212
1213 CKML == {
1214 Indices -> {Index[Generation],Index[Generation]},
1215 TensorClass -> CKML,
1216 Unitary -> True,
1217 Value -> {CKML[1,1] -> cL12*cL13, CKML[1,2] -> sL12*cL13, CKML[1,3] -> sL13 Exp[-I delL13],
1218 CKML[2,1] -> -sL12*cL23-cL12*sL23*sL13 Exp[I delL13], CKML[2,2] -> cL12*cL23-sL12*sL23*sL13 Exp[I delL13], CKML[2,3] -> sL23*cL13,
1219 CKML[3,1] -> sL12*sL23-cL12*cL23*sL13 Exp[I delL13], CKML[3,2] -> -cL12*sL23-sL12*cL23*sL13 Exp[I delL13], CKML[3,3] -> cL23*cL13},
1220 TeX -> Superscript[VL,CKM],
1221 Description -> "Left CKM-Matrix"},
1222 CKMR == {
1223 Indices -> {Index[Generation],Index[Generation]},
1224 TensorClass -> CKMR,
1225 Unitary -> True,
1226 Value -> {CKMR[1,1] -> cR12*cR13, CKMR[1,2] -> sR12*cL13, CKMR[1,3] -> sR13 Exp[-I delR13],
1227 CKMR[2,1] -> -sR12*cR23-cR12*sR23*sR13 Exp[I delR13], CKMR[2,2] -> cR12*cR23-sR12*sR23*sR13 Exp[I delR13], CKMR[2,3] -> sR23*cR13,
1228 CKMR[3,1] -> sR12*sR23-cR12*cR23*sR13 Exp[I delR13], CKMR[3,2] -> -cR12*sR23-sR12*cR23*sR13 Exp[I delR13], CKMR[3,3] -> cR23*cR13},
1229 TeX -> Superscript[VR,CKM],
1230 Description -> "Right CKM-Matrix"},
1231(*
1232 MNSTB == {
1233 ParameterType -> Internal,
1234 Indices -> {Index[Generation],Index[Generation]},
1235 Unitary -> True,
1236 Value -> {MNSTB[1,1] -> Sqrt[2/3.] , MNSTB[1,2] -> 1/ Sqrt[3.], MNSTB[1,3] -> 0,
1237 MNSTB[2,1] -> -1/ Sqrt[6.], MNSTB[2,2] -> 1/ Sqrt[3.], MNSTB[2,3] -> 1/ Sqrt[2.],
1238 MNSTB[3,1] -> 1/ Sqrt[6.], MNSTB[3,2] -> -1/ Sqrt[3.], MNSTB[3,3] -> 1/ Sqrt[2.]},
1239 TeX -> Subscript[Superscript[UL,MNS],TB],
1240 Description -> "TriBimaximal MNS-Matrix"},
1241*)
1242 MNSL == {
1243 ParameterType -> Internal,
1244 Indices -> {Index[Generation],Index[Generation]},
1245 Unitary -> True,
1246 Value -> { MNSL[1,1] -> ccL12*ccL13, MNSL[1,2] -> ssL12*cL13, MNSL[1,3] -> ssL13 Exp[-I delLCP],
1247 MNSL[2,1] -> -ssL12*ccL23-ccL12*ssL23*ssL13 Exp[I delLCP], MNSL[2,2] -> ccL12*ccL23-ssL12*ssL23*ssL13 Exp[I delLCP], MNSL[2,3] -> ssL23*ccL13,
1248 MNSL[3,1] -> ssL12*ssL23-ccL12*ccL23*ssL13 Exp[I delLCP], MNSL[3,2] -> -ccL12*ssL23-ssL12*ccL23*ssL13 Exp[I delLCP], MNSL[3,3] -> ccL23*ccL13},
1249 TeX -> Superscript[UL,MNS],
1250 Description -> "Left MNS-Matrix"},
1251 MNSR == {
1252 ParameterType -> Internal,
1253 Indices -> {Index[Generation],Index[Generation]},
1254 Unitary -> True,
1255 Value -> { MNSR[1,1] -> ccR12*ccR13, MNSR[1,2] -> ssR12*cR13, MNSR[1,3] -> ssR13 Exp[-I delRCP],
1256 MNSR[2,1] -> -ssR12*ccR23-ccR12*ssR23*ssR13 Exp[I delRCP], MNSR[2,2] -> ccR12*ccR23-ssR12*ssR23*ssR13 Exp[I delRCP], MNSR[2,3] -> ssR23*ccR13,
1257 MNSR[3,1] -> ssR12*ssR23-ccR12*ccR23*ssR13 Exp[I delRCP], MNSR[3,2] -> -ccR12*ssR23-ssR12*ccR23*ssR13 Exp[I delRCP], MNSR[3,3] -> ccR23*ccR13},
1258 TeX -> Superscript[UR,MNS],
1259 Description -> "Right MNS-Matrix"},
1260 aEW == {
1261 ParameterType -> Internal,
1262 Value -> 1/aEWM1,
1263 InteractionOrder -> {QED,2},
1264 TeX -> Subscript[\[Alpha],EW],
1265 Description -> "Electroweak coupling constant" },
1266 ee == {
1267 ParameterType -> Internal,
1268 Value -> Sqrt[4Pi aEW],
1269 InteractionOrder -> {QED,1},
1270 TeX -> e,
1271 Description -> "Electric coupling constant" },
1272 sw == {
1273 ParameterType -> Internal,
1274 Value -> Sqrt[sw2],
1275 InteractionOrder -> {QED,0},
1276 TeX -> Subscript[s,w],
1277 Description -> "Sine of the Weinberg angle" },
1278 cw == {
1279 ParameterType -> Internal,
1280 Value -> Sqrt[1-sw2],
1281 InteractionOrder -> {QED,0},
1282 TeX -> Subscript[c,w],
1283 Description -> "Cosine of the Weinberg angle" },
1284 gL == {
1285 ParameterType -> Internal,
1286 Definitions -> {gL->ee/sw},
1287 InteractionOrder -> {QED,1},
1288 TeX -> Subscript[g,L],
1289 Description -> "Weak coupling constant at the Z pole" },
1290 gBL == {
1291 ParameterType -> Internal,
1292 Definitions -> {gBL->gL*gR*sw/Sqrt[cw^2*gR^2-sw^2*gL^2]},
1293 InteractionOrder -> {QED,1},
1294 TeX -> Subscript[g,BL],
1295 Description -> "B-L coupling constant at the Z pole" },
1296 gY == {
1297 ParameterType -> Internal,
1298 Definitions -> {gY->gR*gBL/Sqrt[gR^2+gBL^2]},
1299 InteractionOrder -> {QED,1},
1300 TeX -> Subscript[g,Y],
1301 Description -> "Hypercharge coupling constant at the Z pole" },
1302 gs == {
1303 ParameterType -> Internal,
1304 Value -> Sqrt[4Pi aS],
1305 InteractionOrder -> {QCD,1},
1306 TeX -> Subscript[g,s],
1307 ParameterName -> G,
1308 Description -> "Strong coupling constant at the Z pole" },
1309 vev == {
1310 ParameterType -> Internal,
1311 Value -> 2*MW*sw/ee,
1312 InteractionOrder -> {QED,-1},
1313 TeX -> v,
1314 Description -> "EW SSB scale (v)" },
1315 vp == {
1316 ParameterType -> Internal,
1317 Value -> Sqrt[(16 MZ^2 (MW^2-MZ^2) (1+tb^2)^2 + 4 gR^2 MZ^2 (1+tb^2) (MZ^2 (-1+sw^2) (1+tb^2) + MW^2 (1+(1-2 sw^2) tb^2)) vev^2 + gR^4 MW^2 tb^4 vev^4)/(gR^4 (MW^2+MZ^2 (-1+sw^2)) (1+tb^2)^2 vev^2)],
1318 InteractionOrder -> {QED,-1},
1319 TeX -> v',
1320 Description -> "Weak Right and BL SSB scale (v')" },
1321 MWp == {
1322 ParameterType -> Internal,
1323 Value -> gR*vp/2,
1324 InteractionOrder -> {QED,0},
1325 TeX -> Subscript[M,W'],
1326 Description -> "W' mass" },
1327 sb == {
1328 ParameterType -> Internal,
1329 Value -> tb/Sqrt[1+tb^2],
1330 InteractionOrder -> {QED,0},
1331 TeX -> Subscript[s,\[Beta]],
1332 Description -> "Sine beta" },
1333 cb == {
1334 ParameterType -> Internal,
1335 Value -> 1/Sqrt[1+tb^2],
1336 InteractionOrder -> {QED,0},
1337 TeX -> Subscript[c,\[Beta]],
1338 Description -> "Cosine beta" },
1339 k == {
1340 ParameterType -> Internal,
1341 Value -> vev*sb,
1342 InteractionOrder -> {QED,-1},
1343 Description -> "Higgs vacuum expectation value" },
1344 vL == {
1345 ParameterType -> Internal,
1346 Value -> vev*cb,
1347 InteractionOrder -> {QED,-1},
1348 TeX -> Subscript[v,L],
1349 Description -> "Higgs vacuum expectation value" },
1350 vR == {
1351 ParameterType -> Internal,
1352 Value -> Sqrt[vp^2-k^2],
1353 InteractionOrder -> {QED,-1},
1354 TeX -> Subscript[v,R],
1355 Description -> "Higgs vacuum expectation value" },
1356 lm4 == {
1357 ParameterType -> Internal,
1358 Value -> lm3-(k mu3)/(Sqrt[2]vL vR),
1359 InteractionOrder -> {QED,2},
1360 TeX -> Subscript[\[Lambda],4],
1361 Description -> "Higgs quartic coupling" },
1362 tz == {
1363 ParameterType -> Internal,
1364 Value -> k/vR,
1365 InteractionOrder -> {QED,0},
1366 TeX -> Subscript[t,\[Zeta]],
1367 Description -> "vevs ratio: Tan zeta" },
1368 sz == {
1369 ParameterType -> Internal,
1370 Value -> tz/Sqrt[1+tz^2],
1371 InteractionOrder -> {QED,0},
1372 TeX -> Subscript[s,\[Zeta]],
1373 Description -> "Sine zeta" },
1374 cz == {
1375 ParameterType -> Internal,
1376 Value -> 1/Sqrt[1+tz^2],
1377 InteractionOrder -> {QED,0},
1378 TeX -> Subscript[c,\[Zeta]],
1379 Description -> "Cosine zeta" },
1380 MLL == {
1381 ParameterType -> Internal,
1382 Value -> MW^2/(1-sw2),
1383 InteractionOrder -> {QED,0},
1384 Description -> "ZZ'-Mixing" },
1385 MLR == {
1386 ParameterType -> Internal,
1387 Value -> -(MW^2(cw^2 gR^2 sb^2-sw2 gL^2))/(gL(1-sw2)Sqrt[cw^2 gR^2-sw2 gL^2]),
1388 InteractionOrder -> {QED,0},
1389 Description -> "ZZ'-Mixing" },
1390 MRR == {
1391 ParameterType -> Internal,
1392 Value -> (1/4)(cp^2 gR^2 k^2 + gBL^2 sp^2 vL^2 + (cp gR + gBL sp)^2 vR^2 ),
1393 InteractionOrder -> {QED,0},
1394 Description -> "ZZ'-Mixing" },
1395 cth == {
1396 ParameterType -> Internal,
1397 Value -> 1/Sqrt[1+(MRR-MLL+Sqrt[4MLR^2+(MRR-MLL)^2])^2/(4MLR^2)],
1398 InteractionOrder -> {QED,0},
1399 TeX -> Subscript[s,\[Theta]],
1400 Description -> "ZZ'-Mixing: Sine theta" },
1401 sth == {
1402 ParameterType -> Internal,
1403 Value -> (MRR-MLL+Sqrt[4MLR^2+(MRR-MLL)^2])/(2MLR Sqrt[1+(MRR-MLL+Sqrt[4MLR^2+(MRR-MLL)^2])^2/(4MLR^2)]),
1404 InteractionOrder -> {QED,0},
1405 TeX -> Subscript[c,\[Theta]],
1406 Description -> "ZZ'-Mixing: Cosine theta" },
1407 sp == {
1408 ParameterType -> Internal,
1409 Value -> gBL/Sqrt[gBL^2+gR^2](*gY/gR*),
1410 InteractionOrder -> {QED,0},
1411 TeX -> Subscript[s,\[Phi]],
1412 Description -> "ZZ'-Mixing: Sine phi" },
1413 cp == {
1414 ParameterType -> Internal,
1415 Value -> gR/Sqrt[gBL^2+gR^2](*gY/gBL*),
1416 InteractionOrder -> {QED,0},
1417 TeX -> Subscript[c,\[Phi]],
1418 Description -> "ZZ'-Mixing: Cosine phi" },
1419 mhp1 == {
1420 ParameterType -> Internal,
1421 Value -> Sqrt[(k vL (al3-al2)-(mu3 vR)/Sqrt[2])(vev^2/(k vL))],
1422 InteractionOrder -> {QED,0},
1423 TeX -> Subscript[m,hp1],
1424 Description -> "Charged Higgs mass" },
1425 mhm2 == {
1426 ParameterType -> Internal,
1427 Value -> Sqrt[(k vR (al3-al2)-(mu3 vL)/Sqrt[2])(vp^2/(k vR))],
1428 InteractionOrder -> {QED,0},
1429 TeX -> Subscript[m,hp2],
1430 Description -> "Charged Higgs mass" },
1431 ma1 == {
1432 ParameterType -> Internal,
1433 Value -> Sqrt[2k^2 lm2+(vL^2+vR^2)(al3-al2)-(vL vR mu3)/(Sqrt[2]k)],
1434 InteractionOrder -> {QED,0},
1435 TeX -> Subscript[m,A1],
1436 Description -> "Pseudoscalar Higgs mass" },
1437 ma2 == {
1438 ParameterType -> Internal,
1439 Value -> Sqrt[-((vL^2 vR^2+k^2(vL^2+vR^2))mu3)/(Sqrt[2]k vL vR)],
1440 InteractionOrder -> {QED,0},
1441 TeX -> Subscript[m,A2],
1442 Description -> "Pseudoscalar Higgs mass" },
1443 mh1 == {
1444 ParameterType -> Internal,
1445 Value -> ma1,
1446 InteractionOrder -> {QED,0},
1447 TeX -> Subscript[m,h1],
1448 Description -> "Higgs Mass" },
1449 a0 == {
1450 ParameterType -> Internal,
1451 Value -> 2(vL^2+vR^2)lm3-((vL^2 vR^2+k^2(vL^2+vR^2))mu3)/(Sqrt[2]vL vR k),
1452 InteractionOrder -> {QED,0},
1453 Description -> "Higgs mass parameter" },
1454 a1 == {
1455 ParameterType -> Internal,
1456 Value -> 2k^2,
1457 InteractionOrder -> {QED,-2},
1458 Description -> "Higgs mass parameter" },
1459 b0 == {
1460 ParameterType -> Internal,
1461 Value -> ((vL^2+vR^2)(4k^3vL vR(al1+al2)^2+Sqrt[2]vL^2 vR^2lm3 mu3)+Sqrt[2]k^2(4vL^2vR^2(al1+al2)+(vL^2-vR^2)^2lm3)mu3)/(k vL vR),
1462 InteractionOrder -> {QED,0},
1463 Description -> "Higgs mass parameter" },
1464 b1 == {
1465 ParameterType -> Internal,
1466 Value -> (vL^2+vR^2)(4k^3vL vR(-lm3)+Sqrt[2]k^4mu3)/(k vL vR),
1467 InteractionOrder -> {QED,-2},
1468 Description -> "Higgs mass parameter" },
1469 c0 == {
1470 ParameterType -> Internal,
1471 Value -> 2Sqrt[2]k^3(vL^2-vR^2)^2(al1+al2)^2mu3/(vL vR),
1472 InteractionOrder -> {QED,0},
1473 Description -> "Higgs mass parameter" },
1474 c1 == {
1475 ParameterType -> Internal,
1476 Value -> -2Sqrt[2]k^3(vL^2-vR^2)^2 lm3 mu3/(vL vR),
1477 InteractionOrder -> {QED,-2},
1478 Description -> "Higgs mass parameter" },
1479 lm1 == {
1480 ParameterType -> Internal,
1481 Value -> -(-mh0^6+a0 mh0^4+b0 mh0^2+c0)/(a1 mh0^4+b1 mh0^2+c1),
1482 InteractionOrder -> {QED,2},
1483 TeX -> Subscript[\[Lambda],1],
1484 Description -> "Higgs quartic coupling" },
1485 mu1 == {
1486 ParameterType -> Internal,
1487 Value -> (vL^2+vR^2)(al1+al2)+k^2 lm1+(vL vR mu3)/(Sqrt[2]k),
1488 InteractionOrder -> {QED,0},
1489 TeX -> Superscript[Subscript[\[Mu],1],2],
1490 Description -> "Higgs quadratic coupling" },
1491 mu2 == {
1492 ParameterType -> Internal,
1493 Value -> k^2(al1+al2)+(vL^2+vR^2)lm3,
1494 InteractionOrder -> {QED,0},
1495 TeX -> Superscript[Subscript[\[Mu],2],2],
1496 Description -> "Higgs quadratic coupling" },
1497 MAZZp == {
1498 ParameterType -> Internal,
1499 Indices -> {Index[ngauge],Index[ngauge]},
1500 Value -> {MAZZp[1,1]->(1/4)gL^2(vL^2+k^2), MAZZp[1,2]->-(1/4)gL gR k^2, MAZZp[1,3]->-(1/4)gL gBL vL^2,
1501 MAZZp[2,1]->MAZZp[1,2], MAZZp[2,2]->(1/4)gR^2(vR^2+k^2), MAZZp[2,3]->-(1/4)gR gBL vR^2,
1502 MAZZp[3,1]->MAZZp[1,3], MAZZp[3,2]->MAZZp[2,3], MAZZp[3,3]->(1/4)gBL^2(vL^2+vR^2)},
1503 InteractionOrder -> {QED,0},
1504 TeX -> MG0,
1505 Description -> "Neutral Gauge Bosons Mass Matrix" },
1506 TAZZp == {
1507 ParameterType -> Internal,
1508 Indices -> {Index[ngauge],Index[ngauge]},
1509 Value -> {TAZZp[1,1]-> cw, TAZZp[1,2]->sw, TAZZp[1,3]->0,
1510 TAZZp[2,1]->-sp sw, TAZZp[2,2]->sp cw, TAZZp[2,3]->cp,
1511 TAZZp[3,1]->-cp sw, TAZZp[3,2]->cp cw, TAZZp[3,3]->-sp, Conjugate[TAZZp[i_,j_]]->TAZZp[i,j] },
1512 InteractionOrder -> {QED,0},
1513 TeX -> TG0,
1514 Description -> "Neutral Gauge Bosons Rotation Matrix" },
1515 MZp == {
1516 ParameterType -> Internal,
1517 Value -> Sqrt[Tr[Table[MAZZp[i,j],{i,3},{j,3}]]-MZ^2],
1518 InteractionOrder -> {QED,0},
1519 TeX -> Subscript[M,Z'],
1520 Description -> "Z'-Mass" },
1521(*
1522 MLL == {
1523 ParameterType -> Internal,
1524 Value -> (Transpose[Table[TAZZp[i,j],{i,3},{j,3}]].Table[MAZZp[i,j],{i,3},{j,3}].Table[TAZZp[i,j],{i,3},{j,3}])[[1,1]],
1525 InteractionOrder -> {QED,0},
1526 Description -> "ZZ'-Mixing" },
1527 MLR == {
1528 ParameterType -> Internal,
1529 Value -> (Transpose[Table[TAZZp[i,j],{i,3},{j,3}]].Table[MAZZp[i,j],{i,3},{j,3}].Table[TAZZp[i,j],{i,3},{j,3}])[[1,3]],
1530 InteractionOrder -> {QED,0},
1531 Description -> "ZZ'-Mixing" },
1532 MRR == {
1533 ParameterType -> Internal,
1534 Value -> (Transpose[Table[TAZZp[i,j],{i,3},{j,3}]].Table[MAZZp[i,j],{i,3},{j,3}].Table[TAZZp[i,j],{i,3},{j,3}])[[3,3]],
1535 InteractionOrder -> {QED,0},
1536 Description -> "ZZ'-Mixing" },
1537 tth2 == {
1538 ParameterType -> Internal,
1539 Value -> 2MLR/(MLL-MRR),
1540 InteractionOrder -> {QED,0},
1541 TeX -> Subscript[t,2\[Theta]],
1542 Description -> "ZZ'-Mixing: Tan 2 theta" },
1543 tth == {
1544 ParameterType -> Internal,
1545 Value -> (-1+Sqrt[1+tth2^2])/tth2,
1546 InteractionOrder -> {QED,0},
1547 TeX -> Subscript[t,\[Theta]],
1548 Description -> "ZZ'-Mixing: Tan theta" },
1549 cth == {
1550 ParameterType -> Internal,
1551 Value -> 1/Sqrt[1+(MRR-MLL+Sqrt[4MLR^2+(MRR-MLL)^2])^2/(4MLR^2)],
1552 Value -> 1/Sqrt[1+tth^2],
1553 InteractionOrder -> {QED,0},
1554 TeX -> Subscript[s,\[Theta]],
1555 Description -> "ZZ'-Mixing: Sine theta" },
1556 sth == {
1557 ParameterType -> Internal,
1558 Value -> (MRR-MLL+Sqrt[4MLR^2+(MRR-MLL)^2])/(2MLR Sqrt[1+(MRR-MLL+Sqrt[4MLR^2+(MRR-MLL)^2])^2/(4MLR^2)]),
1559 Value -> tth/Sqrt[1+tth^2],
1560 InteractionOrder -> {QED,0},
1561 TeX -> Subscript[c,\[Theta]],
1562 Description -> "ZZ'-Mixing: Cosine theta" },
1563*)
1564 MH == {
1565 ParameterType -> Internal,
1566 Indices -> {Index[higgs],Index[higgs]},
1567 Value -> {MH[1,1]->2k^2 lm1-(vL vR mu3)/(Sqrt[2]k), MH[1,2]->2k vL(al1+al2)+(vR mu3)/(Sqrt[2]),
1568 MH[1,3]->2k vR(al1+al2)+(vL mu3)/(Sqrt[2]),
1569 MH[2,1]->MH[1,2], MH[2,2]->2vL^2 lm3-(k vR mu3)/(Sqrt[2]vL), MH[2,3]->2vL vR lm3-(k mu3)/(Sqrt[2]),
1570 MH[3,1]->MH[1,3], MH[3,2]->MH[2,3], MH[3,3]->2vR^2 lm3-(k vL mu3)/(Sqrt[2]vR)},
1571 InteractionOrder -> {QED,0},
1572 TeX -> Mh,
1573 Description -> "CP-even Higgs Mass Matrix" },
1574 trmh == {
1575 ParameterType -> Internal,
1576 Value -> Tr[Table[MH[i,j],{i,3},{j,3}]],
1577 InteractionOrder -> {QED,0},
1578 TeX -> TrMH,
1579 Description -> "CP-even Higgs mass matrix trace" },
1580 detmh == {
1581 ParameterType -> Internal,
1582 Value -> Det[Table[MH[i,j],{i,3},{j,3}]],
1583 InteractionOrder -> {QED,0},
1584 TeX -> DetMH,
1585 Description -> "CP-even Higgs mass matrix determinant" },
1586 mh2 == {
1587 ParameterType -> Internal,
1588 Value -> Sqrt[(1/2)(trmh-mh0^2-Sqrt[(trmh-mh0^2)^2-4detmh/mh0^2])],
1589 InteractionOrder -> {QED,0},
1590 TeX -> Subscript[m,h2],
1591 Description -> "Higgs Mass" },
1592 mh3 == {
1593 ParameterType -> Internal,
1594 (* Value -> Sqrt[(1/2)(trmh-mh0^2+Sqrt[(trmh-mh0^2)^2-4detmh/mh0^2])],*)
1595 Value -> Sqrt[trmh-mh0^2-mh2^2],
1596 InteractionOrder -> {QED,0},
1597 TeX -> Subscript[m,h3],
1598 Description -> "Higgs Mass" },
1599 yl == {
1600 ParameterType -> Internal,
1601 Indices -> {Index[Generation],Index[Generation]},
1602 Definitions -> {yl[i_?NumericQ, j_?NumericQ] :> 0 /; (i =!= j)},
1603 Value -> {yl[1,1] -> Sqrt[2]yme / (vev*sb), yl[2,2] -> Sqrt[2]ymm / (vev*sb), yl[3,3] -> Sqrt[2]ymtau / (vev*sb)},
1604 InteractionOrder -> {QED,1},
1605 ParameterName -> {yl[1,1] -> ye, yl[2,2] -> ym, yl[3,3] -> ytau},
1606 TeX -> Superscript[y,l],
1607 Description -> "Lepton Yukawa couplings" },
1608 yn == {
1609 ParameterType -> Internal,
1610 Indices -> {Index[Generation],Index[Generation]},
1611 Definitions -> {yn[i_?NumericQ, j_?NumericQ] :> 0 /; (i =!= j)},
1612 Value -> {yn[1,1] -> Sqrt[2]ymne / vR, yn[2,2] -> Sqrt[2]ymnm / vR, yn[3,3] -> Sqrt[2]ymntau / vR},
1613 InteractionOrder -> {QED,1},
1614 ParameterName -> {yn[1,1] -> yne, yn[2,2] -> ynm, yn[3,3] -> yntau},
1615 TeX -> Superscript[y,n],
1616 Description -> "Scotino Yukawa couplings" },
1617 ynu == {
1618 ParameterType -> Internal,
1619 Indices -> {Index[Generation],Index[Generation]},
1620 Definitions -> {ynu[i_?NumericQ, j_?NumericQ] :> 0 /; (i =!= j)},
1621 Value -> {ynu[1,1] -> Sqrt[2]ymnue / (vev*cb), ynu[2,2] -> Sqrt[2]ymnum / (vev*cb), ynu[3,3] -> Sqrt[2]ymnutau / (vev*cb)},
1622 InteractionOrder -> {QED,1},
1623 ParameterName -> {ynu[1,1] -> ynue, ynu[2,2] -> ynum, ynu[3,3] -> ynutau},
1624 TeX -> Superscript[y,\[Nu]],
1625 Description -> "Neutrino Yukawa couplings" },
1626 yu == {
1627 ParameterType -> Internal,
1628 Indices -> {Index[Generation],Index[Generation]},
1629 Definitions -> {yu[i_?NumericQ, j_?NumericQ] :> 0 /; (i =!= j)},
1630 Value -> {yu[1,1] -> Sqrt[2]ymup/(vev*sb), yu[2,2] -> Sqrt[2]ymc/(vev*sb), yu[3,3] -> Sqrt[2]ymt/(vev*sb)},
1631 InteractionOrder -> {QED,1},
1632 ParameterName -> {yu[1,1] -> yup, yu[2,2] -> yc, yu[3,3] -> yt},
1633 TeX -> Superscript[y,u],
1634 Description -> "Up-type Yukawa couplings" },
1635 yd == {
1636 ParameterType -> Internal,
1637 Indices -> {Index[Generation],Index[Generation]},
1638 Definitions -> {yd[i_?NumericQ, j_?NumericQ] :> 0 /; (i =!= j)},
1639 Value -> {yd[1,1] -> Sqrt[2]ymdo/(vev*cb), yd[2,2] -> Sqrt[2]yms/(vev*cb), yd[3,3] -> Sqrt[2]ymb/(vev*cb)},
1640 InteractionOrder -> {QED,1},
1641 ParameterName -> {yd[1,1] -> ydo, yd[2,2] -> ys, yd[3,3] -> yb},
1642 TeX -> Superscript[y,d],
1643 Description -> "Down-type Yukawa couplings" },
1644 ydd == {
1645 ParameterType -> Internal,
1646 Indices -> {Index[Generation],Index[Generation]},
1647 Definitions -> {ydd[i_?NumericQ, j_?NumericQ] :> 0 /; (i =!= j)},
1648 Value -> {ydd[1,1] -> Sqrt[2]ymddo/vR, ydd[2,2] -> Sqrt[2]ymds/vR, ydd[3,3] -> Sqrt[2]ymdb/vR},
1649 InteractionOrder -> {QED,1},
1650 ParameterName -> {ydd[1,1] -> yddo, ydd[2,2] -> yds, ydd[3,3] -> ydb},
1651 TeX -> Superscript[y,d'],
1652 Description -> "Exotic-Down-type Yukawa couplings" },
1653
1654 (* Higgs Functions Definitions for Neutral Higgs Mixing *)
1655
1656 XH == {
1657 ParameterType -> Internal,
1658 Indices -> {Index[higgs]},
1659 Value -> {XH[1] -> mh0^2, XH[2] -> mh2^2, XH[3] -> mh3^2 },
1660 InteractionOrder -> {QED,0},
1661 TeX -> xh,
1662 Description -> "Squared CP-even Higgs Masses" },
1663 DH == {
1664 ParameterType -> Internal,
1665 Indices -> {Index[higgs],Index[higgs]},
1666 Value -> {DH[i_,k_]:>Det[Drop[Table[MH[ii,jj]-KroneckerDelta[ii,jj] XH[k],{ii,3},{jj,3}],{i},{1}]]},
1667 InteractionOrder -> {QED,0},
1668 TeX -> dh,
1669 Description -> "Higgs mass matrix minors determinants" },
1670 FH == {
1671 ParameterType -> Internal,
1672 Indices -> {Index[higgs],Index[higgs]},
1673 Value -> {FH[i_,k_]:>(-1)^(3+i)DH[i,k]/DH[3,k]},
1674 InteractionOrder -> {QED,0},
1675 TeX -> fh,
1676 Description -> "Higgs Mixing functions" },
1677 FVH == {
1678 ParameterType -> Internal,
1679 Indices -> {Index[higgs],Index[higgs],Index[higgs]},
1680 Value -> {FVH[kk_?NumericQ, ii_?NumericQ, jj_?NumericQ] :> If[kk==1,FH[ii,jj],If[jj<kk,0,FVH[kk-1,ii,jj]-FVH[kk-1,ii,kk-1] Sum[FVH[1,ll,jj]FVH[kk-1,ll,kk-1],{ll,3}]/Sum[FVH[kk-1,ll,kk-1]^2,{ll,3}]]] },
1681 InteractionOrder -> {QED,0},
1682 TeX -> fvh,
1683 Description -> "Higgs Mixing recursion functions" },
1684 VH == {
1685 ParameterType -> Internal,
1686 Indices -> {Index[higgs],Index[higgs]},
1687 Value -> {VH[jj_?NumericQ, kk_?NumericQ] :> FVH[kk,jj,kk]/Sqrt[Sum[FVH[kk,ll,kk]^2,{ll,3}]]},
1688 InteractionOrder -> {QED,0},
1689 ParameterName -> {VH[1,1] -> THP0, VH[1,2] -> THP2, VH[1,3] -> THP3,
1690 VH[2,1] -> THL0, VH[2,2] -> THL2, VH[2,3] -> THL3,
1691 VH[3,1] -> THR0, VH[3,2] -> THR2, VH[3,3] -> THR3 },
1692 TeX -> vh,
1693 Description -> "CP-even Higgs mixing" },
1694
1695(* Definitions for the effective hYY and hgg Lagrangian *)
1696
1697 MP == {
1698 ParameterType -> Internal,
1699 Indices -> {Index[bjorkeneff]},
1700 Value -> {MP[1]->MT, MP[2]->MW, MP[3]->MWp, MP[4]->mhp1, MP[5]->mhm2},
1701 InteractionOrder -> {QED,0},
1702 TeX -> mp,
1703 Description -> "Loop particle masses" },
1704 XP == {
1705 ParameterType -> Internal,
1706 Indices -> {Index[bjorkeneff]},
1707 Value -> {XP[j_]:>mh0^2/(4MP[j]^2)},
1708 InteractionOrder -> {QED,0},
1709 TeX -> xp,
1710 Description -> "Bjorken parameter" },
1711 FP == {
1712 ParameterType -> Internal,
1713 Indices -> {Index[bjorkeneff]},
1714 Value -> {FP[1]:>Module[{},2(xf+(xf-1)ArcSin[Sqrt[xf]]^2) xf^(-2)/.xf->XP[1]],
1715 FP[2]:>Module[{},-(2xv^2+3xv+3(2xv-1)ArcSin[Sqrt[xv]]^2) xv^(-2)/.xv->XP[2]],
1716 FP[3]:>Module[{},-(2xv^2+3xv+3(2xv-1)ArcSin[Sqrt[xv]]^2) xv^(-2)/.xv->XP[3]],
1717 FP[4]:>Module[{},-(xs-ArcSin[Sqrt[xs]]^2) xs^(-2)/.xs->XP[4]],
1718 FP[5]:>Module[{},-(xs-ArcSin[Sqrt[xs]]^2) xs^(-2)/.xs->XP[5]]},
1719 InteractionOrder -> {QED,0},
1720 TeX -> fp,
1721 Description -> "Loop functions" },
1722 mu11 == {
1723 ParameterType -> Internal,
1724 Value -> THP0(lm1 k cb^2+(al3-al2)vL sb cb+(al1+al3)k sb^2)+THL0(lm3 vL sb^2+(al3-al2)k sb cb+(al1+al3)vL cb^2)+THR0(lm4 vR sb^2-Sqrt[2]mu3 sb cb+(al1+al2)vR cb^2),
1725 TeX -> Subscript[\[Mu],11],
1726 Description -> "Neutral-Charged Higgs Coupling" },
1727 mu22 == {
1728 ParameterType -> Internal,
1729 Value -> THP0(lm1 k cz^2+(al3-al2)vR sz cz+(al1+al3)k sz^2)+THR0(lm3 vR sz^2+(al3-al2)k sz cz+(al1+al3)vR cz^2)+THL0(lm4 vL sz^2-Sqrt[2]mu3 sz cz+(al1+al2)vL cz^2),
1730 TeX -> Subscript[\[Mu],22],
1731 Description -> "Neutral-Charged Higgs Coupling" },
1732 LP == {
1733 ParameterType -> Internal,
1734 Indices -> {Index[bjorkeneff]},
1735 Value -> {LP[1]:>Module[{}, (3)(2/3)^2 (THP0/sb) ff/.ff->FP[1]],
1736 LP[2]:>Module[{}, (THP0 sb+THL0 cb) fv/.fv->FP[2]],
1737 LP[3]:>Module[{}, (vev/vp)(THP0 sz+THR0 cz) fv/.fv->FP[3]],
1738 LP[4]:>Module[{}, (-mu11 vev/MP[4]^2) fs/.fs->FP[4]],
1739 LP[5]:>Module[{}, (-mu22 vev/MP[5]^2) fs/.fs->FP[5]]},
1740 InteractionOrder -> {QED,0},
1741 TeX -> Subscript[\[Lambda],hpp],
1742 Description -> "loop-induced higgs Couplings" },
1743 Ghaa == {
1744 ParameterType -> Internal,
1745 Value -> ( 1 ) (1/2)Sqrt[Gf Sqrt[2]](aEW/(4Pi)) Sum[LP[i],{i,5}],
1746 InteractionOrder -> {HIW,1},
1747 TeX -> Subscript[G,h\[Gamma]\[Gamma]],
1748 Description -> "loop-induced h-aa Coupling" },
1749 Ghgg == {
1750 ParameterType -> Internal,
1751 Value -> (3/8) (1/2)Sqrt[Gf Sqrt[2]] (aS/(4Pi)) LP[1],
1752 InteractionOrder -> {HIG,1},
1753 TeX -> Subscript[G,hgg],
1754 Description -> "loop-induced h-gg Coupling" }
1755};
1756