← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/short_ML_SMQCD_LoopInduced/gg_hh/loop_matrix.f

merged with 2.3 rev286

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/short_ML_SMQCD_LoopInduced/gg_hh/loop_matrix.f
73
73
C     and helicities for the point in phase space P(0:3,NEXTERNAL)
74
74
C     and external lines W(0:6,NEXTERNAL)
75
75
C     
76
 
C     Process: g g > h h QED=2 QCD=2 [ virt = QCD ]
 
76
C     Process: g g > h h QED=2 QCD=2 [ sqrvirt = QCD ]
77
77
C     
78
78
      IMPLICIT NONE
79
79
C     
311
311
        ML_INIT = .FALSE.
312
312
      ENDIF
313
313
 
 
314
C     Setup the file paths
 
315
      CALL JOINPATH(MLPATH,PARAMFNAME,PARAMFN)
 
316
      CALL JOINPATH(MLPATH,PROC_PREFIX,TMP)
 
317
      CALL JOINPATH(TMP,HELCONFIGFNAME,HELCONFIGFN)
 
318
      CALL JOINPATH(TMP,LOOPFILTERFNAME,LOOPFILTERFN)
 
319
      CALL JOINPATH(TMP,COLORNUMFNAME,COLORNUMFN)
 
320
      CALL JOINPATH(TMP,COLORDENOMFNAME,COLORDENOMFN)
 
321
      CALL JOINPATH(TMP,HELFILTERFNAME,HELFILTERFN)
 
322
 
 
323
      OPEN(1, FILE=COLORNUMFN, ERR=104, STATUS='OLD',          
 
324
     $  ACTION='READ')
 
325
      DO I=1,NLOOPAMPS
 
326
        READ(1,*,END=105) (CF_N(I,J),J=1,NLOOPAMPS)
 
327
      ENDDO
 
328
      GOTO 105
 
329
 104  CONTINUE
 
330
      STOP 'Color factors could not be initialized from file ML5_0_Col'
 
331
     $ //'orNumFactors.dat. File not found'
 
332
 105  CONTINUE
 
333
      CLOSE(1)
 
334
      OPEN(1, FILE=COLORDENOMFN, ERR=106, STATUS='OLD',          
 
335
     $  ACTION='READ')
 
336
      DO I=1,NLOOPAMPS
 
337
        READ(1,*,END=107) (CF_D(I,J),J=1,NLOOPAMPS)
 
338
      ENDDO
 
339
      GOTO 107
 
340
 106  CONTINUE
 
341
      STOP 'Color factors could not be initialized from file ML5_0_Col'
 
342
     $ //'orDenomFactors.dat. File not found'
 
343
 107  CONTINUE
 
344
      CLOSE(1)
 
345
      OPEN(1, FILE=HELCONFIGFN, ERR=108, STATUS='OLD',                
 
346
     $   ACTION='READ')
 
347
      DO H=1,NCOMB
 
348
        READ(1,*,END=109) (HELC(I,H),I=1,NEXTERNAL)
 
349
      ENDDO
 
350
      GOTO 109
 
351
 108  CONTINUE
 
352
      STOP 'Color helictiy configurations could not be initialize'
 
353
     $ //'d from file ML5_0_HelConfigs.dat. File not found'
 
354
 109  CONTINUE
 
355
      CLOSE(1)
 
356
      IF(BOOTANDSTOP) THEN
 
357
        WRITE(*,*) 'Stopped by user request.'
 
358
        STOP
 
359
      ENDIF
 
360
 
314
361
      IF(NTRY.EQ.0) THEN
315
 
 
316
 
C       Setup the file paths
317
 
        CALL JOINPATH(MLPATH,PARAMFNAME,PARAMFN)
318
 
        CALL JOINPATH(MLPATH,PROC_PREFIX,TMP)
319
 
        CALL JOINPATH(TMP,HELCONFIGFNAME,HELCONFIGFN)
320
 
        CALL JOINPATH(TMP,LOOPFILTERFNAME,LOOPFILTERFN)
321
 
        CALL JOINPATH(TMP,COLORNUMFNAME,COLORNUMFN)
322
 
        CALL JOINPATH(TMP,COLORDENOMFNAME,COLORDENOMFN)
