← Back to branch summary

~maddevelopers/mg5amcnlo/new_clustering

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/short_ML_SMQCD_default/ddx_ttx/born_matrix.f

  • Committer: Rikkert Frederix
  • Date: 2021-09-09 15:51:40 UTC
  • mfrom: (78.75.502 3.2.1)
  • Revision ID: frederix@physik.uzh.ch-20210909155140-rg6umfq68h6h47cf
merge with 3.2.1

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/short_ML_SMQCD_default/ddx_ttx/born_matrix.f
87
87
      INTEGER USERHEL
88
88
      COMMON/ML5_0_HELUSERCHOICE/USERHEL
89
89
      DATA USERHEL/-1/
 
90
      LOGICAL HELRESET
 
91
      COMMON/ML5_0_HELRESET/HELRESET
 
92
      DATA HELRESET/.TRUE./
90
93
 
91
94
      DATA (NHEL(I,   1),I=1,4) / 1,-1,-1, 1/
92
95
      DATA (NHEL(I,   2),I=1,4) / 1,-1,-1,-1/
116
119
      LOGICAL ML5_0_IS_BORN_HEL_SELECTED
117
120
 
118
121
C     ----------
 
122
C     Check if helreset mode is on
 
123
C     ---------
 
124
      IF (HELRESET) THEN
 
125
        NTRY = 0
 
126
        DO I=1,NCOMB
 
127
          GOODHEL(I) = .FALSE.
 
128
        ENDDO
 
129
        HELRESET = .FALSE.
 
130
      ENDIF
 
131
 
 
132
C     ----------
119
133
C     BEGIN CODE
120
134
C     ----------
121
135
      IF(USERHEL.EQ.-1) NTRY=NTRY+1
148
162
              CYCLE
149
163
            ENDIF
150
164
            T=ML5_0_MATRIX(P ,NHEL(1,IHEL),JC(1))
151
 
            IF(POLARIZATIONS(0,0).EQ.-1.OR.ML5_0_IS_BORN_HEL_SELECTED(I
152
 
     $HEL)) THEN
 
165
            IF(POLARIZATIONS(0,0).EQ.
 
166
     $       -1.OR.ML5_0_IS_BORN_HEL_SELECTED(IHEL)) THEN
153
167
              ANS=ANS+T
154
168
            ENDIF
155
169
            IF (T .NE. 0D0 .AND. .NOT.    GOODHEL(IHEL)) THEN
178
192
C     By the MadGraph5_aMC@NLO Development Team
179
193
C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
180
194
C     
181
 
C     Returns amplitude squared summed/avg over colors
 
195
C     Returns amplitude squared -- no average over initial
 
196
C      state/symmetry factor
182
197
C     for the point with external lines W(0:6,NEXTERNAL)
183
198
C     
184
199
C     Process: d d~ > t t~ QCD<=2 QED=0 [ virt = QCD ]
207
222
C     
208
223
      INTEGER I,J
209
224
      COMPLEX*16 ZTEMP
210
 
      REAL*8 DENOM(NCOLOR), CF(NCOLOR,NCOLOR)
211
 
      COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR)
 
225
      REAL*8 CF(NCOLOR,NCOLOR)
 
226
      COMPLEX*16 AMP(NGRAPHS), JAMP(NCOLOR), TMP_JAMP(0)
212
227
      COMPLEX*16 W(20,NWAVEFUNCS)
213
228
      COMPLEX*16 DUM0,DUM1
214
229
      DATA DUM0, DUM1/(0D0, 0D0), (1D0, 0D0)/
220
235
C     
221
236
C     COLOR DATA
222
237
C     
223
 
      DATA DENOM(1)/1/
224
 
      DATA (CF(I,  1),I=  1,  2) /    9,    3/
 
238
      DATA (CF(I,  1),I=  1,  2) /9.000000000000000D+00
 
239
     $ ,3.000000000000000D+00/
225
240
C     1 T(2,1) T(3,4)
226
 
      DATA DENOM(2)/1/
227
 
      DATA (CF(I,  2),I=  1,  2) /    3,    9/
 
241
      DATA (CF(I,  2),I=  1,  2) /3.000000000000000D+00
 
242
     $ ,9.000000000000000D+00/
228
243
C     1 T(2,4) T(3,1)
229
244
C     ----------
230
245
C     BEGIN CODE
245
260
        DO J = 1, NCOLOR
246
261
          ZTEMP = ZTEMP + CF(J,I)*JAMP(J)
247
262
        ENDDO
248
 
        ML5_0_MATRIX = ML5_0_MATRIX+ZTEMP*DCONJG(JAMP(I))/DENOM(I)
 
263
        ML5_0_MATRIX = ML5_0_MATRIX+ZTEMP*DCONJG(JAMP(I))
249
264
      ENDDO
250
265
 
251
266
      END