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: g u > t t~ u WEIGHTED=3 [ QCD ]
11
C Process: g c > t t~ c WEIGHTED=3 [ QCD ]
12
C Process: g d > t t~ d WEIGHTED=3 [ QCD ]
13
C Process: g s > t t~ s WEIGHTED=3 [ QCD ]
15
C ****************************************************
23
INCLUDE 'nexternal.inc'
25
PARAMETER (CONV=389379660D0) !CONV TO PICOBARNS
30
DOUBLE PRECISION PP(0:3,NEXTERNAL), LUM
35
INTEGER I, ICROSS,ITYPE,LP
36
DOUBLE PRECISION P1(0:3,NEXTERNAL)
38
DOUBLE PRECISION D2,U2,S2,C2
39
DOUBLE PRECISION XPQ(-7:7)
44
DOUBLE PRECISION ALPHAS2,REWGT,PDG2PDF
50
DOUBLE PRECISION PD(0:MAXPROC)
51
COMMON /SUBPROC/ PD, IPROC
55
COMMON/CMIRROR/IMIRROR
61
DATA D2,U2,S2,C2/4*1D0/
70
IF (IMIRROR.EQ.2) THEN
71
IF (ABS(LPP(2)) .GE. 1) THEN
73
G1=PDG2PDF(ABS(LPP(2)),0*LP,XBK(2),DSQRT(Q2FACT(2)))
75
IF (ABS(LPP(1)) .GE. 1) THEN
77
D2=PDG2PDF(ABS(LPP(1)),1*LP,XBK(1),DSQRT(Q2FACT(1)))
78
U2=PDG2PDF(ABS(LPP(1)),2*LP,XBK(1),DSQRT(Q2FACT(1)))
79
S2=PDG2PDF(ABS(LPP(1)),3*LP,XBK(1),DSQRT(Q2FACT(1)))
80
C2=PDG2PDF(ABS(LPP(1)),4*LP,XBK(1),DSQRT(Q2FACT(1)))
84
IPROC=IPROC+1 ! g u > t t~ u
86
IPROC=IPROC+1 ! g c > t t~ c
88
IPROC=IPROC+1 ! g d > t t~ d
90
IPROC=IPROC+1 ! g s > t t~ s
93
IF (ABS(LPP(1)) .GE. 1) THEN
95
G1=PDG2PDF(ABS(LPP(1)),0*LP,XBK(1),DSQRT(Q2FACT(1)))
97
IF (ABS(LPP(2)) .GE. 1) THEN
99
D2=PDG2PDF(ABS(LPP(2)),1*LP,XBK(2),DSQRT(Q2FACT(2)))
100
U2=PDG2PDF(ABS(LPP(2)),2*LP,XBK(2),DSQRT(Q2FACT(2)))
101
S2=PDG2PDF(ABS(LPP(2)),3*LP,XBK(2),DSQRT(Q2FACT(2)))
102
C2=PDG2PDF(ABS(LPP(2)),4*LP,XBK(2),DSQRT(Q2FACT(2)))
106
IPROC=IPROC+1 ! g u > t t~ u
108
IPROC=IPROC+1 ! g c > t t~ c
110
IPROC=IPROC+1 ! g d > t t~ d
112
IPROC=IPROC+1 ! g s > t t~ s
116
LUM = LUM + PD(I) * CONV