BSMCharacterisation: BSMC_v2.2.fr

File BSMC_v2.2.fr, 67.8 KB (added by BenjF, 3 months ago)

BSMC FR model, version 2.2

Line 
1(*****************************************************************************)
2(*  This is the FeynRules model file for the BSM Characterisation            *)
3(*****************************************************************************)
4
5(* ************************** *)
6(* *****  Information   ***** *)
7(* ************************** *)
8
9M$ModelName = "BSMC";
10
11M$Information = {
12 Authors      -> {"B. Fuks", "K. Mawatari"},
13 Version      -> "2.2",
14 Date         -> "16.07.2018",
15 Institutions -> {"LPTHE Paris / Sorbonne U.", "Osaka U."},
16 Emails       -> {"fuks@lpthe.jussieu.fr", "kentarou.mawatari@het.phys.sci.osaka-u.ac.jp"},
17 URLs         -> "http://feynrules.irmp.ucl.ac.be/wiki/BSMCharacterisation",
18 References   -> {"A. Falkowski, B. Fuks, K. Mawatari, K. Mimasu, F. Riva and V. Sanz, arXiv:1508.05895"}
19};
20
21
22(* ************************** *)
23(* *****  Change  log   ***** *)
24(* ************************** *)
25
26(* 14.07.2015 v1.0 - release version (based on the Higgs basis note (Jun/14)).                    *)
27(* 27.07.2015 v1.1 - modifying the flavor structure.                                              *)
28(* 07.10.2015 v1.2 - the dipole operators now match the HEFT note.                                *)
29(* 04.03.2016 v2.0 - version corresponding to the Rosetta paper arXiv:1508.05895.                 *)
30(* 31.05.2017 v2.1 - change of the signs of the SM Hgg, HaZ and Haa couplings (thanks S. Di Vita) *)
31(* 16.07.2018 v2.2 - Bug in the dipole operators  + normalization of one operator                 *)
32(*                   (thanks G. Durieux, Marc Riembau and Thibaud Vantalon)                       *)
33
34
35FeynmanGauge = False;
36
37(* ************************** *)
38(* **  Interaction  orders ** *)
39(* ************************** *)
40
41M$InteractionOrderLimit = { {QNP, 2} };
42M$InteractionOrderHierarchy = { {QNP, 1} };
43
44
45(* ************************** *)
46(* *****   Parameters   ***** *)
47(* ************************** *)
48M$Parameters = {
49
50  (* Quadratic terms *)
51  dM  == { TeX -> "\[Delta]m", ParameterType -> External, Value -> 0.1, BlockName -> BCxmass,  OrderBlock -> 1, Description -> "W-mass corrections"},
52
53  (* Lvertex parameters *)
54  RdGLzu ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdGLzu[i_,j_]->0.1}, BlockName ->   BCxdGLzu, Description -> "Z-u vertex interactions (real part, L)"},
55  IdGLzu ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdGLzu[i_,j_]->0.0}, BlockName -> IMBCxdGLzu, Description -> "Z-u vertex interactions (im part, L)"  },
56  RdGLzd ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdGLzd[i_,j_]->0.1}, BlockName ->   BCxdGLzd, Description -> "Z-d vertex interactions (real part, L)"},
57  IdGLzd ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdGLzd[i_,j_]->0.0}, BlockName -> IMBCxdGLzd, Description -> "Z-d vertex interactions (im part, L)"  },
58  RdGLze ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdGLze[i_,j_]->0.1}, BlockName ->   BCxdGLze, Description -> "Z-l vertex interactions (real part, L)"},
59  IdGLze ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdGLze[i_,j_]->0.0}, BlockName -> IMBCxdGLze, Description -> "Z-l vertex interactions (im part, L)"},
60  RdGLzv ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdGLzv[i_,j_]->0.1}, BlockName ->   BCxdGLzv, Description -> "Z-v vertex interactions (real part, L)"},
61  IdGLzv ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdGLzv[i_,j_]->0.0}, BlockName -> IMBCxdGLzv, Description -> "Z-v vertex interactions (im part, L)"},
62  RdGRzu ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdGRzu[i_,j_]->0.1}, BlockName ->   BCxdGRzu, Description -> "Z-u vertex interactions (real part, R)"},
63  IdGRzu ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdGRzu[i_,j_]->0.0}, BlockName -> IMBCxdGRzu, Description -> "Z-u vertex interactions (im part, R)"  },
64  RdGRzd ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdGRzd[i_,j_]->0.1}, BlockName ->   BCxdGRzd, Description -> "Z-d vertex interactions (real part, R)"},
65  IdGRzd ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdGRzd[i_,j_]->0.0}, BlockName -> IMBCxdGRzd, Description -> "Z-d vertex interactions (im part, R)"  },
66  RdGRze ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdGRze[i_,j_]->0.1}, BlockName ->   BCxdGRze, Description -> "Z-l vertex interactions (real part, R)"},
67  IdGRze ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdGRze[i_,j_]->0.0}, BlockName -> IMBCxdGRze, Description -> "Z-l vertex interactions (im part, R)"},
68  RdGLwq ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdGLwq[i_,j_]->0.1}, BlockName ->   BCxdGLWq, Description -> "W-q vertex interactions (real part, L)"},
69  IdGLwq ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdGLwq[i_,j_]->0.0}, BlockName -> IMBCxdGLWq, Description -> "W-q vertex interactions (im part, L)"  },
70  RdGRwq ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdGRwq[i_,j_]->0.1}, BlockName ->   BCxdGRWq, Description -> "W-q vertex interactions (real part, R)"},
71  IdGRwq ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdGRwq[i_,j_]->0.0}, BlockName -> IMBCxdGRWq, Description -> "W-q vertex interactions (im part, R)"  },
72  RdGLwl ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdGLwl[i_,j_]->0.1}, BlockName ->   BCxdGLWl, Description -> "W-l vertex interactions (real part, L)"},
73  IdGLwl ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdGLwl[i_,j_]->0.0}, BlockName -> IMBCxdGLWl, Description -> "W-l vertex interactions (im part, L)"  },
74
75  dGLzu ==  { TeX -> Subsuperscript["\[Delta]g", L, Zu], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
76    Value -> {dGLzu[i_,j_]->RdGLzu[i,j] + I IdGLzu[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "Z-u vertex interactions (L)"},
77  dGLzd ==  { TeX -> Subsuperscript["\[Delta]g", L, Zd], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
78    Value -> {dGLzd[i_,j_]->RdGLzd[i,j] + I IdGLzd[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "Z-d vertex interactions (L)"},
79  dGLze ==  { TeX -> Subsuperscript["\[Delta]g", L, Zl], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
80    Value -> {dGLze[i_,j_]->RdGLze[i,j] + I IdGLze[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "Z-l vertex interactions (L)"},
81  dGLzv ==  { TeX -> Subsuperscript["\[Delta]g", L, Zv], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
82    Value -> {dGLzv[i_,j_]->RdGLzv[i,j] + I IdGLzv[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "Z-v vertex interactions (L)"},
83  dGRzu ==  { TeX -> Subsuperscript["\[Delta]g", R, Zu], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
84    Value -> {dGRzu[i_,j_]->RdGRzu[i,j] + I IdGRzu[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "Z-u vertex interactions (R)"},
85  dGRzd ==  { TeX -> Subsuperscript["\[Delta]g", R, Zd], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
86    Value -> {dGRzd[i_,j_]->RdGRzd[i,j] + I IdGRzd[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "Z-d vertex interactions (R)"},
87  dGRze ==  { TeX -> Subsuperscript["\[Delta]g", R, Zl], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
88    Value -> {dGRze[i_,j_]->RdGRze[i,j] + I IdGRze[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "Z-l vertex interactions (R)"},
89  dGLwq ==  { TeX -> Subsuperscript["\[Delta]g", L, Wq], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
90    Value -> {dGLwq[i_,j_]->RdGLwq[i,j] + I IdGLwq[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "W-q vertex interactions (L)"},
91  dGRwq ==  { TeX -> Subsuperscript["\[Delta]g", R, Wq], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
92    Value -> {dGRwq[i_,j_]->RdGRwq[i,j] + I IdGRwq[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "W-q vertex interactions (R)"},
93  dGLwl ==  { TeX -> Subsuperscript["\[Delta]g", L, Wl], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
94    Value -> {dGLwl[i_,j_]->RdGLwl[i,j] + I IdGLwl[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "W-l vertex interactions (L)"},
95
96  (* Ldip parameters *)
97  Rdgu ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdgu[i_,j_]->0.1}, BlockName ->   BCxdgu, Description -> "G-u dipole interactions (real part)"},
98  Idgu ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idgu[i_,j_]->0.0}, BlockName -> IMBCxdgu, Description -> "G-u dipole interactions (im part)"  },
99  Rdgd ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdgd[i_,j_]->0.1}, BlockName ->   BCxdgd, Description -> "G-d dipole interactions (real part)"},
100  Idgd ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idgd[i_,j_]->0.0}, BlockName -> IMBCxdgd, Description -> "G-d dipole interactions (im part)"  },
101  Rdau ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdau[i_,j_]->0.1}, BlockName ->   BCxdau, Description -> "A-u dipole interactions (real part)"},
102  Idau ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idau[i_,j_]->0.0}, BlockName -> IMBCxdau, Description -> "A-u dipole interactions (im part)"  },
103  Rdad ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdad[i_,j_]->0.1}, BlockName ->   BCxdad, Description -> "A-d dipole interactions (real part)"},
104  Idad ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idad[i_,j_]->0.0}, BlockName -> IMBCxdad, Description -> "A-d dipole interactions (im part)"  },
105  Rdae ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdae[i_,j_]->0.1}, BlockName ->   BCxdae, Description -> "A-e dipole interactions (real part)"},
106  Idae ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idae[i_,j_]->0.0}, BlockName -> IMBCxdae, Description -> "A-e dipole interactions (im part)"  },
107  Rdzu ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdzu[i_,j_]->0.1}, BlockName ->   BCxdzu, Description -> "Z-u dipole interactions (real part)"},
108  Idzu ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idzu[i_,j_]->0.0}, BlockName -> IMBCxdzu, Description -> "Z-u dipole interactions (im part)"  },
109  Rdzd ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdzd[i_,j_]->0.1}, BlockName ->   BCxdzd, Description -> "Z-d dipole interactions (real part)"},
110  Idzd ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idzd[i_,j_]->0.0}, BlockName -> IMBCxdzd, Description -> "Z-d dipole interactions (im part)"  },
111  Rdze ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdze[i_,j_]->0.1}, BlockName ->   BCxdze, Description -> "Z-e dipole interactions (real part)"},
112  Idze ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idze[i_,j_]->0.0}, BlockName -> IMBCxdze, Description -> "Z-e dipole interactions (im part)"  },
113  Rdwd ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdwd[i_,j_]->0.1}, BlockName ->   BCxdwd, Description -> "W-d dipole interactions (real part)"},
114  Idwd ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idwd[i_,j_]->0.0}, BlockName -> IMBCxdwd, Description -> "W-d dipole interactions (im part)"  },
115  Rdwu ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdwu[i_,j_]->0.1}, BlockName ->   BCxdwu, Description -> "W-u dipole interactions (real part)"},
116  Idwu ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idwu[i_,j_]->0.0}, BlockName -> IMBCxdwu, Description -> "W-u dipole interactions (im part)"  },
117  Rdwl ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdwl[i_,j_]->0.1}, BlockName ->   BCxdwl, Description -> "W-l dipole interactions (real part)"},
118  Idwl ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idwl[i_,j_]->0.0}, BlockName -> IMBCxdwl, Description -> "W-l dipole interactions (im part)"  },
119  Rtdgu == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rtdgu[i_,j_]->0.1}, BlockName ->   BCxtdgu, Description -> "G-u dipole interactions (real part)"},
120  Itdgu == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Itdgu[i_,j_]->0.0}, BlockName -> IMBCxtdgu, Description -> "G-u dipole interactions (im part, CP-odd)  "},
121  Rtdgd == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rtdgd[i_,j_]->0.1}, BlockName ->   BCxtdgd, Description -> "G-d dipole interactions (real part, CP-odd)"},
122  Itdgd == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Itdgd[i_,j_]->0.0}, BlockName -> IMBCxtdgd, Description -> "G-d dipole interactions (im part, CP-odd)"  },
123  Rtdau == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rtdau[i_,j_]->0.1}, BlockName ->   BCxtdau, Description -> "A-u dipole interactions (real part, CP-odd)"},
124  Itdau == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Itdau[i_,j_]->0.0}, BlockName -> IMBCxtdau, Description -> "A-u dipole interactions (im part, CP-odd)"  },
125  Rtdad == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rtdad[i_,j_]->0.1}, BlockName ->   BCxtdad, Description -> "A-d dipole interactions (real part, CP-odd)"},
126  Itdad == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Itdad[i_,j_]->0.0}, BlockName -> IMBCxtdad, Description -> "A-d dipole interactions (im part, CP-odd)"  },
127  Rtdae == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rtdae[i_,j_]->0.1}, BlockName ->   BCxtdae, Description -> "A-e dipole interactions (real part, CP-odd)"},
128  Itdae == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Itdae[i_,j_]->0.0}, BlockName -> IMBCxtdae, Description -> "A-e dipole interactions (im part, CP-odd)"  },
129  Rtdzu == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rtdzu[i_,j_]->0.1}, BlockName ->   BCxtdzu, Description -> "Z-u dipole interactions (real part, CP-odd)"},
130  Itdzu == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Itdzu[i_,j_]->0.0}, BlockName -> IMBCxtdzu, Description -> "Z-u dipole interactions (im part, CP-odd)"  },
131  Rtdzd == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rtdzd[i_,j_]->0.1}, BlockName ->   BCxtdzd, Description -> "Z-d dipole interactions (real part, CP-odd)"},
132  Itdzd == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Itdzd[i_,j_]->0.0}, BlockName -> IMBCxtdzd, Description -> "Z-d dipole interactions (im part, CP-odd)"  },
133  Rtdze == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rtdze[i_,j_]->0.1}, BlockName ->   BCxtdze, Description -> "Z-e dipole interactions (real part, CP-odd)"},
134  Itdze == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Itdze[i_,j_]->0.0}, BlockName -> IMBCxtdze, Description -> "Z-e dipole interactions (im part, CP-odd)"  },
135
136  dgu ==  { TeX -> Subscript["d", gu], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
137    Value -> {dgu[i_,j_]->Rdgu[i,j] + I Idgu[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "G-u dipole interactions"},
138  dgd ==  { TeX -> Subscript["d", gd], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
139    Value -> {dgd[i_,j_]->Rdgd[i,j] + I Idgd[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "G-d dipole interactions"},
140  dau ==  { TeX -> Subscript["d", au], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
141    Value -> {dau[i_,j_]->Rdau[i,j] + I Idau[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "A-u dipole interactions"},
142  dad ==  { TeX -> Subscript["d", ad], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
143    Value -> {dad[i_,j_]->Rdad[i,j] + I Idad[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "A-d dipole interactions"},
144  dae ==  { TeX -> Subscript["d", ae], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
145    Value -> {dae[i_,j_]->Rdae[i,j] + I Idae[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "A-e dipole interactions"},
146  dzu ==  { TeX -> Subscript["d", zu], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
147    Value -> {dzu[i_,j_]->Rdzu[i,j] + I Idzu[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "Z-u dipole interactions"},
148  dzd ==  { TeX -> Subscript["d", zd], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
149    Value -> {dzd[i_,j_]->Rdzd[i,j] + I Idzd[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "Z-d dipole interactions"},
150  dze ==  { TeX -> Subscript["d", ze], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
151    Value -> {dze[i_,j_]->Rdze[i,j] + I Idze[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "Z-e dipole interactions"},
152  dwu ==  { TeX -> Subscript["d", wu], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
153    Value -> {dwu[i_,j_]->Rdwu[i,j] + I Idwu[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "w-u dipole interactions"},
154  dwd ==  { TeX -> Subscript["d", wd], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
155    Value -> {dwd[i_,j_]->Rdwd[i,j] + I Idwd[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "w-d dipole interactions"},
156  dwl ==  { TeX -> Subscript["d", wl], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
157    Value -> {dwl[i_,j_]->Rdwl[i,j] + I Idwl[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "W-l dipole interactions"},
158  tdgu ==  { TeX -> Subscript[OverTilde["d"], gu], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
159    Value -> {tdgu[i_,j_]->Rtdgu[i,j] + I Itdgu[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "G-u dipole interactions (CP-odd)"},
160  tdgd ==  { TeX -> Subscript[OverTilde["d"], gd], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
161    Value -> {tdgd[i_,j_]->Rtdgd[i,j] + I Itdgd[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "G-d dipole interactions (CP-odd)"},
162  tdau ==  { TeX -> Subscript[OverTilde["d"], au], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
163    Value -> {tdau[i_,j_]->Rtdau[i,j] + I Itdau[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "A-u dipole interactions (CP-odd)"},
164  tdad ==  { TeX -> Subscript[OverTilde["d"], ad], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
165    Value -> {tdad[i_,j_]->Rtdad[i,j] + I Itdad[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "A-d dipole interactions (CP-odd)"},
166  tdae ==  { TeX -> Subscript[OverTilde["d"], ae], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
167    Value -> {tdae[i_,j_]->Rtdae[i,j] + I Itdae[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "A-e dipole interactions (CP-odd)"},
168  tdzu ==  { TeX -> Subscript[OverTilde["d"], zu], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
169    Value -> {tdzu[i_,j_]->Rtdzu[i,j] + I Itdzu[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "Z-u dipole interactions (CP-odd)"},
170  tdzd ==  { TeX -> Subscript[OverTilde["d"], zd], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
171    Value -> {tdzd[i_,j_]->Rtdzd[i,j] + I Itdzd[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "Z-d dipole interactions (CP-odd)"},
172  tdze ==  { TeX -> Subscript[OverTilde["d"], ze], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
173    Value -> {tdze[i_,j_]->Rtdze[i,j] + I Itdze[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "Z-e dipole interactions (CP-odd)"},
174
175 (* Ltgc parameters *)
176  dKa  == { TeX -> Subscript["\[Delta]\[Kappa]",a], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxtgc, OrderBlock -> 1, Description -> "fsA-W-W interactions"},
177  dKz  == { TeX -> Subscript["\[Delta]\[Kappa]",z], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxtgc, OrderBlock -> 2, Description -> "fsZ-W-W interactions"},
178  dG1z == { TeX -> Subscript["\[Delta]g","1z"],     ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxtgc, OrderBlock -> 3, Description -> "fsW-W-Z interactions"},
179  La   == { TeX -> Subscript["\[Lambda]",a],        ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxtgc, OrderBlock -> 4, Description -> "fsW-fsW-fsA interactions"},
180  Lz   == { TeX -> Subscript["\[Lambda]",z],        ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxtgc, OrderBlock -> 5, Description -> "fsW-fsW-fsZ interactions"},
181  C3g  == { TeX -> Subscript[c,"3G"],               ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxtgc, OrderBlock -> 6, Description -> "fsG-fsG-fsG interactions"},
182  tKa == { TeX -> Subscript[OverTilde["\[Kappa]"],a],  ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxtgc, OrderBlock -> 7,  Description -> "fsA-W-W interactions (CP-odd)"},
183  tKz == { TeX -> Subscript[OverTilde["\[Kappa]"],z],  ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxtgc, OrderBlock -> 8,  Description -> "fsZ-W-W interactions (CP-odd)"},
184  tLa == { TeX -> Subscript[OverTilde["\[Lambda]"],a], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxtgc, OrderBlock -> 9,  Description -> "fsW-fsW-fsA interactions (CP-odd)"},
185  tLz == { TeX -> Subscript[OverTilde["\[Lambda]"],z], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxtgc, OrderBlock -> 10, Description -> "fsW-fsW-fsZ interactions (CP-odd)"},
186  tC3g== { TeX -> Subscript[OverTilde["c"],"3g"],      ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxtgc, OrderBlock -> 11, Description -> "fsG-fsG-fsG interactions (CP-odd)"},
187
188 (* Lqgc parameters *)
189  dGw4   == { TeX -> Subscript["\[Delta]g",w4],   ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxqgc, OrderBlock -> 1, Description -> "W^4 interactions"},
190  dGw2z2 == { TeX -> Subscript["\[Delta]g",w2z2], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxqgc, OrderBlock -> 2, Description -> "W^2-Z^2 interactions"},
191  dGw2za == { TeX -> Subscript["\[Delta]g",w2za], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxqgc, OrderBlock -> 3, Description -> "W^2-Z-A interactions"},
192  Lw4    == { TeX -> Subscript["\[Lambda]",w4],   ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxqgc, OrderBlock -> 4, Description -> "fsW^2-W^2 interactions"},
193  Lw2z2  == { TeX -> Subscript["\[Lambda]",w2z2], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxqgc, OrderBlock -> 5, Description -> "fsW-fsZ-W-Z interactions"},
194  Lw2a2  == { TeX -> Subscript["\[Lambda]",w2a2], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxqgc, OrderBlock -> 6, Description -> "fsW-fsA-W-A interactions"},
195  Lw2az  == { TeX -> Subscript["\[Lambda]",w2az], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxqgc, OrderBlock -> 7, Description -> "fsW-fsA-W-Z interactions"},
196  Lw2za  == { TeX -> Subscript["\[Lambda]",w2za], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxqgc, OrderBlock -> 8, Description -> "fsW-fsZ-W-A interactions"},
197  C4g    == { TeX -> Subscript[c,"4G"],           ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxqgc, OrderBlock -> 9, Description -> "fsG-fsG-G-G interactions"},
198  tLw4   == { TeX -> Subscript[OverTilde["\[Lambda]"],w4],   ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxqgc, OrderBlock -> 10, Description -> "fsW^2-W^2 interactions (CP-odd)"},
199  tLw2z2 == { TeX -> Subscript[OverTilde["\[Lambda]"],w2z2], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxqgc, OrderBlock -> 11, Description -> "fsW-fsZ-W-Z interactions (CP-odd)"},
200  tLw2a2 == { TeX -> Subscript[OverTilde["\[Lambda]"],w2a2], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxqgc, OrderBlock -> 12, Description -> "fsW-fsA-W-A interactions (CP-odd)"},
201  tLw2az == { TeX -> Subscript[OverTilde["\[Lambda]"],w2az], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxqgc, OrderBlock -> 13, Description -> "fsW-fsA-W-Z interactions (CP-odd)"},
202  tLw2za == { TeX -> Subscript[OverTilde["\[Lambda]"],w2za], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxqgc, OrderBlock -> 14, Description -> "fsW-fsZ-W-A interactions (CP-odd)"},
203  tC4g   == { TeX -> Subscript[OverTilde[c],"4G"],           ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxqgc, OrderBlock -> 15, Description -> "fsG-fsG-G-G interactions (CP-odd)"},
204
205  (* Lh parameters *)
206  dCw  == { TeX -> Subscript["\[Delta]c",w],   ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxh, OrderBlock -> 1, Description -> "W-W-h interactions"},
207  dCz  == { TeX -> Subscript["\[Delta]c",z],   ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxh, OrderBlock -> 2, Description -> "Z-Z-h interactions"},
208  Cww  == { TeX -> Subscript[c,ww],            ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxh, OrderBlock -> 3, Description -> "fsW-fsW-h interactions"},
209  Cgg  == { TeX -> Subscript[c,gg],            ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxh, OrderBlock -> 4, Description -> "fsG-fsG-h interactions"},
210  Caa  == { TeX -> Subscript[c,aa],            ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxh, OrderBlock -> 5, Description -> "fsA-fsA-h interactions"},
211  Cza  == { TeX -> Subscript[c,za],            ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxh, OrderBlock -> 6, Description -> "fsZ-fsA-h interactions"},
212  Czz  == { TeX -> Subscript[c,zz],            ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxh, OrderBlock -> 7, Description -> "fsZ-fsZ-h interactions"},
213  Cwbx == { TeX -> Subscript[c,w,Square[""]],  ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxh, OrderBlock -> 8, Description -> "fsW-W-h interactions"},
214  Czbx == { TeX -> Subscript[c,z,Square[""]],  ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxh, OrderBlock -> 9, Description -> "fsZ-Z-h interactions"},
215  Cabx == { TeX -> Subscript[c,a,Square[""]],  ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxh, OrderBlock -> 10,Description -> "fsA-Z-h interactions"},
216  tCww  == { TeX -> Subscript[OverTilde[c],ww], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxh, OrderBlock -> 11, Description -> "fsW-fsW-h interactions (CP-odd)"},
217  tCgg  == { TeX -> Subscript[OverTilde[c],gg], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxh, OrderBlock -> 12, Description -> "fsG-fsG-h interactions (CP-odd)"},
218  tCaa  == { TeX -> Subscript[OverTilde[c],aa], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxh, OrderBlock -> 13, Description -> "fsA-fsA-h interactions (CP-odd)"},
219  tCza  == { TeX -> Subscript[OverTilde[c],za], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxh, OrderBlock -> 14, Description -> "fsZ-fsA-h interactions (CP-odd)"},
220  tCzz  == { TeX -> Subscript[OverTilde[c],zz], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxh, OrderBlock -> 15, Description -> "fsZ-fsZ-h interactions (CP-odd)"},
221
222  dYu == { TeX -> Subscript["\[Delta]Y",u], ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {dYu[i_,j_]->0.1}, InteractionOrder -> {QNP, 1},
223    BlockName -> BCxdYu, Description -> "Norm of Huu couplings"},
224  dYd == { TeX -> Subscript["\[Delta]Y",d], ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {dYd[i_,j_]->0.1}, InteractionOrder -> {QNP, 1},
225    BlockName -> BCxdYd, Description -> "Norm of Hdd couplings"},
226  dYe == { TeX -> Subscript["\[Delta]Y",e], ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {dYe[i_,j_]->0.1}, InteractionOrder -> {QNP, 1},
227    BlockName -> BCxdYe, Description -> "Norm of Hee couplings"},
228  Su  == { TeX -> Subscript[S,"\[Phi]u"], ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Su[i_,j_]->0.707107}, BlockName -> BCxSu, Description -> "Sin(phase) of Huu couplings"},
229  Sd  == { TeX -> Subscript[S,"\[Phi]d"], ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Sd[i_,j_]->0.707107}, BlockName -> BCxSd, Description -> "Sin(phase) of Hdd couplings"},
230  Se  == { TeX -> Subscript[S,"\[Phi]e"], ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Se[i_,j_]->0.707107}, BlockName -> BCxSe, Description -> "Sin(phase) of Hee couplings"},
231  Cu  == { TeX -> Subscript[C,"\[Phi]u"], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, Value -> {Cu[i_,j_]->Sqrt[1-Su[i,j]^2]}, Description -> "Cos(phase) of Huu couplings"},
232  Cd  == { TeX -> Subscript[C,"\[Phi]d"], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, Value -> {Cd[i_,j_]->Sqrt[1-Sd[i,j]^2]}, Description -> "Cos(phase) of Huu couplings"},
233  Ce  == { TeX -> Subscript[C,"\[Phi]e"], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, Value -> {Ce[i_,j_]->Sqrt[1-Se[i,j]^2]}, Description -> "Cos(phase) of Huu couplings"},
234
235  (* Lhvff parameters *)
236  RdGLhzu ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdGLhzu[i_,j_]->0.1}, BlockName ->   BCxdGLhzu, Description -> "h-Z-u vertex interactions (real part, L)"},
237  IdGLhzu ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdGLhzu[i_,j_]->0.0}, BlockName -> IMBCxdGLhzu, Description -> "h-Z-u vertex interactions (im part, L)"  },
238  RdGLhzd ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdGLhzd[i_,j_]->0.1}, BlockName ->   BCxdGLhzd, Description -> "h-Z-d vertex interactions (real part, L)"},
239  IdGLhzd ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdGLhzd[i_,j_]->0.0}, BlockName -> IMBCxdGLhzd, Description -> "h-Z-d vertex interactions (im part, L)"  },
240  RdGLhze ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdGLhze[i_,j_]->0.1}, BlockName ->   BCxdGLhze, Description -> "h-Z-l vertex interactions (real part, L)"},
241  IdGLhze ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdGLhze[i_,j_]->0.0}, BlockName -> IMBCxdGLhze, Description -> "h-Z-l vertex interactions (im part, L)"},
242  RdGLhzv ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdGLhzv[i_,j_]->0.1}, BlockName ->   BCxdGLhzv, Description -> "h-Z-v vertex interactions (real part, L)"},
243  IdGLhzv ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdGLhzv[i_,j_]->0.0}, BlockName -> IMBCxdGLhzv, Description -> "h-Z-v vertex interactions (im part, L)"},
244  RdGRhzu ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdGRhzu[i_,j_]->0.1}, BlockName ->   BCxdGRhzu, Description -> "h-Z-u vertex interactions (real part, R)"},
245  IdGRhzu ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdGRhzu[i_,j_]->0.0}, BlockName -> IMBCxdGRhzu, Description -> "h-Z-u vertex interactions (im part, R)"  },
246  RdGRhzd ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdGRhzd[i_,j_]->0.1}, BlockName ->   BCxdGRhzd, Description -> "h-Z-d vertex interactions (real part, R)"},
247  IdGRhzd ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdGRhzd[i_,j_]->0.0}, BlockName -> IMBCxdGRhzd, Description -> "h-Z-d vertex interactions (im part, R)"  },
248  RdGRhze ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdGRhze[i_,j_]->0.1}, BlockName ->   BCxdGRhze, Description -> "h-Z-l vertex interactions (real part, R)"},
249  IdGRhze ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdGRhze[i_,j_]->0.0}, BlockName -> IMBCxdGRhze, Description -> "h-Z-l vertex interactions (im part, R)"},
250  RdGLhwq ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdGLhwq[i_,j_]->0.1}, BlockName ->   BCxdGLhWq, Description -> "h-W-q vertex interactions (real part, L)"},
251  IdGLhwq ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdGLhwq[i_,j_]->0.0}, BlockName -> IMBCxdGLhWq, Description -> "h-W-q vertex interactions (im part, L)"  },
252  RdGRhwq ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdGRhwq[i_,j_]->0.1}, BlockName ->   BCxdGRhWq, Description -> "h-W-q vertex interactions (real part, R)"},
253  IdGRhwq ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdGRhwq[i_,j_]->0.0}, BlockName -> IMBCxdGRhWq, Description -> "h-W-q vertex interactions (im part, R)"  },
254  RdGLhwl ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdGLhwl[i_,j_]->0.1}, BlockName ->   BCxdGLhWl, Description -> "h-W-l vertex interactions (real part, L)"},
255  IdGLhwl ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdGLhwl[i_,j_]->0.0}, BlockName -> IMBCxdGLhWl, Description -> "h-W-l vertex interactions (im part, L)"  },
256
257  dGLhzu ==  { TeX -> Subsuperscript["\[Delta]g", L, hZu], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
258    Value -> {dGLhzu[i_,j_]->RdGLhzu[i,j] + I IdGLhzu[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-Z-u vertex interactions (L)"},
259  dGLhzd ==  { TeX -> Subsuperscript["\[Delta]g", L, hZd], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
260    Value -> {dGLhzd[i_,j_]->RdGLhzd[i,j] + I IdGLhzd[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-Z-d vertex interactions (L)"},
261  dGLhze ==  { TeX -> Subsuperscript["\[Delta]g", L, hZl], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
262    Value -> {dGLhze[i_,j_]->RdGLhze[i,j] + I IdGLhze[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-Z-l vertex interactions (L)"},
263  dGLhzv ==  { TeX -> Subsuperscript["\[Delta]g", L, hZv], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
264    Value -> {dGLhzv[i_,j_]->RdGLhzv[i,j] + I IdGLhzv[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-Z-v vertex interactions (L)"},
265  dGRhzu ==  { TeX -> Subsuperscript["\[Delta]g", R, hZu], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
266    Value -> {dGRhzu[i_,j_]->RdGRhzu[i,j] + I IdGRhzu[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-Z-u vertex interactions (R)"},
267  dGRhzd ==  { TeX -> Subsuperscript["\[Delta]g", R, hZd], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
268    Value -> {dGRhzd[i_,j_]->RdGRhzd[i,j] + I IdGRhzd[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-Z-d vertex interactions (R)"},
269  dGRhze ==  { TeX -> Subsuperscript["\[Delta]g", R, hZl], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
270    Value -> {dGRhze[i_,j_]->RdGRhze[i,j] + I IdGRhze[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-Z-l vertex interactions (R)"},
271  dGLhwq ==  { TeX -> Subsuperscript["\[Delta]g", L, hWq], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
272    Value -> {dGLhwq[i_,j_]->RdGLhwq[i,j] + I IdGLhwq[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-W-q vertex interactions (L)"},
273  dGRhwq ==  { TeX -> Subsuperscript["\[Delta]g", R, hWq], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
274    Value -> {dGRhwq[i_,j_]->RdGRhwq[i,j] + I IdGRhwq[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-W-q vertex interactions (R)"},
275  dGLhwl ==  { TeX -> Subsuperscript["\[Delta]g", L, hWl], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
276    Value -> {dGLhwl[i_,j_]->RdGLhwl[i,j] + I IdGLhwl[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-W-l vertex interactions (L)"},
277
278
279  (* Lhdvff parameters *)
280  Rdhgu ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdhgu[i_,j_]->0.1}, BlockName ->   BCxdhgu, Description -> "h-G-u dipole interactions (real part)"},
281  Idhgu ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idhgu[i_,j_]->0.0}, BlockName -> IMBCxdhgu, Description -> "h-G-u dipole interactions (im part)"  },
282  Rdhgd ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdhgd[i_,j_]->0.1}, BlockName ->   BCxdhgd, Description -> "h-G-d dipole interactions (real part)"},
283  Idhgd ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idhgd[i_,j_]->0.0}, BlockName -> IMBCxdhgd, Description -> "h-G-d dipole interactions (im part)"  },
284  Rdhau ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdhau[i_,j_]->0.1}, BlockName ->   BCxdhau, Description -> "h-A-u dipole interactions (real part)"},
285  Idhau ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idhau[i_,j_]->0.0}, BlockName -> IMBCxdhau, Description -> "h-A-u dipole interactions (im part)"  },
286  Rdhad ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdhad[i_,j_]->0.1}, BlockName ->   BCxdhad, Description -> "h-A-d dipole interactions (real part)"},
287  Idhad ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idhad[i_,j_]->0.0}, BlockName -> IMBCxdhad, Description -> "h-A-d dipole interactions (im part)"  },
288  Rdhae ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdhae[i_,j_]->0.1}, BlockName ->   BCxdhae, Description -> "h-A-e dipole interactions (real part)"},
289  Idhae ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idhae[i_,j_]->0.0}, BlockName -> IMBCxdhae, Description -> "h-A-e dipole interactions (im part)"  },
290  Rdhzu ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdhzu[i_,j_]->0.1}, BlockName ->   BCxdhzu, Description -> "h-Z-u dipole interactions (real part)"},
291  Idhzu ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idhzu[i_,j_]->0.0}, BlockName -> IMBCxdhzu, Description -> "h-Z-u dipole interactions (im part)"  },
292  Rdhzd ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdhzd[i_,j_]->0.1}, BlockName ->   BCxdhzd, Description -> "h-Z-d dipole interactions (real part)"},
293  Idhzd ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idhzd[i_,j_]->0.0}, BlockName -> IMBCxdhzd, Description -> "h-Z-d dipole interactions (im part)"  },
294  Rdhze ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdhze[i_,j_]->0.1}, BlockName ->   BCxdhze, Description -> "h-Z-e dipole interactions (real part)"},
295  Idhze ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idhze[i_,j_]->0.0}, BlockName -> IMBCxdhze, Description -> "h-Z-e dipole interactions (im part)"  },
296  Rdhwd ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdhwd[i_,j_]->0.1}, BlockName ->   BCxdhwd, Description -> "h-W-d dipole interactions (real part)"},
297  Idhwd ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idhwd[i_,j_]->0.0}, BlockName -> IMBCxdhwd, Description -> "h-W-d dipole interactions (im part)"  },
298  Rdhwu ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdhwu[i_,j_]->0.1}, BlockName ->   BCxdhwu, Description -> "h-W-u dipole interactions (real part)"},
299  Idhwu ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idhwu[i_,j_]->0.0}, BlockName -> IMBCxdhwu, Description -> "h-W-u dipole interactions (im part)"  },
300  Rdhwl ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rdhwl[i_,j_]->0.1}, BlockName ->   BCxdhwl, Description -> "h-W-l dipole interactions (real part)"},
301  Idhwl ==  { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Idhwl[i_,j_]->0.0}, BlockName -> IMBCxdhwl, Description -> "h-W-l dipole interactions (im part)"  },
302  Rtdhgu == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rtdhgu[i_,j_]->0.1}, BlockName ->   BCxtdhgu, Description -> "h-G-u dipole interactions (real part)"},
303  Itdhgu == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Itdhgu[i_,j_]->0.0}, BlockName -> IMBCxtdhgu, Description -> "h-G-u dipole interactions (im part, CP-odd)  "},
304  Rtdhgd == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rtdhgd[i_,j_]->0.1}, BlockName ->   BCxtdhgd, Description -> "h-G-d dipole interactions (real part, CP-odd)"},
305  Itdhgd == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Itdhgd[i_,j_]->0.0}, BlockName -> IMBCxtdhgd, Description -> "h-G-d dipole interactions (im part, CP-odd)"  },
306  Rtdhau == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rtdhau[i_,j_]->0.1}, BlockName ->   BCxtdhau, Description -> "h-A-u dipole interactions (real part, CP-odd)"},
307  Itdhau == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Itdhau[i_,j_]->0.0}, BlockName -> IMBCxtdhau, Description -> "h-A-u dipole interactions (im part, CP-odd)"  },
308  Rtdhad == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rtdhad[i_,j_]->0.1}, BlockName ->   BCxtdhad, Description -> "h-A-d dipole interactions (real part, CP-odd)"},
309  Itdhad == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Itdhad[i_,j_]->0.0}, BlockName -> IMBCxtdhad, Description -> "h-A-d dipole interactions (im part, CP-odd)"  },
310  Rtdhae == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rtdhae[i_,j_]->0.1}, BlockName ->   BCxtdhae, Description -> "h-A-e dipole interactions (real part, CP-odd)"},
311  Itdhae == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Itdhae[i_,j_]->0.0}, BlockName -> IMBCxtdhae, Description -> "h-A-e dipole interactions (im part, CP-odd)"  },
312  Rtdhzu == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rtdhzu[i_,j_]->0.1}, BlockName ->   BCxtdhzu, Description -> "h-Z-u dipole interactions (real part, CP-odd)"},
313  Itdhzu == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Itdhzu[i_,j_]->0.0}, BlockName -> IMBCxtdhzu, Description -> "h-Z-u dipole interactions (im part, CP-odd)"  },
314  Rtdhzd == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rtdhzd[i_,j_]->0.1}, BlockName ->   BCxtdhzd, Description -> "h-Z-d dipole interactions (real part, CP-odd)"},
315  Itdhzd == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Itdhzd[i_,j_]->0.0}, BlockName -> IMBCxtdhzd, Description -> "h-Z-d dipole interactions (im part, CP-odd)"  },
316  Rtdhze == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Rtdhze[i_,j_]->0.1}, BlockName ->   BCxtdhze, Description -> "h-Z-e dipole interactions (real part, CP-odd)"},
317  Itdhze == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {Itdhze[i_,j_]->0.0}, BlockName -> IMBCxtdhze, Description -> "h-Z-e dipole interactions (im part, CP-odd)"  },
318
319  dhgu ==  { TeX -> Subscript["d", hgu], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
320    Value -> {dhgu[i_,j_]->Rdhgu[i,j] + I Idhgu[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-G-u dipole interactions"},
321  dhgd ==  { TeX -> Subscript["d", hgd], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
322    Value -> {dhgd[i_,j_]->Rdhgd[i,j] + I Idhgd[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-G-d dipole interactions"},
323  dhau ==  { TeX -> Subscript["d", hau], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
324    Value -> {dhau[i_,j_]->Rdhau[i,j] + I Idhau[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-A-u dipole interactions"},
325  dhad ==  { TeX -> Subscript["d", had], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
326    Value -> {dhad[i_,j_]->Rdhad[i,j] + I Idhad[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-A-d dipole interactions"},
327  dhae ==  { TeX -> Subscript["d", hae], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
328    Value -> {dhae[i_,j_]->Rdhae[i,j] + I Idhae[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-A-e dipole interactions"},
329  dhzu ==  { TeX -> Subscript["d", hzu], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
330    Value -> {dhzu[i_,j_]->Rdhzu[i,j] + I Idhzu[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-Z-u dipole interactions"},
331  dhzd ==  { TeX -> Subscript["d", hzd], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
332    Value -> {dhzd[i_,j_]->Rdhzd[i,j] + I Idhzd[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-Z-d dipole interactions"},
333  dhze ==  { TeX -> Subscript["d", hze], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
334    Value -> {dhze[i_,j_]->Rdhze[i,j] + I Idhze[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-Z-e dipole interactions"},
335  dhwu ==  { TeX -> Subscript["d", hwu], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
336    Value -> {dhwu[i_,j_]->Rdhwu[i,j] + I Idhwu[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-w-u dipole interactions"},
337  dhwd ==  { TeX -> Subscript["d", hwd], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
338    Value -> {dhwd[i_,j_]->Rdhwd[i,j] + I Idhwd[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-w-d dipole interactions"},
339  dhwl ==  { TeX -> Subscript["d", hwl], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
340    Value -> {dhwl[i_,j_]->Rdhwl[i,j] + I Idhwl[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-W-l dipole interactions"},
341
342  tdhgu ==  { TeX -> Subscript[OverTilde["d"], hgu], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
343    Value -> {tdhgu[i_,j_]->Rtdhgu[i,j] + I Itdhgu[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-G-u dipole interactions (CP-odd)"},
344  tdhgd ==  { TeX -> Subscript[OverTilde["d"], hgd], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
345    Value -> {tdhgd[i_,j_]->Rtdhgd[i,j] + I Itdhgd[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-G-d dipole interactions (CP-odd)"},
346  tdhau ==  { TeX -> Subscript[OverTilde["d"], hau], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
347    Value -> {tdhau[i_,j_]->Rtdhau[i,j] + I Itdhau[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-A-u dipole interactions (CP-odd)"},
348  tdhad ==  { TeX -> Subscript[OverTilde["d"], had], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
349    Value -> {tdhad[i_,j_]->Rtdhad[i,j] + I Itdhad[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-A-d dipole interactions (CP-odd)"},
350  tdhae ==  { TeX -> Subscript[OverTilde["d"], hae], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
351    Value -> {tdhae[i_,j_]->Rtdhae[i,j] + I Itdhae[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-A-e dipole interactions (CP-odd)"},
352  tdhzu ==  { TeX -> Subscript[OverTilde["d"], hzu], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
353    Value -> {tdhzu[i_,j_]->Rtdhzu[i,j] + I Itdhzu[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-Z-u dipole interactions (CP-odd)"},
354  tdhzd ==  { TeX -> Subscript[OverTilde["d"], hzd], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
355    Value -> {tdhzd[i_,j_]->Rtdhzd[i,j] + I Itdhzd[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-Z-d dipole interactions (CP-odd)"},
356  tdhze ==  { TeX -> Subscript[OverTilde["d"], hze], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter->True,
357    Value -> {tdhze[i_,j_]->Rtdhze[i,j] + I Itdhze[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "h-Z-e dipole interactions (CP-odd)"},
358
359  (* Lhself parameters *)
360  dL3 == { TeX -> Subscript["\[Delta]\[Lambda]",3], ParameterType -> External, Value -> 0.1, InteractionOrder -> {{QED,1}, {QNP, 1}}, BlockName -> BCxhself,  OrderBlock -> 1, Description -> "Trilinear Higgs self-coupling"},
361  dL4 == { TeX -> Subscript["\[Delta]\[Lambda]",4], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxhself,  OrderBlock -> 2, Description -> "Quartic Higgs self-coupling"},
362
363  (* Lhh parameters *)
364  dCw2  == { TeX -> Subsuperscript["\[Delta]c",w,"(2)"],    ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxhh, OrderBlock -> 1, Description -> "W-W-h-h interactions"},
365  dCz2  == { TeX -> Subsuperscript["\[Delta]c",z,"(2)"],    ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxhh, OrderBlock -> 2, Description -> "Z-Z-h-h interactions"},
366  Cww2  == { TeX -> Subsuperscript[c,ww,"(2)"],             ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxhh, OrderBlock -> 3, Description -> "fsW-fsW-h-h interactions"},
367  Cgg2  == { TeX -> Subsuperscript[c,gg,"(2)"],             ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxhh, OrderBlock -> 4, Description -> "fsG-fsG-h-h interactions"},
368  Caa2  == { TeX -> Subsuperscript[c,aa,"(2)"],             ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxhh, OrderBlock -> 5, Description -> "fsA-fsA-h-h interactions"},
369  Cza2  == { TeX -> Subsuperscript[c,za,"(2)"],             ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxhh, OrderBlock -> 6, Description -> "fsZ-fsA-h-h interactions"},
370  Czz2  == { TeX -> Subsuperscript[c,zz,"(2)"],             ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxhh, OrderBlock -> 7, Description -> "fsZ-fsZ-h-h interactions"},
371  Cwbx2 == { TeX -> Subsuperscript[c,w,Square[""],"(2)"],   ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxhh, OrderBlock -> 8, Description -> "fsW-W-h-h interactions"},
372  Czbx2 == { TeX -> Subsuperscript[c,z,Square[""],"(2)"],   ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxhh, OrderBlock -> 9, Description -> "fsZ-Z-h-h interactions"},
373  Cabx2 == { TeX -> Subsuperscript[c,a,Square[""],"(2)"],   ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxhh, OrderBlock -> 10,Description -> "fsA-Z-h-h interactions"},
374  tCww2  == { TeX -> Subsuperscript[OverTilde[c],ww,"(2)"], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxhh, OrderBlock -> 11, Description -> "fsW-fsW-h-h interactions (CP-odd)"},
375  tCgg2  == { TeX -> Subsuperscript[OverTilde[c],gg,"(2)"], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxhh, OrderBlock -> 12, Description -> "fsG-fsG-h-h interactions (CP-odd)"},
376  tCaa2  == { TeX -> Subsuperscript[OverTilde[c],aa,"(2)"], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxhh, OrderBlock -> 13, Description -> "fsA-fsA-h-h interactions (CP-odd)"},
377  tCza2  == { TeX -> Subsuperscript[OverTilde[c],za,"(2)"], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxhh, OrderBlock -> 14, Description -> "fsZ-fsA-h-h interactions (CP-odd)"},
378  tCzz2  == { TeX -> Subsuperscript[OverTilde[c],zz,"(2)"], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCxhh, OrderBlock -> 15, Description -> "fsZ-fsZ-h-h interactions (CP-odd)"},
379
380  RdYu2 == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdYu2[i_,j_]->0.1}, BlockName ->   BCxY2u, Description -> "HHuu couplings (real part)"},
381  RdYd2 == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdYd2[i_,j_]->0.1}, BlockName ->   BCxY2d, Description -> "HHdd couplings (real part)"},
382  RdYe2 == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {RdYe2[i_,j_]->0.1}, BlockName ->   BCxY2e, Description -> "HHee couplings (real part)"},
383  IdYu2 == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdYu2[i_,j_]->0.0}, BlockName -> IMBCxY2u, Description -> "HHuu couplings (im part)"},
384  IdYd2 == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdYd2[i_,j_]->0.0}, BlockName -> IMBCxY2d, Description -> "HHdd couplings (im part)"},
385  IdYe2 == { ParameterType -> External, Indices -> {Index[Generation], Index[Generation]}, Value -> {IdYe2[i_,j_]->0.0}, BlockName -> IMBCxY2e, Description -> "HHee couplings (im part)"},
386
387  dYu2 == { TeX -> Subsuperscript["\[Delta]Y",u,"(2)"], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter -> True,
388    Value -> {dYu2[i_,j_]->RdYu2[i,j]+I*IdYu2[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "HHuu couplings"},
389  dYd2 == { TeX -> Subsuperscript["\[Delta]Y",d,"(2)"], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter -> True,
390    Value -> {dYd2[i_,j_]->RdYd2[i,j]+I*IdYd2[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "HHdd couplings"},
391  dYe2 == { TeX -> Subsuperscript["\[Delta]Y",e,"(2)"], ParameterType -> Internal, Indices -> {Index[Generation], Index[Generation]}, ComplexParameter -> True,
392    Value -> {dYe2[i_,j_]->RdYe2[i,j]+I*IdYe2[i,j]}, InteractionOrder -> {QNP, 1}, Description -> "HHee couplings"},
393
394  (* Four-fermion interactions *)
395  cll1122 == { TeX -> Subscript[c,ll,1122],    ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCx4f,  OrderBlock -> 1, Description -> "Four fermion interaction strength"},
396  cpuu3333== { TeX -> Subscript["c'",uu,3333], ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCx4f,  OrderBlock -> 2, Description -> "Four fermion interaction strength"},
397  cll1221 == { TeX -> Subscript[c,ll,1221],    ParameterType -> External, Value -> 0.1, InteractionOrder -> {QNP, 1}, BlockName -> BCx4f,  OrderBlock -> 3, Description -> "Four fermion interaction strength"},
398
399  (* loop-induced SM couplings *)
400  gHgg == { ParameterType -> Internal,
401            Value -> 1/(4*Pi)/(3*Pi),
402            TeX -> Subscript[g,Hgg],
403            Description -> "loop-induced SM Hgg coupling"},                     
404  gHaa == { ParameterType -> Internal,
405            Value -> -1/(4*Pi)/(Pi)*(47/18),
406            TeX -> Subscript[g,Haa],
407            Description -> "loop-induced SM Haa coupling"},             
408  gHza == { ParameterType -> Internal,
409            Value -> -Sqrt[ee^2/(4*Pi)*Gf*MZ^2/(8*Sqrt[2]*Pi)]*(94*cw^2-13)/(9*Pi)/ee^2,
410            TeX -> Subscript[g,Hza],
411            Description -> "loop-induced SM Hza coupling"}
412
413};
414
415
416(* ************************** *)
417(* *****   Lagrangian   ***** *)
418(* ************************** *)
419
420(* eq. (3.3) in the Higgs basis note *)
421LMass:= 2 dM gw^2 vev^2/4 W[mu] Wbar[mu];
422
423(* eq. (3.4) *)
424Lvertex := Block[{LZ,LW},
425  LZ =  Sqrt[gw^2+g1^2] Z[mu] ( uqbar.Ga[mu].dGLzu.ProjM.uq + dqbar.Ga[mu].dGLzd.ProjM.dq + lbar.Ga[mu].dGLze.ProjM.l + vlbar.Ga[mu].dGLzv.ProjM.vl +
426                                uqbar.Ga[mu].dGRzu.ProjP.uq + dqbar.Ga[mu].dGRzd.ProjP.dq + lbar.Ga[mu].dGRze.ProjP.l );
427  LW= gw/Sqrt[2] W[mu] (uqbar.Ga[mu].dGLwq.ProjM.dq + uqbar.Ga[mu].dGRwq.ProjP.dq + vlbar.Ga[mu].dGLwl.ProjM.l);
428  LZ+LW+HC[LW]
429];
430
431(* eq. (3.5) for the dipole terms *)
432Ldip:=Block[{LG,LA,LZ,LW},
433  LG = gs FS[G,mu,nu,a] (uqbar.Sig[mu,nu].T[a].dgu.uq + dqbar.Sig[mu,nu].T[a].dgd.dq) +
434       gs Dual[FS][G,mu,nu,a] (uqbar.Sig[mu,nu].T[a].tdgu.uq + dqbar.Sig[mu,nu].T[a].tdgd.dq);
435  LA = ee FS[A,mu,nu] (uqbar.Sig[mu,nu].dau.uq + dqbar.Sig[mu,nu].dad.dq + lbar.Sig[mu,nu].dae.l) +
436       ee Dual[FS][A,mu,nu] (uqbar.Sig[mu,nu].tdau.uq + dqbar.Sig[mu,nu].tdad.dq + lbar.Sig[mu,nu].tdae.l);
437  LZ = Sqrt[gw^2+g1^2] FS[Z,mu,nu] (uqbar.Sig[mu,nu].dzu.uq + dqbar.Sig[mu,nu].dzd.dq + lbar.Sig[mu,nu].dze.l) +
438       Sqrt[gw^2+g1^2] Dual[FS][Z,mu,nu] (uqbar.Sig[mu,nu].tdzu.uq + dqbar.Sig[mu,nu].tdzd.dq + lbar.Sig[mu,nu].tdze.l);
439  LW = Sqrt[2] gw (FS[Wbar,mu,nu] dqbar.Sig[mu,nu].dwu.ProjP.uq + FS[W,mu,nu] uqbar.Sig[mu,nu].dwd.ProjP.dq + FS[W,mu,nu] vlbar.Sig[mu,nu].dwl.ProjP.l);
440  -1/(4 vev) (LG+LA+LZ+LW+HC[LW])
441];
442
443(* eqs. (3.6) and (3.7) for the gauge boson self-interactions *)
444Ltgc := I ee ( dKa FS[A,mu,nu] W[mu] Wbar[nu]  + tKa Dual[FS][A,mu,nu] W[mu] Wbar[nu]  ) +
445        I gw cw ( dG1z (FS[W,mu,nu] Wbar[mu] - FS[Wbar,mu,nu] W[mu] ) Z[nu] + dKz FS[Z,mu,nu] W[mu] Wbar[nu]  + tKz Dual[FS][Z,mu,nu] W[mu] Wbar[nu] ) +
446        I ee/MW^2 ( La FS[W,mu,nu] FS[Wbar,nu,rho] FS[A,rho,mu] + tLa FS[W,mu,nu] FS[Wbar,nu,rho] Dual[FS][A,rho,mu] ) +
447        I gw cw/MW^2 ( Lz  FS[W,mu,nu] FS[Wbar,nu,rho] FS[Z,rho,mu] + tLz FS[W,mu,nu] FS[Wbar,nu,rho] Dual[FS][Z,rho,mu] ) +
448        gs^3/vev^2 f[a,b,c] ( C3g FS[G,mu,nu,a] FS[G,nu,rho,b] FS[G,rho,mu,c] + tC3g Dual[FS][G,mu,nu,a] FS[G,nu,rho,b] FS[G,rho,mu,c] );
449
450Lqgc := gw^2/2    dGw4    ( W[mu] W[mu] Wbar[nu] Wbar[nu] - W[mu] Wbar[mu] W[nu] Wbar[nu] ) +
451        gw^2 cw^2 dGw2z2  ( W[mu] Z[mu] Wbar[nu] Z[nu]    - W[mu] Wbar[mu] Z[nu] Z[nu]    ) +
452        gw^2 cw sw dGw2za ( W[mu] Z[mu] Wbar[nu] A[nu]    + W[mu] A[mu] Wbar[nu] Z[nu] - 2 W[mu] Wbar[mu] Z[nu] A[nu] ) -
453        gw^2/2/MW^2    Lw4   ( FS[W,mu,nu] FS[Wbar,nu,rho] - FS[Wbar,mu,nu] FS[W,nu,rho] ) ( W[mu] Wbar[rho] - Wbar[mu] W[rho] ) -
454        gw^2 cw^2/MW^2 Lw2z2 ( W[mu] ( FS[Z,mu,nu] FS[Wbar,nu,rho] - FS[Wbar,mu,nu] FS[Z,nu,rho] ) Z[rho] + Wbar[mu] ( FS[Z,mu,nu] FS[W,nu,rho] - FS[W,mu,nu] FS[Z,nu,rho] ) Z[rho] ) -
455        ee^2/MW^2      Lw2a2 ( W[mu] ( FS[A,mu,nu] FS[Wbar,nu,rho] - FS[Wbar,mu,nu] FS[A,nu,rho] ) A[rho] + Wbar[mu] ( FS[A,mu,nu] FS[W,nu,rho] - FS[W,mu,nu] FS[A,nu,rho] ) A[rho] ) -
456        ee gw cw/MW^2  Lw2az ( W[mu] ( FS[A,mu,nu] FS[Wbar,nu,rho] - FS[Wbar,mu,nu] FS[A,nu,rho] ) Z[rho] + Wbar[mu] ( FS[A,mu,nu] FS[W,nu,rho] - FS[W,mu,nu] FS[A,nu,rho] ) Z[rho] ) -
457        ee gw cw/MW^2  Lw2za ( W[mu] ( FS[Z,mu,nu] FS[Wbar,nu,rho] - FS[Wbar,mu,nu] FS[Z,nu,rho] ) A[rho] + Wbar[mu] ( FS[Z,mu,nu] FS[W,nu,rho] - FS[W,mu,nu] FS[Z,nu,rho] ) A[rho] ) +
458        3 gs^3/vev^2   C4g f[a,b,c] f[c,d,e] FS[G,mu,nu,a] FS[G,nu,rho,b] G[rho,d] G[mu,e] -
459        gw^2/2/MW^2    tLw4   ( FS[W,mu,nu] Dual[FS][Wbar,nu,rho] - FS[Wbar,mu,nu] Dual[FS][W,nu,rho] ) ( W[mu] Wbar[rho] - Wbar[mu] W[rho] ) -
460        gw^2 cw^2/MW^2 tLw2z2 ( W[mu] ( FS[Z,mu,nu] Dual[FS][Wbar,nu,rho] - FS[Wbar,mu,nu] Dual[FS][Z,nu,rho] ) Z[rho] + Wbar[mu] ( FS[Z,mu,nu] FS[W,nu,rho] - FS[W,mu,nu] FS[Z,nu,rho] ) Z[rho] ) -
461        ee^2/MW^2      tLw2a2 ( W[mu] ( FS[A,mu,nu] Dual[FS][Wbar,nu,rho] - FS[Wbar,mu,nu] Dual[FS][A,nu,rho] ) A[rho] + Wbar[mu] ( FS[A,mu,nu] FS[W,nu,rho] - FS[W,mu,nu] FS[A,nu,rho] ) A[rho] ) -
462        ee gw cw/MW^2  tLw2az ( W[mu] ( FS[A,mu,nu] Dual[FS][Wbar,nu,rho] - FS[Wbar,mu,nu] Dual[FS][A,nu,rho] ) Z[rho] + Wbar[mu] ( FS[A,mu,nu] FS[W,nu,rho] - FS[W,mu,nu] FS[A,nu,rho] ) Z[rho] ) -
463        ee gw cw/MW^2  tLw2za ( W[mu] ( FS[Z,mu,nu] Dual[FS][Wbar,nu,rho] - FS[Wbar,mu,nu] Dual[FS][Z,nu,rho] ) A[rho] + Wbar[mu] ( FS[Z,mu,nu] FS[W,nu,rho] - FS[W,mu,nu] FS[Z,nu,rho] ) A[rho] ) +
464        3 gs^3/vev^2   tC4g f[a,b,c] f[c,d,e] Dual[FS][G,mu,nu,a] FS[G,nu,rho,b] G[rho,d] G[mu,e];
465
466(* eq.(3.8) for the 3-point single Higgs interactions *)
467Lh := Block[{Lvv,Lff},
468  Lvv = 2 dCw MW^2 Wbar[mu] W[mu] + dCz MZ^2 Z[mu] Z[mu] +
469           gs^2/4     gHgg FS[G,mu,nu,a]  FS[G,mu,nu,a] +
470           gs^2/4      Cgg FS[G,mu,nu,a]  FS[G,mu,nu,a] + gs^2/4      tCgg FS[G,mu,nu,a]  Dual[FS][G,mu,nu,a] +
471           gw^2/2      Cww FS[Wbar,mu,nu] FS[W,mu,nu]   + gw^2/2      tCww FS[Wbar,mu,nu] Dual[FS][W,mu,nu]  +
472           ee^2/4     gHaa FS[A,mu,nu]    FS[A,mu,nu]   +
473           ee^2/4      Caa FS[A,mu,nu]    FS[A,mu,nu]   + ee^2/4      tCaa FS[A,mu,nu]    Dual[FS][A,mu,nu] +
474           ee^2/2     gHza FS[Z,mu,nu]    FS[A,mu,nu]   +
475           ee gw/2/cw  Cza FS[Z,mu,nu]    FS[A,mu,nu]   + ee gw/2/cw  tCza FS[Z,mu,nu]    Dual[FS][A,mu,nu] +
476           gw^2/4/cw^2 Czz FS[Z,mu,nu]    FS[Z,mu,nu]   + gw^2/4/cw^2 tCzz FS[Z,mu,nu]    Dual[FS][Z,mu,nu] +
477           gw^2 Cwbx (Wbar[mu] del[FS[W,mu,nu],nu] + W[mu] del[FS[Wbar,mu,nu],nu]) + g1 gw Cabx Z[mu] del[FS[A,mu,nu],nu] + gw^2 Czbx Z[mu] del[FS[Z,mu,nu],nu];
478  Lff = - ExpandIndices[Plus@@Flatten[Table[Sqrt[Mu[ff] Mu[ffp]] dYu[ff,ffp] (Cu[ff,ffp] IndexDelta[sp1,sp2]-I*Su[ff,ffp] Ga[5,sp1,sp2]) uqbar[sp1,ff,cc].uq[sp2,ffp,cc], {ffp,1,3},{ff,1,3}]]] -
479          ExpandIndices[Plus@@Flatten[Table[Sqrt[Md[ff] Md[ffp]] dYd[ff,ffp] (Cd[ff,ffp] IndexDelta[sp1,sp2]-I*Sd[ff,ffp] Ga[5,sp1,sp2]) dqbar[sp1,ff,cc].dq[sp2,ffp,cc], {ffp,1,3},{ff,1,3}]]] -
480          ExpandIndices[Plus@@Flatten[Table[Sqrt[Ml[ff] Ml[ffp]] dYe[ff,ffp] (Ce[ff,ffp] IndexDelta[sp1,sp2]-I*Se[ff,ffp] Ga[5,sp1,sp2]) lbar[sp1,ff].l[sp2,ffp], {ffp,1,3},{ff,1,3}]]];
481  H/vev (Lvv + Lff)
482];
483
484(* eq (3.9) for the hvff interactions *)
485Lhvff := Block[{LW,LZ},
486  LW = Sqrt[2] gw H/vev W[mu] (uqbar.Ga[mu].dGLhwq.ProjM.dq + uqbar.Ga[mu].dGRhwq.ProjP.dq + vlbar.Ga[mu].dGLhwl.ProjM.l);
487  LZ = 2 Sqrt[gw^2+g1^2] H/vev Z[mu] ( uqbar.Ga[mu].dGLhzu.ProjM.uq + dqbar.Ga[mu].dGLhzd.ProjM.dq + lbar.Ga[mu].dGLhze.ProjM.l + vlbar.Ga[mu].dGLhzv.ProjM.vl +
488                                       uqbar.Ga[mu].dGRhzu.ProjP.uq + dqbar.Ga[mu].dGRhzd.ProjP.dq + lbar.Ga[mu].dGRhze.ProjP.l );
489  LZ + LW + HC[LW]
490];
491
492(* eq. (3.11) for the hvff dipole terms *)
493Lhdvff:=Block[{LG,LA,LZ,LW},
494  LG = gs FS[G,mu,nu,a] (uqbar.Sig[mu,nu].T[a].dhgu.uq + dqbar.Sig[mu,nu].T[a].dhgd.dq) +
495       gs Dual[FS][G,mu,nu,a] (uqbar.Sig[mu,nu].T[a].tdhgu.uq + dqbar.Sig[mu,nu].T[a].tdhgd.dq);
496  LA = ee FS[A,mu,nu] (uqbar.Sig[mu,nu].dhau.uq + dqbar.Sig[mu,nu].dhad.dq + lbar.Sig[mu,nu].dhae.l) +
497       ee Dual[FS][A,mu,nu] (uqbar.Sig[mu,nu].tdhau.uq + dqbar.Sig[mu,nu].tdhad.dq + lbar.Sig[mu,nu].tdhae.l);
498  LZ = Sqrt[gw^2+g1^2] FS[Z,mu,nu] (uqbar.Sig[mu,nu].dhzu.uq + dqbar.Sig[mu,nu].dhzd.dq + lbar.Sig[mu,nu].dhze.l) +
499       Sqrt[gw^2+g1^2] Dual[FS][Z,mu,nu] (uqbar.Sig[mu,nu].tdhzu.uq + dqbar.Sig[mu,nu].tdhzd.dq + lbar.Sig[mu,nu].tdhze.l);
500  LW = Sqrt[2] gw (FS[Wbar,mu,nu] dqbar.Sig[mu,nu].dhwu.ProjP.uq +  FS[W,mu,nu] uqbar.Sig[mu,nu].dhwd.ProjP.dq + FS[W,mu,nu] vlbar.Sig[mu,nu].dhwl.ProjP.l);
501  -H/(4 vev^2) (LG+LA+LZ+LW + HC[LW])
502];
503
504(* eq.(3.12)  for the Higgs self interactions *)
505Lhself := -dL3 vev H^3 - dL4 H^4;
506
507(* eq (3.13) for the diHiggs couplings *)
508Lhh := Block[{Lvv,Lff},
509  Lvv = vev^2 dCw2 gw^2 Wbar[mu] W[mu] + dCz2 vev^2 (g1^2+gw^2)/2 Z[mu] Z[mu] +
510         gs^2/8        Cgg2 FS[G,mu,nu,a]  FS[G,mu,nu,a] + gs^2/8          tCgg2 FS[G,mu,nu,a]  Dual[FS][G,mu,nu,a] +
511         gw^2/4        Cww2 FS[Wbar,mu,nu] FS[W,mu,nu]   + gw^2/4          tCww2 FS[Wbar,mu,nu] Dual[FS][W,mu,nu]   +
512         ee^2/8        Caa2 FS[A,mu,nu]    FS[A,mu,nu]   + ee^2/8          tCaa2 FS[A,mu,nu]    Dual[FS][A,mu,nu]   +
513         gw g1/4       Cza2 FS[Z,mu,nu]    FS[A,mu,nu]   + gw g1/4         tCza2 FS[Z,mu,nu]    Dual[FS][A,mu,nu]   +
514         (gw^2+g1^2)/8 Czz2 FS[Z,mu,nu]    FS[Z,mu,nu]   + (gw^2 + g1^2)/8 tCzz2 FS[Z,mu,nu]    Dual[FS][Z,mu,nu]   -
515         gw^2 Cwbx2/2 (Wbar[mu] del[FS[W,mu,nu],nu] + W[mu] del[FS[Wbar,mu,nu],nu]) -
516         gw^2 Cabx2/2 Z[mu] del[FS[A,mu,nu],nu] -
517         gw^2 Czbx2/2 Z[mu] del[FS[Z,mu,nu],nu];
518  Lff = -ExpandIndices[Plus@@Flatten[Table[Sqrt[Mu[ff] Mu[ffp]]/2 dYu2[ff,ffp] ProjM[sp1,sp2] uqbar[sp1,ff,cc].uq[sp2,ffp,cc], {ffp,1,3},{ff,1,3}]]] -
519         ExpandIndices[Plus@@Flatten[Table[Sqrt[Md[ff] Md[ffp]]/2 dYd2[ff,ffp] ProjM[sp1,sp2] dqbar[sp1,ff,cc].dq[sp2,ffp,cc], {ffp,1,3},{ff,1,3}]]] -
520         ExpandIndices[Plus@@Flatten[Table[Sqrt[Ml[ff] Ml[ffp]]/2 dYe2[ff,ffp] ProjM[sp1,sp2] lbar[sp1,ff].l[sp2,ffp], {ffp,1,3},{ff,1,3}]]];
521  H^2/vev^2 ( Lvv + Lff + HC[Lff] )
522];
523
524(* LOther not included *)
525
526(* The sum of all contributions *)
527LBSMC:= {LMass, Lvertex, Ldip, Ltgc, Lqgc, Lh, Lhvff, Lhdvff, Lhself, Lhh};