← Back to branch summary

~maddevelopers/mg5amcnlo/3.0.2-alpha0

« back to all changes in this revision

Viewing changes to HELAS/uioxxx.F

Added Template and HELAS into bzr

Show diffs side-by-side

added added

removed removed

Lines of Context:
HELAS/uioxxx.F
 
1
      subroutine uioxxx(fi,fo,gt,fmass,tmass,twidth , uio)
 
2
c
 
3
c This subroutine computes an off-shell tensor current from
 
4
c the fermion-anti-fermion-tensor boson coupling.
 
5
c
 
6
c input:
 
7
c       complex fi(6)          : flow-in  fermion                   |fi>
 
8
c       complex fo(6)          : flow-out fermion                   <fo|
 
9
c       complex gt             : coupling constant       gtf=-1/Lambda/4
 
10
c       real    fmass          : mass  of input fermion
 
11
c       real    tmass          : mass  of output tensor T
 
12
c       real    twidth         : width of output tensor T
 
13
c
 
14
c output:
 
15
c       complex uio(18)        : tensor current       j^mu^nu(<fo|T|fi>)
 
16
c
 
17
c- by Q.Li - OCT. 2006
 
18
c- Added massless tensor - P. de Aquino - Oct. 2009 
 
19
c
 
20
      implicit none
 
21
      double complex fi(6), fo(6), gt, uio(18)
 
22
      double precision fmass, tmass, twidth
 
23
 
 
24
      double complex yio(6,4)
 
25
      integer i,j
 
26
      double precision pi(4),po(4),km(4),kp(4)
 
27
      double precision MET(4,4)
 
28
      double complex cone,cZero, d, tt1,tt2,tt3
 
29
      double precision rZero, rTwo
 
30
      double precision KT2,K1KT,K2KT
 
31
      parameter( rZero = 0.0d0, rTwo = 2.0d0 )
 
32
      parameter( cZero = ( 0.0d0, 0.0d0 ), cone=(0.0d0,1.0d0) )
 
33
 
 
34
      
 
35
      yio(5,1) = -fi(5)+fo(5)
 
36
      yio(6,1) = -fi(6)+fo(6)
 
37
      
 
38
      pi(1) = dreal(fi(5))
 
39
      pi(2) = dreal(fi(6))
 
40
      pi(3) = dimag(fi(6))
 
41
      pi(4) = dimag(fi(5))
 
42
      
 
43
      po(1) = dreal(fo(5))
 
44
      po(2) = dreal(fo(6))
 
45
      po(3) = dimag(fo(6))
 
46
      po(4) = dimag(fo(5))
 
47
      
 
48
      km(1) = dreal(yio(5,1))
 
49
      km(2) = dreal(yio(6,1))
 
50
      km(3) = dimag(yio(6,1))
 
51
      km(4) = dimag(yio(5,1))
 
52
 
 
53
      kp(1) = po(1)+pi(1)
 
54
      kp(2) = po(2)+pi(2)
 
55
      kp(3) = po(3)+pi(3)
 
56
      kp(4) = po(4)+pi(4)
 
57
      
 
58
      do i=1,4
 
59
         do j=1,4
 
60
            MET(i,j) = 0.0d0
 
61
         enddo 
 
62
      enddo
 
63
      MET(1,1) =  1.0d0
 
64
      MET(2,2) = -1.0d0
 
65
      MET(3,3) = -1.0d0
 
66
      MET(4,4) = -1.0d0
 
67
      
 
68
      K1KT = pi(1)*km(1)-pi(2)*km(2)-pi(3)*km(3)-pi(4)*km(4)
 
69
      K2KT = po(1)*km(1)-po(2)*km(2)-po(3)*km(3)-po(4)*km(4)
 
70
      KT2 = km(1)**2-km(2)**2-km(3)**2-km(4)**2
 
71
      
 
72
      if ( tmass.eq.rZero ) then
 
73
         d = - gt/dcmplx( KT2, rZero )
 
74
 
 
75
         tt1 = fi(3)*(-4*fmass*fo(3) + fo(2)*(-kp(2) - cone*kp(3)) 
 
76
     &+ fo(1)*(kp(1) - kp(4)))  
 
77
     &+ fi(2)*(-4*fmass*fo(2) + fo(3)*(kp(2) - cone*kp(3)) 
 
78
     &+ fo(4)*(kp(1) - kp(4)))  
 
