← Back to branch summary

~maddevelopers/mg5amcnlo/2.0.2_onshellsubtract_newresh

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%matrix_2.f

  • Committer: Marco Zaro
  • Date: 2014-03-31 09:10:24 UTC
  • mfrom: (285.1.4 2.0.2_onshellsubtract)
  • Revision ID: marco.zaro@gmail.com-20140331091024-gkn2xz7jzq6t5bap
merged with latest 2.0.2_onshellsubtract

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%matrix_2.f
1
 
      SUBROUTINE SMATRIX_2(P,ANS)
2
 
C     
3
 
C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
4
 
C     By the MadGraph5_aMC@NLO Development Team
5
 
C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
6
 
C     
7
 
C     Returns amplitude squared summed/avg over colors
8
 
C     and helicities
9
 
C     for the point in phase space P(0:3,NEXTERNAL)
10
 
C     
11
 
C     Process: g u~ > t t~ u~ WEIGHTED=3 [ QCD ]
12
 
C     Process: g c~ > t t~ c~ WEIGHTED=3 [ QCD ]
13
 
C     Process: g d~ > t t~ d~ WEIGHTED=3 [ QCD ]
14
 
C     Process: g s~ > t t~ s~ WEIGHTED=3 [ QCD ]
15
 
C     
16
 
      IMPLICIT NONE
17
 
C     
18
 
C     CONSTANTS
19
 
C     
20
 
      INCLUDE 'nexternal.inc'
21
 
      INTEGER     NCOMB
22
 
      PARAMETER ( NCOMB=32)
23
 
C     
24
 
C     ARGUMENTS 
25
 
C     
26
 
      REAL*8 P(0:3,NEXTERNAL),ANS
27
 
C     
28
 
C     LOCAL VARIABLES 
29
 
C     
30
 
      INTEGER IHEL,IDEN,I,T_IDENT(NCOMB)
31
 
      REAL*8 MATRIX_2
32
 
      REAL*8 T,T_SAVE(NCOMB)
33
 
      SAVE T_SAVE,T_IDENT
34
 
      INTEGER NHEL(NEXTERNAL,NCOMB)
35
 
      DATA (NHEL(I,   1),I=1,5) /-1,-1,-1,-1,-1/
36
 
      DATA (NHEL(I,   2),I=1,5) /-1,-1,-1,-1, 1/
37
 
      DATA (NHEL(I,   3),I=1,5) /-1,-1,-1, 1,-1/
38
 
      DATA (NHEL(I,   4),I=1,5) /-1,-1,-1, 1, 1/
39
 
      DATA (NHEL(I,   5),I=1,5) /-1,-1, 1,-1,-1/
40
 
      DATA (NHEL(I,   6),I=1,5) /-1,-1, 1,-1, 1/
41
 
      DATA (NHEL(I,   7),I=1,5) /-1,-1, 1, 1,-1/
42
 
      DATA (NHEL(I,   8),I=1,5) /-1,-1, 1, 1, 1/
43
 
      DATA (NHEL(I,   9),I=1,5) /-1, 1,-1,-1,-1/
44
 
      DATA (NHEL(I,  10),I=1,5) /-1, 1,-1,-1, 1/
45
 
      DATA (NHEL(I,  11),I=1,5) /-1, 1,-1, 1,-1/
46
 
      DATA (NHEL(I,  12),I=1,5) /-1, 1,-1, 1, 1/
47
 
      DATA (NHEL(I,  13),I=1,5) /-1, 1, 1,-1,-1/
48
 
      DATA (NHEL(I,  14),I=1,5) /-1, 1, 1,-1, 1/
49
 
      DATA (NHEL(I,  15),I=1,5) /-1, 1, 1, 1,-1/
50
 
      DATA (NHEL(I,  16),I=1,5) /-1, 1, 1, 1, 1/
51
 
      DATA (NHEL(I,  17),I=1,5) / 1,-1,-1,-1,-1/
52
 
      DATA (NHEL(I,  18),I=1,5) / 1,-1,-1,-1, 1/
53
 
      DATA (NHEL(I,  19),I=1,5) / 1,-1,-1, 1,-1/
54
 
      DATA (NHEL(I,  20),I=1,5) / 1,-1,-1, 1, 1/
55
 
      DATA (NHEL(I,  21),I=1,5) / 1,-1, 1,-1,-1/
56
 
      DATA (NHEL(I,  22),I=1,5) / 1,-1, 1,-1, 1/
57
 
      DATA (NHEL(I,  23),I=1,5) / 1,-1, 1, 1,-1/
58
 
      DATA (NHEL(I,  24),I=1,5) / 1,-1, 1, 1, 1/
59
 
      DATA (NHEL(I,  25),I=1,5) / 1, 1,-1,-1,-1/
