← Back to branch summary

~maddevelopers/mg5amcnlo/3.0.2-alpha0

« back to all changes in this revision

Viewing changes to HELAS/ttsxxx.F

Added Template and HELAS into bzr

Show diffs side-by-side

added added

removed removed

Lines of Context:
HELAS/ttsxxx.F
 
1
      subroutine ttsxxx(tc1,tc2,sc,gc,vertex)
 
2
c
 
3
c- by RF - Feb. 2006 
 
4
c
 
5
c This subroutine computes an amplitude of the tts coupling.
 
6
c
 
7
c     input:
 
8
c          complex tc1(18)           : Incoming tensor particle
 
9
c          complex tc2(18)           : Incoming tensor particle
 
10
c          complex sc(3)             : Incoming scalar particle (Higgs)
 
11
c          complex gc(2)             : coupling constant: gt(1) scalar
 
12
c                                                         gt(2) not used
 
13
c
 
14
c     output:
 
15
c          complex vertex            : amplitude for a tts vertex
 
16
c
 
17
 
 
18
      implicit none
 
19
c--   dimension of the current set to arbitrary length
 
20
c      INTEGER DIM
 
21
c      PARAMETER(DIM=18)
 
22
      include "dimension.inc"
 
23
      double complex tc1(DIM),tc2(DIM),sc(DIM)
 
24
      double complex vertex, gc(2)
 
25
 
 
26
 
 
27
c Take the inner product between the tensor particles
 
28
c and multiply it with the scalar particle and the coupling constant.
 
29
c Note that the tensor particle is antisymmetric, thus all diagonal terms
 
30
c are zero.
 
31
 
 
32
      if (gc(1).NE.(0D0,0D0)) then
 
33
 
 
34
      vertex = - gc(1)*sc(1)* (
 
35
     &                 !      + tc1( 1) * tc2( 1)
 
36
     &                       - tc1( 2) * tc2( 2)
 
37
     &                       - tc1( 3) * tc2( 3)
 
38
     &                       - tc1( 4) * tc2( 4)
 
39
 
 
40
     &                       - tc1( 5) * tc2( 5)
 
41
     &                 !      + tc1( 6) * tc2( 6)
 
42
     &                       + tc1( 7) * tc2( 7)
 
43
     &                       + tc1( 8) * tc2( 8)
 
44
 
 
45
     &                       - tc1( 9) * tc2( 9)
 
46
     &                       + tc1(10) * tc2(10)
 
47
     &                 !      + tc1(11) * tc2(11)
 
48
     &                       + tc1(12) * tc2(12)
 
49
 
 
50
     &                       - tc1(13) * tc2(13)
 
51
     &                       + tc1(14) * tc2(14)
 
52
     &                       + tc1(15) * tc2(15)
 
53
     &                 !      + tc1(16) * tc2(16)
 
54
     &                                           )
 
55
 
 
56
 
 
57
      else
 
58
      vertex = (0D0,0D0)
 
59
      endif      
 
60
 
 
61
 
 
62
      return
 
63
      end