← 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/gg_wmtbx/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/gg_wmtbx/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
 
1195
1197
C         Counter-term amplitude(s) for loop diagram number 122
1196
1198
          CALL VVV1_0(W(1,1,H),W(1,2,H),W(1,14,H),R2_3GG,AMPL(1,236))
1197
1199
C         Amplitude(s) for UVCT diagram with ID 135
1198
 
          CALL FFV1_0(W(1,5,H),W(1,7,H),W(1,6,H),GC_5,AMPL(1,237))
1199
 
          AMPL(1,237)=AMPL(1,237)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
1200
 
     $     +1.0D0*UVWFCT_T_0+1.0D0*UVWFCT_B_0)
 
1200
          CALL FFV1_0(W(1,5,H),W(1,7,H),W(1,6,H),GC_5,AMPL(2,237))
 
1201
          AMPL(2,237)=AMPL(2,237)*(4.0D0*UVWFCT_G_1_1EPS+2.0D0
 
1202
     $     *UVWFCT_B_0_1EPS)
1201
1203
C         Amplitude(s) for UVCT diagram with ID 136
1202
 
          CALL FFV1_0(W(1,5,H),W(1,7,H),W(1,6,H),GC_5,AMPL(2,238))
1203
 
          AMPL(2,238)=AMPL(2,238)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
1204
 
     $     *UVWFCT_G_1_1EPS)
 
1204
          CALL FFV1_0(W(1,5,H),W(1,7,H),W(1,6,H),GC_5,AMPL(1,238))
 
1205
          AMPL(1,238)=AMPL(1,238)*(2.0D0*UVWFCT_G_2+1.0D0*UVWFCT_T_0
 
1206
     $     +1.0D0*UVWFCT_B_0+2.0D0*UVWFCT_G_1)
1205
1207
C         Amplitude(s) for UVCT diagram with ID 137
1206
 
          CALL FFV1_0(W(1,8,H),W(1,4,H),W(1,6,H),GC_5,AMPL(1,239))
1207
 
          AMPL(1,239)=AMPL(1,239)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
1208
 
     $     +1.0D0*UVWFCT_T_0+1.0D0*UVWFCT_B_0)
 
1208
          CALL FFV1_0(W(1,8,H),W(1,4,H),W(1,6,H),GC_5,AMPL(2,239))
 
1209
          AMPL(2,239)=AMPL(2,239)*(4.0D0*UVWFCT_G_1_1EPS+2.0D0
 
1210
     $     *UVWFCT_B_0_1EPS)
1209
1211
C         Amplitude(s) for UVCT diagram with ID 138
1210
 
          CALL FFV1_0(W(1,8,H),W(1,4,H),W(1,6,H),GC_5,AMPL(2,240))
1211
 
          AMPL(2,240)=AMPL(2,240)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
1212
 
     $     *UVWFCT_G_1_1EPS)
 
1212
          CALL FFV1_0(W(1,8,H),W(1,4,H),W(1,6,H),GC_5,AMPL(1,240))
 
1213
          AMPL(1,240)=AMPL(1,240)*(2.0D0*UVWFCT_G_2+1.0D0*UVWFCT_T_0
 
1214
     $     +1.0D0*UVWFCT_B_0+2.0D0*UVWFCT_G_1)
1213
1215
C         Amplitude(s) for UVCT diagram with ID 139
1214
 
          CALL FFV2_0(W(1,10,H),W(1,9,H),W(1,3,H),GC_11,AMPL(1,241))
1215
 
          AMPL(1,241)=AMPL(1,241)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
1216
 
     $     +1.0D0*UVWFCT_T_0+1.0D0*UVWFCT_B_0)
 
1216
          CALL FFV2_0(W(1,10,H),W(1,9,H),W(1,3,H),GC_11,AMPL(2,241))
 
1217
          AMPL(2,241)=AMPL(2,241)*(4.0D0*UVWFCT_G_1_1EPS+2.0D0
 
1218
     $     *UVWFCT_B_0_1EPS)
1217
1219
C         Amplitude(s) for UVCT diagram with ID 140
1218
 
          CALL FFV2_0(W(1,10,H),W(1,9,H),W(1,3,H),GC_11,AMPL(2,242))
1219
 
          AMPL(2,242)=AMPL(2,242)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
1220
 
     $     *UVWFCT_G_1_1EPS)
 
1220
          CALL FFV2_0(W(1,10,H),W(1,9,H),W(1,3,H),GC_11,AMPL(1,242))
 
1221
          AMPL(1,242)=AMPL(1,242)*(2.0D0*UVWFCT_G_2+1.0D0*UVWFCT_T_0
 
1222
     $     +1.0D0*UVWFCT_B_0+2.0D0*UVWFCT_G_1)
1221
1223
C         Amplitude(s) for UVCT diagram with ID 141
1222
 
          CALL FFV1_0(W(1,8,H),W(1,9,H),W(1,2,H),GC_5,AMPL(1,243))
1223
 
          AMPL(1,243)=AMPL(1,243)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
