← Back to branch summary

~maddevelopers/mg5amcnlo/2.7.1.3

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/FKS/gg_ttx/born_hel.f

  • Committer: olivier Mattelaer
  • Date: 2016-05-12 11:00:18 UTC
  • mfrom: (262.1.150 2.3.4)
  • Revision ID: olivier.mattelaer@uclouvain.be-20160512110018-sevb79f0wm4g8mpp
pass to 2.4.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/FKS/gg_ttx/born_hel.f
1
 
      SUBROUTINE SBORN_HEL(P1,ANS)
2
 
C     
3
 
C     Generated by MadGraph5_aMC@NLO v. 2.1.1.alpha, 2014-02-22
4
 
C     By the MadGraph5_aMC@NLO Development Team
5
 
C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
6
 
C     
7
 
C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
8
 
C     AND HELICITIES
9
 
C     FOR THE POINT IN PHASE SPACE P1(0:3,NEXTERNAL-1)
10
 
C     
11
 
C     Process: g g > t t~ WEIGHTED=2 [ QCD ]
12
 
C     
13
 
      IMPLICIT NONE
14
 
C     
15
 
C     CONSTANTS
16
 
C     
17
 
      INCLUDE 'nexternal.inc'
18
 
      INCLUDE 'born_nhel.inc'
19
 
      INTEGER     NCOMB
20
 
      PARAMETER ( NCOMB=  16 )
21
 
      INTEGER    THEL
22
 
      PARAMETER (THEL=NCOMB*8)
23
 
      INTEGER NGRAPHS
24
 
      PARAMETER (NGRAPHS = 3)
25
 
C     
26
 
C     ARGUMENTS 
27
 
C     
28
 
      REAL*8 P1(0:3,NEXTERNAL-1),ANS
29
 
C     
30
 
C     LOCAL VARIABLES 
31
 
C     
32
 
      INTEGER IHEL,IDEN,J
33
 
      REAL*8 BORN_HEL
34
 
      INTEGER IDEN_VALUES(8)
35
 
      DATA IDEN_VALUES /256, 256, 256, 256, 256, 256, 256, 256/
36
 
C     
37
 
C     GLOBAL VARIABLES
38
 
C     
39
 
      LOGICAL GOODHEL(NCOMB,8)
40
 
      COMMON /C_GOODHEL/ GOODHEL
41
 
      DOUBLE PRECISION SAVEMOM(NEXTERNAL-1,2)
42
 
      COMMON/TO_SAVEMOM/SAVEMOM
43
 
      LOGICAL CALCULATEDBORN
44
 
      COMMON/CCALCULATEDBORN/CALCULATEDBORN
45
 
      INTEGER NFKSPROCESS
46
 
      COMMON/C_NFKSPROCESS/NFKSPROCESS
47
 
      DOUBLE PRECISION WGT_HEL(MAX_BHEL)
48
 
      COMMON/C_BORN_HEL/WGT_HEL
49
 
 
50
 
C     ----------
51
 
C     BEGIN CODE
52
 
C     ----------
53
 
      IDEN=IDEN_VALUES(NFKSPROCESS)
54
 
      IF (CALCULATEDBORN) THEN
55
 
        DO J=1,NEXTERNAL-1
56
 
          IF (SAVEMOM(J,1).NE.P1(0,J) .OR. SAVEMOM(J,2).NE.P1(3
57
 
     $     ,J)) THEN
58
 
            CALCULATEDBORN=.FALSE.
59
 
            WRITE (*,*) 'momenta not the same in Born_hel'
60
 
            STOP
61
 
          ENDIF
62
 
        ENDDO
63
 
      ELSE
64
 
        WRITE(*,*) 'Error in born_hel: should be called only with
65
 
     $    calculatedborn = true'
66
 
        STOP
67
 
      ENDIF
68
 
      ANS = 0D0
69
 
      DO IHEL=1,NCOMB
70
 
        WGT_HEL(IHEL)=0D0
71
 
        IF (GOODHEL(IHEL,NFKSPROCESS)) THEN
