← Back to branch summary

~maddevelopers/mg5amcnlo/new_clustering

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_wprod_fksew/%SubProcesses%P0_dxu_veep%V0_dxu_veep%mp_helas_calls_ampb_1.f

  • Committer: Rikkert Frederix
  • Date: 2021-09-09 15:51:40 UTC
  • mfrom: (78.75.502 3.2.1)
  • Revision ID: frederix@physik.uzh.ch-20210909155140-rg6umfq68h6h47cf
merge with 3.2.1

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/IOExportFKSTest/test_wprod_fksew/%SubProcesses%P0_dxu_veep%V0_dxu_veep%mp_helas_calls_ampb_1.f
 
1
      SUBROUTINE MP_HELAS_CALLS_AMPB_1(P,NHEL,H,IC)
 
2
C     
 
3
      USE POLYNOMIAL_CONSTANTS
 
4
      IMPLICIT NONE
 
5
C     
 
6
C     CONSTANTS
 
7
C     
 
8
      INTEGER    NEXTERNAL
 
9
      PARAMETER (NEXTERNAL=4)
 
10
      INTEGER    NCOMB
 
11
      PARAMETER (NCOMB=16)
 
12
 
 
13
      INTEGER NBORNAMPS
 
14
      PARAMETER (NBORNAMPS=1)
 
15
      INTEGER    NLOOPS, NLOOPGROUPS, NCTAMPS
 
16
      PARAMETER (NLOOPS=35, NLOOPGROUPS=25, NCTAMPS=15)
 
17
      INTEGER    NLOOPAMPS
 
18
      PARAMETER (NLOOPAMPS=50)
 
19
      INTEGER    NWAVEFUNCS,NLOOPWAVEFUNCS
 
20
      PARAMETER (NWAVEFUNCS=6,NLOOPWAVEFUNCS=73)
 
21
      REAL*16     ZERO
 
22
      PARAMETER (ZERO=0.0E0_16)
 
23
      COMPLEX*32     IZERO
 
24
      PARAMETER (IZERO=CMPLX(0.0E0_16,0.0E0_16,KIND=16))
 
25
C     These are constants related to the split orders
 
26
      INTEGER    NSO, NSQUAREDSO, NAMPSO
 
27
      PARAMETER (NSO=2, NSQUAREDSO=1, NAMPSO=2)
 
28
C     
 
29
C     ARGUMENTS
 
30
C     
 
31
      REAL*16 P(0:3,NEXTERNAL)
 
32
      INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
 
33
      INTEGER H
 
34
C     
 
35
C     LOCAL VARIABLES
 
36
C     
 
37
      INTEGER I,J,K
 
