← Back to branch summary

~maddevelopers/mg5amcnlo/3.0.2-alpha0

« back to all changes in this revision

Viewing changes to HELAS/uttBxx.F

Added Template and HELAS into bzr

Show diffs side-by-side

added added

removed removed

Lines of Context:
HELAS/uttBxx.F
 
1
      subroutine uttbxx(wt1,wt2,gt,tmass,twidth,wt)
 
2
c-------------------CP3  2009.10-----------------
 
3
c
 
4
c This subroutine computes an off-shell tensor current from 
 
5
c the coupling of two non-propagating tensor bosons.
 
6
c
 
7
c input:
 
8
c       complex wt1(18)           : input non-propagating tensor                                                       t
 
9
c       complex wt2(18)         :  input non-propagating tensor                      T
 
10
c       complex gt             : coupling constant         gt=-1/Lambda
 
11
c output:
 
12
c       complex wt(18)         :  tensor 
 
13
 
 
14
      implicit none
 
15
      double complex wt(18), wt1(18), gt, wt2(18)
 
16
      double precision pT(4),pT2,tmass,twidth
 
17
      double complex T(6,4),T1(6,4),T2(6,4),d
 
18
      double precision MT(4,4)
 
19
      integer i, j,m,n
 
20
      double complex cZero
 
21
      double precision rZero, rTwo
 
22
      parameter( rZero = 0.0d0, rTwo = 2.0d0 )
 
23
      parameter( cZero = ( 0.0d0, 0.0d0 ) )
 
24
 
 
25
 
 
26
      T2(1,1) = wt2(1)
 
27
      T2(1,2) = wt2(2)
 
28
      T2(1,3) = wt2(3)
 
29
      T2(1,4) = wt2(4)
 
30
      T2(2,1) = wt2(5)
 
31
      T2(2,2) = wt2(6)
 
32
      T2(2,3) = wt2(7)
 
33
      T2(2,4) = wt2(8)
 
34
      T2(3,1) = wt2(9)
 
35
      T2(3,2) = wt2(10)
 
36
      T2(3,3) = wt2(11)
 
37
      T2(3,4) = wt2(12)
 
38
      T2(4,1) = wt2(13)
 
39
      T2(4,2) = wt2(14)
 
40
      T2(4,3) = wt2(15)
 
41
      T2(4,4) = wt2(16)
 
42
      T2(5,1) = wt2(17)
 
43
      T2(6,1) = wt2(18)
 
44
 
 
45
      T1(1,1) = wt1(1)
 
46
      T1(1,2) = wt1(2)
 
47
      T1(1,3) = wt1(3)
 
48
      T1(1,4) = wt1(4)
 
49
      T1(2,1) = wt1(5)
 
50
      T1(2,2) = wt1(6)
 
51
      T1(2,3) = wt1(7)
 
52
      T1(2,4) = wt1(8)
 
53
      T1(3,1) = wt1(9)
 
54
      T1(3,2) = wt1(10)
 
55
      T1(3,3) = wt1(11)
 
56
      T1(3,4) = wt1(12)
 
57
      T1(4,1) = wt1(13)
 
58
      T1(4,2) = wt1(14)
 
59
      T1(4,3) = wt1(15)
 
60
      T1(4,4) = wt1(16)
 
61
      T1(5,1) = wt1(17)
 
62
      T1(6,1) = wt1(18)
 
63
 
 
64
      T(5,1) = T2(5,1)+T1(5,1)
 
65
      T(6,1) = T2(6,1)+T1(6,1)
 
66
 
 
67
 
 
68
      do i=1,4
 
69
         do j=1,4
 
70
            MT(i,j) = 0.0d0
 
71
         enddo 
 
72
      enddo
 
73
      MT(1,1) =  1.0d0
 
74
      MT(2,2) = -1.0d0
 
75
      MT(3,3) = -1.0d0
 
76
      MT(4,4) = -1.0d0
 
77
      
 
78
    
 
79
        
 
80
         do m=1,4
 
81
           do n=1,4
 
82
 
 
83
          T(m,n)= -((2*MT(1,m)*MT(1,n)
 
84
     -   - (2*MT(m,n))/3.)*T1(1,2)*T2(1,2)) + 
 
85
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(1,2)*T2(1,2) - 
 
86
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(1,2)*T2(1,2) - 
 
87
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(1,2)*T2(1,2) + 
 
88
     -  2*(MT(2,n)*MT(3,m) + MT(2,m)*MT(3,n))*T1(1,3)*T2(1,2) + 
 
89
     -  2*(MT(2,n)*MT(4,m) + MT(2,m)*MT(4,n))*T1(1,4)*T2(1,2) + 
 
90
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(2,1)*T2(1,2) - 
 
91
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(2,1)*T2(1,2) + 
 
92
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(2,1)*T2(1,2) + 
 
93
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(2,1)*T2(1,2) - 
 
94
     -  2*(MT(1,n)*MT(3,m) + MT(1,m)*MT(3,n))*T1(2,3)*T2(1,2) - 
 
95
     -  2*(MT(1,n)*MT(4,m) + MT(1,m)*MT(4,n))*T1(2,4)*T2(1,2) - 
 
96
     -  2*(MT(2,n)*MT(3,m) + MT(2,m)*MT(3,n))*T1(3,1)*T2(1,2) + 
 
97
     -  2*(MT(1,n)*MT(3,m) + MT(1,m)*MT(3,n))*T1(3,2)*T2(1,2) - 
 
98
     -  2*(MT(2,n)*MT(4,m) + MT(2,m)*MT(4,n))*T1(4,1)*T2(1,2) + 
 
99
     -  2*(MT(1,n)*MT(4,m) + MT(1,m)*MT(4,n))*T1(4,2)*T2(1,2) + 
 
100
     -  2*(MT(2,n)*MT(3,m) + MT(2,m)*MT(3,n))*T1(1,2)*T2(1,3) - 
 
101
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(1,3)*T2(1,3) - 
 
102
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(1,3)*T2(1,3) + 
 
103
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(1,3)*T2(1,3) - 
 
104
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(1,3)*T2(1,3) + 
 
105
     -  2*(MT(3,n)*MT(4,m) + MT(3,m)*MT(4,n))*T1(1,4)*T2(1,3) - 
 
106
     -  2*(MT(2,n)*MT(3,m) + MT(2,m)*MT(3,n))*T1(2,1)*T2(1,3) + 
 
107
     -  2*(MT(1,n)*MT(2,m) + MT(1,m)*MT(2,n))*T1(2,3)*T2(1,3) + 
 
108
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(3,1)*T2(1,3) + 
 
109
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(3,1)*T2(1,3) - 
 
110
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(3,1)*T2(1,3) + 
 
111
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(3,1)*T2(1,3) - 
 
112
     -  2*(MT(1,n)*MT(2,m) + MT(1,m)*MT(2,n))*T1(3,2)*T2(1,3) - 
 
113
     -  2*(MT(1,n)*MT(4,m) + MT(1,m)*MT(4,n))*T1(3,4)*T2(1,3) - 
 
114
     -  2*(MT(3,n)*MT(4,m) + MT(3,m)*MT(4,n))*T1(4,1)*T2(1,3) + 
 
115
     -  2*(MT(1,n)*MT(4,m) + MT(1,m)*MT(4,n))*T1(4,3)*T2(1,3) + 
 
116
     -  2*(MT(2,n)*MT(4,m) + MT(2,m)*MT(4,n))*T1(1,2)*T2(1,4) + 
 
117
     -  2*(MT(3,n)*MT(4,m) + MT(3,m)*MT(4,n))*T1(1,3)*T2(1,4) - 
 
118
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(1,4)*T2(1,4) - 
 
119
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(1,4)*T2(1,4) - 
 
120
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(1,4)*T2(1,4) + 
 
121
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(1,4)*T2(1,4) - 
 
122
     -  2*(MT(2,n)*MT(4,m) + MT(2,m)*MT(4,n))*T1(2,1)*T2(1,4) + 
 
123
     -  2*(MT(1,n)*MT(2,m) + MT(1,m)*MT(2,n))*T1(2,4)*T2(1,4) - 
 
124
     -  2*(MT(3,n)*MT(4,m) + MT(3,m)*MT(4,n))*T1(3,1)*T2(1,4) + 
 
125
     -  2*(MT(1,n)*MT(3,m) + MT(1,m)*MT(3,n))*T1(3,4)*T2(1,4) + 
 
126
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(4,1)*T2(1,4) + 
 
127
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(4,1)*T2(1,4) + 
 
128
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(4,1)*T2(1,4) - 
 
129
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(4,1)*T2(1,4) - 
 
130
     -  2*(MT(1,n)*MT(2,m) + MT(1,m)*MT(2,n))*T1(4,2)*T2(1,4) - 
 
131
     -  2*(MT(1,n)*MT(3,m) + MT(1,m)*MT(3,n))*T1(4,3)*T2(1,4) + 
 
132
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(1,2)*T2(2,1) - 
 
133
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(1,2)*T2(2,1) + 
 
134
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(1,2)*T2(2,1) + 
 
135
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(1,2)*T2(2,1) - 
 
