← Back to branch summary

~maddevelopers/mg5amcnlo/FKS_EW_granny-run_py8_at_evtgen

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_ttx%b_sf_008.f

  • Committer: Marco Zaro
  • Date: 2017-03-08 10:31:32 UTC
  • Revision ID: marco.zaro@gmail.com-20170308103132-vsnzjovp9f90gc4y
fixing tests

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_ttx%b_sf_008.f
1
 
      SUBROUTINE SB_SF_008(P1,ANS)
 
1
      SUBROUTINE SB_SF_008(P,ANS_SUMMED)
 
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     
 
8
C     Return the sum of the split orders which are required in
 
9
C      orders.inc (BORN_ORDERS)
 
10
C     Also the values needed for the counterterms are stored in the
 
11
C      C_BORN_CNT common block
 
12
C     
 
13
C     
 
14
C     Process: g g > t t~ [ LOonly = QED QCD ] QCD^2=6 QED^2=0
 
15
C     spectators: 4 4 
 
16
 
 
17
C     
 
18
C     
 
19
C     CONSTANTS
 
20
C     
 
21
      IMPLICIT NONE
 
22
      INCLUDE 'nexternal.inc'
 
23
      INTEGER NSQAMPSO
 
24
      PARAMETER (NSQAMPSO=1)
 
25
C     
 
26
C     ARGUMENTS 
 
27
C     
 
28
      REAL*8 P(0:3,NEXTERNAL), ANS_SUMMED
 
29
C     
 
30
C     VARIABLES
 
31
C     
 
32
      INTEGER I,J
 
33
      INCLUDE 'orders.inc'
 
34
      REAL*8 ANS(0:NSQAMPSO)
 
35
      LOGICAL KEEP_ORDER_CNT(NSPLITORDERS, NSQAMPSO)
 
36
      COMMON /C_KEEP_ORDER_CNT/ KEEP_ORDER_CNT
 
37
      INTEGER AMP_ORDERS(NSPLITORDERS)
 
38
C     
 
39
C     FUNCTIONS
 
40
C     
 
41
      INTEGER GETORDPOWFROMINDEX_B
 
42
      INTEGER ORDERS_TO_AMP_SPLIT_POS
 
43
C     
 
44
C     BEGIN CODE
 
45
C     
 
46
      CALL SB_SF_008_SPLITORDERS(P,ANS)
 
47
C     color-linked borns are called for QCD-type emissions
 
48
      ANS_SUMMED = 0D0
 
49
 
 
50
C     reset the amp_split_cnt array
 
51
      DO I = 1, AMP_SPLIT_SIZE
 
52
        DO J = 1, NSPLITORDERS
 
53
          AMP_SPLIT_CNT(I,1,J) = 0D0
 
54
          AMP_SPLIT_CNT(I,2,J) = 0D0
 
55
        ENDDO
 
56
      ENDDO
 
57
 
 
58
      DO I = 1, NSQAMPSO
 
59
        IF (KEEP_ORDER_CNT(QCD_POS, I)) THEN
 
60
          ANS_SUMMED = ANS_SUMMED + ANS(I)
 
61
          DO J = 1, NSPLITORDERS
 
62
            AMP_ORDERS(J) = GETORDPOWFROMINDEX_B(J, I)
 
63
C           take into account the fact that this is for QCD
 
64
            IF (J.EQ.QCD_POS) AMP_ORDERS(J) = AMP_ORDERS(J) + 2
 
65
          ENDDO
 
66
          AMP_SPLIT_CNT(ORDERS_TO_AMP_SPLIT_POS(AMP_ORDERS),1,QCD_POS)
 
67
     $      = ANS(I)
 
68
        ENDIF
 
69
      ENDDO
 
70
 
 
71
      RETURN
 
72
      END
 
73
 
 
74
 
 
75
      SUBROUTINE SB_SF_008_SPLITORDERS(P1,ANS)
