1
subroutine vvsshx(v1,v2,sc1,sc2,g1, vertex)
6
c This subroutine computes an amplitude of the vector-vector-Higgs-Higgs
10
c complex v1(6) : first vector
11
c complex v2(6) : second vector
12
c complex sc1(3) : first scalar
13
c complex sc2(3) : second scalar
14
c complex g1(2) : first coupling constant
17
c complex vertex : amplitude gamma(v1,v2,s,s)
20
c-- dimension of the current set to arbitrary length
23
include "dimension.inc"
24
double complex v1(DIM),v2(DIM),sc1(DIM),sc2(DIM),vertex
25
double complex vertex1,vertex2,g1(2)
26
double complex v12,p2v1,p1v2,v13,v14,v23,v24,v34
27
double precision p12,p13,p14,p23,p24,p34
28
double precision p1(0:3),p2(0:3)
44
if (g1(1).NE.(0D0,0D0)) then
46
v12 = v1(1)*v2(1) - v1(2)*v2(2) - v1(3)*v2(3) - v1(4)*v2(4)
47
p12 = p1(0)*p2(0) - p1(1)*p2(1) - p1(2)*p2(2) - p1(3)*p2(3)
48
p2v1 = v1(1)*p2(0) - v1(2)*p2(1) - v1(3)*p2(2) - v1(4)*p2(3)
49
p1v2 = p1(0)*v2(1) - p1(1)*v2(2) - p1(2)*v2(3) - p1(3)*v2(4)
51
vertex1 = g1(1) *(v12*p12 - p2v1*p1v2)
54
if (g1(2).NE.(0D0,0D0)) then
55
p12 = p1(0)*p2(1) - p1(1)*p2(0)
56
p13 = p1(0)*p2(2) - p1(2)*p2(0)
57
p14 = p1(0)*p2(3) - p1(3)*p2(0)
58
p23 = p1(1)*p2(2) - p1(2)*p2(1)
59
p24 = p1(1)*p2(3) - p1(3)*p2(1)
60
p34 = p1(2)*p2(3) - p1(3)*p2(2)
62
v12 = v1(1)*v2(2) - v1(2)*v2(1)
63
v13 = v1(1)*v2(3) - v1(3)*v2(1)
64
v14 = v1(1)*v2(4) - v1(4)*v2(1)
65
v23 = v1(2)*v2(3) - v1(3)*v2(2)
66
v24 = v1(2)*v2(4) - v1(4)*v2(2)
67
v34 = v1(3)*v2(4) - v1(4)*v2(3)
69
vertex2 = - g1(2)*( v12*p34 - v13*p24 + v14*p23
70
& +v23*p14 - v24*p13 + v34*p12 )
73
vertex = sc1(1)*sc2(1)*(vertex1 + vertex2)