136
     -  2*(MT(2,n)*MT(3,m) + MT(2,m)*MT(3,n))*T1(1,3)*T2(2,1) - 
 
137
     -  2*(MT(2,n)*MT(4,m) + MT(2,m)*MT(4,n))*T1(1,4)*T2(2,1) - 
 
138
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(2,1)*T2(2,1) + 
 
139
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(2,1)*T2(2,1) - 
 
140
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(2,1)*T2(2,1) - 
 
141
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(2,1)*T2(2,1) + 
 
142
     -  2*(MT(1,n)*MT(3,m) + MT(1,m)*MT(3,n))*T1(2,3)*T2(2,1) + 
 
143
     -  2*(MT(1,n)*MT(4,m) + MT(1,m)*MT(4,n))*T1(2,4)*T2(2,1) + 
 
144
     -  2*(MT(2,n)*MT(3,m) + MT(2,m)*MT(3,n))*T1(3,1)*T2(2,1) - 
 
145
     -  2*(MT(1,n)*MT(3,m) + MT(1,m)*MT(3,n))*T1(3,2)*T2(2,1) + 
 
146
     -  2*(MT(2,n)*MT(4,m) + MT(2,m)*MT(4,n))*T1(4,1)*T2(2,1) - 
 
147
     -  2*(MT(1,n)*MT(4,m) + MT(1,m)*MT(4,n))*T1(4,2)*T2(2,1) - 
 
148
     -  2*(MT(1,n)*MT(3,m) + MT(1,m)*MT(3,n))*T1(1,2)*T2(2,3) + 
 
149
     -  2*(MT(1,n)*MT(2,m) + MT(1,m)*MT(2,n))*T1(1,3)*T2(2,3) + 
 
150
     -  2*(MT(1,n)*MT(3,m) + MT(1,m)*MT(3,n))*T1(2,1)*T2(2,3) - 
 
151
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(2,3)*T2(2,3) - 
 
152
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(2,3)*T2(2,3) - 
 
153
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(2,3)*T2(2,3) + 
 
154
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(2,3)*T2(2,3) - 
 
155
     -  2*(MT(3,n)*MT(4,m) + MT(3,m)*MT(4,n))*T1(2,4)*T2(2,3) - 
 
156
     -  2*(MT(1,n)*MT(2,m) + MT(1,m)*MT(2,n))*T1(3,1)*T2(2,3) + 
 
157
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(3,2)*T2(2,3) + 
 
158
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(3,2)*T2(2,3) + 
 
159
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(3,2)*T2(2,3) - 
 
160
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(3,2)*T2(2,3) + 
 
161
     -  2*(MT(2,n)*MT(4,m) + MT(2,m)*MT(4,n))*T1(3,4)*T2(2,3) + 
 
162
     -  2*(MT(3,n)*MT(4,m) + MT(3,m)*MT(4,n))*T1(4,2)*T2(2,3) - 
 
163
     -  2*(MT(2,n)*MT(4,m) + MT(2,m)*MT(4,n))*T1(4,3)*T2(2,3) - 
 
164
     -  2*(MT(1,n)*MT(4,m) + MT(1,m)*MT(4,n))*T1(1,2)*T2(2,4) + 
 
165
     -  2*(MT(1,n)*MT(2,m) + MT(1,m)*MT(2,n))*T1(1,4)*T2(2,4) + 
 
166
     -  2*(MT(1,n)*MT(4,m) + MT(1,m)*MT(4,n))*T1(2,1)*T2(2,4) - 
 
167
     -  2*(MT(3,n)*MT(4,m) + MT(3,m)*MT(4,n))*T1(2,3)*T2(2,4) - 
 
168
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(2,4)*T2(2,4) - 
 
169
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(2,4)*T2(2,4) + 
 
170
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(2,4)*T2(2,4) - 
 
171
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(2,4)*T2(2,4) + 
 
172
     -  2*(MT(3,n)*MT(4,m) + MT(3,m)*MT(4,n))*T1(3,2)*T2(2,4) - 
 
173
     -  2*(MT(2,n)*MT(3,m) + MT(2,m)*MT(3,n))*T1(3,4)*T2(2,4) - 
 
174
     -  2*(MT(1,n)*MT(2,m) + MT(1,m)*MT(2,n))*T1(4,1)*T2(2,4) + 
 
175
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(4,2)*T2(2,4) + 
 
176
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(4,2)*T2(2,4) - 
 
177
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(4,2)*T2(2,4) + 
 
178
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(4,2)*T2(2,4) + 
 
179
     -  2*(MT(2,n)*MT(3,m) + MT(2,m)*MT(3,n))*T1(4,3)*T2(2,4) - 
 
180
     -  2*(MT(2,n)*MT(3,m) + MT(2,m)*MT(3,n))*T1(1,2)*T2(3,1) + 
 
181
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(1,3)*T2(3,1) + 
 
182
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(1,3)*T2(3,1) - 
 
183
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(1,3)*T2(3,1) + 
 
184
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(1,3)*T2(3,1) - 
 
185
     -  2*(MT(3,n)*MT(4,m) + MT(3,m)*MT(4,n))*T1(1,4)*T2(3,1) + 
 
186
     -  2*(MT(2,n)*MT(3,m) + MT(2,m)*MT(3,n))*T1(2,1)*T2(3,1) - 
 
187
     -  2*(MT(1,n)*MT(2,m) + MT(1,m)*MT(2,n))*T1(2,3)*T2(3,1) - 
 
188
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(3,1)*T2(3,1) - 
 
189
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(3,1)*T2(3,1) + 
 
190
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(3,1)*T2(3,1) - 
 
191
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(3,1)*T2(3,1) + 
 
192
     -  2*(MT(1,n)*MT(2,m) + MT(1,m)*MT(2,n))*T1(3,2)*T2(3,1) + 
 
193
     -  2*(MT(1,n)*MT(4,m) + MT(1,m)*MT(4,n))*T1(3,4)*T2(3,1) + 
 
194
     -  2*(MT(3,n)*MT(4,m) + MT(3,m)*MT(4,n))*T1(4,1)*T2(3,1) - 
 
195
     -  2*(MT(1,n)*MT(4,m) + MT(1,m)*MT(4,n))*T1(4,3)*T2(3,1) + 
 
196
     -  2*(MT(1,n)*MT(3,m) + MT(1,m)*MT(3,n))*T1(1,2)*T2(3,2) - 
 
197
     -  2*(MT(1,n)*MT(2,m) + MT(1,m)*MT(2,n))*T1(1,3)*T2(3,2) - 
 
198
     -  2*(MT(1,n)*MT(3,m) + MT(1,m)*MT(3,n))*T1(2,1)*T2(3,2) + 
 
199
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(2,3)*T2(3,2) + 
 
200
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(2,3)*T2(3,2) + 
 
201
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(2,3)*T2(3,2) - 
 
202
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(2,3)*T2(3,2) + 
 
203
     -  2*(MT(3,n)*MT(4,m) + MT(3,m)*MT(4,n))*T1(2,4)*T2(3,2) + 
 
204
     -  2*(MT(1,n)*MT(2,m) + MT(1,m)*MT(2,n))*T1(3,1)*T2(3,2) - 
 
205
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(3,2)*T2(3,2) - 
 
206
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(3,2)*T2(3,2) - 
 
207
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(3,2)*T2(3,2) + 
 
208
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(3,2)*T2(3,2) - 
 
209
     -  2*(MT(2,n)*MT(4,m) + MT(2,m)*MT(4,n))*T1(3,4)*T2(3,2) - 
 
210
     -  2*(MT(3,n)*MT(4,m) + MT(3,m)*MT(4,n))*T1(4,2)*T2(3,2) + 
 
211
     -  2*(MT(2,n)*MT(4,m) + MT(2,m)*MT(4,n))*T1(4,3)*T2(3,2) - 
 
212
     -  2*(MT(1,n)*MT(4,m) + MT(1,m)*MT(4,n))*T1(1,3)*T2(3,4) + 
 
213
     -  2*(MT(1,n)*MT(3,m) + MT(1,m)*MT(3,n))*T1(1,4)*T2(3,4) + 
 
214
     -  2*(MT(2,n)*MT(4,m) + MT(2,m)*MT(4,n))*T1(2,3)*T2(3,4) - 
 
215
     -  2*(MT(2,n)*MT(3,m) + MT(2,m)*MT(3,n))*T1(2,4)*T2(3,4) + 
 
216
     -  2*(MT(1,n)*MT(4,m) + MT(1,m)*MT(4,n))*T1(3,1)*T2(3,4) - 
 
217
     -  2*(MT(2,n)*MT(4,m) + MT(2,m)*MT(4,n))*T1(3,2)*T2(3,4) - 
 
218
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(3,4)*T2(3,4) + 
 
219
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(3,4)*T2(3,4) - 
 
220
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(3,4)*T2(3,4) - 
 
221
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(3,4)*T2(3,4) - 
 
222
     -  2*(MT(1,n)*MT(3,m) + MT(1,m)*MT(3,n))*T1(4,1)*T2(3,4) + 
 
223
     -  2*(MT(2,n)*MT(3,m) + MT(2,m)*MT(3,n))*T1(4,2)*T2(3,4) + 
 
