← 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_default/gg_ttx/loop_matrix.f

mwrge with 2.7.3

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/short_ML_SMQCD_default/gg_ttx/loop_matrix.f
700
700
 
701
701
      VALIDH=-1
702
702
      DO H=1,NCOMB
703
 
        IF ((HELPICKED.EQ.H).OR.((HELPICKED.EQ.-1).AND.(CHECKPHASE.OR.(
704
 
     $.NOT.HELDOUBLECHECKED).OR.GOODHEL(H)))) THEN
 
703
        IF ((HELPICKED.EQ.H).OR.((HELPICKED.EQ.-1)
 
704
     $   .AND.(CHECKPHASE.OR.(.NOT.HELDOUBLECHECKED).OR.GOODHEL(H))))
 
705
     $    THEN
705
706
 
706
707
C         Handle the possible requirement of specific polarizations
707
 
          IF ((.NOT.CHECKPHASE).AND.HELDOUBLECHECKED.AND.POLARIZATIONS(
708
 
     $0,0).EQ.0.AND.(.NOT.ML5_0_IS_HEL_SELECTED(H))) THEN
 
708
          IF ((.NOT.CHECKPHASE)
 
709
     $     .AND.HELDOUBLECHECKED.AND.POLARIZATIONS(0,0)
 
710
     $     .EQ.0.AND.(.NOT.ML5_0_IS_HEL_SELECTED(H))) THEN
709
711
            CYCLE
710
712
          ENDIF
711
713
 
884
886
C         Counter-term amplitude(s) for loop diagram number 33
885
887
          CALL VVV1_0(W(1,1,H),W(1,2,H),W(1,8,H),R2_3GG,AMPL(1,79))
886
888
C         Amplitude(s) for UVCT diagram with ID 40
887
 
          CALL FFV1_0(W(1,4,H),W(1,3,H),W(1,5,H),GC_5,AMPL(1,80))
888
 
          AMPL(1,80)=AMPL(1,80)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
889
 
     $     +2.0D0*UVWFCT_T_0)
 
889
          CALL FFV1_0(W(1,4,H),W(1,3,H),W(1,5,H),GC_5,AMPL(2,80))
 
890
          AMPL(2,80)=AMPL(2,80)*(4.0D0*UVWFCT_G_1_1EPS+2.0D0
 
891
     $     *UVWFCT_B_0_1EPS)
890
892
C         Amplitude(s) for UVCT diagram with ID 41
891
 
          CALL FFV1_0(W(1,4,H),W(1,3,H),W(1,5,H),GC_5,AMPL(2,81))
892
 
          AMPL(2,81)=AMPL(2,81)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
893
 
     $     *UVWFCT_G_1_1EPS)
 
893
          CALL FFV1_0(W(1,4,H),W(1,3,H),W(1,5,H),GC_5,AMPL(1,81))
 
894
          AMPL(1,81)=AMPL(1,81)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_T_0
 
895
     $     +2.0D0*UVWFCT_G_1)
894
896
C         Amplitude(s) for UVCT diagram with ID 42
895
 
          CALL FFV1_0(W(1,4,H),W(1,6,H),W(1,2,H),GC_5,AMPL(1,82))
896
 
          AMPL(1,82)=AMPL(1,82)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
897
 
     $     +2.0D0*UVWFCT_T_0)
 
897
          CALL FFV1_0(W(1,4,H),W(1,6,H),W(1,2,H),GC_5,AMPL(2,82))
 
898
          AMPL(2,82)=AMPL(2,82)*(4.0D0*UVWFCT_G_1_1EPS+2.0D0
 
899
     $     *UVWFCT_B_0_1EPS)
898
900
C         Amplitude(s) for UVCT diagram with ID 43
899
 
          CALL FFV1_0(W(1,4,H),W(1,6,H),W(1,2,H),GC_5,AMPL(2,83))
900
 
          AMPL(2,83)=AMPL(2,83)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
901
 
     $     *UVWFCT_G_1_1EPS)
 
901
          CALL FFV1_0(W(1,4,H),W(1,6,H),W(1,2,H),GC_5,AMPL(1,83))
 
902
          AMPL(1,83)=AMPL(1,83)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_T_0
 
903
     $     +2.0D0*UVWFCT_G_1)
