ALRM: ALRM_LO.fr

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