2
C ****************************************************
4
C Generated by MadGraph5_aMC@NLO v. %(version)s, %(date)s
5
C By the MadGraph5_aMC@NLO Development Team
6
C Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
7
C RETURNS PARTON LUMINOSITIES FOR MADFKS
10
C Process: d g > t t~ d [ real = QED QCD ] QCD^2<=4 QED^2<=2
11
C Process: s g > t t~ s [ real = QED QCD ] QCD^2<=4 QED^2<=2
13
C ****************************************************
21
INCLUDE 'nexternal.inc'
23
PARAMETER (CONV=389379660D0) !CONV TO PICOBARNS
34
DOUBLE PRECISION D1,S1
40
DOUBLE PRECISION PDG2PDF
46
DOUBLE PRECISION PD(0:MAXPROC)
47
COMMON /SUBPROC/ PD, IPROC
51
COMMON/CMIRROR/IMIRROR
53
C STUFF FOR DRESSED EE COLLISIONS
56
DOUBLE PRECISION EE_COMP_PROD
57
DOUBLE PRECISION DUMMY_COMPONENTS(N_EE)
58
DOUBLE PRECISION D1_COMPONENTS(N_EE),S1_COMPONENTS(N_EE)
59
DOUBLE PRECISION G2_COMPONENTS(N_EE)
66
CHARACTER*7 PDLABEL,EPA_LABEL
68
COMMON/TO_PDF/LHAID,PDLABEL,EPA_LABEL
83
IF (ABS(LPP(1)) .GE. 1) THEN
84
D1=PDG2PDF(LPP(1),1,1,XBK(1),DSQRT(Q2FACT(1)))
85
IF ((ABS(LPP(1)).EQ.4.OR.ABS(LPP(1)).EQ.3)
86
$ .AND.PDLABEL.NE.'none') D1_COMPONENTS(1:N_EE) =
87
$ EE_COMPONENTS(1:N_EE)
88
S1=PDG2PDF(LPP(1),3,1,XBK(1),DSQRT(Q2FACT(1)))
89
IF ((ABS(LPP(1)).EQ.4.OR.ABS(LPP(1)).EQ.3)
90
$ .AND.PDLABEL.NE.'none') S1_COMPONENTS(1:N_EE) =
91
$ EE_COMPONENTS(1:N_EE)
93
IF (ABS(LPP(2)) .GE. 1) THEN
94
G2=PDG2PDF(LPP(2),0,2,XBK(2),DSQRT(Q2FACT(2)))
95
IF ((ABS(LPP(2)).EQ.4.OR.ABS(LPP(2)).EQ.3)
96
$ .AND.PDLABEL.NE.'none') G2_COMPONENTS(1:N_EE) =
97
$ EE_COMPONENTS(1:N_EE)
101
IPROC=IPROC+1 ! d g > t t~ d
103
IF (ABS(LPP(1)).EQ.ABS(LPP(2)).AND. (ABS(LPP(1))
104
$ .EQ.3.OR.ABS(LPP(1)).EQ.4).AND.PDLABEL.NE.'none')PD(IPROC)
105
$ =EE_COMP_PROD(D1_COMPONENTS,G2_COMPONENTS)
106
IPROC=IPROC+1 ! s g > t t~ s
108
IF (ABS(LPP(1)).EQ.ABS(LPP(2)).AND. (ABS(LPP(1))
109
$ .EQ.3.OR.ABS(LPP(1)).EQ.4).AND.PDLABEL.NE.'none')PD(IPROC)
110
$ =EE_COMP_PROD(S1_COMPONENTS,G2_COMPONENTS)
112
IF (NINCOMING.EQ.2) THEN
113
LUM = LUM + PD(I) * CONV