1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
|
subroutine jvsshx(vc,sc1,sc2,g1,xm,xw , jvs)
c
c- by RF - Mar. 2006
c
c
c This subroutine computes an off-shell vector current from the vector-
c vector-Higgs-Higgs effective coupling.
c
implicit none
c-- dimension of the current set to arbitrary length
INTEGER DIM
PARAMETER(DIM=18)
c include "dimension.inc"
double complex vc(DIM),sc1(DIM),sc2(DIM),jvs(DIM)
double complex jvs1(DIM),jvs2(DIM),dg,qvc
double precision qp1,p12,p13,p14,p23,p24,p34
double precision p1(4),q(4),xm,xw,q2,vm2
double complex g1(2)
jvs(5) = vc(5)+sc1(2)+sc2(2)
jvs(6) = vc(6)+sc1(3)+sc2(3)
p1(1) = dble( vc(5))
p1(2) = dble( vc(6))
p1(3) = dimag(vc(6))
p1(4) = dimag(vc(5))
q(1) = -dble( jvs(5))
q(2) = -dble( jvs(6))
q(3) = -dimag(jvs(6))
q(4) = -dimag(jvs(5))
q2 = q(1)**2 - q(2)**2 - q(3)**2 - q(4)**2
jvs1(1) = (0D0,0D0)
jvs1(2) = (0D0,0D0)
jvs1(3) = (0D0,0D0)
jvs1(4) = (0D0,0D0)
jvs2(1) = (0D0,0D0)
jvs2(2) = (0D0,0D0)
jvs2(3) = (0D0,0D0)
jvs2(4) = (0D0,0D0)
dg = sc1(1)*sc2(1) /q2
if (g1(1).NE.(0D0,0D0)) then
qvc = vc(1)*q(1) - vc(2)*q(2) - vc(3)*q(3) - vc(4)*q(4)
qp1 = q(1)*p1(1) - q(2)*p1(2) - q(3)*p1(3) - q(4)*p1(4)
jvs1(1) = g1(1)* (vc(1)*qp1 - p1(1)*qvc)
jvs1(2) = g1(1)* (vc(2)*qp1 - p1(2)*qvc)
jvs1(3) = g1(1)* (vc(3)*qp1 - p1(3)*qvc)
jvs1(4) = g1(1)* (vc(4)*qp1 - p1(4)*qvc)
endif
if (g1(2).NE.(0D0,0D0)) then
p12 = p1(1)*q(2) - p1(2)*q(1)
p13 = p1(1)*q(3) - p1(3)*q(1)
p14 = p1(1)*q(4) - p1(4)*q(1)
p23 = p1(2)*q(3) - p1(3)*q(2)
p24 = p1(2)*q(4) - p1(4)*q(2)
p34 = p1(3)*q(4) - p1(4)*q(3)
jvs2(1)= - g1(2)* (-vc(2)*p34 +vc(3)*p24 -vc(4)*p23)
jvs2(2)= g1(2)* ( vc(1)*p34 -vc(3)*p14 +vc(4)*p13)
jvs2(3)= g1(2)* (-vc(1)*p24 +vc(2)*p14 -vc(4)*p12)
jvs2(4)= g1(2)* ( vc(1)*p23 -vc(2)*p13 +vc(3)*p12)
endif
jvs(1) = dg * (jvs1(1) + jvs2(1))
jvs(2) = dg * (jvs1(2) + jvs2(2))
jvs(3) = dg * (jvs1(3) + jvs2(3))
jvs(4) = dg * (jvs1(4) + jvs2(4))
return
end
|