902
904
C         Amplitude(s) for UVCT diagram with ID 44
903
 
          CALL FFV1_0(W(1,7,H),W(1,3,H),W(1,2,H),GC_5,AMPL(1,84))
904
 
          AMPL(1,84)=AMPL(1,84)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
905
 
     $     +2.0D0*UVWFCT_T_0)
 
905
          CALL FFV1_0(W(1,7,H),W(1,3,H),W(1,2,H),GC_5,AMPL(2,84))
 
906
          AMPL(2,84)=AMPL(2,84)*(4.0D0*UVWFCT_G_1_1EPS+2.0D0
 
907
     $     *UVWFCT_B_0_1EPS)
906
908
C         Amplitude(s) for UVCT diagram with ID 45
907
 
          CALL FFV1_0(W(1,7,H),W(1,3,H),W(1,2,H),GC_5,AMPL(2,85))
908
 
          AMPL(2,85)=AMPL(2,85)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
909
 
     $     *UVWFCT_G_1_1EPS)
 
909
          CALL FFV1_0(W(1,7,H),W(1,3,H),W(1,2,H),GC_5,AMPL(1,85))
 
910
          AMPL(1,85)=AMPL(1,85)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_T_0
 
911
     $     +2.0D0*UVWFCT_G_1)
910
912
 300      CONTINUE
911
913
 
912
914
 
1213
1215
            ENDDO
1214
1216
          ENDIF
1215
1217
        ELSEIF (.NOT.HELDOUBLECHECKED)THEN
1216
 
          IF ((.NOT.GOODHEL(HELPICKED)).AND.(.NOT.ML5_0_ISZERO(ABS(ANS(
1217
 
     $1))+ABS(ANS(2))+ABS(ANS(3)),REF/DBLE(NCOMB),-1))) THEN
 
1218
          IF ((.NOT.GOODHEL(HELPICKED))
 
1219
     $     .AND.(.NOT.ML5_0_ISZERO(ABS(ANS(1))+ABS(ANS(2))+ABS(ANS(3))
 
1220
     $     ,REF/DBLE(NCOMB),-1))) THEN
1218
1221
            WRITE(*,*) '##W15 Helicity filter could not be'
1219
1222
     $       //' successfully double checked.'
1220
1223
            WRITE(*,*) '##One reason for this is that you have changed'
1305
1308
 
1306
1309
        CTMODE=BASIC_CT_MODE
1307
1310
 
1308
 
        IF(.NOT.EVAL_DONE(3).AND. ((DOING_QP_EVALS.AND.NROTATIONS_QP.GE
1309
 
     $.1).OR.((.NOT.DOING_QP_EVALS).AND.NROTATIONS_DP.GE.1)) ) THEN
 
1311
        IF(.NOT.EVAL_DONE(3).AND.
 
1312
     $    ((DOING_QP_EVALS.AND.NROTATIONS_QP.GE.1)
 
1313
     $   .OR.((.NOT.DOING_QP_EVALS).AND.NROTATIONS_DP.GE.1)) ) THEN
1310
1314
          EVAL_DONE(3)=.TRUE.
1311
1315
          CALL ML5_0_ROTATE_PS(PS,P,1)
1312
1316
          IF (DOING_QP_EVALS) CALL ML5_0_MP_ROTATE_PS(MP_PS,MP_P,1)
1313
1317
          GOTO 200
1314
1318
        ENDIF
1315
1319
 
1316
 
        IF(.NOT.EVAL_DONE(4).AND. ((DOING_QP_EVALS.AND.NROTATIONS_QP.GE
1317
 
     $.2).OR.((.NOT.DOING_QP_EVALS).AND.NROTATIONS_DP.GE.2)) ) THEN
 
1320
        IF(.NOT.EVAL_DONE(4).AND.
 
1321
     $    ((DOING_QP_EVALS.AND.NROTATIONS_QP.GE.2)
 
1322
     $   .OR.((.NOT.DOING_QP_EVALS).AND.NROTATIONS_DP.GE.2)) ) THEN
1318
1323
          EVAL_DONE(4)=.TRUE.
1319
1324
          CALL ML5_0_ROTATE_PS(PS,P,2)
1320
1325
          IF (DOING_QP_EVALS) CALL ML5_0_MP_ROTATE_PS(MP_PS,MP_P,2)