← Back to branch summary

~maddevelopers/mg5amcnlo/2.6.3_rwgt

« back to all changes in this revision

Viewing changes to tests/input_files/loop_smgrav/particles.py

  • Committer: olivier-mattelaer
  • Date: 2017-08-16 21:28:30 UTC
  • mfrom: (272.1.60 2.5.6)
  • Revision ID: olivier-mattelaer-20170816212830-sfl8jfv8pmpm5l2j
pass to 2.6.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/loop_smgrav/particles.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
 
 
6
from __future__ import division
 
7
from object_library import all_particles, Particle
 
8
import parameters as Param
 
9
import CT_parameters as CTParam
 
10
 
 
11
# ======================================================================
 
12
# QCD particles
 
13
# ======================================================================
 
14
 
 
15
d = Particle(pdg_code = 1,
 
16
             name = 'd',
 
17
             antiname = 'd~',
 
18
             spin = 2,
 
19
             color = 3,
 
20
             mass = Param.ZERO,
 
21
             width = Param.ZERO,
 
22
             texname = 'd',
 
23
             antitexname = 'd',
 
24
             charge = -1/3,
 
25
             LeptonNumber = 0,
 
26
             GhostNumber = 0)
 
27
 
 
28
d__tilde__ = d.anti()
 
29
 
 
30
u = Particle(pdg_code = 2,
 
31
             name = 'u',
 
32
             antiname = 'u~',
 
33
             spin = 2,
 
34
             color = 3,
 
35
             mass = Param.ZERO,
 
36
             width = Param.ZERO,
 
37
             texname = 'u',
 
38
             antitexname = 'u',
 
39
             charge = 2/3,
 
40
             LeptonNumber = 0,
 
41
             GhostNumber = 0)
 
42
 
 
43
u__tilde__ = u.anti()
 
44
 
 
45
s = Particle(pdg_code = 3,
 
46
             name = 's',
 
47
             antiname = 's~',
 
48
             spin = 2,
 
49
             color = 3,
 
50
             mass = Param.ZERO,
 
51
             width = Param.ZERO,
 
52
             texname = 's',
 
53
             antitexname = 's',
 
54
             charge = -1/3,
 
55
             LeptonNumber = 0,
 
56
             GhostNumber = 0)
 
57
 
 
58
s__tilde__ = s.anti()
 
59
 
 
60
c = Particle(pdg_code = 4,
 
61
             name = 'c',
 
62
             antiname = 'c~',
 
63
             spin = 2,
 
64
             color = 3,
 
65
             mass = Param.MC,
 
66
             width = Param.ZERO,
 
67
             texname = 'c',
 
68
             antitexname = 'c',
 
69
             charge = 2/3,
 
70
             LeptonNumber = 0,
 
71
             GhostNumber = 0)
 
72
 
 
73
c__tilde__ = c.anti()
 
74
 
 
75
b = Particle(pdg_code = 5,
 
76
             name = 'b',
 
77
             antiname = 'b~',
 
78
             spin = 2,
 
79
             color = 3,
 
80
             mass = Param.MB,
 
81
             width = Param.ZERO,
 
82
             texname = 'b',
 
83
             antitexname = 'b',
 
84
             charge = -1/3,
 
85
             LeptonNumber = 0,
 
86
             GhostNumber = 0)
 
87
 
 
88
b__tilde__ = b.anti()
 
89
 
 
90
t = Particle(pdg_code = 6,
 
91
             name = 't',
 
92
             antiname = 't~',
 
93
             spin = 2,
 
94
             color = 3,
 
95
             mass = Param.MT,
 
96
             width = Param.WT,
 
97
             texname = 't',
 
98
             antitexname = 't',
 
99
             charge = 2/3,
 
100
             LeptonNumber = 0,
 
101
             GhostNumber = 0)
 
102
 
 
103
t__tilde__ = t.anti()
 
104
 
 
105
G = Particle(pdg_code = 21,
 
106
             name = 'G',
 
107
             antiname = 'G',
 
108
             spin = 3,
 
109
             color = 8,
 
110
             mass = Param.ZERO,
 
111
             width = Param.ZERO,
 
112
             texname = 'G',
 
113
             antitexname = 'G',
 
114
             charge = 0,
 
115
             LeptonNumber = 0,
 
116
             GhostNumber = 0)
 
117
 
 
118
gh = Particle(pdg_code = 82,
 
119
             name = 'gh',
 
120
             antiname = 'gh~',
 
121
             spin = -1,
 
122
             color = 8,
 
123
             mass = Param.ZERO,
 
124
             width = Param.ZERO,
 
125
             texname = 'gh',
 
126
             antitexname = 'gh',
 
127
             charge = 0,
 
128
             LeptonNumber = 0,
 
129
             line='dashed',
 
130
             GhostNumber = 1)
 
