← Back to branch summary

~alifson/chiralityflow/trunk

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uxu_ttx%b_sf_001.f

  • Committer: andrew.lifson at lu
  • Date: 2021-09-01 15:34:39 UTC
  • Revision ID: andrew.lifson@thep.lu.se-20210901153439-7fasjhav4cp4m88r
testing a new repository of a madgraph folder

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_uxu_ttx%b_sf_001.f
 
1
      SUBROUTINE SB_SF_001(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~ [ LOonly = QCD QED ] QCD^2<=6 QED^2<=0
 
15
C     Process: c~ c > t t~ [ LOonly = QCD QED ] QCD^2<=6 QED^2<=0
 
16
C     Process: d~ d > t t~ [ LOonly = QCD QED ] QCD^2<=6 QED^2<=0
 
17
C     Process: s~ s > t t~ [ LOonly = QCD QED ] QCD^2<=6 QED^2<=0
 
18
C     spectators: 1 2 
 
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_001_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_001_SPLITORDERS(P1,ANS)
 
89
C     
 
90
C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
 
91
C     By the MadGraph5_aMC@NLO Development Team
 
92
C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
 
93
C     
 
94
C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
 
95
C     AND HELICITIES
 
96
C     FOR THE POINT IN PHASE SPACE P(0:3,NEXTERNAL-1)
 
97
C     
 
98
C     Process: u~ u > t t~ [ LOonly = QCD QED ] QCD^2<=6 QED^2<=0
 
99
C     Process: c~ c > t t~ [ LOonly = QCD QED ] QCD^2<=6 QED^2<=0
 
100
C     Process: d~ d > t t~ [ LOonly = QCD QED ] QCD^2<=6 QED^2<=0
 
101
C     Process: s~ s > t t~ [ LOonly = QCD QED ] QCD^2<=6 QED^2<=0
 
102
C     spectators: 1 2 
 
103
 
 
104
C     
 
105
      IMPLICIT NONE
 
106
C     
 
107
C     CONSTANTS
 
108
C     
 
109
      INCLUDE 'nexternal.inc'
 
110
      INTEGER     NCOMB
 
111
      PARAMETER ( NCOMB=  16 )
 
112
      INTEGER NSQAMPSO
 
113
      PARAMETER (NSQAMPSO=1)
 
114
      INTEGER    THEL
 
115
      PARAMETER (THEL=NCOMB*0)
 
116
      INTEGER NGRAPHS
 
117
      PARAMETER (NGRAPHS=   1)
 
118
C     
 
119
C     ARGUMENTS 
 
120
C     
 
121
      REAL*8 P1(0:3,NEXTERNAL-1),ANS(0:NSQAMPSO)
 
122
C     
 
123
C     LOCAL VARIABLES 
 
124
C     
 
125
      INTEGER IHEL,IDEN,I,J
 
126
      DOUBLE PRECISION T(NSQAMPSO)
 
127
      INTEGER IDEN_VALUES(1)
 
128
      DATA IDEN_VALUES / 36 /
 
129
C     
 
130
C     GLOBAL VARIABLES
 
131
C     
 
132
      LOGICAL GOODHEL(NCOMB,0)
 
133
      COMMON /C_GOODHEL/ GOODHEL
 
134
      DOUBLE PRECISION SAVEMOM(NEXTERNAL-1,2)
 
135
      COMMON/TO_SAVEMOM/SAVEMOM
 
136
      LOGICAL CALCULATEDBORN
 
137
      COMMON/CCALCULATEDBORN/CALCULATEDBORN
 
138
      INTEGER NFKSPROCESS
 
139
      COMMON/C_NFKSPROCESS/NFKSPROCESS
 
140
C     ----------
 
141
C     BEGIN CODE
 
142
C     ----------
 
143
      IDEN=IDEN_VALUES(NFKSPROCESS)
 
144
      IF (CALCULATEDBORN) THEN
 
145
        DO J=1,NEXTERNAL-1
 
146
          IF (SAVEMOM(J,1).NE.P1(0,J) .OR. SAVEMOM(J,2).NE.P1(3,J))
 
147
     $      THEN
 
148
            CALCULATEDBORN=.FALSE.
 
149
            WRITE(*,*) 'Error in sb_sf: momenta not the same in the'
 
150
     $       //' born'
 
151
            STOP
 
152
          ENDIF
 
153
        ENDDO
 
154
      ELSE
 
155
        WRITE(*,*) 'Error in sb_sf: color_linked borns should be'
 
156
     $   //' called only with calculatedborn = true'
 
157
        STOP
 
158
      ENDIF
 
159
      DO I=0,NSQAMPSO
 
160
        ANS(I) = 0D0
 
161
      ENDDO
 
162
      DO IHEL=1,NCOMB
 
163
        IF (GOODHEL(IHEL,NFKSPROCESS)) THEN
 
164
          CALL B_SF_001(IHEL,T)
 
165
          DO I=1,NSQAMPSO
 
166
            ANS(I)=ANS(I)+T(I)
 
167
          ENDDO
 
168
        ENDIF
 
169
      ENDDO
 
170
      DO I=1,NSQAMPSO
 
171
        ANS(I)=ANS(I)/DBLE(IDEN)
 
172
        ANS(0)=ANS(0)+ANS(I)
 
173
      ENDDO
 
174
      END
 
175
 
 
176
 
 
177
      SUBROUTINE B_SF_001(HELL,ANS)
 
178
C     
 
179
C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
 
180
C     By the MadGraph5_aMC@NLO Development Team
 
181
C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
 
182
C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
 
183
C     FOR THE POINT WITH EXTERNAL LINES W(0:6,NEXTERNAL-1)
 
184
 
 
185
C     Process: u~ u > t t~ [ LOonly = QCD QED ] QCD^2<=6 QED^2<=0
 
186
C     Process: c~ c > t t~ [ LOonly = QCD QED ] QCD^2<=6 QED^2<=0
 
187
C     Process: d~ d > t t~ [ LOonly = QCD QED ] QCD^2<=6 QED^2<=0
 
188
C     Process: s~ s > t t~ [ LOonly = QCD QED ] QCD^2<=6 QED^2<=0
 
189
C     spectators: 1 2 
 
190
 
 
191
C     
 
192
      IMPLICIT NONE
 
193
C     
 
194
C     CONSTANTS
 
195
C     
 
196
      INTEGER NAMPSO, NSQAMPSO
 
197
      PARAMETER (NAMPSO=1, NSQAMPSO=1)
 
198
      INTEGER     NGRAPHS
 
199
      PARAMETER ( NGRAPHS = 1 )
 
200
      INTEGER NCOLOR1, NCOLOR2
 
201
      PARAMETER (NCOLOR1=2, NCOLOR2=2)
 
202
      REAL*8     ZERO
 
203
      PARAMETER (ZERO=0D0)
 
204
      COMPLEX*16 IMAG1
 
205
      PARAMETER (IMAG1 = (0D0,1D0))
 
206
      INCLUDE 'nexternal.inc'
 
207
      INCLUDE 'born_nhel.inc'
 
208
C     
 
209
C     ARGUMENTS 
 
210
C     
 
211
      INTEGER HELL
 
212
      REAL*8 ANS(NSQAMPSO)
 
213
C     
 
214
C     LOCAL VARIABLES 
 
215
C     
 
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)
 
