← Back to branch summary

~maddevelopers/mg5amcnlo/MS_with_decayevents

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%b_sf_003.f

  • Committer: olivier Mattelaer
  • Date: 2014-03-31 16:58:22 UTC
  • mfrom: (252.1.59 2.1.1)
  • Revision ID: olivier.mattelaer@uclouvain.be-20140331165822-4eb1040586ic1y1v
pass to 2.1.1

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%b_sf_003.f
 
1
      SUBROUTINE SB_SF_003(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 4 
 
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_003
 
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_003(P1,IHEL)
 
74
        ENDIF
 
75
      ENDDO
 
76
      ANS=ANS/DBLE(IDEN)
 
77
      END
 
78
 
 
79
 
 
80
      REAL*8 FUNCTION B_SF_003(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 4 
 
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_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)
 
155
        ENDDO
 
156
        B_SF_003 =B_SF_003+ZTEMP*DCONJG(JAMP1(I))/DENOM(I)
 
157
      ENDDO
 
158
      END
 
159
 
 
160
 
 
161