← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/long_ML_SMQCD_default/dux_mumvmxg/loop_matrix.f

mwrge with 2.7.3

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/long_ML_SMQCD_default/dux_mumvmxg/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
 
789
791
     $     ,23))
790
792
          CALL FFV1_0(W(1,10,H),W(1,2,H),W(1,5,H),R2_GQQ,AMPL(1,24))
791
793
C         Amplitude(s) for UVCT diagram with ID 14
792
 
          CALL FFV2_0(W(1,6,H),W(1,2,H),W(1,7,H),GC_11,AMPL(1,25))
793
 
          AMPL(1,25)=AMPL(1,25)*(1.0D0*UVWFCT_G_2+1.0D0*UVWFCT_G_1)
 
794
          CALL FFV2_0(W(1,6,H),W(1,2,H),W(1,7,H),GC_11,AMPL(2,25))
 
795
          AMPL(2,25)=AMPL(2,25)*(2.0D0*UVWFCT_G_1_1EPS)
794
796
C         Amplitude(s) for UVCT diagram with ID 15
795
 
          CALL FFV2_0(W(1,6,H),W(1,2,H),W(1,7,H),GC_11,AMPL(2,26))
796
 
          AMPL(2,26)=AMPL(2,26)*(2.0D0*UVWFCT_G_1_1EPS)
 
797
          CALL FFV2_0(W(1,6,H),W(1,2,H),W(1,7,H),GC_11,AMPL(1,26))
 
798
          AMPL(1,26)=AMPL(1,26)*(1.0D0*UVWFCT_G_2+1.0D0*UVWFCT_G_1)
797
799
C         Amplitude(s) for UVCT diagram with ID 16
798
 
          CALL FFV2_0(W(1,1,H),W(1,8,H),W(1,7,H),GC_11,AMPL(1,27))
799
 
          AMPL(1,27)=AMPL(1,27)*(1.0D0*UVWFCT_G_2+1.0D0*UVWFCT_G_1)
 
800
          CALL FFV2_0(W(1,1,H),W(1,8,H),W(1,7,H),GC_11,AMPL(2,27))
 
801
          AMPL(2,27)=AMPL(2,27)*(2.0D0*UVWFCT_G_1_1EPS)
800
802
C         Amplitude(s) for UVCT diagram with ID 17
801
 
          CALL FFV2_0(W(1,1,H),W(1,8,H),W(1,7,H),GC_11,AMPL(2,28))
802
 
          AMPL(2,28)=AMPL(2,28)*(2.0D0*UVWFCT_G_1_1EPS)
 
803
          CALL FFV2_0(W(1,1,H),W(1,8,H),W(1,7,H),GC_11,AMPL(1,28))
 
804
          AMPL(1,28)=AMPL(1,28)*(1.0D0*UVWFCT_G_2+1.0D0*UVWFCT_G_1)
803
805
 300      CONTINUE
804
806
 
805
807
 
960
962
            ENDDO
961
963
          ENDIF
962
964
        ELSEIF (.NOT.HELDOUBLECHECKED)THEN
963
 
          IF ((.NOT.GOODHEL(HELPICKED)).AND.(.NOT.ML5_0_ISZERO(ABS(ANS(
964
 
     $1))+ABS(ANS(2))+ABS(ANS(3)),REF/DBLE(NCOMB),-1))) THEN
 
965
          IF ((.NOT.GOODHEL(HELPICKED))
 
966
     $     .AND.(.NOT.ML5_0_ISZERO(ABS(ANS(1))+ABS(ANS(2))+ABS(ANS(3))
 
967
     $     ,REF/DBLE(NCOMB),-1))) THEN
965
968
            WRITE(*,*) '##W15 Helicity filter could not be'
966
969
     $       //' successfully double checked.'
967
970
            WRITE(*,*) '##One reason for this is that you have changed'
1052
1055
 
1053
1056
        CTMODE=BASIC_CT_MODE
1054
1057
 
1055
 
        IF(.NOT.EVAL_DONE(3).AND. ((DOING_QP_EVALS.AND.NROTATIONS_QP.GE
1056
 
     $.1).OR.((.NOT.DOING_QP_EVALS).AND.NROTATIONS_DP.GE.1)) ) THEN
 
1058
        IF(.NOT.EVAL_DONE(3).AND.
 
1059
     $    ((DOING_QP_EVALS.AND.NROTATIONS_QP.GE.1)
 
1060
     $   .OR.((.NOT.DOING_QP_EVALS).AND.NROTATIONS_DP.GE.1)) ) THEN
1057
1061
          EVAL_DONE(3)=.TRUE.
1058
1062
          CALL ML5_0_ROTATE_PS(PS,P,1)
1059
1063
          IF (DOING_QP_EVALS) CALL ML5_0_MP_ROTATE_PS(MP_PS,MP_P,1)
1060
1064
          GOTO 200
1061
1065
        ENDIF
1062
1066
 
1063
 
        IF(.NOT.EVAL_DONE(4).AND. ((DOING_QP_EVALS.AND.NROTATIONS_QP.GE
1064
 
     $.2).OR.((.NOT.DOING_QP_EVALS).AND.NROTATIONS_DP.GE.2)) ) THEN
 
1067
        IF(.NOT.EVAL_DONE(4).AND.
 
1068
     $    ((DOING_QP_EVALS.AND.NROTATIONS_QP.GE.2)
 
1069
     $   .OR.((.NOT.DOING_QP_EVALS).AND.NROTATIONS_DP.GE.2)) ) THEN
1065
1070
          EVAL_DONE(4)=.TRUE.
1066
1071
          CALL ML5_0_ROTATE_PS(PS,P,2)
1067
1072
          IF (DOING_QP_EVALS) CALL ML5_0_MP_ROTATE_PS(MP_PS,MP_P,2)