← 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_uux_ttx%b_sf_003.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_uux_ttx%b_sf_003.f
1
 
      SUBROUTINE SB_SF_003(P1,ANS)
 
1
      SUBROUTINE SB_SF_003(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 = QED QCD ] QCD^2=6 QED^2=0
 
15
C     Process: c c~ > t t~ [ LOonly = QED QCD ] QCD^2=6 QED^2=0
 
16
C     Process: d d~ > t t~ [ LOonly = QED QCD ] QCD^2=6 QED^2=0
 
17
C     Process: s s~ > t t~ [ LOonly = QED QCD ] QCD^2=6 QED^2=0
 
18
C     spectators: 1 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
C     
 
42
C     FUNCTIONS
 
43
C     
 
44
      INTEGER GETORDPOWFROMINDEX_B
 
45
      INTEGER ORDERS_TO_AMP_SPLIT_POS
 
46
C     
 
47
C     BEGIN CODE
 
48
C     
 
49
      CALL SB_SF_003_SPLITORDERS(P,ANS)
 
50
C     color-linked borns are called for QCD-type emissions
 
51
      ANS_SUMMED = 0D0
 
52
 
 
53
C     reset the amp_split_cnt array
 
54
      DO I = 1, AMP_SPLIT_SIZE
 
55
        DO J = 1, NSPLITORDERS
 
56
          AMP_SPLIT_CNT(I,1,J) = 0D0
 
57
          AMP_SPLIT_CNT(I,2,J) = 0D0
 
58
        ENDDO
 
59
      ENDDO
 
60
 
 
61
      DO I = 1, NSQAMPSO
 
62
        IF (KEEP_ORDER_CNT(QCD_POS, I)) THEN
 
63
          ANS_SUMMED = ANS_SUMMED + ANS(I)
 
64
          DO J = 1, NSPLITORDERS
 
65
            AMP_ORDERS(J) = GETORDPOWFROMINDEX_B(J, I)
 
66
C           take into account the fact that this is for QCD
 
67
            IF (J.EQ.QCD_POS) AMP_ORDERS(J) = AMP_ORDERS(J) + 2
 
68
          ENDDO
 
69
          AMP_SPLIT_CNT(ORDERS_TO_AMP_SPLIT_POS(AMP_ORDERS),1,QCD_POS)
 
70
     $      = ANS(I)
 
71
        ENDIF
 
72
      ENDDO
 
73
 
 
74
      RETURN
 
75
      END
 
76
 
 
77
 
 
78
      SUBROUTINE SB_SF_003_SPLITORDERS(P1,ANS)
2
79
C     
3
80
C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
4
81
C     By the MadGraph5_aMC@NLO Development Team
8
85
C     AND HELICITIES
9
86
C     FOR THE POINT IN PHASE SPACE P(0:3,NEXTERNAL-1)
10
87
C     
11
 
C     Process: u u~ > t t~ WEIGHTED<=2 [ LOonly = QCD ]
12
 
C     Process: c c~ > t t~ WEIGHTED<=2 [ LOonly = QCD ]
13
 
C     Process: d d~ > t t~ WEIGHTED<=2 [ LOonly = QCD ]
14
 
C     Process: s s~ > t t~ WEIGHTED<=2 [ LOonly = QCD ]
 
88
C     Process: u u~ > t t~ [ LOonly = QED QCD ] QCD^2=6 QED^2=0
 
89
C     Process: c c~ > t t~ [ LOonly = QED QCD ] QCD^2=6 QED^2=0
 
90
C     Process: d d~ > t t~ [ LOonly = QED QCD ] QCD^2=6 QED^2=0
 
91
C     Process: s s~ > t t~ [ LOonly = QED QCD ] QCD^2=6 QED^2=0
15
92
C     spectators: 1 4 
16
93
 
17
94
C     
22
99
      INCLUDE 'nexternal.inc'
23
100
      INTEGER     NCOMB
24
101
      PARAMETER ( NCOMB=  16 )
 
102
      INTEGER NSQAMPSO
 
103
      PARAMETER (NSQAMPSO=1)
25
104
      INTEGER    THEL
26
105
      PARAMETER (THEL=NCOMB*0)
27
106
      INTEGER NGRAPHS
29
108
C     
30
109
C     ARGUMENTS 
31
110
C     
32
 
      REAL*8 P1(0:3,NEXTERNAL-1),ANS
 
111
      REAL*8 P1(0:3,NEXTERNAL-1),ANS(0:NSQAMPSO)
33
112
C     
34
113
C     LOCAL VARIABLES 
35
114
C     
36
 
      INTEGER IHEL,IDEN,J
37
 
      REAL*8 B_SF_003
 
115
      INTEGER IHEL,IDEN,I,J
 
116
      DOUBLE PRECISION T(NSQAMPSO)
38
117
      INTEGER IDEN_VALUES(1)
39
118
      DATA IDEN_VALUES / 36 /
40
119
C     
67
146
     $   //' called only with calculatedborn = true'
68
147
        STOP
69
148
      ENDIF
70
 
      ANS = 0D0
 
149
      DO I=0,NSQAMPSO
 
150
        ANS(I) = 0D0
 
151
      ENDDO
71
152
      DO IHEL=1,NCOMB
72
153
        IF (GOODHEL(IHEL,NFKSPROCESS)) THEN
73
 
          ANS=ANS+B_SF_003(P1,IHEL)
 
154
          CALL B_SF_003(IHEL,T)
 
155
          DO I=1,NSQAMPSO
 
