← Back to branch summary

~maddevelopers/mg5amcnlo/2.5.3

« back to all changes in this revision

Viewing changes to HELAS/utsaxx.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/utsaxx.F
1
 
      subroutine utsaxx(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(18)           : Incoming tensor particle
9
 
c          complex sc(3)             : Incoming scalar particle
10
 
c          complex gt(2)             : coupling constant: gt(1) scalar
11
 
c                                                         gt(2) not used
12
 
c
13
 
c     output:
14
 
c          complex jts               : off-shell tensor current
15
 
c
16
 
c     not used:
17
 
c          xm, xw
18
 
c
19
 
 
20
 
      implicit none
21
 
c--   dimension of the current set to arbitrary length
22
 
      INTEGER DIM
23
 
      PARAMETER(DIM=18)
24
 
c      include "dimension.inc"
25
 
      double complex tc1(DIM),jts(DIM),sc(DIM), gt(2)
26
 
      double precision xm, xw
27
 
 
28
 
c The outgoing tensor current is the same as the incoming multiplied by the
29
 
c coupling constant and the scalar particle. Note that the diagonal tensor
30
 
c terms are always zero because the tensor particle is anti-symmetric. The 
31
 
c tensor particle does not propagate, thus no multiplication with the tensor 
32
 
c propagator.
33
 
 
34
 
      jts(17) = sc(2) + tc1(17)
35
 
      jts(18) = sc(3) + tc1(18)
36
 
 
37
 
 
38
 
      if (gt(1).NE.(0D0,0D0)) then
39
 
 
40
 
         jts( 1) = (0D0,0D0)        !-gt(1) * sc(1) * tc1( 1)
41
 
         jts( 2) =  -gt(1) * sc(1) * tc1( 2)
42
 
         jts( 3) =  -gt(1) * sc(1) * tc1( 3)
43
 
         jts( 4) =  -gt(1) * sc(1) * tc1( 4)
44
 
 
45
 
         jts( 5) =  -gt(1) * sc(1) * tc1( 5)
46
 
         jts( 6) = (0D0,0D0)        !-gt(1) * sc(1) * tc1( 6)
47
 
         jts( 7) =  -gt(1) * sc(1) * tc1( 7)
48
 
         jts( 8) =  -gt(1) * sc(1) * tc1( 8)
49
 
 
50
 
         jts( 9) =  -gt(1) * sc(1) * tc1( 9)
51
 
         jts(10) =  -gt(1) * sc(1) * tc1(10)
52
 
         jts(11) = (0D0,0D0)        !-gt(1) * sc(1) * tc1(11)
53
 
         jts(12) =  -gt(1) * sc(1) * tc1(12)
54
 
         
55
 
         jts(13) =  -gt(1) * sc(1) * tc1(13)
56
 
         jts(14) =  -gt(1) * sc(1) * tc1(14)
57
 
         jts(15) =  -gt(1) * sc(1) * tc1(15)
58
 
         jts(16) = (0D0,0D0)        !-gt(1) * sc(1) * tc1(16)
59
 
 
60
 
      else
61
 
         jts( 1)=(0D0,0D0)
62
 
         jts( 2)=(0D0,0D0)
63
 
         jts( 3)=(0D0,0D0)
64
 
         jts( 4)=(0D0,0D0)
65
 
         jts( 5)=(0D0,0D0)
66
 
         jts( 6)=(0D0,0D0)
67
 
         jts( 7)=(0D0,0D0)
68
 
         jts( 8)=(0D0,0D0)
69
 
         jts( 9)=(0D0,0D0)
70
 
         jts(10)=(0D0,0D0)
71
 
         jts(11)=(0D0,0D0)
72
 
         jts(12)=(0D0,0D0)
73
 
         jts(13)=(0D0,0D0)
74
 
         jts(14)=(0D0,0D0)
75
 
         jts(15)=(0D0,0D0)
76
 
         jts(16)=(0D0,0D0)
77
 
      endif
78
 
 
79
 
 
80
 
      return
81
 
      end