224
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(4,3)*T2(3,4) - 
 
225
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(4,3)*T2(3,4) + 
 
226
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(4,3)*T2(3,4) + 
 
227
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(4,3)*T2(3,4) - 
 
228
     -  2*(MT(2,n)*MT(4,m) + MT(2,m)*MT(4,n))*T1(1,2)*T2(4,1) - 
 
229
     -  2*(MT(3,n)*MT(4,m) + MT(3,m)*MT(4,n))*T1(1,3)*T2(4,1) + 
 
230
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(1,4)*T2(4,1) + 
 
231
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(1,4)*T2(4,1) + 
 
232
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(1,4)*T2(4,1) - 
 
233
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(1,4)*T2(4,1) + 
 
234
     -  2*(MT(2,n)*MT(4,m) + MT(2,m)*MT(4,n))*T1(2,1)*T2(4,1) - 
 
235
     -  2*(MT(1,n)*MT(2,m) + MT(1,m)*MT(2,n))*T1(2,4)*T2(4,1) + 
 
236
     -  2*(MT(3,n)*MT(4,m) + MT(3,m)*MT(4,n))*T1(3,1)*T2(4,1) - 
 
237
     -  2*(MT(1,n)*MT(3,m) + MT(1,m)*MT(3,n))*T1(3,4)*T2(4,1) - 
 
238
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(4,1)*T2(4,1) - 
 
239
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(4,1)*T2(4,1) - 
 
240
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(4,1)*T2(4,1) + 
 
241
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(4,1)*T2(4,1) + 
 
242
     -  2*(MT(1,n)*MT(2,m) + MT(1,m)*MT(2,n))*T1(4,2)*T2(4,1) + 
 
243
     -  2*(MT(1,n)*MT(3,m) + MT(1,m)*MT(3,n))*T1(4,3)*T2(4,1) + 
 
244
     -  2*(MT(1,n)*MT(4,m) + MT(1,m)*MT(4,n))*T1(1,2)*T2(4,2) - 
 
245
     -  2*(MT(1,n)*MT(2,m) + MT(1,m)*MT(2,n))*T1(1,4)*T2(4,2) - 
 
246
     -  2*(MT(1,n)*MT(4,m) + MT(1,m)*MT(4,n))*T1(2,1)*T2(4,2) + 
 
247
     -  2*(MT(3,n)*MT(4,m) + MT(3,m)*MT(4,n))*T1(2,3)*T2(4,2) + 
 
248
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(2,4)*T2(4,2) + 
 
249
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(2,4)*T2(4,2) - 
 
250
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(2,4)*T2(4,2) + 
 
251
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(2,4)*T2(4,2) - 
 
252
     -  2*(MT(3,n)*MT(4,m) + MT(3,m)*MT(4,n))*T1(3,2)*T2(4,2) + 
 
253
     -  2*(MT(2,n)*MT(3,m) + MT(2,m)*MT(3,n))*T1(3,4)*T2(4,2) + 
 
254
     -  2*(MT(1,n)*MT(2,m) + MT(1,m)*MT(2,n))*T1(4,1)*T2(4,2) - 
 
255
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(4,2)*T2(4,2) - 
 
256
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(4,2)*T2(4,2) + 
 
257
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(4,2)*T2(4,2) - 
 
258
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(4,2)*T2(4,2) - 
 
259
     -  2*(MT(2,n)*MT(3,m) + MT(2,m)*MT(3,n))*T1(4,3)*T2(4,2) + 
 
260
     -  2*(MT(1,n)*MT(4,m) + MT(1,m)*MT(4,n))*T1(1,3)*T2(4,3) - 
 
261
     -  2*(MT(1,n)*MT(3,m) + MT(1,m)*MT(3,n))*T1(1,4)*T2(4,3) - 
 
262
     -  2*(MT(2,n)*MT(4,m) + MT(2,m)*MT(4,n))*T1(2,3)*T2(4,3) + 
 
263
     -  2*(MT(2,n)*MT(3,m) + MT(2,m)*MT(3,n))*T1(2,4)*T2(4,3) - 
 
264
     -  2*(MT(1,n)*MT(4,m) + MT(1,m)*MT(4,n))*T1(3,1)*T2(4,3) + 
 
265
     -  2*(MT(2,n)*MT(4,m) + MT(2,m)*MT(4,n))*T1(3,2)*T2(4,3) + 
 
266
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(3,4)*T2(4,3) - 
 
267
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(3,4)*T2(4,3) + 
 
268
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(3,4)*T2(4,3) + 
 
269
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(3,4)*T2(4,3) + 
 
270
     -  2*(MT(1,n)*MT(3,m) + MT(1,m)*MT(3,n))*T1(4,1)*T2(4,3) - 
 
271
     -  2*(MT(2,n)*MT(3,m) + MT(2,m)*MT(3,n))*T1(4,2)*T2(4,3) - 
 
272
     -  (2*MT(1,m)*MT(1,n) - (2*MT(m,n))/3.)*T1(4,3)*T2(4,3) + 
 
273
     -  (2*MT(2,m)*MT(2,n) + (2*MT(m,n))/3.)*T1(4,3)*T2(4,3) - 
 
274
     -  (2*MT(3,m)*MT(3,n) + (2*MT(m,n))/3.)*T1(4,3)*T2(4,3) - 
 
275
     -  (2*MT(4,m)*MT(4,n) + (2*MT(m,n))/3.)*T1(4,3)*T2(4,3)
 
276
           enddo
 
277
         enddo
 
278
       
 
279
       pT(1) = dreal(T(5,1))
 
280
       pT(2) = dreal(T(6,1))
 
281
       pT(3) = dimag(T(6,1))
 
282
       pT(4) = dimag(T(5,1))
 
283
       pT2 = pT(1)**2-pT(2)**2-pT(3)**2-pT(4)**2
 
284
 
 
285
 
 
286
         do m=1,4
 
287
           do n=1,4
 
288
          T(m,n)=T(m,n)
 
