← Back to branch summary

~maddevelopers/mg5amcnlo/2.7.1.3

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_udx_wp%V0_udx_wp%helas_calls_ampb_1.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/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_udx_wp%V0_udx_wp%helas_calls_ampb_1.f
 
1
      SUBROUTINE HELAS_CALLS_AMPB_1(P,NHEL,H,IC)
 
2
C     
 
3
C     Modules
 
4
C     
 
5
      USE POLYNOMIAL_CONSTANTS
 
6
C     
 
7
      IMPLICIT NONE
 
8
C     
 
9
C     CONSTANTS
 
10
C     
 
11
      INTEGER    NEXTERNAL
 
12
      PARAMETER (NEXTERNAL=3)
 
13
      INTEGER    NCOMB
 
14
      PARAMETER (NCOMB=12)
 
15
      INTEGER NBORNAMPS
 
16
      PARAMETER (NBORNAMPS=1)
 
17
      INTEGER    NLOOPS, NLOOPGROUPS, NCTAMPS
 
18
      PARAMETER (NLOOPS=1, NLOOPGROUPS=1, NCTAMPS=1)
 
19
      INTEGER    NLOOPAMPS
 
20
      PARAMETER (NLOOPAMPS=2)
 
21
      INTEGER    NWAVEFUNCS,NLOOPWAVEFUNCS
 
22
      PARAMETER (NWAVEFUNCS=3,NLOOPWAVEFUNCS=3)
 
23
      REAL*8     ZERO
 
24
      PARAMETER (ZERO=0D0)
 
25
      REAL*16     MP__ZERO
 
26
      PARAMETER (MP__ZERO=0.0E0_16)
 
27
C     These are constants related to the split orders
 
28
      INTEGER    NSO, NSQUAREDSO, NAMPSO
 
29
      PARAMETER (NSO=1, NSQUAREDSO=1, NAMPSO=2)
 
30
C     
 
31
C     ARGUMENTS
 
32
C     
 
33
      REAL*8 P(0:3,NEXTERNAL)
 
34
      INTEGER NHEL(NEXTERNAL), IC(NEXTERNAL)
 
35
      INTEGER H
 
36
C     
 
37
C     LOCAL VARIABLES
 
38
C     
 
39
      INTEGER I,J,K
 
40
      COMPLEX*16 COEFS(MAXLWFSIZE,0:VERTEXMAXCOEFS-1,MAXLWFSIZE)
 
41
 
 
42
      LOGICAL DUMMYFALSE
 
43
      DATA DUMMYFALSE/.FALSE./
 
44
C     
 
45
C     GLOBAL VARIABLES
 
46
C     
 
47
      INCLUDE 'coupl.inc'
 
48
      INCLUDE 'mp_coupl.inc'
 
49
 
 
50
      INTEGER HELOFFSET
 
51
      INTEGER GOODHEL(NCOMB)
 
52
      LOGICAL GOODAMP(NSQUAREDSO,NLOOPGROUPS)
 
53
      COMMON/FILTERS/GOODAMP,GOODHEL,HELOFFSET
 
54
 
 
55
      LOGICAL CHECKPHASE
 
56
      LOGICAL HELDOUBLECHECKED
 
57
      COMMON/INIT/CHECKPHASE, HELDOUBLECHECKED
 
58
 
 
59
      INTEGER SQSO_TARGET
 
60
      COMMON/SOCHOICE/SQSO_TARGET
 
61
 
 
62
      LOGICAL UVCT_REQ_SO_DONE,MP_UVCT_REQ_SO_DONE,CT_REQ_SO_DONE
 
63
     $ ,MP_CT_REQ_SO_DONE,LOOP_REQ_SO_DONE,MP_LOOP_REQ_SO_DONE
 
64
     $ ,CTCALL_REQ_SO_DONE,FILTER_SO
 
65
      COMMON/SO_REQS/UVCT_REQ_SO_DONE,MP_UVCT_REQ_SO_DONE
 
66
     $ ,CT_REQ_SO_DONE,MP_CT_REQ_SO_DONE,LOOP_REQ_SO_DONE
 
67
     $ ,MP_LOOP_REQ_SO_DONE,CTCALL_REQ_SO_DONE,FILTER_SO
 
68
 
 
69
      INTEGER I_SO
 
70
      COMMON/I_SO/I_SO
 
71
      INTEGER I_LIB
 
72
      COMMON/I_LIB/I_LIB
 
73
 
 
74
      COMPLEX*16 AMP(NBORNAMPS)
 
75
      COMMON/AMPS/AMP
 
76
      COMPLEX*16 W(20,NWAVEFUNCS)
 
77
      COMMON/W/W
 
78
 
 
79
      COMPLEX*16 WL(MAXLWFSIZE,0:LOOPMAXCOEFS-1,MAXLWFSIZE
 
80
     $ ,0:NLOOPWAVEFUNCS)
 
81
      COMPLEX*16 PL(0:3,0:NLOOPWAVEFUNCS)
 
82
      COMMON/WL/WL,PL
 
83
 
 
84
      COMPLEX*16 AMPL(3,NCTAMPS)
 
85
      COMMON/AMPL/AMPL
 
86
 
 
87
C     
 
88
C     ----------
 
89
C     BEGIN CODE
 
90
C     ----------
 
91
 
 
92
C     The target squared split order contribution is already reached
 
93
C      if true.
 
94
      IF (FILTER_SO.AND.CT_REQ_SO_DONE) THEN
 
95
        GOTO 1001
 
96
      ENDIF
 
97
 
 
98
      CALL IXXXXX(P(0,1),ZERO,NHEL(1),+1*IC(1),W(1,1))
 
99
      CALL OXXXXX(P(0,2),ZERO,NHEL(2),-1*IC(2),W(1,2))
 
100
      CALL VXXXXX(P(0,3),MDL_MW,NHEL(3),+1*IC(3),W(1,3))
 
101
C     Amplitude(s) for born diagram with ID 1
 
102
      CALL FFV2_0(W(1,1),W(1,2),W(1,3),GC_47,AMP(1))
 
103
C     Counter-term amplitude(s) for loop diagram number 2
 
104
      CALL FFV2_0(W(1,1),W(1,2),W(1,3),R2_SXCW,AMPL(1,1))
 
105
C     At this point, all CT amps needed for (QCD=2), i.e. of split
 
106
C      order ID=1, are computed.
 
107
      IF(FILTER_SO.AND.SQSO_TARGET.EQ.1) GOTO 2000
 
108
 
 
109
      GOTO 1001
 
110
 2000 CONTINUE
 
111
      CT_REQ_SO_DONE=.TRUE.
 
112
 1001 CONTINUE
 
113
      END
 
114