79
     &+ fi(4)*(-4*fmass*fo(4) + fo(1)*(-kp(2) + cone*kp(3)) 
 
80
     &+ fo(2)*(kp(1) + kp(4)))  
 
81
     &+ fi(1)*(-4*fmass*fo(1) + fo(4)*(kp(2) + cone*kp(3)) 
 
82
     &+ fo(3)*(kp(1) + kp(4)))
 
83
      
 
84
         do i = 1,4
 
85
            do j=1,4
 
86
               yio(i,j) = tt1*MET(i,j)
 
87
            end do
 
88
         enddo
 
89
        
 
90
         yio(1,1) = yio(1,1)
 
91
     &+ 2*fi(3)*fo(1)*kp(1) + 2*fi(4)*fo(2)*kp(1)  
 
92
     &+ 2*fi(1)*fo(3)*kp(1) + 2*fi(2)*fo(4)*kp(1)
 
93
        
 
94
         yio(2,2) = yio(2,2)
 
95
     &+ 2*fi(4)*fo(1)*kp(2) + 2*fi(3)*fo(2)*kp(2)  
 
96
     &- 2*fi(2)*fo(3)*kp(2) - 2*fi(1)*fo(4)*kp(2)
 
97
   
 
98
         yio(3,3) = yio(3,3)
 
99
     &- 2*cone*fi(4)*fo(1)*kp(3) + 2*cone*fi(3)*fo(2)*kp(3) 
 
100
     &+ 2*cone*fi(2)*fo(3)*kp(3) - 2*cone*fi(1)*fo(4)*kp(3)
 
101
        
 
102
         yio(4,4) = yio(4,4)
 
103
     &+ 2*fi(3)*fo(1)*kp(4) - 2*fi(4)*fo(2)*kp(4) 
 
104
     &- 2*fi(1)*fo(3)*kp(4) + 2*fi(2)*fo(4)*kp(4)       
 
105
   
 
106
         yio(1,2) = yio(1,2)
 
107
     &+ fi(3)*(fo(2)*kp(1) + fo(1)*kp(2)) 
 
108
     &+ fi(4)*(fo(1)*kp(1) + fo(2)*kp(2)) 
 
109
     &+ fi(1)*(-(fo(4)*kp(1))+ fo(3)*kp(2)) 
 
110
     &+ fi(2)*(-(fo(3)*kp(1)) + fo(4)*kp(2))
 
111
   
 
112
         yio(1,3) = yio(1,3)
 
113
     &+ fi(3)*(cone*fo(2)*kp(1) + fo(1)*kp(3)) 
 
114
     &+ fi(4)*(-(cone*fo(1)*kp(1)) + fo(2)*kp(3)) 
 
115
     &+ fi(1)*(-(cone*fo(4)*kp(1)) + fo(3)*kp(3)) 
 
116
     &+ fi(2)*(cone*fo(3)*kp(1) + fo(4)*kp(3))
 
117
     
 
118
         yio(1,4) = yio(1,4)
 
119
     &+ fi(4)*fo(2)*(-kp(1) + kp(4)) 
 
120
     &+ fi(1)*fo(3)*(-kp(1) + kp(4)) 
 
121
     &+ fi(3)*fo(1)*(kp(1) + kp(4))
 
122
     &+ fi(2)*fo(4)*(kp(1) + kp(4))
 
123
     
 
124
         yio(2,3) = yio(2,3)
 
125
     &+ fi(1)*fo(4)*(-(cone*kp(2)) - kp(3)) 
 
126
     &+ fi(2)*fo(3)*(cone*kp(2) - kp(3))  
 
127
     &+ fi(4)*fo(1)*(-(cone*kp(2)) + kp(3)) 
 
128
     &+ fi(3)*fo(2)*(cone*kp(2) + kp(3))
 
129
   
 
130
         yio(2,4) = yio(2,4)
 
131
     &+ fi(4)*(-(fo(2)*kp(2)) + fo(1)*kp(4)) 
 
132
     &+ fi(3)*(fo(1)*kp(2) + fo(2)*kp(4))  
 
133
     &+ fi(2)*(fo(4)*kp(2) - fo(3)*kp(4)) 
 
134
     &+ fi(1)*(-(fo(3)*kp(2)) - fo(4)*kp(4))
 