1224
 
     $     +1.0D0*UVWFCT_T_0+1.0D0*UVWFCT_B_0)
 
1224
          CALL FFV1_0(W(1,8,H),W(1,9,H),W(1,2,H),GC_5,AMPL(2,243))
 
1225
          AMPL(2,243)=AMPL(2,243)*(4.0D0*UVWFCT_G_1_1EPS+2.0D0
 
1226
     $     *UVWFCT_B_0_1EPS)
1225
1227
C         Amplitude(s) for UVCT diagram with ID 142
1226
 
          CALL FFV1_0(W(1,8,H),W(1,9,H),W(1,2,H),GC_5,AMPL(2,244))
1227
 
          AMPL(2,244)=AMPL(2,244)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
1228
 
     $     *UVWFCT_G_1_1EPS)
 
1228
          CALL FFV1_0(W(1,8,H),W(1,9,H),W(1,2,H),GC_5,AMPL(1,244))
 
1229
          AMPL(1,244)=AMPL(1,244)*(2.0D0*UVWFCT_G_2+1.0D0*UVWFCT_T_0
 
1230
     $     +1.0D0*UVWFCT_B_0+2.0D0*UVWFCT_G_1)
1229
1231
C         Amplitude(s) for UVCT diagram with ID 143
1230
 
          CALL FFV2_0(W(1,11,H),W(1,12,H),W(1,3,H),GC_11,AMPL(1,245))
1231
 
          AMPL(1,245)=AMPL(1,245)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
1232
 
     $     +1.0D0*UVWFCT_T_0+1.0D0*UVWFCT_B_0)
 
1232
          CALL FFV2_0(W(1,11,H),W(1,12,H),W(1,3,H),GC_11,AMPL(2,245))
 
1233
          AMPL(2,245)=AMPL(2,245)*(4.0D0*UVWFCT_G_1_1EPS+2.0D0
 
1234
     $     *UVWFCT_B_0_1EPS)
1233
1235
C         Amplitude(s) for UVCT diagram with ID 144
1234
 
          CALL FFV2_0(W(1,11,H),W(1,12,H),W(1,3,H),GC_11,AMPL(2,246))
1235
 
          AMPL(2,246)=AMPL(2,246)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
1236
 
     $     *UVWFCT_G_1_1EPS)
 
1236
          CALL FFV2_0(W(1,11,H),W(1,12,H),W(1,3,H),GC_11,AMPL(1,246))
 
1237
          AMPL(1,246)=AMPL(1,246)*(2.0D0*UVWFCT_G_2+1.0D0*UVWFCT_T_0
 
1238
     $     +1.0D0*UVWFCT_B_0+2.0D0*UVWFCT_G_1)
1237
1239
C         Amplitude(s) for UVCT diagram with ID 145
1238
 
          CALL FFV1_0(W(1,11,H),W(1,7,H),W(1,2,H),GC_5,AMPL(1,247))
1239
 
          AMPL(1,247)=AMPL(1,247)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
1240
 
     $     +1.0D0*UVWFCT_T_0+1.0D0*UVWFCT_B_0)
 
1240
          CALL FFV1_0(W(1,11,H),W(1,7,H),W(1,2,H),GC_5,AMPL(2,247))
 
1241
          AMPL(2,247)=AMPL(2,247)*(4.0D0*UVWFCT_G_1_1EPS+2.0D0
 
1242
     $     *UVWFCT_B_0_1EPS)
1241
1243
C         Amplitude(s) for UVCT diagram with ID 146
1242
 
          CALL FFV1_0(W(1,11,H),W(1,7,H),W(1,2,H),GC_5,AMPL(2,248))
1243
 
          AMPL(2,248)=AMPL(2,248)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
1244
 
     $     *UVWFCT_G_1_1EPS)
 
1244
          CALL FFV1_0(W(1,11,H),W(1,7,H),W(1,2,H),GC_5,AMPL(1,248))
 
1245
          AMPL(1,248)=AMPL(1,248)*(2.0D0*UVWFCT_G_2+1.0D0*UVWFCT_T_0
 
1246
     $     +1.0D0*UVWFCT_B_0+2.0D0*UVWFCT_G_1)
1245
1247
C         Amplitude(s) for UVCT diagram with ID 147
1246
 
          CALL FFV1_0(W(1,8,H),W(1,12,H),W(1,1,H),GC_5,AMPL(1,249))
1247
 
          AMPL(1,249)=AMPL(1,249)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
1248
 
     $     +1.0D0*UVWFCT_T_0+1.0D0*UVWFCT_B_0)
 
1248
          CALL FFV1_0(W(1,8,H),W(1,12,H),W(1,1,H),GC_5,AMPL(2,249))
 
1249
          AMPL(2,249)=AMPL(2,249)*(4.0D0*UVWFCT_G_1_1EPS+2.0D0
 
1250
     $     *UVWFCT_B_0_1EPS)
1249
1251
C         Amplitude(s) for UVCT diagram with ID 148
1250
 
          CALL FFV1_0(W(1,8,H),W(1,12,H),W(1,1,H),GC_5,AMPL(2,250))