60
 
      DATA (NHEL(I,  26),I=1,5) / 1, 1,-1,-1, 1/
61
 
      DATA (NHEL(I,  27),I=1,5) / 1, 1,-1, 1,-1/
62
 
      DATA (NHEL(I,  28),I=1,5) / 1, 1,-1, 1, 1/
63
 
      DATA (NHEL(I,  29),I=1,5) / 1, 1, 1,-1,-1/
64
 
      DATA (NHEL(I,  30),I=1,5) / 1, 1, 1,-1, 1/
65
 
      DATA (NHEL(I,  31),I=1,5) / 1, 1, 1, 1,-1/
66
 
      DATA (NHEL(I,  32),I=1,5) / 1, 1, 1, 1, 1/
67
 
      LOGICAL GOODHEL(NCOMB)
68
 
      DATA GOODHEL/NCOMB*.FALSE./
69
 
      INTEGER NTRY
70
 
      DATA NTRY/0/
71
 
      DATA IDEN/96/
72
 
C     ----------
73
 
C     BEGIN CODE
74
 
C     ----------
75
 
      NTRY=NTRY+1
76
 
      ANS = 0D0
77
 
      DO IHEL=1,NCOMB
78
 
        IF (GOODHEL(IHEL) .OR. NTRY .LT. 2) THEN
79
 
          IF (NTRY.LT.2) THEN
80
 
C           for the first ps-point, check for helicities that give
81
 
C           identical matrix elements
82
 
            T=MATRIX_2(P ,NHEL(1,IHEL))
83
 
            T_SAVE(IHEL)=T
84
 
            T_IDENT(IHEL)=-1
85
 
            DO I=1,IHEL-1
86
 
              IF (T.EQ.0D0) EXIT
87
 
              IF (T_SAVE(I).EQ.0D0) CYCLE
88
 
              IF (ABS(T/T_SAVE(I)-1D0) .LT. 1D-12) THEN
89
 
C               WRITE (*,*) 'FOUND IDENTICAL',T,IHEL,T_SAVE(I),I
90
 
                T_IDENT(IHEL) = I
91
 
              ENDIF
92
 
            ENDDO
93
 
          ELSE
94
 
            IF (T_IDENT(IHEL).GT.0) THEN
95
 
C             if two helicity states are identical, dont recompute
96
 
              T=T_SAVE(T_IDENT(IHEL))
97
 
              T_SAVE(IHEL)=T
98
 
            ELSE
99
 
              T=MATRIX_2(P ,NHEL(1,IHEL))
100
 
              T_SAVE(IHEL)=T
101
 
            ENDIF
102
 
          ENDIF
103
 
C         add to the sum of helicities
104
 
          ANS=ANS+T
105
 
          IF (T .NE. 0D0 .AND. .NOT. GOODHEL(IHEL)) THEN
106
 
            GOODHEL(IHEL)=.TRUE.
107
 
          ENDIF
108
 
        ENDIF
109
 
      ENDDO
110
 
      ANS=ANS/DBLE(IDEN)
111
 
      END
112
 
 
113
 
 
114
 
      REAL*8 FUNCTION MATRIX_2(P,NHEL)
115
 
C     
116
 
C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
117
 
C     By the MadGraph5_aMC@NLO Development Team
118
 
C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
119
 
C     
120
 
C     Returns amplitude squared summed/avg over colors
121
 
C     for the point with external lines W(0:6,NEXTERNAL)
122
 
C     
123
 
C     Process: g u~ > t t~ u~ WEIGHTED=3 [ QCD ]
124
 
C     Process: g c~ > t t~ c~ WEIGHTED=3 [ QCD ]
125
 
C     Process: g d~ > t t~ d~ WEIGHTED=3 [ QCD ]
126
 
C     Process: g s~ > t t~ s~ WEIGHTED=3 [ QCD ]
127
 
C     
128
 
      IMPLICIT NONE
129
 
C     
130
 
C     CONSTANTS
131
 
C     
132
 
      INTEGER    NGRAPHS
133
 
      PARAMETER (NGRAPHS=5)
134
 
      INTEGER    NWAVEFUNCS, NCOLOR
135
 
      PARAMETER (NWAVEFUNCS=8, NCOLOR=4)
136
 
      REAL*8     ZERO
137
 
      PARAMETER (ZERO=0D0)
138
 
      COMPLEX*16 IMAG1
139
 
      PARAMETER (IMAG1=(0D0,1D0))
140
 
      INCLUDE 'nexternal.inc'
141
 
      INCLUDE 'coupl.inc'
142
 
C     
143
 
C     ARGUMENTS 
144
 
C     
145
 
      REAL*8 P(0:3,NEXTERNAL)
