← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/MadLoop_output_from_the_interface/TIR_output/%ggttx_IOTest%SubProcesses%P0_gg_ttx%born_matrix.f

  • Committer: olivier Mattelaer
  • Date: 2015-03-05 00:14:16 UTC
  • mfrom: (258.1.9 2.3)
  • mto: (258.8.1 2.3)
  • mto: This revision was merged to the branch mainline in revision 259.
  • Revision ID: olivier.mattelaer@uclouvain.be-20150305001416-y9mzeykfzwnl9t0j
partial merge

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/MadLoop_output_from_the_interface/TIR_output/%ggttx_IOTest%SubProcesses%P0_gg_ttx%born_matrix.f
67
67
C     and helicities
68
68
C     for the point in phase space P(0:3,NEXTERNAL)
69
69
C     
70
 
C     Process: g g > t t~ [ QCD ]
 
70
C     Process: g g > t t~ [ virt = QCD ]
71
71
C     
72
72
      IMPLICIT NONE
73
73
C     
98
98
      LOGICAL GOODHEL(NCOMB)
99
99
      DATA NTRY/0/
100
100
      DATA GOODHEL/NCOMB*.FALSE./
101
 
      DATA (NHEL(I,   1),I=1,4) /-1,-1,-1,-1/
102
 
      DATA (NHEL(I,   2),I=1,4) /-1,-1,-1, 1/
103
 
      DATA (NHEL(I,   3),I=1,4) /-1,-1, 1,-1/
104
 
      DATA (NHEL(I,   4),I=1,4) /-1,-1, 1, 1/
105
 
      DATA (NHEL(I,   5),I=1,4) /-1, 1,-1,-1/
106
 
      DATA (NHEL(I,   6),I=1,4) /-1, 1,-1, 1/
107
 
      DATA (NHEL(I,   7),I=1,4) /-1, 1, 1,-1/
108
 
      DATA (NHEL(I,   8),I=1,4) /-1, 1, 1, 1/
109
 
      DATA (NHEL(I,   9),I=1,4) / 1,-1,-1,-1/
110
 
      DATA (NHEL(I,  10),I=1,4) / 1,-1,-1, 1/
111
 
      DATA (NHEL(I,  11),I=1,4) / 1,-1, 1,-1/
112
 
      DATA (NHEL(I,  12),I=1,4) / 1,-1, 1, 1/
113
 
      DATA (NHEL(I,  13),I=1,4) / 1, 1,-1,-1/
114
 
      DATA (NHEL(I,  14),I=1,4) / 1, 1,-1, 1/
115
 
      DATA (NHEL(I,  15),I=1,4) / 1, 1, 1,-1/
116
 
      DATA (NHEL(I,  16),I=1,4) / 1, 1, 1, 1/
 
101
      DATA (NHEL(I,   1),I=1,4) /-1,-1,-1, 1/
 
102
      DATA (NHEL(I,   2),I=1,4) /-1,-1,-1,-1/
 
103
      DATA (NHEL(I,   3),I=1,4) /-1,-1, 1, 1/
 
104
      DATA (NHEL(I,   4),I=1,4) /-1,-1, 1,-1/
 
105
      DATA (NHEL(I,   5),I=1,4) /-1, 1,-1, 1/
 
106
      DATA (NHEL(I,   6),I=1,4) /-1, 1,-1,-1/
 
107
      DATA (NHEL(I,   7),I=1,4) /-1, 1, 1, 1/
 
108
      DATA (NHEL(I,   8),I=1,4) /-1, 1, 1,-1/
 
109
      DATA (NHEL(I,   9),I=1,4) / 1,-1,-1, 1/
 
110
      DATA (NHEL(I,  10),I=1,4) / 1,-1,-1,-1/
 
111
      DATA (NHEL(I,  11),I=1,4) / 1,-1, 1, 1/
 
112
      DATA (NHEL(I,  12),I=1,4) / 1,-1, 1,-1/
 
