← Back to branch summary

~maddevelopers/mg5amcnlo/3.0.2-alpha0

« back to all changes in this revision

Viewing changes to HELAS/hsstxx.F

Added Template and HELAS into bzr

Show diffs side-by-side

added added

removed removed

Lines of Context:
HELAS/hsstxx.F
 
1
      subroutine hsstxx(tc1,sc,gt,xm,xw,jts)
 
2
c
 
3
c- by RF - Feb. 2006 
 
4
c
 
5
c This subroutine computes an off-shell tensor current from the tts coupling.
 
6
c
 
7
c     input:
 
8
c          complex tc1               : Incoming tensor particle
 
9
c          complex sc                : Incoming scalar particle (Higgs)
 
10
c          real    gt                : coupling constant for the tts vertex
 
11
c
 
12
c     output:
 
13
c          complex jts               : off-shell tensor current
 
14
c
 
15
c     not used:
 
16
c          xm, xw
 
17
c
 
18
 
 
19
      implicit none
 
20
c--   dimension of the current set to arbitrary length
 
21
      INTEGER DIM
 
22
      PARAMETER(DIM=18)
 
23
c      include 'dim.inc'
 
24
      double complex tc1(DIM),jts(DIM),sc(DIM)
 
25
      double precision gt, xm, xw
 
26
 
 
27
c The outgoing tensor current is the same as the incoming multiplied by the
 
28
c coupling constant and the scalar particle.
 
29
c Note that the diagonal tensor terms are always zero because
 
30
c the tensor particle is anti-symmetric.
 
31
 
 
32
      jts( 1) = 0 !gt * sc(1) * tc1( 1)
 
33
      jts( 2) =  gt * sc(1) * tc1( 2)
 
34
      jts( 3) =  gt * sc(1) * tc1( 3)
 
35
      jts( 4) =  gt * sc(1) * tc1( 4)
 
36
 
 
37
      jts( 5) =  gt * sc(1) * tc1( 5)
 
38
      jts( 6) = 0 !gt * sc(1) * tc1( 6)
 
39
      jts( 7) =  gt * sc(1) * tc1( 7)
 
40
      jts( 8) =  gt * sc(1) * tc1( 8)
 
41
 
 
42
      jts( 9) =  gt * sc(1) * tc1( 9)
 
43
      jts(10) =  gt * sc(1) * tc1(10)
 
44
      jts(11) = 0 !gt * sc(1) * tc1(11)
 
45
      jts(12) =  gt * sc(1) * tc1(12)
 
46
 
 
47
      jts(13) =  gt * sc(1) * tc1(13)
 
48
      jts(14) =  gt * sc(1) * tc1(14)
 
49
      jts(15) =  gt * sc(1) * tc1(15)
 
50
      jts(16) = 0 !gt * sc(1) * tc1(16)
 
51
 
 
52
      jts(17) = sc(2) + tc1(17)
 
53
      jts(18) = sc(3) + tc1(18)
 
54
 
 
55
      return
 
56
      end