289
     - +((-2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) + 
 
290
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(1,2)*
 
291
     -   T2(1,2) + (2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
292
     -        pT(2)*pT(m)*MT(2,n)) - 
 
293
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(1,2)*T2(1,2)
 
294
     -   + (-2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) + 
 
295
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(1,2)*T2(1,2)
 
296
     -   + (-2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) + 
 
297
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(1,2)*T2(1,2)
 
298
     -   + 2*(-(pT(3)*pT(n)*MT(2,m)) - pT(3)*pT(m)*MT(2,n) - 
 
299
     -     pT(2)*pT(n)*MT(3,m) - pT(2)*pT(m)*MT(3,n) + 
 
300
     -     (2*pT(2)*pT(3)*MT(m,n))/3.)*T1(1,3)*T2(1,2) + 
 
301
     -  2*(-(pT(4)*pT(n)*MT(2,m)) - pT(4)*pT(m)*MT(2,n) - 
 
302
     -     pT(2)*pT(n)*MT(4,m) - pT(2)*pT(m)*MT(4,n) + 
 
303
     -     (2*pT(2)*pT(4)*MT(m,n))/3.)*T1(1,4)*T2(1,2) + 
 
304
     -  (2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) - 
 
305
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(2,1)*
 
306
     -   T2(1,2) + (-2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
307
     -        pT(2)*pT(m)*MT(2,n)) + 
 
308
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(2,1)*T2(1,2)
 
309
     -   + (2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) - 
 
310
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(2,1)*T2(1,2)
 
311
     -   + (2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) - 
 
312
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(2,1)*T2(1,2)
 
313
     -   - 2*(-(pT(3)*pT(n)*MT(1,m)) - pT(3)*pT(m)*MT(1,n) - 
 
314
     -     pT(1)*pT(n)*MT(3,m) - pT(1)*pT(m)*MT(3,n) + 
 
315
     -     (2*pT(1)*pT(3)*MT(m,n))/3.)*T1(2,3)*T2(1,2) - 
 
316
     -  2*(-(pT(4)*pT(n)*MT(1,m)) - pT(4)*pT(m)*MT(1,n) - 
 
317
     -     pT(1)*pT(n)*MT(4,m) - pT(1)*pT(m)*MT(4,n) + 
 
318
     -     (2*pT(1)*pT(4)*MT(m,n))/3.)*T1(2,4)*T2(1,2) - 
 
319
     -  2*(-(pT(3)*pT(n)*MT(2,m)) - pT(3)*pT(m)*MT(2,n) - 
 
320
     -     pT(2)*pT(n)*MT(3,m) - pT(2)*pT(m)*MT(3,n) + 
 
321
     -     (2*pT(2)*pT(3)*MT(m,n))/3.)*T1(3,1)*T2(1,2) + 
 
322
     -  2*(-(pT(3)*pT(n)*MT(1,m)) - pT(3)*pT(m)*MT(1,n) - 
 
323
     -     pT(1)*pT(n)*MT(3,m) - pT(1)*pT(m)*MT(3,n) + 
 
324
     -     (2*pT(1)*pT(3)*MT(m,n))/3.)*T1(3,2)*T2(1,2) - 
 
325
     -  2*(-(pT(4)*pT(n)*MT(2,m)) - pT(4)*pT(m)*MT(2,n) - 
 
326
     -     pT(2)*pT(n)*MT(4,m) - pT(2)*pT(m)*MT(4,n) + 
 
327
     -     (2*pT(2)*pT(4)*MT(m,n))/3.)*T1(4,1)*T2(1,2) + 
 
328
     -  2*(-(pT(4)*pT(n)*MT(1,m)) - pT(4)*pT(m)*MT(1,n) - 
 
329
     -     pT(1)*pT(n)*MT(4,m) - pT(1)*pT(m)*MT(4,n) + 
 
330
     -     (2*pT(1)*pT(4)*MT(m,n))/3.)*T1(4,2)*T2(1,2) + 
 
331
     -  2*(-(pT(3)*pT(n)*MT(2,m)) - pT(3)*pT(m)*MT(2,n) - 
 
332
     -     pT(2)*pT(n)*MT(3,m) - pT(2)*pT(m)*MT(3,n) + 
 
333
     -     (2*pT(2)*pT(3)*MT(m,n))/3.)*T1(1,2)*T2(1,3) + 
 
334
     -  (-2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) + 
 
335
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(1,3)*
 
336
     -   T2(1,3) + (-2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
337
     -        pT(2)*pT(m)*MT(2,n)) + 
 
338
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(1,3)*T2(1,3)
 
339
     -   + (2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) - 
 
340
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(1,3)*T2(1,3)
 
341
     -   + (-2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) + 
 
342
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(1,3)*T2(1,3)
 
343
     -   + 2*(-(pT(4)*pT(n)*MT(3,m)) - pT(4)*pT(m)*MT(3,n) - 
 
344
     -     pT(3)*pT(n)*MT(4,m) - pT(3)*pT(m)*MT(4,n) + 
 
345
     -     (2*pT(3)*pT(4)*MT(m,n))/3.)*T1(1,4)*T2(1,3) - 
 
346
     -  2*(-(pT(3)*pT(n)*MT(2,m)) - pT(3)*pT(m)*MT(2,n) - 
 
347
     -     pT(2)*pT(n)*MT(3,m) - pT(2)*pT(m)*MT(3,n) + 
 
348
     -     (2*pT(2)*pT(3)*MT(m,n))/3.)*T1(2,1)*T2(1,3) + 
 
349
     -  2*(-(pT(2)*pT(n)*MT(1,m)) - pT(2)*pT(m)*MT(1,n) - 
 
350
     -     pT(1)*pT(n)*MT(2,m) - pT(1)*pT(m)*MT(2,n) + 
 
351
     -     (2*pT(1)*pT(2)*MT(m,n))/3.)*T1(2,3)*T2(1,3) + 
 
352
     -  (2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) - 
 
353
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(3,1)*
 
354
     -   T2(1,3) + (2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
355
     -        pT(2)*pT(m)*MT(2,n)) - 
 
356
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(3,1)*T2(1,3)
 
357
     -   + (-2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) + 
 
358
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(3,1)*T2(1,3)
 
359
     -   + (2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) - 
 
360
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(3,1)*T2(1,3)
 
361
     -   - 2*(-(pT(2)*pT(n)*MT(1,m)) - pT(2)*pT(m)*MT(1,n) - 
 
362
     -     pT(1)*pT(n)*MT(2,m) - pT(1)*pT(m)*MT(2,n) + 
 
363
     -     (2*pT(1)*pT(2)*MT(m,n))/3.)*T1(3,2)*T2(1,3) - 
 
364
     -  2*(-(pT(4)*pT(n)*MT(1,m)) - pT(4)*pT(m)*MT(1,n) - 
 
365
     -     pT(1)*pT(n)*MT(4,m) - pT(1)*pT(m)*MT(4,n) + 
 
366
     -     (2*pT(1)*pT(4)*MT(m,n))/3.)*T1(3,4)*T2(1,3) - 
 
367
     -  2*(-(pT(4)*pT(n)*MT(3,m)) - pT(4)*pT(m)*MT(3,n) - 
 
368
     -     pT(3)*pT(n)*MT(4,m) - pT(3)*pT(m)*MT(4,n) + 
 
369
     -     (2*pT(3)*pT(4)*MT(m,n))/3.)*T1(4,1)*T2(1,3) + 
 
370
     -  2*(-(pT(4)*pT(n)*MT(1,m)) - pT(4)*pT(m)*MT(1,n) - 
 
371
     -     pT(1)*pT(n)*MT(4,m) - pT(1)*pT(m)*MT(4,n) + 
 
372
     -     (2*pT(1)*pT(4)*MT(m,n))/3.)*T1(4,3)*T2(1,3) + 
 
373
     -  2*(-(pT(4)*pT(n)*MT(2,m)) - pT(4)*pT(m)*MT(2,n) - 
 
374
     -     pT(2)*pT(n)*MT(4,m) - pT(2)*pT(m)*MT(4,n) + 
 
375
     -     (2*pT(2)*pT(4)*MT(m,n))/3.)*T1(1,2)*T2(1,4) + 
 
376
     -  2*(-(pT(4)*pT(n)*MT(3,m)) - pT(4)*pT(m)*MT(3,n) - 
 
377
     -     pT(3)*pT(n)*MT(4,m) - pT(3)*pT(m)*MT(4,n) + 
 
378
     -     (2*pT(3)*pT(4)*MT(m,n))/3.)*T1(1,3)*T2(1,4) + 
 
379
     -  (-2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) + 
 
380
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(1,4)*
 
381
     -   T2(1,4) + (-2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
382
     -        pT(2)*pT(m)*MT(2,n)) + 
 
383
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(1,4)*T2(1,4)
 
384
     -   + (-2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) + 
 
385
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(1,4)*T2(1,4)
 
386
     -   + (2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) - 
 
387
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(1,4)*T2(1,4)
 
388
     -   - 2*(-(pT(4)*pT(n)*MT(2,m)) - pT(4)*pT(m)*MT(2,n) - 
 
389
     -     pT(2)*pT(n)*MT(4,m) - pT(2)*pT(m)*MT(4,n) + 
 
390
     -     (2*pT(2)*pT(4)*MT(m,n))/3.)*T1(2,1)*T2(1,4) + 
 
391
     -  2*(-(pT(2)*pT(n)*MT(1,m)) - pT(2)*pT(m)*MT(1,n) - 
 
392
     -     pT(1)*pT(n)*MT(2,m) - pT(1)*pT(m)*MT(2,n) + 
 
393
     -     (2*pT(1)*pT(2)*MT(m,n))/3.)*T1(2,4)*T2(1,4) - 
 
394
     -  2*(-(pT(4)*pT(n)*MT(3,m)) - pT(4)*pT(m)*MT(3,n) - 
 
395
     -     pT(3)*pT(n)*MT(4,m) - pT(3)*pT(m)*MT(4,n) + 
 
396
     -     (2*pT(3)*pT(4)*MT(m,n))/3.)*T1(3,1)*T2(1,4) + 
 
397
     -  2*(-(pT(3)*pT(n)*MT(1,m)) - pT(3)*pT(m)*MT(1,n) - 
 
398
     -     pT(1)*pT(n)*MT(3,m) - pT(1)*pT(m)*MT(3,n) + 
 
399
     -     (2*pT(1)*pT(3)*MT(m,n))/3.)*T1(3,4)*T2(1,4) + 
 
400
     -  (2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) - 
 
401
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(4,1)*
 
402
     -   T2(1,4) + (2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
403
     -        pT(2)*pT(m)*MT(2,n)) - 
 
404
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(4,1)*T2(1,4)
 
405
     -   + (2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) - 
 
406
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(4,1)*T2(1,4)
 
407
     -   + (-2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) + 
 
408
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(4,1)*T2(1,4)
 
409
     -   - 2*(-(pT(2)*pT(n)*MT(1,m)) - pT(2)*pT(m)*MT(1,n) - 
 
410
     -     pT(1)*pT(n)*MT(2,m) - pT(1)*pT(m)*MT(2,n) + 
 
411
     -     (2*pT(1)*pT(2)*MT(m,n))/3.)*T1(4,2)*T2(1,4) - 
 
412
     -  2*(-(pT(3)*pT(n)*MT(1,m)) - pT(3)*pT(m)*MT(1,n) - 
 
413
     -     pT(1)*pT(n)*MT(3,m) - pT(1)*pT(m)*MT(3,n) + 
 
414
     -     (2*pT(1)*pT(3)*MT(m,n))/3.)*T1(4,3)*T2(1,4) + 
 
415
     -  (2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) - 
 
416
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(1,2)*
 
417
     -   T2(2,1) + (-2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
418
     -        pT(2)*pT(m)*MT(2,n)) + 
 
419
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(1,2)*T2(2,1)
 
420
     -   + (2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) - 
 
421
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(1,2)*T2(2,1)
 
422
     -   + (2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) - 
 
423
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(1,2)*T2(2,1)
 
424
     -   - 2*(-(pT(3)*pT(n)*MT(2,m)) - pT(3)*pT(m)*MT(2,n) - 
 
425
     -     pT(2)*pT(n)*MT(3,m) - pT(2)*pT(m)*MT(3,n) + 
 
426
     -     (2*pT(2)*pT(3)*MT(m,n))/3.)*T1(1,3)*T2(2,1) - 
 
427
     -  2*(-(pT(4)*pT(n)*MT(2,m)) - pT(4)*pT(m)*MT(2,n) - 
 
428
     -     pT(2)*pT(n)*MT(4,m) - pT(2)*pT(m)*MT(4,n) + 
 
429
     -     (2*pT(2)*pT(4)*MT(m,n))/3.)*T1(1,4)*T2(2,1) + 
 
430
     -  (-2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) + 
 
431
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(2,1)*
 
432
     -   T2(2,1) + (2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
433
     -        pT(2)*pT(m)*MT(2,n)) - 
 
434
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(2,1)*T2(2,1)
 
435
     -   + (-2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) + 
 
436
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(2,1)*T2(2,1)
 
437
     -   + (-2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) + 
 
438
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(2,1)*T2(2,1)
 
439
     -   + 2*(-(pT(3)*pT(n)*MT(1,m)) - pT(3)*pT(m)*MT(1,n) - 
 
440
     -     pT(1)*pT(n)*MT(3,m) - pT(1)*pT(m)*MT(3,n) + 
 
441
     -     (2*pT(1)*pT(3)*MT(m,n))/3.)*T1(2,3)*T2(2,1) + 
 
442
     -  2*(-(pT(4)*pT(n)*MT(1,m)) - pT(4)*pT(m)*MT(1,n) - 
 
443
     -     pT(1)*pT(n)*MT(4,m) - pT(1)*pT(m)*MT(4,n) + 
 
444
     -     (2*pT(1)*pT(4)*MT(m,n))/3.)*T1(2,4)*T2(2,1) + 
 
445
     -  2*(-(pT(3)*pT(n)*MT(2,m)) - pT(3)*pT(m)*MT(2,n) - 
 
446
     -     pT(2)*pT(n)*MT(3,m) - pT(2)*pT(m)*MT(3,n) + 
 
447
     -     (2*pT(2)*pT(3)*MT(m,n))/3.)*T1(3,1)*T2(2,1) - 
 
448
     -  2*(-(pT(3)*pT(n)*MT(1,m)) - pT(3)*pT(m)*MT(1,n) - 
 
449
     -     pT(1)*pT(n)*MT(3,m) - pT(1)*pT(m)*MT(3,n) + 
 
450
     -     (2*pT(1)*pT(3)*MT(m,n))/3.)*T1(3,2)*T2(2,1) + 
 
451
     -  2*(-(pT(4)*pT(n)*MT(2,m)) - pT(4)*pT(m)*MT(2,n) - 
 
452
     -     pT(2)*pT(n)*MT(4,m) - pT(2)*pT(m)*MT(4,n) + 
 
453
     -     (2*pT(2)*pT(4)*MT(m,n))/3.)*T1(4,1)*T2(2,1) - 
 
454
     -  2*(-(pT(4)*pT(n)*MT(1,m)) - pT(4)*pT(m)*MT(1,n) - 
 
455
     -     pT(1)*pT(n)*MT(4,m) - pT(1)*pT(m)*MT(4,n) + 
 
456
     -     (2*pT(1)*pT(4)*MT(m,n))/3.)*T1(4,2)*T2(2,1) - 
 
457
     -  2*(-(pT(3)*pT(n)*MT(1,m)) - pT(3)*pT(m)*MT(1,n) - 
 
458
     -     pT(1)*pT(n)*MT(3,m) - pT(1)*pT(m)*MT(3,n) + 
 
459
     -     (2*pT(1)*pT(3)*MT(m,n))/3.)*T1(1,2)*T2(2,3) + 
 
460
     -  2*(-(pT(2)*pT(n)*MT(1,m)) - pT(2)*pT(m)*MT(1,n) - 
 
461
     -     pT(1)*pT(n)*MT(2,m) - pT(1)*pT(m)*MT(2,n) + 
 
462
     -     (2*pT(1)*pT(2)*MT(m,n))/3.)*T1(1,3)*T2(2,3) + 
 
463
     -  2*(-(pT(3)*pT(n)*MT(1,m)) - pT(3)*pT(m)*MT(1,n) - 
 
464
     -     pT(1)*pT(n)*MT(3,m) - pT(1)*pT(m)*MT(3,n) + 
 
465
     -     (2*pT(1)*pT(3)*MT(m,n))/3.)*T1(2,1)*T2(2,3) + 
 
466
     -  (-2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) + 
 
467
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(2,3)*
 
468
     -   T2(2,3) + (-2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
469
     -        pT(2)*pT(m)*MT(2,n)) + 
 
470
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(2,3)*T2(2,3)
 
471
     -   + (-2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) + 
 
472
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(2,3)*T2(2,3)
 
473
     -   + (2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) - 
 
474
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(2,3)*T2(2,3)
 
475
     -   - 2*(-(pT(4)*pT(n)*MT(3,m)) - pT(4)*pT(m)*MT(3,n) - 
 
476
     -     pT(3)*pT(n)*MT(4,m) - pT(3)*pT(m)*MT(4,n) + 
 
477
     -     (2*pT(3)*pT(4)*MT(m,n))/3.)*T1(2,4)*T2(2,3) - 
 
478
     -  2*(-(pT(2)*pT(n)*MT(1,m)) - pT(2)*pT(m)*MT(1,n) - 
 
479
     -     pT(1)*pT(n)*MT(2,m) - pT(1)*pT(m)*MT(2,n) + 
 
480
     -     (2*pT(1)*pT(2)*MT(m,n))/3.)*T1(3,1)*T2(2,3) + 
 
481
     -  (2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) - 
 
482
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(3,2)*
 
483
     -   T2(2,3) + (2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
484
     -        pT(2)*pT(m)*MT(2,n)) - 
 
485
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(3,2)*T2(2,3)
 
486
     -   + (2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) - 
 
487
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(3,2)*T2(2,3)
 
488
     -   + (-2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) + 
 
489
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(3,2)*T2(2,3)
 
490
     -   + 2*(-(pT(4)*pT(n)*MT(2,m)) - pT(4)*pT(m)*MT(2,n) - 
 
491
     -     pT(2)*pT(n)*MT(4,m) - pT(2)*pT(m)*MT(4,n) + 
 
492
     -     (2*pT(2)*pT(4)*MT(m,n))/3.)*T1(3,4)*T2(2,3) + 
 
493
     -  2*(-(pT(4)*pT(n)*MT(3,m)) - pT(4)*pT(m)*MT(3,n) - 
 
494
     -     pT(3)*pT(n)*MT(4,m) - pT(3)*pT(m)*MT(4,n) + 
 
495
     -     (2*pT(3)*pT(4)*MT(m,n))/3.)*T1(4,2)*T2(2,3) - 
 
496
     -  2*(-(pT(4)*pT(n)*MT(2,m)) - pT(4)*pT(m)*MT(2,n) - 
 
497
     -     pT(2)*pT(n)*MT(4,m) - pT(2)*pT(m)*MT(4,n) + 
 
498
     -     (2*pT(2)*pT(4)*MT(m,n))/3.)*T1(4,3)*T2(2,3) - 
 
499
     -  2*(-(pT(4)*pT(n)*MT(1,m)) - pT(4)*pT(m)*MT(1,n) - 
 
500
     -     pT(1)*pT(n)*MT(4,m) - pT(1)*pT(m)*MT(4,n) + 
 
501
     -     (2*pT(1)*pT(4)*MT(m,n))/3.)*T1(1,2)*T2(2,4) + 
 
502
     -  2*(-(pT(2)*pT(n)*MT(1,m)) - pT(2)*pT(m)*MT(1,n) - 
 
503
     -     pT(1)*pT(n)*MT(2,m) - pT(1)*pT(m)*MT(2,n) + 
 
504
     -     (2*pT(1)*pT(2)*MT(m,n))/3.)*T1(1,4)*T2(2,4) + 
 
505
     -  2*(-(pT(4)*pT(n)*MT(1,m)) - pT(4)*pT(m)*MT(1,n) - 
 
506
     -     pT(1)*pT(n)*MT(4,m) - pT(1)*pT(m)*MT(4,n) + 
 
507
     -     (2*pT(1)*pT(4)*MT(m,n))/3.)*T1(2,1)*T2(2,4) - 
 
508
     -  2*(-(pT(4)*pT(n)*MT(3,m)) - pT(4)*pT(m)*MT(3,n) - 
 
509
     -     pT(3)*pT(n)*MT(4,m) - pT(3)*pT(m)*MT(4,n) + 
 
510
     -     (2*pT(3)*pT(4)*MT(m,n))/3.)*T1(2,3)*T2(2,4) + 
 
511
     -  (-2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) + 
 
512
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(2,4)*
 
513
     -   T2(2,4) + (-2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
514
     -        pT(2)*pT(m)*MT(2,n)) + 
 
515
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(2,4)*T2(2,4)
 
516
     -   + (2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) - 
 
517
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(2,4)*T2(2,4)
 
518
     -   + (-2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) + 
 
519
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(2,4)*T2(2,4)
 
520
     -   + 2*(-(pT(4)*pT(n)*MT(3,m)) - pT(4)*pT(m)*MT(3,n) - 
 
521
     -     pT(3)*pT(n)*MT(4,m) - pT(3)*pT(m)*MT(4,n) + 
 
522
     -     (2*pT(3)*pT(4)*MT(m,n))/3.)*T1(3,2)*T2(2,4) - 
 
523
     -  2*(-(pT(3)*pT(n)*MT(2,m)) - pT(3)*pT(m)*MT(2,n) - 
 
524
     -     pT(2)*pT(n)*MT(3,m) - pT(2)*pT(m)*MT(3,n) + 
 
525
     -     (2*pT(2)*pT(3)*MT(m,n))/3.)*T1(3,4)*T2(2,4) - 
 
526
     -  2*(-(pT(2)*pT(n)*MT(1,m)) - pT(2)*pT(m)*MT(1,n) - 
 
527
     -     pT(1)*pT(n)*MT(2,m) - pT(1)*pT(m)*MT(2,n) + 
 
528
     -     (2*pT(1)*pT(2)*MT(m,n))/3.)*T1(4,1)*T2(2,4) + 
 
529
     -  (2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) - 
 
530
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(4,2)*
 
531
     -   T2(2,4) + (2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
532
     -        pT(2)*pT(m)*MT(2,n)) - 
 
533
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(4,2)*T2(2,4)
 
534
     -   + (-2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) + 
 
535
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(4,2)*T2(2,4)
 
536
     -   + (2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) - 
 
537
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(4,2)*T2(2,4)
 
538
     -   + 2*(-(pT(3)*pT(n)*MT(2,m)) - pT(3)*pT(m)*MT(2,n) - 
 
539
     -     pT(2)*pT(n)*MT(3,m) - pT(2)*pT(m)*MT(3,n) + 
 
540
     -     (2*pT(2)*pT(3)*MT(m,n))/3.)*T1(4,3)*T2(2,4) - 
 
541
     -  2*(-(pT(3)*pT(n)*MT(2,m)) - pT(3)*pT(m)*MT(2,n) - 
 
542
     -     pT(2)*pT(n)*MT(3,m) - pT(2)*pT(m)*MT(3,n) + 
 
543
     -     (2*pT(2)*pT(3)*MT(m,n))/3.)*T1(1,2)*T2(3,1) + 
 
544
     -  (2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) - 
 
545
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(1,3)*
 
546
     -   T2(3,1) + (2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
547
     -        pT(2)*pT(m)*MT(2,n)) - 
 
548
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(1,3)*T2(3,1)
 
549
     -   + (-2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) + 
 
550
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(1,3)*T2(3,1)
 
551
     -   + (2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) - 
 
552
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(1,3)*T2(3,1)
 
553
     -   - 2*(-(pT(4)*pT(n)*MT(3,m)) - pT(4)*pT(m)*MT(3,n) - 
 
554
     -     pT(3)*pT(n)*MT(4,m) - pT(3)*pT(m)*MT(4,n) + 
 
555
     -     (2*pT(3)*pT(4)*MT(m,n))/3.)*T1(1,4)*T2(3,1) + 
 
556
     -  2*(-(pT(3)*pT(n)*MT(2,m)) - pT(3)*pT(m)*MT(2,n) - 
 
557
     -     pT(2)*pT(n)*MT(3,m) - pT(2)*pT(m)*MT(3,n) + 
 
558
     -     (2*pT(2)*pT(3)*MT(m,n))/3.)*T1(2,1)*T2(3,1) - 
 
559
     -  2*(-(pT(2)*pT(n)*MT(1,m)) - pT(2)*pT(m)*MT(1,n) - 
 
560
     -     pT(1)*pT(n)*MT(2,m) - pT(1)*pT(m)*MT(2,n) + 
 
561
     -     (2*pT(1)*pT(2)*MT(m,n))/3.)*T1(2,3)*T2(3,1) + 
 
562
     -  (-2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) + 
 
563
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(3,1)*
 
564
     -   T2(3,1) + (-2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
565
     -        pT(2)*pT(m)*MT(2,n)) + 
 
566
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(3,1)*T2(3,1)
 
567
     -   + (2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) - 
 
568
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(3,1)*T2(3,1)
 
569
     -   + (-2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) + 
 
570
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(3,1)*T2(3,1)
 
571
     -   + 2*(-(pT(2)*pT(n)*MT(1,m)) - pT(2)*pT(m)*MT(1,n) - 
 
572
     -     pT(1)*pT(n)*MT(2,m) - pT(1)*pT(m)*MT(2,n) + 
 
573
     -     (2*pT(1)*pT(2)*MT(m,n))/3.)*T1(3,2)*T2(3,1) + 
 
574
     -  2*(-(pT(4)*pT(n)*MT(1,m)) - pT(4)*pT(m)*MT(1,n) - 
 
575
     -     pT(1)*pT(n)*MT(4,m) - pT(1)*pT(m)*MT(4,n) + 
 
576
     -     (2*pT(1)*pT(4)*MT(m,n))/3.)*T1(3,4)*T2(3,1) + 
 
577
     -  2*(-(pT(4)*pT(n)*MT(3,m)) - pT(4)*pT(m)*MT(3,n) - 
 
578
     -     pT(3)*pT(n)*MT(4,m) - pT(3)*pT(m)*MT(4,n) + 
 
579
     -     (2*pT(3)*pT(4)*MT(m,n))/3.)*T1(4,1)*T2(3,1) - 
 
580
     -  2*(-(pT(4)*pT(n)*MT(1,m)) - pT(4)*pT(m)*MT(1,n) - 
 
581
     -     pT(1)*pT(n)*MT(4,m) - pT(1)*pT(m)*MT(4,n) + 
 
582
     -     (2*pT(1)*pT(4)*MT(m,n))/3.)*T1(4,3)*T2(3,1) + 
 
583
     -  2*(-(pT(3)*pT(n)*MT(1,m)) - pT(3)*pT(m)*MT(1,n) - 
 
584
     -     pT(1)*pT(n)*MT(3,m) - pT(1)*pT(m)*MT(3,n) + 
 
585
     -     (2*pT(1)*pT(3)*MT(m,n))/3.)*T1(1,2)*T2(3,2) - 
 
586
     -  2*(-(pT(2)*pT(n)*MT(1,m)) - pT(2)*pT(m)*MT(1,n) - 
 
587
     -     pT(1)*pT(n)*MT(2,m) - pT(1)*pT(m)*MT(2,n) + 
 
588
     -     (2*pT(1)*pT(2)*MT(m,n))/3.)*T1(1,3)*T2(3,2) - 
 
589
     -  2*(-(pT(3)*pT(n)*MT(1,m)) - pT(3)*pT(m)*MT(1,n) - 
 
590
     -     pT(1)*pT(n)*MT(3,m) - pT(1)*pT(m)*MT(3,n) + 
 
591
     -     (2*pT(1)*pT(3)*MT(m,n))/3.)*T1(2,1)*T2(3,2) + 
 
592
     -  (2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) - 
 
593
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(2,3)*
 
594
     -   T2(3,2) + (2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
595
     -        pT(2)*pT(m)*MT(2,n)) - 
 
596
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(2,3)*T2(3,2)
 
597
     -   + (2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) - 
 
598
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(2,3)*T2(3,2)
 
599
     -   + (-2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) + 
 
600
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(2,3)*T2(3,2)
 
601
     -   + 2*(-(pT(4)*pT(n)*MT(3,m)) - pT(4)*pT(m)*MT(3,n) - 
 
602
     -     pT(3)*pT(n)*MT(4,m) - pT(3)*pT(m)*MT(4,n) + 
 
603
     -     (2*pT(3)*pT(4)*MT(m,n))/3.)*T1(2,4)*T2(3,2) + 
 
604
     -  2*(-(pT(2)*pT(n)*MT(1,m)) - pT(2)*pT(m)*MT(1,n) - 
 
605
     -     pT(1)*pT(n)*MT(2,m) - pT(1)*pT(m)*MT(2,n) + 
 
606
     -     (2*pT(1)*pT(2)*MT(m,n))/3.)*T1(3,1)*T2(3,2) + 
 
607
     -  (-2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) + 
 
608
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(3,2)*
 
609
     -   T2(3,2) + (-2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
610
     -        pT(2)*pT(m)*MT(2,n)) + 
 
611
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(3,2)*T2(3,2)
 
612
     -   + (-2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) + 
 
613
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(3,2)*T2(3,2)
 
614
     -   + (2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) - 
 
615
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(3,2)*T2(3,2)
 
616
     -   - 2*(-(pT(4)*pT(n)*MT(2,m)) - pT(4)*pT(m)*MT(2,n) - 
 
617
     -     pT(2)*pT(n)*MT(4,m) - pT(2)*pT(m)*MT(4,n) + 
 
618
     -     (2*pT(2)*pT(4)*MT(m,n))/3.)*T1(3,4)*T2(3,2) - 
 
619
     -  2*(-(pT(4)*pT(n)*MT(3,m)) - pT(4)*pT(m)*MT(3,n) - 
 
620
     -     pT(3)*pT(n)*MT(4,m) - pT(3)*pT(m)*MT(4,n) + 
 
621
     -     (2*pT(3)*pT(4)*MT(m,n))/3.)*T1(4,2)*T2(3,2) + 
 
622
     -  2*(-(pT(4)*pT(n)*MT(2,m)) - pT(4)*pT(m)*MT(2,n) - 
 
623
     -     pT(2)*pT(n)*MT(4,m) - pT(2)*pT(m)*MT(4,n) + 
 
624
     -     (2*pT(2)*pT(4)*MT(m,n))/3.)*T1(4,3)*T2(3,2) - 
 
625
     -  2*(-(pT(4)*pT(n)*MT(1,m)) - pT(4)*pT(m)*MT(1,n) - 
 
626
     -     pT(1)*pT(n)*MT(4,m) - pT(1)*pT(m)*MT(4,n) + 
 
627
     -     (2*pT(1)*pT(4)*MT(m,n))/3.)*T1(1,3)*T2(3,4) + 
 
628
     -  2*(-(pT(3)*pT(n)*MT(1,m)) - pT(3)*pT(m)*MT(1,n) - 
 
629
     -     pT(1)*pT(n)*MT(3,m) - pT(1)*pT(m)*MT(3,n) + 
 
630
     -     (2*pT(1)*pT(3)*MT(m,n))/3.)*T1(1,4)*T2(3,4) + 
 
631
     -  2*(-(pT(4)*pT(n)*MT(2,m)) - pT(4)*pT(m)*MT(2,n) - 
 
632
     -     pT(2)*pT(n)*MT(4,m) - pT(2)*pT(m)*MT(4,n) + 
 
633
     -     (2*pT(2)*pT(4)*MT(m,n))/3.)*T1(2,3)*T2(3,4) - 
 
634
     -  2*(-(pT(3)*pT(n)*MT(2,m)) - pT(3)*pT(m)*MT(2,n) - 
 
635
     -     pT(2)*pT(n)*MT(3,m) - pT(2)*pT(m)*MT(3,n) + 
 
636
     -     (2*pT(2)*pT(3)*MT(m,n))/3.)*T1(2,4)*T2(3,4) + 
 
637
     -  2*(-(pT(4)*pT(n)*MT(1,m)) - pT(4)*pT(m)*MT(1,n) - 
 
638
     -     pT(1)*pT(n)*MT(4,m) - pT(1)*pT(m)*MT(4,n) + 
 
639
     -     (2*pT(1)*pT(4)*MT(m,n))/3.)*T1(3,1)*T2(3,4) - 
 
640
     -  2*(-(pT(4)*pT(n)*MT(2,m)) - pT(4)*pT(m)*MT(2,n) - 
 
641
     -     pT(2)*pT(n)*MT(4,m) - pT(2)*pT(m)*MT(4,n) + 
 
642
     -     (2*pT(2)*pT(4)*MT(m,n))/3.)*T1(3,2)*T2(3,4) + 
 
643
     -  (-2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) + 
 
644
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(3,4)*
 
645
     -   T2(3,4) + (2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
646
     -        pT(2)*pT(m)*MT(2,n)) - 
 
647
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(3,4)*T2(3,4)
 
648
     -   + (-2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) + 
 
649
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(3,4)*T2(3,4)
 
650
     -   + (-2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) + 
 
651
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(3,4)*T2(3,4)
 
652
     -   - 2*(-(pT(3)*pT(n)*MT(1,m)) - pT(3)*pT(m)*MT(1,n) - 
 
653
     -     pT(1)*pT(n)*MT(3,m) - pT(1)*pT(m)*MT(3,n) + 
 
654
     -     (2*pT(1)*pT(3)*MT(m,n))/3.)*T1(4,1)*T2(3,4) + 
 
655
     -  2*(-(pT(3)*pT(n)*MT(2,m)) - pT(3)*pT(m)*MT(2,n) - 
 
656
     -     pT(2)*pT(n)*MT(3,m) - pT(2)*pT(m)*MT(3,n) + 
 
657
     -     (2*pT(2)*pT(3)*MT(m,n))/3.)*T1(4,2)*T2(3,4) + 
 
658
     -  (2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) - 
 
659
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(4,3)*
 
660
     -   T2(3,4) + (-2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
661
     -        pT(2)*pT(m)*MT(2,n)) + 
 
662
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(4,3)*T2(3,4)
 
663
     -   + (2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) - 
 
664
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(4,3)*T2(3,4)
 
665
     -   + (2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) - 
 
666
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(4,3)*T2(3,4)
 
667
     -   - 2*(-(pT(4)*pT(n)*MT(2,m)) - pT(4)*pT(m)*MT(2,n) - 
 
668
     -     pT(2)*pT(n)*MT(4,m) - pT(2)*pT(m)*MT(4,n) + 
 
669
     -     (2*pT(2)*pT(4)*MT(m,n))/3.)*T1(1,2)*T2(4,1) - 
 
670
     -  2*(-(pT(4)*pT(n)*MT(3,m)) - pT(4)*pT(m)*MT(3,n) - 
 
671
     -     pT(3)*pT(n)*MT(4,m) - pT(3)*pT(m)*MT(4,n) + 
 
672
     -     (2*pT(3)*pT(4)*MT(m,n))/3.)*T1(1,3)*T2(4,1) + 
 
673
     -  (2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) - 
 
674
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(1,4)*
 
675
     -   T2(4,1) + (2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
676
     -        pT(2)*pT(m)*MT(2,n)) - 
 
677
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(1,4)*T2(4,1)
 
678
     -   + (2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) - 
 
679
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(1,4)*T2(4,1)
 
680
     -   + (-2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) + 
 
681
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(1,4)*T2(4,1)
 
682
     -   + 2*(-(pT(4)*pT(n)*MT(2,m)) - pT(4)*pT(m)*MT(2,n) - 
 
683
     -     pT(2)*pT(n)*MT(4,m) - pT(2)*pT(m)*MT(4,n) + 
 
684
     -     (2*pT(2)*pT(4)*MT(m,n))/3.)*T1(2,1)*T2(4,1) - 
 
685
     -  2*(-(pT(2)*pT(n)*MT(1,m)) - pT(2)*pT(m)*MT(1,n) - 
 
686
     -     pT(1)*pT(n)*MT(2,m) - pT(1)*pT(m)*MT(2,n) + 
 
687
     -     (2*pT(1)*pT(2)*MT(m,n))/3.)*T1(2,4)*T2(4,1) + 
 
688
     -  2*(-(pT(4)*pT(n)*MT(3,m)) - pT(4)*pT(m)*MT(3,n) - 
 
689
     -     pT(3)*pT(n)*MT(4,m) - pT(3)*pT(m)*MT(4,n) + 
 
690
     -     (2*pT(3)*pT(4)*MT(m,n))/3.)*T1(3,1)*T2(4,1) - 
 
691
     -  2*(-(pT(3)*pT(n)*MT(1,m)) - pT(3)*pT(m)*MT(1,n) - 
 
692
     -     pT(1)*pT(n)*MT(3,m) - pT(1)*pT(m)*MT(3,n) + 
 
693
     -     (2*pT(1)*pT(3)*MT(m,n))/3.)*T1(3,4)*T2(4,1) + 
 
694
     -  (-2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) + 
 
695
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(4,1)*
 
696
     -   T2(4,1) + (-2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
697
     -        pT(2)*pT(m)*MT(2,n)) + 
 
698
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(4,1)*T2(4,1)
 
699
     -   + (-2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) + 
 
700
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(4,1)*T2(4,1)
 
701
     -   + (2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) - 
 
702
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(4,1)*T2(4,1)
 
703
     -   + 2*(-(pT(2)*pT(n)*MT(1,m)) - pT(2)*pT(m)*MT(1,n) - 
 
704
     -     pT(1)*pT(n)*MT(2,m) - pT(1)*pT(m)*MT(2,n) + 
 
705
     -     (2*pT(1)*pT(2)*MT(m,n))/3.)*T1(4,2)*T2(4,1) + 
 
706
     -  2*(-(pT(3)*pT(n)*MT(1,m)) - pT(3)*pT(m)*MT(1,n) - 
 
707
     -     pT(1)*pT(n)*MT(3,m) - pT(1)*pT(m)*MT(3,n) + 
 
708
     -     (2*pT(1)*pT(3)*MT(m,n))/3.)*T1(4,3)*T2(4,1) + 
 
709
     -  2*(-(pT(4)*pT(n)*MT(1,m)) - pT(4)*pT(m)*MT(1,n) - 
 
710
     -     pT(1)*pT(n)*MT(4,m) - pT(1)*pT(m)*MT(4,n) + 
 
711
     -     (2*pT(1)*pT(4)*MT(m,n))/3.)*T1(1,2)*T2(4,2) - 
 
712
     -  2*(-(pT(2)*pT(n)*MT(1,m)) - pT(2)*pT(m)*MT(1,n) - 
 
713
     -     pT(1)*pT(n)*MT(2,m) - pT(1)*pT(m)*MT(2,n) + 
 
714
     -     (2*pT(1)*pT(2)*MT(m,n))/3.)*T1(1,4)*T2(4,2) - 
 
715
     -  2*(-(pT(4)*pT(n)*MT(1,m)) - pT(4)*pT(m)*MT(1,n) - 
 
716
     -     pT(1)*pT(n)*MT(4,m) - pT(1)*pT(m)*MT(4,n) + 
 
717
     -     (2*pT(1)*pT(4)*MT(m,n))/3.)*T1(2,1)*T2(4,2) + 
 
718
     -  2*(-(pT(4)*pT(n)*MT(3,m)) - pT(4)*pT(m)*MT(3,n) - 
 
719
     -     pT(3)*pT(n)*MT(4,m) - pT(3)*pT(m)*MT(4,n) + 
 
720
     -     (2*pT(3)*pT(4)*MT(m,n))/3.)*T1(2,3)*T2(4,2) + 
 
721
     -  (2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) - 
 
722
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(2,4)*
 
723
     -   T2(4,2) + (2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
724
     -        pT(2)*pT(m)*MT(2,n)) - 
 
725
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(2,4)*T2(4,2)
 
726
     -   + (-2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) + 
 
727
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(2,4)*T2(4,2)
 
728
     -   + (2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) - 
 
729
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(2,4)*T2(4,2)
 
730
     -   - 2*(-(pT(4)*pT(n)*MT(3,m)) - pT(4)*pT(m)*MT(3,n) - 
 
731
     -     pT(3)*pT(n)*MT(4,m) - pT(3)*pT(m)*MT(4,n) + 
 
732
     -     (2*pT(3)*pT(4)*MT(m,n))/3.)*T1(3,2)*T2(4,2) + 
 
733
     -  2*(-(pT(3)*pT(n)*MT(2,m)) - pT(3)*pT(m)*MT(2,n) - 
 
734
     -     pT(2)*pT(n)*MT(3,m) - pT(2)*pT(m)*MT(3,n) + 
 
735
     -     (2*pT(2)*pT(3)*MT(m,n))/3.)*T1(3,4)*T2(4,2) + 
 
736
     -  2*(-(pT(2)*pT(n)*MT(1,m)) - pT(2)*pT(m)*MT(1,n) - 
 
737
     -     pT(1)*pT(n)*MT(2,m) - pT(1)*pT(m)*MT(2,n) + 
 
738
     -     (2*pT(1)*pT(2)*MT(m,n))/3.)*T1(4,1)*T2(4,2) + 
 
739
     -  (-2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) + 
 
740
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(4,2)*
 
741
     -   T2(4,2) + (-2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
742
     -        pT(2)*pT(m)*MT(2,n)) + 
 
743
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(4,2)*T2(4,2)
 
744
     -   + (2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) - 
 
745
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(4,2)*T2(4,2)
 
746
     -   + (-2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) + 
 
747
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(4,2)*T2(4,2)
 
748
     -   - 2*(-(pT(3)*pT(n)*MT(2,m)) - pT(3)*pT(m)*MT(2,n) - 
 
749
     -     pT(2)*pT(n)*MT(3,m) - pT(2)*pT(m)*MT(3,n) + 
 
750
     -     (2*pT(2)*pT(3)*MT(m,n))/3.)*T1(4,3)*T2(4,2) + 
 
751
     -  2*(-(pT(4)*pT(n)*MT(1,m)) - pT(4)*pT(m)*MT(1,n) - 
 
752
     -     pT(1)*pT(n)*MT(4,m) - pT(1)*pT(m)*MT(4,n) + 
 
753
     -     (2*pT(1)*pT(4)*MT(m,n))/3.)*T1(1,3)*T2(4,3) - 
 
754
     -  2*(-(pT(3)*pT(n)*MT(1,m)) - pT(3)*pT(m)*MT(1,n) - 
 
755
     -     pT(1)*pT(n)*MT(3,m) - pT(1)*pT(m)*MT(3,n) + 
 
756
     -     (2*pT(1)*pT(3)*MT(m,n))/3.)*T1(1,4)*T2(4,3) - 
 
757
     -  2*(-(pT(4)*pT(n)*MT(2,m)) - pT(4)*pT(m)*MT(2,n) - 
 
758
     -     pT(2)*pT(n)*MT(4,m) - pT(2)*pT(m)*MT(4,n) + 
 
759
     -     (2*pT(2)*pT(4)*MT(m,n))/3.)*T1(2,3)*T2(4,3) + 
 
760
     -  2*(-(pT(3)*pT(n)*MT(2,m)) - pT(3)*pT(m)*MT(2,n) - 
 
761
     -     pT(2)*pT(n)*MT(3,m) - pT(2)*pT(m)*MT(3,n) + 
 
762
     -     (2*pT(2)*pT(3)*MT(m,n))/3.)*T1(2,4)*T2(4,3) - 
 
763
     -  2*(-(pT(4)*pT(n)*MT(1,m)) - pT(4)*pT(m)*MT(1,n) - 
 
764
     -     pT(1)*pT(n)*MT(4,m) - pT(1)*pT(m)*MT(4,n) + 
 
765
     -     (2*pT(1)*pT(4)*MT(m,n))/3.)*T1(3,1)*T2(4,3) + 
 
766
     -  2*(-(pT(4)*pT(n)*MT(2,m)) - pT(4)*pT(m)*MT(2,n) - 
 
767
     -     pT(2)*pT(n)*MT(4,m) - pT(2)*pT(m)*MT(4,n) + 
 
768
     -     (2*pT(2)*pT(4)*MT(m,n))/3.)*T1(3,2)*T2(4,3) + 
 
769
     -  (2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) - 
 
770
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(3,4)*
 
771
     -   T2(4,3) + (-2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
772
     -        pT(2)*pT(m)*MT(2,n)) + 
 
773
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(3,4)*T2(4,3)
 
774
     -   + (2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) - 
 
775
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(3,4)*T2(4,3)
 
776
     -   + (2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) - 
 
777
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)*T1(3,4)*T2(4,3)
 
