← Back to branch summary

~maddevelopers/mg5amcnlo/MS_with_decayevents

« back to all changes in this revision

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

  • Committer: olivier Mattelaer
  • Date: 2014-03-31 16:58:22 UTC
  • mfrom: (252.1.59 2.1.1)
  • Revision ID: olivier.mattelaer@uclouvain.be-20140331165822-4eb1040586ic1y1v
pass to 2.1.1

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uxu_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 CX1,SX1,UX1,DX1
 
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 CX1,SX1,UX1,DX1/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
          CX1=PDG2PDF(ABS(LPP(2)),-4*LP,XBK(2),DSQRT(Q2FACT(2)))
 
74
          SX1=PDG2PDF(ABS(LPP(2)),-3*LP,XBK(2),DSQRT(Q2FACT(2)))
 
75
          UX1=PDG2PDF(ABS(LPP(2)),-2*LP,XBK(2),DSQRT(Q2FACT(2)))
 
76
          DX1=PDG2PDF(ABS(LPP(2)),-1*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) = UX1*G2
 
86
        IPROC=IPROC+1  ! c~ g > t t~ c~
 
87
        PD(IPROC) = CX1*G2
 
88
        IPROC=IPROC+1  ! d~ g > t t~ d~
 
89
        PD(IPROC) = DX1*G2
 
90
        IPROC=IPROC+1  ! s~ g > t t~ s~
 
91
        PD(IPROC) = SX1*G2
 
92
      ELSE
 
93
        IF (ABS(LPP(1)) .GE. 1) THEN
 
94
          LP=SIGN(1,LPP(1))
 
95
          CX1=PDG2PDF(ABS(LPP(1)),-4*LP,XBK(1),DSQRT(Q2FACT(1)))
 
96
          SX1=PDG2PDF(ABS(LPP(1)),-3*LP,XBK(1),DSQRT(Q2FACT(1)))
 
97
          UX1=PDG2PDF(ABS(LPP(1)),-2*LP,XBK(1),DSQRT(Q2FACT(1)))
 
98
          DX1=PDG2PDF(ABS(LPP(1)),-1*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) = UX1*G2
 
108
        IPROC=IPROC+1  ! c~ g > t t~ c~
 
109
        PD(IPROC) = CX1*G2
 
110
        IPROC=IPROC+1  ! d~ g > t t~ d~
 
111
        PD(IPROC) = DX1*G2
 
112
        IPROC=IPROC+1  ! s~ g > t t~ s~
 
113
        PD(IPROC) = SX1*G2
 
114
      ENDIF
 
115
      DO I=1,IPROC
 
116
        LUM = LUM + PD(I) * CONV
 
117
      ENDDO
 
118
      RETURN
 
119
      END
 
120