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: u~ a > t t~ u~ [ real = QED QCD ] QCD^2<=4 QED^2<=2
11
C Process: c~ a > t t~ c~ [ 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 CX1,UX1
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 CX1_COMPONENTS(N_EE),UX1_COMPONENTS(N_EE)
59
DOUBLE PRECISION A2_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
CX1=PDG2PDF(LPP(1),-4,1,XBK(1),DSQRT(Q2FACT(1)))
85
IF ((ABS(LPP(1)).EQ.4.OR.ABS(LPP(1)).EQ.3)
86
$ .AND.PDLABEL.NE.'none') CX1_COMPONENTS(1:N_EE) =
87
$ EE_COMPONENTS(1:N_EE)
88
UX1=PDG2PDF(LPP(1),-2,1,XBK(1),DSQRT(Q2FACT(1)))
89
IF ((ABS(LPP(1)).EQ.4.OR.ABS(LPP(1)).EQ.3)
90
$ .AND.PDLABEL.NE.'none') UX1_COMPONENTS(1:N_EE) =
91
$ EE_COMPONENTS(1:N_EE)
93
IF (ABS(LPP(2)) .GE. 1) THEN
94
A2=PDG2PDF(LPP(2),7,2,XBK(2),DSQRT(Q2FACT(2)))
95
IF ((ABS(LPP(2)).EQ.4.OR.ABS(LPP(2)).EQ.3)
96
$ .AND.PDLABEL.NE.'none') A2_COMPONENTS(1:N_EE) =
97
$ EE_COMPONENTS(1:N_EE)
101
IPROC=IPROC+1 ! u~ a > t t~ u~
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(UX1_COMPONENTS,A2_COMPONENTS)
106
IPROC=IPROC+1 ! c~ a > t t~ c~
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(CX1_COMPONENTS,A2_COMPONENTS)
112
IF (NINCOMING.EQ.2) THEN
113
LUM = LUM + PD(I) * CONV