331: 331_1.fr

File 331_1.fr, 60.7 KB (added by zhangdongming, 3 years ago)

The 331 model file for beta=+1/\sqrt(3).

Line 
1(***************************************************************************************************************)
2(******      This is the FeynRules mod-file for the 331 model where beta equals to 1/Sqrt[3] ******)
3(******                                                                                                   ******)
4(******     Authors: Dongming Zhang          ******)
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 = "331_1 Model";
16
17M$Information = {
18  Authors      -> {"Dongming Zhang"},
19  Version      -> "1.0.0",
20  Date         -> "22. 01. 2014",
21  Institutions -> {"Peking University"},
22  Emails       -> {"zhangdongming@pku.edu.cn"},
23  URLs         -> "http://feynrules.irmp.ucl.ac.be/wiki/331"
24};
25
26FeynmanGauge = True;
27
28(* ************************** *)
29(* *****  Change  log   ***** *)
30(* ************************** *)
31
32(* ************************** *)
33(* *****      vevs      ***** *)
34(* ************************** *)
35M$vevs = { {Rho[2],v}, {Phi[1],v2}, {Chi[3],v3} };
36
37(* ************************** *)
38(* *****  Gauge groups  ***** *)
39(* ************************** *)
40M$GaugeGroups = {
41  U1X  == {
42    Abelian          -> True, 
43    CouplingConstant -> gx,
44    GaugeBoson       -> K,
45    Charge           -> X
46  },
47  SU3L == {
48    Abelian           -> False,
49    CouplingConstant  -> gw,
50    GaugeBoson        -> Wi,
51    StructureConstant -> x,
52    Representations   -> {Ta,SU3T},
53    Definitions       -> {Ta[a_,b_,c_]->Gellmann[a,b,c]/2,FSU3L[i_,j_,k_]:> I x[i,j,k]},
54    SymmetricTensor   -> dSUN
55  },
56  ASU3L == {
57    Abelian           -> False,
58    CouplingConstant  -> gw,
59    GaugeBoson        -> WWi,
60    StructureConstant -> x,
61    Representations   -> {Tb,ASU3T},
62    Definitions       -> {Tb[a_,b_,c_]->-Gellmann[a,c,b]/2,FSU3L[i_,j_,k_]:> I x[i,j,k]},
63    SymmetricTensor   -> dSUN
64  },
65  SU3C == {
66    Abelian           -> False,
67    CouplingConstant  -> gs,
68    GaugeBoson        -> G,
69    StructureConstant -> f,
70    Representations   -> {T,Colour},
71    SymmetricTensor   -> dSUN
72  }
73};
74
75(* ************************** *)
76(* *** Gellmann matrices  *** *)
77(* ************************** *)
78
79Table[Gellmann[i, j, k] = 0, {i, 1, 8}, {j, 1, 3}, {k, 1, 3}] //
80  Flatten;
81
82Gellmann[1] = {{0, 1, 0}, {1, 0, 0}, {0, 0, 0}};
83Gellmann[2] = {{0, -I, 0}, {I, 0, 0}, {0, 0, 0}};
84Gellmann[3] = {{1, 0, 0}, {0, -1, 0}, {0, 0, 0}};
85Gellmann[4] = {{0, 0, 1}, {0, 0, 0}, {1, 0, 0}};
86Gellmann[5] = {{0, 0, -I}, {0, 0, 0}, {I, 0, 0}};
87Gellmann[6] = {{0, 0, 0}, {0, 0, 1}, {0, 1, 0}};
88Gellmann[7] = {{0, 0, 0}, {0, 0, -I}, {0, I, 0}};
89Gellmann[8] = 1/Sqrt[3] {{1, 0, 0}, {0, 1, 0}, {0, 0, -2}};
90
91Gellmann[1, 1, 2] = 1; Gellmann[1, 2, 1] = 1;
92Gellmann[2, 1, 2] = -I; Gellmann[2, 2, 1] = I;
93Gellmann[3, 1, 1] = 1; Gellmann[3, 2, 2] = -1;
94Gellmann[4, 1, 3] = 1; Gellmann[4, 3, 1] = 1;
95Gellmann[5, 1, 3] = -I; Gellmann[5, 3, 1] = I;
96Gellmann[6, 2, 3] = 1; Gellmann[6, 3, 2] = 1;
97Gellmann[7, 2, 3] = -I; Gellmann[7, 3, 2] = I;
98Gellmann[8, 1, 1] = 1/Sqrt[3]; Gellmann[8, 2, 2] = 1/Sqrt[3];
99Gellmann[8, 3, 3] = -2/Sqrt[3];
100
101
102Gellmann[i_Integer, j_Integer, k_Integer] := Gellmann[i][[j, k]];
103Gellmann[xx___, Index[_, i_Integer], yy___] := Gellmann[xx, i, yy];
104
105Gellmann /:
106  Gellmann[i1_, i2_, i3_?(Not[NumericQ[#]] &)] Gellmann[j1_, i3_,
107    j3_] :=     
108  Gellmann[i1, i2, 1] Gellmann[j1, 1, j3] +
109   Gellmann[i1, i2, 2] Gellmann[j1, 2, j3] +
110   Gellmann[i1, i2, 3] Gellmann[j1, 3, j3];
111
112Table[x[i, j, k] = 0, {i, 1, 8}, {j, 1, 8}, {k, 1, 8}] // Flatten;
113x[1, 2, 3] = 1; x[2, 3, 1] = 1; x[3, 1, 2] = 1;
114x[2, 1, 3] = -1; x[1, 3, 2] = -1; x[3, 2, 1] = -1;
115x[1, 5, 6] = -1/2; x[3, 6, 7] = -1/2; x[1, 7, 4] = -1/2;
116x[2, 6, 4] = -1/2; x[2, 7, 5] = -1/2; x[3, 5, 4] = -1/2;
117x[6, 1, 5] = -1/2; x[7, 3, 6] = -1/2; x[4, 1, 7] = -1/2;
118x[4, 2, 6] = -1/2; x[5, 2, 7] = -1/2; x[4, 3, 5] = -1/2;
119x[5, 6, 1] = -1/2; x[6, 7, 3] = -1/2; x[7, 4, 1] = -1/2;
120x[6, 4, 2] = -1/2; x[7, 5, 2] = -1/2; x[5, 4, 3] = -1/2;
121x[1, 6, 5] = 1/2; x[3, 7, 6] = 1/2; x[1, 4, 7] = 1/2;
122x[2, 4, 6] = 1/2; x[2, 5, 7] = 1/2; x[3, 4, 5] = 1/2;
123x[5, 1, 6] = 1/2; x[7, 6, 3] = 1/2; x[4, 7, 1] = 1/2;
124x[4, 6, 2] = 1/2; x[5, 7, 2] = 1/2; x[4, 5, 3] = 1/2;
125x[6, 5, 1] = 1/2; x[6, 3, 7] = 1/2; x[7, 1, 4] = 1/2;
126x[6, 2, 4] = 1/2; x[7, 2, 5] = 1/2; x[5, 3, 4] = 1/2;
127x[4, 5, 8] = Sqrt[3]/2; x[6, 7, 8] = Sqrt[3]/2;
128x[8, 4, 5] = Sqrt[3]/2; x[8, 6, 7] = Sqrt[3]/2;
129x[5, 8, 4] = Sqrt[3]/2; x[7, 8, 6] = Sqrt[3]/2;
130x[4, 8, 5] = -Sqrt[3]/2; x[6, 8, 7] = -Sqrt[3]/2;
131x[5, 4, 8] = -Sqrt[3]/2;
132x[7, 6, 8] = -Sqrt[3]/2; x[8, 7, 6] = -Sqrt[3]/2;
133x[8, 5, 4] = -Sqrt[3]/2;
134
135x /: x[ii___, Except[Index[___] | _?NumericQ, jj_], kk___] f_[aa___,
136    Index[name_, jj_], cc___] :=
137  x[ii, Index[name, jj], kk] f[aa, Index[name, jj], cc];
138x /: x[ii___, Except[Index[___] | _?NumericQ, jj_],
139    kk___] f_[aa___, Index[name_, jj_], cc___][ind___] :=
140  x[ii, Index[name, jj], kk] f[aa, Index[name, jj], cc][ind];
141x /: x[ii___, Except[Index[___] | _?NumericQ, jj_], kk___] f_[aa___,
142    g_[xx___, Index[name_, jj_], yy___], cc___] :=
143  x[ii, Index[name, jj], kk] f[aa, g[xx, Index[name, jj], yy], cc];
144x /: x[ii___, Except[Index[___] | _?NumericQ, jj_],
145    kk___] f_[aa___, g_[xx___, Index[name_, jj_], yy___], cc___][
146    ind___] :=
147  x[ii, Index[name, jj], kk] f[aa, g[xx, Index[name, jj], yy], cc][
148    ind];
149
150x[ii___, Except[_Index | _Done[Index] | _FV,
151    jj_?(Not[NumericQ[#]] &)], kk___, Index[name_, ll_], mm___] :=
152  x[ii, Index[name, jj], kk, Index[name, ll], mm];
153x[ii___, Index[name_, ll_], kk___,
154   Except[_Index | _Done[Index] | _FV, jj_?(Not[NumericQ[#]] &)],
155   mm___] := x[ii, Index[name, ll], kk, Index[name, jj], mm];
156
157x /: x[i_?((Not[NumericQ[#]] && Not[MatchQ[#, Index[_, _?NumericQ]]] &&
158          Not[MatchQ[#, Done[Index][_, _?NumericQ]]]) &), j_, k_] x[
159    i_, m_, n_] :=
160  x[1, j, k] x[1, m, n] + x[2, j, k] x[2, m, n] +
161   x[3, j, k] x[3, m, n] + x[4, j, k] x[4, m, n] +
162   x[5, j, k] x[5, m, n] + x[6, j, k] x[6, m, n] +
163   x[7, j, k] x[7, m, n] + x[8, j, k] x[8, m, n];
164
165x /: x[i_?((Not[NumericQ[#]] && Not[MatchQ[#, Index[_, _?NumericQ]]] &&
166          Not[MatchQ[#, Done[Index][_, _?NumericQ]]]) &), j_, k_] x[
167    m_, n_, i_] := x[i, j, k] x[i, m, n];
168x /: x[i_?((Not[NumericQ[#]] && Not[MatchQ[#, Index[_, _?NumericQ]]] &&
169          Not[MatchQ[#, Done[Index][_, _?NumericQ]]]) &), j_, k_] x[
170    m_, i_, n_] := x[i, j, k] x[i, n, m];
171x /: x[j_,
172    i_?((Not[NumericQ[#]] && Not[MatchQ[#, Index[_, _?NumericQ]]] &&
173         Not[MatchQ[#, Done[Index][_, _?NumericQ]]]) &), k_] x[m_, i_,
174     n_] := x[i, k, j] x[i, n, m];
175x /: x[j_,
176    i_?((Not[NumericQ[#]] && Not[MatchQ[#, Index[_, _?NumericQ]]] &&
177         Not[MatchQ[#, Done[Index][_, _?NumericQ]]]) &), k_] x[m_, n_,
178     i_] := x[i, k, j] x[i, m, n];
179x /: x[j_, k_,
180    i_?((Not[NumericQ[#]] && Not[MatchQ[#, Index[_, _?NumericQ]]] &&
181         Not[MatchQ[#, Done[Index][_, _?NumericQ]]]) &)] x[m_, n_,
182    i_] := x[i, j, k] x[i, m, n];
183
184x /: x[___, i_, ___, j_, ___] FV[a_, i_] FV[a_, j_] := 0;
185x /: x[___, i_, ___, j_, ___] del[del[_, i_], j_] := 0;
186x /: x[___, i_, ___, j_, ___] del[del[_, j_], i_] := 0;
187
188x[xx___, Index[name_, i_?NumericQ], yy___] := x[xx, i, yy];
189
190(* ************************** *)
191(* *****    Indices     ***** *)
192(* ************************** *)
193
194IndexRange[Index[SU3W      ]] = Unfold[Range[8]];
195IndexRange[Index[ASU3W     ]] = Unfold[Range[8]];
196IndexRange[Index[ASU3T     ]] = Unfold[Range[3]];
197IndexRange[Index[SU3T      ]] = Unfold[Range[3]];
198IndexRange[Index[Gluon     ]] = NoUnfold[Range[8]];
199IndexRange[Index[Colour    ]] = NoUnfold[Range[3]];
200IndexRange[Index[Generation1]] = Range[3,3];
201IndexRange[Index[Generation2]] = Range[2];
202IndexRange[Index[Generation]] = Range[3];
203
204IndexStyle[SU3W,       j];
205IndexStyle[ASU3W,      o];
206IndexStyle[ASU3T,      p];
207IndexStyle[SU3T,       k];
208IndexStyle[Gluon,      a];
209IndexStyle[Colour,     m];
210IndexStyle[Generation1, r];
211IndexStyle[Generation2, g];
212IndexStyle[Generation, f];
213
214
215(* ************************** *)
216(* *** Interaction orders *** *)
217(* ***  (as used by mg5)  *** *)
218(* ************************** *)
219
220M$InteractionOrderHierarchy = {
221  {QCD, 1},
222  {QED, 2}
223};
224
225
226(* ************************** *)
227(* **** Particle classes **** *)
228(* ************************** *)
229M$ClassesDescription = {
230
231(* Gauge bosons: physical vector fields *)
232  V[1] == {
233    ClassName       -> A,
234    SelfConjugate   -> True, 
235    Mass            -> 0, 
236    Width           -> 0, 
237    ParticleName    -> "a",
238    PDG             -> 22,
239    PropagatorLabel -> "a",
240    PropagatorType  -> W,
241    PropagatorArrow -> None,
242    FullName        -> "Photon"
243  },
244  V[2] == {
245    ClassName       -> Z,
246    SelfConjugate   -> True,
247    Mass            -> {MZ, 91.1876},
248    Width           -> {WZ, 2.4952},
249    ParticleName    -> "Z",
250    PDG             -> 23,
251    PropagatorLabel -> "Z",
252    PropagatorType  -> Sine,
253    PropagatorArrow -> None,
254    FullName        -> "Z"
255  },
256  V[3] == {
257    ClassName       -> ZP,
258    SelfConjugate   -> True,
259    Mass            -> {MZP, 4000},
260    Width           -> {WZP, 10},
261    ParticleName    -> "ZP", 
262    PropagatorLabel -> "ZP",
263    PropagatorType  -> Sine,
264    PropagatorArrow -> None,
265    FullName        -> "ZP"
266  },
267  V[4] == {
268    ClassName        -> W,
269    SelfConjugate    -> False,
270    Mass             -> {MW, 80.385},
271    Width            -> {WW, 2.085},
272    ParticleName     -> "W+",
273    AntiParticleName -> "W-",
274    QuantumNumbers   -> {Q -> 1},
275    PDG              -> 24,
276    PropagatorLabel  -> "W",
277    PropagatorType   -> Sine,
278    PropagatorArrow  -> Forward,
279    FullName         -> "W"
280  },
281  V[5] == {
282    ClassName        -> YY,
283    SelfConjugate    -> False,
284    Mass             -> {MY, Internal},
285    Width            -> {WY, 10},
286    ParticleName     -> "~Y-",
287    AntiParticleName -> "~Y+",
288    QuantumNumbers   -> {Q -> -1},
289    PropagatorLabel  -> "YY",
290    PropagatorType   -> Sine,
291    PropagatorArrow  -> Forward,
292    FullName         -> "YY"
293  },
294  V[6] == {
295    ClassName        -> V,
296    SelfConjugate    -> False,
297    Mass             -> {MV, Internal},
298    Width            -> {WV, 10},
299    ParticleName     -> "~V",
300    AntiParticleName -> "~V~",
301    QuantumNumbers   -> {Q -> 0},
302    PropagatorLabel  -> "V",
303    PropagatorType   -> Sine,
304    PropagatorArrow  -> None,
305    FullName         -> "V"
306  },
307  V[7] == {
308    ClassName        -> G,
309    SelfConjugate    -> True,
310    Indices          -> {Index[Gluon]},
311    Mass             -> 0,
312    Width            -> 0,
313    ParticleName     -> "g",
314    PDG              -> 21,
315    PropagatorLabel  -> "G",
316    PropagatorType   -> C,
317    PropagatorArrow  -> None,
318    FullName         -> "G"
319  },
320
321(* Ghosts: related to physical gauge bosons *)
322  U[1] == {
323    ClassName       -> ghA,
324    SelfConjugate   -> False,
325    Ghost           -> A,
326    QuantumNumbers  -> {GhostNumber -> 1},
327    Mass            -> 0,
328    Width           -> 0,
329    PropagatorLabel -> "uA",
330    PropagatorType  -> GhostDash,
331    PropagatorArrow -> Forward
332  },
333  U[2] == {
334    ClassName       -> ghZ,
335    SelfConjugate   -> False,
336    Ghost           -> Z,
337    QuantumNumbers  -> {GhostNumber -> 1},
338    Mass            -> {MZ,91.1876}, 
339    Width           -> {WZ, 2.4952},
340    PropagatorLabel -> "uZ",
341    PropagatorType  -> GhostDash,
342    PropagatorArrow -> Forward
343  },
344  U[3] == {
345    ClassName       -> ghZP,
346    SelfConjugate   -> False,
347    Ghost           -> ZP,
348    QuantumNumbers  -> {GhostNumber -> 1},
349    Mass            -> {MZP, 4000}, 
350    Width           -> {WZP, 10},
351    PropagatorLabel -> "uZP",
352    PropagatorType  -> GhostDash,
353    PropagatorArrow -> Forward
354  },
355  U[41] == {
356    ClassName       -> ghWp,
357    SelfConjugate   -> False,
358    Ghost           -> W,
359    QuantumNumbers  -> {GhostNumber -> 1, Q -> 1},
360    Mass            -> {MW, 80.385},
361    Width           -> {WW, 2.085},
362    PropagatorLabel -> "uWp",
363    PropagatorType  -> GhostDash,
364    PropagatorArrow -> Forward
365  },
366  U[42] == {
367    ClassName       -> ghWm,
368    SelfConjugate   -> False,
369    Ghost           -> Wbar,
370    QuantumNumbers  -> {GhostNumber -> 1, Q -> -1},
371    Mass            -> {MW, 80.385},
372    Width           -> {WW, 2.085},
373    PropagatorLabel -> "uWm",
374    PropagatorType  -> GhostDash,
375    PropagatorArrow -> Forward
376  },
377  U[51] == {
378    ClassName       -> ghYp,
379    SelfConjugate   -> False,
380    Ghost           -> YY,
381    QuantumNumbers  -> {GhostNumber -> 1, Q -> -1},
382    Mass            -> {MY, Internal},
383    Width           -> {WY, 10},
384    PropagatorLabel -> "uYp",
385    PropagatorType  -> GhostDash,
386    PropagatorArrow -> Forward
387  },
388  U[52] == {
389    ClassName       -> ghYm,
390    SelfConjugate   -> False,
391    Ghost           -> YYbar,
392    QuantumNumbers  -> {GhostNumber -> 1, Q -> 1},
393    Mass            -> {MY, Internal},
394    Width           -> {WY, 10},
395    PropagatorLabel -> "uYm",
396    PropagatorType  -> GhostDash,
397    PropagatorArrow -> Forward
398  },
399  U[61] == {
400    ClassName       -> ghVp,
401    SelfConjugate   -> False,
402    Ghost           -> V,
403    QuantumNumbers  -> {GhostNumber -> 1, Q -> 0},
404    Mass            -> {MV,Internal},
405    Width           -> {WV, 10},
406    PropagatorLabel -> "uV",
407    PropagatorType  -> GhostDash,
408    PropagatorArrow -> Forward
409  },
410  U[62] == {
411    ClassName       -> ghVm,
412    SelfConjugate   -> False,
413    Ghost           -> Vbar,
414    QuantumNumbers  -> {GhostNumber -> 1, Q -> 0},
415    Mass            -> {MV,Internal},
416    Width           -> {WV, 10},
417    PropagatorLabel -> "uV~",
418    PropagatorType  -> GhostDash,
419    PropagatorArrow -> Forward
420  },
421  U[7] == {
422    ClassName       -> ghG,
423    SelfConjugate   -> False,
424    Indices         -> {Index[Gluon]},
425    Ghost           -> G,
426    QuantumNumbers  ->{GhostNumber -> 1},
427    Mass            -> 0,
428    Width           -> 0,
429    PropagatorLabel -> "uG",
430    PropagatorType  -> GhostDash,
431    PropagatorArrow -> Forward
432  },
433
434(* Gauge bosons: unphysical vector fields *)
435  V[12] == {
436    ClassName     -> K,
437    Unphysical    -> True,
438    SelfConjugate -> True,
439    Definitions   -> { K[mu_] -> c3 (cz ZP[mu]-sz Z[mu]) + s3 (-sw (cz Z[mu]+sz ZP[mu]) + cw A[mu])}
440  },
441  V[13] == {
442    ClassName     -> Wi,
443    Unphysical    -> True,
444    SelfConjugate -> True,
445    Indices       -> {Index[SU3W]},
446    FlavorIndex   -> SU3W,
447    Definitions   -> {
448Wi[mu_,1] -> (Wbar[mu]+W[mu])/Sqrt[2], Wi[mu_,2] -> (Wbar[mu]-W[mu])/(I*Sqrt[2]),
449Wi[mu_,4] -> (YYbar[mu]+YY[mu])/Sqrt[2], Wi[mu_,5] -> (YY[mu]-YYbar[mu])/(I*Sqrt[2]),
450Wi[mu_,6] -> (Vbar[mu]+V[mu])/Sqrt[2], Wi[mu_,7] -> (V[mu]-Vbar[mu])/(I*Sqrt[2]),
451Wi[mu_,3] -> cw (cz Z[mu]+sz ZP[mu]) + sw A[mu], Wi[mu_,8] -> -s3 (cz ZP[mu]-sz Z[mu]) + c3 (-sw (cz Z[mu]+sz ZP[mu]) +cw A[mu])}
452  },
453  V[14] == {
454    ClassName     -> WWi,
455    Unphysical    -> True,
456    SelfConjugate -> True,
457    Indices       -> {Index[ASU3W]},
458    FlavorIndex   -> ASU3W,
459    Definitions   -> {
460WWi[mu_,1] -> (Wbar[mu]+W[mu])/Sqrt[2], WWi[mu_,2] -> (Wbar[mu]-W[mu])/(I*Sqrt[2]),
461WWi[mu_,4] -> (YYbar[mu]+YY[mu])/Sqrt[2], WWi[mu_,5] -> (YY[mu]-YYbar[mu])/(I*Sqrt[2]),
462WWi[mu_,6] -> (Vbar[mu]+V[mu])/Sqrt[2], WWi[mu_,7] -> (V[mu]-Vbar[mu])/(I*Sqrt[2]),
463WWi[mu_,3] -> cw (cz Z[mu]+sz ZP[mu]) + sw A[mu], WWi[mu_,8] -> -s3 (cz ZP[mu]-sz Z[mu])+c3 (-sw (cz Z[mu]+sz ZP[mu])+cw A[mu])}
464  },
465
466(* Ghosts: related to unphysical gauge bosons *)
467  U[12] == {
468    ClassName     -> ghK,
469    Unphysical    -> True,
470    SelfConjugate -> False,
471    Ghost         -> K,
472    Definitions   -> { ghK -> -c3 (cz ghZP-sz ghZ) + s3 (-sw (cz ghZ+sz ghZP) + cw ghA)}
473  },
474  U[13] == {
475    ClassName     -> ghWi,
476    Unphysical    -> True,
477    SelfConjugate -> False,
478    Ghost         -> Wi,
479    Indices       -> {Index[SU3W]},
480    FlavorIndex   -> SU3W,
481    Definitions   -> { ghWi[1] -> (ghWp+ghWm)/Sqrt[2], ghWi[2] -> (ghWm-ghWp)/(I*Sqrt[2]), ghWi[4] -> (ghYp+ghYm)/Sqrt[2], ghWi[5] -> (ghYp-ghYm)/(I*Sqrt[2]), ghWi[6] -> (ghVp+ghVm)/Sqrt[2], ghWi[7] -> (ghVp-ghVm)/(I*Sqrt[2]),  ghWi[3] -> cw (cz ghZ+sz ghZP)+sw ghA, ghWi[8] ->  -s3 (cz ghZP-sz ghZ)+c3 (-sw (cz ghZ+sz ghZP) + cw ghA)}
482  },
483
484(* Fermions: physical fields *)
485  F[1] == {
486    ClassName        -> vl,
487    ClassMembers     -> {ve,vm,vt},
488    Indices          -> {Index[Generation]},
489    FlavorIndex      -> Generation,
490    SelfConjugate    -> False,
491    Mass             -> 0,
492    Width            -> 0,
493    QuantumNumbers   -> {LeptonNumber -> 1},
494    PropagatorLabel  -> {"v", "ve", "vm", "vt"} ,
495    PropagatorType   -> S,
496    PropagatorArrow  -> Forward,
497    PDG              -> {12,14,16},
498    ParticleName     -> {"ve","vm","vt"},
499    AntiParticleName -> {"ve~","vm~","vt~"},
500    FullName         -> {"Electron-neutrino", "Mu-neutrino", "Tau-neutrino"}
501  },
502  F[2] == {
503    ClassName        -> l,
504    ClassMembers     -> {e, mu, ta},
505    Indices          -> {Index[Generation]},
506    FlavorIndex      -> Generation,
507    SelfConjugate    -> False,
508    Mass             -> {Ml, {Me,5.11*^-4}, {MMU,0.10566}, {MTA,1.777}},
509    Width            -> 0,
510    QuantumNumbers   -> {Q -> -1, LeptonNumber -> 1},
511    PropagatorLabel  -> {"l", "e", "mu", "ta"},
512    PropagatorType   -> Straight,
513    PropagatorArrow  -> Forward,
514    PDG              -> {11, 13, 15},
515    ParticleName     -> {"e-", "mu-", "ta-"},
516    AntiParticleName -> {"e+", "mu+", "ta+"},
517    FullName         -> {"Electron", "Muon", "Tau"}
518  },
519  F[3] == {
520    ClassName        -> EE,
521    ClassMembers     -> {Ee, Emu, Eta},
522    Indices          -> {Index[Generation]},
523    FlavorIndex      -> Generation,
524    SelfConjugate    -> False,
525    Mass             -> {ME, {MEE,1*^3}, {MEMU,1*^3}, {META,1*^3}},
526    Width            -> {WE, {WEE,10},{WEMU,10},{WETA,10}},
527    QuantumNumbers   -> {LeptonNumber -> 1},
528    PropagatorLabel  -> {"E", "Ee", "Emu", "Eta"},
529    PropagatorType   -> Straight,
530    PropagatorArrow  -> Forward,
531    ParticleName     -> {"~Ee", "~Emu", "~Eta"},
532    AntiParticleName -> {"~Ee~", "~Emu~", "~Eta~"},
533    FullName         -> {"HElectron", "HMuon", "HTau"}
534  },
535  F[4] == {
536    ClassName        -> uq,
537    ClassMembers     -> {u, c, t},
538    Indices          -> {Index[Generation], Index[Colour]},
539    FlavorIndex      -> Generation,
540    SelfConjugate    -> False,
541    Mass             -> {Mu, {MU, 2.55*^-3}, {MC,1.27}, {MT,172}},
542    Width            -> {0, 0, {WT,1.50833649}},
543    QuantumNumbers   -> {Q -> 2/3},
544    PropagatorLabel  -> {"uq", "u", "c", "t"},
545    PropagatorType   -> Straight,
546    PropagatorArrow  -> Forward,
547    PDG              -> {2, 4, 6},
548    ParticleName     -> {"u",  "c",  "t" },
549    AntiParticleName -> {"u~", "c~", "t~"},
550    FullName         -> {"u-quark", "c-quark", "t-quark"}
551  },
552  F[5] == {
553    ClassName        -> dq,
554    ClassMembers     -> {d, s, b},
555    Indices          -> {Index[Generation], Index[Colour]},
556    FlavorIndex      -> Generation,
557    SelfConjugate    -> False,
558    Mass             -> {Md, {MD,5.04*^-3}, {MS,0.101}, {MB,4.7}},
559    Width            -> 0,
560    QuantumNumbers   -> {Q -> -1/3},
561    PropagatorLabel  -> {"dq", "d", "s", "b"},
562    PropagatorType   -> Straight,
563    PropagatorArrow  -> Forward,
564    PDG              -> {1,3,5},
565    ParticleName     -> {"d",  "s",  "b" },
566    AntiParticleName -> {"d~", "s~", "b~"},
567    FullName         -> {"d-quark", "s-quark", "b-quark"}
568  },
569   F[6] == {
570    ClassName        -> Jq12,
571    ClassMembers     -> {Jd, Js},
572    Indices          -> {Index[Generation2], Index[Colour]},
573    FlavorIndex      -> Generation2,
574    SelfConjugate    -> False,
575    Mass             -> {MJ12, {MJD,1*^3}, {MJS,1*^3}},
576    Width            -> {WJQ12,{WJD,10},{WJS,10}},
577    QuantumNumbers   -> {Q -> -1/3},
578    PropagatorLabel  -> {"Jq12", "Jd", "Js"},
579    PropagatorType   -> Straight,
580    PropagatorArrow  -> Forward,
581    ParticleName     -> {"~Jd",  "~Js"},
582    AntiParticleName -> {"~Jd~", "~Js~"},
583    FullName         -> {"Jd-quark", "Js-quark"}
584  },
585  F[7] == {
586    ClassName        -> Jt,
587    Indices          -> {Index[Colour]},
588    SelfConjugate    -> False,
589    Mass              -> {MJT,1*^3},
590    Width            -> 10,
591    QuantumNumbers   -> {Q -> 2/3},
592    PropagatorLabel  -> "Jt",
593    PropagatorType   -> Straight,
594    PropagatorArrow  -> Forward,
595    ParticleName     -> "~Jt",
596    AntiParticleName -> "~Jt~",
597    FullName         -> "Jt-quark"
598  },
599
600
601
602(* Fermions: unphysical fields *)
603   F[11] == {
604    ClassName      -> LL,
605    Unphysical     -> True,
606    Indices        -> {Index[ASU3T], Index[Generation]},
607    FlavorIndex    -> ASU3T,
608    SelfConjugate  -> False,
609    QuantumNumbers -> {X -> -1/3},
610    Definitions    -> {
611      LL[sp1_,1,ff_] :> Module[{sp2}, ProjM[sp1,sp2] l[sp2,ff]],
612      LL[sp1_,2,ff_] :> -Module[{sp2}, ProjM[sp1,sp2] vl[sp2,ff]],
613      LL[sp1_,3,ff_] :> Module[{sp2}, ProjM[sp1,sp2] EE[sp2,ff]]}
614  },
615  F[12] == {
616    ClassName      -> lR,
617    Unphysical     -> True,
618    Indices        -> {Index[Generation]},
619    FlavorIndex    -> Generation,
620    SelfConjugate  -> False,
621    QuantumNumbers -> {X -> -1},
622    Definitions    -> { lR[sp1_,ff_] :> Module[{sp2}, ProjP[sp1,sp2] l[sp2,ff]] }
623  },
624  F[13] == {
625    ClassName      -> EER,
626    Unphysical     -> True,
627    Indices        -> {Index[Generation]},
628    FlavorIndex    -> Generation,
629    SelfConjugate  -> False,
630    QuantumNumbers -> {X -> 0},
631    Definitions    -> { EER[sp1_,ff_] :> Module[{sp2}, ProjP[sp1,sp2] EE[sp2,ff]] }
632  },
633  F[14] == {
634    ClassName      -> QL12,
635    Unphysical     -> True,
636    Indices        -> {Index[SU3T], Index[Generation2], Index[Colour]},
637    FlavorIndex    -> SU3T,
638    SelfConjugate  -> False,
639    QuantumNumbers -> {X -> 0},
640    Definitions    -> {
641      QL12[sp1_,1,1,cc_] :> Module[{sp2,ff2}, RU[1,ff2]  ProjM[sp1,sp2] uq[sp2,ff2,cc]],
642      QL12[sp1_,1,2,cc_] :> Module[{sp2,ff2}, RU[2,ff2]  ProjM[sp1,sp2] uq[sp2,ff2,cc]],
643      QL12[sp1_,2,1,cc_] :> Module[{sp2,ff2}, CKM[1,ff2] ProjM[sp1,sp2] dq[sp2,ff2,cc]],
644      QL12[sp1_,2,2,cc_] :> Module[{sp2,ff2}, CKM[2,ff2] ProjM[sp1,sp2] dq[sp2,ff2,cc]],
645      QL12[sp1_,3,ff12_,cc_] :> Module[{sp2}, ProjM[sp1,sp2] Jq12[sp2,ff12,cc]]}
646  },
647 F[15] == {
648    ClassName      -> QL3,
649    Unphysical     -> True,
650    Indices        -> {Index[ASU3T], Index[Generation1], Index[Colour]},
651    FlavorIndex    -> ASU3T,
652    SelfConjugate  -> False,
653    QuantumNumbers -> {X -> 1/3},
654    Definitions    -> {
655      QL3[sp1_,1,3,cc_] :> Module[{sp2,ff2}, CKM[3,ff2] ProjM[sp1,sp2] dq[sp2,ff2,cc]],
656      QL3[sp1_,2,3,cc_] :>-Module[{sp2,ff2}, RU[3,ff2]  ProjM[sp1,sp2] uq[sp2,ff2,cc]],
657      QL3[sp1_,3,3,cc_] :> Module[{sp2}, ProjM[sp1,sp2] Jt[sp2,cc]]}
658  },
659  F[16] == {
660    ClassName      -> QL,
661    Unphysical     -> True,
662    Indices        -> {Index[SU3T], Index[Generation], Index[Colour]},
663    FlavorIndex    -> SU3T,
664    SelfConjugate  -> False,
665    QuantumNumbers -> {X -> 0, X -> 0, X -> 1/3},
666    Definitions    -> {
667      QL[sp1_,1,1,cc_] :> QL12[sp1,1,1,cc],
668      QL[sp1_,1,2,cc_] :> QL12[sp1,1,2,cc],
669      QL[sp1_,1,3,cc_] :> QL3[sp1,1,3,cc],
670      QL[sp1_,2,1,cc_] :> QL12[sp1,2,1,cc],
671      QL[sp1_,2,2,cc_] :> QL12[sp1,2,2,cc],
672      QL[sp1_,2,3,cc_] :> QL3[sp1,2,3,cc],
673      QL[sp1_,3,1,cc_] :> QL12[sp1,3,1,cc],
674      QL[sp1_,3,2,cc_] :> QL12[sp1,3,2,cc],
675      QL[sp1_,3,3,cc_] :> QL3[sp1,3,3,cc]}
676  },
677  F[17] == {
678    ClassName      -> uR,
679    Unphysical     -> True,
680    Indices        -> {Index[Generation], Index[Colour]},
681    FlavorIndex    -> Generation,
682    SelfConjugate  -> False,
683    QuantumNumbers -> {X -> 2/3},
684    Definitions    -> { uR[sp1_,ff_,cc_] :> Module[{sp2}, ProjP[sp1,sp2] uq[sp2,ff,cc]] }
685  },
686  F[18] == {
687    ClassName      -> dR,
688    Unphysical     -> True,
689    Indices        -> {Index[Generation], Index[Colour]},
690    FlavorIndex    -> Generation,
691    SelfConjugate  -> False,
692    QuantumNumbers -> {X -> -1/3},
693    Definitions    -> { dR[sp1_,ff_,cc_] :> Module[{sp2}, ProjP[sp1,sp2] dq[sp2,ff,cc]] }
694  },
695   F[19] == {
696    ClassName      -> JR12,
697    Unphysical     -> True,
698    Indices        -> {Index[Generation2], Index[Colour]},
699    FlavorIndex    -> Generation2,
700    SelfConjugate  -> False,
701    QuantumNumbers -> {X -> -1/3},
702    Definitions    -> { JR12[sp1_,ff_,cc_] :> Module[{sp2}, ProjP[sp1,sp2] Jq12[sp2,ff,cc]] }
703  },
704   F[20] == {
705    ClassName      -> JR3,
706    Unphysical     -> True,
707    Indices        -> {Index[Generation1], Index[Colour]},
708    FlavorIndex    -> Generation1,
709    SelfConjugate  -> False,
710    QuantumNumbers -> {X -> 2/3},
711    Definitions    -> { JR3[sp1_,ff_,cc_] :> Module[{sp2}, ProjP[sp1,sp2] Jt[sp2,cc]] }
712  },
713   F[21] == {
714    ClassName      -> JR,
715    Unphysical     -> True,
716    Indices        -> {Index[Generation], Index[Colour]},
717    FlavorIndex    -> Generation,
718    SelfConjugate  -> False,
719    QuantumNumbers -> {X -> -1/3, X -> -1/3, X -> 2/3},
720    Definitions    -> { JR[sp1_,1,cc_] :> JR12[sp1,1,cc],
721                        JR[sp1_,2,cc_] :> JR12[sp1,2,cc],
722                        JR[sp1_,3,cc_] :> JR3[sp1,3,cc] }
723  },
724 
725(* Higgs: physical scalars  *)
726  S[1] == {
727    ClassName       -> h,
728    SelfConjugate   -> True,
729    Mass            -> {Mh,125},
730    Width           -> {Wh,0.00407},
731    PropagatorLabel -> "h",
732    PropagatorType  -> D,
733    PropagatorArrow -> None,
734    PDG             -> 25,
735    ParticleName    -> "h",
736    FullName        -> "h"
737  },
738  S[2] == {
739    ClassName       -> H2,
740    SelfConjugate   -> True,
741    Mass            -> {MH2,Internal},
742    Width           -> {WH2,10},
743    PropagatorLabel -> "H2",
744    PropagatorType  -> D,
745    PropagatorArrow -> None,
746    ParticleName    -> "H2",
747    FullName        -> "H2"
748  },
749  S[3] == {
750    ClassName       -> H3,
751    SelfConjugate   -> True,
752    Mass            -> {MH3,Internal},
753    Width           -> {WH3,10},
754    PropagatorLabel -> "H3",
755    PropagatorType  -> D,
756    PropagatorArrow -> None,
757    ParticleName    -> "H3",
758    FullName        -> "H3"
759  },
760  S[4] == {
761    ClassName       -> H0,
762    SelfConjugate   -> True,
763    Mass            -> {MH0,Internal},
764    Width           -> {WH0,10},
765    PropagatorLabel -> "H0",
766    PropagatorType  -> D,
767    PropagatorArrow -> None,
768    ParticleName    -> "H0",
769    FullName        -> "H0"
770  },
771  S[5] == {
772    ClassName       -> HW,
773    SelfConjugate   -> False,
774    Mass            -> {MHW,Internal},
775    Width           -> {WHW,10},
776    ParticleName    -> "HW+",
777    AntiParticleName -> "HW-",
778    QuantumNumbers  -> {Q -> 1},
779    PropagatorLabel -> "HW",
780    PropagatorType  -> D,
781    PropagatorArrow -> Forward,
782    FullName        -> "HW"
783  },
784  S[6] == {
785    ClassName       -> HY,
786    SelfConjugate   -> False,
787    Mass            -> {MHY,Internal},
788    Width           -> {WHY,10},
789    ParticleName    -> "~HY-",
790    AntiParticleName -> "~HY+",
791    QuantumNumbers  -> {Q -> -1},
792    PropagatorLabel -> "HY",
793    PropagatorType  -> D,
794    PropagatorArrow -> Forward,
795    FullName        -> "HY"
796  },
797  S[71] == {
798    ClassName       -> HV,
799    SelfConjugate   -> False,
800    Mass            -> {MHV,Internal},
801    Width           -> {WHV,10},
802    ParticleName    -> "~HV",
803    AntiParticleName -> "~HV~",
804    QuantumNumbers  -> {Q -> 0},
805    PropagatorLabel -> "HV",
806    PropagatorType  -> D,
807    PropagatorArrow -> None,
808    FullName        -> "HV"
809  },
810
811
812(* Higgs: physical scalars  *)
813  S[8] == {
814    ClassName       -> GZ,
815    SelfConjugate   -> True,
816    Goldstone       -> Z,
817    Mass            -> {MZ, 91.1876},
818    Width           -> {WZ, 2.4952},
819    PropagatorLabel -> "GZ",
820    PropagatorType  -> D,
821    PropagatorArrow -> None,
822    PDG             -> 250,
823    ParticleName    -> "GZ",
824    FullName        -> "GZ"
825  },
826  S[9] == {
827    ClassName       -> GZP,
828    SelfConjugate   -> True,
829    Goldstone       -> ZP,
830    Mass            -> {MZP, 4000},
831    Width           -> {WZP, 10},
832    PropagatorLabel -> "GZP",
833    PropagatorType  -> D,
834    PropagatorArrow -> None,
835    ParticleName    -> "GZP",
836    FullName        -> "GZP"
837  },
838  S[10] == {
839    ClassName        -> GW,
840    SelfConjugate    -> False,
841    Goldstone        -> W,
842    Mass             -> {MW, 80.385},
843    QuantumNumbers   -> {Q -> 1},
844    Width            -> {WW, 2.085},
845    PropagatorLabel  -> "GW",
846    PropagatorType   -> D,
847    PropagatorArrow  -> None,
848    PDG              -> 251,
849    ParticleName     -> "GW+",
850    AntiParticleName -> "GW-",
851    FullName         -> "GW"
852  },
853   S[11] == {
854    ClassName        -> GY,
855    SelfConjugate    -> False,
856    Goldstone        -> YY,
857    Mass             -> {MY, Internal},
858    QuantumNumbers   -> {Q -> -1},
859    Width            -> {WY, 10},
860    PropagatorLabel  -> "GY",
861    PropagatorType   -> D,
862    PropagatorArrow  -> None,
863    ParticleName     -> "GY-",
864    AntiParticleName -> "GY+",
865    FullName         -> "GY"
866  },
867   S[12] == {
868    ClassName        -> GV,
869    SelfConjugate    -> False,
870    Goldstone        -> V,
871    Mass             -> {MV, Internal},
872    QuantumNumbers  -> {Q -> 0},
873    Width            -> {WV, 10},
874    PropagatorLabel  -> "GV",
875    PropagatorType   -> D,
876    PropagatorArrow  -> None,
877    ParticleName     -> "GV",
878    AntiParticleName -> "GV~",
879    FullName         -> "GV"
880  },
881
882
883(* Higgs: unphysical scalars  *)
884  S[13] == {
885    ClassName      -> Rho,
886    Unphysical     -> True,
887    Indices        -> {Index[SU3T]},
888    FlavorIndex    -> SU3T,
889    SelfConjugate  -> False,
890    QuantumNumbers -> {X -> 1/3},
891    Definitions    -> { Rho[1] -> -I (HW svv2+GW cvv2), Rho[2] -> (v + UH11 h + UH12 H2 + UH13 H3 + I (Uh11 H0 + Uh12 GZ + Uh13 GZP))/Sqrt[2], Rho[3] -> -I (HV svv3+GV cvv3)}
892  },
893  S[14] == {
894    ClassName      -> Phi,
895    Unphysical     -> True,
896    Indices        -> {Index[SU3T]},
897    FlavorIndex    -> SU3T,
898    SelfConjugate  -> False,
899    QuantumNumbers -> {X -> -2/3},
900    Definitions    -> { Phi[1] -> (v2 + UH21 h + UH22 H2 + UH23 H3 + I (Uh21 H0  + Uh23 GZP))/Sqrt[2], Phi[2] -> -I(GWbar svv2 -HWbar cvv2), Phi[3] -> -I(HY sv2v3+GY cv2v3) }
901  },
902   S[15] == {
903    ClassName      -> Chi,
904    Unphysical     -> True,
905    Indices        -> {Index[SU3T]},
906    FlavorIndex    -> SU3T,
907    SelfConjugate  -> False,
908    QuantumNumbers -> {X -> 1/3},
909    Definitions    -> { Chi[1] -> -I (GYbar sv2v3-HYbar cv2v3), Chi[2] -> -I (GVbar svv3-HVbar cvv3), Chi[3] -> (v3 + UH31 h + UH32 H2 + UH33 H3 + I (Uh31 H0 + Uh32 GZ + Uh33 GZP))/Sqrt[2] }
910  },
911   S[16] == {
912    ClassName      -> Su,
913    Unphysical     -> True,
914    Indices        -> {Index[Generation], Index[Generation],Index[SU3T]},
915    FlavorIndex    -> SU3T,
916    SelfConjugate  -> False,
917    Definitions    -> {Su[1,1,kk_] -> Phi[kk],  Su[1,2,kk_] -> 0, Su[1,3,kk_] -> 0, Su[2,1,kk_] -> 0, Su[2,2,kk_] -> Phi[kk], Su[2,3,kk_] -> 0, Su[3,1,kk_] -> 0, Su[3,2,kk_] -> 0, Su[3,3,kk_] -> -Rhobar[kk]}
918  },
919   S[17] == {
920    ClassName      -> Sd,
921    Unphysical     -> True,
922    Indices        -> {Index[Generation], Index[Generation],Index[SU3T]},
923    FlavorIndex    -> SU3T,
924    SelfConjugate  -> False,
925    Definitions    -> {Sd[1,1,kk_] -> Rho[kk],  Sd[1,2,kk_] -> 0, Sd[1,3,kk_] -> 0, Sd[2,1,kk_] -> 0, Sd[2,2,kk_] -> Rho[kk], Sd[2,3,kk_] -> 0, Sd[3,1,kk_] -> 0, Sd[3,2,kk_] -> 0, Sd[3,3,kk_] -> Phibar[kk]}
926  },
927   S[18] == {
928    ClassName      -> SJ,
929    Unphysical     -> True,
930    Indices        -> {Index[Generation], Index[Generation],Index[SU3T]},
931    FlavorIndex    -> SU3T,
932    SelfConjugate  -> False,
933    Definitions    -> {SJ[1,1,kk_] -> Chi[kk],  SJ[1,2,kk_] -> 0, SJ[1,3,kk_] -> 0, SJ[2,1,kk_] -> 0, SJ[2,2,kk_] -> Chi[kk], SJ[2,3,kk_] -> 0, SJ[3,1,kk_] -> 0, SJ[3,2,kk_] -> 0, SJ[3,3,kk_] -> Chibar[kk]}
934  }
935};
936 
937
938(* ************************** *)
939(* *****     Gauge      ***** *)
940(* *****   Parameters   ***** *)
941(* *****   (FeynArts)   ***** *)
942(* ************************** *)
943
944GaugeXi[ V[1]  ] = GaugeXi[A];
945GaugeXi[ V[2]  ] = GaugeXi[Z];
946GaugeXi[ V[3]  ] = GaugeXi[ZP];
947GaugeXi[ V[4]  ] = GaugeXi[W];
948GaugeXi[ V[5]  ] = GaugeXi[YY];
949GaugeXi[ V[6]  ] = GaugeXi[V];
950GaugeXi[ V[7]  ] = GaugeXi[G];
951GaugeXi[ S[1]  ] = 1;
952GaugeXi[ S[2]  ] = 1;
953GaugeXi[ S[3]  ] = 1;
954GaugeXi[ S[4]  ] = 1;
955GaugeXi[ S[5]  ] = 1;
956GaugeXi[ S[6]  ] = 1;
957GaugeXi[ S[7]  ] = 1;
958GaugeXi[ S[8]  ] = GaugeXi[Z];
959GaugeXi[ S[9]  ] = GaugeXi[ZP];
960GaugeXi[ S[10]  ] = GaugeXi[W];
961GaugeXi[ S[11]  ] = GaugeXi[YY];
962GaugeXi[ S[12]  ] = GaugeXi[V];
963GaugeXi[ U[1]  ] = GaugeXi[A];
964GaugeXi[ U[2]  ] = GaugeXi[Z];
965GaugeXi[ U[3]  ] = GaugeXi[ZP];
966GaugeXi[ U[41] ] = GaugeXi[W];
967GaugeXi[ U[42] ] = GaugeXi[W];
968GaugeXi[ U[51] ] = GaugeXi[YY];
969GaugeXi[ U[52] ] = GaugeXi[YY];
970GaugeXi[ U[61] ] = GaugeXi[V];
971GaugeXi[ U[62] ] = GaugeXi[V];
972GaugeXi[ U[7]  ] = GaugeXi[G];
973
974
975(* ************************** *)
976(* *****   Parameters   ***** *)
977(* ************************** *)
978M$Parameters = {
979
980  (* External parameters *)
981  aEWM1 == {
982    ParameterType    -> External,
983    BlockName        -> SMINPUTS,
984    OrderBlock       -> 1,
985    Value            -> 127.9,
986    InteractionOrder -> {QED,-2},
987    Description      -> "Inverse of the EW coupling constant at the Z pole"
988  },
989  Gf == {
990    ParameterType    -> External,
991    BlockName        -> SMINPUTS,
992    OrderBlock       -> 2,
993    Value            -> 1.16637*^-5,
994    InteractionOrder -> {QED,2},
995    TeX              -> Subscript[G,f],
996    Description      -> "Fermi constant"
997  },
998  aS    == {
999    ParameterType    -> External,
1000    BlockName        -> SMINPUTS,
1001    OrderBlock       -> 3,
1002    Value            -> 0.1184,
1003    InteractionOrder -> {QCD,2},
1004    TeX              -> Subscript[\[Alpha],s],
1005    Description      -> "Strong coupling constant at the Z pole"
1006  },
1007  ymdo == {
1008    ParameterType -> External,
1009    BlockName     -> YUKAWA,
1010    OrderBlock    -> 1,
1011    Value         -> 5.04*^-3,
1012    Description   -> "Down Yukawa mass"
1013  },
1014  ymup == {
1015    ParameterType -> External,
1016    BlockName     -> YUKAWA,
1017    OrderBlock    -> 2,
1018    Value         -> 2.55*^-3,
1019    Description   -> "Up Yukawa mass"
1020  },
1021  yms == {
1022    ParameterType -> External,
1023    BlockName     -> YUKAWA,
1024    OrderBlock    -> 3,
1025    Value         -> 0.101,
1026    Description   -> "Strange Yukawa mass"
1027  },
1028  ymc == {
1029    ParameterType -> External,
1030    BlockName     -> YUKAWA,
1031    OrderBlock    -> 4,
1032    Value         -> 1.27,
1033    Description   -> "Charm Yukawa mass"
1034  },
1035  ymb == {
1036    ParameterType -> External,
1037    BlockName     -> YUKAWA,
1038    OrderBlock    -> 5,
1039    Value         -> 4.7,
1040    Description   -> "Bottom Yukawa mass"
1041  },
1042  ymt == {
1043    ParameterType -> External,
1044    BlockName     -> YUKAWA,
1045    OrderBlock    -> 6,
1046    Value         -> 172,
1047    Description   -> "Top Yukawa mass"
1048  },
1049  ymD == {
1050    ParameterType -> External,
1051    BlockName     -> YUKAWA,
1052    OrderBlock    -> 7,
1053    Value         -> 1*^3,
1054    Description   -> "Heavy Down Yukawa mass"
1055  },
1056  ymS == {
1057    ParameterType -> External,
1058    BlockName     -> YUKAWA,
1059    OrderBlock    -> 8,
1060    Value         -> 1*^3,
1061    Description   -> "Heavy Strange Yukawa mass"
1062  },
1063  ymT == {
1064    ParameterType -> External,
1065    BlockName     -> YUKAWA,
1066    OrderBlock    -> 9,
1067    Value         -> 1*^3,
1068    Description   -> "Heavy Top Yukawa mass"
1069  },
1070  yme == {
1071    ParameterType -> External,
1072    BlockName     -> YUKAWA,
1073    OrderBlock    -> 11,
1074    Value         -> 5.11*^-4,
1075    Description   -> "Electron Yukawa mass"
1076  },
1077  ymm == {
1078    ParameterType -> External,
1079    BlockName     -> YUKAWA,
1080    OrderBlock    -> 13,
1081    Value         -> 0.10566,
1082    Description   -> "Muon Yukawa mass"
1083  },
1084  ymtau == {
1085    ParameterType -> External,
1086    BlockName     -> YUKAWA,
1087    OrderBlock    -> 15,
1088    Value         -> 1.777,
1089    Description   -> "Tau Yukawa mass"
1090  },
1091  ymEe == {
1092    ParameterType -> External,
1093    BlockName     -> YUKAWA,
1094    OrderBlock    -> 16,
1095    Value         -> 1*^3,
1096    Description   -> "Heavy Electron Yukawa mass"
1097  },
1098  ymEm == {
1099    ParameterType -> External,
1100    BlockName     -> YUKAWA,
1101    OrderBlock    -> 17,
1102    Value         -> 1*^3,
1103    Description   -> "Heavy Muon Yukawa mass"
1104  },
1105  ymEtau == {
1106    ParameterType -> External,
1107    BlockName     -> YUKAWA,
1108    OrderBlock    -> 18,
1109    Value         -> 1*^3,
1110    Description   -> "Heavy Tau Yukawa mass"
1111  },
1112  cabi == {
1113    ParameterType -> External,
1114    BlockName     -> CKMBLOCK,
1115    OrderBlock    -> 1,
1116    Value         -> 0.227736,
1117    TeX           -> Subscript[\[Theta], c],
1118    Description   -> "Cabibbo angle"
1119  },
1120  v2 == {
1121    ParameterType    -> External,
1122    BlockName     -> OTHERS,
1123    OrderBlock    -> 21,
1124    Value            -> 174,
1125    InteractionOrder -> {QED,-1},
1126    Description      -> "phi vacuum expectation value"
1127  },
1128  lam2 == {
1129    ParameterType    -> External,
1130    BlockName     -> OTHERS,
1131    OrderBlock    -> 2,
1132    Value            -> -1,
1133    InteractionOrder -> {QED, 2},
1134    TeX              -> Subscript[\[Lambda], 2],
1135    Description      -> "phi quartic coupling"
1136  },
1137  lam3 == {
1138    ParameterType    -> External,
1139    BlockName     -> OTHERS,
1140    OrderBlock    -> 3,
1141    Value            -> -1,
1142    InteractionOrder -> {QED, 2},
1143    TeX              -> Subscript[\[Lambda], 3],
1144    Description      -> "Chi quartic coupling"
1145  },
1146  lam12 == {
1147    ParameterType    -> External,
1148    BlockName     -> OTHERS,
1149    OrderBlock    -> 4,
1150    Value            -> -1,
1151    InteractionOrder -> {QED, 2},
1152    TeX              -> Subscript[\[Lambda], 12],
1153    Description      -> "Rho Rho and phi phi quartic coupling"
1154  },
1155  lam13 == {
1156    ParameterType    -> External,
1157    BlockName     -> OTHERS,
1158    OrderBlock    -> 5,
1159    Value            -> -1,
1160    InteractionOrder -> {QED, 2},
1161    TeX              -> Subscript[\[Lambda], 13],
1162    Description      -> "Rho Rho and Chi Chi quartic coupling"
1163  },
1164  lam23 == {
1165    ParameterType    -> External,
1166    BlockName     -> OTHERS,
1167    OrderBlock    -> 6,
1168    Value            -> -1,
1169    InteractionOrder -> {QED, 2},
1170    TeX              -> Subscript[\[Lambda], 23],
1171    Description      -> "Chi Chi and phi phi quartic coupling"
1172  },
1173  lam12P == {
1174    ParameterType    -> External,
1175    BlockName     -> OTHERS,
1176    OrderBlock    -> 7,
1177    Value            -> -1,
1178    InteractionOrder -> {QED, 2},
1179    TeX              -> Subscript[\[Lambda]', 12],
1180    Description      -> "Rho phi and phi Rho quartic coupling"
1181  },
1182  lam13P == {
1183    ParameterType    -> External,
1184    BlockName     -> OTHERS,
1185    OrderBlock    -> 8,
1186    Value            -> -1,
1187    InteractionOrder -> {QED, 2},
1188    TeX              -> Subscript[\[Lambda]', 13],
1189    Description      -> "Rho Chi and Chi Rho quartic coupling"
1190  },
1191  lam23P == {
1192    ParameterType    -> External,
1193    BlockName     -> OTHERS,
1194    OrderBlock    -> 9,
1195    Value            -> -1,
1196    InteractionOrder -> {QED, 2},
1197    TeX              -> Subscript[\[Lambda]', 23],
1198    Description      -> "phi Chi and Chi phi quartic coupling"
1199  },
1200  UH11 == {
1201    ParameterType -> External,
1202    BlockName     -> OTHERS,
1203    OrderBlock    -> 10,
1204    Value         -> -0.716449,
1205    TeX           -> Subsuperscript[U,11,H],
1206    Description   -> "11 rotation-matrix element of Realscalar"
1207  },
1208  UH12 == {
1209    ParameterType -> External,
1210    BlockName     -> OTHERS,
1211    OrderBlock    -> 11,
1212    Value         -> 0.69729,
1213    TeX           -> Subsuperscript[U,12,H],
1214    Description   -> "12 rotation-matrix element of Realscalar"
1215  },
1216  UH13 == {
1217    ParameterType -> External,
1218    BlockName     -> OTHERS,
1219    OrderBlock    -> 12,
1220    Value         -> 0.0220899,
1221    TeX           -> Subsuperscript[U,13,H],
1222    Description   -> "13 rotation-matrix element of Realscalar"
1223  },
1224  UH21 == {
1225    ParameterType -> External,
1226    BlockName     -> OTHERS,
1227    OrderBlock    -> 13,
1228    Value         -> -0.69693,
1229    TeX           -> Subsuperscript[U,21,H],
1230    Description   -> "21 rotation-matrix element of Realscalar"
1231  },
1232  UH22 == {
1233    ParameterType -> External,
1234    BlockName     -> OTHERS,
1235    OrderBlock    -> 14,
1236    Value         -> -0.716789,
1237    TeX           -> Subsuperscript[U,22,H],
1238    Description   -> "22 rotation-matrix element of Realscalar"
1239  },
1240  UH23 == {
1241    ParameterType -> External,
1242    BlockName     -> OTHERS,
1243    OrderBlock    -> 15,
1244    Value         -> 0.0224204,
1245    TeX           -> Subsuperscript[U,23,H],
1246    Description   -> "23 rotation-matrix element of Realscalar"
1247  },
1248  UH31 == {
1249    ParameterType -> External,
1250    BlockName     -> OTHERS,
1251    OrderBlock    -> 16,
1252    Value         -> 0.0314673,
1253    TeX           -> Subsuperscript[U,31,H],
1254    Description   -> "31 rotation-matrix element of Realscalar"
1255  },
1256  UH32 == {
1257    ParameterType -> External,
1258    BlockName     -> OTHERS,
1259    OrderBlock    -> 17,
1260    Value         -> 0.000667955,
1261    TeX           -> Subsuperscript[U,32,H],
1262    Description   -> "32 rotation-matrix element of Realscalar"
1263  },
1264  UH33 == {
1265    ParameterType -> External,
1266    BlockName     -> OTHERS,
1267    OrderBlock    -> 18,
1268    Value         -> 0.999505,
1269    TeX           -> Subsuperscript[U,33,H],
1270    Description   -> "33 rotation-matrix element of Realscalar"
1271  },
1272  tz == {
1273    ParameterType    -> External,
1274    BlockName     -> OTHERS,
1275    OrderBlock    -> 19,
1276    Definitions      -> {tz -> 0},
1277    Description      -> "Tan of z zp mixing angle"
1278  },
1279  lam1 == {
1280    ParameterType    -> External,
1281    BlockName     -> OTHERS,
1282    OrderBlock    -> 20,
1283    Value            -> 1.03833,
1284    InteractionOrder -> {QED, 2},
1285    TeX              -> Subscript[\[Lambda], 1],
1286    Description      -> "Rho quartic coupling"
1287  },
1288  lamWS == {
1289    ParameterType    -> External,
1290    BlockName     -> WOLFENSTEIN,
1291    OrderBlock    -> 1,
1292    Value            -> 0.2253,
1293    TeX              -> \[Lambda],
1294    Description      -> "Wolfenstein variable lam"
1295  },
1296  AWS == {
1297    ParameterType    -> External,
1298    BlockName     -> WOLFENSTEIN,
1299    OrderBlock    -> 2,
1300    Value            -> 0.808,
1301    TeX              -> A,
1302    Description      -> "Wolfenstein variable A"
1303  },
1304  rhoWS == {
1305    ParameterType    -> External,
1306    BlockName     -> WOLFENSTEIN,
1307    OrderBlock    -> 3,
1308    Value            -> 0.132,
1309    TeX              -> \[Rho],
1310    Description      -> "Wolfenstein variable rho"
1311  },
1312  etaWS == {
1313    ParameterType    -> External,
1314    BlockName     -> WOLFENSTEIN,
1315    OrderBlock    -> 4,
1316    Value            -> 0.341,
1317    TeX              -> \[Eta],
1318    Description      -> "Wolfenstein variable eta"
1319  },
1320
1321
1322  (* Internal Parameters *)
1323  beta == {
1324    ParameterType    -> Internal,
1325    Definitions      -> {beta->1/Sqrt[3]},
1326    TeX              -> \[Beta],
1327    Description      -> "beta"
1328  },
1329  aEW == {
1330    ParameterType    -> Internal,
1331    Value            -> 1/aEWM1,
1332    InteractionOrder -> {QED,2},
1333    TeX              -> Subscript[\[Alpha], EW],
1334    Description      -> "Electroweak coupling contant"
1335  },
1336  ee == {
1337    ParameterType    -> Internal,
1338    Value            -> Sqrt[4 Pi aEW],
1339    InteractionOrder -> {QED,1},
1340    TeX              -> e, 
1341    Description      -> "Electric coupling constant"
1342  },
1343  cw == {
1344    ParameterType -> Internal,
1345    Value         -> MW/MZ,
1346    TeX           -> Subscript[c,w],
1347    Description   -> "Cosine of the Weinberg angle"
1348  },
1349  sw == {
1350    ParameterType -> Internal,
1351    Value         -> Sqrt[1-cw^2],
1352    TeX           -> Subscript[s,w],
1353    Description   -> "Sine of the Weinberg angle"
1354  },
1355  c3 == {
1356    ParameterType -> Internal,
1357    Definitions   -> {c3->beta sw/cw},
1358    TeX           -> Subscript[c,3],
1359    Description   -> "Cosine of the 331 angle"
1360  },
1361  s3 == {
1362    ParameterType -> Internal,
1363    Definitions   -> {s3->Sqrt[1-(1+beta^2)*sw^2]/cw},
1364    TeX           -> Subscript[s,3],
1365    Description   -> "Sine of the 331 angle"
1366  },
1367  gw == {
1368    ParameterType    -> Internal,
1369    Definitions      -> {gw->ee/sw},
1370    InteractionOrder -> {QED,1}, 
1371    TeX              -> Subscript[g,w],
1372    Description      -> "Weak coupling constant at the Z pole"
1373  },
1374  gx == {
1375    ParameterType    -> Internal,
1376    Definitions      -> {gx->gw sw/Sqrt[1-(1+beta^2)*sw^2]},
1377    InteractionOrder -> {QED,1}, 
1378    TeX              -> Subscript[g,x],
1379    Description      -> "U(1)X coupling constant at the Z pole"
1380  },
1381  gs == {
1382    ParameterType    -> Internal,
1383    Value            -> Sqrt[4 Pi aS],
1384    InteractionOrder -> {QCD,1}, 
1385    TeX              -> Subscript[g,s],
1386    ParameterName    -> G,
1387    Description      -> "Strong coupling constant at the Z pole"
1388  },
1389  v == {
1390    ParameterType -> Internal,
1391    Value         -> Sqrt[(Sqrt[2]Gf)^-1-v2^2],
1392    InteractionOrder -> {QED,-1},
1393    Description   -> "Rho vaccum expectation value"
1394  },
1395  v3 == {
1396    ParameterType -> Internal,
1397    Value         -> MZP Sqrt[3 - 4 sw^2]/(gw cw),
1398    InteractionOrder -> {QED,-1},
1399    Description   -> "Chi vaccum expectation value"
1400  },
1401  cz == {
1402    ParameterType    -> Internal,
1403    Definitions      -> {cz->1/Sqrt[1+tz^2]},
1404    Description      -> "Cosin of z zp mixing angle"
1405  },
1406  sz == {
1407    ParameterType    -> Internal,
1408    Definitions      -> {sz->-tz/Sqrt[1+tz^2]},
1409    Description      -> "Sin of z zp mixing angle"
1410  },
1411  svv2 == {
1412    ParameterType -> Internal,
1413    Value         -> 1/Sqrt[1+(v/v2)^2],
1414    TeX           -> Subscript[s,vv2],
1415    Description   -> "Sine of the vv2 angle"
1416  },
1417  cvv2 == {
1418    ParameterType -> Internal,
1419    Value         -> Sqrt[1-svv2^2],
1420    TeX           -> Subscript[c,vv2],
1421    Description   -> "Cosine of the vv2 angle"
1422  },
1423  svv3 == {
1424    ParameterType -> Internal,
1425    Value         -> 1/Sqrt[1+(v/v3)^2],
1426    TeX           -> Subscript[s,vv3],
1427    Description   -> "Sine of the vv3 angle"
1428  },
1429  cvv3 == {
1430    ParameterType -> Internal,
1431    Value         -> Sqrt[1-svv3^2],
1432    TeX           -> Subscript[c,vv3],
1433    Description   -> "Cosine of the vv3 angle"
1434  },
1435  sv2v3 == {
1436    ParameterType -> Internal,
1437    Value         -> 1/Sqrt[1+(v2/v3)^2],
1438    TeX           -> Subscript[s,v2v3],
1439    Description   -> "Sine of the v2v3 angle"
1440  },   
1441  cv2v3 == {
1442    ParameterType -> Internal,
1443    Value         -> Sqrt[1-sv2v3^2],
1444    TeX           -> Subscript[c,v2v3],
1445    Description   -> "Cosine of the v2v3 angle"
1446  },
1447  ff == {
1448    ParameterType    -> Internal,
1449    Value            -> -v3/(2Sqrt[2]),
1450    InteractionOrder -> {QED,1},
1451    TeX              -> f,
1452    Description      -> "Coefficient of the cubic piece of the Higgs potential"
1453  },
1454  mu1 == {
1455    ParameterType -> Internal,
1456    Definitions   -> {mu1->Sqrt[-lam1*v^2-lam12*v2^2/2-lam13*v3^2/2+ff*v3*v2/v]},
1457    Description   -> "Coefficient of the quadratic piece of the Rho potential"
1458  },
1459  mu2 == {
1460    ParameterType -> Internal,
1461    Definitions   -> {mu2->Sqrt[-lam2*v2^2-lam12*v^2/2-lam23*v3^2/2+ff*v3*v/v2]},
1462    Description   -> "Coefficient of the quadratic piece of the phi potential"
1463  },
1464  mu3 == {
1465    ParameterType -> Internal,
1466    Definitions   -> {mu3->Sqrt[-lam3*v3^2-lam13*v^2/2-lam23*v2^2/2+ff*v*v2/v3]},
1467    Description   -> "Coefficient of the quadratic piece of the Chi potential"
1468  },
1469  yl == {
1470    ParameterType    -> Internal,
1471    Indices          -> {Index[Generation], Index[Generation]},
1472    Definitions      -> {yl[i_?NumericQ, j_?NumericQ] :> 0  /; (i =!= j)},
1473    Value            -> {yl[1,1] -> Sqrt[2] yme / v2, yl[2,2] -> Sqrt[2] ymm / v2, yl[3,3] -> Sqrt[2] ymtau / v2},
1474    InteractionOrder -> {QED, 1},
1475    ParameterName    -> {yl[1,1] -> ye, yl[2,2] -> ym, yl[3,3] -> ytau},
1476    TeX              -> Superscript[y, l],
1477    Description      -> "Lepton Yukawa couplings"
1478  },
1479  yE == {
1480    ParameterType    -> Internal,
1481    Indices          -> {Index[Generation], Index[Generation]},
1482    Definitions      -> {yE[i_?NumericQ, j_?NumericQ] :> 0  /; (i =!= j)},
1483    Value            -> {yE[1,1] -> Sqrt[2] ymEe / v3, yE[2,2] -> Sqrt[2] ymEm / v3, yE[3,3] -> Sqrt[2] ymEtau / v3},
1484    InteractionOrder -> {QED, 1},
1485    ParameterName    -> {yE[1,1] -> yEe, yl[2,2] -> yEm, yl[3,3] -> yEtau},
1486    TeX              -> Superscript[y, E],
1487    Description      -> "Heavy Lepton Yukawa couplings"
1488  },
1489  yu == {
1490    ParameterType    -> Internal,
1491    Indices          -> {Index[Generation], Index[Generation]},
1492    Definitions      -> {yu[i_?NumericQ, j_?NumericQ] :> 0  /; (i =!= j)},
1493    Value            -> {yu[1,1] -> Sqrt[2] ymup/v2, yu[2,2] -> Sqrt[2] ymc/v2, yu[3,3] -> Sqrt[2] ymt/v2},
1494    InteractionOrder -> {QED, 1},
1495    ParameterName    -> {yu[1,1] -> yup, yu[2,2] -> yc, yu[3,3] -> yt},
1496    TeX              -> Superscript[y, u],
1497    Description      -> "Up-type Yukawa couplings"
1498  },
1499  yd == {
1500    ParameterType    -> Internal,
1501    Indices          -> {Index[Generation], Index[Generation]},
1502    Definitions      -> {yd[i_?NumericQ, j_?NumericQ] :> 0  /; (i =!= j)},
1503    Value            -> {yd[1,1] -> Sqrt[2] ymdo/v, yd[2,2] -> Sqrt[2] yms/v, yd[3,3] -> Sqrt[2] ymb/v},
1504    InteractionOrder -> {QED, 1},
1505    ParameterName    -> {yd[1,1] -> ydo, yd[2,2] -> ys, yd[3,3] -> yb},
1506    TeX              -> Superscript[y, d],
1507    Description      -> "Down-type Yukawa couplings"
1508  },
1509  yJ == {
1510    ParameterType    -> Internal,
1511    Indices          -> {Index[Generation], Index[Generation]},
1512    Definitions      -> {yJ[i_?NumericQ, j_?NumericQ] :> 0  /; (i =!= j)},
1513    Value            -> {yJ[1,1] -> Sqrt[2] ymD/v3, yJ[2,2] -> Sqrt[2] ymS/v3, yJ[3,3] -> Sqrt[2] ymT/v3},
1514    InteractionOrder -> {QED, 1},
1515    ParameterName    -> {yJ[1,1] -> yD, yJ[2,2] -> yS, yJ[3,3] -> yT},
1516    TeX              -> Superscript[y, J],
1517    Description      -> "Heavy-type Yukawa couplings"
1518  },
1519  MY == {
1520    ParameterType -> Internal,
1521    Value         -> 1/2 gw Sqrt[v3^2+v2^2],
1522    TeX           -> Subscript[M,Y],
1523    Description   -> "YY mass"
1524  },
1525  MV == {
1526    ParameterType -> Internal,
1527    Value         -> 1/2 gw Sqrt[v3^2+v^2],
1528    TeX           -> Subscript[M,V],
1529    Description   -> "V mass"
1530  },
1531  MH2 == {
1532    ParameterType -> Internal,
1533    Value         -> Sqrt[2] Sqrt[ff UH22 UH32 v-lam1 UH12^2 v^2+ff UH12 UH32 v2-lam12 UH12 UH22 v v2-lam2 UH22^2 v2^2-(ff UH32^2 v v2)/(2 v3)+ff UH12 UH22 v3-lam13 UH12 UH32 v v3-(ff UH22^2 v v3)/(2 v2)-lam23 UH22 UH32 v2 v3-(ff UH12^2 v2 v3)/(2 v)-lam3 UH32^2 v3^2],
1534    TeX           -> Subscript[M,H2],
1535    Description   -> "H2 mass"
1536  },
1537  MH3 == {
1538    ParameterType -> Internal,
1539    Value         -> Sqrt[2] Sqrt[ff UH23 UH33 v-lam1 UH13^2 v^2+ff UH13 UH33 v2-lam12 UH13 UH23 v v2-lam2 UH23^2 v2^2-(ff UH33^2 v v2)/(2 v3)+ff UH13 UH23 v3-lam13 UH13 UH33 v v3-(ff UH23^2 v v3)/(2 v2)-lam23 UH23 UH33 v2 v3-(ff UH13^2 v2 v3)/(2 v)-lam3 UH33^2 v3^2],
1540    TeX           -> Subscript[M,H3],
1541    Description   -> "H3 mass"
1542  },
1543  Uh11 == {
1544    ParameterType -> Internal,
1545    Value         -> v3/(v Sqrt[1 + (1/v^2 + 1/v2^2) v3^2]),
1546    TeX           -> Subsuperscript[U, 11, h],
1547    Description   -> "11 rotation-matrix element of Pseudoscalar"
1548  },
1549  Uh12 == {
1550    ParameterType -> Internal,
1551    Value         -> -v/(Sqrt[1 + v^2/v3^2] v3),
1552    TeX           -> Subsuperscript[U, 12, h],
1553    Description   -> "12 rotation-matrix element of Pseudoscalar"
1554  },
1555  Uh13 == {
1556    ParameterType -> Internal,
1557    Value         -> -(v v2 v3^2 (v^2 + v3^2))/(Sqrt[
1558 v2^2 (v^2 + v3^2)^2] Sqrt[(v^2 + v3^2) (v2^2 v3^2 +
1559    v^2 (v2^2 + v3^2))]),
1560    TeX           -> Subsuperscript[U, 13, h],
1561    Description   -> "13 rotation-matrix element of Pseudoscalar"
1562  },
1563  Uh21 == {
1564    ParameterType -> Internal,
1565    Value         -> v3/(v2 Sqrt[1 + (1/v^2 + 1/v2^2) v3^2]),
1566    TeX           -> Subsuperscript[U, 21, h],
1567    Description   -> "21 rotation-matrix element of Pseudoscalar"
1568  },
1569  Uh23 == {
1570    ParameterType -> Internal,
1571    Value         -> Sqrt[
1572 v2^2 (v^2 + v3^2)^2]/Sqrt[(v^2 + v3^2) (v2^2 v3^2 +
1573   v^2 (v2^2 + v3^2))],
1574    TeX           -> Subsuperscript[U, 23, h],
1575    Description   -> "23 rotation-matrix element of Pseudoscalar"
1576  },
1577  Uh31 == {
1578    ParameterType -> Internal,
1579    Value         -> 1/Sqrt[1 + (1/v^2 + 1/v2^2) v3^2],
1580    TeX           -> Subsuperscript[U, 31, h],
1581    Description   -> "31 rotation-matrix element of Pseudoscalar"
1582  },
1583  Uh32 == {
1584    ParameterType -> Internal,
1585    Value         -> 1/Sqrt[1 + v^2/v3^2],
1586    TeX           -> Subsuperscript[U, 32, h],
1587    Description   -> "32 rotation-matrix element of Pseudoscalar"
1588  },
1589  Uh33 == {
1590    ParameterType -> Internal,
1591    Value         -> -(v^2 v2 v3 (v^2 + v3^2))/(Sqrt[
1592 v2^2 (v^2 + v3^2)^2] Sqrt[(v^2 + v3^2) (v2^2 v3^2 +
1593    v^2 (v2^2 + v3^2))]),
1594    TeX           -> Subsuperscript[U, 33, h],
1595    Description   -> "33 rotation-matrix element of Pseudoscalar"
1596  },
1597  MH0 == {
1598    ParameterType -> Internal,
1599    Value         -> Sqrt[2] Sqrt[-ff Uh21 Uh31 v-ff Uh11 Uh31 v2-(ff Uh31^2 v v2)/(2 v3)-ff Uh11 Uh21 v3-(ff Uh21^2 v v3)/(2 v2)-(ff Uh11^2 v2 v3)/(2 v)],
1600    TeX           -> Subscript[M,H0],
1601    Description   -> "H0 mass"
1602  },
1603  MHW == {
1604    ParameterType -> Internal,
1605    Value         -> Sqrt[ (v^2+v2^2) ((-ff v3)/(v v2)-lam12P/2)],
1606    TeX           -> Subscript[M,HW],
1607    Description   -> "HW mass"
1608  },
1609  MHY == {
1610    ParameterType -> Internal,
1611    Value         -> Sqrt[ (v2^2+v3^2) ((-ff v)/(v2 v3)-lam23P/2)],
1612    TeX           -> Subscript[M,HY],
1613    Description   -> "HY mass"
1614  },
1615  MHV == {
1616    ParameterType -> Internal,
1617    Value         -> Sqrt[ (v^2+v3^2) ((-ff v2)/(v v3)-lam13P/2)],
1618    TeX           -> Subscript[M,HV],
1619    Description   -> "HV mass"
1620  },
1621
1622
1623
1624
1625(* N. B. : only Cabibbo mixing! *)
1626  CKM == {
1627    ParameterType -> Internal,
1628    Indices       -> {Index[Generation], Index[Generation]},
1629    Unitary       -> True,
1630    Value         -> {CKM[1,1] -> 1-lamWS^2/2,  CKM[1,2] -> lamWS, CKM[1,3] -> AWS*lamWS^3*(rhoWS-I*etaWS), CKM[2,1] -> -lamWS, CKM[2,2] -> 1-lamWS^2/2, CKM[2,3] -> AWS*lamWS^2, CKM[3,1] -> AWS*lamWS^3*(1-rhoWS-I*etaWS), CKM[3,2] -> -AWS*lamWS^2, CKM[3,3] -> 1},
1631    TeX         -> Superscript[V,CKM],
1632    Description -> "CKM-Matrix"
1633  },
1634  RU == {
1635    ParameterType -> Internal,
1636    Indices       -> {Index[Generation], Index[Generation]},
1637    Unitary       -> True,
1638    Definitions         -> {RU[1,1] -> 1,  RU[1,2] -> 0, RU[1,3] -> 0, RU[2,1] -> 0, RU[2,2] -> 1, RU[2,3] -> 0, RU[3,1] -> 0, RU[3,2] -> 0, RU[3,3] -> 1},
1639    TeX         -> Superscript[R,u],
1640    Description -> "RU-Matrix"
1641  }
1642
1643};
1644
1645(* ************************** *)
1646(* *****   Lagrangian   ***** *)
1647(* ************************** *)
1648
1649LGauge := Block[{mu,nu,ii,aa},
1650
1651  ExpandIndices[-1/4 FS[K,mu,nu] FS[K,mu,nu] - 1/4 FS[Wi,mu,nu,ii]FS[Wi,mu,nu,ii] - 1/4 FS[G,mu,nu,aa] FS[G,mu,nu,aa], FlavorExpand->SU3W]];
1652
1653
1654LFermions := Block[{mu,fermi},
1655 
1656  fermi=ExpandIndices[I*(
1657QL12bar.Ga[mu].DC[QL12, mu] + QL3bar.Ga[mu].DC[QL3, mu] + LLbar.Ga[mu].DC[LL, mu] + uRbar.Ga[mu].DC[uR, mu] + dRbar.Ga[mu].DC[dR, mu] + JR12bar.Ga[mu].DC[JR12,mu] + JR3bar.Ga[mu].DC[JR3,mu] + lRbar.Ga[mu].DC[lR, mu] + EERbar.Ga[mu].DC[EER, mu]),
1658  FlavorExpand->{SU3W,SU3T,ASU3W,ASU3T}];
1659  fermi = ExpandIndices[fermi]/.{Conjugate[RU[a_,b_]] RU[a_,a_] -> 1 - Conjugate[RU[2,3]] RU[2,1]- Conjugate[RU[3,3]] RU[3,1]}];
1660
1661LHiggs := Block[{ii,mu, feynmangaugerules},
1662  feynmangaugerules = If[Not[FeynmanGauge], {GZ|GZP|GW|GWbar|GY|GYbar|GV|GVP ->0}, {}];
1663 
1664  ExpandIndices[DC[Rhobar[ii],mu] DC[Rho[ii],mu] + DC[Phibar[ii],mu] DC[Phi[ii],mu] + DC[Chibar[ii],mu] DC[Chi[ii],mu] + mu1^2 Rhobar[ii] Rho[ii] + lam1 Rhobar[ii] Rho[ii] Rhobar[jj] Rho[jj] + mu2^2 Phibar[ii] Phi[ii] + lam2 Phibar[ii] Phi[ii] Phibar[jj] Phi[jj] + mu3^2 Chibar[ii] Chi[ii] + lam3 Chibar[ii] Chi[ii] Chibar[jj] Chi[jj] + lam12 Rhobar[ii] Rho[ii] Phibar[jj] Phi[jj] + lam13 Rhobar[ii] Rho[ii] Chibar[jj] Chi[jj] + lam23 Phibar[ii] Phi[ii] Chibar[jj] Chi[jj] + lam12P Rhobar[ii] Phi[ii] Phibar[jj] Rho[jj] + lam13P Rhobar[ii] Chi[ii] Chibar[jj] Rho[jj] + lam23P Phibar[ii] Chi[ii] Chibar[jj] Phi[jj] + Sqrt[2] ff Eps[ii,jj,kk] Rho[ii] Phi[jj] Chi[kk] + Sqrt[2] ff HC[Eps[ii,jj,kk] Rho[ii] Phi[jj] Chi[kk]], FlavorExpand->{SU3T,SU3W}]/.feynmangaugerules
1665 ];
1666
1667LYukawa := Block[{sp,ii,cc,ff1,ff3,ff4,yuk,feynmangaugerules},
1668  feynmangaugerules = If[Not[FeynmanGauge], {GZ|GZP|GW|GWbar|GY|GYbar|GV|GVP ->0}, {}];
1669 
1670  yuk = ExpandIndices[
1671   -yl[ff1, ff3] LLbar[sp, ii, ff1].lR [sp, ff3] Phibar[ii]
1672   -yE[ff1, ff3] LLbar[sp, ii, ff1].EER [sp, ff3] Chibar[ii]
1673   -yu[ff1, ff3] RU[ff2, ff1]  QLbar[sp, ii, ff4, cc].uR [sp, ff3, cc] Su[ff4, ff1, ii]
1674   -yd[ff1, ff3] CKM[ff2, ff1] QLbar[sp, ii, ff4, cc].dR [sp, ff3, cc] Sd[ff4, ff2, ii]
1675   -yJ[ff1, ff3] QLbar[sp, ii, ff2, cc].JR[sp, ff3, cc] SJ[ff2, ff1, ii],
1676FlavorExpand -> {SU3T}];
1677  yuk = ExpandIndices[yuk] /.{CKM[a_, b_] Conjugate[CKM[a_, c_]] -> 1/3 IndexDelta[b, c],CKM[b_, a_] Conjugate[CKM[c_, a_]] -> 1/3 IndexDelta[b, c]};
1678  yuk+HC[yuk]/.feynmangaugerules
1679 ];
1680
1681LGhost := Block[{LGh1,LGhw,LGhs,LGhhiggs,LGhrho,LGhphi,LGhchi,mu, generators,gh,ghbar,Vectorize,rho1,rho2,rho3,rho4,phi1,phi2,phi3,phi4,chi1,chi2,chi3,chi4,togoldstones,rho,rho0,phi,phi0,chi,chi0},
1682  (* Pure gauge piece *)       
1683  LGh1 = -ghKbar.del[DC[ghK,mu],mu];
1684  LGhw = -ghWibar.del[DC[ghWi,mu],mu];
1685  LGhs = -ghGbar.del[DC[ghG,mu],mu];
1686
1687  (* Scalar pieces: see Peskin pages 739-742 *)
1688  (* rho1, rho2, rho3 and rho4 are the real degrees of freedom of HW and GW *)
1689  (* Vectorize transforms a triplet in a vector in the phi-basis, i.e. the basis of real degrees of freedom *)
1690  gh    = {ghK, ghWi[1], ghWi[2], ghWi[3], ghWi[4], ghWi[5], ghWi[6], ghWi[7], ghWi[8]};
1691  ghbar = {ghKbar, ghWibar[1], ghWibar[2], ghWibar[3], ghWibar[4], ghWibar[5], ghWibar[6], ghWibar[7], ghWibar[8]};
1692  generators = {-I/2 gx IdentityMatrix[3], -I/2 gw Gellmann[1], -I/2 gw Gellmann[2], -I/2 gw Gellmann[3], -I/2 gw Gellmann[4], -I/2 gw Gellmann[5], -I/2 gw Gellmann[6], -I/2 gw Gellmann[7], -I/2 gw Gellmann[8]};
1693  rho = Expand[{(-I rho1 - rho2)/Sqrt[2], (UH11 h + UH12 H2 + UH13 H3 + I (Uh11 H0 + Uh12 GZ + Uh13 GZP))/Sqrt[2], (-I rho3 - rho4)/Sqrt[2] } ];
1694  rho0 = {0, v/Sqrt[2], 0};
1695  phi = Expand[{(UH21 h + UH22 H2 + UH23 H3 + I (Uh21 H0  + Uh23 GZP))/Sqrt[2], (-I phi1 - phi2)/Sqrt[2], (-I phi3 - phi4)/Sqrt[2] } ];
1696  phi0 = {v2/Sqrt[2], 0, 0};
1697  chi = Expand[{(-I chi1 - chi2)/Sqrt[2], (-I chi3 - chi4)/Sqrt[2], (UH31 h + UH32 H2 + UH33 H3 + I (Uh31 H0 + Uh32 GZ + Uh33 GZP))/Sqrt[2] } ];
1698  chi0 = {0, 0, v3/Sqrt[2]};
1699  Vectorize[{a_, b_, c_}]:= 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]]}/.{Im[_]->0, Re[num_]->num}];
1700  togoldstones := {
1701rho1 -> (HW svv2 + GW cvv2 + HWbar svv2 + GWbar cvv2)/Sqrt[2],
1702rho2 -> (- HW svv2 - GW cvv2  + HWbar svv2 + GWbar cvv2)/(I Sqrt[2]),
1703rho3 -> (GV cvv3 + HV svv3 + GVbar cvv3 + HVbar svv3)/Sqrt[2],
1704rho4 -> (- GV cvv3 - HV svv3  + GVbar cvv3 + HVbar svv3)/(I Sqrt[2]),
1705phi1 -> (- HW cvv2 + GW svv2 - HWbar cvv2 + GWbar svv2)/Sqrt[2],
1706phi2 -> (HWbar cvv2 - GWbar svv2 - HW cvv2 + GW svv2)/(I Sqrt[2]),
1707phi3 -> (GY cv2v3 + HY sv2v3 + GYbar cv2v3 + HYbar sv2v3)/Sqrt[2],
1708phi4 -> (- GY cv2v3 - HY sv2v3 + GYbar cv2v3 + HYbar sv2v3)/(I Sqrt[2]),
1709chi1 -> (- HY cv2v3 + GY sv2v3 - HYbar cv2v3 + GYbar sv2v3)/Sqrt[2],
1710chi2 -> (HYbar cv2v3 - GYbar sv2v3 - HY cv2v3 + GY sv2v3)/(I Sqrt[2]),
1711chi3 -> (- HV cvv3 + GV svv3 - HVbar cvv3 + GVbar svv3)/Sqrt[2],
1712chi4 -> (HVbar cvv3 - GVbar svv3 - HV cvv3 + GV svv3)/(I Sqrt[2])};
1713  LGhrho=
1714Plus@@Flatten[Table[-ghbar[[kkk]].gh[[lll]] Vectorize[generators[[kkk]].(rho0)].Vectorize[generators[[lll]].(rho + rho0)],{kkk,9},{lll,9}]];
1715  LGhphi=
1716Plus@@Flatten[Table[-ghbar[[kkkk]].gh[[llll]] Vectorize[generators[[kkkk]].(phi0)].Vectorize[generators[[llll]].(phi + phi0)],{kkkk,9},{llll,9}]];
1717  LGhchi=
1718Plus@@Flatten[Table[-ghbar[[kkkkk]].gh[[lllll]] Vectorize[generators[[kkkkk]].(chi0)].Vectorize[generators[[lllll]].(chi + chi0)],{kkkkk,9},{lllll,9}]];
1719
1720  LGhhiggs=LGhrho+ LGhphi+ LGhchi /.togoldstones;
1721
1722ExpandIndices[ LGhs + If[FeynmanGauge, LGhw +LGh1 + LGhhiggs ,0], FlavorExpand->SU3W]];
1723
1724LThree:= LHiggs+LFermions+LYukawa+LGhost+LGauge;