778
     -   + 2*(-(pT(3)*pT(n)*MT(1,m)) - pT(3)*pT(m)*MT(1,n) - 
 
779
     -     pT(1)*pT(n)*MT(3,m) - pT(1)*pT(m)*MT(3,n) + 
 
780
     -     (2*pT(1)*pT(3)*MT(m,n))/3.)*T1(4,1)*T2(4,3) - 
 
781
     -  2*(-(pT(3)*pT(n)*MT(2,m)) - pT(3)*pT(m)*MT(2,n) - 
 
782
     -     pT(2)*pT(n)*MT(3,m) - pT(2)*pT(m)*MT(3,n) + 
 
783
     -     (2*pT(2)*pT(3)*MT(m,n))/3.)*T1(4,2)*T2(4,3) + 
 
784
     -  (-2*(-(pT(1)*pT(n)*MT(1,m)) - pT(1)*pT(m)*MT(1,n)) + 
 
785
     -     (2*(-(pT(m)*pT(n)) - pT(1)**2*MT(m,n)))/3.)*T1(4,3)*
 
786
     -   T2(4,3) + (2*(-(pT(2)*pT(n)*MT(2,m)) - 
 
787
     -        pT(2)*pT(m)*MT(2,n)) - 
 
788
     -     (2*(pT(m)*pT(n) - pT(2)**2*MT(m,n)))/3.)*T1(4,3)*T2(4,3)
 
