← Back to branch summary

~maddevelopers/mg5amcnlo/3.0.2-alpha0

« back to all changes in this revision

Viewing changes to HELAS/tttxxx.F

Added Template and HELAS into bzr

Show diffs side-by-side

added added

removed removed

Lines of Context:
HELAS/tttxxx.F
 
1
      subroutine tttxxx(tt,ta,tb,gt,vertex)
 
2
c
 
3
c This subroutine computes an amplitude of the three-point coupling of
 
4
c two gauge bosons and a tensor boson.
 
5
c
 
6
c input:
 
7
c       complex ta(18)          : first  tensor                           ta
 
8
c       complex tb(18)          : second tensor                       tb
 
9
c       complex tt(18)           : input  tensor                         tt
 
10
c       complex gt             : coupling constant         gtv=-1/Lambda
 
11
c output:
 
12
c       complex vertex         : amplitude                gamma(tt,ta,tb)
 
13
c
 
14
      implicit none
 
15
      double complex ta(18), tb(18), tt(18), gt, vertex
 
16
      double complex T1(6,4),T2(6,4),T(6,4)
 
17
 
 
18
      integer i, j
 
19
 
 
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
      T(1,1) = tt(1)
 
27
      T(1,2) = tt(2)
 
28
      T(1,3) = tt(3)
 
29
      T(1,4) = tt(4)
 
30
      T(2,1) = tt(5)
 
31
      T(2,2) = tt(6)
 
32
      T(2,3) = tt(7)
 
33
      T(2,4) = tt(8)
 
34
      T(3,1) = tt(9)
 
35
      T(3,2) = tt(10)
 
36
      T(3,3) = tt(11)
 
37
      T(3,4) = tt(12)
 
38
      T(4,1) = tt(13)
 
39
      T(4,2) = tt(14)
 
40
      T(4,3) = tt(15)
 
41
      T(4,4) = tt(16)
 
42
      T(5,1) = tt(17)
 
43
      T(6,1) = tt(18)
 
44
 
 
45
      T1(1,1) = ta(1)
 
46
      T1(1,2) = ta(2)
 
47
      T1(1,3) = ta(3)
 
48
      T1(1,4) = ta(4)
 
49
      T1(2,1) = ta(5)
 
50
      T1(2,2) = ta(6)
 
51
      T1(2,3) = ta(7)
 
52
      T1(2,4) = ta(8)
 
53
      T1(3,1) = ta(9)
 
54
      T1(3,2) = ta(10)
 
55
      T1(3,3) = ta(11)
 
56
      T1(3,4) = ta(12)
 
57
      T1(4,1) = ta(13)
 
58
      T1(4,2) = ta(14)
 
59
      T1(4,3) = ta(15)
 
60
      T1(4,4) = ta(16)
 
61
      T1(5,1) = ta(17)
 
62
      T1(6,1) = ta(18)
 
63
 
 
64
      T2(1,1) = tb(1)
 
65
      T2(1,2) = tb(2)
 
66
      T2(1,3) = tb(3)
 
67
      T2(1,4) = tb(4)
 
68
      T2(2,1) = tb(5)
 
69
      T2(2,2) = tb(6)
 
70
      T2(2,3) = tb(7)
 
71
      T2(2,4) = tb(8)
 
72
      T2(3,1) = tb(9)
 
73
      T2(3,2) = tb(10)
 
74
      T2(3,3) = tb(11)
 
75
      T2(3,4) = tb(12)
 
76
      T2(4,1) = tb(13)
 
77
      T2(4,2) = tb(14)
 
78
      T2(4,3) = tb(15)
 
79
      T2(4,4) = tb(16)
 
80
      T2(5,1) = tb(17)
 
81
      T2(6,1) = tb(18)
 