146
 
      INTEGER NHEL(NEXTERNAL)
147
 
C     
148
 
C     LOCAL VARIABLES 
149
 
C     
150
 
      INTEGER I,J
151
 
      INTEGER IC(NEXTERNAL)
152
 
      DATA IC /NEXTERNAL*1/
153
 
      REAL*8 DENOM(NCOLOR), CF(NCOLOR,NCOLOR)
154
 
      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP(NCOLOR), W(8,NWAVEFUNCS)
155
 
C     
156
 
C     COLOR DATA
157
 
C     
158
 
      DATA DENOM(1)/1/
159
 
      DATA (CF(I,  1),I=  1,  4) /   12,    4,    4,    0/
160
 
C     1 T(1,2,4) T(3,5)
161
 
      DATA DENOM(2)/1/
162
 
      DATA (CF(I,  2),I=  1,  4) /    4,   12,    0,    4/
163
 
C     1 T(1,2,5) T(3,4)
164
 
      DATA DENOM(3)/1/
165
 
      DATA (CF(I,  3),I=  1,  4) /    4,    0,   12,    4/
166
 
C     1 T(1,3,4) T(2,5)
167
 
      DATA DENOM(4)/1/
168
 
      DATA (CF(I,  4),I=  1,  4) /    0,    4,    4,   12/
169
 
C     1 T(1,3,5) T(2,4)
170
 
C     ----------
171
 
C     BEGIN CODE
172
 
C     ----------
173
 
      CALL VXXXXX(P(0,1),ZERO,NHEL(1),-1*IC(1),W(1,1))
174
 
      CALL OXXXXX(P(0,2),ZERO,NHEL(2),-1*IC(2),W(1,2))
175
 
      CALL OXXXXX(P(0,3),MDL_MT,NHEL(3),+1*IC(3),W(1,3))
176
 
      CALL IXXXXX(P(0,4),MDL_MT,NHEL(4),-1*IC(4),W(1,4))
177
 
      CALL IXXXXX(P(0,5),ZERO,NHEL(5),-1*IC(5),W(1,5))
178
 
      CALL FFV1_1(W(1,2),W(1,1),GC_11,ZERO,ZERO,W(1,6))
179
 
      CALL FFV1P0_3(W(1,4),W(1,3),GC_11,ZERO,ZERO,W(1,7))
180
 
C     Amplitude(s) for diagram number 1
181
 
      CALL FFV1_0(W(1,5),W(1,6),W(1,7),GC_11,AMP(1))
182
 
      CALL FFV1_1(W(1,3),W(1,1),GC_11,MDL_MT,MDL_WT,W(1,6))
183
 
      CALL FFV1P0_3(W(1,5),W(1,2),GC_11,ZERO,ZERO,W(1,8))
184
 
C     Amplitude(s) for diagram number 2
185
 
      CALL FFV1_0(W(1,4),W(1,6),W(1,8),GC_11,AMP(2))
186
 
      CALL FFV1_2(W(1,4),W(1,1),GC_11,MDL_MT,MDL_WT,W(1,6))
187
 
C     Amplitude(s) for diagram number 3
188
 
      CALL FFV1_0(W(1,6),W(1,3),W(1,8),GC_11,AMP(3))
189
 
      CALL FFV1_2(W(1,5),W(1,1),GC_11,ZERO,ZERO,W(1,6))
190
 
C     Amplitude(s) for diagram number 4
191
 
      CALL FFV1_0(W(1,6),W(1,2),W(1,7),GC_11,AMP(4))
192
 
C     Amplitude(s) for diagram number 5
193
 
      CALL VVV1_0(W(1,1),W(1,8),W(1,7),GC_10,AMP(5))
194
 
      JAMP(1)=+1D0/2D0*(+AMP(1)+AMP(3)+IMAG1*AMP(5))
195
 
      JAMP(2)=+1D0/2D0*(-1D0/3D0*AMP(1)-1D0/3D0*AMP(4))
196
 
      JAMP(3)=+1D0/2D0*(-1D0/3D0*AMP(2)-1D0/3D0*AMP(3))
197
 
      JAMP(4)=+1D0/2D0*(+AMP(2)+AMP(4)-IMAG1*AMP(5))
198
 
      MATRIX_2 = 0.D0
199
 
      DO I = 1, NCOLOR
200
 
        ZTEMP = (0.D0,0.D0)
201
 
        DO J = 1, NCOLOR
202
 
          ZTEMP = ZTEMP + CF(J,I)*JAMP(J)
203
 
        ENDDO
204
 
        MATRIX_2 = MATRIX_2+ZTEMP*DCONJG(JAMP(I))/DENOM(I)
205
 
      ENDDO
206
 
      END
207