38
      COMPLEX*32 COEFS(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
 
39
C     
 
40
C     GLOBAL VARIABLES
 
41
C     
 
42
      INCLUDE 'mp_coupl_same_name.inc'
 
43
 
 
44
      INTEGER GOODHEL(NCOMB)
 
45
      LOGICAL GOODAMP(NSQUAREDSO,NLOOPGROUPS)
 
46
      COMMON/FILTERS/GOODAMP,GOODHEL
 
47
 
 
48
      INTEGER SQSO_TARGET
 
49
      COMMON/SOCHOICE/SQSO_TARGET
 
50
 
 
51
      LOGICAL UVCT_REQ_SO_DONE,MP_UVCT_REQ_SO_DONE,CT_REQ_SO_DONE
 
52
     $ ,MP_CT_REQ_SO_DONE,LOOP_REQ_SO_DONE,MP_LOOP_REQ_SO_DONE
 
53
     $ ,CTCALL_REQ_SO_DONE,FILTER_SO
 
54
      COMMON/SO_REQS/UVCT_REQ_SO_DONE,MP_UVCT_REQ_SO_DONE
 
55
     $ ,CT_REQ_SO_DONE,MP_CT_REQ_SO_DONE,LOOP_REQ_SO_DONE
 
56
     $ ,MP_LOOP_REQ_SO_DONE,CTCALL_REQ_SO_DONE,FILTER_SO
 
57
 
 
58
      COMPLEX*32 AMP(NBORNAMPS)
 
59
      COMMON/MP_AMPS/AMP
 
60
      COMPLEX*32 W(20,NWAVEFUNCS)
 
61
      COMMON/MP_W/W
 
62
 
 
63
      COMPLEX*32 WL(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE,
 
64
     $ -1:NLOOPWAVEFUNCS)
 
65
      COMPLEX*32 PL(0:3,-1:NLOOPWAVEFUNCS)
 
66
      COMMON/MP_WL/WL,PL
 
67
 
 
68
      COMPLEX*32 AMPL(3,NCTAMPS)
 
69
      COMMON/MP_AMPL/AMPL
 
70
 
 
71
C     
 
72
C     ----------
 
73
C     BEGIN CODE
 
74
C     ----------
 
75
 
 
76
C     The target squared split order contribution is already reached
 
77
C      if true.
 
78
      IF (FILTER_SO.AND.MP_CT_REQ_SO_DONE) THEN
 
79
        GOTO 1001
 
80
      ENDIF
 
81
 
 
82
      CALL MP_OXXXXX(P(0,1),ZERO,NHEL(1),-1*IC(1),W(1,1))
 
83
      CALL MP_IXXXXX(P(0,2),ZERO,NHEL(2),+1*IC(2),W(1,2))
 
84
      CALL MP_OXXXXX(P(0,3),ZERO,NHEL(3),+1*IC(3),W(1,3))
 
85
      CALL MP_IXXXXX(P(0,4),ZERO,NHEL(4),-1*IC(4),W(1,4))
 
86
      CALL MP_FFV2P0_3(W(1,2),W(1,1),GC_124,CMPLX(CMASS_MDL_MW,KIND=16)
 
87
     $ ,W(1,5))
 
88
C     Amplitude(s) for born diagram with ID 1
 
89
      CALL MP_FFV2_0(W(1,4),W(1,3),W(1,5),GC_124,AMP(1))
 
90
      CALL MP_FFV2P0_3(W(1,4),W(1,3),GC_124,CMPLX(CMASS_MDL_MW,KIND=16)
 
91
     $ ,W(1,6))
 
92
C     Counter-term amplitude(s) for loop diagram number 2
 
93
      CALL MP_L_WMWPMASS2_L_WMWPMASS4_0(W(1,6),W(1,5)
 
94
     $ ,C_UVWMWPMASS2EW_1EPS,C_UVWMWPMASS1EW_1EPS,AMPL(2,1))
 
95
      CALL MP_L_WMWPMASS2_L_WMWPMASS4_0(W(1,6),W(1,5),C_UVWMWPMASS2EW
 
96
     $ ,C_UVWMWPMASS1EW,AMPL(1,2))
 
97
C     Counter-term amplitude(s) for loop diagram number 4
 
98
      CALL MP_L_VEXVEA21_0(W(1,4),W(1,3),W(1,5),C_UVEPVEWM1EW_1EPS
 
99
     $ ,AMPL(2,3))
 
100
      CALL MP_L_VEXVEA21_0(W(1,4),W(1,3),W(1,5),C_UVEPVEWM1EW,AMPL(1,4)
 
101
     $ )
 
102
C     Counter-term amplitude(s) for loop diagram number 5
 
103
      CALL MP_L_VEXVEA21_0(W(1,2),W(1,1),W(1,6),C_UVCXSWP1EW_1EPS
 
104
     $ ,AMPL(2,5))
 
105
      CALL MP_L_VEXVEA21_0(W(1,2),W(1,1),W(1,6),C_UVCXSWP1EW,AMPL(1,6))
 
106
C     Counter-term amplitude(s) for loop diagram number 14
 
107
      CALL MP_FFV2_0(W(1,4),W(1,3),W(1,5),R2_VLW,AMPL(1,7))
 
108
C     Counter-term amplitude(s) for loop diagram number 16
 
109
      CALL MP_FFV2_0(W(1,2),W(1,1),W(1,6),R2_BXTW2CP,AMPL(1,8))
 
110
C     Counter-term amplitude(s) for loop diagram number 23
 
111
      CALL MP_R2_GG_1_R2_GG_2_R2_GG_3_0(W(1,5),W(1,6),R2_WWBOSON1
 
112
     $ ,R2_WWBOSON2,R2_WWBOSON3,AMPL(1,9))
 
113
C     Counter-term amplitude(s) for loop diagram number 34
 
114
      CALL MP_R2_GG_1_0(W(1,5),W(1,6),R2_WWCS1,AMPL(1,10))
 
115
      CALL MP_R2_GG_1_0(W(1,5),W(1,6),R2_WWCS1,AMPL(1,11))
 
116
C     Counter-term amplitude(s) for loop diagram number 35
 
117
      CALL MP_R2_GG_1_R2_GG_3_0(W(1,5),W(1,6),R2_WWCS1,R2_WWTB3,AMPL(1
 
118
     $ ,12))
 
119
C     Counter-term amplitude(s) for loop diagram number 36
 
120
      CALL MP_R2_GG_1_0(W(1,5),W(1,6),R2_WWL,AMPL(1,13))
 
121
      CALL MP_R2_GG_1_0(W(1,5),W(1,6),R2_WWL,AMPL(1,14))
 
122
      CALL MP_R2_GG_1_0(W(1,5),W(1,6),R2_WWL,AMPL(1,15))
 
123
C     At this point, all CT amps needed for (QCD=0 QED=6), i.e. of
 
124
C      split order ID=1, are computed.
 
125
      IF(FILTER_SO.AND.SQSO_TARGET.EQ.1) GOTO 2000
 
126
 
 
127
      GOTO 1001
 
128
 2000 CONTINUE
 
129
      MP_CT_REQ_SO_DONE=.TRUE.
 
130
 1001 CONTINUE
 
131
      END
 
132