← Back to branch summary

~maddevelopers/mg5amcnlo/3.0.2-alpha0

« back to all changes in this revision

Viewing changes to HELAS/vvtaxx.F

Added Template and HELAS into bzr

Show diffs side-by-side

added added

removed removed

Lines of Context:
HELAS/vvtaxx.F
 
1
      subroutine vvtaxx(ga,gb,tc,g,xm, vertex)
 
2
c
 
3
c- by RF - Feb. 2006
 
4
c
 
5
c This subroutine computes the portion of the amplitude of the four-point 
 
6
c coupling of 2 massless color octet gauge bosons
 
7
C with an effective color octect antisymmetric tensor.
 
8
c                                                                       
 
9
c input:
 
10
c     complex ga(6)                       : first vector  (gluon)
 
11
c     complex gb(6)                       : second vector (gluon)
 
12
c     complex tc(18)                      : tensor current
 
13
c     real    g                           : coupling constant
 
14
c 
 
15
c output:
 
16
c     complex vertex                      : amplitude
 
17
c
 
18
 
 
19
      implicit none
 
20
 
 
21
c dimension of the current set to arbitrary length
 
22
c      integer DIM
 
23
c      parameter (DIM=18)
 
24
      include "dimension.inc"
 
25
      double complex ga(DIM),gb(DIM),tc(DIM)
 
26
      double precision xm,xw,g
 
27
      double precision sqrTwo
 
28
      parameter( sqrTwo=1.41421356237309514547462185873882845044d0 )
 
29
 
 
30
      double complex dvertx, vertex
 
31
 
 
32
 
 
33
      dvertx = ! + tc( 1)*( ga(1)*gb(1) - ga(1)*gb(1) )  Always zero
 
34
     &          - tc( 2)*( ga(1)*gb(2) - ga(2)*gb(1) ) 
 
35
     &          - tc( 3)*( ga(1)*gb(3) - ga(3)*gb(1) ) 
 
36
     &          - tc( 4)*( ga(1)*gb(4) - ga(4)*gb(1) ) 
 
37
     
 
38
     &          - tc( 5)*( ga(2)*gb(1) - ga(1)*gb(2) ) 
 
39
     &         ! + tc( 6)*( ga(2)*gb(2) - ga(2)*gb(2) )  Always zero
 
40
     &          + tc( 7)*( ga(2)*gb(3) - ga(3)*gb(2) ) 
 
41
     &          + tc( 8)*( ga(2)*gb(4) - ga(4)*gb(2) ) 
 
42
     
 
43
     &          - tc( 9)*( ga(3)*gb(1) - ga(1)*gb(3) ) 
 
44
     &          + tc(10)*( ga(3)*gb(2) - ga(2)*gb(3) ) 
 
45
     &         ! + tc(11)*( ga(3)*gb(3) - ga(3)*gb(3) )  Always zero
 
46
     &          + tc(12)*( ga(3)*gb(4) - ga(4)*gb(3) ) 
 
47
     
 
48
     &          - tc(13)*( ga(4)*gb(1) - ga(1)*gb(4) ) 
 
49
     &          + tc(14)*( ga(4)*gb(2) - ga(2)*gb(4) ) 
 
50
     &          + tc(15)*( ga(4)*gb(3) - ga(3)*gb(4) ) 
 
51
     &         ! + tc(16)*( ga(4)*gb(4) - ga(4)*gb(4) )  Always zero
 
52
 
 
53
      vertex = g * dvertx /sqrTwo
 
54
 
 
55
 
 
56
      return
 
57
      end