← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to models/loop_qcd_qed_sm_Gmu__CMS__/lorentz.py

1. Merged with branch 2.3.1_CMS

Show diffs side-by-side

added added

removed removed

Lines of Context:
models/loop_qcd_qed_sm_Gmu__CMS__/lorentz.py
1
 
# This file was automatically created by FeynRules $Revision: 999 $
2
 
# Mathematica version: 7.0 for Linux x86 (64-bit) (February 18, 2009)
3
 
# Date: Mon 30 Jan 2012 19:57:04
4
 
 
5
 
 
6
 
from object_library import all_lorentz, Lorentz
7
 
 
8
 
from function_library import complexconjugate, re, im, csc, sec, acsc, asec
9
 
 
10
 
###################################
11
 
# CounterTerms Lorentz structures #
12
 
###################################
13
 
 
14
 
R2_GG_1 = Lorentz(name = 'R2_GG_1',
15
 
               spins = [ 3, 3 ],
16
 
               structure = 'P(-1,1)*P(-1,1)*Metric(1,2)')
17
 
 
18
 
R2_GG_2 = Lorentz(name = 'R2_GG_2',
19
 
               spins = [ 3, 3 ],
20
 
               structure = 'P(1,1)*P(2,1)')
21
 
 
22
 
R2_GG_3 = Lorentz(name = 'R2_GG_3',
23
 
               spins = [ 3, 3 ],
24
 
               structure = 'Metric(1,2)')
25
 
 
26
 
R2_QQ_1 = Lorentz(name = 'R2_QQ_1',
27
 
               spins = [ 2, 2 ],
28
 
               structure = 'P(-1,1)*Gamma(-1,2,1)')
29
 
 
30
 
R2_QQ_2 = Lorentz(name = 'R2_QQ_2',
31
 
               spins = [ 2, 2 ],
32
 
               structure = 'Identity(1,2)')
33
 
 
34
 
R2_QQ_3 = Lorentz(name = 'R2_QQ_3',
35
 
               spins = [ 2, 2 ],
36
 
               structure = 'P(-1,1)*Gamma(-1,2,-2)*ProjP(-2,1)')
37
 
 
38
 
R2_QQ_4 = Lorentz(name = 'R2_QQ_4',
39
 
                spins = [ 2, 2 ],
40
 
                structure = 'P(-1,1)*Gamma(-1,2,-2)*ProjM(-2,1)')
41
 
 
42
 
R2_SS_1 = Lorentz(name = 'R2_SS_1',
43
 
                  spins = [ 1, 1 ],
44
 
                  structure = '1')
45
 
 
46
 
R2_SS_2 = Lorentz(name = 'R2_SS_2',
47
 
                  spins = [ 1, 1 ],
48
 
                  structure = 'P(-1,1)*P(-1,1)')
49
 
 
50
 
GHGHG = Lorentz(name = 'GHGHG',
51
 
                 spins = [ 1, 1, 3 ],
52
 
                structure = 'P(3,1)')
53
 
 
54
 
#=============================================================================================
55
 
#  4-gluon R2 vertex
56
 
#=============================================================================================
57
 
 
58
 
 
59
 
R2_4G_1234 = Lorentz(name = 'R2_4G_1234',
60
 
                spins = [ 3, 3, 3, 3 ],
61
 
                structure = 'Metric(1,2)*Metric(3,4)')
62
 
 
63
 
R2_4G_1324 = Lorentz(name = 'R2_4G_1324',
64
 
                spins = [ 3, 3, 3, 3 ],
65
 
                structure = 'Metric(1,3)*Metric(2,4)')
66
 
 
67
 
R2_4G_1423 = Lorentz(name = 'R2_4G_1423',
68
 
                spins = [ 3, 3, 3, 3 ],
69
 
                structure = 'Metric(1,4)*Metric(2,3)')
70
 
 
71
 
# From FeynRules
72
 
 
73
 
