~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_udx_wp%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:
 
1
      SUBROUTINE SBORN_HEL(P1,ANS)
 
2
C     
 
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
 
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: u d~ > w+ WEIGHTED<=2 [ all = QCD ]
 
12
C     Process: c s~ > w+ WEIGHTED<=2 [ all = QCD ]
 
13
C     
 
14
      IMPLICIT NONE
 
15
C     
 
16
C     CONSTANTS
 
17
C     
 
18
      INCLUDE 'nexternal.inc'
 
19
      INCLUDE 'born_nhel.inc'
 
20
      INTEGER     NCOMB
 
21
      PARAMETER ( NCOMB=  12 )
 
22
      INTEGER    THEL
 
23
      PARAMETER (THEL=NCOMB*4)
 
24
      INTEGER NGRAPHS
 
25
      PARAMETER (NGRAPHS = 1)
 
26
C     
 
27
C     ARGUMENTS 
 
28
C     
 
29
      REAL*8 P1(0:3,NEXTERNAL-1),ANS
 
30
C     
 
31
C     LOCAL VARIABLES 
 
32
C     
 
33
      INTEGER IHEL,IDEN,J
 
34
      REAL*8 BORN_HEL
 
35
      INTEGER IDEN_VALUES(4)
 
36
      DATA IDEN_VALUES /36, 36, 36, 36/
 
37
C     
 
38
C     GLOBAL VARIABLES
 
39
C     
 
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
 
46
      INTEGER NFKSPROCESS
 
47
      COMMON/C_NFKSPROCESS/NFKSPROCESS
 
48
      DOUBLE PRECISION WGT_HEL(MAX_BHEL)
 
49
      COMMON/C_BORN_HEL/WGT_HEL
 
50
 
 
51
C     ----------
 
52
C     BEGIN CODE
 
53
C     ----------
 
54
      IDEN=IDEN_VALUES(NFKSPROCESS)
 
55
      IF (CALCULATEDBORN) THEN
 
56
        DO J=1,NEXTERNAL-1
 
57
          IF (SAVEMOM(J,1).NE.P1(0,J) .OR. SAVEMOM(J,2).NE.P1(3,J))
 
58
     $      THEN
 
59
            CALCULATEDBORN=.FALSE.
 
60
            WRITE (*,*) 'momenta not the same in Born_hel'
 
61
            STOP
 
62
          ENDIF
 
63
        ENDDO
 
64
      ELSE
 
65
        WRITE(*,*) 'Error in born_hel: should be called only with'
 
66
     $   //' calculatedborn = true'
 
67
        STOP
 
68
      ENDIF
 
69
      ANS = 0D0
 
70
      DO IHEL=1,NCOMB
 
71
        WGT_HEL(IHEL)=0D0
 
72
        IF (GOODHEL(IHEL,NFKSPROCESS)) THEN
 
73
          WGT_HEL(IHEL)=BORN_HEL(P1,IHEL)/DBLE(IDEN)
 
74
          ANS=ANS+WGT_HEL(IHEL)
 
75
        ENDIF
 
76
      ENDDO
 
77
      END
 
78
 
 
79
 
 
80
      REAL*8 FUNCTION BORN_HEL(P,HELL)
 
81
C     
 
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)
 
87
 
 
88
C     Process: u d~ > w+ WEIGHTED<=2 [ all = QCD ]
 
89
C     Process: c s~ > w+ WEIGHTED<=2 [ all = QCD ]
 
90
C     
 
91
      IMPLICIT NONE
 
92
C     
 
93
C     CONSTANTS
 
94
C     
 
95
      INTEGER     NGRAPHS
 
96
      PARAMETER ( NGRAPHS = 1 )
 
97
      INTEGER    NCOLOR
 
98
      PARAMETER (NCOLOR=1)
 
99
      REAL*8     ZERO
 
100
      PARAMETER (ZERO=0D0)
 
101
      COMPLEX*16 IMAG1
 
102
      PARAMETER (IMAG1 = (0D0,1D0))
 
103
      INCLUDE 'nexternal.inc'
 
104
      INCLUDE 'born_nhel.inc'
 
105
C     
 
106
C     ARGUMENTS 
 
107
C     
 
108
      REAL*8 P(0:3,NEXTERNAL-1)
 
109
      INTEGER HELL
 
110
C     
 
111
C     LOCAL VARIABLES 
 
112
C     
 
113
      INTEGER I,J
 
114
      REAL*8 DENOM(NCOLOR), CF(NCOLOR,NCOLOR)
 
115
      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP(NCOLOR)
 
116
C     
 
117
C     GLOBAL VARIABLES
 
118
C     
 
119
      DOUBLE COMPLEX SAVEAMP(NGRAPHS,MAX_BHEL)
 
120
      COMMON/TO_SAVEAMP/SAVEAMP
 
121
      LOGICAL CALCULATEDBORN
 
122
      COMMON/CCALCULATEDBORN/CALCULATEDBORN
 
123
C     
 
124
C     COLOR DATA
 
125
C     
 
126
      DATA DENOM(1)/1/
 
127
      DATA (CF(I,  1),I=  1,  1) /    3/
 
128
C     1 T(2,1)
 
129
C     ----------
 
130
C     BEGIN CODE
 
131
C     ----------
 
132
      IF (.NOT. CALCULATEDBORN) THEN
 
133
        WRITE(*,*) 'Error in born_hel.f: this should be called only'
 
134
     $   //' with calculatedborn = true'
 
135
        STOP
 
136
      ELSEIF (CALCULATEDBORN) THEN
 
137
        DO I=1,NGRAPHS
 
138
          AMP(I)=SAVEAMP(I,HELL)
 
139
        ENDDO
 
140
      ENDIF
 
141
      JAMP(1)=+AMP(1)
 
142
      BORN_HEL = 0.D0
 
143
      DO I = 1, NCOLOR
 
144
        ZTEMP = (0.D0,0.D0)
 
145
        DO J = 1, NCOLOR
 
146
          ZTEMP = ZTEMP + CF(J,I)*JAMP(J)
 
147
        ENDDO
 
148
        BORN_HEL =BORN_HEL+ZTEMP*DCONJG(JAMP(I))/DENOM(I)
 
149
      ENDDO
 
150
      END
 
151