2
76
C     
3
77
C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
4
78
C     By the MadGraph5_aMC@NLO Development Team
8
82
C     AND HELICITIES
9
83
C     FOR THE POINT IN PHASE SPACE P(0:3,NEXTERNAL-1)
10
84
C     
11
 
C     Process: g g > t t~ WEIGHTED<=2 [ LOonly = QCD ]
 
85
C     Process: g g > t t~ [ LOonly = QED QCD ] QCD^2=6 QED^2=0
12
86
C     spectators: 4 4 
13
87
 
14
88
C     
19
93
      INCLUDE 'nexternal.inc'
20
94
      INTEGER     NCOMB
21
95
      PARAMETER ( NCOMB=  16 )
 
96
      INTEGER NSQAMPSO
 
97
      PARAMETER (NSQAMPSO=1)
22
98
      INTEGER    THEL
23
99
      PARAMETER (THEL=NCOMB*0)
24
100
      INTEGER NGRAPHS
26
102
C     
27
103
C     ARGUMENTS 
28
104
C     
29
 
      REAL*8 P1(0:3,NEXTERNAL-1),ANS
 
105
      REAL*8 P1(0:3,NEXTERNAL-1),ANS(0:NSQAMPSO)
30
106
C     
31
107
C     LOCAL VARIABLES 
32
108
C     
33
 
      INTEGER IHEL,IDEN,J
34
 
      REAL*8 B_SF_008
 
109
      INTEGER IHEL,IDEN,I,J
 
110
      DOUBLE PRECISION T(NSQAMPSO)
35
111
      INTEGER IDEN_VALUES(1)
36
112
      DATA IDEN_VALUES / 256 /
37
113
C     
64
140
     $   //' called only with calculatedborn = true'
65
141
        STOP
66
142
      ENDIF
67
 
      ANS = 0D0
 
143
      DO I=0,NSQAMPSO
 
144
        ANS(I) = 0D0
 
145
      ENDDO
68
146
      DO IHEL=1,NCOMB
69
147
        IF (GOODHEL(IHEL,NFKSPROCESS)) THEN
70
 
          ANS=ANS+B_SF_008(P1,IHEL)
 
148
          CALL B_SF_008(IHEL,T)
 
149
          DO I=1,NSQAMPSO
 
150
            ANS(I)=ANS(I)+T(I)
 
151
          ENDDO
71
152
        ENDIF
72
153
      ENDDO
73
 
      ANS=ANS/DBLE(IDEN)
 
154
      DO I=1,NSQAMPSO
 
155
        ANS(I)=ANS(I)/DBLE(IDEN)
 
156
        ANS(0)=ANS(0)+ANS(I)
 
157
      ENDDO
74
158
      END
75
159
 
76
160
 
77
 
      REAL*8 FUNCTION B_SF_008(P,HELL)
 
161
      SUBROUTINE B_SF_008(HELL,ANS)
78
162
C     
79
163
C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
80
164
C     By the MadGraph5_aMC@NLO Development Team
82
166
C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
83
167
C     FOR THE POINT WITH EXTERNAL LINES W(0:6,NEXTERNAL-1)
84
168
 
85
 
C     Process: g g > t t~ WEIGHTED<=2 [ LOonly = QCD ]
 
169
C     Process: g g > t t~ [ LOonly = QED QCD ] QCD^2=6 QED^2=0
86
170
C     spectators: 4 4 
87
171
 
88
172
C     
90
174
C     
91
175
C     CONSTANTS
92
176
C     
 
177
      INTEGER NAMPSO, NSQAMPSO
 
178
      PARAMETER (NAMPSO=1, NSQAMPSO=1)
93
179
      INTEGER     NGRAPHS
94
180
      PARAMETER ( NGRAPHS = 3 )
95
181
      INTEGER NCOLOR1, NCOLOR2