156
            ANS(I)=ANS(I)+T(I)
 
157
          ENDDO
74
158
        ENDIF
75
159
      ENDDO
76
 
      ANS=ANS/DBLE(IDEN)
 
160
      DO I=1,NSQAMPSO
 
161
        ANS(I)=ANS(I)/DBLE(IDEN)
 
162
        ANS(0)=ANS(0)+ANS(I)
 
163
      ENDDO
77
164
      END
78
165
 
79
166
 
80
 
      REAL*8 FUNCTION B_SF_003(P,HELL)
 
167
      SUBROUTINE B_SF_003(HELL,ANS)
81
168
C     
82
169
C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
83
170
C     By the MadGraph5_aMC@NLO Development Team
85
172
C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
86
173
C     FOR THE POINT WITH EXTERNAL LINES W(0:6,NEXTERNAL-1)
87
174
 
88
 
C     Process: u u~ > t t~ WEIGHTED<=2 [ LOonly = QCD ]
89
 
C     Process: c c~ > t t~ WEIGHTED<=2 [ LOonly = QCD ]
90
 
C     Process: d d~ > t t~ WEIGHTED<=2 [ LOonly = QCD ]
91
 
C     Process: s s~ > t t~ WEIGHTED<=2 [ LOonly = QCD ]
 
175
C     Process: u u~ > t t~ [ LOonly = QED QCD ] QCD^2=6 QED^2=0
 
176
C     Process: c c~ > t t~ [ LOonly = QED QCD ] QCD^2=6 QED^2=0
 
177
C     Process: d d~ > t t~ [ LOonly = QED QCD ] QCD^2=6 QED^2=0
 
178
C     Process: s s~ > t t~ [ LOonly = QED QCD ] QCD^2=6 QED^2=0
92
179
C     spectators: 1 4 
93
180
 
94
181
C     
96
183
C     
97
184
C     CONSTANTS
98
185
C     
 
186
      INTEGER NAMPSO, NSQAMPSO
 
187
      PARAMETER (NAMPSO=1, NSQAMPSO=1)
99
188
      INTEGER     NGRAPHS
100
189
      PARAMETER ( NGRAPHS = 1 )
101
190
      INTEGER NCOLOR1, NCOLOR2
109
198
C     
110
199
C     ARGUMENTS 
111
200
C     
112
 
      REAL*8 P(0:3,NEXTERNAL-1)
113
201
      INTEGER HELL
 
202
      REAL*8 ANS(NSQAMPSO)
114
203
C     
115
204
C     LOCAL VARIABLES 
116
205
C     
117
 
      INTEGER I,J
 
206
      INTEGER I,J,M,N
118
207
      REAL*8 DENOM(NCOLOR1), CF(NCOLOR2,NCOLOR1)
119
 
      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP1(NCOLOR1), JAMP2(NCOLOR2)
 
208
      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP1(NCOLOR1,NAMPSO),
 
209
     $  JAMP2(NCOLOR2,NAMPSO)
120
210
C     
121
211
C     GLOBAL VARIABLES
122
212
C     
125
215
      LOGICAL CALCULATEDBORN
126
216
      COMMON/CCALCULATEDBORN/CALCULATEDBORN
127
217
C     
 
218
C     FUNCTION
 
219
C     
 
220
      INTEGER SQSOINDEXB
 
221
C     
128
222
C     COLOR DATA
129
223
C     
130
224
      DATA DENOM(1)/1/
143
237
          AMP(I)=SAVEAMP(I,HELL)
144
238
        ENDDO
145
239
      ENDIF
146
 
      JAMP1(1)=+1D0/2D0*(+1D0/3D0*AMP(1))
147
 
      JAMP1(2)=+1D0/2D0*(-AMP(1))
148
 
      JAMP2(1)=+1D0/4D0*(-AMP(1)-1D0/9D0*AMP(1))
149
 
      JAMP2(2)=+1D0/2D0*(+1D0/3D0*AMP(1))
150
 
      B_SF_003 = 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)
 
240
C     JAMP1s contributing to orders QCD=2 QED=0
 
241
      JAMP1(1,1)=+1D0/2D0*(+1D0/3D0*AMP(1))
 
242
      JAMP1(2,1)=+1D0/2D0*(-AMP(1))
 
243
C     JAMP2s contributing to orders QCD=2 QED=0
 
244
      JAMP2(1,1)=+1D0/4D0*(-AMP(1)-1D0/9D0*AMP(1))
 
245
      JAMP2(2,1)=+1D0/2D0*(+1D0/3D0*AMP(1))
 
246
      DO I = 1, NSQAMPSO
 
247
        ANS(I) = 0D0
 
248
      ENDDO
 
249
      DO M = 1, NAMPSO
 
250
        DO I = 1, NCOLOR1
 
251
          ZTEMP = (0.D0,0.D0)
 
252
          DO J = 1, NCOLOR2
 
253
            ZTEMP = ZTEMP + CF(J,I)*JAMP2(J,M)
 
254
          ENDDO
 
255
          DO N = 1, NAMPSO
 
256
            ANS(SQSOINDEXB(M,N))=ANS(SQSOINDEXB(M,N))+ZTEMP
 
257
     $       *DCONJG(JAMP1(I,N))/DENOM(I)
 
258
          ENDDO
155
259
        ENDDO
156
 
        B_SF_003 =B_SF_003+ZTEMP*DCONJG(JAMP1(I))/DENOM(I)
157
260
      ENDDO
158
261
      END
159
262