← Back to branch summary

~maddevelopers/mg5amcnlo/new_clustering

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%b_sf_008.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_pptt_fksreal/%SubProcesses%P0_uux_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: u u~ > t t~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
 
15
C     Process: c c~ > t t~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
 
16
C     Process: d d~ > t t~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
 
17
C     Process: s s~ > t t~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
 
18
C     spectators: 4 4 
 
19
 
 
20
C     
 
21
C     
 
22
C     CONSTANTS
 
23
C     
 
24
      IMPLICIT NONE
 
25
      INCLUDE 'nexternal.inc'
 
26
      INTEGER NSQAMPSO
 
27
      PARAMETER (NSQAMPSO=1)
 
28
C     
 
29
C     ARGUMENTS 
 
30
C     
 
31
      REAL*8 P(0:3,NEXTERNAL), ANS_SUMMED
 
32
C     
 
33
C     VARIABLES
 
34
C     
 
35
      INTEGER I,J
 
36
      INCLUDE 'orders.inc'
 
37
      REAL*8 ANS(0:NSQAMPSO)
 
38
      LOGICAL KEEP_ORDER_CNT(NSPLITORDERS, NSQAMPSO)
 
39
      COMMON /C_KEEP_ORDER_CNT/ KEEP_ORDER_CNT
 
40
      INTEGER AMP_ORDERS(NSPLITORDERS)
 
41
      DOUBLE PRECISION TINY
 
42
      PARAMETER (TINY = 1D-12)
 
43
      DOUBLE PRECISION MAX_VAL
 
44
C     
 
45
C     FUNCTIONS
 
46
C     
 
47
      INTEGER GETORDPOWFROMINDEX_B
 
48
      INTEGER ORDERS_TO_AMP_SPLIT_POS
 
49
C     
 
50
C     BEGIN CODE
 
51
C     
 
52
      CALL SB_SF_008_SPLITORDERS(P,ANS)
 
53
C     color-linked borns are called for QCD-type emissions
 
54
      ANS_SUMMED = 0D0
 
55
      MAX_VAL = 0D0
 
56
 
 
57
C     reset the amp_split_cnt array
 
58
      AMP_SPLIT_CNT(1:AMP_SPLIT_SIZE,1:2,1:NSPLITORDERS) = DCMPLX(0D0
 
59
     $ ,0D0)
 
60
 
 
61
 
 
62
      DO I = 1, NSQAMPSO
 
63
        MAX_VAL = MAX(MAX_VAL, ABS(ANS(I)))
 
64
      ENDDO
 
65
 
 
66
      DO I = 1, NSQAMPSO
 
67
        IF (KEEP_ORDER_CNT(QCD_POS, I)) THEN
 
68
          ANS_SUMMED = ANS_SUMMED + ANS(I)
 
69
          DO J = 1, NSPLITORDERS
 
70
            AMP_ORDERS(J) = GETORDPOWFROMINDEX_B(J, I)
 
71
C           take into account the fact that this is for QCD
 
72
            IF (J.EQ.QCD_POS) AMP_ORDERS(J) = AMP_ORDERS(J) + 2
 
73
          ENDDO
 
74
            !amp_split_cnt(orders_to_amp_split_pos(amp_orders),1,qcd_pos) = ans(I)
 
75
          IF(ABS(ANS(I)).GT.MAX_VAL*TINY)
 
76
     $      AMP_SPLIT_CNT(ORDERS_TO_AMP_SPLIT_POS(AMP_ORDERS),1
 
77
     $     ,QCD_POS) = ANS(I)
 
78
        ENDIF
 
79
      ENDDO
 
80
 
 
81
C     this is to avoid fake non-zero contributions 
 
82
      IF (ABS(ANS_SUMMED).LT.MAX_VAL*TINY) ANS_SUMMED=0D0
 
83
 
 
84
      RETURN
 
85
      END
 
86
 
 
87
 
 
88
      SUBROUTINE SB_SF_008_SPLITORDERS(P1,ANS)
2
89
C     
3
90
C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
4
91
C     By the MadGraph5_aMC@NLO Development Team
8
95
C     AND HELICITIES
9
96
C     FOR THE POINT IN PHASE SPACE P(0:3,NEXTERNAL-1)
10
97
C     
11
 
