~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uux_ttx%born.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:
8
8
C     AND HELICITIES
9
9
C     FOR THE POINT IN PHASE SPACE P1(0:3,NEXTERNAL-1)
10
10
C     
11
 
C     Process: u u~ > t t~ WEIGHTED=2 [ LOonly = QCD ]
12
 
C     Process: c c~ > t t~ WEIGHTED=2 [ LOonly = QCD ]
13
 
C     Process: d d~ > t t~ WEIGHTED=2 [ LOonly = QCD ]
14
 
C     Process: s s~ > t t~ WEIGHTED=2 [ LOonly = QCD ]
 
11
C     Process: u u~ > t t~ WEIGHTED<=2 [ LOonly = QCD ]
 
12
C     Process: c c~ > t t~ WEIGHTED<=2 [ LOonly = QCD ]
 
13
C     Process: d d~ > t t~ WEIGHTED<=2 [ LOonly = QCD ]
 
14
C     Process: s s~ > t t~ WEIGHTED<=2 [ LOonly = QCD ]
15
15
C     
16
16
      IMPLICIT NONE
17
17
C     
79
79
      COMMON/CCALCULATEDBORN/CALCULATEDBORN
80
80
      INTEGER NFKSPROCESS
81
81
      COMMON/C_NFKSPROCESS/NFKSPROCESS
 
82
      DOUBLE PRECISION       WGT_ME_BORN,WGT_ME_REAL
 
83
      COMMON /C_WGT_ME_TREE/ WGT_ME_BORN,WGT_ME_REAL
82
84
C     ----------
83
85
C     BEGIN CODE
84
86
C     ----------
100
102
      ENDDO
101
103
      IF (CALCULATEDBORN) THEN
102
104
        DO J=1,NEXTERNAL-1
103
 
          IF (SAVEMOM(J,1).NE.P1(0,J) .OR. SAVEMOM(J,2).NE.P1(3
104
 
     $     ,J)) THEN
 
105
          IF (SAVEMOM(J,1).NE.P1(0,J) .OR. SAVEMOM(J,2).NE.P1(3,J))
 
106
     $      THEN
105
107
            CALCULATEDBORN=.FALSE.
106
108
            WRITE (*,*) 'momenta not the same in Born'
107
109
            STOP
124
126
      HEL_FAC=1D0
125
127
      DO IHEL=1,NCOMB
126
128
        IF (NHEL(GLU_IJ,IHEL).EQ.NHEL(GLU_IJ,1)) THEN
127
 
          IF ((GOODHEL(IHEL,NFKSPROCESS) .OR. GOODHEL(IHEL+SKIP(NFKSPRO
128
 
     $     CESS),NFKSPROCESS) .OR. NTRY(NFKSPROCESS) .LT. 2) ) THEN
 
129
          IF ((GOODHEL(IHEL,NFKSPROCESS) .OR. GOODHEL(IHEL
 
130
     $     +SKIP(NFKSPROCESS),NFKSPROCESS) .OR. NTRY(NFKSPROCESS) .LT.
 
131
     $      2) ) THEN
129
132
            ANS(1)=ANS(1)+BORN(P1,NHEL(1,IHEL),IHEL,BORNTILDE,BORNS)
130
133
            ANS(2)=ANS(2)+BORNTILDE
131
 
            IF ( BORNS(1).NE.0D0 .AND. .NOT. GOODHEL(IHEL,NFKSPROCESS
132
 
     $       ) ) THEN
 
134
            IF ( BORNS(1).NE.0D0 .AND. .NOT. GOODHEL(IHEL,NFKSPROCESS)
 
135
     $        ) THEN
133
136
              GOODHEL(IHEL,NFKSPROCESS)=.TRUE.
134
137
            ENDIF
135
 
            IF ( BORNS(2).NE.0D0 .AND. .NOT. GOODHEL(IHEL+SKIP(NFKSPROC
136
 
     $       ESS),NFKSPROCESS) ) THEN
 
138
            IF ( BORNS(2).NE.0D0 .AND. .NOT. GOODHEL(IHEL
 
139
     $       +SKIP(NFKSPROCESS),NFKSPROCESS) ) THEN
137
140
              GOODHEL(IHEL+SKIP(NFKSPROCESS),NFKSPROCESS)=.TRUE.
138
141
            ENDIF
139
142
          ENDIF
141
144
      ENDDO
142
145
      ANS(1)=ANS(1)/DBLE(IDEN)
143
146
      ANS(2)=ANS(2)/DBLE(IDEN)
 
147
      WGT_ME_BORN=DBLE(ANS(1))
144
148
      CALCULATEDBORN=.TRUE.
145
149
      END
146
150
 
153
157
C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
154
158
C     FOR THE POINT WITH EXTERNAL LINES W(0:6,NEXTERNAL-1)
155
159
 
156
 
C     Process: u u~ > t t~ WEIGHTED=2 [ LOonly = QCD ]
157
 
C     Process: c c~ > t t~ WEIGHTED=2 [ LOonly = QCD ]
158
 
C     Process: d d~ > t t~ WEIGHTED=2 [ LOonly = QCD ]
159
 
C     Process: s s~ > t t~ WEIGHTED=2 [ LOonly = QCD ]
 
160
C     Process: u u~ > t t~ WEIGHTED<=2 [ LOonly = QCD ]
 
161
C     Process: c c~ > t t~ WEIGHTED<=2 [ LOonly = QCD ]
 
162
C     Process: d d~ > t t~ WEIGHTED<=2 [ LOonly = QCD ]
 
163
C     Process: s s~ > t t~ WEIGHTED<=2 [ LOonly = QCD ]
160
164
C     
161
165
      IMPLICIT NONE
162
166
C     
188
192
      PARAMETER (NMO=NEXTERNAL-1)
189
193
      DATA IC /NMO*1/
190
194
      REAL*8 DENOM(NCOLOR), CF(NCOLOR,NCOLOR)
191
 
      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP(NCOLOR), W(8,NWAVEFUNCS)
192
 
     $ , JAMPH(2, NCOLOR)
 
195
      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP(NCOLOR), W(8,NWAVEFUNCS),
 
196
     $  JAMPH(2, NCOLOR)
193
197
C     
194
198
C     GLOBAL VARIABLES
195
199
C