← Back to branch summary

~maddevelopers/mg5amcnlo/2.0.2_onshellsubtract_newresh

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%parton_lum_4.f

  • Committer: Marco Zaro
  • Date: 2015-04-09 16:12:12 UTC
  • mfrom: (250.22.93 2.3)
  • Revision ID: marco.zaro@gmail.com-20150409161212-t7umsvk7p4z983rt
merged with 2.3 r345

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_gg_ttx%parton_lum_4.f
 
1
      SUBROUTINE DLUM_4(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: g d~ > t t~ d~ WEIGHTED=3 [ real = QCD ]
 
11
C     Process: g s~ > t t~ s~ WEIGHTED=3 [ real = QCD ]
 
12
C     Process: g u~ > t t~ u~ WEIGHTED=3 [ real = QCD ]
 
13
C     Process: g c~ > t t~ c~ WEIGHTED=3 [ real = 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 G1
 
38
      DOUBLE PRECISION CX2,SX2,UX2,DX2
 
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 G1/1*1D0/
 
61
      DATA CX2,SX2,UX2,DX2/4*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
          G1=PDG2PDF(ABS(LPP(2)),0*LP,XBK(2),DSQRT(Q2FACT(2)))
 
74
        ENDIF
 
75
        IF (ABS(LPP(1)) .GE. 1) THEN
 
76
          LP=SIGN(1,LPP(1))
 
77
          CX2=PDG2PDF(ABS(LPP(1)),-4*LP,XBK(1),DSQRT(Q2FACT(1)))
 
78
          SX2=PDG2PDF(ABS(LPP(1)),-3*LP,XBK(1),DSQRT(Q2FACT(1)))
 
79
          UX2=PDG2PDF(ABS(LPP(1)),-2*LP,XBK(1),DSQRT(Q2FACT(1)))
 
80
          DX2=PDG2PDF(ABS(LPP(1)),-1*LP,XBK(1),DSQRT(Q2FACT(1)))
 
81
        ENDIF
 
82
        PD(0) = 0D0
 
83
        IPROC = 0
 
84
        IPROC=IPROC+1  ! g d~ > t t~ d~
 
85
        PD(IPROC) = G1*DX2
 
86
        IPROC=IPROC+1  ! g s~ > t t~ s~
 
87
        PD(IPROC) = G1*SX2
 
88
        IPROC=IPROC+1  ! g u~ > t t~ u~
 
89
        PD(IPROC) = G1*UX2
 
90
        IPROC=IPROC+1  ! g c~ > t t~ c~
 
91
        PD(IPROC) = G1*CX2
 
92
      ELSE
 
93
        IF (ABS(LPP(1)) .GE. 1) THEN
 
94
          LP=SIGN(1,LPP(1))
 
95
          G1=PDG2PDF(ABS(LPP(1)),0*LP,XBK(1),DSQRT(Q2FACT(1)))
 
96
        ENDIF
 
97
        IF (ABS(LPP(2)) .GE. 1) THEN
 
98
          LP=SIGN(1,LPP(2))
 
99
          CX2=PDG2PDF(ABS(LPP(2)),-4*LP,XBK(2),DSQRT(Q2FACT(2)))
 
100
          SX2=PDG2PDF(ABS(LPP(2)),-3*LP,XBK(2),DSQRT(Q2FACT(2)))
 
101
          UX2=PDG2PDF(ABS(LPP(2)),-2*LP,XBK(2),DSQRT(Q2FACT(2)))
 
102
          DX2=PDG2PDF(ABS(LPP(2)),-1*LP,XBK(2),DSQRT(Q2FACT(2)))
 
103
        ENDIF
 
104
        PD(0) = 0D0
 
105
        IPROC = 0
 
106
        IPROC=IPROC+1  ! g d~ > t t~ d~
 
107
        PD(IPROC) = G1*DX2
 
108
        IPROC=IPROC+1  ! g s~ > t t~ s~
 
109
        PD(IPROC) = G1*SX2
 
110
        IPROC=IPROC+1  ! g u~ > t t~ u~
 
111
        PD(IPROC) = G1*UX2
 
112
        IPROC=IPROC+1  ! g c~ > t t~ c~
 
113
        PD(IPROC) = G1*CX2
 
114
      ENDIF
 
115
      DO I=1,IPROC
 
116
        LUM = LUM + PD(I) * CONV
 
117
      ENDDO
 
118
      RETURN
 
119
      END
 
120