~maddevelopers/mg5amcnlo/2.7.3-lepcoll

« back to all changes in this revision

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

• browse files at revision 308
• compare with another revision
• download diff
• download tarball
• view history from revision 308

Show diffs side-by-side

added

removed

tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%b_sf_006.f

1

 

      SUBROUTINE SB_SF_006(P1,ANS)

 

1

      SUBROUTINE SB_SF_006(P,ANS_SUMMED)

 

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     

 

8

C     Return the sum of the split orders which are required in

 

9

C      orders.inc (BORN_ORDERS)

 

10

C     Also the values needed for the counterterms are stored in the

 

11

C      C_BORN_CNT common block

 

12

C     

 

13

C     

 

14

C     Process: g g > t t~ [ real = QCD QED ] QCD^2<=6 QED^2<=0

 

15

C     spectators: 3 3 

 

16

 

 

17

C     

 

18

C     

 

19

C     CONSTANTS

 

20

C     

 

21

      IMPLICIT NONE

 

22

      INCLUDE 'nexternal.inc'

 

23

      INTEGER NSQAMPSO

 

24

      PARAMETER (NSQAMPSO=1)

 

25

C     

 

26

C     ARGUMENTS 

 

27

C     

 

28

      REAL*8 P(0:3,NEXTERNAL), ANS_SUMMED

 

29

C     

 

30

C     VARIABLES

 

31

C     

 

32

      INTEGER I,J

 

33

      INCLUDE 'orders.inc'

 

34

      REAL*8 ANS(0:NSQAMPSO)

 

35

      LOGICAL KEEP_ORDER_CNT(NSPLITORDERS, NSQAMPSO)

 

36

      COMMON /C_KEEP_ORDER_CNT/ KEEP_ORDER_CNT

 

37

      INTEGER AMP_ORDERS(NSPLITORDERS)

 

38

      DOUBLE PRECISION TINY

 

39

      PARAMETER (TINY = 1D-12)

 

40

      DOUBLE PRECISION MAX_VAL

 

41

C     

 

42

C     FUNCTIONS

 

43

C     

 

44

      INTEGER GETORDPOWFROMINDEX_B

 

45

      INTEGER ORDERS_TO_AMP_SPLIT_POS

 

46

C     

 

47

C     BEGIN CODE

 

48

C     

 

49

      CALL SB_SF_006_SPLITORDERS(P,ANS)

 

50

C     color-linked borns are called for QCD-type emissions

 

51

      ANS_SUMMED = 0D0

 

52

      MAX_VAL = 0D0

 

53

 

 

54

C     reset the amp_split_cnt array

 

55

      AMP_SPLIT_CNT(1:AMP_SPLIT_SIZE,1:2,1:NSPLITORDERS) = DCMPLX(0D0

 

56

     $ ,0D0)

 

57

 

 

58

 

 

59

      DO I = 1, NSQAMPSO

 

60

        MAX_VAL = MAX(MAX_VAL, ABS(ANS(I)))

 

61

      ENDDO

 

62

 

 

63

      DO I = 1, NSQAMPSO

 

64

        IF (KEEP_ORDER_CNT(QCD_POS, I)) THEN

 

65

          ANS_SUMMED = ANS_SUMMED + ANS(I)

 

66

          DO J = 1, NSPLITORDERS

 

67

            AMP_ORDERS(J) = GETORDPOWFROMINDEX_B(J, I)

 

68

C           take into account the fact that this is for QCD

 

69

            IF (J.EQ.QCD_POS) AMP_ORDERS(J) = AMP_ORDERS(J) + 2

 

70

          ENDDO

 

71

            !amp_split_cnt(orders_to_amp_split_pos(amp_orders),1,qcd_pos) = ans(I)

 

72

          IF(ABS(ANS(I)).GT.MAX_VAL*TINY)

 

73

     $      AMP_SPLIT_CNT(ORDERS_TO_AMP_SPLIT_POS(AMP_ORDERS),1

 

74

     $     ,QCD_POS) = ANS(I)

 

75

        ENDIF

 

76

      ENDDO

 

77

 

 

78

C     this is to avoid fake non-zero contributions 

 

79

      IF (ABS(ANS_SUMMED).LT.MAX_VAL*TINY) ANS_SUMMED=0D0

 