82
 
 
83
 
 
84
      vertex =-gt*(T(1,1)*T1(1,2)*T2(1,2) - T(2,2)*T1(1,2)*T2(1,2) + 
 
85
     -  T(3,3)*T1(1,2)*T2(1,2) + T(4,4)*T1(1,2)*T2(1,2) - 
 
86
     -  T(2,3)*T1(1,3)*T2(1,2) - T(3,2)*T1(1,3)*T2(1,2) - 
 
87
     -  T(2,4)*T1(1,4)*T2(1,2) - T(4,2)*T1(1,4)*T2(1,2) - 
 
88
     -  T(1,1)*T1(2,1)*T2(1,2) + T(2,2)*T1(2,1)*T2(1,2) - 
 
89
     -  T(3,3)*T1(2,1)*T2(1,2) - T(4,4)*T1(2,1)*T2(1,2) + 
 
90
     -  T(1,3)*T1(2,3)*T2(1,2) + T(3,1)*T1(2,3)*T2(1,2) + 
 
91
     -  T(1,4)*T1(2,4)*T2(1,2) + T(4,1)*T1(2,4)*T2(1,2) + 
 
92
     -  T(2,3)*T1(3,1)*T2(1,2) + T(3,2)*T1(3,1)*T2(1,2) - 
 
93
     -  T(1,3)*T1(3,2)*T2(1,2) - T(3,1)*T1(3,2)*T2(1,2) + 
 
94
     -  T(2,4)*T1(4,1)*T2(1,2) + T(4,2)*T1(4,1)*T2(1,2) - 
 
95
     -  T(1,4)*T1(4,2)*T2(1,2) - T(4,1)*T1(4,2)*T2(1,2) - 
 
96
     -  T(2,3)*T1(1,2)*T2(1,3) - T(3,2)*T1(1,2)*T2(1,3) + 
 
97
     -  T(1,1)*T1(1,3)*T2(1,3) + T(2,2)*T1(1,3)*T2(1,3) - 
 
98
     -  T(3,3)*T1(1,3)*T2(1,3) + T(4,4)*T1(1,3)*T2(1,3) - 
 
99
     -  T(3,4)*T1(1,4)*T2(1,3) - T(4,3)*T1(1,4)*T2(1,3) + 
 
100
     -  T(2,3)*T1(2,1)*T2(1,3) + T(3,2)*T1(2,1)*T2(1,3) - 
 
101
     -  T(1,2)*T1(2,3)*T2(1,3) - T(2,1)*T1(2,3)*T2(1,3) - 
 
102
     -  T(1,1)*T1(3,1)*T2(1,3) - T(2,2)*T1(3,1)*T2(1,3) + 
 
103
     -  T(3,3)*T1(3,1)*T2(1,3) - T(4,4)*T1(3,1)*T2(1,3) + 
 
104
     -  T(1,2)*T1(3,2)*T2(1,3) + T(2,1)*T1(3,2)*T2(1,3) + 
 
105
     -  T(1,4)*T1(3,4)*T2(1,3) + T(4,1)*T1(3,4)*T2(1,3) + 
 
106
     -  T(3,4)*T1(4,1)*T2(1,3) + T(4,3)*T1(4,1)*T2(1,3) - 
 
107
     -  T(1,4)*T1(4,3)*T2(1,3) - T(4,1)*T1(4,3)*T2(1,3) - 
 
108
     -  T(2,4)*T1(1,2)*T2(1,4) - T(4,2)*T1(1,2)*T2(1,4) - 
 
109
     -  T(3,4)*T1(1,3)*T2(1,4) - T(4,3)*T1(1,3)*T2(1,4) + 
 
110
     -  T(1,1)*T1(1,4)*T2(1,4) + T(2,2)*T1(1,4)*T2(1,4) + 
 
111
     -  T(3,3)*T1(1,4)*T2(1,4) - T(4,4)*T1(1,4)*T2(1,4) + 
 
112
     -  T(2,4)*T1(2,1)*T2(1,4) + T(4,2)*T1(2,1)*T2(1,4) - 
 
113
     -  T(1,2)*T1(2,4)*T2(1,4) - T(2,1)*T1(2,4)*T2(1,4) + 
 
114
     -  T(3,4)*T1(3,1)*T2(1,4) + T(4,3)*T1(3,1)*T2(1,4) - 
 
115
     -  T(1,3)*T1(3,4)*T2(1,4) - T(3,1)*T1(3,4)*T2(1,4) - 
 
116
     -  T(1,1)*T1(4,1)*T2(1,4) - T(2,2)*T1(4,1)*T2(1,4) - 
 
117
     -  T(3,3)*T1(4,1)*T2(1,4) + T(4,4)*T1(4,1)*T2(1,4) + 
 
118
     -  T(1,2)*T1(4,2)*T2(1,4) + T(2,1)*T1(4,2)*T2(1,4) + 
 
119
     -  T(1,3)*T1(4,3)*T2(1,4) + T(3,1)*T1(4,3)*T2(1,4) - 
 
120
     -  T(1,1)*T1(1,2)*T2(2,1) + T(2,2)*T1(1,2)*T2(2,1) - 
 
121
     -  T(3,3)*T1(1,2)*T2(2,1) - T(4,4)*T1(1,2)*T2(2,1) + 
 
122
     -  T(2,3)*T1(1,3)*T2(2,1) + T(3,2)*T1(1,3)*T2(2,1) + 
 
123
     -  T(2,4)*T1(1,4)*T2(2,1) + T(4,2)*T1(1,4)*T2(2,1) + 
 
124
     -  T(1,1)*T1(2,1)*T2(2,1) - T(2,2)*T1(2,1)*T2(2,1) + 
 
125
     -  T(3,3)*T1(2,1)*T2(2,1) + T(4,4)*T1(2,1)*T2(2,1) - 
 
126
     -  T(1,3)*T1(2,3)*T2(2,1) - T(3,1)*T1(2,3)*T2(2,1) - 
 
127
     -  T(1,4)*T1(2,4)*T2(2,1) - T(4,1)*T1(2,4)*T2(2,1) - 
 
128
     -  T(2,3)*T1(3,1)*T2(2,1) - T(3,2)*T1(3,1)*T2(2,1) + 
 
129
     -  T(1,3)*T1(3,2)*T2(2,1) + T(3,1)*T1(3,2)*T2(2,1) - 
 
130
     -  T(2,4)*T1(4,1)*T2(2,1) - T(4,2)*T1(4,1)*T2(2,1) + 
 
131
     -  T(1,4)*T1(4,2)*T2(2,1) + T(4,1)*T1(4,2)*T2(2,1) + 
 
132
     -  T(1,3)*T1(1,2)*T2(2,3) + T(3,1)*T1(1,2)*T2(2,3) - 
 
133
     -  T(1,2)*T1(1,3)*T2(2,3) - T(2,1)*T1(1,3)*T2(2,3) - 
 
134
     -  T(1,3)*T1(2,1)*T2(2,3) - T(3,1)*T1(2,1)*T2(2,3) + 
 
135
     -  T(1,1)*T1(2,3)*T2(2,3) + T(2,2)*T1(2,3)*T2(2,3) + 
 
136
     -  T(3,3)*T1(2,3)*T2(2,3) - T(4,4)*T1(2,3)*T2(2,3) + 
 
137
     -  T(3,4)*T1(2,4)*T2(2,3) + T(4,3)*T1(2,4)*T2(2,3) + 
 
138
     -  T(1,2)*T1(3,1)*T2(2,3) + T(2,1)*T1(3,1)*T2(2,3) - 
 
139
     -  T(1,1)*T1(3,2)*T2(2,3) - T(2,2)*T1(3,2)*T2(2,3) - 
 
140
     -  T(3,3)*T1(3,2)*T2(2,3) + T(4,4)*T1(3,2)*T2(2,3) - 
 
141
     -  T(2,4)*T1(3,4)*T2(2,3) - T(4,2)*T1(3,4)*T2(2,3) - 
 
142
     -  T(3,4)*T1(4,2)*T2(2,3) - T(4,3)*T1(4,2)*T2(2,3) + 
 
143
     -  T(2,4)*T1(4,3)*T2(2,3) + T(4,2)*T1(4,3)*T2(2,3) + 
 
144
     -  T(1,4)*T1(1,2)*T2(2,4) + T(4,1)*T1(1,2)*T2(2,4) - 
 
145
     -  T(1,2)*T1(1,4)*T2(2,4) - T(2,1)*T1(1,4)*T2(2,4) - 
 
146
     -  T(1,4)*T1(2,1)*T2(2,4) - T(4,1)*T1(2,1)*T2(2,4) + 
 
147
     -  T(3,4)*T1(2,3)*T2(2,4) + T(4,3)*T1(2,3)*T2(2,4) + 
 
148
     -  T(1,1)*T1(2,4)*T2(2,4) + T(2,2)*T1(2,4)*T2(2,4) - 
 
149
     -  T(3,3)*T1(2,4)*T2(2,4) + T(4,4)*T1(2,4)*T2(2,4) - 
 
150
     -  T(3,4)*T1(3,2)*T2(2,4) - T(4,3)*T1(3,2)*T2(2,4) + 
 
151
     -  T(2,3)*T1(3,4)*T2(2,4) + T(3,2)*T1(3,4)*T2(2,4) + 
 
152
     -  T(1,2)*T1(4,1)*T2(2,4) + T(2,1)*T1(4,1)*T2(2,4) - 
 
153
     -  T(1,1)*T1(4,2)*T2(2,4) - T(2,2)*T1(4,2)*T2(2,4) + 
 
154
     -  T(3,3)*T1(4,2)*T2(2,4) - T(4,4)*T1(4,2)*T2(2,4) - 
 
155
     -  T(2,3)*T1(4,3)*T2(2,4) - T(3,2)*T1(4,3)*T2(2,4) + 
 
156
     -  T(2,3)*T1(1,2)*T2(3,1) + T(3,2)*T1(1,2)*T2(3,1) - 
 
157
     -  T(1,1)*T1(1,3)*T2(3,1) - T(2,2)*T1(1,3)*T2(3,1) + 
 
158
     -  T(3,3)*T1(1,3)*T2(3,1) - T(4,4)*T1(1,3)*T2(3,1) + 
 
159
     -  T(3,4)*T1(1,4)*T2(3,1) + T(4,3)*T1(1,4)*T2(3,1) - 
 
160
     -  T(2,3)*T1(2,1)*T2(3,1) - T(3,2)*T1(2,1)*T2(3,1) + 
 
161
     -  T(1,2)*T1(2,3)*T2(3,1) + T(2,1)*T1(2,3)*T2(3,1) + 
 
162
     -  T(1,1)*T1(3,1)*T2(3,1) + T(2,2)*T1(3,1)*T2(3,1) - 
 
163
     -  T(3,3)*T1(3,1)*T2(3,1) + T(4,4)*T1(3,1)*T2(3,1) - 
 
164
     -  T(1,2)*T1(3,2)*T2(3,1) - T(2,1)*T1(3,2)*T2(3,1) - 
 
165
     -  T(1,4)*T1(3,4)*T2(3,1) - T(4,1)*T1(3,4)*T2(3,1) - 
 
166
     -  T(3,4)*T1(4,1)*T2(3,1) - T(4,3)*T1(4,1)*T2(3,1) + 
 
167
     -  T(1,4)*T1(4,3)*T2(3,1) + T(4,1)*T1(4,3)*T2(3,1) - 
 
168
     -  T(1,3)*T1(1,2)*T2(3,2) - T(3,1)*T1(1,2)*T2(3,2) + 
 
169
     -  T(1,2)*T1(1,3)*T2(3,2) + T(2,1)*T1(1,3)*T2(3,2) + 
 
170
     -  T(1,3)*T1(2,1)*T2(3,2) + T(3,1)*T1(2,1)*T2(3,2) - 
 
171
     -  T(1,1)*T1(2,3)*T2(3,2) - T(2,2)*T1(2,3)*T2(3,2) - 
 
172
     -  T(3,3)*T1(2,3)*T2(3,2) + T(4,4)*T1(2,3)*T2(3,2) - 
 
173
     -  T(3,4)*T1(2,4)*T2(3,2) - T(4,3)*T1(2,4)*T2(3,2) - 
 
174
     -  T(1,2)*T1(3,1)*T2(3,2) - T(2,1)*T1(3,1)*T2(3,2) + 
 
175
     -  T(1,1)*T1(3,2)*T2(3,2) + T(2,2)*T1(3,2)*T2(3,2) + 
 
176
     -  T(3,3)*T1(3,2)*T2(3,2) - T(4,4)*T1(3,2)*T2(3,2) + 
 
177
     -  T(2,4)*T1(3,4)*T2(3,2) + T(4,2)*T1(3,4)*T2(3,2) + 
 
178
     -  T(3,4)*T1(4,2)*T2(3,2) + T(4,3)*T1(4,2)*T2(3,2) - 
 
179
     -  T(2,4)*T1(4,3)*T2(3,2) - T(4,2)*T1(4,3)*T2(3,2) + 
 
180
     -  T(1,4)*T1(1,3)*T2(3,4) + T(4,1)*T1(1,3)*T2(3,4) - 
 
181
     -  T(1,3)*T1(1,4)*T2(3,4) - T(3,1)*T1(1,4)*T2(3,4) - 
 
182
     -  T(2,4)*T1(2,3)*T2(3,4) - T(4,2)*T1(2,3)*T2(3,4) + 
 
183
     -  T(2,3)*T1(2,4)*T2(3,4) + T(3,2)*T1(2,4)*T2(3,4) - 
 
184
     -  T(1,4)*T1(3,1)*T2(3,4) - T(4,1)*T1(3,1)*T2(3,4) + 
 
185
     -  T(2,4)*T1(3,2)*T2(3,4) + T(4,2)*T1(3,2)*T2(3,4) + 
 
186
     -  T(1,1)*T1(3,4)*T2(3,4) - T(2,2)*T1(3,4)*T2(3,4) + 
 
187
     -  T(3,3)*T1(3,4)*T2(3,4) + T(4,4)*T1(3,4)*T2(3,4) + 
 
188
     -  T(1,3)*T1(4,1)*T2(3,4) + T(3,1)*T1(4,1)*T2(3,4) - 
 
189
     -  T(2,3)*T1(4,2)*T2(3,4) - T(3,2)*T1(4,2)*T2(3,4) - 
 
190
     -  T(1,1)*T1(4,3)*T2(3,4) + T(2,2)*T1(4,3)*T2(3,4) - 
 
191
     -  T(3,3)*T1(4,3)*T2(3,4) - T(4,4)*T1(4,3)*T2(3,4) + 
 
192
     -  T(2,4)*T1(1,2)*T2(4,1) + T(4,2)*T1(1,2)*T2(4,1) + 
 
193
     -  T(3,4)*T1(1,3)*T2(4,1) + T(4,3)*T1(1,3)*T2(4,1) - 
 
194
     -  T(1,1)*T1(1,4)*T2(4,1) - T(2,2)*T1(1,4)*T2(4,1) - 
 
195
     -  T(3,3)*T1(1,4)*T2(4,1) + T(4,4)*T1(1,4)*T2(4,1) - 
 
196
     -  T(2,4)*T1(2,1)*T2(4,1) - T(4,2)*T1(2,1)*T2(4,1) + 
 
197
     -  T(1,2)*T1(2,4)*T2(4,1) + T(2,1)*T1(2,4)*T2(4,1) - 
 
198
     -  T(3,4)*T1(3,1)*T2(4,1) - T(4,3)*T1(3,1)*T2(4,1) + 
 
199
     -  T(1,3)*T1(3,4)*T2(4,1) + T(3,1)*T1(3,4)*T2(4,1) + 
 
200
     -  T(1,1)*T1(4,1)*T2(4,1) + T(2,2)*T1(4,1)*T2(4,1) + 
 
201
     -  T(3,3)*T1(4,1)*T2(4,1) - T(4,4)*T1(4,1)*T2(4,1) - 
 
202
     -  T(1,2)*T1(4,2)*T2(4,1) - T(2,1)*T1(4,2)*T2(4,1) - 
 
203
     -  T(1,3)*T1(4,3)*T2(4,1) - T(3,1)*T1(4,3)*T2(4,1) - 
 
204
     -  T(1,4)*T1(1,2)*T2(4,2) - T(4,1)*T1(1,2)*T2(4,2) + 
 
205
     -  T(1,2)*T1(1,4)*T2(4,2) + T(2,1)*T1(1,4)*T2(4,2) + 
 
206
     -  T(1,4)*T1(2,1)*T2(4,2) + T(4,1)*T1(2,1)*T2(4,2) - 
 
207
     -  T(3,4)*T1(2,3)*T2(4,2) - T(4,3)*T1(2,3)*T2(4,2) - 
 
208
     -  T(1,1)*T1(2,4)*T2(4,2) - T(2,2)*T1(2,4)*T2(4,2) + 
 
209
     -  T(3,3)*T1(2,4)*T2(4,2) - T(4,4)*T1(2,4)*T2(4,2) + 
 
210
     -  T(3,4)*T1(3,2)*T2(4,2) + T(4,3)*T1(3,2)*T2(4,2) - 
 
211
     -  T(2,3)*T1(3,4)*T2(4,2) - T(3,2)*T1(3,4)*T2(4,2) - 
 
212
     -  T(1,2)*T1(4,1)*T2(4,2) - T(2,1)*T1(4,1)*T2(4,2) + 
 
213
     -  T(1,1)*T1(4,2)*T2(4,2) + T(2,2)*T1(4,2)*T2(4,2) - 
 
214
     -  T(3,3)*T1(4,2)*T2(4,2) + T(4,4)*T1(4,2)*T2(4,2) + 
 
215
     -  T(2,3)*T1(4,3)*T2(4,2) + T(3,2)*T1(4,3)*T2(4,2) - 
 
216
     -  T(1,4)*T1(1,3)*T2(4,3) - T(4,1)*T1(1,3)*T2(4,3) + 
 
217
     -  T(1,3)*T1(1,4)*T2(4,3) + T(3,1)*T1(1,4)*T2(4,3) + 
 
218
     -  T(2,4)*T1(2,3)*T2(4,3) + T(4,2)*T1(2,3)*T2(4,3) - 
 
219
     -  T(2,3)*T1(2,4)*T2(4,3) - T(3,2)*T1(2,4)*T2(4,3) + 
 
220
     -  T(1,4)*T1(3,1)*T2(4,3) + T(4,1)*T1(3,1)*T2(4,3) - 
 
221
     -  T(2,4)*T1(3,2)*T2(4,3) - T(4,2)*T1(3,2)*T2(4,3) - 
 
222
     -  T(1,1)*T1(3,4)*T2(4,3) + T(2,2)*T1(3,4)*T2(4,3) - 
 
223
     -  T(3,3)*T1(3,4)*T2(4,3) - T(4,4)*T1(3,4)*T2(4,3) - 
 
224
     -  T(1,3)*T1(4,1)*T2(4,3) - T(3,1)*T1(4,1)*T2(4,3) + 
 
225
     -  T(2,3)*T1(4,2)*T2(4,3) + T(3,2)*T1(4,2)*T2(4,3) + 
 
226
     -  T(1,1)*T1(4,3)*T2(4,3) - T(2,2)*T1(4,3)*T2(4,3) + 
 
227
     -  T(3,3)*T1(4,3)*T2(4,3) + T(4,4)*T1(4,3)*T2(4,3))
 
228
 
 
229
       return
 
230
       end
 
231