135
     
 
136
         yio(3,4) = yio(3,4)
 
137
     &+ fi(4)*(-(fo(2)*kp(3)) - cone*fo(1)*kp(4)) 
 
138
     &+ fi(3)*(fo(1)*kp(3) + cone*fo(2)*kp(4)) 
 
139
     &+ fi(2)*(fo(4)*kp(3) + cone*fo(3)*kp(4)) 
 
140
     &+ fi(1)*(-(fo(3)*kp(3)) - cone*fo(4)*kp(4))
 
141
   
 
142
         yio(2,1) = yio(1,2)
 
143
         yio(3,1) = yio(1,3)
 
144
         yio(4,1) = yio(1,4)         
 
145
         yio(3,2) = yio(2,3)
 
146
         yio(4,2) = yio(2,4)
 
147
         yio(4,3) = yio(3,4)
 
148
        
 
149
         do i = 1,4
 
150
            do j=1,4
 
151
               yio(i,j) = yio(i,j)*d
 
152
            end do
 
153
         enddo
 
154
   
 
155
      else if (tmass.gt.rZero) then
 
156
           d = - gt/dcmplx( KT2-tmass**2, tmass*twidth )
 
157
            
 
158
           tt1 = fi(3)*(((-8.0d0*fmass)/3.d0 
 
159
     &  + (8.0d0*fmass*KT2)/(3.d0*tmass**2))*fo(3)
 
160
     & +  fo(2)*((4.0d0*(K1KT + K2KT)
 
161
     &            *(-km(2) - cone*km(3)))/(3.d0*tmass**2)
 
162
     & - (4.0d0*KT2*(-kp(2) - cone*kp(3)))/(3.d0*tmass**2)) + 
 
163
     &    fo(1)*((4.0d0*(K1KT + K2KT)*(km(1) - km(4)))/(3.d0*tmass**2) 
 
164
     & - (4.0d0*KT2*(kp(1) - kp(4)))/(3.d0*tmass**2))) + 
 
165
     &    fi(2)*(((-8.0d0*fmass)/3.d0 
 
166
     &        + (8.0d0*fmass*KT2)/(3.d0*tmass**2))*fo(2)  
 
167
     & +  fo(3)*((4.0d0*(K1KT + K2KT)*(km(2) 
 
168
     &          - cone*km(3)))/(3.d0*tmass**2) 
 
169
     &- (4.0d0*KT2*(kp(2) - cone*kp(3)))/(3.d0*tmass**2)) + 
 
170
     &    fo(4)*((4.0d0*(K1KT + K2KT)*(km(1) - km(4)))/(3.d0*tmass**2)
 
171
     & - (4.0d0*KT2*(kp(1) - kp(4)))/(3.d0*tmass**2))) + 
 
172
     &    fi(4)*(((-8.0d0*fmass)/3.d0 
 
173
     &+ (8.0d0*fmass*KT2)/(3.d0*tmass**2))*fo(4)  
 
174
     & +  fo(1)*((4.0d0*(K1KT + K2KT)
 
175
     &       *(-km(2) + cone*km(3)))/(3.d0*tmass**2)
 
176
     & - (4.0d0*KT2*(-kp(2) + cone*kp(3)))/(3.d0*tmass**2)) + 
 
177
     &    fo(2)*((4.0d0*(K1KT + K2KT)*(km(1) + km(4)))/(3.d0*tmass**2)
 
178
     & - (4.0d0*KT2*(kp(1) + kp(4)))/(3.d0*tmass**2))) + 
 
179
     &    fi(1)*(((-8.0d0*fmass)/3.d0 
 
180
     &+ (8.0d0*fmass*KT2)/(3.d0*tmass**2))*fo(1) 
 
181
     & +  fo(4)*((4.0d0*(K1KT + K2KT)*(km(2)
 
182
     &      + cone*km(3)))/(3.d0*tmass**2) 
 
183
     &- (4.0d0*KT2*(kp(2) + cone*kp(3)))/(3.d0*tmass**2)) + 
 
184
     &    fo(3)*((4.0d0*(K1KT + K2KT)
 
185
     &*(km(1) + km(4)))/(3.d0*tmass**2)
 
186
     & - (4.0d0*KT2*(kp(1) + kp(4)))/(3.d0*tmass**2)))
 
