1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
|
SUBROUTINE DLUM_2(LUM)
C ****************************************************
C
C Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
C By the MadGraph5_aMC@NLO Development Team
C Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
C RETURNS PARTON LUMINOSITIES FOR MADFKS
C
C
C Process: g u > t t~ u WEIGHTED=3 [ real = QCD ]
C Process: g c > t t~ c WEIGHTED=3 [ real = QCD ]
C Process: g d > t t~ d WEIGHTED=3 [ real = QCD ]
C Process: g s > t t~ s WEIGHTED=3 [ real = QCD ]
C
C ****************************************************
C
IMPLICIT NONE
C
C CONSTANTS
C
C
INCLUDE 'genps.inc'
INCLUDE 'nexternal.inc'
DOUBLE PRECISION CONV
PARAMETER (CONV=389379660D0) !CONV TO PICOBARNS
C
C ARGUMENTS
C
C
DOUBLE PRECISION PP(0:3,NEXTERNAL), LUM
C
C LOCAL VARIABLES
C
C
INTEGER I, ICROSS,ITYPE,LP
DOUBLE PRECISION P1(0:3,NEXTERNAL)
DOUBLE PRECISION G1
DOUBLE PRECISION D2,U2,S2,C2
DOUBLE PRECISION XPQ(-7:7)
C
C EXTERNAL FUNCTIONS
C
C
DOUBLE PRECISION ALPHAS2,REWGT,PDG2PDF
C
C GLOBAL VARIABLES
C
C
INTEGER IPROC
DOUBLE PRECISION PD(0:MAXPROC)
COMMON /SUBPROC/ PD, IPROC
INCLUDE 'coupl.inc'
INCLUDE 'run.inc'
INTEGER IMIRROR
COMMON/CMIRROR/IMIRROR
C
C DATA
C
C
DATA G1/1*1D0/
DATA D2,U2,S2,C2/4*1D0/
DATA ICROSS/1/
C ----------
C
C BEGIN CODE
C
C ----------
C
LUM = 0D0
IF (IMIRROR.EQ.2) THEN
IF (ABS(LPP(2)) .GE. 1) THEN
LP=SIGN(1,LPP(2))
G1=PDG2PDF(ABS(LPP(2)),0*LP,XBK(2),DSQRT(Q2FACT(2)))
ENDIF
IF (ABS(LPP(1)) .GE. 1) THEN
LP=SIGN(1,LPP(1))
D2=PDG2PDF(ABS(LPP(1)),1*LP,XBK(1),DSQRT(Q2FACT(1)))
U2=PDG2PDF(ABS(LPP(1)),2*LP,XBK(1),DSQRT(Q2FACT(1)))
S2=PDG2PDF(ABS(LPP(1)),3*LP,XBK(1),DSQRT(Q2FACT(1)))
C2=PDG2PDF(ABS(LPP(1)),4*LP,XBK(1),DSQRT(Q2FACT(1)))
ENDIF
PD(0) = 0D0
IPROC = 0
IPROC=IPROC+1 ! g u > t t~ u
PD(IPROC) = G1*U2
IPROC=IPROC+1 ! g c > t t~ c
PD(IPROC) = G1*C2
IPROC=IPROC+1 ! g d > t t~ d
PD(IPROC) = G1*D2
IPROC=IPROC+1 ! g s > t t~ s
PD(IPROC) = G1*S2
ELSE
IF (ABS(LPP(1)) .GE. 1) THEN
LP=SIGN(1,LPP(1))
G1=PDG2PDF(ABS(LPP(1)),0*LP,XBK(1),DSQRT(Q2FACT(1)))
ENDIF
IF (ABS(LPP(2)) .GE. 1) THEN
LP=SIGN(1,LPP(2))
D2=PDG2PDF(ABS(LPP(2)),1*LP,XBK(2),DSQRT(Q2FACT(2)))
U2=PDG2PDF(ABS(LPP(2)),2*LP,XBK(2),DSQRT(Q2FACT(2)))
S2=PDG2PDF(ABS(LPP(2)),3*LP,XBK(2),DSQRT(Q2FACT(2)))
C2=PDG2PDF(ABS(LPP(2)),4*LP,XBK(2),DSQRT(Q2FACT(2)))
ENDIF
PD(0) = 0D0
IPROC = 0
IPROC=IPROC+1 ! g u > t t~ u
PD(IPROC) = G1*U2
IPROC=IPROC+1 ! g c > t t~ c
PD(IPROC) = G1*C2
IPROC=IPROC+1 ! g d > t t~ d
PD(IPROC) = G1*D2
IPROC=IPROC+1 ! g s > t t~ s
PD(IPROC) = G1*S2
ENDIF
DO I=1,IPROC
LUM = LUM + PD(I) * CONV
ENDDO
RETURN
END
|