R2RGA_VVVV10 = Lorentz(name = 'R2RGA_VVVV10',
74
 
                       spins = [ 3, 3, 3, 3 ],
75
 
                       structure = 'Metric(1,4)*Metric(2,3) + Metric(1,3)*Metric(2,4)\
76
 
                       + Metric(1,2)*Metric(3,4)')
77
 
 
78
 
R2RGA_VVVV2 = Lorentz(name = 'R2RGA_VVVV2',
79
 
                      spins = [ 3, 3, 3, 3 ],
80
 
                      structure = 'Metric(1,4)*Metric(2,3)')
81
 
 
82
 
R2RGA_VVVV3 = Lorentz(name = 'R2RGA_VVVV3',
83
 
                      spins = [ 3, 3, 3, 3 ],
84
 
                      structure = 'Metric(1,3)*Metric(2,4)')
85
 
 
86
 
R2RGA_VVVV5 = Lorentz(name = 'R2RGA_VVVV5',
87
 
                      spins = [ 3, 3, 3, 3 ],
88
 
                      structure = 'Metric(1,2)*Metric(3,4)')
89
 
 
90
 
#=============================================================================================
91
 
 
92
 
R2_GGZ = Lorentz(name = 'R2_GGZ',
93
 
                 spins = [ 3, 3, 3 ],
94
 
                 structure = 'Epsilon(3,1,2,-1)*P(-1,2)-Epsilon(3,1,2,-1)*P(-1,1)') 
95
 
 
96
 
R2_GGVV = Lorentz(name = 'R2_GGVV',
97
 
                 spins = [ 3, 3, 3, 3 ],
98
 
                 structure = 'Metric(1,2)*Metric(3,4)+Metric(1,3)*Metric(2,4)+Metric(1,4)*Metric(2,3)')
99
 
 
100
 
R2_GGHH = Lorentz(name = 'R2_GGHH',
101
 
                 spins = [ 3, 3, 1, 1 ],
102
 
                 structure = 'Metric(1,2)')
103
 
 
104
 
R2_GGGVa = Lorentz(name = 'R2_GGGVa',
105
 
                 spins = [ 3, 3, 3, 3 ],
106
 
                 structure = 'Epsilon(4,1,2,3)')
107
 
 
108
 
R2_VVVV1 = Lorentz(name = 'R2_VVVV1',
109
 
                spins = [ 3, 3, 3, 3 ],
110
 
                structure = 'Metric(1,2)*Metric(3,4)+Metric(1,3)*Metric(2,4)+Metric(1,4)*Metric(2,3)')
111
 
 
112
 
R2_VVVV2 = Lorentz(name = 'R2_VVVV2',
113
 
                   spins = [ 3, 3, 3, 3 ],
114
 
                   structure = 'Metric(1,2)*Metric(3,4)')
115
 
 
116
 
R2_VVVV3 = Lorentz(name = 'R2_VVVV3',
117
 
                   spins = [ 3, 3, 3, 3 ],
118
 
                   structure = 'Metric(1,3)*Metric(2,4)+Metric(1,4)*Metric(2,3)')
119
 
 
120
 
###################
121
 
# Base structures #
122
 
###################
123
 
 
124
 
 
125
 
UUS1 = Lorentz(name = 'UUS1',
126
 
               spins = [ 1, 1, 1 ],
127
 
               structure = '1')
128
 
 
129
 
UUV1 = Lorentz(name = 'UUV1',
130
 
               spins = [ 1, 1, 3 ],
131
 
               structure = 'P(3,2) + P(3,3)')
132
 
 
133
 
SSS1 = Lorentz(name = 'SSS1',
134
 
               spins = [ 1, 1, 1 ],
135
 
               structure = '1')
136
 
 
137
 
FFS1 = Lorentz(name = 'FFS1',
138
 
               spins = [ 2, 2, 1 ],
139
 
               structure = 'ProjM(2,1)')
140
 
 
141
 
FFS2 = Lorentz(name = 'FFS2',
142
 
               spins = [ 2, 2, 1 ],
143
 
               structure = 'ProjM(2,1) - ProjP(2,1)')
144
 
 
145
 
FFS3 = Lorentz(name = 'FFS3',
146
 
               spins = [ 2, 2, 1 ],
147
 
               structure = 'ProjP(2,1)')
148
 
 
149
 
FFS4 = Lorentz(name = 'FFS4',
150
 
               spins = [ 2, 2, 1 ],
151
 
               structure = 'ProjM(2,1) + ProjP(2,1)')
152
 
 
153
 
FFS5 = Lorentz(name = 'FFS5',
154
 
               spins = [ 2, 2, 1 ],
155
 
               structure = 'Identity(2,1)')
156
 
 
157
 
FFV1 = Lorentz(name = 'FFV1',
158
 
               spins = [ 2, 2, 3 ],
159
 
               structure = 'Gamma(3,2,1)')
160
 
 
161
 
FFV2 = Lorentz(name = 'FFV2',
162
 
               spins = [ 2, 2, 3 ],
163
 
               structure = 'Gamma(3,2,-1)*ProjM(-1,1)')
164
 
 
165
 
FFV3 = Lorentz(name = 'FFV3',
166
 
               spins = [ 2, 2, 3 ],
167
 
               structure = 'Gamma(3,2,-1)*ProjM(-1,1) - 2*Gamma(3,2,-1)*ProjP(-1,1)')
168
 
 
169
 
FFV4 = Lorentz(name = 'FFV4',
170
 
               spins = [ 2, 2, 3 ],
171
 
               structure = 'Gamma(3,2,-1)*ProjM(-1,1) + 2*Gamma(3,2,-1)*ProjP(-1,1)')
172
 
 
173
 
FFV5 = Lorentz(name = 'FFV5',
174
 
               spins = [ 2, 2, 3 ],
175
 
               structure = 'Gamma(3,2,-1)*ProjM(-1,1) + 4*Gamma(3,2,-1)*ProjP(-1,1)')
176
 
 
177
 
FFV6 = Lorentz(name = 'FFV6',
178
 
               spins = [ 2, 2, 3 ],
179
 
               structure = 'Gamma(3,2,-1)*ProjP(-1,1)')
180
 
 
181
 
VSS1 = Lorentz(name = 'VSS1',
182
 
               spins = [ 3, 1, 1 ],
183
 
               structure = 'P(1,2) - P(1,3)')
184
 
 
185
 
VVS1 = Lorentz(name = 'VVS1',
186
 
               spins = [ 3, 3, 1 ],
187
 
               structure = 'Metric(1,2)')
188
 
 
189
 
VVV1 = Lorentz(name = 'VVV1',
190
 
               spins = [ 3, 3, 3 ],
191
 
               structure = 'P(3,1)*Metric(1,2) - P(3,2)*Metric(1,2) - P(2,1)*Metric(1,3) + P(2,3)*Metric(1,3) + P(1,2)*Metric(2,3) - P(1,3)*Metric(2,3)')
192
 
 
193
 
SSSS1 = Lorentz(name = 'SSSS1',
194
 
                spins = [ 1, 1, 1, 1 ],
195
 
                structure = '1')
196
 
 
197
 
VVSS1 = Lorentz(name = 'VVSS1',
198
 
                spins = [ 3, 3, 1, 1 ],
199
 
                structure = 'Metric(1,2)')
200
 
 
201
 
VVVV1 = Lorentz(name = 'VVVV1',
202
 
                spins = [ 3, 3, 3, 3 ],
203
 
                structure = 'Metric(1,4)*Metric(2,3) - Metric(1,3)*Metric(2,4)')
204
 
 
205
 
VVVV2 = Lorentz(name = 'VVVV2',
206
 
                spins = [ 3, 3, 3, 3 ],
207
 
                structure = 'Metric(1,4)*Metric(2,3) + Metric(1,3)*Metric(2,4) - 2*Metric(1,2)*Metric(3,4)')
208
 
 
209
 
VVVV3 = Lorentz(name = 'VVVV3',
210
 
                spins = [ 3, 3, 3, 3 ],
211
 
                structure = 'Metric(1,4)*Metric(2,3) - Metric(1,2)*Metric(3,4)')
212
 
 
213
 
VVVV4 = Lorentz(name = 'VVVV4',
214
 
                spins = [ 3, 3, 3, 3 ],
215
 
                structure = 'Metric(1,3)*Metric(2,4) - Metric(1,2)*Metric(3,4)')
216
 
 
217
 
VVVV5 = Lorentz(name = 'VVVV5',
218
 
                spins = [ 3, 3, 3, 3 ],
219
 
                structure = 'Metric(1,4)*Metric(2,3) - (Metric(1,3)*Metric(2,4))/2. - (Metric(1,2)*Metric(3,4))/2.')
220
 
 
221
 
 
222
 
##############################################
223
 
# UV CounterTerms Lorentz structures for QED #
224
 
# Generate by WriteUFO automatically         # 
225
 
##############################################
226
 
 
227
 
l_WmWpMass1 = Lorentz(name = 'l_WmWpMass1',
228
 
                      spins = [ 3, 3 ],
229
 
                      structure = '-Metric(1,2)*P(-1,1)*P(-1,1)')
230
 
 
231
 
 
232
 
l_WmWpMass2 = Lorentz(name = 'l_WmWpMass2',
233
 
                      spins = [ 3, 3 ],
234
 
                      structure = 'Metric(1,2)')
235
 
 
236
 
 
237
 
l_WmWpMass3 = Lorentz(name = 'l_WmWpMass3',
238
 
                      spins = [ 3, 3 ],
239
 
                      structure = '-P(1,1)*P(2,1)')
240
 
 
241
 
 
242
 
l_GpWmMass4 = Lorentz(name = 'l_GpWmMass4',
243
 
                      spins = [ 1, 3 ],
244
 
                      structure = 'P(2,1)')
245
 
 
246
 
 
247
 
l_GpWmMass5 = Lorentz(name = 'l_GpWmMass5',
248
 
                      spins = [ 1, 3 ],
249
 
                      structure = 'P(2,2)')
250
 
 
251
 
 
252
 
l_HMass6 = Lorentz(name = 'l_HMass6',
253
 
                   spins = [ 1, 1 ],
254
 
                   structure = '-P(-1,1)*P(-1,1)')
255
 
 
256
 
 
257
 
l_HMass7 = Lorentz(name = 'l_HMass7',
258
 
                   spins = [ 1, 1 ],
259
 
                   structure = '1')
260
 
 
261
 
 
262
 
l_vevexMass8 = Lorentz(name = 'l_vevexMass8',
263
 
                       spins = [ 2, 2 ],
264
 
                       structure = 'P(-1,1)*Gamma(-1,2,-2)*ProjM(-2,1)')
265
 
 
266
 
 
267
 
l_vevexMass9 = Lorentz(name = 'l_vevexMass9',
268
 
                       spins = [ 2, 2 ],
269
 
                       structure = 'P(-1,2)*Gamma(-1,2,-2)*ProjP(-2,1)')
270
 
 
271
 
 
272
 
l_vevexMass10 = Lorentz(name = 'l_vevexMass10',
273
 
                        spins = [ 2, 2 ],
274
 
                        structure = 'ProjM(2,1)')
275
 
 
276
 
 
277
 
l_vevexMass11 = Lorentz(name = 'l_vevexMass11',
278
 
                        spins = [ 2, 2 ],
279
 
                        structure = 'ProjP(2,1)')
280
 
 
281
 
 
282
 
l_WpWpWmWm12 = Lorentz(name = 'l_WpWpWmWm12',
283
 
                       spins = [ 3, 3, 3, 3 ],
284
 
                       structure = 'Metric(1,2)*Metric(3,4)')
285
 
 
286
 
 
287
 
l_WpWpWmWm13 = Lorentz(name = 'l_WpWpWmWm13',
288
 
                       spins = [ 3, 3, 3, 3 ],
289
 
                       structure = 'Metric(1,3)*Metric(2,4)')
290
 
 
291
 
 
292
 
l_WpWpWmWm14 = Lorentz(name = 'l_WpWpWmWm14',
293
 
                       spins = [ 3, 3, 3, 3 ],
294
 
                       structure = 'Metric(1,4)*Metric(2,3)')
295
 
 
296
 
 
297
 
l_AWpWm15 = Lorentz(name = 'l_AWpWm15',
298
 
                    spins = [ 3, 3, 3 ],
299
 
                    structure = 'Metric(1,2)*(P(3,2)-P(3,1))+Metric(2,3)*(P(1,3)-P(1,2))+Metric(3,1)*(P(2,1)-P(2,3))')
300
 
 
301
 
 
302
 
l_HHHH16 = Lorentz(name = 'l_HHHH16',
303
 
                   spins = [ 1, 1, 1, 1 ],
304
 
                   structure = '1')
305
 
 
306
 
 
307
 
l_HHH17 = Lorentz(name = 'l_HHH17',
308
 
                  spins = [ 1, 1, 1 ],
309
 
                  structure = '1')
310
 
 
311
 
 
312
 
l_HHWmWp18 = Lorentz(name = 'l_HHWmWp18',
313
 
                     spins = [ 1, 1, 3, 3 ],
314
 
                     structure = 'Metric(3,4)')
315
 
 
316
 
 
317
 
l_G0HA19 = Lorentz(name = 'l_G0HA19',
318
 
                   spins = [ 1, 1, 3 ],
319
 
                   structure = 'P(3,1)-P(3,2)')
320
 
 
321
 
 
322
 
l_HWpWm20 = Lorentz(name = 'l_HWpWm20',
323
 
                    spins = [ 1, 3, 3 ],
324
 
                    structure = 'Metric(2,3)')
325
 
 
326
 
 
327
 
l_vexveA21 = Lorentz(name = 'l_vexveA21',
328
 
                     spins = [ 2, 2, 3 ],
329
 
                     structure = 'Gamma(3,2,-1)*ProjM(-1,1)')
330
 
 
331
 
 
332
 
l_vexveA22 = Lorentz(name = 'l_vexveA22',
333
 
                     spins = [ 2, 2, 3 ],
334
 
                     structure = 'Gamma(3,2,-1)*ProjP(-1,1)')
335
 
 
336
 
 
337
 
l_epemH23 = Lorentz(name = 'l_epemH23',
338
 
                    spins = [ 2, 2, 1 ],
339
 
                    structure = 'ProjM(2,1)')
340
 
 
341
 
 
342
 
l_epemH24 = Lorentz(name = 'l_epemH24',
343
 
                    spins = [ 2, 2, 1 ],
344
 
                    structure = 'ProjP(2,1)')
345
 
 
346
 
 
347
 
l_umumxA25 = Lorentz(name = 'l_umumxA25',
348
 
                     spins = [ -1, -1, 3 ],
349
 
                     structure = 'P(3,1)')
350
 
 
351
 
 
352
 
l_umumxA26 = Lorentz(name = 'l_umumxA26',
353
 
                     spins = [ -1, -1, 3 ],
354
 
                     structure = 'P(3,2)')
355
 
 
356
 
 
357
 
l_HuZuZx27 = Lorentz(name = 'l_HuZuZx27',
358
 
                     spins = [ 1, -1, -1 ],
359
 
                     structure = '1')
360
 
 
361