113
      DATA (NHEL(I,  13),I=1,4) / 1, 1,-1, 1/
 
114
      DATA (NHEL(I,  14),I=1,4) / 1, 1,-1,-1/
 
115
      DATA (NHEL(I,  15),I=1,4) / 1, 1, 1, 1/
 
116
      DATA (NHEL(I,  16),I=1,4) / 1, 1, 1,-1/
117
117
      DATA IDEN/256/
118
118
C     
119
119
C     GLOBAL VARIABLES
201
201
C     Returns amplitude squared summed/avg over colors
202
202
C     for the point with external lines W(0:6,NEXTERNAL)
203
203
C     
204
 
C     Process: g g > t t~ [ QCD ]
 
204
C     Process: g g > t t~ [ virt = QCD ]
205
205
C     
206
206
      IMPLICIT NONE
207
207
C     
317
317
      INTEGER I, SQORDERS(NSO)
318
318
      INTEGER AMPSPLITORDERS(NAMPSO,NSO)
319
319
      DATA (AMPSPLITORDERS(  1,I),I=  1,  1) /    2/
 
320
      COMMON/ML5_0_AMPSPLITORDERS/AMPSPLITORDERS
320
321
C     
321
322
C     FUNCTION
322
323
C     
352
353
      INTEGER I,J
353
354
      INTEGER SQSPLITORDERS(NSQSO,NSO)
354
355
      DATA (SQSPLITORDERS(  1,I),I=  1,  1) /    4/
 
356
      COMMON/ML5_0_SQPLITORDERS/SQPLITORDERS
355
357
C     
356
358
C     BEGIN CODE
357
359
C     
364
366
 1009   CONTINUE
365
367
      ENDDO
366
368
 
367
 
      WRITE(*,*) 'ERROR:: Stopping function ML5_0_SOINDEX_FOR_SQUARED_'
368
 
     $ //'ORDERS'
 
369
      WRITE(*,*) 'ERROR:: Stopping in function'
 
370
      WRITE(*,*) 'ML5_0_SOINDEX_FOR_SQUARED_ORDERS'
369
371
      WRITE(*,*) 'Could not find squared orders ',(ORDERS(I),I=1,NSO)
370
372
      STOP
371
373
 
386
388
 
387
389
      END
388
390
 
 
391
C     This is the inverse subroutine of SOINDEX_FOR_SQUARED_ORDERS.
 
392
C      Not directly useful, but provided nonetheless.
 
393
      SUBROUTINE ML5_0_GET_SQUARED_ORDERS_FOR_SOINDEX(SOINDEX,ORDERS)
 
394
C     
 
395
C     This functions returns the orders identified by the squared
 
396
C      split order index in argument. Order values correspond to
 
397
C      following list of couplings (and in this order):
 
398
C     ['QCD']
 
399
C     
 
400
C     CONSTANTS
 
401
C     
 
402
      INTEGER    NSO, NSQSO
 
403
      PARAMETER (NSO=1, NSQSO=1)
 
404
C     
 
405
C     ARGUMENTS
 
406
C     
 
407
      INTEGER SOINDEX, ORDERS(NSO)
 
408
C     
 
409
C     LOCAL VARIABLES
 
410
C     
 
411
      INTEGER I
 
412
      INTEGER SQPLITORDERS(NSQSO,NSO)
 
413
      COMMON/ML5_0_SQPLITORDERS/SQPLITORDERS
 
414
C     
 
415
C     BEGIN CODE
 
416
C     
 
417
      IF (SOINDEX.GT.0.AND.SOINDEX.LE.NSQSO) THEN
 
418
        DO I=1,NSO
 
419
          ORDERS(I) =  SQPLITORDERS(SOINDEX,I)
 
420
        ENDDO
 
421
        RETURN
 
422
      ENDIF
 
423
 
 
424
      WRITE(*,*) 'ERROR:: Stopping function ML5_0_GET_SQUARED_ORDERS_F'
 
