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
|
SUBROUTINE DLUM_4(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 d~ > t t~ d~ [ real = QED QCD ] QCD^2<=4 QED^2<=2
C Process: g s~ > t t~ s~ [ real = QED QCD ] QCD^2<=4 QED^2<=2
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 LUM
C
C LOCAL VARIABLES
C
C
INTEGER I, ICROSS,LP
DOUBLE PRECISION G1
DOUBLE PRECISION SX2,DX2
C
C EXTERNAL FUNCTIONS
C
C
DOUBLE PRECISION 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 SX2,DX2/2*1D0/
DATA ICROSS/1/
C ----------
C
C BEGIN CODE
C
C ----------
C
LUM = 0D0
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))
SX2=PDG2PDF(ABS(LPP(2)),-3*LP,XBK(2),DSQRT(Q2FACT(2)))
DX2=PDG2PDF(ABS(LPP(2)),-1*LP,XBK(2),DSQRT(Q2FACT(2)))
ENDIF
PD(0) = 0D0
IPROC = 0
IPROC=IPROC+1 ! g d~ > t t~ d~
PD(IPROC) = G1*DX2
IPROC=IPROC+1 ! g s~ > t t~ s~
PD(IPROC) = G1*SX2
DO I=1,IPROC
IF (NINCOMING.EQ.2) THEN
LUM = LUM + PD(I) * CONV
ELSE
LUM = LUM + PD(I)
ENDIF
ENDDO
RETURN
END
|