← Back to branch summary

~maddevelopers/mg5amcnlo/3.0.4

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%matrix_4.f

  • Committer: olivier-mattelaer
  • Date: 2021-03-25 09:12:41 UTC
  • mfrom: (78.443.30 2.9.3)
  • Revision ID: olivier-mattelaer-20210325091241-b9d2xyz0wwgc1jvs
merge with 2.9.3 -> forcing DP for CF computationn

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%matrix_4.f
9
9
C      orders.inc (NLO_ORDERS)
10
10
C     
11
11
C     
12
 
C     Process: g d~ > t t~ d~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
13
 
C     Process: g s~ > t t~ s~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
14
 
C     Process: g u~ > t t~ u~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
15
 
C     Process: g c~ > t t~ c~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
 
12
C     Process: g d~ > t t~ d~ [ real = QCD QED ] QCD^2<=6 QED^2<=0
 
13
C     Process: g s~ > t t~ s~ [ real = QCD QED ] QCD^2<=6 QED^2<=0
 
14
C     Process: g u~ > t t~ u~ [ real = QCD QED ] QCD^2<=6 QED^2<=0
 
15
C     Process: g c~ > t t~ c~ [ real = QCD QED ] QCD^2<=6 QED^2<=0
16
16
C     
17
17
C     
18
18
C     CONSTANTS
111
111
C     and helicities
112
112
C     for the point in phase space P(0:3,NEXTERNAL)
113
113
C     
114
 
C     Process: g d~ > t t~ d~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
115
 
C     Process: g s~ > t t~ s~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
116
 
C     Process: g u~ > t t~ u~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
117
 
C     Process: g c~ > t t~ c~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
 
114
C     Process: g d~ > t t~ d~ [ real = QCD QED ] QCD^2<=6 QED^2<=0
 
115
C     Process: g s~ > t t~ s~ [ real = QCD QED ] QCD^2<=6 QED^2<=0
 
116
C     Process: g u~ > t t~ u~ [ real = QCD QED ] QCD^2<=6 QED^2<=0
 
117
C     Process: g c~ > t t~ c~ [ real = QCD QED ] QCD^2<=6 QED^2<=0
118
118
C     
119
119
      IMPLICIT NONE
120
120
C     
238
238
C     Returns amplitude squared summed/avg over colors
239
239
C     for the point with external lines W(0:6,NEXTERNAL)
240
240
C     
241
 
C     Process: g d~ > t t~ d~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
242
 
C     Process: g s~ > t t~ s~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
243
 
C     Process: g u~ > t t~ u~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
244
 
C     Process: g c~ > t t~ c~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
 
241
C     Process: g d~ > t t~ d~ [ real = QCD QED ] QCD^2<=6 QED^2<=0
 
242
C     Process: g s~ > t t~ s~ [ real = QCD QED ] QCD^2<=6 QED^2<=0
 
243
C     Process: g u~ > t t~ u~ [ real = QCD QED ] QCD^2<=6 QED^2<=0
 
244
C     Process: g c~ > t t~ c~ [ real = QCD QED ] QCD^2<=6 QED^2<=0
245
245
C     
246
246
      IMPLICIT NONE
247
247
C     
271
271
      INTEGER I,J,M,N
272
272
      INTEGER IC(NEXTERNAL)
273
273
      DATA IC /NEXTERNAL*1/
274
 
      REAL*8 DENOM(NCOLOR), CF(NCOLOR,NCOLOR)
 
274
      REAL*8  CF(NCOLOR,NCOLOR)
275
275
      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP(NCOLOR,NAMPSO), W(8
276
276
     $ ,NWAVEFUNCS)
 
277
      COMPLEX*16 TMP_JAMP(0)
277
278
C     
278
279
C     FUNCTION
279
280
C     
281
282
C     
282
283
C     COLOR DATA
283
284
C     
284
 
      DATA DENOM(1)/1/