C     Process: u u~ > t t~ WEIGHTED<=2 [ real = QCD ]
12
 
C     Process: c c~ > t t~ WEIGHTED<=2 [ real = QCD ]
13
 
C     Process: d d~ > t t~ WEIGHTED<=2 [ real = QCD ]
14
 
C     Process: s s~ > t t~ WEIGHTED<=2 [ real = QCD ]
 
98
C     Process: u u~ > t t~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
 
99
C     Process: c c~ > t t~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
 
100
C     Process: d d~ > t t~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
 
101
C     Process: s s~ > t t~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
15
102
C     spectators: 4 4 
16
103
 
17
104
C     
22
109
      INCLUDE 'nexternal.inc'
23
110
      INTEGER     NCOMB
24
111
      PARAMETER ( NCOMB=  16 )
 
112
      INTEGER NSQAMPSO
 
113
      PARAMETER (NSQAMPSO=1)
25
114
      INTEGER    THEL
26
115
      PARAMETER (THEL=NCOMB*6)
27
116
      INTEGER NGRAPHS
29
118
C     
30
119
C     ARGUMENTS 
31
120
C     
32
 
      REAL*8 P1(0:3,NEXTERNAL-1),ANS
 
121
      REAL*8 P1(0:3,NEXTERNAL-1),ANS(0:NSQAMPSO)
33
122
C     
34
123
C     LOCAL VARIABLES 
35
124
C     
36
 
      INTEGER IHEL,IDEN,J
37
 
      REAL*8 B_SF_008
 
125
      INTEGER IHEL,IDEN,I,J
 
126
      DOUBLE PRECISION T(NSQAMPSO)
38
127
      INTEGER IDEN_VALUES(6)
39
128
      DATA IDEN_VALUES /36, 36, 36, 36, 36, 36/
40
129
C     
67
156
     $   //' called only with calculatedborn = true'
68
157
        STOP
69
158
      ENDIF
70
 
      ANS = 0D0
 
159
      DO I=0,NSQAMPSO
 
160
        ANS(I) = 0D0
 
161
      ENDDO
71
162
      DO IHEL=1,NCOMB
72
163
        IF (GOODHEL(IHEL,NFKSPROCESS)) THEN
73
 
          ANS=ANS+B_SF_008(P1,IHEL)
 
164
          CALL B_SF_008(IHEL,T)
 
165
          DO I=1,NSQAMPSO
 
166
            ANS(I)=ANS(I)+T(I)
 
167
          ENDDO
74
168
        ENDIF
75
169
      ENDDO
76
 
      ANS=ANS/DBLE(IDEN)
 
170
      DO I=1,NSQAMPSO
 
171
        ANS(I)=ANS(I)/DBLE(IDEN)
 
172
        ANS(0)=ANS(0)+ANS(I)
 
173
      ENDDO
77
174
      END
78
175
 
79
176
 
80
 
      REAL*8 FUNCTION B_SF_008(P,HELL)
 
177
      SUBROUTINE B_SF_008(HELL,ANS)
81
178
C     
82
179
C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
83
180
C     By the MadGraph5_aMC@NLO Development Team
85
182
C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
86
183
C     FOR THE POINT WITH EXTERNAL LINES W(0:6,NEXTERNAL-1)
87
184
 
88
 
C     Process: u u~ > t t~ WEIGHTED<=2 [ real = QCD ]
89
 
C     Process: c c~ > t t~ WEIGHTED<=2 [ real = QCD ]
90
 
C     Process: d d~ > t t~ WEIGHTED<=2 [ real = QCD ]
91
 
C     Process: s s~ > t t~ WEIGHTED<=2 [ real = QCD ]
 
185
C     Process: u u~ > t t~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
 
186
C     Process: c c~ > t t~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
 
187
C     Process: d d~ > t t~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
 
188
C     Process: s s~ > t t~ [ real = QED QCD ] QCD^2<=6 QED^2<=0
92
189
C     spectators: 4 4 
93
190
 
94
191
C     
96
193
C     
97
194
C     CONSTANTS
98
195
C     
 
196
      INTEGER NAMPSO, NSQAMPSO
 
197
      PARAMETER (NAMPSO=1, NSQAMPSO=1)
99
198
      INTEGER     NGRAPHS
100
199
      PARAMETER ( NGRAPHS = 1 )