187
 
 
188
      
 
189
           tt2 = (fi(3)*(((-16.0d0*KT2*fmass)/(3.0D0*tmass**4)
 
190
     &   - (16.d0*fmass)/(3.d0*tmass**2))*fo(3) + 
 
191
     &       fo(2)*((8.d0*(K1KT + K2KT)*(-km(2)
 
192
     & - cone*km(3)))/(3.d0*tmass**4) 
 
193
     &+ (8.d0*KT2*(-kp(2) - cone*kp(3)))/(3.0D0*tmass**4) + 
 
194
     &          (4.d0*(-kp(2) - cone*kp(3)))/tmass**2) + 
 
195
     &       fo(1)*((8.d0*(K1KT + K2KT)*(km(1)
 
196
     & - km(4)))/(3.d0*tmass**4) 
 
197
     &+ (8.d0*KT2*(kp(1) - kp(4)))/(3.d0*tmass**4) + 
 
198
     &          (4.d0*(kp(1) - kp(4)))/tmass**2)) + fi(2)*
 
199
     &     (((-16.d0*KT2*fmass)/(3.d0*tmass**4) 
 
200
     &- (16.d0*fmass)/(3.d0*tmass**2))*fo(2) + 
 
201
     &       fo(3)*((8.d0*(K1KT + K2KT)*(km(2)
 
202
     & - cone*km(3)))/(3.d0*tmass**4) 
 
203
     &+ (8.d0*KT2*(kp(2) - cone*kp(3)))/(3.d0*tmass**4) + 
 
204
     &          (4.d0*(kp(2) - cone*kp(3)))/tmass**2) + 
 
205
     &       fo(4)*((8.d0*(K1KT + K2KT)*(km(1) - km(4)))/(3.d0*tmass**4)
 
206
     & + (8.d0*KT2*(kp(1) - kp(4)))/(3.d0*tmass**4) + 
 
207
     &          (4.d0*(kp(1) - kp(4)))/tmass**2)) + fi(4)*
 
208
     &     (((-16.d0*KT2*fmass)/(3.d0*tmass**4) 
 
209
     &- (16.d0*fmass)/(3.d0*tmass**2))*fo(4) + 
 
210
     &       fo(1)*((8.d0*(K1KT + K2KT)*(-km(2) 
 
211
     &+ cone*km(3)))/(3.d0*tmass**4)
 
212
     & + (8.d0*KT2*(-kp(2) + cone*kp(3)))/(3.d0*tmass**4) + 
 
213
     &          (4.d0*(-kp(2) + cone*kp(3)))/tmass**2) + 
 
214
     &       fo(2)*((8.d0*(K1KT + K2KT)*(km(1) 
 
215
     &+ km(4)))/(3.d0*tmass**4)
 
216
     & + (8.d0*KT2*(kp(1) + kp(4)))/(3.d0*tmass**4) + 
 
217
     &          (4.d0*(kp(1) + kp(4)))/tmass**2)) + fi(1)*
 
218
     &     (((-16.d0*KT2*fmass)/(3.d0*tmass**4) 
 
219
     &- (16.d0*fmass)/(3.d0*tmass**2))*fo(1) + 
 
220
     &       fo(4)*((8.d0*(K1KT + K2KT)*(km(2)
 
221
     & + cone*km(3)))/(3.d0*tmass**4)
 
222
     & + (8.d0*KT2*(kp(2) + cone*kp(3)))/(3.d0*tmass**4) + 
 
223
     &          (4.d0*(kp(2) + cone*kp(3)))/tmass**2) + 
 
224
     &       fo(3)*((8.d0*(K1KT + K2KT)
 
225
     &*(km(1) + km(4)))/(3.d0*tmass**4)
 
226
     & + (8.d0*KT2*(kp(1) + kp(4)))/(3.d0*tmass**4) + 
 
227
     &          (4.d0*(kp(1) + kp(4)))/tmass**2)))
 
228
   
 
229
           tt3 = fi(3)*((-2.0d0*fo(2)*(-km(2) - cone*km(3)))/tmass**2
 
230
     &   - (2.0d0*fo(1)*(km(1) - km(4)))/tmass**2) + 
 
231
     &    fi(2)*((-2.0d0*fo(3)*(km(2) - cone*km(3)))/tmass**2 
 
232
     &- (2.0d0*fo(4)*(km(1) - km(4)))/tmass**2) + 
 