72
 
          WGT_HEL(IHEL)=BORN_HEL(P1,IHEL)/DBLE(IDEN)
73
 
          ANS=ANS+WGT_HEL(IHEL)
74
 
        ENDIF
75
 
      ENDDO
76
 
      END
77
 
 
78
 
 
79
 
      REAL*8 FUNCTION BORN_HEL(P,HELL)
80
 
C     
81
 
C     Generated by MadGraph5_aMC@NLO v. 2.1.1.alpha, 2014-02-22
82
 
C     By the MadGraph5_aMC@NLO Development Team
83
 
C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
84
 
C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
85
 
C     FOR THE POINT WITH EXTERNAL LINES W(0:6,NEXTERNAL-1)
86
 
 
87
 
C     Process: g g > t t~ WEIGHTED=2 [ QCD ]
88
 
C     
89
 
      IMPLICIT NONE
90
 
C     
91
 
C     CONSTANTS
92
 
C     
93
 
      INTEGER     NGRAPHS
94
 
      PARAMETER ( NGRAPHS = 3 )
95
 
      INTEGER    NCOLOR
96
 
      PARAMETER (NCOLOR=2)
97
 
      REAL*8     ZERO
98
 
      PARAMETER (ZERO=0D0)
99
 
      COMPLEX*16 IMAG1
100
 
      PARAMETER (IMAG1 = (0D0,1D0))
101
 
      INCLUDE 'nexternal.inc'
102
 
      INCLUDE 'born_nhel.inc'
103
 
C     
104
 
C     ARGUMENTS 
105
 
C     
106
 
      REAL*8 P(0:3,NEXTERNAL-1)
107
 
      INTEGER HELL
108
 
C     
109
 
C     LOCAL VARIABLES 
110
 
C     
111
 
      INTEGER I,J
112
 
      REAL*8 DENOM(NCOLOR), CF(NCOLOR,NCOLOR)
113
 
      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP(NCOLOR)
114
 
C     
115
 
C     GLOBAL VARIABLES
116
 
C     
117
 
      DOUBLE COMPLEX SAVEAMP(NGRAPHS,MAX_BHEL)
118
 
      COMMON/TO_SAVEAMP/SAVEAMP
119
 
      LOGICAL CALCULATEDBORN
120
 
      COMMON/CCALCULATEDBORN/CALCULATEDBORN
121
 
C     
122
 
C     COLOR DATA
123
 
C     
124
 
      DATA DENOM(1)/3/
125
 
      DATA (CF(I,  1),I=  1,  2) /   16,   -2/
126
 
C     1 T(1,2,3,4)
127
 
      DATA DENOM(2)/3/
128
 
      DATA (CF(I,  2),I=  1,  2) /   -2,   16/
129
 
C     1 T(2,1,3,4)
130
 
C     ----------
131
 
C     BEGIN CODE
132
 
C     ----------
133
 
      IF (.NOT. CALCULATEDBORN) THEN
134
 
        WRITE(*,*) 'Error in born_hel.f: this should be called only
135
 
     $    with calculatedborn = true'
136
 
        STOP
137
 
      ELSEIF (CALCULATEDBORN) THEN
138
 
        DO I=1,NGRAPHS
139
 
          AMP(I)=SAVEAMP(I,HELL)
140
 
        ENDDO
141
 
      ENDIF
142
 
      JAMP(1)=+IMAG1*AMP(1)-AMP(2)
143
 
      JAMP(2)=-IMAG1*AMP(1)-AMP(3)
144
 
      BORN_HEL = 0.D0
145
 
      DO I = 1, NCOLOR
146
 
        ZTEMP = (0.D0,0.D0)
147
 
        DO J = 1, NCOLOR
148
 
          ZTEMP = ZTEMP + CF(J,I)*JAMP(J)
149
 
        ENDDO
150
 
        BORN_HEL =BORN_HEL+ZTEMP*DCONJG(JAMP(I))/DENOM(I)
151
 
      ENDDO
152
 
      END
153