-
Notifications
You must be signed in to change notification settings - Fork 0
/
lorentz.py
199 lines (146 loc) · 9.51 KB
/
lorentz.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
# This file was automatically created by FeynRules $Revision: 1058$
# Mathematica version: 8.0 for Mac OS X x86 (64-bit) (November 6, 2010)
# Date: Fri 6 Jul 2012 14:58:26
from object_library import all_lorentz, Lorentz
from function_library import complexconjugate, re, im, csc, sec, acsc, asec
UUS1 = Lorentz(name = 'UUS1',
spins = [ -1, -1, 1 ],
structure = '1')
UUV1 = Lorentz(name = 'UUV1',
spins = [ -1, -1, 3 ],
structure = 'P(3,2) + P(3,3)')
SSS1 = Lorentz(name = 'SSS1',
spins = [ 1, 1, 1 ],
structure = '1')
FFS1 = Lorentz(name = 'FFS1',
spins = [ 2, 2, 1 ],
structure = 'Gamma5(2,1)')
FFS2 = Lorentz(name = 'FFS2',
spins = [ 2, 2, 1 ],
structure = 'Identity(2,1)')
FFS3 = Lorentz(name = 'FFS3',
spins = [ 2, 2, 1 ],
structure = 'ProjM(2,1)')
FFS4 = Lorentz(name = 'FFS4',
spins = [ 2, 2, 1 ],
structure = 'ProjP(2,1)')
FFV1 = Lorentz(name = 'FFV1',
spins = [ 2, 2, 3 ],
structure = 'Gamma(3,2,1)')
FFV2 = Lorentz(name = 'FFV2',
spins = [ 2, 2, 3 ],
structure = 'Gamma(3,2,-1)*ProjM(-1,1)')
FFV3 = Lorentz(name = 'FFV3',
spins = [ 2, 2, 3 ],
structure = 'Gamma(3,2,-1)*ProjM(-1,1) - 2*Gamma(3,2,-1)*ProjP(-1,1)')
FFV4 = Lorentz(name = 'FFV4',
spins = [ 2, 2, 3 ],
structure = 'Gamma(3,2,-1)*ProjM(-1,1) + 2*Gamma(3,2,-1)*ProjP(-1,1)')
FFV5 = Lorentz(name = 'FFV5',
spins = [ 2, 2, 3 ],
structure = 'Gamma(3,2,-1)*ProjM(-1,1) + 4*Gamma(3,2,-1)*ProjP(-1,1)')
VSS1 = Lorentz(name = 'VSS1',
spins = [ 3, 1, 1 ],
structure = 'P(1,2) - P(1,3)')
VVS1 = Lorentz(name = 'VVS1',
spins = [ 3, 3, 1 ],
structure = 'Metric(1,2)')
VVV1 = Lorentz(name = 'VVV1',
spins = [ 3, 3, 3 ],
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)')
SSSS1 = Lorentz(name = 'SSSS1',
spins = [ 1, 1, 1, 1 ],
structure = '1')
FFFF1 = Lorentz(name = 'FFFF1',
spins = [ 2, 2, 2, 2 ],
structure = 'Gamma(-1,2,1)*Gamma(-1,4,3)')
FFFF2 = Lorentz(name = 'FFFF2',
spins = [ 2, 2, 2, 2 ],
structure = 'Identity(1,2)*Identity(3,4)')
FFFF3 = Lorentz(name = 'FFFF3',
spins = [ 2, 2, 2, 2 ],
structure = 'Gamma(-2,-4,3)*Gamma(-2,-3,1)*Gamma(-1,2,-3)*Gamma(-1,4,-4) - Gamma(-1,-2,3)*Gamma(-1,4,-2)*Identity(1,2) - Gamma(-1,-2,1)*Gamma(-1,2,-2)*Identity(3,4) + 4*Identity(1,2)*Identity(3,4)')
FFFF4 = Lorentz(name = 'FFFF4',
spins = [ 2, 2, 2, 2 ],
structure = 'Identity(3,4)*ProjM(2,1) - Identity(3,4)*ProjP(2,1)')
FFFF5 = Lorentz(name = 'FFFF5',
spins = [ 2, 2, 2, 2 ],
structure = 'Identity(1,2)*ProjM(4,3) - Identity(1,2)*ProjP(4,3)')
FFFF6 = Lorentz(name = 'FFFF6',
spins = [ 2, 2, 2, 2 ],
structure = 'ProjM(2,1)*ProjM(4,3) - ProjM(4,3)*ProjP(2,1) - ProjM(2,1)*ProjP(4,3) + ProjP(2,1)*ProjP(4,3)')
FFFF7 = Lorentz(name = 'FFFF7',
spins = [ 2, 2, 2, 2 ],
structure = 'Gamma(-1,2,-2)*Gamma(-1,4,3)*ProjM(-2,1) - Gamma(-1,2,-2)*Gamma(-1,4,3)*ProjP(-2,1)')
FFFF8 = Lorentz(name = 'FFFF8',
spins = [ 2, 2, 2, 2 ],
structure = 'Gamma(-1,2,1)*Gamma(-1,4,-2)*ProjM(-2,3) - Gamma(-1,2,1)*Gamma(-1,4,-2)*ProjP(-2,3)')
FFFF9 = Lorentz(name = 'FFFF9',
spins = [ 2, 2, 2, 2 ],
structure = 'Gamma(-1,2,-2)*Gamma(-1,4,-3)*ProjM(-3,3)*ProjM(-2,1) - Gamma(-1,2,-3)*Gamma(-1,4,-2)*ProjM(-2,3)*ProjP(-3,1) - Gamma(-1,2,-2)*Gamma(-1,4,-3)*ProjM(-2,1)*ProjP(-3,3) + Gamma(-1,2,-2)*Gamma(-1,4,-3)*ProjP(-3,3)*ProjP(-2,1)')
FFFF10 = Lorentz(name = 'FFFF10',
spins = [ 2, 2, 2, 2 ],
structure = '-(Gamma(-2,-5,-4)*Gamma(-2,-3,3)*Gamma(-1,2,-5)*Gamma(-1,4,-3)*ProjM(-4,1)) + Gamma(-1,-3,-2)*Gamma(-1,2,-3)*Identity(3,4)*ProjM(-2,1) + Gamma(-1,-2,3)*Gamma(-1,4,-2)*ProjM(2,1) - 4*Identity(3,4)*ProjM(2,1) + Gamma(-2,-5,-4)*Gamma(-2,-3,3)*Gamma(-1,2,-5)*Gamma(-1,4,-3)*ProjP(-4,1) - Gamma(-1,-3,-2)*Gamma(-1,2,-3)*Identity(3,4)*ProjP(-2,1) - Gamma(-1,-2,3)*Gamma(-1,4,-2)*ProjP(2,1) + 4*Identity(3,4)*ProjP(2,1)')
FFFF11 = Lorentz(name = 'FFFF11',
spins = [ 2, 2, 2, 2 ],
structure = '-(Gamma(-2,-5,-4)*Gamma(-2,-3,1)*Gamma(-1,2,-3)*Gamma(-1,4,-5)*ProjM(-4,3)) + Gamma(-1,-3,-2)*Gamma(-1,4,-3)*Identity(1,2)*ProjM(-2,3) + Gamma(-1,-2,1)*Gamma(-1,2,-2)*ProjM(4,3) - 4*Identity(1,2)*ProjM(4,3) + Gamma(-2,-5,-4)*Gamma(-2,-3,1)*Gamma(-1,2,-3)*Gamma(-1,4,-5)*ProjP(-4,3) - Gamma(-1,-3,-2)*Gamma(-1,4,-3)*Identity(1,2)*ProjP(-2,3) - Gamma(-1,-2,1)*Gamma(-1,2,-2)*ProjP(4,3) + 4*Identity(1,2)*ProjP(4,3)')
FFVV1 = Lorentz(name = 'FFVV1',
spins = [ 2, 2, 3, 3 ],
structure = '-(Epsilon(3,4,-1,-2)*P(-2,4)*P(-1,3)*Identity(2,1)) + Epsilon(3,4,-1,-2)*P(-2,3)*P(-1,4)*Identity(2,1)')
FFVV2 = Lorentz(name = 'FFVV2',
spins = [ 2, 2, 3, 3 ],
structure = 'P(3,4)*P(4,3)*Identity(2,1) - P(-1,3)*P(-1,4)*Identity(2,1)*Metric(3,4)')
FFVV3 = Lorentz(name = 'FFVV3',
spins = [ 2, 2, 3, 3 ],
structure = '-(Epsilon(3,4,-1,-2)*P(-2,4)*P(-1,3)*ProjM(2,1)) + Epsilon(3,4,-1,-2)*P(-2,3)*P(-1,4)*ProjM(2,1) + Epsilon(3,4,-1,-2)*P(-2,4)*P(-1,3)*ProjP(2,1) - Epsilon(3,4,-1,-2)*P(-2,3)*P(-1,4)*ProjP(2,1)')
FFVV4 = Lorentz(name = 'FFVV4',
spins = [ 2, 2, 3, 3 ],
structure = 'P(3,4)*P(4,3)*ProjM(2,1) - P(-1,3)*P(-1,4)*Metric(3,4)*ProjM(2,1) - P(3,4)*P(4,3)*ProjP(2,1) + P(-1,3)*P(-1,4)*Metric(3,4)*ProjP(2,1)')
VVSS1 = Lorentz(name = 'VVSS1',
spins = [ 3, 3, 1, 1 ],
structure = 'Metric(1,2)')
VVVV1 = Lorentz(name = 'VVVV1',
spins = [ 3, 3, 3, 3 ],
structure = 'Metric(1,4)*Metric(2,3) - Metric(1,3)*Metric(2,4)')
VVVV2 = Lorentz(name = 'VVVV2',
spins = [ 3, 3, 3, 3 ],
structure = 'Metric(1,4)*Metric(2,3) + Metric(1,3)*Metric(2,4) - 2*Metric(1,2)*Metric(3,4)')
VVVV3 = Lorentz(name = 'VVVV3',
spins = [ 3, 3, 3, 3 ],
structure = 'Metric(1,4)*Metric(2,3) - Metric(1,2)*Metric(3,4)')
VVVV4 = Lorentz(name = 'VVVV4',
spins = [ 3, 3, 3, 3 ],
structure = 'Metric(1,3)*Metric(2,4) - Metric(1,2)*Metric(3,4)')
VVVV5 = Lorentz(name = 'VVVV5',
spins = [ 3, 3, 3, 3 ],
structure = 'Metric(1,4)*Metric(2,3) - (Metric(1,3)*Metric(2,4))/2. - (Metric(1,2)*Metric(3,4))/2.')
FFVVV1 = Lorentz(name = 'FFVVV1',
spins = [ 2, 2, 3, 3, 3 ],
structure = 'Epsilon(3,4,5,-1)*P(-1,3)*Identity(2,1) + Epsilon(3,4,5,-1)*P(-1,4)*Identity(2,1) + Epsilon(3,4,5,-1)*P(-1,5)*Identity(2,1)')
FFVVV2 = Lorentz(name = 'FFVVV2',
spins = [ 2, 2, 3, 3, 3 ],
structure = 'P(5,3)*Identity(2,1)*Metric(3,4) - P(5,4)*Identity(2,1)*Metric(3,4) - P(4,3)*Identity(2,1)*Metric(3,5) + P(4,5)*Identity(2,1)*Metric(3,5) + P(3,4)*Identity(2,1)*Metric(4,5) - P(3,5)*Identity(2,1)*Metric(4,5)')
FFVVV3 = Lorentz(name = 'FFVVV3',
spins = [ 2, 2, 3, 3, 3 ],
structure = 'Epsilon(3,4,5,-1)*P(-1,3)*ProjM(2,1) + Epsilon(3,4,5,-1)*P(-1,4)*ProjM(2,1) + Epsilon(3,4,5,-1)*P(-1,5)*ProjM(2,1) - Epsilon(3,4,5,-1)*P(-1,3)*ProjP(2,1) - Epsilon(3,4,5,-1)*P(-1,4)*ProjP(2,1) - Epsilon(3,4,5,-1)*P(-1,5)*ProjP(2,1)')
FFVVV4 = Lorentz(name = 'FFVVV4',
spins = [ 2, 2, 3, 3, 3 ],
structure = 'P(5,3)*Metric(3,4)*ProjM(2,1) - P(5,4)*Metric(3,4)*ProjM(2,1) - P(4,3)*Metric(3,5)*ProjM(2,1) + P(4,5)*Metric(3,5)*ProjM(2,1) + P(3,4)*Metric(4,5)*ProjM(2,1) - P(3,5)*Metric(4,5)*ProjM(2,1) - P(5,3)*Metric(3,4)*ProjP(2,1) + P(5,4)*Metric(3,4)*ProjP(2,1) + P(4,3)*Metric(3,5)*ProjP(2,1) - P(4,5)*Metric(3,5)*ProjP(2,1) - P(3,4)*Metric(4,5)*ProjP(2,1) + P(3,5)*Metric(4,5)*ProjP(2,1)')
FFVVVV1 = Lorentz(name = 'FFVVVV1',
spins = [ 2, 2, 3, 3, 3, 3 ],
structure = 'Identity(2,1)*Metric(3,6)*Metric(4,5) - Identity(2,1)*Metric(3,5)*Metric(4,6)')
FFVVVV2 = Lorentz(name = 'FFVVVV2',
spins = [ 2, 2, 3, 3, 3, 3 ],
structure = 'Identity(2,1)*Metric(3,6)*Metric(4,5) - Identity(2,1)*Metric(3,4)*Metric(5,6)')
FFVVVV3 = Lorentz(name = 'FFVVVV3',
spins = [ 2, 2, 3, 3, 3, 3 ],
structure = 'Identity(2,1)*Metric(3,5)*Metric(4,6) - Identity(2,1)*Metric(3,4)*Metric(5,6)')
FFVVVV4 = Lorentz(name = 'FFVVVV4',
spins = [ 2, 2, 3, 3, 3, 3 ],
structure = 'Metric(3,6)*Metric(4,5)*ProjM(2,1) - Metric(3,5)*Metric(4,6)*ProjM(2,1) - Metric(3,6)*Metric(4,5)*ProjP(2,1) + Metric(3,5)*Metric(4,6)*ProjP(2,1)')
FFVVVV5 = Lorentz(name = 'FFVVVV5',
spins = [ 2, 2, 3, 3, 3, 3 ],
structure = 'Metric(3,6)*Metric(4,5)*ProjM(2,1) - Metric(3,4)*Metric(5,6)*ProjM(2,1) - Metric(3,6)*Metric(4,5)*ProjP(2,1) + Metric(3,4)*Metric(5,6)*ProjP(2,1)')
FFVVVV6 = Lorentz(name = 'FFVVVV6',
spins = [ 2, 2, 3, 3, 3, 3 ],
structure = 'Metric(3,5)*Metric(4,6)*ProjM(2,1) - Metric(3,4)*Metric(5,6)*ProjM(2,1) - Metric(3,5)*Metric(4,6)*ProjP(2,1) + Metric(3,4)*Metric(5,6)*ProjP(2,1)')