1
SUBROUTINE SBORN_HEL(P1,ANS)
3
C Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
4
C By the MadGraph5_aMC@NLO Development Team
5
C Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
7
C RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
9
C FOR THE POINT IN PHASE SPACE P1(0:3,NEXTERNAL-1)
11
C Process: u d~ > w+ WEIGHTED<=2 [ all = QCD ]
12
C Process: c s~ > w+ WEIGHTED<=2 [ all = QCD ]
18
INCLUDE 'nexternal.inc'
19
INCLUDE 'born_nhel.inc'
21
PARAMETER ( NCOMB= 12 )
23
PARAMETER (THEL=NCOMB*4)
25
PARAMETER (NGRAPHS = 1)
29
REAL*8 P1(0:3,NEXTERNAL-1),ANS
35
INTEGER IDEN_VALUES(4)
36
DATA IDEN_VALUES /36, 36, 36, 36/
40
LOGICAL GOODHEL(NCOMB,4)
41
COMMON /C_GOODHEL/ GOODHEL
42
DOUBLE PRECISION SAVEMOM(NEXTERNAL-1,2)
43
COMMON/TO_SAVEMOM/SAVEMOM
44
LOGICAL CALCULATEDBORN
45
COMMON/CCALCULATEDBORN/CALCULATEDBORN
47
COMMON/C_NFKSPROCESS/NFKSPROCESS
48
DOUBLE PRECISION WGT_HEL(MAX_BHEL)
49
COMMON/C_BORN_HEL/WGT_HEL
54
IDEN=IDEN_VALUES(NFKSPROCESS)
55
IF (CALCULATEDBORN) THEN
57
IF (SAVEMOM(J,1).NE.P1(0,J) .OR. SAVEMOM(J,2).NE.P1(3,J))
59
CALCULATEDBORN=.FALSE.
60
WRITE (*,*) 'momenta not the same in Born_hel'
65
WRITE(*,*) 'Error in born_hel: should be called only with'
66
$ //' calculatedborn = true'
72
IF (GOODHEL(IHEL,NFKSPROCESS)) THEN
73
WGT_HEL(IHEL)=BORN_HEL(P1,IHEL)/DBLE(IDEN)
80
REAL*8 FUNCTION BORN_HEL(P,HELL)
82
C Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
83
C By the MadGraph5_aMC@NLO Development Team
84
C Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
85
C RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
86
C FOR THE POINT WITH EXTERNAL LINES W(0:6,NEXTERNAL-1)
88
C Process: u d~ > w+ WEIGHTED<=2 [ all = QCD ]
89
C Process: c s~ > w+ WEIGHTED<=2 [ all = QCD ]
96
PARAMETER ( NGRAPHS = 1 )
102
PARAMETER (IMAG1 = (0D0,1D0))
103
INCLUDE 'nexternal.inc'
104
INCLUDE 'born_nhel.inc'
108
REAL*8 P(0:3,NEXTERNAL-1)
114
REAL*8 DENOM(NCOLOR), CF(NCOLOR,NCOLOR)
115
COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP(NCOLOR)
119
DOUBLE COMPLEX SAVEAMP(NGRAPHS,MAX_BHEL)
120
COMMON/TO_SAVEAMP/SAVEAMP
121
LOGICAL CALCULATEDBORN
122
COMMON/CCALCULATEDBORN/CALCULATEDBORN
127
DATA (CF(I, 1),I= 1, 1) / 3/
132
IF (.NOT. CALCULATEDBORN) THEN
133
WRITE(*,*) 'Error in born_hel.f: this should be called only'
134
$ //' with calculatedborn = true'
136
ELSEIF (CALCULATEDBORN) THEN
138
AMP(I)=SAVEAMP(I,HELL)
146
ZTEMP = ZTEMP + CF(J,I)*JAMP(J)
148
BORN_HEL =BORN_HEL+ZTEMP*DCONJG(JAMP(I))/DENOM(I)