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
|
SUBROUTINE DLUM_1(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: u~ u > t t~ g [ real = QCD QED ] QCD^2<=4 QED^2<=2
C Process: c~ c > t t~ g [ real = QCD QED ] 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 CX1,UX1
DOUBLE PRECISION U2,C2
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 CX1,UX1/2*1D0/
DATA U2,C2/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))
CX1=PDG2PDF(ABS(LPP(1)),-4*LP,XBK(1),DSQRT(Q2FACT(1)))
UX1=PDG2PDF(ABS(LPP(1)),-2*LP,XBK(1),DSQRT(Q2FACT(1)))
ENDIF
IF (ABS(LPP(2)) .GE. 1) THEN
LP=SIGN(1,LPP(2))
U2=PDG2PDF(ABS(LPP(2)),2*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 ! u~ u > t t~ g
PD(IPROC) = UX1*U2
IPROC=IPROC+1 ! c~ c > t t~ g
PD(IPROC) = CX1*C2
DO I=1,IPROC
IF (NINCOMING.EQ.2) THEN
LUM = LUM + PD(I) * CONV
ELSE
LUM = LUM + PD(I)
ENDIF
ENDDO
RETURN
END
|