233
     &    fi(4)*((-2.0d0*fo(1)*(-km(2) + cone*km(3)))/tmass**2
 
234
     & - (2.0d0*fo(2)*(km(1) + km(4)))/tmass**2) + 
 
235
     &    fi(1)*((-2.0d0*fo(4)*(km(2) + cone*km(3)))/tmass**2 
 
236
     &- (2.0d0*fo(3)*(km(1) + km(4)))/tmass**2) 
 
237
         
 
238
           do i = 1,4
 
239
              do j=1,4
 
240
                 yio(i,j) = tt1*MET(i,j)+tt2*km(i)*km(j)
 
241
     &+tt3*(km(i)*kp(j)+km(j)*kp(i))
 
242
              end do
 
243
           enddo
 
244
        
 
245
           yio(1,1) = yio(1,1)
 
246
     &+fi(3)*fo(1)*((-4.0d0*(K1KT + K2KT)*km(1))/tmass**2 + 4.0d0*kp(1))  
 
247
     & +fi(4)*fo(2)*((-4.d0*(K1KT + K2KT)*km(1))/tmass**2 + 4.0d0*kp(1)) 
 
248
     & +fi(1)*fo(3)*((-4.0d0*(K1KT + K2KT)*km(1))/tmass**2+ 4.0d0*kp(1)) 
 
249
     &+ fi(2)*fo(4)*((-4.0d0*(K1KT + K2KT)*km(1))/tmass**2+ 4.0d0*kp(1))
 
250
        
 
251
           yio(2,2) = yio(2,2)
 
252
     &+fi(2)*fo(3)*((4.0d0*(K1KT + K2KT)*km(2))/tmass**2 - 4.d0*kp(2))
 
253
     &+fi(1)*fo(4)*((4.0d0*(K1KT + K2KT)*km(2))/tmass**2 - 4.d0*kp(2))  
 
254
     &+fi(4)*fo(1)*((-4.d0*(K1KT + K2KT)*km(2))/tmass**2 + 4.d0*kp(2))
 
255
     &+fi(3)*fo(2)*((-4.d0*(K1KT + K2KT)*km(2))/tmass**2 + 4.d0*kp(2))
 
256
      
 
257
           yio(3,3) = yio(3,3)
 
258
     &+fi(4)*fo(1)*((4.0d0*cone*(K1KT + K2KT)*km(3))/tmass**2
 
259
     & - 4.0d0*cone*kp(3)) + 
 
260
     &    fi(1)*fo(4)*((4.0d0*cone*(K1KT + K2KT)*km(3))/tmass**2
 
261
     & - 4.0d0*cone*kp(3)) + 
 
262
     &    fi(3)*fo(2)*((-4.0d0*cone*(K1KT + K2KT)*km(3))/tmass**2
 
263
     & + 4.d0*cone*kp(3)) + 
 
264
     &    fi(2)*fo(3)*((-4.d0*cone*(K1KT + K2KT)*km(3))/tmass**2
 
265
     & + 4.d0*cone*kp(3))       
 
266
        
 
267
           yio(4,4) = yio(4,4)
 
268
     &+fi(4)*fo(2)*((4.d0*(K1KT + K2KT)*km(4))/tmass**2 - 4.d0*kp(4)) 
 
269
     &+fi(1)*fo(3)*((4.d0*(K1KT + K2KT)*km(4))/tmass**2 - 4.d0*kp(4)) 
 
270
     &+fi(3)*fo(1)*((-4.d0*(K1KT + K2KT)*km(4))/tmass**2+ 4.d0*kp(4))
 
271
     &+fi(2)*fo(4)*((-4.d0*(K1KT + K2KT)*km(4))/tmass**2+ 4.d0*kp(4))   
 
272
         
 
273
           yio(1,2) = yio(1,2)
 
274
     &+fi(3)*(fo(2)*((-2.d0*(K1KT + K2KT)*km(1))/tmass**2 + 2.d0*kp(1)) 
 
275
     &+ fo(1)*((-2.d0*(K1KT + K2KT)*km(2))/tmass**2 + 2.d0*kp(2))) + 
 
