← 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_udx_wp%parton_lum_2.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_udx_wp%parton_lum_2.f
 
1
      SUBROUTINE DLUM_2(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~ > w+ u~ WEIGHTED<=3 [ all = QCD ]
 
11
C     Process: g s~ > w+ c~ 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 G1
 
36
      DOUBLE PRECISION SX2,DX2
 
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 G1/1*1D0/
 
59
      DATA SX2,DX2/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
          G1=PDG2PDF(ABS(LPP(2)),0*LP,XBK(2),DSQRT(Q2FACT(2)))
 
72
        ENDIF
 
73
        IF (ABS(LPP(1)) .GE. 1) THEN
 
74
          LP=SIGN(1,LPP(1))
 
75
          SX2=PDG2PDF(ABS(LPP(1)),-3*LP,XBK(1),DSQRT(Q2FACT(1)))
 
76
          DX2=PDG2PDF(ABS(LPP(1)),-1*LP,XBK(1),DSQRT(Q2FACT(1)))
 
77
        ENDIF
 
78
        PD(0) = 0D0
 
79
        IPROC = 0
 
80
        IPROC=IPROC+1  ! g d~ > w+ u~
 
81
        PD(IPROC) = G1*DX2
 
82
        IPROC=IPROC+1  ! g s~ > w+ c~
 
83
        PD(IPROC) = G1*SX2
 
84
      ELSE
 
85
        IF (ABS(LPP(1)) .GE. 1) THEN
 
86
          LP=SIGN(1,LPP(1))
 
87
          G1=PDG2PDF(ABS(LPP(1)),0*LP,XBK(1),DSQRT(Q2FACT(1)))
 
88
        ENDIF
 
89
        IF (ABS(LPP(2)) .GE. 1) THEN
 
90
          LP=SIGN(1,LPP(2))
 
91
          SX2=PDG2PDF(ABS(LPP(2)),-3*LP,XBK(2),DSQRT(Q2FACT(2)))
 
92
          DX2=PDG2PDF(ABS(LPP(2)),-1*LP,XBK(2),DSQRT(Q2FACT(2)))
 
93
        ENDIF
 
94
        PD(0) = 0D0
 
95
        IPROC = 0
 
96
        IPROC=IPROC+1  ! g d~ > w+ u~
 
97
        PD(IPROC) = G1*DX2
 
98
        IPROC=IPROC+1  ! g s~ > w+ c~
 
99
        PD(IPROC) = G1*SX2
 
100
      ENDIF
 
101
      DO I=1,IPROC
 
102
        IF (NINCOMING.EQ.2) THEN
 
103
          LUM = LUM + PD(I) * CONV
 
104
        ELSE
 
105
          LUM = LUM + PD(I)
 
106
        ENDIF
 
107
      ENDDO
 
108
      RETURN
 
109
      END
 
110