← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%parton_lum_3.f

  • Committer: olivier Mattelaer
  • Date: 2014-06-11 19:47:48 UTC
  • mfrom: (253.1.38 2.1.2)
  • mto: (253.1.70 2.1.2)
  • mto: This revision was merged to the branch mainline in revision 254.
  • Revision ID: olivier.mattelaer@uclouvain.be-20140611194748-nhxx3210gudd1qfw
merge with 2.1.2

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%parton_lum_3.f
 
1
      SUBROUTINE DLUM_3(LUM)
 
2
C     ****************************************************            
 
3
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                          
 
8
C        
 
9
C     
 
10
C     Process: u g > t t~ u WEIGHTED=3 [ QCD ]
 
11
C     Process: c g > t t~ c WEIGHTED=3 [ QCD ]
 
12
C     Process: d g > t t~ d WEIGHTED=3 [ QCD ]
 
13
C     Process: s g > t t~ s WEIGHTED=3 [ QCD ]
 
14
C     
 
15
C     ****************************************************            
 
16
C         
 
17
      IMPLICIT NONE
 
18
C     
 
19
C     CONSTANTS                                                       
 
20
C         
 
21
C     
 
22
      INCLUDE 'genps.inc'
 
23
      INCLUDE 'nexternal.inc'
 
24
      DOUBLE PRECISION       CONV
 
25
      PARAMETER (CONV=389379660D0)  !CONV TO PICOBARNS             
 
26
C     
 
27
C     ARGUMENTS                                                       
 
28
C         
 
29
C     
 
30
      DOUBLE PRECISION PP(0:3,NEXTERNAL), LUM
 
31
C     
 
32
C     LOCAL VARIABLES                                                 
 
33
C         
 
34
C     
 
35
      INTEGER I, ICROSS,ITYPE,LP
 
36
      DOUBLE PRECISION P1(0:3,NEXTERNAL)
 
37
      DOUBLE PRECISION D1,U1,S1,C1
 
38
      DOUBLE PRECISION G2
 
39
      DOUBLE PRECISION XPQ(-7:7)
 
40
C     
 
41
C     EXTERNAL FUNCTIONS                                              
 
42
C         
 
43
C     
 
44
      DOUBLE PRECISION ALPHAS2,REWGT,PDG2PDF
 
45
C     
 
46
C     GLOBAL VARIABLES                                                
 
47
C         
 
48
C     
 
49
      INTEGER              IPROC
 
50
      DOUBLE PRECISION PD(0:MAXPROC)
 
51
      COMMON /SUBPROC/ PD, IPROC
 
52
      INCLUDE 'coupl.inc'
 
53
      INCLUDE 'run.inc'
 
54
      INTEGER IMIRROR
 
55
      COMMON/CMIRROR/IMIRROR
 
56
C     
 
57
C     DATA                                                            
 
58
C         
 
59
C     
 
60
      DATA D1,U1,S1,C1/4*1D0/
 
61
      DATA G2/1*1D0/
 
62
      DATA ICROSS/1/
 
63
C     ----------                                                      
 
64
C         
 
65
C     BEGIN CODE                                                      
 
66
C         
 
67
C     ----------                                                      
 
68
C         
 
69
      LUM = 0D0
 
70
      IF (IMIRROR.EQ.2) THEN
 
71
        IF (ABS(LPP(2)) .GE. 1) THEN
 
72
          LP=SIGN(1,LPP(2))
 
73
          D1=PDG2PDF(ABS(LPP(2)),1*LP,XBK(2),DSQRT(Q2FACT(2)))
 
74
          U1=PDG2PDF(ABS(LPP(2)),2*LP,XBK(2),DSQRT(Q2FACT(2)))
 
75
          S1=PDG2PDF(ABS(LPP(2)),3*LP,XBK(2),DSQRT(Q2FACT(2)))
 
76
          C1=PDG2PDF(ABS(LPP(2)),4*LP,XBK(2),DSQRT(Q2FACT(2)))
 
77
        ENDIF
 
78
        IF (ABS(LPP(1)) .GE. 1) THEN
 
79
          LP=SIGN(1,LPP(1))
 
80
          G2=PDG2PDF(ABS(LPP(1)),0*LP,XBK(1),DSQRT(Q2FACT(1)))
 
81
        ENDIF
 
82
        PD(0) = 0D0
 
83
        IPROC = 0
 
84
        IPROC=IPROC+1  ! u g > t t~ u
 
85
        PD(IPROC) = U1*G2
 
86
        IPROC=IPROC+1  ! c g > t t~ c
 
87
        PD(IPROC) = C1*G2
 
88
        IPROC=IPROC+1  ! d g > t t~ d
 
89
        PD(IPROC) = D1*G2
 
90
        IPROC=IPROC+1  ! s g > t t~ s
 
91
        PD(IPROC) = S1*G2
 
92
      ELSE
 
93
        IF (ABS(LPP(1)) .GE. 1) THEN
 
94
          LP=SIGN(1,LPP(1))
 
95
          D1=PDG2PDF(ABS(LPP(1)),1*LP,XBK(1),DSQRT(Q2FACT(1)))
 
96
          U1=PDG2PDF(ABS(LPP(1)),2*LP,XBK(1),DSQRT(Q2FACT(1)))
 
97
          S1=PDG2PDF(ABS(LPP(1)),3*LP,XBK(1),DSQRT(Q2FACT(1)))
 
98
          C1=PDG2PDF(ABS(LPP(1)),4*LP,XBK(1),DSQRT(Q2FACT(1)))
 
99
        ENDIF
 
100
        IF (ABS(LPP(2)) .GE. 1) THEN
 
101
          LP=SIGN(1,LPP(2))
 
102
          G2=PDG2PDF(ABS(LPP(2)),0*LP,XBK(2),DSQRT(Q2FACT(2)))
 
103
        ENDIF
 
104
        PD(0) = 0D0
 
105
        IPROC = 0
 
106
        IPROC=IPROC+1  ! u g > t t~ u
 
107
        PD(IPROC) = U1*G2
 
108
        IPROC=IPROC+1  ! c g > t t~ c
 
109
        PD(IPROC) = C1*G2
 
110
        IPROC=IPROC+1  ! d g > t t~ d
 
111
        PD(IPROC) = D1*G2
 
112
        IPROC=IPROC+1  ! s g > t t~ s
 
113
        PD(IPROC) = S1*G2
 
114
      ENDIF
 
115
      DO I=1,IPROC
 
116
        LUM = LUM + PD(I) * CONV
 
117
      ENDDO
 
118
      RETURN
 
119
      END
 
120