323
 
        CALL JOINPATH(TMP,HELFILTERFNAME,HELFILTERFN)
324
 
 
325
362
        CALL ML5_0_SET_N_EVALS(N_DP_EVAL,N_QP_EVAL)
326
 
        HELDOUBLECHECKED=.NOT.DOUBLECHECKHELICITYFILTER
 
363
        HELDOUBLECHECKED=(.NOT.DOUBLECHECKHELICITYFILTER).OR.(HELICITYF
 
364
     $   ILTERLEVEL.EQ.0)
327
365
        DO J=1,NCOMB
328
366
          DO I=1,NCTAMPS
329
367
            GOODAMP(I,J)=.TRUE.
344
382
        ENDDO
345
383
 101    CONTINUE
346
384
        CLOSE(1)
 
385
        IF (HELICITYFILTERLEVEL.EQ.0) THEN
 
386
          FOUNDHELFILTER=.TRUE.
 
387
          DO J=1,NCOMB
 
388
            GOODHEL(J)=.TRUE.
 
389
          ENDDO
 
390
          GOTO 122
 
391
        ENDIF
347
392
        OPEN(1, FILE=HELFILTERFN, ERR=102, STATUS='OLD',          
348
393
     $    ACTION='READ')
349
394
        READ(1,*,END=103) (GOODHEL(I),I=1,NCOMB)
355
400
        ENDDO
356
401
 103    CONTINUE
357
402
        CLOSE(1)
358
 
        OPEN(1, FILE=COLORNUMFN, ERR=104, STATUS='OLD',          
359
 
     $    ACTION='READ')
360
 
        DO I=1,NLOOPAMPS
361
 
          READ(1,*,END=105) (CF_N(I,J),J=1,NLOOPAMPS)
362
 
        ENDDO
363
 
        GOTO 105
364
 
 104    CONTINUE
365
 
        STOP 'Color factors could not be initialized from fil'
366
 
     $   //'e ML5_0_ColorNumFactors.dat. File not found'
367
 
 105    CONTINUE
368
 
        CLOSE(1)
369
 
        OPEN(1, FILE=COLORDENOMFN, ERR=106, STATUS='OLD',          
370
 
     $    ACTION='READ')
371
 
        DO I=1,NLOOPAMPS
372
 
          READ(1,*,END=107) (CF_D(I,J),J=1,NLOOPAMPS)
373
 
        ENDDO
374
 
        GOTO 107
375
 
 106    CONTINUE
376
 
        STOP 'Color factors could not be initialized from fil'
377
 
     $   //'e ML5_0_ColorDenomFactors.dat. File not found'
378
 
 107    CONTINUE
379
 
        CLOSE(1)
380
 
        OPEN(1, FILE=HELCONFIGFN, ERR=108, STATUS='OLD',              
381
 
     $       ACTION='READ')
382
 
        DO H=1,NCOMB
383
 
          READ(1,*,END=109) (HELC(I,H),I=1,NEXTERNAL)
384
 
        ENDDO
385
 
        GOTO 109
386
 
 108    CONTINUE
387
 
        STOP 'Color helictiy configurations could not be initialize'
388
 
     $   //'d from file ML5_0_HelConfigs.dat. File not found'
389
 
 109    CONTINUE
390
 
        CLOSE(1)
391
 
        IF(BOOTANDSTOP) THEN
392
 
          WRITE(*,*) 'Stopped by user request.'
393
 
          STOP
394
 
        ENDIF
 
403
 122    CONTINUE
395
404
      ENDIF
396
405
 
397
406
      MP_DONE=.FALSE.
599
608
     $     ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16)
600
609
     $     ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(MDL_MB)
601
610
     $     ,CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5
602
 
     $     ,GC_33,MP__GC_33,GC_33,MP__GC_33,4,1,5,AMPL(1,5),S(5))
 
611
     $     ,GC_33,MP__GC_33,GC_33,MP__GC_33,4,1,1,5,AMPL(1,5),S(5))
603
612
C         Loop amplitude for loop diagram with ID 2
604
613
          CALL ML5_0_LOOP_4_4(1,1,2,3,4,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
605
614
     $     ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16)
606
615
     $     ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(MDL_MB)
