~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to models/loop_sm_FR/lorentz.py

merge with 2.0.0beta4 revision 216

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
# This file was automatically created by FeynRules 1.7.167
 
2
# Mathematica version: 8.0 for Mac OS X x86 (64-bit) (November 6, 2010)
 
3
# Date: Tue 7 May 2013 06:54:14
 
4
 
 
5
 
 
6
from object_library import all_lorentz, Lorentz
 
7
 
 
8
from function_library import complexconjugate, re, im, csc, sec, acsc, asec, cot
 
9
try:
 
10
   import form_factors as ForFac 
 
11
except ImportError:
 
12
   pass
 
13
 
 
14
 
 
15
FF1 = Lorentz(name = 'FF1',
 
16
              spins = [ 2, 2 ],
 
17
              structure = 'P(-1,1)*Gamma(-1,2,1)')
 
18
 
 
19
FF2 = Lorentz(name = 'FF2',
 
20
              spins = [ 2, 2 ],
 
21
              structure = 'ProjM(2,1) + ProjP(2,1)')
 
22
 
 
23
FF3 = Lorentz(name = 'FF3',
 
24
              spins = [ 2, 2 ],
 
25
              structure = 'P(-1,1)*Gamma(-1,2,-2)*ProjM(-2,1) + P(-1,1)*Gamma(-1,2,-2)*ProjP(-2,1)')
 
26
 
 
27
FF4 = Lorentz(name = 'FF4',
 
28
              spins = [ 2, 2 ],
 
29
              structure = '-(P(-1,2)*Gamma(-1,2,-2)*ProjM(-2,1)) + P(-1,1)*Gamma(-1,2,-2)*ProjP(-2,1)')
 
30
 
 
31
FF5 = Lorentz(name = 'FF5',
 
32
              spins = [ 2, 2 ],
 
33
              structure = '-(P(-1,1)*Gamma(-1,2,1)) - P(-1,2)*Gamma(-1,2,-2)*ProjM(-2,1) + P(-1,1)*Gamma(-1,2,-2)*ProjP(-2,1)')
 
34
 
 
35
FF6 = Lorentz(name = 'FF6',
 
36
              spins = [ 2, 2 ],
 
37
              structure = 'P(-1,1)*Gamma(-1,2,-2)*ProjM(-2,1) - 2*P(-1,2)*Gamma(-1,2,-2)*ProjM(-2,1) + 3*P(-1,1)*Gamma(-1,2,-2)*ProjP(-2,1)')
 
38
 
 
39
VV1 = Lorentz(name = 'VV1',
 
40
              spins = [ 3, 3 ],
 
41
              structure = 'Metric(1,2)')
 
42
 
 
43
VV2 = Lorentz(name = 'VV2',
 
44
              spins = [ 3, 3 ],
 
45
              structure = 'P(-1,2)**2*Metric(1,2)')
 
46
 
 
47
VV3 = Lorentz(name = 'VV3',
 
48
              spins = [ 3, 3 ],
 
49
              structure = 'P(1,2)*P(2,1) - P(-1,1)*P(-1,2)*Metric(1,2)')
 
50
 
 
51
VV4 = Lorentz(name = 'VV4',
 
52
              spins = [ 3, 3 ],
 
53
              structure = 'P(1,2)*P(2,1) - (P(-1,1)*P(-1,2)*Metric(1,2))/3. - (P(-1,2)**2*Metric(1,2))/3.')
 
54
 
 
55
VV5 = Lorentz(name = 'VV5',
 
56
              spins = [ 3, 3 ],
 
57
              structure = 'P(1,2)*P(2,1) + (P(-1,1)*P(-1,2)*Metric(1,2))/3. + (P(-1,2)**2*Metric(1,2))/3.')
 
58
 
 
59
VV6 = Lorentz(name = 'VV6',
 
60
              spins = [ 3, 3 ],
 
61
              structure = 'P(1,2)*P(2,1) + P(-1,2)**2*Metric(1,2)')
 
62
 
 
63
VV7 = Lorentz(name = 'VV7',
 
64
              spins = [ 3, 3 ],
 
65
              structure = 'P(-1,1)*P(-1,2)*Metric(1,2) + P(-1,2)**2*Metric(1,2)')
 
66
 
 
67
VV8 = Lorentz(name = 'VV8',
 
68
              spins = [ 3, 3 ],
 
69
              structure = 'P(1,2)*P(2,1) + (3*P(-1,2)**2*Metric(1,2))/2.')
 
70
 
 
71
UUS1 = Lorentz(name = 'UUS1',
 
72
               spins = [ -1, -1, 1 ],
 
73
               structure = '1')
 