80

 

 

81

      RETURN

 

82

      END

 

83

 

 

84

 

 

85

      SUBROUTINE SB_SF_006_SPLITORDERS(P1,ANS)

2

86

C     

3

87

C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s

4

88

C     By the MadGraph5_aMC@NLO Development Team

8

92

C     AND HELICITIES

9

93

C     FOR THE POINT IN PHASE SPACE P(0:3,NEXTERNAL-1)

10

94

C     

11

 

C     Process: g g > t t~ WEIGHTED<=2 [ real = QCD ]

 

95

C     Process: g g > t t~ [ real = QCD QED ] QCD^2<=6 QED^2<=0

12

96

C     spectators: 3 3 

13

97

 

14

98

C     

19

103

      INCLUDE 'nexternal.inc'

20

104

      INTEGER     NCOMB

21

105

      PARAMETER ( NCOMB=  16 )

 

106

      INTEGER NSQAMPSO

 

107

      PARAMETER (NSQAMPSO=1)

22

108

      INTEGER    THEL

23

109

      PARAMETER (THEL=NCOMB*8)

24

110

      INTEGER NGRAPHS

26

112

C     

27

113

C     ARGUMENTS 

28

114

C     

29

 

      REAL*8 P1(0:3,NEXTERNAL-1),ANS

 

115

      REAL*8 P1(0:3,NEXTERNAL-1),ANS(0:NSQAMPSO)

30

116

C     

31

117

C     LOCAL VARIABLES 

32

118

C     

33

 

      INTEGER IHEL,IDEN,J

34

 

      REAL*8 B_SF_006

 

119

      INTEGER IHEL,IDEN,I,J

 

120

      DOUBLE PRECISION T(NSQAMPSO)

35

121

      INTEGER IDEN_VALUES(8)

36

122

      DATA IDEN_VALUES /256, 256, 256, 256, 256, 256, 256, 256/

37

123

C     

64

150

     $   //' called only with calculatedborn = true'

65

151

        STOP

66

152

      ENDIF

67

 

      ANS = 0D0

 

153

      DO I=0,NSQAMPSO

 

154

        ANS(I) = 0D0

 

155

      ENDDO

68

156

      DO IHEL=1,NCOMB

69

157

        IF (GOODHEL(IHEL,NFKSPROCESS)) THEN

70

 

          ANS=ANS+B_SF_006(P1,IHEL)

 

158

          CALL B_SF_006(IHEL,T)

 

159

          DO I=1,NSQAMPSO

 

160

            ANS(I)=ANS(I)+T(I)

 

161

          ENDDO

71

162

        ENDIF

72

163

      ENDDO

73

 

      ANS=ANS/DBLE(IDEN)

 

164

      DO I=1,NSQAMPSO

 

165

        ANS(I)=ANS(I)/DBLE(IDEN)

 

166

        ANS(0)=ANS(0)+ANS(I)

 

167

      ENDDO

74

168

      END

75

169

 

76

170

 

77

 

      REAL*8 FUNCTION B_SF_006(P,HELL)

 

171

      SUBROUTINE B_SF_006(HELL,ANS)

78

172

C     

79

173

C     Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s

80

174

C     By the MadGraph5_aMC@NLO Development Team

82

176

C     RETURNS AMPLITUDE SQUARED SUMMED/AVG OVER COLORS

83

177

C     FOR THE POINT WITH EXTERNAL LINES W(0:6,NEXTERNAL-1)

84

178

 

85

 

C     Process: g g > t t~ WEIGHTED<=2 [ real = QCD ]

 

179

C     Process: g g > t t~ [ real = QCD QED ] QCD^2<=6 QED^2<=0

86

180

C     spectators: 3 3 

87

181

 

88

182

C     

90

184

C     

91

185

C     CONSTANTS

92

186

C     

 

187

      INTEGER NAMPSO, NSQAMPSO

 

188

      PARAMETER (NAMPSO=1, NSQAMPSO=1)

93

189

      INTEGER     NGRAPHS

94

190

      PARAMETER ( NGRAPHS = 3 )

95

191

      INTEGER NCOLOR1, NCOLOR2

103

199

C     

104

200

C     ARGUMENTS 

105

201

C     