103
189
C     
104
190
C     ARGUMENTS 
105
191
C     
106
 
      REAL*8 P(0:3,NEXTERNAL-1)
107
192
      INTEGER HELL
 
193
      REAL*8 ANS(NSQAMPSO)
108
194
C     
109
195
C     LOCAL VARIABLES 
110
196
C     
111
 
      INTEGER I,J
 
197
      INTEGER I,J,M,N
112
198
      REAL*8 DENOM(NCOLOR1), CF(NCOLOR2,NCOLOR1)
113
 
      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP1(NCOLOR1), JAMP2(NCOLOR2)
 
199
      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP1(NCOLOR1,NAMPSO),
 
200
     $  JAMP2(NCOLOR2,NAMPSO)
114
201
C     
115
202
C     GLOBAL VARIABLES
116
203
C     
119
206
      LOGICAL CALCULATEDBORN
120
207
      COMMON/CCALCULATEDBORN/CALCULATEDBORN
121
208
C     
 
209
C     FUNCTION
 
210
C     
 
211
      INTEGER SQSOINDEXB
 
212
C     
122
213
C     COLOR DATA
123
214
C     
124
215
      DATA DENOM(1)/3/
137
228
          AMP(I)=SAVEAMP(I,HELL)
138
229
        ENDDO
139
230
      ENDIF
140
 
      JAMP1(1)=+IMAG1*AMP(1)-AMP(2)
141
 
      JAMP1(2)=-IMAG1*AMP(1)-AMP(3)
142
 
      JAMP2(1)=+1D0/4D0*(+3D0*IMAG1*AMP(1)-1D0/3D0*IMAG1*AMP(1)-3D0
 
231
C     JAMP1s contributing to orders QCD=2 QED=0
 
232
      JAMP1(1,1)=+IMAG1*AMP(1)-AMP(2)
 
233
      JAMP1(2,1)=-IMAG1*AMP(1)-AMP(3)
 
234
C     JAMP2s contributing to orders QCD=2 QED=0
 
235
      JAMP2(1,1)=+1D0/4D0*(+3D0*IMAG1*AMP(1)-1D0/3D0*IMAG1*AMP(1)-3D0
143
236
     $ *AMP(2)+1D0/3D0*AMP(2))
144
 
      JAMP2(2)=+1D0/4D0*(-3D0*IMAG1*AMP(1)+1D0/3D0*IMAG1*AMP(1)-3D0
 
237
      JAMP2(2,1)=+1D0/4D0*(-3D0*IMAG1*AMP(1)+1D0/3D0*IMAG1*AMP(1)-3D0
145
238
     $ *AMP(3)+1D0/3D0*AMP(3))
146
 
      B_SF_008 = 0.D0
147
 
      DO I = 1, NCOLOR1
148
 
        ZTEMP = (0.D0,0.D0)
149
 
        DO J = 1, NCOLOR2
150
 
          ZTEMP = ZTEMP + CF(J,I)*JAMP2(J)
 
239
      DO I = 1, NSQAMPSO
 
240
        ANS(I) = 0D0
 
241
      ENDDO
 
242
      DO M = 1, NAMPSO
 
243
        DO I = 1, NCOLOR1
 
244
          ZTEMP = (0.D0,0.D0)
 
245
          DO J = 1, NCOLOR2
 
246
            ZTEMP = ZTEMP + CF(J,I)*JAMP2(J,M)
 
247
          ENDDO
 
248
          DO N = 1, NAMPSO
 
249
            ANS(SQSOINDEXB(M,N))=ANS(SQSOINDEXB(M,N))+ZTEMP
 
250
     $       *DCONJG(JAMP1(I,N))/DENOM(I)
 
251
          ENDDO
151
252
        ENDDO
152
 
        B_SF_008 =B_SF_008+ZTEMP*DCONJG(JAMP1(I))/DENOM(I)
153
253
      ENDDO
154
254
      END
155
255