← Back to branch summary

~maddevelopers/mg5amcnlo/2.5.3

« back to all changes in this revision

Viewing changes to HELAS/ggggtx.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/ggggtx.F
1
 
      subroutine ggggtx(va,vb,vc,vd,tc,gc,gt , vertex)
2
 
c      
3
 
c This subroutine computes the portion of the amplitude of the five-point 
4
 
c coupling of a tensor boson with 4 massless color octet gauge bosons
5
 
c (gluons) corresponding to the color structure f^{a,b,e} f{c,d,e}. 
6
 
c
7
 
c To obtain the complete amplitude, this coupling must be called three
8
 
c times (once for each color structure) with the following permutations:
9
 
c     call ggggtx(va,vb,vc,vd,tc,gc,gt , vertex1)
10
 
c     call ggggtx(va,vc,vd,vb,tc,gc,gt , vertex2)
11
 
c     call ggggtx(va,vd,vb,vc,tc,gc,gt , vertex3)
12
 
c corresponding to
13
 
c       f^{a,b,e} f^{c,d,e}
14
 
c       f^{a,c,e} f^{d,b,e}
15
 
c       f^{a,d,e} f^{b,c,e}
16
 
c                                                                       
17
 
c input:                                                                
18
 
c       complex va(6)          : boson with adjoint color index a     va
19
 
c       complex vb(6)          : boson with adjoint color index b     vb
20
 
c       complex vc(6)          : boson with adjoint color index c     vc
21
 
c       complex vd(6)          : boson with adjoint color index d     vd
22
 
c       complex tc(18)         : input tensor                          T
23
 
c       real    gc             : coupling constant                    gs
24
 
c       complex gt             : coupling constant         gtv=-1/Lambda
25
 
c
26
 
c output:
27
 
c       complex vertex         : amplitude          gamma(va,vb,vc,vd,T)
28
 
c
29
 
c- by Q.Li - OCT. 2006
30
 
c
31
 
      implicit none
32
 
      double complex va(6), vb(6), vc(6), vd(6), tc(18), gt, vertex
33
 
      double precision gc
34
 
 
35
 
      double complex vab,vac,vad,vbc,vbd,vcd,ft(6,4),dvertx
36
 
      double complex T00, T12, T13, T14, T23, T24, T34  
37
 
      double complex TV24,TV23,TV14,TV13
38
 
      double complex cZero
39
 
      double precision rZero, rTwo
40
 
      parameter( rZero = 0.0d0, rTwo = 2.0d0 )
41
 
      parameter( cZero = ( 0.0d0, 0.0d0 ) )
42
 
 
43
 
      
44
 
      ft(1,1) = tc(1)
45
 
      ft(1,2) = tc(2)
46
 
      ft(1,3) = tc(3)
47
 
      ft(1,4) = tc(4)
48
 
      ft(2,1) = tc(5)
49
 
      ft(2,2) = tc(6)
50
 
      ft(2,3) = tc(7)
51
 
      ft(2,4) = tc(8)
52
 
      ft(3,1) = tc(9)
53
 
      ft(3,2) = tc(10)
54
 
      ft(3,3) = tc(11)
55
 
      ft(3,4) = tc(12)
56
 
      ft(4,1) = tc(13)
57
 
      ft(4,2) = tc(14)
58
 
      ft(4,3) = tc(15)
59
 
      ft(4,4) = tc(16)
60
 
      ft(5,1) = tc(17)
61
 
      ft(6,1) = tc(18)
62
 
 
63
 
      T00 = ft(1,1)-ft(2,2)-ft(3,3)-ft(4,4)
64
 
      T12 = ft(1,2) + ft(2,1)
65
 
      T13 = ft(1,3) + ft(3,1)
66
 
      T14 = ft(1,4) + ft(4,1)
67
 
      T23 = ft(2,3) + ft(3,2)
68
 
      T24 = ft(2,4) + ft(4,2)
69
 
      T34 = ft(3,4) + ft(4,3)
70
 
 
71
 
      TV14 = rtwo*(ft(1,1)*va(1)*vd(1)+ft(2,2)*va(2)*vd(2)
72
 
     &+ft(3,3)*va(3)*vd(3)+ft(4,4)*va(4)*vd(4))
