← Back to branch summary

~maddevelopers/mg5amcnlo/2.5.3

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fks_loonly/%SubProcesses%P0_gg_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/IOExportFKSTest/test_pptt_fks_loonly/%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 [ LOonly = 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*0)
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(1)
36
 
      DATA IDEN_VALUES / 256 /
37
 
C     
38
 
C     GLOBAL VARIABLES
39
 
C     
40
 
      LOGICAL GOODHEL(NCOMB,0)
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,J))
55
 
     $      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 [ LOonly = 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