607
616
     $     ,CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5
608
 
     $     ,GC_33,MP__GC_33,GC_33,MP__GC_33,4,1,6,AMPL(1,6),S(6))
 
617
     $     ,GC_33,MP__GC_33,GC_33,MP__GC_33,4,1,1,6,AMPL(1,6),S(6))
609
618
C         Loop amplitude for loop diagram with ID 3
610
619
          CALL ML5_0_LOOP_3_3(2,1,2,5,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
611
620
     $     ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16)
612
621
     $     ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5
613
 
     $     ,GC_5,MP__GC_5,GC_33,MP__GC_33,3,1,7,AMPL(1,7),S(7))
 
622
     $     ,GC_5,MP__GC_5,GC_33,MP__GC_33,3,1,1,7,AMPL(1,7),S(7))
614
623
C         Loop amplitude for loop diagram with ID 4
615
624
          CALL ML5_0_LOOP_3_3(3,1,2,5,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
616
625
     $     ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16)
617
626
     $     ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5
618
 
     $     ,GC_5,MP__GC_5,GC_33,MP__GC_33,3,1,8,AMPL(1,8),S(8))
 
627
     $     ,GC_5,MP__GC_5,GC_33,MP__GC_33,3,1,1,8,AMPL(1,8),S(8))
619
628
C         Loop amplitude for loop diagram with ID 5
620
629
          CALL ML5_0_LOOP_4_4(4,1,2,4,3,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
621
630
     $     ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16)
622
631
     $     ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(MDL_MB)
623
632
     $     ,CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5
624
 
     $     ,GC_33,MP__GC_33,GC_33,MP__GC_33,4,1,9,AMPL(1,9),S(9))
 
633
     $     ,GC_33,MP__GC_33,GC_33,MP__GC_33,4,1,1,9,AMPL(1,9),S(9))
625
634
C         Loop amplitude for loop diagram with ID 6
626
635
          CALL ML5_0_LOOP_4_4(5,1,3,2,4,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
627
636
     $     ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16)
628
637
     $     ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(MDL_MB)
629
638
     $     ,CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_33,MP__GC_33
630
 
     $     ,GC_5,MP__GC_5,GC_33,MP__GC_33,4,1,10,AMPL(1,10),S(10))
 
639
     $     ,GC_5,MP__GC_5,GC_33,MP__GC_33,4,1,1,10,AMPL(1,10),S(10))
631
640
C         Loop amplitude for loop diagram with ID 7
632
641
          CALL ML5_0_LOOP_4_4(4,1,2,3,4,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
633
642
     $     ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16)
634
643
     $     ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(MDL_MB)
635
644
     $     ,CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5
636
 
     $     ,GC_33,MP__GC_33,GC_33,MP__GC_33,4,1,11,AMPL(1,11),S(11))
 
645
     $     ,GC_33,MP__GC_33,GC_33,MP__GC_33,4,1,1,11,AMPL(1,11),S(11))
637
646
C         Loop amplitude for loop diagram with ID 8
638
647
          CALL ML5_0_LOOP_4_4(6,1,3,2,4,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB
639
648
     $     ,KIND=16),DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16)
640
649
     $     ,DCMPLX(MDL_MB),CMPLX(MP__MDL_MB,KIND=16),DCMPLX(MDL_MB)
641
650
     $     ,CMPLX(MP__MDL_MB,KIND=16),GC_5,MP__GC_5,GC_33,MP__GC_33
642
 
     $     ,GC_5,MP__GC_5,GC_33,MP__GC_33,4,1,12,AMPL(1,12),S(12))
 
651
     $     ,GC_5,MP__GC_5,GC_33,MP__GC_33,4,1,1,12,AMPL(1,12),S(12))
643
652
C         Loop amplitude for loop diagram with ID 9
644
653
          CALL ML5_0_LOOP_4_4(1,1,2,4,3,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
645
654
     $     ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16)
646
655
     $     ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(MDL_MT)
647
656
     $     ,CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5
648
 
     $     ,GC_37,MP__GC_37,GC_37,MP__GC_37,4,1,13,AMPL(1,13),S(13))
 
657
     $     ,GC_37,MP__GC_37,GC_37,MP__GC_37,4,1,1,13,AMPL(1,13),S(13))
