← Back to branch summary

~maddevelopers/mg5amcnlo/2.5.3

« back to all changes in this revision

Viewing changes to HELAS/hvvtxx.F

  • Committer: olivier-mattelaer
  • Date: 2017-03-08 12:31:17 UTC
  • Revision ID: olivier-mattelaer-20170308123117-h0zkqjyh9sihsc61
empty version to have an effective freeze of the code

Show diffs side-by-side

added added

removed removed

Lines of Context:
HELAS/hvvtxx.F
1
 
      subroutine hvvtxx(w1,w2,g,xm,xw,jt)
2
 
c
3
 
c- by RF - Feb. 2006
4
 
c
5
 
c This subroutine computes the portion of the off-shell current
6
 
c for the color-octect tensor t3 in terms of w1 and w2 
7
 
c
8
 
c input:
9
 
c       complex w1(6)          : first  vector                        w1
10
 
c       complex w2(6)          : second vector                        w2
11
 
c       real    g              : first  coupling constant
12
 
c       real    xm             : not used
13
 
c       real    xw             : not used
14
 
c
15
 
c output:
16
 
c       complex jt(18)        : tensor current  j^(mu,nu)(w':w1,w2,w3)
17
 
c
18
 
      implicit none
19
 
 
20
 
c dimension of the current set to arbitrary length
21
 
      integer DIM
22
 
      parameter (DIM=18)
23
 
c      include "dimension.inc"
24
 
      double complex w1(DIM),w2(DIM),jt(DIM)
25
 
 
26
 
      double precision xm,xw,g,s2g
27
 
 
28
 
      double precision sqrTwo
29
 
      parameter( sqrTwo = 1.41421356237309514547462185873882845044d0 )
30
 
 
31
 
      s2g = g * sqrTwo
32
 
      
33
 
      jt( 1) = 0 ! dcmplx( s2g * (w1(1)*w2(1)-w1(1)*w2(1)) )
34
 
      jt( 2) =  dcmplx( s2g * (w1(1)*w2(2)-w1(2)*w2(1)) )
35
 
      jt( 3) =  dcmplx( s2g * (w1(1)*w2(3)-w1(3)*w2(1)) )
36
 
      jt( 4) =  dcmplx( s2g * (w1(1)*w2(4)-w1(4)*w2(1)) )
37
 
 
38
 
      jt( 5) =  dcmplx( s2g * (w1(2)*w2(1)-w1(1)*w2(2)) )
39
 
      jt( 6) = 0 ! dcmplx( s2g * (w1(2)*w2(2)-w1(2)*w2(2)) )
40
 
      jt( 7) =  dcmplx( s2g * (w1(2)*w2(3)-w1(3)*w2(2)) )
41
 
      jt( 8) =  dcmplx( s2g * (w1(2)*w2(4)-w1(4)*w2(2)) )
42
 
 
43
 
      jt( 9) =  dcmplx( s2g * (w1(3)*w2(1)-w1(1)*w2(3)) )
44
 
      jt(10) =  dcmplx( s2g * (w1(3)*w2(2)-w1(2)*w2(3)) )
45
 
      jt(11) = 0 ! dcmplx( s2g * (w1(3)*w2(3)-w1(3)*w2(3)) )
46
 
      jt(12) =  dcmplx( s2g * (w1(3)*w2(4)-w1(4)*w2(3)) )
47
 
 
48
 
      jt(13) =  dcmplx( s2g * (w1(4)*w2(1)-w1(1)*w2(4)) )
49
 
      jt(14) =  dcmplx( s2g * (w1(4)*w2(2)-w1(2)*w2(4)) )
50
 
      jt(15) =  dcmplx( s2g * (w1(4)*w2(3)-w1(3)*w2(4)) )
51
 
      jt(16) = 0 ! dcmplx( s2g * (w1(4)*w2(4)-w1(4)*w2(4)) )
52
 
 
53
 
      jt(17) = w1(5) + w2(5)
54
 
      jt(18) = w1(6) + w2(6)
55
 
 
56
 
      return
57
 
      end