← 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/mp_born_amps_and_wfs.f

  • Committer: olivier Mattelaer
  • Date: 2016-05-12 11:00:18 UTC
  • mfrom: (262.1.150 2.3.4)
  • Revision ID: olivier.mattelaer@uclouvain.be-20160512110018-sevb79f0wm4g8mpp
pass to 2.4.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/short_ML_SMQCD_default/gg_ttx/mp_born_amps_and_wfs.f
7
7
C     Computes all the AMP and WFS in quadruple precision for the 
8
8
C     phase space point P(0:3,NEXTERNAL)
9
9
C     
10
 
C     Process: g g > t t~ QED=0 QCD=2 [ virt = QCD ]
 
10
C     Process: g g > t t~ QED=0 QCD<=2 [ virt = QCD ]
11
11
C     
12
12
      IMPLICIT NONE
13
13
C     
93
93
 
94
94
      DO H=1,NCOMB
95
95
        IF ((HELPICKED.EQ.H).OR.((HELPICKED.EQ.-1).AND.((CHECKPHASE.OR.
96
 
     $   .NOT.HELDOUBLECHECKED).OR.GOODHEL(H)))) THEN
 
96
     $.NOT.HELDOUBLECHECKED).OR.GOODHEL(H)))) THEN
97
97
          DO I=1,NEXTERNAL
98
98
            NHEL(I)=HELC(I,H)
99
99
          ENDDO
187
187
          CALL MP_FFV1_0(W(1,4,H),W(1,6,H),W(1,2,H),UV_GQQG_1EPS
188
188
     $     ,AMPL(2,42))
189
189
          CALL MP_FFV1_0(W(1,4,H),W(1,6,H),W(1,2,H),R2_GQQ,AMPL(1,43))
190
 
          CALL MP_FFV1_1(W(1,3,H),W(1,2,H),GC_5,MDL_MT,MDL_WT,W(1,10
191
 
     $     ,H))
 
190
          CALL MP_FFV1_1(W(1,3,H),W(1,2,H),GC_5,MDL_MT,MDL_WT,W(1,10,H)
 
191
     $     )
192
192
C         Counter-term amplitude(s) for loop diagram number 16
193
193
          CALL MP_R2_QQ_1_R2_QQ_2_0(W(1,7,H),W(1,10,H),R2_QQQ,R2_QQT
194
194
     $     ,AMPL(1,44))
195
195
          CALL MP_R2_QQ_2_0(W(1,7,H),W(1,10,H),UV_TMASS,AMPL(1,45))
196
 
          CALL MP_R2_QQ_2_0(W(1,7,H),W(1,10,H),UV_TMASS_1EPS,AMPL(2
197
 
     $     ,46))
 
196
          CALL MP_R2_QQ_2_0(W(1,7,H),W(1,10,H),UV_TMASS_1EPS,AMPL(2,46)
 
197
     $     )
198
198
C         Counter-term amplitude(s) for loop diagram number 17
199
199
          CALL MP_FFV1_0(W(1,7,H),W(1,3,H),W(1,2,H),UV_GQQQ_1EPS
200
200
     $     ,AMPL(2,47))
222
222
     $     ,AMPL(2,59))
223
223
          CALL MP_FFV1_0(W(1,4,H),W(1,10,H),W(1,1,H),UV_GQQQ_1EPS
224
224
     $     ,AMPL(2,60))
225
 
          CALL MP_FFV1_0(W(1,4,H),W(1,10,H),W(1,1,H),UV_GQQB,AMPL(1
226
 
     $     ,61))
 
225
          CALL MP_FFV1_0(W(1,4,H),W(1,10,H),W(1,1,H),UV_GQQB,AMPL(1,61)
 
226
     $     )
227
227
          CALL MP_FFV1_0(W(1,4,H),W(1,10,H),W(1,1,H),UV_GQQQ_1EPS
228
228
     $     ,AMPL(2,62))
229
 
          CALL MP_FFV1_0(W(1,4,H),W(1,10,H),W(1,1,H),UV_GQQT,AMPL(1
230
 
     $     ,63))
 
229
          CALL MP_FFV1_0(W(1,4,H),W(1,10,H),W(1,1,H),UV_GQQT,AMPL(1,63)
 
230
     $     )
231
231
          CALL MP_FFV1_0(W(1,4,H),W(1,10,H),W(1,1,H),UV_GQQQ_1EPS
232
232
     $     ,AMPL(2,64))
233
233
          CALL MP_FFV1_0(W(1,4,H),W(1,10,H),W(1,1,H),UV_GQQG_1EPS
264
264
     $     +2.0D0*UVWFCT_T_0)
265
265
C         Amplitude(s) for UVCT diagram with ID 41
266
266
          CALL MP_FFV1_0(W(1,4,H),W(1,3,H),W(1,5,H),GC_5,AMPL(2,81))
267
 
          AMPL(2,81)=AMPL(2,81)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0*UVWFCT_G_2
268
 
     $     _1EPS)
 
267
          AMPL(2,81)=AMPL(2,81)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
 
268
     $     *UVWFCT_G_2_1EPS)
269
269
C         Amplitude(s) for UVCT diagram with ID 42
270
270
          CALL MP_FFV1_0(W(1,4,H),W(1,6,H),W(1,2,H),GC_5,AMPL(1,82))
271
271
          AMPL(1,82)=AMPL(1,82)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
272
272
     $     +2.0D0*UVWFCT_T_0)
273
273
C         Amplitude(s) for UVCT diagram with ID 43
274
274
          CALL MP_FFV1_0(W(1,4,H),W(1,6,H),W(1,2,H),GC_5,AMPL(2,83))
275
 
          AMPL(2,83)=AMPL(2,83)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0*UVWFCT_G_2
276
 
     $     _1EPS)
 
275
          AMPL(2,83)=AMPL(2,83)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
 
276
     $     *UVWFCT_G_2_1EPS)
277
277
C         Amplitude(s) for UVCT diagram with ID 44
278
278
          CALL MP_FFV1_0(W(1,7,H),W(1,3,H),W(1,2,H),GC_5,AMPL(1,84))
279
279
          AMPL(1,84)=AMPL(1,84)*(2.0D0*UVWFCT_G_2+2.0D0*UVWFCT_G_1
280
280
     $     +2.0D0*UVWFCT_T_0)
281
281
C         Amplitude(s) for UVCT diagram with ID 45
282
282
          CALL MP_FFV1_0(W(1,7,H),W(1,3,H),W(1,2,H),GC_5,AMPL(2,85))
283
 
          AMPL(2,85)=AMPL(2,85)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0*UVWFCT_G_2
284
 
     $     _1EPS)
 
283
          AMPL(2,85)=AMPL(2,85)*(2.0D0*UVWFCT_B_0_1EPS+4.0D0
 
284
     $     *UVWFCT_G_2_1EPS)
285
285
C         Copy the qp wfs to the dp ones as they are used to setup the
286
286
C          CT calls.
287
287
          DO I=1,NWAVEFUNCS