← Back to branch summary

~maddevelopers/mg5amcnlo/2.5.3

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uxu_ttx%b_sf_002.f

  • Committer: olivier-mattelaer
  • Date: 2017-03-08 12:31:17 UTC
  • Revision ID: olivier-mattelaer-20170308123117-h0zkqjyh9sihsc61
empty version to have an effective freeze of the code

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uxu_ttx%b_sf_002.f
1
 
      SUBROUTINE SB_SF_002(P1,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-1)
10
 
C     
11
 
C     Process: u~ u > t t~ WEIGHTED=2 [ QCD ]
12
 
C     Process: c~ c > t t~ WEIGHTED=2 [ QCD ]
13
 
C     Process: d~ d > t t~ WEIGHTED=2 [ QCD ]
14
 
C     Process: s~ s > t t~ WEIGHTED=2 [ QCD ]
15
 
C     spectators: 1 3 
16
 
 
17
 
C     
18
 
      IMPLICIT NONE
19
 
C     
20
 
C     CONSTANTS
21
 
C     
22
 
      INCLUDE 'nexternal.inc'
23
 
      INTEGER     NCOMB
24
 
      PARAMETER ( NCOMB=  16 )
25
 
      INTEGER    THEL
26
 
      PARAMETER (THEL=NCOMB*6)
27
 
      INTEGER NGRAPHS
28
 
      PARAMETER (NGRAPHS=   1)
29
 
C     
30
 
C     ARGUMENTS 
31
 
C     
32
 
      REAL*8 P1(0:3,NEXTERNAL-1),ANS
33
 
C     
34
 
C     LOCAL VARIABLES 
35
 
C     
36
 
      INTEGER IHEL,IDEN,J
37
 
      REAL*8 B_SF_002
38
 
      INTEGER IDEN_VALUES(6)
39
 
      DATA IDEN_VALUES /36, 36, 36, 36, 36, 36/
40
 
C     
41
 
C     GLOBAL VARIABLES
42
 
C     
43
 
      LOGICAL GOODHEL(NCOMB,6)
44
 
      COMMON /C_GOODHEL/ GOODHEL
45
 
      DOUBLE PRECISION SAVEMOM(NEXTERNAL-1,2)
46
 
      COMMON/TO_SAVEMOM/SAVEMOM
47
 
      LOGICAL CALCULATEDBORN
48
 
      COMMON/CCALCULATEDBORN/CALCULATEDBORN
49
 
      INTEGER NFKSPROCESS
50
 
      COMMON/C_NFKSPROCESS/NFKSPROCESS
51
 
C     ----------
52
 
C     BEGIN CODE
53
 
C     ----------
54
 
      IDEN=IDEN_VALUES(NFKSPROCESS)
55
 
      IF (CALCULATEDBORN) THEN
56
 
        DO J=1,NEXTERNAL-1
57
 
          IF (SAVEMOM(J,1).NE.P1(0,J) .OR. SAVEMOM(J,2).NE.P1(3
58
 
     $     ,J)) THEN
59
 
            CALCULATEDBORN=.FALSE.
60
 
            WRITE(*,*) 'Error in sb_sf: momenta not the same in the
61
 
     $        born'
62
 
            STOP
63
 
          ENDIF
64
 
        ENDDO
65
 
      ELSE
66
 
        WRITE(*,*) 'Error in sb_sf: color_linked borns should be
67
 
     $    called only with calculatedborn = true'
68
 
        STOP
69
 
      ENDIF
70
 
      ANS = 0D0
71
 
      DO IHEL=1,NCOMB
72
 
        IF (GOODHEL(IHEL,NFKSPROCESS)) THEN
73
 
          ANS=ANS+B_SF_002(P1,IHEL)
74
 
        ENDIF
75
 
      ENDDO
76
 
      ANS=ANS/DBLE(IDEN)
77
 
      END
78
 
 
79
 
 
80
 
      REAL*8 FUNCTION B_SF_002(P,HELL)
81
 
C     
82
 
C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
83
 
C     By the MadGraph5_aMC@NLO Development Team
84
 
C     Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
85
 
C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS
86
 
C     FOR THE POINT WITH EXTERNAL LINES W(0:6,NEXTERNAL-1)
87
 
 
88
 
C     Process: u~ u > t t~ WEIGHTED=2 [ QCD ]
89
 
C     Process: c~ c > t t~ WEIGHTED=2 [ QCD ]
90
 
C     Process: d~ d > t t~ WEIGHTED=2 [ QCD ]
91
 
C     Process: s~ s > t t~ WEIGHTED=2 [ QCD ]
92
 
C     spectators: 1 3 
93
 
 
94
 
C     
95
 
      IMPLICIT NONE
96
 
C     
97
 
C     CONSTANTS
98
 
C     
99
 
      INTEGER     NGRAPHS
100
 
      PARAMETER ( NGRAPHS = 1 )
101
 
      INTEGER NCOLOR1, NCOLOR2
102
 
      PARAMETER (NCOLOR1=2, NCOLOR2=2)
103
 
      REAL*8     ZERO
104
 
      PARAMETER (ZERO=0D0)
105
 
      COMPLEX*16 IMAG1
106
 
      PARAMETER (IMAG1 = (0D0,1D0))
107
 
      INCLUDE 'nexternal.inc'
108
 
      INCLUDE 'born_nhel.inc'
109
 
C     
110
 
C     ARGUMENTS 
111
 
C     
112
 
      REAL*8 P(0:3,NEXTERNAL-1)
113
 
      INTEGER HELL
114
 
C     
115
 
C     LOCAL VARIABLES 
116
 
C     
117
 
      INTEGER I,J
118
 
      REAL*8 DENOM(NCOLOR1), CF(NCOLOR2,NCOLOR1)
119
 
      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP1(NCOLOR1), JAMP2(NCOLOR2)
120
 
C     
121
 
C     GLOBAL VARIABLES
122
 
C     
123
 
      DOUBLE COMPLEX SAVEAMP(NGRAPHS,MAX_BHEL)
124
 
      COMMON/TO_SAVEAMP/SAVEAMP
125
 
      LOGICAL CALCULATEDBORN
126
 
      COMMON/CCALCULATEDBORN/CALCULATEDBORN
127
 
C     
128
 
C     COLOR DATA
129
 
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/
134
 
C     ----------
135
 
C     BEGIN CODE
136
 
C     ----------
137
 
      IF (.NOT. CALCULATEDBORN) THEN
138
 
        WRITE(*,*) 'Error in b_sf: color_linked borns should be called
139
 
     $    only with calculatedborn = true'
140
 
        STOP
141
 
      ELSEIF (CALCULATEDBORN) THEN
142
 
        DO I=1,NGRAPHS
143
 
          AMP(I)=SAVEAMP(I,HELL)
144
 
        ENDDO
145
 
      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_002 = 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)
155
 
        ENDDO
156
 
        B_SF_002 =B_SF_002+ZTEMP*DCONJG(JAMP1(I))/DENOM(I)
157
 
      ENDDO
158
 
      END
159
 
 
160
 
 
161