~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

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

mwrge with 2.7.3

Show diffs side-by-side

added added

removed removed

Lines of Context:
102
102
      CALL MP_UPDATE_AS_PARAM()
103
103
 
104
104
      DO H=1,NCOMB
105
 
        IF ((HELPICKED.EQ.H).OR.((HELPICKED.EQ.-1).AND.((CHECKPHASE.OR.
106
 
     $.NOT.HELDOUBLECHECKED).OR.GOODHEL(H)))) THEN
 
105
        IF ((HELPICKED.EQ.H).OR.((HELPICKED.EQ.-1)
 
106
     $   .AND.((CHECKPHASE.OR..NOT.HELDOUBLECHECKED).OR.GOODHEL(H))))
 
107
     $    THEN
107
108
C         Handle the possible requirement of specific polarizations
108
 
          IF ((.NOT.CHECKPHASE).AND.HELDOUBLECHECKED.AND.POLARIZATIONS(
109
 
     $0,0).EQ.0.AND.(.NOT.ML5_0_IS_HEL_SELECTED(H))) THEN
 
109
          IF ((.NOT.CHECKPHASE)
 
110
     $     .AND.HELDOUBLECHECKED.AND.POLARIZATIONS(0,0)
 
111
     $     .EQ.0.AND.(.NOT.ML5_0_IS_HEL_SELECTED(H))) THEN
110
112
            CYCLE
111
113
          ENDIF
112
114
          DO I=1,NEXTERNAL
174
176
     $     ,AMPL(2,23))
175
177
          CALL MP_FFV1_0(W(1,10,H),W(1,2,H),W(1,5,H),R2_GQQ,AMPL(1,24))
176
178
C         Amplitude(s) for UVCT diagram with ID 14
177
 
          CALL MP_FFV2_0(W(1,6,H),W(1,2,H),W(1,7,H),GC_11,AMPL(1,25))
178
 
          AMPL(1,25)=AMPL(1,25)*(1.0D0*UVWFCT_G_2+1.0D0*UVWFCT_G_1)
 
179
          CALL MP_FFV2_0(W(1,6,H),W(1,2,H),W(1,7,H),GC_11,AMPL(2,25))
 
180
          AMPL(2,25)=AMPL(2,25)*(2.0D0*UVWFCT_G_1_1EPS)
179
181
C         Amplitude(s) for UVCT diagram with ID 15
180
 
          CALL MP_FFV2_0(W(1,6,H),W(1,2,H),W(1,7,H),GC_11,AMPL(2,26))
181
 
          AMPL(2,26)=AMPL(2,26)*(2.0D0*UVWFCT_G_1_1EPS)
 
182
          CALL MP_FFV2_0(W(1,6,H),W(1,2,H),W(1,7,H),GC_11,AMPL(1,26))
 
183
          AMPL(1,26)=AMPL(1,26)*(1.0D0*UVWFCT_G_2+1.0D0*UVWFCT_G_1)
182
184
C         Amplitude(s) for UVCT diagram with ID 16
183
 
          CALL MP_FFV2_0(W(1,1,H),W(1,8,H),W(1,7,H),GC_11,AMPL(1,27))
184
 
          AMPL(1,27)=AMPL(1,27)*(1.0D0*UVWFCT_G_2+1.0D0*UVWFCT_G_1)
 
185
          CALL MP_FFV2_0(W(1,1,H),W(1,8,H),W(1,7,H),GC_11,AMPL(2,27))
 
186
          AMPL(2,27)=AMPL(2,27)*(2.0D0*UVWFCT_G_1_1EPS)
185
187
C         Amplitude(s) for UVCT diagram with ID 17
186
 
          CALL MP_FFV2_0(W(1,1,H),W(1,8,H),W(1,7,H),GC_11,AMPL(2,28))
187
 
          AMPL(2,28)=AMPL(2,28)*(2.0D0*UVWFCT_G_1_1EPS)
 
188
          CALL MP_FFV2_0(W(1,1,H),W(1,8,H),W(1,7,H),GC_11,AMPL(1,28))
 
189
          AMPL(1,28)=AMPL(1,28)*(1.0D0*UVWFCT_G_2+1.0D0*UVWFCT_G_1)
188
190
C         Copy the qp wfs to the dp ones as they are used to setup the
189
191
C          CT calls.
190
192
          DO I=1,NWAVEFUNCS