221
C     
 
222
C     GLOBAL VARIABLES
 
223
C     
 
224
      DOUBLE COMPLEX SAVEAMP(NGRAPHS,MAX_BHEL)
 
225
      COMMON/TO_SAVEAMP/SAVEAMP
 
226
      LOGICAL CALCULATEDBORN
 
227
      COMMON/CCALCULATEDBORN/CALCULATEDBORN
 
228
C     
 
229
C     FUNCTION
 
230
C     
 
231
      INTEGER SQSOINDEXB
 
232
C     
 
233
C     COLOR DATA
 
234
C     
 
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/
 
239
C     ----------
 
240
C     BEGIN CODE
 
241
C     ----------
 
242
      IF (.NOT. CALCULATEDBORN) THEN
 
243
        WRITE(*,*) 'Error in b_sf: color_linked borns should be called'
 
244
     $   //' only with calculatedborn = true'
 
245
        STOP
 
246
      ELSEIF (CALCULATEDBORN) THEN
 
247
        DO I=1,NGRAPHS
 
248
          AMP(I)=SAVEAMP(I,HELL)
 
249
        ENDDO
 
250
      ENDIF
 
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) = (-2.777777777777778D-02)*AMP(1)
 
256
      JAMP2(2,1) = (8.333333333333333D-02)*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
 
270
        ENDDO
 
271
      ENDDO
 
272
      END
 
273
 
 
274
 
 
275