← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

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

  • Committer: olivier Mattelaer
  • Date: 2014-06-11 19:47:48 UTC
  • mfrom: (253.1.38 2.1.2)
  • mto: (253.1.70 2.1.2)
  • mto: This revision was merged to the branch mainline in revision 254.
  • Revision ID: olivier.mattelaer@uclouvain.be-20140611194748-nhxx3210gudd1qfw
merge with 2.1.2

Show diffs side-by-side

added added

removed removed

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