276
     &  fi(4)*(fo(1)*((-2.d0*(K1KT + K2KT)*km(1))/tmass**2 + 2.d0*kp(1))
 
277
     & + fo(2)*((-2.d0*(K1KT + K2KT)*km(2))/tmass**2 + 2.d0*kp(2))) + 
 
278
     &  fi(1)*(fo(4)*((2.d0*(K1KT + K2KT)*km(1))/tmass**2 - 2.d0*kp(1))
 
279
     & + fo(3)*((-2.d0*(K1KT + K2KT)*km(2))/tmass**2 + 2.d0*kp(2))) + 
 
280
     & fi(2)*(fo(3)*((2.d0*(K1KT + K2KT)*km(1))/tmass**2 - 2.d0*kp(1))
 
281
     & + fo(4)*((-2.d0*(K1KT + K2KT)*km(2))/tmass**2 + 2.d0*kp(2)))
 
282
      
 
283
           yio(1,3) = yio(1,3)
 
284
     &+fi(3)*(fo(2)*((-2.0d0*cone*(K1KT + K2KT)*km(1))/tmass**2
 
285
     & + 2.d0*cone*kp(1)) + 
 
286
     &    fo(1)*((-2.d0*(K1KT + K2KT)*km(3))/tmass**2 +2.d0*kp(3))) + 
 
287
     &    fi(4)*(fo(1)*((2.d0*cone*(K1KT + K2KT)*km(1))/tmass**2
 
288
     & - 2.d0*cone*kp(1)) + 
 
289
     &    fo(2)*((-2.d0*(K1KT + K2KT)*km(3))/tmass**2 +2.d0*kp(3))) + 
 
290
     &    fi(1)*(fo(4)*((2.d0*cone*(K1KT + K2KT)*km(1))/tmass**2 
 
291
     &- 2.d0*cone*kp(1)) + 
 
292
     &    fo(3)*((-2.d0*(K1KT + K2KT)*km(3))/tmass**2 +2.d0*kp(3))) + 
 
293
     &    fi(2)*(fo(3)*((-2.d0*cone*(K1KT + K2KT)*km(1))/tmass**2 
 
294
     &+ 2.d0*cone*kp(1)) + 
 
295
     &    fo(4)*((-2.d0*(K1KT + K2KT)*km(3))/tmass**2 + 2.d0*kp(3)))
 
296
        
 
297
           yio(1,4) = yio(1,4)
 
298
     &+fi(4)*fo(2)*((-2.d0*(K1KT + K2KT)*(-km(1) + km(4)))/tmass**2
 
299
     & - 2.d0*kp(1) + 2.d0*kp(4)) + 
 
300
     & fi(1)*fo(3)*((-2.d0*(K1KT + K2KT)*(-km(1) + km(4)))/tmass**2
 
301
     & - 2.d0*kp(1) + 2.d0*kp(4)) + 
 
302
     & fi(3)*fo(1)*((-2.d0*(K1KT + K2KT)*(km(1) + km(4)))/tmass**2 
 
303
     &+ 2.d0*kp(1) + 2.d0*kp(4)) + 
 
304
     & fi(2)*fo(4)*((-2.d0*(K1KT + K2KT)*(km(1) + km(4)))/tmass**2
 
305
     & + 2.d0*kp(1) + 2.d0*kp(4))
 
306
      
 
307
           yio(2,3) = yio(2,3)
 
308
     &+fi(1)*fo(4)*((-2.0d0*(K1KT + K2KT)*(-(cone*km(2))
 
309
     & - km(3)))/tmass**2 - 2.0d0*cone*kp(2) - 2.d0*kp(3)) + 
 
310
     &    fi(2)*fo(3)*((-2.d0*(K1KT + K2KT)*(cone*km(2)
 
311
     & - km(3)))/tmass**2 + 2.d0*cone*kp(2) - 2.d0*kp(3)) + 
 
312
     &    fi(4)*fo(1)*((-2.d0*(K1KT + K2KT)*(-(cone*km(2))
 
313
     & + km(3)))/tmass**2 - 2.d0*cone*kp(2) + 2.d0*kp(3)) + 
 
314
     &    fi(3)*fo(2)*((-2.d0*(K1KT + K2KT)*(cone*km(2)
 
315
     & + km(3)))/tmass**2 + 2.d0*cone*kp(2) + 2.d0*kp(3))
 
316
      
 
317
      
 
318
           yio(2,4) = yio(2,4)
 