74
 
 
75
UUV1 = Lorentz(name = 'UUV1',
 
76
               spins = [ -1, -1, 3 ],
 
77
               structure = 'P(3,2) + P(3,3)')
 
78
 
 
79
SSS1 = Lorentz(name = 'SSS1',
 
80
               spins = [ 1, 1, 1 ],
 
81
               structure = '1')
 
82
 
 
83
FFS1 = Lorentz(name = 'FFS1',
 
84
               spins = [ 2, 2, 1 ],
 
85
               structure = 'Gamma5(2,1)')
 
86
 
 
87
FFS2 = Lorentz(name = 'FFS2',
 
88
               spins = [ 2, 2, 1 ],
 
89
               structure = 'Identity(2,1)')
 
90
 
 
91
FFS3 = Lorentz(name = 'FFS3',
 
92
               spins = [ 2, 2, 1 ],
 
93
               structure = 'ProjM(2,1)')
 
94
 
 
95
FFS4 = Lorentz(name = 'FFS4',
 
96
               spins = [ 2, 2, 1 ],
 
97
               structure = 'ProjM(2,1) - ProjP(2,1)')
 
98
 
 
99
FFS5 = Lorentz(name = 'FFS5',
 
100
               spins = [ 2, 2, 1 ],
 
101
               structure = 'ProjP(2,1)')
 
102
 
 
103
FFS6 = Lorentz(name = 'FFS6',
 
104
               spins = [ 2, 2, 1 ],
 
105
               structure = 'ProjM(2,1) + ProjP(2,1)')
 
106
 
 
107
FFV1 = Lorentz(name = 'FFV1',
 
108
               spins = [ 2, 2, 3 ],
 
109
               structure = 'Gamma(3,2,1)')
 
110
 
 
111
FFV2 = Lorentz(name = 'FFV2',
 
112
               spins = [ 2, 2, 3 ],
 
113
               structure = 'Gamma(3,2,-1)*ProjM(-1,1)')
 
114
 
 
115
FFV3 = Lorentz(name = 'FFV3',
 
116
               spins = [ 2, 2, 3 ],
 
117
               structure = 'Gamma(3,2,-1)*ProjM(-1,1) - 2*Gamma(3,2,-1)*ProjP(-1,1)')
 
118
 
 
119
FFV4 = Lorentz(name = 'FFV4',
 
120
               spins = [ 2, 2, 3 ],
 
121
               structure = 'Gamma(3,2,1) - Gamma(3,2,-1)*ProjM(-1,1) - Gamma(3,2,-1)*ProjP(-1,1)')
 
122
 
 
123
FFV5 = Lorentz(name = 'FFV5',
 
124
               spins = [ 2, 2, 3 ],
 
125
               structure = 'Gamma(3,2,1) - (4*Gamma(3,2,-1)*ProjM(-1,1))/13. - (4*Gamma(3,2,-1)*ProjP(-1,1))/13.')
 
126
 
 
127
FFV6 = Lorentz(name = 'FFV6',
 
128
               spins = [ 2, 2, 3 ],
 
129
               structure = 'Gamma(3,2,1) + (8*Gamma(3,2,-1)*ProjM(-1,1))/13. + (8*Gamma(3,2,-1)*ProjP(-1,1))/13.')
 
130
 
 
131
FFV7 = Lorentz(name = 'FFV7',
 
132
               spins = [ 2, 2, 3 ],
 
133
               structure = 'Gamma(3,2,-1)*ProjM(-1,1) + Gamma(3,2,-1)*ProjP(-1,1)')
 
134
 
 
135
FFV8 = Lorentz(name = 'FFV8',
 
136
               spins = [ 2, 2, 3 ],
 
137
               structure = 'Gamma(3,2,-1)*ProjM(-1,1) + 2*Gamma(3,2,-1)*ProjP(-1,1)')
 
138
 
 
139
FFV9 = Lorentz(name = 'FFV9',
 
140
               spins = [ 2, 2, 3 ],
 
141
               structure = 'Gamma(3,2,1) + 2*Gamma(3,2,-1)*ProjM(-1,1) + 2*Gamma(3,2,-1)*ProjP(-1,1)')
 
142
 
 
143
FFV10 = Lorentz(name = 'FFV10',
 
144
                spins = [ 2, 2, 3 ],
 
145
                structure = 'Gamma(3,2,-1)*ProjM(-1,1) + 4*Gamma(3,2,-1)*ProjP(-1,1)')
 