131
 
 
132
gh__tilde__ = gh.anti()
 
133
 
 
134
# ======================================================================
 
135
# Other particles 
 
136
# ======================================================================
 
137
 
 
138
ve = Particle(pdg_code = 12,
 
139
              name = 've',
 
140
              antiname = 've~',
 
141
              spin = 2,
 
142
              color = 1,
 
143
              mass = Param.ZERO,
 
144
              width = Param.ZERO,
 
145
              texname = 've',
 
146
              antitexname = 've',
 
147
              charge = 0,
 
148
              LeptonNumber = 1,
 
149
              GhostNumber = 0)
 
150
 
 
151
ve__tilde__ = ve.anti()
 
152
 
 
153
vm = Particle(pdg_code = 14,
 
154
              name = 'vm',
 
155
              antiname = 'vm~',
 
156
              spin = 2,
 
157
              color = 1,
 
158
              mass = Param.ZERO,
 
159
              width = Param.ZERO,
 
160
              texname = 'vm',
 
161
              antitexname = 'vm',
 
162
              charge = 0,
 
163
              LeptonNumber = 1,
 
164
              GhostNumber = 0)
 
165
 
 
166
vm__tilde__ = vm.anti()
 
167
 
 
168
vt = Particle(pdg_code = 16,
 
169
              name = 'vt',
 
170
              antiname = 'vt~',
 
171
              spin = 2,
 
172
              color = 1,
 
173
              mass = Param.ZERO,
 
174
              width = Param.ZERO,
 
175
              texname = 'vt',
 
176
              antitexname = 'vt',
 
177
              charge = 0,
 
178
              LeptonNumber = 1,
 
179
              GhostNumber = 0)
 
180
 
 
181
vt__tilde__ = vt.anti()
 
182
 
 
183
A = Particle(pdg_code = 22,
 
184
             name = 'A',
 
185
             antiname = 'A',
 
186
             spin = 3,
 
187
             color = 1,
 
188
             mass = Param.ZERO,
 
189
             width = Param.ZERO,
 
190
             texname = 'A',
 
191
             antitexname = 'A',
 
192
             charge = 0,
 
193
             LeptonNumber = 0,
 
194
             GhostNumber = 0)
 
195
 
 
196
Z = Particle(pdg_code = 23,
 
197
             name = 'Z',
 
198
             antiname = 'Z',
 
199
             spin = 3,
 
200
             color = 1,
 
201
             mass = Param.MZ,
 
202
             width = Param.WZ,
 
203
             texname = 'Z',
 
204
             antitexname = 'Z',
 
205
             charge = 0,
 
206
             LeptonNumber = 0,
 
207
             GhostNumber = 0)
 
208
 
 
209
W__plus__ = Particle(pdg_code = 24,
 
210
                     name = 'W+',
 
211
                     antiname = 'W-',
 
212
                     spin = 3,
 
213
                     color = 1,
 
214
                     mass = Param.MW,
 
215
                     width = Param.WW,
 
216
                     texname = 'W+',
 
217
                     antitexname = 'W+',
 
218
                     charge = 1,
 
219
                     LeptonNumber = 0,
 
220
                     GhostNumber = 0)
 
221
 
 
222
W__minus__ = W__plus__.anti()
 
223
 
 
224
H = Particle(pdg_code = 25,
 
225
             name = 'H',
 
226
             antiname = 'H',
 
227
             spin = 1,
 
228
             color = 1,
 
229
             mass = Param.MH,
 
230
             width = Param.WH,
 
231
             texname = '\\phi',
 
232
             antitexname = '\\phi',
 
233
             charge = 0,
 
234
             LeptonNumber = 0,
 
235
             GhostNumber = 0)
 
236
 
 
237
e__minus__ = Particle(pdg_code = 11,
 
238
                      name = 'e-',
 
239
                      antiname = 'e+',
 
240
                      spin = 2,
 
241
                      color = 1,
 
242
                      mass = Param.Me,
 
243
                      width = Param.ZERO,
 
244
                      texname = 'e-',
 
245
                      antitexname = 'e-',
 
246
                      charge = -1,
 
247
                      LeptonNumber = 1,
 
248
                      GhostNumber = 0)
 
249
 
 
250
e__plus__ = e__minus__.anti()
 