319
     &+fi(2)*(fo(4)*((-2.d0*(K1KT + K2KT)*km(2))/tmass**2+2.d0*kp(2)) + 
 
320
     &       fo(3)*((2.d0*(K1KT + K2KT)*km(4))/tmass**2 - 2.d0*kp(4))) + 
 
321
     &  fi(1)*(fo(3)*((2.d0*(K1KT + K2KT)*km(2))/tmass**2 - 2.d0*kp(2))
 
322
     & + fo(4)*((2.d0*(K1KT + K2KT)*km(4))/tmass**2 - 2.d0*kp(4))) + 
 
323
     &  fi(4)*(fo(2)*((2.d0*(K1KT + K2KT)*km(2))/tmass**2 - 2.d0*kp(2))
 
324
     & + fo(1)*((-2.d0*(K1KT + K2KT)*km(4))/tmass**2 + 2.d0*kp(4))) + 
 
325
     & fi(3)*(fo(1)*((-2.d0*(K1KT + K2KT)*km(2))/tmass**2 + 2.d0*kp(2))
 
326
     & + fo(2)*((-2.d0*(K1KT + K2KT)*km(4))/tmass**2 + 2.d0*kp(4)))
 
327
         
 
328
           yio(3,4) = yio(3,4)
 
329
     &+fi(4)*(fo(2)*((2.d0*(K1KT + K2KT)*km(3))/tmass**2- 2.d0*kp(3)) + 
 
330
     &       fo(1)*((2.d0*cone*(K1KT + K2KT)*km(4))/tmass**2
 
331
     & - 2.d0*cone*kp(4))) + 
 
332
     &    fi(1)*(fo(3)*((2.d0*(K1KT + K2KT)*km(3))/tmass**2
 
333
     & - 2.d0*kp(3)) + 
 
334
     &       fo(4)*((2.d0*cone*(K1KT + K2KT)*km(4))/tmass**2
 
335
     & - 2.d0*cone*kp(4))) + 
 
336
     &    fi(3)*(fo(1)*((-2.d0*(K1KT + K2KT)*km(3))/tmass**2
 
337
     & + 2.d0*kp(3)) + 
 
338
     &       fo(2)*((-2.d0*cone*(K1KT + K2KT)*km(4))/tmass**2
 
339
     & + 2.d0*cone*kp(4))) + 
 
340
     &    fi(2)*(fo(4)*((-2.d0*(K1KT + K2KT)*km(3))/tmass**2
 
341
     & + 2.d0*kp(3)) + 
 
342
     &       fo(3)*((-2.d0*cone*(K1KT + K2KT)*km(4))/tmass**2
 
343
     &+ 2.d0*cone*kp(4)))
 
344
      
 
345
           yio(2,1) = yio(1,2)
 
346
           yio(3,1) = yio(1,3)
 
347
           yio(4,1) = yio(1,4)
 
348
            
 
349
           yio(3,2) = yio(2,3)
 
350
           yio(4,2) = yio(2,4)
 
351
           yio(4,3) = yio(3,4)
 
352
           
 
353
           do i = 1,4
 
354
              do j=1,4
 
355
                 yio(i,j) = yio(i,j)*d/2.0d0
 
356
              end do
 
357
           enddo            
 
358
                        
 
359
         else 
 
360
            write(*,*) 'invalid tensor mass'
 
361
            stop
 
362
      end if
 
363
 
 
364
      uio(1) = yio(1,1)
 
365
      uio(2) = yio(1,2)
 
366
      uio(3) = yio(1,3)
 
367
      uio(4) = yio(1,4)
 
368
      uio(5) = yio(2,1)
 
369
      uio(6) = yio(2,2)
 
370
      uio(7) = yio(2,3)
 
371
      uio(8) = yio(2,4)
 
372
      uio(9) = yio(3,1)
 
373
      uio(10) = yio(3,2)
 
374
      uio(11) = yio(3,3)
 
375
      uio(12) = yio(3,4)
 
376
      uio(13) = yio(4,1)
 
377
      uio(14) = yio(4,2)
 
378
      uio(15) = yio(4,3)
 
379
      uio(16) = yio(4,4)
 
380
      uio(17) = yio(5,1)
 
381
      uio(18) = yio(6,1)
 
382
 
 
383
      return
 
384
      end