« back to all changes in this revision
Viewing changes to tests/loop_smgrav/CT_couplings.py
- browse files at revision 282.1.20
- compare with another revision
- download diff
- download tarball
- view history from revision 282.1.20
- Committer: olivier-mattelaer
- Date: 2019-09-13 14:04:58 UTC
- mto: This revision was merged to the branch mainline in revision 283.
- Revision ID: olivier-mattelaer-20190913140458-5ilmavms0mv5dhu0
change some default of the run_card and fix some tests
- files removed:
- tests/loop_smgrav
- files modified:
- UpdateNotes.txt
added
removed
tests/loop_smgrav/CT_couplings.py
1
# This file was automatically created by FeynRules $Revision: 535 $
2
# Mathematica version: 7.0 for Mac OS X x86 (64-bit) (November 11, 2008)
3
# Date: Fri 18 Mar 2011 18:40:51
4
5
from object_library import all_couplings, Coupling
6
7
from function_library import complexconjugate, re, im, csc, sec, acsc, asec
8
9
################
10
# R2 couplings #
11
################
12
13
# ========= #
14
# Pure QCD #
15
# ========= #
16
17
R2_3Gq = Coupling(name = 'R2_3Gq',
18
value = '2.0*G**3/(48.0*cmath.pi**2)',
19
order = {'QCD':3})
20
21
R2_3Gg = Coupling(name = 'R2_3Gg',
22
value = 'Ncol*G**3/(48.0*cmath.pi**2)*(7.0/4.0+lhv)',
23
order = {'QCD':3})
24
25
#=============================================================================================
26
# 4-gluon R2 couplings
27
#=============================================================================================
28
29
# Gluon contribution to it
30
31
GC_4GR2_Gluon_delta5 = Coupling(name = 'GC_4GR2_Gluon_delta5',
32
value = '-4.0*complex(0,1)*RGR2*(2.0*lhv+5.0)',
33
order = {'QCD':4})
34
35
GC_4GR2_Gluon_delta7 = Coupling(name = 'GC_4GR2_Gluon_delta7',
36
value = '2.0*complex(0,1)*RGR2*(2.0*lhv+7.0)',
37
order = {'QCD':4})
38
39
GC_4GR2_2Struct = Coupling(name = 'GC_4GR2_2Struct',
40
value = '2.0*complex(0,1)*RGR2*Ncol*(lhv+3.0)',
41
order = {'QCD':4})
42
43
GC_4GR2_4Struct = Coupling(name = 'GC_4GR2_4Struct',
44
value = '-complex(0,1)*RGR2*Ncol*(4.0*lhv+11.0)',
45
order = {'QCD':4})
46
47
# Fermion contribution to it
48
49
GC_4GR2_Fermion_delta5 = Coupling(name = 'GC_4GR2_Fermion_delta5',
50
value = '(2.0/Ncol)*5.0*complex(0,1)*RGR2',
51
order = {'QCD':4})
52
53
GC_4GR2_Fermion_delta11 = Coupling(name = 'GC_4GR2_Fermion_delta11',
54
value = '-(2.0/Ncol)*11.0*complex(0,1)*RGR2',
55
order = {'QCD':4})
56
57
GC_4GR2_5Struct = Coupling(name = 'GC_4GR2_5Struct',
58
value = '5.0*complex(0,1)*RGR2',
59
order = {'QCD':4})
60
61
GC_4GR2_11Struct = Coupling(name = 'GC_4GR2_11Struct',
62
value = '-11.0*complex(0,1)*RGR2',
63
order = {'QCD':4})
64
65
# From the auto UFO from FR
66
67
R2GC_137_43 = Coupling(name = 'R2GC_137_43',
68
value = '(complex(0,1)*G**4)/(192.*cmath.pi**2)',
69
order = {'QCD':4})
70
71
R2GC_137_44 = Coupling(name = 'R2GC_137_44',
72
value = '-(complex(0,1)*G**4)/(64.*cmath.pi**2)',
73
order = {'QCD':4})
74
75
R2GC_138_45 = Coupling(name = 'R2GC_138_45',
76
value = '-(complex(0,1)*G**4)/(192.*cmath.pi**2)',
77
order = {'QCD':4})
78
79
R2GC_138_46 = Coupling(name = 'R2GC_138_46',
80
value = '(complex(0,1)*G**4)/(64.*cmath.pi**2)',
81
order = {'QCD':4})
82
83
R2GC_139_47 = Coupling(name = 'R2GC_139_47',
84
value = '-(complex(0,1)*G**4)/(48.*cmath.pi**2)',
85
order = {'QCD':4})
86
87
R2GC_140_48 = Coupling(name = 'R2GC_140_48',
88
value = '(complex(0,1)*G**4)/(288.*cmath.pi**2)',
89
order = {'QCD':4})
90
91
R2GC_140_49 = Coupling(name = 'R2GC_140_49',
92
value = '-(complex(0,1)*G**4)/(32.*cmath.pi**2)',
93
order = {'QCD':4})
94
95
R2GC_141_50 = Coupling(name = 'R2GC_141_50',
96
value = '-(complex(0,1)*G**4)/(16.*cmath.pi**2)',
97
order = {'QCD':4})
98
99
R2GC_141_51 = Coupling(name = 'R2GC_141_51',
100
value = '(-7*complex(0,1)*G**4)/(32.*cmath.pi**2)',
101
order = {'QCD':4})
102
103
R2GC_142_52 = Coupling(name = 'R2GC_142_52',
104
value = '(complex(0,1)*G**4)/(16.*cmath.pi**2)',
105
order = {'QCD':4})
106
107
R2GC_142_53 = Coupling(name = 'R2GC_142_53',
108
value = '(7*complex(0,1)*G**4)/(32.*cmath.pi**2)',
109
order = {'QCD':4})
110
111
R2GC_143_54 = Coupling(name = 'R2GC_143_54',
112
value = '(11*complex(0,1)*G**4)/(192.*cmath.pi**2)',
113
order = {'QCD':4})
114
115
R2GC_143_55 = Coupling(name = 'R2GC_143_55',
116
value = '(15*complex(0,1)*G**4)/(64.*cmath.pi**2)',
117
order = {'QCD':4})
118
119
R2GC_144_56 = Coupling(name = 'R2GC_144_56',
120
value = '(-3*complex(0,1)*G**4)/(64.*cmath.pi**2)',
121
order = {'QCD':4})
122
123
R2GC_144_57 = Coupling(name = 'R2GC_144_57',
124
value = '(-17*complex(0,1)*G**4)/(64.*cmath.pi**2)',
125
order = {'QCD':4})
126
127
R2GC_145_58 = Coupling(name = 'R2GC_145_58',
128
value = '-G**4/(192.*cmath.pi**2)',
129
order = {'QCD':4})
130
131
R2GC_145_59 = Coupling(name = 'R2GC_145_59',
132
value = 'G**4/(64.*cmath.pi**2)',
133
order = {'QCD':4})
134
135
#=============================================================================================
136
137
R2_GQQ = Coupling(name = 'R2_GQQ',
138
value = '-complex(0,1)*G**3/(16.0*cmath.pi**2)*((Ncol**2-1)/(2.0*Ncol))*(1.0+lhv)',
139
order = {'QCD':3})
140
141
R2_GGq = Coupling(name = 'R2_GGq',
142
value = '(2.0)*complex(0,1)*G**2/(48.0*cmath.pi**2)',
143
order = {'QCD':2})
144
145
R2_GGb = Coupling(name = 'R2_GGb',
146
value = '(2.0)*complex(0,1)*G**2*(-6.0*MB**2)/(48.0*cmath.pi**2)',
147
order = {'QCD':2})
148
149
R2_GGc = Coupling(name = 'R2_GGc',
150
value = '(2.0)*complex(0,1)*G**2*(-6.0*MC**2)/(48.0*cmath.pi**2)',
151
order = {'QCD':2})
152
153
R2_GGt = Coupling(name = 'R2_GGt',
154
value = '(2.0)*complex(0,1)*G**2*(-6.0*MT**2)/(48.0*cmath.pi**2)',
155
order = {'QCD':2})
156
157
R2_GGg_1 = Coupling(name = 'R2_GGg_1',
158
value = '(2.0)*complex(0,1)*G**2*Ncol/(48.0*cmath.pi**2)*(1.0/2.0+lhv)',
159
order = {'QCD':2})
160
161
R2_GGg_2 = Coupling(name = 'R2_GGg_2',
162
value = '-(2.0)*complex(0,1)*G**2*Ncol/(48.0*cmath.pi**2)*lhv',
163
order = {'QCD':2})
164
165
R2_QQq = Coupling(name = 'R2_QQq',
166
value = 'lhv*complex(0,1)*G**2*(Ncol**2-1)/(32.0*cmath.pi**2*Ncol)',
167
order = {'QCD':2})
168
169
R2_QQc = Coupling(name = 'R2_QQc',
170
value = 'lhv*complex(0,1)*G**2*(Ncol**2-1)*(2.0*MC)/(32.0*cmath.pi**2*Ncol)',
171
order = {'QCD':2})
172
173
R2_QQb = Coupling(name = 'R2_QQb',
174
value = 'lhv*complex(0,1)*G**2*(Ncol**2-1)*(2.0*MB)/(32.0*cmath.pi**2*Ncol)',
175
order = {'QCD':2})
176
177
R2_QQt = Coupling(name = 'R2_QQt',
178
value = 'lhv*complex(0,1)*G**2*(Ncol**2-1)*(2.0*MT)/(32.0*cmath.pi**2*Ncol)',
179
order = {'QCD':2})
180
181
# ============== #
182
# Mixed QCD-QED #
183
# ============== #
184
185
# Quark couplings to A and Z
186
187
R2_DDA = Coupling(name = 'R2_DDA',
188
value = '(-(ee*complex(0,1))/3.0)*R2MixedFactor',
189
order = {'QCD':2,'QED':1})
190
191
R2_UUA = Coupling(name = 'R2_UUA',
192
value = '(2.0*(ee*complex(0,1))/3.0)*R2MixedFactor',
193
order = {'QCD':2,'QED':1})
194
195
R2_DDZ_V2 = Coupling(name = 'R2_DDZ_V2',
196
value = '(-(cw*ee*complex(0,1))/(2.0*sw))*R2MixedFactor',
197
order = {'QCD':2,'QED':1})
198
199
R2_DDZ_V3 = Coupling(name = 'R2_DDZ_V3',
200
value = '(-(ee*complex(0,1)*sw)/(6.0*cw))*R2MixedFactor',
201
order = {'QCD':2,'QED':1})
202
203
R2_UUZ_V2 = Coupling(name = 'R2_UUZ_V2',
204
value = '((cw*ee*complex(0,1))/(2.*sw))*R2MixedFactor',
205
order = {'QCD':2,'QED':1})
206
207
R2_UUZ_V5 = Coupling(name = 'R2_UUZ_V5',
208
value = '(-(ee*complex(0,1)*sw)/(6.*cw))*R2MixedFactor',
209
order = {'QCD':2,'QED':1})
210
211
# Quark couplings to W, most general CKM
212
213
R2_dxuW = Coupling(name = 'R2_dxuW',
214
value = '((CKM11*ee*complex(0,1))/(sw*cmath.sqrt(2)))*R2MixedFactor',
215
order = {'QCD':2,'QED':1})
216
217
R2_dxcW = Coupling(name = 'R2_dxcW',
218
value = '((CKM21*ee*complex(0,1))/(sw*cmath.sqrt(2)))*R2MixedFactor',
219
order = {'QCD':2,'QED':1})
220
221
R2_dxtW = Coupling(name = 'R2_dxtW',
222
value = '((CKM31*ee*complex(0,1))/(sw*cmath.sqrt(2)))*R2MixedFactor',
223
order = {'QCD':2,'QED':1})
224
225
R2_sxuW = Coupling(name = 'R2_sxuW',
226
value = '((CKM12*ee*complex(0,1))/(sw*cmath.sqrt(2)))*R2MixedFactor',
227
order = {'QCD':2,'QED':1})
228
229
R2_sxcW = Coupling(name = 'R2_sxcW',
230
value = '((CKM22*ee*complex(0,1))/(sw*cmath.sqrt(2)))*R2MixedFactor',
231
order = {'QCD':2,'QED':1})
232
233
R2_sxtW = Coupling(name = 'R2_sxtW',
234
value = '((CKM32*ee*complex(0,1))/(sw*cmath.sqrt(2)))*R2MixedFactor',
235
order = {'QCD':2,'QED':1})
236
237
R2_bxuW = Coupling(name = 'R2_bxuW',
238
value = '((CKM13*ee*complex(0,1))/(sw*cmath.sqrt(2)))*R2MixedFactor',
239
order = {'QCD':2,'QED':1})
240
241
R2_bxcW = Coupling(name = 'R2_bxcW',
242
value = '((CKM23*ee*complex(0,1))/(sw*cmath.sqrt(2)))*R2MixedFactor',
243
order = {'QCD':2,'QED':1})
244
245
R2_bxtW = Coupling(name = 'R2_bxtW',
246
value = '((CKM33*ee*complex(0,1))/(sw*cmath.sqrt(2)))*R2MixedFactor',
247
order = {'QCD':2,'QED':1})
248
249
R2_uxdW = Coupling(name = 'R2_uxdW',
250
value = '((ee*complex(0,1)*complexconjugate(CKM11))/(sw*cmath.sqrt(2)))*R2MixedFactor',
251
order = {'QCD':2,'QED':1})
252
253
R2_cxdW = Coupling(name = 'R2_cxdW',
254
value = '((ee*complex(0,1)*complexconjugate(CKM21))/(sw*cmath.sqrt(2)))*R2MixedFactor',
255
order = {'QCD':2,'QED':1})
256
257
R2_txdW = Coupling(name = 'R2_txdW',
258
value = '((ee*complex(0,1)*complexconjugate(CKM31))/(sw*cmath.sqrt(2)))*R2MixedFactor',
259
order = {'QCD':2,'QED':1})
260
261
R2_uxsW = Coupling(name = 'R2_uxsW',
262
value = '((ee*complex(0,1)*complexconjugate(CKM12))/(sw*cmath.sqrt(2)))*R2MixedFactor',
263
order = {'QCD':2,'QED':1})
264
265
R2_cxsW = Coupling(name = 'R2_cxsW',
266
value = '((ee*complex(0,1)*complexconjugate(CKM22))/(sw*cmath.sqrt(2)))*R2MixedFactor',
267
order = {'QCD':2,'QED':1})
268
269
R2_txsW = Coupling(name = 'R2_txsW',
270
value = '((ee*complex(0,1)*complexconjugate(CKM32))/(sw*cmath.sqrt(2)))*R2MixedFactor',
271
order = {'QCD':2,'QED':1})
272
273
R2_uxbW = Coupling(name = 'R2_uxbW',
274
value = '((ee*complex(0,1)*complexconjugate(CKM13))/(sw*cmath.sqrt(2)))*R2MixedFactor',
275
order = {'QCD':2,'QED':1})
276
277
R2_cxbW = Coupling(name = 'R2_cxbW',
278
value = '((ee*complex(0,1)*complexconjugate(CKM23))/(sw*cmath.sqrt(2)))*R2MixedFactor',
279
order = {'QCD':2,'QED':1})
280
281
R2_txbW = Coupling(name = 'R2_txbW',
282
value = '((ee*complex(0,1)*complexconjugate(CKM33))/(sw*cmath.sqrt(2)))*R2MixedFactor',
283
order = {'QCD':2,'QED':1})
284
285
# R2 for the Higgs couplings to massive quarks
286
287
R2_bbH = Coupling(name = 'R2_bbH',
288
value = '(-((complex(0,1)*yb)/cmath.sqrt(2)))*(2.0*R2MixedFactor)',
289
order = {'QCD':2,'QED':1})
290
291
R2_ttH = Coupling(name = 'R2_ttH',
292
value = '(-((complex(0,1)*yt)/cmath.sqrt(2)))*(2.0*R2MixedFactor)',
293
order = {'QCD':2,'QED':1})
294
295
R2_ccH = Coupling(name = 'R2_ccH',
296
value = '(-((complex(0,1)*yc)/cmath.sqrt (2)))*(2.0*R2MixedFactor)',
297
order = {'QCD':2,'QED':1})
298
299
# R2 interactions non proportional to the SM
300
301
# R2 for the Higgs interactions
302
R2_GGHc = Coupling(name = 'R2_GGHc',
303
value = '4.0*(-((complex(0,1)*yc)/cmath.sqrt(2)))*(1.0/2.0)*(G**2/(8.0*cmath.pi**2))*MC',
304
order = {'QCD':2,'QED':1})
305
306
R2_GGHb = Coupling(name = 'R2_GGHb',
307
value = '4.0*(-((complex(0,1)*yb)/cmath.sqrt(2)))*(1.0/2.0)*(G**2/(8.0*cmath.pi**2))*MB',
308
order = {'QCD':2,'QED':1})
309
310
R2_GGHt = Coupling(name = 'R2_GGHt',
311
value = '4.0*(-((complex(0,1)*yt)/cmath.sqrt(2)))*(1.0/2.0)*(G**2/(8.0*cmath.pi**2))*MT',
312
order = {'QCD':2,'QED':1})
313
314
# R2 for the weak vector bosons interaction with gluons
315
316
R2_GGZup = Coupling(name = 'R2_GGZup',
317
value = {0:'AxialZUp*(1.0/2.0)*(G**2/(12.0*cmath.pi**2))'},
318
order = {'QCD':2,'QED':1})
319
320
R2_GGZdown = Coupling(name = 'R2_GGZdown',
321
value = {0:'AxialZDown*(1.0/2.0)*(G**2/(12.0*cmath.pi**2))'},
322
order = {'QCD':2,'QED':1})
323
324
# EDIT VH
325
# There is a factor four added here to all the GGVV R2 couplings below without proper understanding of it
326
# I have not investigate too much about it though.
327
# END EDIT VH
328
329
R2_GGZAup = Coupling(name = 'R2_GGZAup',
330
value = {0:'4.0*(-VectorAUp*VectorZUp)*(1.0/2.0)*(-(complex(0,1)*G**2)/(24.0*cmath.pi**2))'},
331
order = {'QCD':2,'QED':2})
332
333
R2_GGZAdown = Coupling(name = 'R2_GGZAdown',
334
value = {0:'4.0*(-VectorADown*VectorZDown)*(1.0/2.0)*(-(complex(0,1)*G**2)/(24.0*cmath.pi**2))'},
335
order = {'QCD':2,'QED':2})
336
337
R2_GGZZdown = Coupling(name = 'R2_GGZZdown',
338
value = {0:'4.0*(-VectorZDown*VectorZDown-AxialZDown*AxialZDown)*(1.0/2.0)*(-(complex(0,1)*G**2)/(24.0*cmath.pi**2))'},
339
order = {'QCD':2,'QED':2})
340
341
R2_GGZZup = Coupling(name = 'R2_GGZZup',
342
value = {0:'4.0*(-VectorZUp*VectorZUp-AxialZUp*AxialZUp)*(1.0/2.0)*(-(complex(0,1)*G**2)/(24.0*cmath.pi**2))'},
343
order = {'QCD':2,'QED':2})
344
345
R2_GGAAdown = Coupling(name = 'R2_GGAAdown',
346
value = {0:'4.0*(-VectorADown*VectorADown)*(1.0/2.0)*(-(complex(0,1)*G**2)/(24.0*cmath.pi**2))'},
347
order = {'QCD':2,'QED':2})
348
349
R2_GGAAup = Coupling(name = 'R2_GGAAup',
350
value = {0:'4.0*(-VectorAUp*VectorAUp)*(1.0/2.0)*(-(complex(0,1)*G**2)/(24.0*cmath.pi**2))'},
351
order = {'QCD':2,'QED':2})
352
353
R2_GGWWud = Coupling(name = 'R2_GGWWud',
354
value = {0:'4.0*(CKM11*complexconjugate(CKM11))*(-VectorWmDxU*VectorWpUxD-AxialWmDxU*AxialWpUxD)*(1.0/2.0)*(-(complex(0,1)*G**2)/(24.0*cmath.pi**2))'},
355
order = {'QCD':2,'QED':2})
356
357
R2_GGWWus = Coupling(name = 'R2_GGWWus',
358
value = {0:'4.0*(CKM12*complexconjugate(CKM12))*(-VectorWmDxU*VectorWpUxD-AxialWmDxU*AxialWpUxD)*(1.0/2.0)*(-(complex(0,1)*G**2)/(24.0*cmath.pi**2))'},
359
order = {'QCD':2,'QED':2})
360
361
R2_GGWWub = Coupling(name = 'R2_GGWWub',
362
value = {0:'4.0*(CKM13*complexconjugate(CKM13))*(-VectorWmDxU*VectorWpUxD-AxialWmDxU*AxialWpUxD)*(1.0/2.0)*(-(complex(0,1)*G**2)/(24.0*cmath.pi**2))'},
363
order = {'QCD':2,'QED':2})
364
365
R2_GGWWcd = Coupling(name = 'R2_GGWWcd',
366
value = {0:'4.0*(CKM21*complexconjugate(CKM21))*(-VectorWmDxU*VectorWpUxD-AxialWmDxU*AxialWpUxD)*(1.0/2.0)*(-(complex(0,1)*G**2)/(24.0*cmath.pi**2))'},
367
order = {'QCD':2,'QED':2})
368
369
R2_GGWWcs = Coupling(name = 'R2_GGWWcs',
370
value = {0:'4.0*(CKM22*complexconjugate(CKM22))*(-VectorWmDxU*VectorWpUxD-AxialWmDxU*AxialWpUxD)*(1.0/2.0)*(-(complex(0,1)*G**2)/(24.0*cmath.pi**2))'},
371
order = {'QCD':2,'QED':2})
372
373
R2_GGWWcb = Coupling(name = 'R2_GGWWcb',
374
value = {0:'4.0*(CKM23*complexconjugate(CKM23))*(-VectorWmDxU*VectorWpUxD-AxialWmDxU*AxialWpUxD)*(1.0/2.0)*(-(complex(0,1)*G**2)/(24.0*cmath.pi**2))'},
375
order = {'QCD':2,'QED':2})
376
377
R2_GGWWtd = Coupling(name = 'R2_GGWWtd',
378
value = {0:'4.0*(CKM31*complexconjugate(CKM31))*(-VectorWmDxU*VectorWpUxD-AxialWmDxU*AxialWpUxD)*(1.0/2.0)*(-(complex(0,1)*G**2)/(24.0*cmath.pi**2))'},
379
order = {'QCD':2,'QED':2})
380
381
R2_GGWWts = Coupling(name = 'R2_GGWWts',
382
value = {0:'4.0*(CKM32*complexconjugate(CKM32))*(-VectorWmDxU*VectorWpUxD-AxialWmDxU*AxialWpUxD)*(1.0/2.0)*(-(complex(0,1)*G**2)/(24.0*cmath.pi**2))'},
383
order = {'QCD':2,'QED':2})
384
385
R2_GGWWtb = Coupling(name = 'R2_GGWWtb',
386
value = {0:'4.0*(CKM33*complexconjugate(CKM33))*(-VectorWmDxU*VectorWpUxD-AxialWmDxU*AxialWpUxD)*(1.0/2.0)*(-(complex(0,1)*G**2)/(24.0*cmath.pi**2))'},
387
order = {'QCD':2,'QED':2})
388
389
R2_GGHHc = Coupling(name = 'R2_GGHHc',
390
value = {0:'(-yc**2/2.0)*(1.0/2.0)*((complex(0,1)*G**2)/(8.0*cmath.pi**2))'},
391
order = {'QCD':2,'QED':2})
392
393
# EDIT VH
394
# There is a factor +4 added here to the two couplings below without proper understanding of it
395
# I have not investigate too much about it though.
396
# END EDIT VH
397
R2_GGHHb = Coupling(name = 'R2_GGHHb',
398
value = {0:'4.0*(-yb**2/2.0)*(1.0/2.0)*((complex(0,1)*G**2)/(8.0*cmath.pi**2))'},
399
order = {'QCD':2,'QED':2})
400
401
R2_GGHHt = Coupling(name = 'R2_GGHHt',
402
value = {0:'4.0*(-yt**2/2.0)*(1.0/2.0)*((complex(0,1)*G**2)/(8.0*cmath.pi**2))'},
403
order = {'QCD':2,'QED':2})
404
405
# EDIT VH
406
# There is a minus sign added here to the coupling below without proper understanding of it
407
# I have not investigate too much about it though.
408
# END EDIT VH
409
R2_GGGZvecUp = Coupling(name = 'R2_GGGZvecUp',
410
value = {0:'complex(0,1)*VectorZUp*(-1.0/2.0)*(-G**3/(24.0*cmath.pi**2))'},
411
order = {'QCD':3,'QED':1})
412
413
# EDIT VH
414
# There is a minus sign added here to the coupling below without proper understanding of it
415
# I have not investigate too much about it though.
416
# END EDIT VH
417
R2_GGGZvecDown = Coupling(name = 'R2_GGGZvecDown',
418
value = {0:'complex(0,1)*VectorZDown*(-1.0/2.0)*(-G**3/(24.0*cmath.pi**2))'},
419
order = {'QCD':3,'QED':1})
420
421
R2_GGGZaxialUp = Coupling(name = 'R2_GGGZaxialUp',
422
value = {0:'complex(0,1)*AxialZUp*(9.0/2.0)*(-G**3/(24.0*cmath.pi**2))'},
423
order = {'QCD':3,'QED':1})
424
425
R2_GGGZaxialDown = Coupling(name = 'R2_GGGZaxialDown',
426
value = {0:'complex(0,1)*AxialZDown*(9.0/2.0)*(-G**3/(24.0*cmath.pi**2))'},
427
order = {'QCD':3,'QED':1})
428
429
# EDIT VH
430
# There is a minus signe added here to the coupling below without proper understanding of it
431
# I have not investigate too much about it though.
432
# END EDIT VH
433
R2_GGGAvecUp = Coupling(name = 'R2_GGGAvecUp',
434
value = {0:'complex(0,1)*VectorAUp*(-1.0/2.0)*(-G**3/(24.0*cmath.pi**2))'},
435
order = {'QCD':3,'QED':1})
436
# EDIT VH
437
# There is a minus signe added here to the coupling below without proper understanding of it
438
# I have not investigate too much about it though.
439
# END EDIT VH
440
R2_GGGAvecDown = Coupling(name = 'R2_GGGAvecDown',
441
value = {0:'complex(0,1)*VectorADown*(-1.0/2.0)*(-G**3/(24.0*cmath.pi**2))'},
442
order = {'QCD':3,'QED':1})
443
444
################
445
# UV couplings #
446
################
447
448
# ========= #
449
# Pure QCD #
450
# ========= #
451
452
UV_3Gg = Coupling(name = 'UV_3Gg',
453
value = '-G_UVg*G',
454
order = {'QCD':3})
455
456
UV_3Gq = Coupling(name = 'UV_3Gq',
457
value = '-G_UVq*G',
458
order = {'QCD':3})
459
460
UV_3Gc = Coupling(name = 'UV_3Gc',
461
value = '-G_UVc*G',
462
order = {'QCD':3})
463
464
UV_3Gb = Coupling(name = 'UV_3Gb',
465
value = '-G_UVb*G',
466
order = {'QCD':3})
467
468
UV_3Gt = Coupling(name = 'UV_3Gt',
469
value = '-G_UVt*G',
470
order = {'QCD':3})
471
472
UV_4Gg = Coupling(name = 'UV_4Gg',
473
value = '2.0*complex(0,1)*G_UVg*(G**2)',
474
order = {'QCD':4})
475
476
UV_4Gq = Coupling(name = 'UV_4Gq',
477
value = '2.0*complex(0,1)*G_UVq*(G**2)',
478
order = {'QCD':4})
479
480
UV_4Gc = Coupling(name = 'UV_4Gc',
481
value = '2.0*complex(0,1)*G_UVc*(G**2)',
482
order = {'QCD':4})
483
484
UV_4Gb = Coupling(name = 'UV_4Gb',
485
value = '2.0*complex(0,1)*G_UVb*(G**2)',
486
order = {'QCD':4})
487
488
UV_4Gt = Coupling(name = 'UV_4Ggt',
489
value = '2.0*complex(0,1)*G_UVt*(G**2)',
490
order = {'QCD':4})
491
492
UV_GQQg = Coupling(name = 'UV_GQQg',
493
value = 'complex(0,1)*G_UVg*G',
494
order = {'QCD':3})
495
496
UV_GQQq = Coupling(name = 'UV_GQQq',
497
value = 'complex(0,1)*G_UVq*G',
498
order = {'QCD':3})
499
500
UV_GQQc = Coupling(name = 'UV_GQQc',
501
value = 'complex(0,1)*G_UVc*G',
502
order = {'QCD':3})
503
504
UV_GQQb = Coupling(name = 'UV_GQQb',
505
value = 'complex(0,1)*G_UVb*G',
506
order = {'QCD':3})
507
508
UV_GQQt = Coupling(name = 'UV_GQQt',
509
value = 'complex(0,1)*G_UVt*G',
510
order = {'QCD':3})
511
512
UV_cMass = Coupling(name = 'UV_cMass',
513
value = 'cMass_UV',
514
order = {'QCD':2})
515
516
UV_bMass = Coupling(name = 'UV_bMass',
517
value = 'bMass_UV',
518
order = {'QCD':2})
519
520
UV_tMass = Coupling(name = 'UV_tMass',
521
value = 'tMass_UV',
522
order = {'QCD':2})
523
524
# ============== #
525
# Mixed QCD-QED #
526
# ============== #
527
528
UV_Hcc = Coupling(name = 'UV_Hcc',
529
value = '-((complex(0,1)*yc)/cmath.sqrt(2))*UV_yuk_c',
530
order = {'QED':1,'QCD':2})
531
532
UV_Htt = Coupling(name = 'UV_Htt',
533
value = '-((complex(0,1)*yt)/cmath.sqrt(2))*UV_yuk_t',
534
order = {'QED':1,'QCD':2})
535
536
UV_Hbb = Coupling(name = 'UV_Hbb',
537
value = '-((complex(0,1)*yb)/cmath.sqrt(2))*UV_yuk_b',
538
order = {'QED':1,'QCD':2})
539
540
# ============================== #
541
# Goldstone R2 CT couplings #
542
# ============================== #
543
544
R2_GGGpGm_ub = Coupling(name = 'R2_GGGpGm_ub',
545
value = 'complex(0,1)*R2_GGGpGm_factor*(Vector_ubGp*Vector_ubGm-Axial_ubGp*Axial_ubGm)',
546
order = {'QED':2,'QCD':2})
547
548
R2_GGGpGm_cd = Coupling(name = 'R2_GGGpGm_cd',
549
value = 'complex(0,1)*R2_GGGpGm_factor*(Vector_cdGp*Vector_cdGm-Axial_cdGp*Axial_cdGm)',
550
order = {'QED':2,'QCD':2})
551
552
R2_GGGpGm_cs = Coupling(name = 'R2_GGGpGm_cs',
553
value = 'complex(0,1)*R2_GGGpGm_factor*(Vector_csGp*Vector_csGm-Axial_csGp*Axial_csGm)',
554
order = {'QED':2,'QCD':2})
555
556
R2_GGGpGm_cb = Coupling(name = 'R2_GGGpGm_cb',
557
value = 'complex(0,1)*R2_GGGpGm_factor*(Vector_cbGp*Vector_cbGm-Axial_cbGp*Axial_cbGm)',
558
order = {'QED':2,'QCD':2})
559
560
R2_GGGpGm_td = Coupling(name = 'R2_GGGpGm_td',
561
value = 'complex(0,1)*R2_GGGpGm_factor*(Vector_tdGp*Vector_tdGm-Axial_tdGp*Axial_tdGm)',
562
order = {'QED':2,'QCD':2})
563
564
R2_GGGpGm_ts = Coupling(name = 'R2_GGGpGm_ts',
565
value = 'complex(0,1)*R2_GGGpGm_factor*(Vector_tsGp*Vector_tsGm-Axial_tsGp*Axial_tsGm)',
566
order = {'QED':2,'QCD':2})
567
568
R2_GGGpGm_tb = Coupling(name = 'R2_GGGpGm_tb',
569
value = 'complex(0,1)*R2_GGGpGm_factor*(Vector_tbGp*Vector_tbGm-Axial_tbGp*Axial_tbGm)',
570
order = {'QED':2,'QCD':2})
571
572
R2_GGG0G0_c = Coupling(name = 'R2_GGG0G0_c',
573
value = 'complex(0,1)*R2_GGG0G0_factor*(-(1.0/2.0)*yc**2)',
574
order = {'QED':2,'QCD':2})
575
576
R2_GGG0G0_b = Coupling(name = 'R2_GGG0G0_b',
577
value = 'complex(0,1)*R2_GGG0G0_factor*(-(1.0/2.0)*yb**2)',
578
order = {'QED':2,'QCD':2})
579
580
R2_GGG0G0_t = Coupling(name = 'R2_GGG0G0_t',
581
value = 'complex(0,1)*R2_GGG0G0_factor*(-(1.0/2.0)*yt**2)',
582
order = {'QED':2,'QCD':2})
583
584
GC_R2_1015 = Coupling(name = 'GC_R2_1015',
585
value = 'I1x33*(R2MixedFactor*2.0)',
586
order = {'QED':1,'QCD':2})
587
588
GC_R2_1021 = Coupling(name = 'GC_R2_1021',
589
value = '-I2x33*(R2MixedFactor*2.0)',
590
order = {'QED':1,'QCD':2})
591
592
GC_R2_1027 = Coupling(name = 'GC_R2_1027',
593
value = 'I3x33*(R2MixedFactor*2.0)',
594
order = {'QED':1,'QCD':2})
595
596
GC_R2_1030 = Coupling(name = 'GC_R2_1030',
597
value = '-I4x33*(R2MixedFactor*2.0)',
598
order = {'QED':1,'QCD':2})
599
600
GC_R2_1082 = Coupling(name = 'GC_R2_1082',
601
value = '-(yb/cmath.sqrt(2))*(R2MixedFactor*2.0)',
602
order = {'QED':1,'QCD':2})
603
604
GC_R2_1085 = Coupling(name = 'GC_R2_1085',
605
value = '(yc/cmath.sqrt(2))*(R2MixedFactor*2.0)',
606
order = {'QED':1,'QCD':2})
607
608
GC_R2_1095 = Coupling(name = 'GC_R2_1095',
609
value = '(yt/cmath.sqrt(2))*(R2MixedFactor*2.0)',
610
order = {'QED':1,'QCD':2})
611
612
GC_R2_1018 = Coupling(name = 'GC_R2_1018',
613
value = '(-I2x22)*(R2MixedFactor*2.0)',
614
order = {'QED':1,'QCD':2})
615
616
GC_R2_1019 = Coupling(name = 'GC_R2_1019',
617
value = '(-I2x23)*(R2MixedFactor*2.0)',
618
order = {'QED':1,'QCD':2})
619
620
GC_R2_1013 = Coupling(name = 'GC_R2_1013',
621
value = '(I1x31)*(R2MixedFactor*2.0)',
622
order = {'QED':1,'QCD':2})
623
624
GC_R2_1014 = Coupling(name = 'GC_R2_1014',
625
value = '(I1x32)*(R2MixedFactor*2.0)',
626
order = {'QED':1,'QCD':2})
627
628
GC_R2_1020 = Coupling(name = 'GC_R2_1020',
629
value = '(-I2x32)*(R2MixedFactor*2.0)',
630
order = {'QED':1,'QCD':2})
631
632
GC_R2_1028 = Coupling(name = 'GC_R2_1028',
633
value = '(-I4x13)*(R2MixedFactor*2.0)',
634
order = {'QED':1,'QCD':2})
635
636
GC_R2_1022 = Coupling(name = 'GC_R2_1022',
637
value = '(I3x21)*(R2MixedFactor*2.0)',
638
order = {'QED':1,'QCD':2})
639
640
GC_R2_1023 = Coupling(name = 'GC_R2_1023',
641
value = '(I3x22)*(R2MixedFactor*2.0)',
642
order = {'QED':1,'QCD':2})
643
644
GC_R2_1024 = Coupling(name = 'GC_R2_1024',
645
value = '(I3x23)*(R2MixedFactor*2.0)',
646
order = {'QED':1,'QCD':2})
647
648
GC_R2_1029 = Coupling(name = 'GC_R2_1029',
649
value = '(-I4x23)*(R2MixedFactor*2.0)',
650
order = {'QED':1,'QCD':2})
651
652
GC_R2_1025 = Coupling(name = 'GC_R2_1025',
653
value = '(I3x31)*(R2MixedFactor*2.0)',
654
order = {'QED':1,'QCD':2})
655
656
GC_R2_1026 = Coupling(name = 'GC_R2_1026',
657
value = '(I3x32)*(R2MixedFactor*2.0)',
658
order = {'QED':1,'QCD':2})
659
660
# ============================== #
661
# Goldstone UV CT couplings #
662
# ============================== #
663
664
GC_UV_1015 = Coupling(name = 'GC_UV_1015',
665
value = 'I1x33*UV_yuk_b',
666
order = {'QED':1,'QCD':2})
667
668
GC_UV_1021 = Coupling(name = 'GC_UV_1021',
669
value = '-I2x33*UV_yuk_t',
670
order = {'QED':1,'QCD':2})
671
672
GC_UV_1027 = Coupling(name = 'GC_UV_1027',
673
value = 'I3x33*UV_yuk_t',
674
order = {'QED':1,'QCD':2})
675
676
GC_UV_1030 = Coupling(name = 'GC_UV_1030',
677
value = '-I4x33*UV_yuk_b',
678
order = {'QED':1,'QCD':2})
679
680
GC_UV_1082 = Coupling(name = 'GC_UV_1082',
681
value = '-(yb/cmath.sqrt(2))*(UV_yuk_b)',
682
order = {'QED':1,'QCD':2})
683
684
GC_UV_1085 = Coupling(name = 'GC_UV_1085',
685
value = '(yc/cmath.sqrt(2))*(UV_yuk_c)',
686
order = {'QED':1,'QCD':2})
687
688
GC_UV_1095 = Coupling(name = 'GC_UV_1095',
689
value = '(yt/cmath.sqrt(2))*(UV_yuk_t)',
690
order = {'QED':1,'QCD':2})
691
692
GC_UV_1018 = Coupling(name = 'GC_UV_1018',
693
value = '(-I2x22)*(UV_yuk_c)',
694
order = {'QED':1,'QCD':2})
695
696
GC_UV_1019 = Coupling(name = 'GC_UV_1019',
697
value = '(-I2x23)*(UV_yuk_t)',
698
order = {'QED':1,'QCD':2})
699
700
GC_UV_1013 = Coupling(name = 'GC_UV_1013',
701
value = '(I1x31)*(UV_yuk_b)',
702
order = {'QED':1,'QCD':2})
703
704
GC_UV_1014 = Coupling(name = 'GC_UV_1014',
705
value = '(I1x32)*(UV_yuk_b)',
706
order = {'QED':1,'QCD':2})
707
708
GC_UV_1020 = Coupling(name = 'GC_UV_1020',
709
value = '(-I2x32)*(UV_yuk_c)',
710
order = {'QED':1,'QCD':2})
711
712
GC_UV_1028 = Coupling(name = 'GC_UV_1028',
713
value = '(-I4x13)*(UV_yuk_b)',
714
order = {'QED':1,'QCD':2})
715
716
GC_UV_1022 = Coupling(name = 'GC_UV_1022',
717
value = '(I3x21)*(UV_yuk_c)',
718
order = {'QED':1,'QCD':2})
719
720
GC_UV_1023 = Coupling(name = 'GC_UV_1023',
721
value = '(I3x22)*(UV_yuk_c)',
722
order = {'QED':1,'QCD':2})
723
724
GC_UV_1024 = Coupling(name = 'GC_UV_1024',
725
value = '(I3x23)*(UV_yuk_c)',
726
order = {'QED':1,'QCD':2})
727
728
GC_UV_1029 = Coupling(name = 'GC_UV_1029',
729
value = '(-I4x23)*(UV_yuk_b)',
730
order = {'QED':1,'QCD':2})
731
732
GC_UV_1025 = Coupling(name = 'GC_UV_1025',
733
value = '(I3x31)*(UV_yuk_t)',
734
order = {'QED':1,'QCD':2})
735
736
GC_UV_1026 = Coupling(name = 'GC_UV_1026',
737
value = '(I3x32)*(UV_yuk_t)',
738
order = {'QED':1,'QCD':2})
Loggerhead is a web-based interface for Breezy
Version: 2.0.1