73
 
 
74
 
      TV13 = rtwo*(ft(1,1)*va(1)*vc(1)+ft(2,2)*va(2)*vc(2)
75
 
     &+ft(3,3)*va(3)*vc(3)+ft(4,4)*va(4)*vc(4))
76
 
 
77
 
      TV24 = rtwo*(ft(1,1)*vb(1)*vd(1)+ft(2,2)*vb(2)*vd(2)
78
 
     &+ft(3,3)*vb(3)*vd(3)+ft(4,4)*vb(4)*vd(4))
79
 
 
80
 
      TV23 = rtwo*(ft(1,1)*vb(1)*vc(1)+ft(2,2)*vb(2)*vc(2)
81
 
     &+ft(3,3)*vb(3)*vc(3)+ft(4,4)*vb(4)*vc(4))
82
 
 
83
 
        
84
 
      TV14 = TV14- T12*(va(1)*vd(2) + va(2)*vd(1))
85
 
     &          - T13*(va(1)*vd(3) + va(3)*vd(1))
86
 
     &          - T14*(va(1)*vd(4) + va(4)*vd(1))
87
 
     &          + T23*(va(2)*vd(3) + va(3)*vd(2))
88
 
     &          + T24*(va(2)*vd(4) + va(4)*vd(2))
89
 
     &          + T34*(va(3)*vd(4) + va(4)*vd(3))
90
 
      
91
 
      TV13 = TV13 - T12*(va(1)*vc(2) + va(2)*vc(1))
92
 
     &          - T13*(va(1)*vc(3) + va(3)*vc(1))
93
 
     &          - T14*(va(1)*vc(4) + va(4)*vc(1))
94
 
     &          + T23*(va(2)*vc(3) + va(3)*vc(2))
95
 
     &          + T24*(va(2)*vc(4) + va(4)*vc(2))
96
 
     &          + T34*(va(3)*vc(4) + va(4)*vc(3))
97
 
 
98
 
      TV24 = TV24 - T12*(vb(1)*vd(2) + vb(2)*vd(1))
99
 
     &          - T13*(vb(1)*vd(3) + vb(3)*vd(1))
100
 
     &          - T14*(vb(1)*vd(4) + vb(4)*vd(1))
101
 
     &          + T23*(vb(2)*vd(3) + vb(3)*vd(2))
102
 
     &          + T24*(vb(2)*vd(4) + vb(4)*vd(2))
103
 
     &          + T34*(vb(3)*vd(4) + vb(4)*vd(3))
104
 
 
105
 
      TV23 = TV23 - T12*(vb(1)*vc(2) + vb(2)*vc(1))
106
 
     &          - T13*(vb(1)*vc(3) + vb(3)*vc(1))
107
 
     &          - T14*(vb(1)*vc(4) + vb(4)*vc(1))
108
 
     &          + T23*(vb(2)*vc(3) + vb(3)*vc(2))
109
 
     &          + T24*(vb(2)*vc(4) + vb(4)*vc(2))
110
 
     &          + T34*(vb(3)*vc(4) + vb(4)*vc(3))
111
 
        
112
 
 
113
 
      vab = va(1)*vb(1)-va(2)*vb(2)-va(3)*vb(3)-va(4)*vb(4)
114
 
      vac = va(1)*vc(1)-va(2)*vc(2)-va(3)*vc(3)-va(4)*vc(4)
115
 
      vad = va(1)*vd(1)-va(2)*vd(2)-va(3)*vd(3)-va(4)*vd(4)
116
 
      vbc = vb(1)*vc(1)-vb(2)*vc(2)-vb(3)*vc(3)-vb(4)*vc(4)
117
 
      vbd = vb(1)*vd(1)-vb(2)*vd(2)-vb(3)*vd(3)-vb(4)*vd(4)
118
 
      vcd = vc(1)*vd(1)-vc(2)*vd(2)-vc(3)*vd(3)-vc(4)*vd(4)
119
 
 
120
 
      dvertx = -TV13*vbd-TV24*vac+TV23*vad+TV14*vbc 
121
 
     &+vbd*vac*T00-vad*vbc*T00
122
 
 
123
 
      vertex = -dvertx * gc*gc*gt
124
 
 
125
 
      return
126
 
      end