← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_ppw_fksall/%SubProcesses%P0_dxu_wp%parton_lum_1.f

  • Committer: olivier Mattelaer
  • Date: 2016-05-12 11:00:18 UTC
  • mfrom: (262.1.150 2.3.4)
  • Revision ID: olivier.mattelaer@uclouvain.be-20160512110018-sevb79f0wm4g8mpp
pass to 2.4.0

Show diffs side-by-side

added added

removed removed

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