425
     $ //'OR_SOINDEX'
 
426
      WRITE(*,*) 'Could not find squared orders index ',SOINDEX
 
427
      STOP
 
428
 
 
429
      END SUBROUTINE
 
430
 
 
431
C     This is the inverse subroutine of getting amplitude SO orders.
 
432
C      Not directly useful, but provided nonetheless.
 
433
      SUBROUTINE ML5_0_GET_ORDERS_FOR_AMPSOINDEX(SOINDEX,ORDERS)
 
434
C     
 
435
C     This functions returns the orders identified by the split order
 
436
C      index in argument. Order values correspond to following list of
 
437
C      couplings (and in this order):
 
438
C     ['QCD']
 
439
C     
 
440
C     CONSTANTS
 
441
C     
 
442
      INTEGER    NSO, NAMPSO
 
443
      PARAMETER (NSO=1, NAMPSO=1)
 
444
C     
 
445
C     ARGUMENTS
 
446
C     
 
447
      INTEGER SOINDEX, ORDERS(NSO)
 
448
C     
 
449
C     LOCAL VARIABLES
 
450
C     
 
451
      INTEGER I
 
452
      INTEGER AMPSPLITORDERS(NAMPSO,NSO)
 
453
      COMMON/ML5_0_AMPSPLITORDERS/AMPSPLITORDERS
 
454
C     
 
455
C     BEGIN CODE
 
456
C     
 
457
      IF (SOINDEX.GT.0.AND.SOINDEX.LE.NAMPSO) THEN
 
458
        DO I=1,NSO
 
459
          ORDERS(I) =  AMPSPLITORDERS(SOINDEX,I)
 
460
        ENDDO
 
461
        RETURN
 
462
      ENDIF
 
463
 
 
464
      WRITE(*,*) 'ERROR:: Stopping function ML5_0_GET_ORDERS_FOR_AMPSO'
 
465
     $ //'INDEX'
 
466
      WRITE(*,*) 'Could not find amplitude split orders index ',SOINDEX
 
467
      STOP
 
468
 
 
469
      END SUBROUTINE
 
470
 
 
471
C     This function is not directly useful, but included for completene
 
472
C     ss
 
473
      INTEGER FUNCTION ML5_0_SOINDEX_FOR_AMPORDERS(ORDERS)
 
474
C     
 
475
C     This functions returns the integer index identifying the
 
476
C      amplitude split orders passed in argument which correspond to
 
477
C      the values of the following list of couplings (and in this
 
478
C      order):
 
479
C     ['QCD']
 
480
C     
 
481
C     CONSTANTS
 
482
C     
 
483
      INTEGER    NSO, NAMPSO
 
484
      PARAMETER (NSO=1, NAMPSO=1)
 
485
C     
 
486
C     ARGUMENTS
 
487
C     
 
488
      INTEGER ORDERS(NSO)
 
489
C     
 
490
C     LOCAL VARIABLES
 
491
C     
 
492
      INTEGER I,J
 
493
      INTEGER AMPSPLITORDERS(NAMPSO,NSO)
 
494
      COMMON/ML5_0_AMPSPLITORDERS/AMPSPLITORDERS
 
495
C     
 
496
C     BEGIN CODE
 
497
C     
 
498
      DO I=1,NAMPSO
 
499
        DO J=1,NSO
 
500
          IF (ORDERS(J).NE.AMPSPLITORDERS(I,J)) GOTO 1009
 
501
        ENDDO
 
502
        ML5_0_SOINDEX_FOR_AMPORDERS = I
 
503
        RETURN
 
504
 1009   CONTINUE
 
505
      ENDDO
 
506
 
 
507
      WRITE(*,*) 'ERROR:: Stopping function ML5_0_SOINDEX_FOR_AMPORDER'
 
508
     $ //'S'
 
509
      WRITE(*,*) 'Could not find squared orders ',(ORDERS(I),I=1,NSO)
 
510
      STOP
 
511
 
 
512
      END
 
513