285
 
      DATA (CF(I,  1),I=  1,  4) /   12,    4,    4,    0/
 
285
      DATA (CF(I,  1),I=  1,  4) /1.200000000000000D+01
 
286
     $ ,4.000000000000000D+00,4.000000000000000D+00,0.000000000000000D
 
287
     $ +00/
286
288
C     1 T(1,2,4) T(3,5)
287
 
      DATA DENOM(2)/1/
288
 
      DATA (CF(I,  2),I=  1,  4) /    4,   12,    0,    4/
 
289
      DATA (CF(I,  2),I=  1,  4) /4.000000000000000D+00
 
290
     $ ,1.200000000000000D+01,0.000000000000000D+00,4.000000000000000D
 
291
     $ +00/
289
292
C     1 T(1,2,5) T(3,4)
290
 
      DATA DENOM(3)/1/
291
 
      DATA (CF(I,  3),I=  1,  4) /    4,    0,   12,    4/
 
293
      DATA (CF(I,  3),I=  1,  4) /4.000000000000000D+00
 
294
     $ ,0.000000000000000D+00,1.200000000000000D+01,4.000000000000000D
 
295
     $ +00/
292
296
C     1 T(1,3,4) T(2,5)
293
 
      DATA DENOM(4)/1/
294
 
      DATA (CF(I,  4),I=  1,  4) /    0,    4,    4,   12/
 
297
      DATA (CF(I,  4),I=  1,  4) /0.000000000000000D+00
 
298
     $ ,4.000000000000000D+00,4.000000000000000D+00,1.200000000000000D
 
299
     $ +01/
295
300
C     1 T(1,3,5) T(2,4)
296
301
C     ----------
297
302
C     BEGIN CODE
318
323
C     Amplitude(s) for diagram number 5
319
324
      CALL VVV1_0(W(1,1),W(1,8),W(1,7),GC_10,AMP(5))
320
325
C     JAMPs contributing to orders QCD=3 QED=0
321
 
      JAMP(1,1)=+1D0/2D0*(+AMP(1)+AMP(3)+IMAG1*AMP(5))
322
 
      JAMP(2,1)=+1D0/2D0*(-1D0/3D0*AMP(1)-1D0/3D0*AMP(4))
323
 
      JAMP(3,1)=+1D0/2D0*(-1D0/3D0*AMP(2)-1D0/3D0*AMP(3))
324
 
      JAMP(4,1)=+1D0/2D0*(+AMP(2)+AMP(4)-IMAG1*AMP(5))
 
326
      JAMP(1,1) = (5.000000000000000D-01)*AMP(1)+(5.000000000000000D
 
327
     $ -01)*AMP(3)+((0.000000000000000D+00,5.000000000000000D-01))
 
328
     $ *AMP(5)
 
329
      JAMP(2,1) = (-1.666666666666667D-01)*AMP(1)+(-1.666666666666667D
 
330
     $ -01)*AMP(4)
 
331
      JAMP(3,1) = (-1.666666666666667D-01)*AMP(2)+(-1.666666666666667D
 
332
     $ -01)*AMP(3)
 
333
      JAMP(4,1) = (5.000000000000000D-01)*AMP(2)+(5.000000000000000D
 
334
     $ -01)*AMP(4)+((0.000000000000000D+00,-5.000000000000000D-01))
 
335
     $ *AMP(5)
325
336
 
326
337
      DO I=0,NSQAMPSO
327
338
        RES(I)=0D0
334
345
          ENDDO
335
346
          DO N = 1, NAMPSO
336
347
            RES(SQSOINDEX4(M,N)) = RES(SQSOINDEX4(M,N)) + ZTEMP
337
 
     $       *DCONJG(JAMP(I,N))/DENOM(I)
 
348
     $       *DCONJG(JAMP(I,N))
338
349
          ENDDO
339
350
        ENDDO
340
351
      ENDDO