101
200
      INTEGER NCOLOR1, NCOLOR2
109
208
C     
110
209
C     ARGUMENTS 
111
210
C     
112
 
      REAL*8 P(0:3,NEXTERNAL-1)
113
211
      INTEGER HELL
 
212
      REAL*8 ANS(NSQAMPSO)
114
213
C     
115
214
C     LOCAL VARIABLES 
116
215
C     
117
 
      INTEGER I,J
118
 
      REAL*8 DENOM(NCOLOR1), CF(NCOLOR2,NCOLOR1)
119
 
      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP1(NCOLOR1), JAMP2(NCOLOR2)
 
216
      INTEGER I,J,M,N
 
217
      REAL*8 CF(NCOLOR2,NCOLOR1)
 
218
      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP1(NCOLOR1,NAMPSO),
 
219
     $  JAMP2(NCOLOR2,NAMPSO)
 
220
      COMPLEX*16 TMP_JAMP(0)
120
221
C     
121
222
C     GLOBAL VARIABLES
122
223
C     
125
226
      LOGICAL CALCULATEDBORN
126
227
      COMMON/CCALCULATEDBORN/CALCULATEDBORN
127
228
C     
 
229
C     FUNCTION
 
230
C     
 
231
      INTEGER SQSOINDEXB
 
232
C     
128
233
C     COLOR DATA
129
234
C     
130
 
      DATA DENOM(1)/1/
131
 
      DATA (CF(I,  1),I=  1,  2) /    9,    3/
132
 
      DATA DENOM(2)/1/
133
 
      DATA (CF(I,  2),I=  1,  2) /    3,    9/
 
235
      DATA (CF(I,  1),I=  1,  2) /9.000000000000000D+00
 
236
     $ ,3.000000000000000D+00/
 
237
      DATA (CF(I,  2),I=  1,  2) /3.000000000000000D+00
 
238
     $ ,9.000000000000000D+00/
134
239
C     ----------
135
240
C     BEGIN CODE
136
241
C     ----------
143
248
          AMP(I)=SAVEAMP(I,HELL)
144
249
        ENDDO
145
250
      ENDIF
146
 
      JAMP1(1)=+1D0/2D0*(+1D0/3D0*AMP(1))
147
 
      JAMP1(2)=+1D0/2D0*(-AMP(1))
148
 
      JAMP2(1)=+1D0/8D0*(+AMP(1)-1D0/9D0*AMP(1))
149
 
      JAMP2(2)=+1D0/8D0*(-3D0*AMP(1)+1D0/3D0*AMP(1))
150
 
      B_SF_008 = 0.D0
151
 
      DO I = 1, NCOLOR1
152
 
        ZTEMP = (0.D0,0.D0)
153
 
        DO J = 1, NCOLOR2
154
 
          ZTEMP = ZTEMP + CF(J,I)*JAMP2(J)
 
251
C     JAMPs contributing to orders QCD=2 QED=0
 
252
      JAMP1(1,1) = (1.666666666666667D-01)*AMP(1)
 
253
      JAMP1(2,1) = (-5.000000000000000D-01)*AMP(1)
 
254
C     JAMPs contributing to orders QCD=2 QED=0
 
255
      JAMP2(1,1) = (1.111111111111111D-01)*AMP(1)
 
256
      JAMP2(2,1) = (-3.333333333333333D-01)*AMP(1)
 
257
      DO I = 1, NSQAMPSO
 
258
        ANS(I) = 0D0
 
259
      ENDDO
 
260
      DO M = 1, NAMPSO
 
261
        DO I = 1, NCOLOR1
 
262
          ZTEMP = (0.D0,0.D0)
 
263
          DO J = 1, NCOLOR2
 
264
            ZTEMP = ZTEMP + CF(J,I)*JAMP2(J,M)
 
265
          ENDDO
 
266
          DO N = 1, NAMPSO
 
267
            ANS(SQSOINDEXB(M,N))=ANS(SQSOINDEXB(M,N))+ZTEMP
 
268
     $       *DCONJG(JAMP1(I,N))
 
269
          ENDDO
155
270
        ENDDO
156
 
        B_SF_008 =B_SF_008+ZTEMP*DCONJG(JAMP1(I))/DENOM(I)
157
271
      ENDDO
158
272
      END
159
273