106

 

      REAL*8 P(0:3,NEXTERNAL-1)

107

202

      INTEGER HELL

 

203

      REAL*8 ANS(NSQAMPSO)

108

204

C     

109

205

C     LOCAL VARIABLES 

110

206

C     

111

 

      INTEGER I,J

112

 

      REAL*8  CF(NCOLOR2,NCOLOR1)

113

 

      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP1(NCOLOR1), JAMP2(NCOLOR2)

 

207

      INTEGER I,J,M,N

 

208

      REAL*8 CF(NCOLOR2,NCOLOR1)

 

209

      COMPLEX*16 ZTEMP, AMP(NGRAPHS), JAMP1(NCOLOR1,NAMPSO),

 

210

     $  JAMP2(NCOLOR2,NAMPSO)

114

211

      COMPLEX*16 TMP_JAMP(0)

115

212

C     

116

213

C     GLOBAL VARIABLES

120

217

      LOGICAL CALCULATEDBORN

121

218

      COMMON/CCALCULATEDBORN/CALCULATEDBORN

122

219

C     

 

220

C     FUNCTION

 

221

C     

 

222

      INTEGER SQSOINDEXB

 

223

C     

123

224

C     COLOR DATA

124

225

C     

125

226

      DATA (CF(I,  1),I=  1,  2) /5.333333333333333D+00,

138

239

          AMP(I)=SAVEAMP(I,HELL)

139

240

        ENDDO

140

241

      ENDIF

141

 

      JAMP1(1) = ((0.000000000000000D+00,1.000000000000000D+00))*AMP(1)

142

 

     $ +(-1.000000000000000D+00)*AMP(2)

143

 

      JAMP1(2) = ((0.000000000000000D+00,-1.000000000000000D+00))

 

242

C     JAMPs contributing to orders QCD=2 QED=0

 

243

      JAMP1(1,1) = ((0.000000000000000D+00,1.000000000000000D+00))

 

244

     $ *AMP(1)+(-1.000000000000000D+00)*AMP(2)

 

245

      JAMP1(2,1) = ((0.000000000000000D+00,-1.000000000000000D+00))

144

246

     $ *AMP(1)+(-1.000000000000000D+00)*AMP(3)

145

 

      JAMP2(1) = ((0.000000000000000D+00,6.666666666666666D-01))*AMP(1)

146

 

     $ +(-6.666666666666666D-01)*AMP(2)

147

 

      JAMP2(2) = ((0.000000000000000D+00,-6.666666666666666D-01))

 

247

C     JAMPs contributing to orders QCD=2 QED=0

 

248

      JAMP2(1,1) = ((0.000000000000000D+00,6.666666666666666D-01))

 

249

     $ *AMP(1)+(-6.666666666666666D-01)*AMP(2)

 

250

      JAMP2(2,1) = ((0.000000000000000D+00,-6.666666666666666D-01))

148

251

     $ *AMP(1)+(-6.666666666666666D-01)*AMP(3)

149

 

      B_SF_006 = 0.D0

150

 

      DO I = 1, NCOLOR1

151

 

        ZTEMP = (0.D0,0.D0)

152

 

        DO J = 1, NCOLOR2

153

 

          ZTEMP = ZTEMP + CF(J,I)*JAMP2(J)

 

252

      DO I = 1, NSQAMPSO

 

253

        ANS(I) = 0D0

 

254

      ENDDO

 

255

      DO M = 1, NAMPSO

 

256

        DO I = 1, NCOLOR1

 

257

          ZTEMP = (0.D0,0.D0)

 

258

          DO J = 1, NCOLOR2

 

259

            ZTEMP = ZTEMP + CF(J,I)*JAMP2(J,M)

 

260

          ENDDO

 

261

          DO N = 1, NAMPSO

 

262

            ANS(SQSOINDEXB(M,N))=ANS(SQSOINDEXB(M,N))+ZTEMP

 

263

     $       *DCONJG(JAMP1(I,N))

 

264

          ENDDO

154

265

        ENDDO

155

 

        B_SF_006 =B_SF_006+ZTEMP*DCONJG(JAMP1(I))

156

266

      ENDDO

157

267

      END

158

268

 

• Older »

Loggerhead is a web-based interface for Breezy
Version: 2.0.1