789
     -   + (-2*(-(pT(3)*pT(n)*MT(3,m)) - pT(3)*pT(m)*MT(3,n)) + 
 
790
     -     (2*(pT(m)*pT(n) - pT(3)**2*MT(m,n)))/3.)*T1(4,3)*T2(4,3)
 
791
     -   + (-2*(-(pT(4)*pT(n)*MT(4,m)) - pT(4)*pT(m)*MT(4,n)) + 
 
792
     -     (2*(pT(m)*pT(n) - pT(4)**2*MT(m,n)))/3.)
 
793
     -    *T1(4,3)*T2(4,3))/tmass*2
 
794
 
 
795
           enddo
 
796
         enddo
 
797
 
 
798
 
 
799
 
 
800
        if ( tmass.gt.rZero ) then
 
801
         d = - 1.0d0/dcmplx( pT2-tmass**2, tmass*twidth )
 
802
        else
 
803
         d = - 1.0d0/dcmplx( pT2, rZero )
 
804
        end if
 
805
 
 
806
         do m=1,4
 
807
           do n=1,4
 
808
          T(m,n)=T(m,n)*gt*d/2.0d0
 
809
           enddo
 
810
         enddo
 
811
c     2.0 factor from propagator convention
 
812
 
 
813
         
 
814
      wt(1) = T(1,1)
 
815
      wt(2) = T(1,2)
 
816
      wt(3) = T(1,3)
 
817
      wt(4) = T(1,4)
 
818
      wt(5) = T(2,1)
 
819
      wt(6) = T(2,2)
 
820
      wt(7) = T(2,3)
 
821
      wt(8) = T(2,4)
 
822
      wt(9) = T(3,1)
 
823
      wt(10) = T(3,2)
 
824
      wt(11) = T(3,3)
 
825
      wt(12) = T(3,4)
 
826
      wt(13) = T(4,1)
 
827
      wt(14) = T(4,2)
 
828
      wt(15) = T(4,3)
 
829
      wt(16) = T(4,4)
 
830
      wt(17) = T(5,1)
 
831
      wt(18) = T(6,1)
 
832
 
 
833
      return
 
834
      end