ALRM: ALRM_LO.fr

File ALRM_LO.fr, 61.3 KB (added by mashry, 3 months ago)

Model File

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