1251
 
          AMPL(2,250)=AMPL(2,250)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
1252
 
     $     *UVWFCT_G_1_1EPS)
 
1252
          CALL FFV1_0(W(1,8,H),W(1,12,H),W(1,1,H),GC_5,AMPL(1,250))
 
1253
          AMPL(1,250)=AMPL(1,250)*(2.0D0*UVWFCT_G_2+1.0D0*UVWFCT_T_0
 
1254
     $     +1.0D0*UVWFCT_B_0+2.0D0*UVWFCT_G_1)
1253
1255
C         Amplitude(s) for UVCT diagram with ID 149
1254
 
          CALL FFV1_0(W(1,10,H),W(1,7,H),W(1,1,H),GC_5,AMPL(1,251))
1255
 
          AMPL(1,251)=AMPL(1,251)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
1256
 
     $     +1.0D0*UVWFCT_T_0+1.0D0*UVWFCT_B_0)
 
1256
          CALL FFV1_0(W(1,10,H),W(1,7,H),W(1,1,H),GC_5,AMPL(2,251))
 
1257
          AMPL(2,251)=AMPL(2,251)*(4.0D0*UVWFCT_G_1_1EPS+2.0D0
 
1258
     $     *UVWFCT_B_0_1EPS)
1257
1259
C         Amplitude(s) for UVCT diagram with ID 150
1258
 
          CALL FFV1_0(W(1,10,H),W(1,7,H),W(1,1,H),GC_5,AMPL(2,252))
1259
 
          AMPL(2,252)=AMPL(2,252)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
1260
 
     $     *UVWFCT_G_1_1EPS)
 
1260
          CALL FFV1_0(W(1,10,H),W(1,7,H),W(1,1,H),GC_5,AMPL(1,252))
 
1261
          AMPL(1,252)=AMPL(1,252)*(2.0D0*UVWFCT_G_2+1.0D0*UVWFCT_T_0
 
1262
     $     +1.0D0*UVWFCT_B_0+2.0D0*UVWFCT_G_1)
1261
1263
 300      CONTINUE
1262
1264
 
1263
1265
 
2076
2078
            ENDDO
2077
2079
          ENDIF
2078
2080
        ELSEIF (.NOT.HELDOUBLECHECKED)THEN
2079
 
          IF ((.NOT.GOODHEL(HELPICKED)).AND.(.NOT.ML5_0_ISZERO(ABS(ANS(
2080
 
     $1))+ABS(ANS(2))+ABS(ANS(3)),REF/DBLE(NCOMB),-1))) THEN
 
2081
          IF ((.NOT.GOODHEL(HELPICKED))
 
2082
     $     .AND.(.NOT.ML5_0_ISZERO(ABS(ANS(1))+ABS(ANS(2))+ABS(ANS(3))
 
2083
     $     ,REF/DBLE(NCOMB),-1))) THEN
2081
2084
            WRITE(*,*) '##W15 Helicity filter could not be'
2082
2085
     $       //' successfully double checked.'
2083
2086
            WRITE(*,*) '##One reason for this is that you have changed'
2168
2171
 
2169
2172
        CTMODE=BASIC_CT_MODE
2170
2173
 
2171
 
        IF(.NOT.EVAL_DONE(3).AND. ((DOING_QP_EVALS.AND.NROTATIONS_QP.GE
2172
 
     $.1).OR.((.NOT.DOING_QP_EVALS).AND.NROTATIONS_DP.GE.1)) ) THEN
 
2174
        IF(.NOT.EVAL_DONE(3).AND.
 
2175
     $    ((DOING_QP_EVALS.AND.NROTATIONS_QP.GE.1)
 
2176
     $   .OR.((.NOT.DOING_QP_EVALS).AND.NROTATIONS_DP.GE.1)) ) THEN
2173
2177
          EVAL_DONE(3)=.TRUE.
2174
2178
          CALL ML5_0_ROTATE_PS(PS,P,1)
2175
2179
          IF (DOING_QP_EVALS) CALL ML5_0_MP_ROTATE_PS(MP_PS,MP_P,1)
2176
2180
          GOTO 200
2177
2181
        ENDIF
2178
2182
 
2179
 
        IF(.NOT.EVAL_DONE(4).AND. ((DOING_QP_EVALS.AND.NROTATIONS_QP.GE
2180
 
     $.2).OR.((.NOT.DOING_QP_EVALS).AND.NROTATIONS_DP.GE.2)) ) THEN
 
2183
        IF(.NOT.EVAL_DONE(4).AND.
 
2184
     $    ((DOING_QP_EVALS.AND.NROTATIONS_QP.GE.2)
 
2185
     $   .OR.((.NOT.DOING_QP_EVALS).AND.NROTATIONS_DP.GE.2)) ) THEN
2181
2186
          EVAL_DONE(4)=.TRUE.
2182
2187
          CALL ML5_0_ROTATE_PS(PS,P,2)
2183
2188
          IF (DOING_QP_EVALS) CALL ML5_0_MP_ROTATE_PS(MP_PS,MP_P,2)