251
 
 
252
m__minus__ = Particle(pdg_code = 13,
 
253
                      name = 'm-',
 
254
                      antiname = 'm+',
 
255
                      spin = 2,
 
256
                      color = 1,
 
257
                      mass = Param.MM,
 
258
                      width = Param.ZERO,
 
259
                      texname = 'm-',
 
260
                      antitexname = 'm-',
 
261
                      charge = -1,
 
262
                      LeptonNumber = 1,
 
263
                      GhostNumber = 0)
 
264
 
 
265
m__plus__ = m__minus__.anti()
 
266
 
 
267
tt__minus__ = Particle(pdg_code = 15,
 
268
                       name = 'tt-',
 
269
                       antiname = 'tt+',
 
270
                       spin = 2,
 
271
                       color = 1,
 
272
                       mass = Param.MTA,
 
273
                       width = Param.WTau,
 
274
                       texname = 'tt-',
 
275
                       antitexname = 'tt-',
 
276
                       charge = -1,
 
277
                       LeptonNumber = 1,
 
278
                       GhostNumber = 0)
 
279
 
 
280
tt__plus__ = tt__minus__.anti()
 
281
 
 
282
 
 
283
y = Particle(pdg_code = 39,
 
284
             name = 'y',
 
285
             antiname = 'y',
 
286
             spin = 5,
 
287
             color = 1,
 
288
             mass = Param.MH,
 
289
             width = Param.WH,
 
290
             texname = 'y',
 
291
             antitexname = 'y',
 
292
             charge = 0,
 
293
             LeptonNumber = 0,
 
294
             GhostNumber = 0)
 
295
 
 
296
 
 
297
# ======================================================================
 
298
# Goldstones 
 
299
# ======================================================================
 
300
 
 
301
G0 = Particle(pdg_code = 250,
 
302
              name = 'G0',
 
303
              antiname = 'G0',
 
304
              spin = 1,
 
305
              color = 1,
 
306
              mass = Param.MZ,
 
307
              width = Param.WZ,
 
308
              texname = 'G0',
 
309
              antitexname = 'G0',
 
310
              GoldstoneBoson = True,
 
311
              charge = 0,
 
312
              GhostNumber = 0,
 
313
              LeptonNumber = 0)
 
314
 
 
315
G__plus__ = Particle(pdg_code = 251,
 
316
                     name = 'G+',
 
317
                     antiname = 'G-',
 
318
                     spin = 1,
 
319
                     color = 1,
 
320
                     mass = Param.MW,
 
321
                     width = Param.WW,
 
322
                     texname = 'G+',
 
323
                     antitexname = 'G-',
 
324
                     GoldstoneBoson = True,
 
325
                     charge = 1,
 
326
                     GhostNumber = 0,
 
327
                     LeptonNumber = 0)
 
328
 
 
329
G__minus__ = G__plus__.anti()
 
330
 
 
331
# Wavefunction renormalization
 
332
 
 
333
b.loop_particles = [[[b,G]]]
 
334
b.counterterm = {(1,0,0):CTParam.bWcft_UV.value}
 
335
 
 
336
c.loop_particles = [[[c,G]]]
 
337
c.counterterm = {(1,0,0):CTParam.cWcft_UV.value}
 
338
 
 
339
t.loop_particles = [[[t,G]]]
 
340
t.counterterm = {(1,0,0):CTParam.tWcft_UV.value}
 
341
 
 
342
G.loop_particles = [[[c]],[[b]],[[t]]]
 
343
G.counterterm = {(1,0,0):CTParam.GWcft_UV_c.value,(1,0,1):CTParam.GWcft_UV_b.value,(1,0,2):CTParam.GWcft_UV_t.value}
 
344
 
 
345
# Set counterterms values
 
346
 
 
347
Param.MB.loop_particles= [[[b,G]]]
 
348
Param.MB.counterterm = {(1,0,0):CTParam.bMass_UV.value}
 
349
 
 
350
Param.MC.loop_particles= [[[c,G]]]
 
351
Param.MC.counterterm = {(1,0,0):CTParam.cMass_UV.value}
 
352
 
 
353
Param.MT.loop_particles= [[[t,G]]]
 
354
Param.MT.counterterm = {(1,0,0):CTParam.tMass_UV.value}
 
355
 
 
356
Param.G.loop_particles = [[[u],[d],[s]],[[c]],[[b]],[[t]],[[G]]],
 
357
Param.G.counterterm = {(1,0,0):CTParam.G_UVq.value,(1,0,1):CTParam.G_UVc.value,(1,0,2):CTParam.G_UVb.value,(1,0,3):CTParam.G_UVt.value,(1,0,4):CTParam.G_UVg.value},