← Back to branch summary

~maddevelopers/mg5amcnlo/2.7.1.3

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_dxu_wp%b_sf_001.f

  • Committer: olivier Mattelaer
  • Date: 2016-05-12 11:00:18 UTC
  • mfrom: (262.1.150 2.3.4)
  • Revision ID: olivier.mattelaer@uclouvain.be-20160512110018-sevb79f0wm4g8mpp
pass to 2.4.0

Show diffs side-by-side

added added

removed removed

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