~maddevelopers/mg5amcnlo/3.0.2-alpha0

« back to all changes in this revision

Viewing changes to HELAS/uvvcxx.F

Added Template and HELAS into bzr

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
      subroutine uvvcxx(v1,v2,gt,vmass,tmass,twidth , uvvh)
 
2
c
 
3
c This subroutine computes an off-shell tensor current from 
 
4
c the two gauge bosons-pseudo tensor boson coupling.
 
5
c
 
6
c input:
 
7
c       complex v1(3)          : first  vector                        v1
 
8
c       complex v2(3)          : second vector                        v2
 
9
c       real    gt             : coupling constant                gtv= gs
 
10
c       real    vmass          : vector boson mass                   m_v
 
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 uvv(18)        : tensor current         j^mu^nu(T:v1,v2)
 
16
c
 
17
      implicit none
 
18
      double complex v1(6), v2(6), uvvh(18)
 
19
      double precision vmass, tmass, twidth,gt
 
20
      integer i,j
 
21
      double complex yvv(6,4)
 
22
      double complex cZero
 
23
      double precision rZero, rTwo
 
24
      parameter( rZero = 0.0d0, rTwo = 2.0d0 )
 
25
      parameter( cZero = ( 0.0d0, 0.0d0 ) )
 
26
 
 
27
      
 
28
      yvv(5,1) = v1(5)+v2(5)
 
29
      yvv(6,1) = v1(6)+v2(6)
 
30
 
 
31
      do i=1,4
 
32
      do j=1,4
 
33
      yvv(i,j)=gt*v1(i)*v2(j)
 
34
      enddo
 
35
      enddo
 
36
 
 
37
      uvvh(1) = yvv(1,1)
 
38
      uvvh(2) = yvv(1,2)
 
39
      uvvh(3) = yvv(1,3)
 
40
      uvvh(4) = yvv(1,4)
 
41
      uvvh(5) = yvv(2,1)
 
42
      uvvh(6) = yvv(2,2)
 
43
      uvvh(7) = yvv(2,3)
 
44
      uvvh(8) = yvv(2,4)
 
45
      uvvh(9) = yvv(3,1)
 
46
      uvvh(10) = yvv(3,2)
 
47
      uvvh(11) = yvv(3,3)
 
48
      uvvh(12) = yvv(3,4)
 
49
      uvvh(13) = yvv(4,1)
 
50
      uvvh(14) = yvv(4,2)
 
51
      uvvh(15) = yvv(4,3)
 
52
      uvvh(16) = yvv(4,4)
 
53
      uvvh(17) = yvv(5,1)
 
54
      uvvh(18) = yvv(6,1)
 
55
 
 
56
      return
 
57
      end