DevelopmentPage/BeyondProjects: RSLorentz.py

File RSLorentz.py, 11.1 KB (added by Will Link, 14 years ago)

Commented out vertices that reader doesn't like

Line 
1from helasamp_object import *
2Lorentz = []
3
4
5Lorentz.append( { 'Name':'Lorentz1',
6 'Particles':[5, 1, 1],
7 'Structure':[P(31,3)*P(32,2) + P(31,2)*P(32,3) - P('a',2)*P('a',3)*Metric(31,32), Metric(31,32), Metric(31,32)]
8 })
9Lorentz.append( { 'Name':'Lorentz4',
10 'Particles':[5, 3, 3],
11 'Structure':[P(31,3)*P(32,2)*Metric(2,3) + P(31,2)*P(32,3)*Metric(2,3) - P(3,2)*P(32,3)*Metric(2,31) - P(3,2)*P(31,3)*Metric(2,32) - P(2,3)*P(32,2)*Metric(3,31) + P('a',2)*P('a',3)*Metric(2,32)*Metric(3,31) - P(2,3)*P(31,2)*Metric(3,32) + P('a',2)*P('a',3)*Metric(2,31)*Metric(3,32) + P(2,3)*P(3,2)*Metric(31,32) - P('a',2)*P('a',3)*Metric(2,3)*Metric(31,32), Metric(2,32)*Metric(3,31) + Metric(2,31)*Metric(3,32) - Metric(2,3)*Metric(31,32)]
12 })
13Lorentz.append( { 'Name':'Lorentz6',
14 'Particles':[3, 3, 3, 5],
15 'Structure':[P(42,2)*Metric(1,41)*Metric(2,3) - P(42,3)*Metric(1,41)*Metric(2,3) + P(41,2)*Metric(1,42)*Metric(2,3) - P(41,3)*Metric(1,42)*Metric(2,3) - P(42,1)*Metric(1,3)*Metric(2,41) + P(42,3)*Metric(1,3)*Metric(2,41) + P(3,1)*Metric(1,42)*Metric(2,41) - P(3,2)*Metric(1,42)*Metric(2,41) - P(41,1)*Metric(1,3)*Metric(2,42) + P(41,3)*Metric(1,3)*Metric(2,42) + P(3,1)*Metric(1,41)*Metric(2,42) - P(3,2)*Metric(1,41)*Metric(2,42) + P(42,1)*Metric(1,2)*Metric(3,41) - P(42,2)*Metric(1,2)*Metric(3,41) - P(2,1)*Metric(1,42)*Metric(3,41) + P(2,3)*Metric(1,42)*Metric(3,41) + P(1,2)*Metric(2,42)*Metric(3,41) - P(1,3)*Metric(2,42)*Metric(3,41) + P(41,1)*Metric(1,2)*Metric(3,42) - P(41,2)*Metric(1,2)*Metric(3,42) - P(2,1)*Metric(1,41)*Metric(3,42) + P(2,3)*Metric(1,41)*Metric(3,42) + P(1,2)*Metric(2,41)*Metric(3,42) - P(1,3)*Metric(2,41)*Metric(3,42) - P(3,1)*Metric(1,2)*Metric(41,42) + P(3,2)*Metric(1,2)*Metric(41,42) + P(2,1)*Metric(1,3)*Metric(41,42) - P(2,3)*Metric(1,3)*Metric(41,42) - P(1,2)*Metric(2,3)*Metric(41,42) + P(1,3)*Metric(2,3)*Metric(41,42)]
16 })
17Lorentz.append( { 'Name':'Lorentz7',
18 'Particles':[3, 3, 3, 3, 5],
19 'Structure':[Metric(1,52)*Metric(2,51)*Metric(3,4) + Metric(1,51)*Metric(2,52)*Metric(3,4) - Metric(1,52)*Metric(2,4)*Metric(3,51) - Metric(1,51)*Metric(2,4)*Metric(3,52) - Metric(1,3)*Metric(2,52)*Metric(4,51) + Metric(1,2)*Metric(3,52)*Metric(4,51) - Metric(1,3)*Metric(2,51)*Metric(4,52) + Metric(1,2)*Metric(3,51)*Metric(4,52) + Metric(1,3)*Metric(2,4)*Metric(51,52) - Metric(1,2)*Metric(3,4)*Metric(51,52), Metric(1,52)*Metric(2,51)*Metric(3,4) + Metric(1,51)*Metric(2,52)*Metric(3,4) - Metric(1,4)*Metric(2,52)*Metric(3,51) - Metric(1,4)*Metric(2,51)*Metric(3,52) - Metric(1,52)*Metric(2,3)*Metric(4,51) + Metric(1,2)*Metric(3,52)*Metric(4,51) - Metric(1,51)*Metric(2,3)*Metric(4,52) + Metric(1,2)*Metric(3,51)*Metric(4,52) + Metric(1,4)*Metric(2,3)*Metric(51,52) - Metric(1,2)*Metric(3,4)*Metric(51,52), Metric(1,52)*Metric(2,4)*Metric(3,51) - Metric(1,4)*Metric(2,52)*Metric(3,51) + Metric(1,51)*Metric(2,4)*Metric(3,52) - Metric(1,4)*Metric(2,51)*Metric(3,52) - Metric(1,52)*Metric(2,3)*Metric(4,51) + Metric(1,3)*Metric(2,52)*Metric(4,51) - Metric(1,51)*Metric(2,3)*Metric(4,52) + Metric(1,3)*Metric(2,51)*Metric(4,52) + Metric(1,4)*Metric(2,3)*Metric(51,52) - Metric(1,3)*Metric(2,4)*Metric(51,52)]
20 })
21Lorentz.append( { 'Name':'Lorentz8',
22 'Particles':[2, 2, 5],
23 'Structure':[P(12,1)*Gamma(11,2,1) - P(12,2)*Gamma(11,2,1) + P(11,1)*Gamma(12,2,1) - P(11,2)*Gamma(12,2,1) - 2*P('a',1)*Gamma('a',2,1)*Metric(11,12) + 2*P('a',2)*Gamma('a',2,1)*Metric(11,12), Metric(11,12)*ProjM(2,1) + Metric(11,12)*ProjP(2,1)]
24 })
25Lorentz.append( { 'Name':'Lorentz15',
26 'Particles':[5, 2, 2],
27 'Structure':[P(12,2)*Gamma(11,3,2) - P(12,3)*Gamma(11,3,2) + P(11,2)*Gamma(12,3,2) - P(11,3)*Gamma(12,3,2) - 2*P('a',2)*Gamma('a',3,2)*Metric(11,12) + 2*P('a',3)*Gamma('a',3,2)*Metric(11,12), Metric(11,12)*ProjM(3,2) + Metric(11,12)*ProjP(3,2)]
28 })
29Lorentz.append( { 'Name':'Lorentz20',
30 'Particles':[5, 1, 1, 1, 1],
31 'Structure':[Metric(51,52)]
32 })
33Lorentz.append( { 'Name':'Lorentz21',
34 'Particles':[5, 1, 1, 1],
35 'Structure':[Metric(41,42)]
36 })
37Lorentz.append( { 'Name':'Lorentz22',
38 'Particles':[2, 2, 3, 5],
39 'Structure':[Gamma(22,2,1)*Metric(3,21) + Gamma(21,2,1)*Metric(3,22) - 2*Gamma(3,2,1)*Metric(21,22)]
40 })
41Lorentz.append( { 'Name':'Lorentz23',
42 'Particles':[3, 5, 2, 2],
43 'Structure':[Gamma(22,4,3)*Metric(1,21) + Gamma(21,4,3)*Metric(1,22) - 2*Gamma(1,4,3)*Metric(21,22)]
44 })
45Lorentz.append( { 'Name':'Lorentz28',
46 'Particles':[3, 5, 3, 3],
47 'Structure':[P(32,3)*Metric(1,31)*Metric(3,4) - P(32,4)*Metric(1,31)*Metric(3,4) + P(31,3)*Metric(1,32)*Metric(3,4) - P(31,4)*Metric(1,32)*Metric(3,4) - P(32,1)*Metric(1,4)*Metric(3,31) + P(32,4)*Metric(1,4)*Metric(3,31) + P(4,1)*Metric(1,32)*Metric(3,31) - P(4,3)*Metric(1,32)*Metric(3,31) - P(31,1)*Metric(1,4)*Metric(3,32) + P(31,4)*Metric(1,4)*Metric(3,32) + P(4,1)*Metric(1,31)*Metric(3,32) - P(4,3)*Metric(1,31)*Metric(3,32) + P(32,1)*Metric(1,3)*Metric(4,31) - P(32,3)*Metric(1,3)*Metric(4,31) - P(3,1)*Metric(1,32)*Metric(4,31) + P(3,4)*Metric(1,32)*Metric(4,31) + P(1,3)*Metric(3,32)*Metric(4,31) - P(1,4)*Metric(3,32)*Metric(4,31) + P(31,1)*Metric(1,3)*Metric(4,32) - P(31,3)*Metric(1,3)*Metric(4,32) - P(3,1)*Metric(1,31)*Metric(4,32) + P(3,4)*Metric(1,31)*Metric(4,32) + P(1,3)*Metric(3,31)*Metric(4,32) - P(1,4)*Metric(3,31)*Metric(4,32) - P(4,1)*Metric(1,3)*Metric(31,32) + P(4,3)*Metric(1,3)*Metric(31,32) + P(3,1)*Metric(1,4)*Metric(31,32) - P(3,4)*Metric(1,4)*Metric(31,32) - P(1,3)*Metric(3,4)*Metric(31,32) + P(1,4)*Metric(3,4)*Metric(31,32)]
48 })
49Lorentz.append( { 'Name':'Lorentz29',
50 'Particles':[5, 1, 1, 3, 3],
51 'Structure':[Metric(4,32)*Metric(5,31) + Metric(4,31)*Metric(5,32) - Metric(4,5)*Metric(31,32)]
52 })
53Lorentz.append( { 'Name':'Lorentz30',
54 'Particles':[5, 1, 3, 3],
55 'Structure':[Metric(3,22)*Metric(4,21) + Metric(3,21)*Metric(4,22) - Metric(3,4)*Metric(21,22)]
56 })
57Lorentz.append( { 'Name':'Lorentz31',
58 'Particles':[3, 3, 5, 3, 3],
59 'Structure':[Metric(1,32)*Metric(2,31)*Metric(4,5) + Metric(1,31)*Metric(2,32)*Metric(4,5) - (Metric(1,32)*Metric(2,5)*Metric(4,31))/2. - (Metric(1,5)*Metric(2,32)*Metric(4,31))/2. - (Metric(1,31)*Metric(2,5)*Metric(4,32))/2. - (Metric(1,5)*Metric(2,31)*Metric(4,32))/2. - (Metric(1,32)*Metric(2,4)*Metric(5,31))/2. - (Metric(1,4)*Metric(2,32)*Metric(5,31))/2. + Metric(1,2)*Metric(4,32)*Metric(5,31) - (Metric(1,31)*Metric(2,4)*Metric(5,32))/2. - (Metric(1,4)*Metric(2,31)*Metric(5,32))/2. + Metric(1,2)*Metric(4,31)*Metric(5,32) + (Metric(1,5)*Metric(2,4)*Metric(31,32))/2. + (Metric(1,4)*Metric(2,5)*Metric(31,32))/2. - Metric(1,2)*Metric(4,5)*Metric(31,32)]
60 })
61Lorentz.append( { 'Name':'Lorentz32',
62 'Particles':[5, 3, 3, 3],
63 'Structure':[P(22,3)*Metric(2,21)*Metric(3,4) - P(22,4)*Metric(2,21)*Metric(3,4) + P(21,3)*Metric(2,22)*Metric(3,4) - P(21,4)*Metric(2,22)*Metric(3,4) - P(22,2)*Metric(2,4)*Metric(3,21) + P(22,4)*Metric(2,4)*Metric(3,21) + P(4,2)*Metric(2,22)*Metric(3,21) - P(4,3)*Metric(2,22)*Metric(3,21) - P(21,2)*Metric(2,4)*Metric(3,22) + P(21,4)*Metric(2,4)*Metric(3,22) + P(4,2)*Metric(2,21)*Metric(3,22) - P(4,3)*Metric(2,21)*Metric(3,22) + P(22,2)*Metric(2,3)*Metric(4,21) - P(22,3)*Metric(2,3)*Metric(4,21) - P(3,2)*Metric(2,22)*Metric(4,21) + P(3,4)*Metric(2,22)*Metric(4,21) + P(2,3)*Metric(3,22)*Metric(4,21) - P(2,4)*Metric(3,22)*Metric(4,21) + P(21,2)*Metric(2,3)*Metric(4,22) - P(21,3)*Metric(2,3)*Metric(4,22) - P(3,2)*Metric(2,21)*Metric(4,22) + P(3,4)*Metric(2,21)*Metric(4,22) + P(2,3)*Metric(3,21)*Metric(4,22) - P(2,4)*Metric(3,21)*Metric(4,22) - P(4,2)*Metric(2,3)*Metric(21,22) + P(4,3)*Metric(2,3)*Metric(21,22) + P(3,2)*Metric(2,4)*Metric(21,22) - P(3,4)*Metric(2,4)*Metric(21,22) - P(2,3)*Metric(3,4)*Metric(21,22) + P(2,4)*Metric(3,4)*Metric(21,22)]
64 })
65Lorentz.append( { 'Name':'Lorentz33',
66 'Particles':[5, 3, 3, 3, 3],
67 'Structure':[Metric(2,12)*Metric(3,11)*Metric(4,5) + Metric(2,11)*Metric(3,12)*Metric(4,5) - (Metric(2,12)*Metric(3,5)*Metric(4,11))/2. - (Metric(2,5)*Metric(3,12)*Metric(4,11))/2. - (Metric(2,11)*Metric(3,5)*Metric(4,12))/2. - (Metric(2,5)*Metric(3,11)*Metric(4,12))/2. - (Metric(2,12)*Metric(3,4)*Metric(5,11))/2. - (Metric(2,4)*Metric(3,12)*Metric(5,11))/2. + Metric(2,3)*Metric(4,12)*Metric(5,11) - (Metric(2,11)*Metric(3,4)*Metric(5,12))/2. - (Metric(2,4)*Metric(3,11)*Metric(5,12))/2. + Metric(2,3)*Metric(4,11)*Metric(5,12) + (Metric(2,5)*Metric(3,4)*Metric(11,12))/2. + (Metric(2,4)*Metric(3,5)*Metric(11,12))/2. - Metric(2,3)*Metric(4,5)*Metric(11,12)]
68 })
69Lorentz.append( { 'Name':'Lorentz34',
70 'Particles':[2, 2, 5, 1],
71 'Structure':[Identity(1,2)*Metric(21,22)]
72 })
73Lorentz.append( { 'Name':'Lorentz35',
74 'Particles':[5, 1, 2, 2],
75 'Structure':[Identity(3,4)*Metric(21,22)]
76 })
77Lorentz.append( { 'Name':'Lorentz44',
78 'Particles':[5, 1, 1, 3, 3],
79 'Structure':[Metric(4,32)*Metric(5,31) + Metric(4,31)*Metric(5,32) - Metric(4,5)*Metric(31,32), Metric(4,32)*Metric(5,31) + Metric(4,31)*Metric(5,32) - Metric(4,5)*Metric(31,32), Metric(4,32)*Metric(5,31) + Metric(4,31)*Metric(5,32) - Metric(4,5)*Metric(31,32)]
80 })
81Lorentz.append( { 'Name':'Lorentz45',
82 'Particles':[5, 1, 3, 3],
83 'Structure':[Metric(3,22)*Metric(4,21) + Metric(3,21)*Metric(4,22) - Metric(3,4)*Metric(21,22), Metric(3,22)*Metric(4,21) + Metric(3,21)*Metric(4,22) - Metric(3,4)*Metric(21,22), Metric(3,22)*Metric(4,21) + Metric(3,21)*Metric(4,22) - Metric(3,4)*Metric(21,22)]
84 })
85Lorentz.append( { 'Name':'Lorentz46',
86 'Particles':[5, 3, 3, 3, 3],
87 'Structure':[Metric(2,12)*Metric(3,11)*Metric(4,5) + Metric(2,11)*Metric(3,12)*Metric(4,5) - (Metric(2,12)*Metric(3,5)*Metric(4,11))/2. - (Metric(2,5)*Metric(3,12)*Metric(4,11))/2. - (Metric(2,11)*Metric(3,5)*Metric(4,12))/2. - (Metric(2,5)*Metric(3,11)*Metric(4,12))/2. - (Metric(2,12)*Metric(3,4)*Metric(5,11))/2. - (Metric(2,4)*Metric(3,12)*Metric(5,11))/2. + Metric(2,3)*Metric(4,12)*Metric(5,11) - (Metric(2,11)*Metric(3,4)*Metric(5,12))/2. - (Metric(2,4)*Metric(3,11)*Metric(5,12))/2. + Metric(2,3)*Metric(4,11)*Metric(5,12) + (Metric(2,5)*Metric(3,4)*Metric(11,12))/2. + (Metric(2,4)*Metric(3,5)*Metric(11,12))/2. - Metric(2,3)*Metric(4,5)*Metric(11,12)]
88 })
89#Lorentz.append( { 'Name':'Lorentz47',
90# 'Particles':[3, 2, 2, 5],
91# 'Structure':[-(Gamma(22,3,'a')*Metric(1,21)*ProjM('a',2))/2. - (Gamma(21,3,'a')*Metric(1,22)*ProjM('a',2))/2. + Gamma(1,3,'a')*Metric(21,22)*ProjM('a',2) - (Gamma(22,3,'a')*Metric(1,21)*ProjP('a',2))/2. - (Gamma(21,3,'a')*Metric(1,22)*ProjP('a',2))/2. + Gamma(1,3,'a')*Metric(21,22)*ProjP('a',2)]
92# })
93#Lorentz.append( { 'Name':'Lorentz49',
94# 'Particles':[3, 5, 2, 2],
95# 'Structure':[-(Gamma(22,4,'a')*Metric(1,21)*ProjM('a',3))/2. - (Gamma(21,4,'a')*Metric(1,22)*ProjM('a',3))/2. + Gamma(1,4,'a')*Metric(21,22)*ProjM('a',3) - (Gamma(22,4,'a')*Metric(1,21)*ProjP('a',3))/2. - (Gamma(21,4,'a')*Metric(1,22)*ProjP('a',3))/2. + Gamma(1,4,'a')*Metric(21,22)*ProjP('a',3)]
96# })
97#Lorentz.append( { 'Name':'Lorentz53',
98# 'Particles':[2, 5, 2, 3],
99# 'Structure':[-(Gamma(12,3,'a')*Metric(4,11)*ProjM('a',1)) - Gamma(11,3,'a')*Metric(4,12)*ProjM('a',1) + Gamma(4,3,'a')*Metric(11,12)*ProjM('a',1)]
100# })
101#Lorentz.append( { 'Name':'Lorentz54',
102# 'Particles':[5, 2, 2, 3],
103# 'Structure':[-(Gamma(12,3,'a')*Metric(4,11)*ProjM('a',2)) - Gamma(11,3,'a')*Metric(4,12)*ProjM('a',2) + Gamma(4,3,'a')*Metric(11,12)*ProjM('a',2)]
104# })