649
658
C         Loop amplitude for loop diagram with ID 10
650
659
          CALL ML5_0_LOOP_4_4(1,1,2,3,4,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
651
660
     $     ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16)
652
661
     $     ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(MDL_MT)
653
662
     $     ,CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5
654
 
     $     ,GC_37,MP__GC_37,GC_37,MP__GC_37,4,1,14,AMPL(1,14),S(14))
 
663
     $     ,GC_37,MP__GC_37,GC_37,MP__GC_37,4,1,1,14,AMPL(1,14),S(14))
655
664
C         Loop amplitude for loop diagram with ID 11
656
665
          CALL ML5_0_LOOP_3_3(2,1,2,5,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
657
666
     $     ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16)
658
667
     $     ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5
659
 
     $     ,GC_5,MP__GC_5,GC_37,MP__GC_37,3,1,15,AMPL(1,15),S(15))
 
668
     $     ,GC_5,MP__GC_5,GC_37,MP__GC_37,3,1,1,15,AMPL(1,15),S(15))
660
669
C         Loop amplitude for loop diagram with ID 12
661
670
          CALL ML5_0_LOOP_3_3(3,1,2,5,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
662
671
     $     ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16)
663
672
     $     ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5
664
 
     $     ,GC_5,MP__GC_5,GC_37,MP__GC_37,3,1,16,AMPL(1,16),S(16))
 
673
     $     ,GC_5,MP__GC_5,GC_37,MP__GC_37,3,1,1,16,AMPL(1,16),S(16))
665
674
C         Loop amplitude for loop diagram with ID 13
666
675
          CALL ML5_0_LOOP_4_4(4,1,2,4,3,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
667
676
     $     ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16)
668
677
     $     ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(MDL_MT)
669
678
     $     ,CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5
670
 
     $     ,GC_37,MP__GC_37,GC_37,MP__GC_37,4,1,17,AMPL(1,17),S(17))
 
679
     $     ,GC_37,MP__GC_37,GC_37,MP__GC_37,4,1,1,17,AMPL(1,17),S(17))
671
680
C         Loop amplitude for loop diagram with ID 14
672
681
          CALL ML5_0_LOOP_4_4(5,1,3,2,4,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
673
682
     $     ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16)
674
683
     $     ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(MDL_MT)
675
684
     $     ,CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_37,MP__GC_37
676
 
     $     ,GC_5,MP__GC_5,GC_37,MP__GC_37,4,1,18,AMPL(1,18),S(18))
 
685
     $     ,GC_5,MP__GC_5,GC_37,MP__GC_37,4,1,1,18,AMPL(1,18),S(18))
677
686
C         Loop amplitude for loop diagram with ID 15
678
687
          CALL ML5_0_LOOP_4_4(4,1,2,3,4,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
679
688
     $     ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16)
680
689
     $     ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(MDL_MT)
681
690
     $     ,CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_5,MP__GC_5
682
 
     $     ,GC_37,MP__GC_37,GC_37,MP__GC_37,4,1,19,AMPL(1,19),S(19))
 
691
     $     ,GC_37,MP__GC_37,GC_37,MP__GC_37,4,1,1,19,AMPL(1,19),S(19))
683
692
C         Loop amplitude for loop diagram with ID 16
684
693
          CALL ML5_0_LOOP_4_4(6,1,3,2,4,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT
685
694
     $     ,KIND=16),DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16)
686
695
     $     ,DCMPLX(MDL_MT),CMPLX(MP__MDL_MT,KIND=16),DCMPLX(MDL_MT)
687
696
     $     ,CMPLX(MP__MDL_MT,KIND=16),GC_5,MP__GC_5,GC_37,MP__GC_37
688
 
     $     ,GC_5,MP__GC_5,GC_37,MP__GC_37,4,1,20,AMPL(1,20),S(20))
 
697
     $     ,GC_5,MP__GC_5,GC_37,MP__GC_37,4,1,1,20,AMPL(1,20),S(20))
689
698
          HELPICKED=HELPICKED_BU
690
699
          DO I=NCTAMPS+1,NLOOPAMPS
691
700
            IF((CTMODERUN.NE.-1).AND..NOT.CHECKPHASE.AND.(.NOT.S(I))