146
 
 
147
VSS1 = Lorentz(name = 'VSS1',
 
148
               spins = [ 3, 1, 1 ],
 
149
               structure = 'P(1,2) - P(1,3)')
 
150
 
 
151
VVS1 = Lorentz(name = 'VVS1',
 
152
               spins = [ 3, 3, 1 ],
 
153
               structure = 'Metric(1,2)')
 
154
 
 
155
VVV1 = Lorentz(name = 'VVV1',
 
156
               spins = [ 3, 3, 3 ],
 
157
               structure = 'P(3,1)*Metric(1,2) - 3*P(3,2)*Metric(1,2) - P(3,3)*Metric(1,2) - P(2,1)*Metric(1,3) + P(2,2)*Metric(1,3) + 3*P(2,3)*Metric(1,3) + 2*P(1,2)*Metric(2,3) - 2*P(1,3)*Metric(2,3)')
 
158
 
 
159
VVV2 = Lorentz(name = 'VVV2',
 
160
               spins = [ 3, 3, 3 ],
 
161
               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)')
 
162
 
 
163
VVV3 = Lorentz(name = 'VVV3',
 
164
               spins = [ 3, 3, 3 ],
 
165
               structure = 'P(3,1)*Metric(1,2) - (11*P(3,2)*Metric(1,2))/17. + (3*P(3,3)*Metric(1,2))/17. - P(2,1)*Metric(1,3) - (3*P(2,2)*Metric(1,3))/17. + (11*P(2,3)*Metric(1,3))/17. + (14*P(1,2)*Metric(2,3))/17. - (14*P(1,3)*Metric(2,3))/17.')
 
166
 
 
167
SSSS1 = Lorentz(name = 'SSSS1',
 
168
                spins = [ 1, 1, 1, 1 ],
 
169
                structure = '1')
 
170
 
 
171
VVSS1 = Lorentz(name = 'VVSS1',
 
172
                spins = [ 3, 3, 1, 1 ],
 
173
                structure = 'Metric(1,2)')
 
174
 
 
175
VVVV1 = Lorentz(name = 'VVVV1',
 
176
                spins = [ 3, 3, 3, 3 ],
 
177
                structure = 'Epsilon(1,2,3,4)')
 
178
 
 
179
VVVV2 = Lorentz(name = 'VVVV2',
 
180
                spins = [ 3, 3, 3, 3 ],
 
181
                structure = 'Metric(1,4)*Metric(2,3)')
 
182
 
 
183
VVVV3 = Lorentz(name = 'VVVV3',
 
184
                spins = [ 3, 3, 3, 3 ],
 
185
                structure = 'Metric(1,3)*Metric(2,4)')
 
186
 
 
187
VVVV4 = Lorentz(name = 'VVVV4',
 
188
                spins = [ 3, 3, 3, 3 ],
 
189
                structure = 'Metric(1,4)*Metric(2,3) - Metric(1,3)*Metric(2,4)')
 
190
 
 
191
VVVV5 = Lorentz(name = 'VVVV5',
 
192
                spins = [ 3, 3, 3, 3 ],
 
193
                structure = 'Metric(1,2)*Metric(3,4)')
 
194
 
 
195
VVVV6 = Lorentz(name = 'VVVV6',
 
196
                spins = [ 3, 3, 3, 3 ],
 
197
                structure = 'Metric(1,4)*Metric(2,3) + Metric(1,3)*Metric(2,4) - 2*Metric(1,2)*Metric(3,4)')
 
198
 
 
199
VVVV7 = Lorentz(name = 'VVVV7',
 
200
                spins = [ 3, 3, 3, 3 ],
 
201
                structure = 'Metric(1,4)*Metric(2,3) - Metric(1,2)*Metric(3,4)')
 
202
 
 
203
VVVV8 = Lorentz(name = 'VVVV8',
 
204
                spins = [ 3, 3, 3, 3 ],
 
205
                structure = 'Metric(1,3)*Metric(2,4) - Metric(1,2)*Metric(3,4)')
 
206
 
 
207
VVVV9 = Lorentz(name = 'VVVV9',
 
208
                spins = [ 3, 3, 3, 3 ],
 
209
                structure = 'Metric(1,4)*Metric(2,3) - (Metric(1,3)*Metric(2,4))/2. - (Metric(1,2)*Metric(3,4))/2.')
 
210
 
 
211
VVVV10 = Lorentz(name = 'VVVV10',
 
212
                 spins = [ 3, 3, 3, 3 ],
 
213
                 structure = 'Metric(1,4)*Metric(2,3) + Metric(1,3)*Metric(2,4) + Metric(1,2)*Metric(3,4)')
 
214