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
121
122
123
124
|
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: d~ u > ve e+ a [ all = QCD QED ] QCD^2<=0 QED^2<=6
C Process: s~ c > ve e+ a [ all = QCD QED ] QCD^2<=0 QED^2<=6
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 SX1,DX1
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 STUFF FOR DRESSED EE COLLISIONS
C
INCLUDE 'eepdf.inc'
DOUBLE PRECISION EE_COMP_PROD
DOUBLE PRECISION DUMMY_COMPONENTS(N_EE)
DOUBLE PRECISION SX1_COMPONENTS(N_EE),DX1_COMPONENTS(N_EE)
DOUBLE PRECISION U2_COMPONENTS(N_EE),C2_COMPONENTS(N_EE)
INTEGER I_EE
C
C
C
C Common blocks
CHARACTER*7 PDLABEL,EPA_LABEL
INTEGER LHAID
COMMON/TO_PDF/LHAID,PDLABEL,EPA_LABEL
C
C DATA
C
C
DATA SX1,DX1/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
SX1=PDG2PDF(LPP(1),-3,1,XBK(1),DSQRT(Q2FACT(1)))
IF ((ABS(LPP(1)).EQ.4.OR.ABS(LPP(1)).EQ.3)
$ .AND.PDLABEL.NE.'none') SX1_COMPONENTS(1:N_EE) =
$ EE_COMPONENTS(1:N_EE)
DX1=PDG2PDF(LPP(1),-1,1,XBK(1),DSQRT(Q2FACT(1)))
IF ((ABS(LPP(1)).EQ.4.OR.ABS(LPP(1)).EQ.3)
$ .AND.PDLABEL.NE.'none') DX1_COMPONENTS(1:N_EE) =
$ EE_COMPONENTS(1:N_EE)
ENDIF
IF (ABS(LPP(2)) .GE. 1) THEN
U2=PDG2PDF(LPP(2),2,2,XBK(2),DSQRT(Q2FACT(2)))
IF ((ABS(LPP(2)).EQ.4.OR.ABS(LPP(2)).EQ.3)
$ .AND.PDLABEL.NE.'none') U2_COMPONENTS(1:N_EE) =
$ EE_COMPONENTS(1:N_EE)
C2=PDG2PDF(LPP(2),4,2,XBK(2),DSQRT(Q2FACT(2)))
IF ((ABS(LPP(2)).EQ.4.OR.ABS(LPP(2)).EQ.3)
$ .AND.PDLABEL.NE.'none') C2_COMPONENTS(1:N_EE) =
$ EE_COMPONENTS(1:N_EE)
ENDIF
PD(0) = 0D0
IPROC = 0
IPROC=IPROC+1 ! d~ u > ve e+ a
PD(IPROC) = DX1*U2
IF (ABS(LPP(1)).EQ.ABS(LPP(2)).AND. (ABS(LPP(1))
$ .EQ.3.OR.ABS(LPP(1)).EQ.4).AND.PDLABEL.NE.'none')PD(IPROC)
$ =EE_COMP_PROD(DX1_COMPONENTS,U2_COMPONENTS)
IPROC=IPROC+1 ! s~ c > ve e+ a
PD(IPROC) = SX1*C2
IF (ABS(LPP(1)).EQ.ABS(LPP(2)).AND. (ABS(LPP(1))
$ .EQ.3.OR.ABS(LPP(1)).EQ.4).AND.PDLABEL.NE.'none')PD(IPROC)
$ =EE_COMP_PROD(SX1_COMPONENTS,C2_COMPONENTS)
DO I=1,IPROC
IF (NINCOMING.EQ.2) THEN
LUM = LUM + PD(I) * CONV
ELSE
LUM = LUM + PD(I)
ENDIF
ENDDO
RETURN
END
|