← Back to branch summary

~maddevelopers/mg5amcnlo/new_clustering

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksrealew/%SubProcesses%P0_gg_ttx%parton_lum_5.f

  • Committer: Rikkert Frederix
  • Date: 2021-09-09 15:51:40 UTC
  • mfrom: (78.75.502 3.2.1)
  • Revision ID: frederix@physik.uzh.ch-20210909155140-rg6umfq68h6h47cf
merge with 3.2.1

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksrealew/%SubProcesses%P0_gg_ttx%parton_lum_5.f
 
1
      SUBROUTINE DLUM_5(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~ [ real = QED QCD ] QCD^2<=4 QED^2<=2
 
11
C     Process: g s~ > t t~ s~ [ real = QED QCD ] QCD^2<=4 QED^2<=2
 
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 LUM
 
29
C     
 
30
C     LOCAL VARIABLES                                                 
 
31
C         
 
32
C     
 
33
      INTEGER I, ICROSS,LP
 
34
      DOUBLE PRECISION G1
 
35
      DOUBLE PRECISION SX2,DX2
 
36
C     
 
37
C     EXTERNAL FUNCTIONS                                              
 
38
C         
 
39
C     
 
40
      DOUBLE PRECISION PDG2PDF
 
41
C     
 
42
C     GLOBAL VARIABLES                                                
 
43
C         
 
44
C     
 
45
      INTEGER              IPROC
 
46
      DOUBLE PRECISION PD(0:MAXPROC)
 
47
      COMMON /SUBPROC/ PD, IPROC
 
48
      INCLUDE 'coupl.inc'
 
49
      INCLUDE 'run.inc'
 
50
      INTEGER IMIRROR
 
51
      COMMON/CMIRROR/IMIRROR
 
52
C     
 
53
C     STUFF FOR DRESSED EE COLLISIONS
 
54
C     
 
55
      INCLUDE 'eepdf.inc'
 
56
      DOUBLE PRECISION EE_COMP_PROD
 
57
      DOUBLE PRECISION DUMMY_COMPONENTS(N_EE)
 
58
      DOUBLE PRECISION G1_COMPONENTS(N_EE)
 
59
      DOUBLE PRECISION SX2_COMPONENTS(N_EE),DX2_COMPONENTS(N_EE)
 
60
 
 
61
      INTEGER I_EE
 
62
C     
 
63
C     
 
64
C     
 
65
C     Common blocks
 
66
      CHARACTER*7         PDLABEL,EPA_LABEL
 
67
      INTEGER       LHAID
 
68
      COMMON/TO_PDF/LHAID,PDLABEL,EPA_LABEL
 
69
C     
 
70
C     DATA                                                            
 
71
C         
 
72
C     
 
73
      DATA G1/1*1D0/
 
74
      DATA SX2,DX2/2*1D0/
 
75
      DATA ICROSS/1/
 
76
C     ----------                                                      
 
77
C         
 
78
C     BEGIN CODE                                                      
 
79
C         
 
80
C     ----------                                                      
 
81
C         
 
82
      LUM = 0D0
 
83
      IF (ABS(LPP(1)) .GE. 1) THEN
 
84
        G1=PDG2PDF(LPP(1),0,1,XBK(1),DSQRT(Q2FACT(1)))
 
85
        IF ((ABS(LPP(1)).EQ.4.OR.ABS(LPP(1)).EQ.3)
 
86
     $   .AND.PDLABEL.NE.'none') G1_COMPONENTS(1:N_EE) =
 
87
     $    EE_COMPONENTS(1:N_EE)
 
88
      ENDIF
 
89
      IF (ABS(LPP(2)) .GE. 1) THEN
 
90
        SX2=PDG2PDF(LPP(2),-3,2,XBK(2),DSQRT(Q2FACT(2)))
 
91
        IF ((ABS(LPP(2)).EQ.4.OR.ABS(LPP(2)).EQ.3)
 
92
     $   .AND.PDLABEL.NE.'none') SX2_COMPONENTS(1:N_EE) =
 
93
     $    EE_COMPONENTS(1:N_EE)
 
94
        DX2=PDG2PDF(LPP(2),-1,2,XBK(2),DSQRT(Q2FACT(2)))
 
95
        IF ((ABS(LPP(2)).EQ.4.OR.ABS(LPP(2)).EQ.3)
 
96
     $   .AND.PDLABEL.NE.'none') DX2_COMPONENTS(1:N_EE) =
 
97
     $    EE_COMPONENTS(1:N_EE)
 
98
      ENDIF
 
99
      PD(0) = 0D0
 
100
      IPROC = 0
 
101
      IPROC=IPROC+1  ! g d~ > t t~ d~
 
102
      PD(IPROC) = G1*DX2
 
103
      IF (ABS(LPP(1)).EQ.ABS(LPP(2)).AND. (ABS(LPP(1))
 
104
     $ .EQ.3.OR.ABS(LPP(1)).EQ.4).AND.PDLABEL.NE.'none')PD(IPROC)
 
105
     $ =EE_COMP_PROD(G1_COMPONENTS,DX2_COMPONENTS)
 
106
      IPROC=IPROC+1  ! g s~ > t t~ s~
 
107
      PD(IPROC) = G1*SX2
 
108
      IF (ABS(LPP(1)).EQ.ABS(LPP(2)).AND. (ABS(LPP(1))
 
109
     $ .EQ.3.OR.ABS(LPP(1)).EQ.4).AND.PDLABEL.NE.'none')PD(IPROC)
 
110
     $ =EE_COMP_PROD(G1_COMPONENTS,SX2_COMPONENTS)
 
111
      DO I=1,IPROC
 
112
        IF (NINCOMING.EQ.2) THEN
 
113
          LUM = LUM + PD(I) * CONV
 
114
        ELSE
 
115
          LUM = LUM + PD(I)
 
116
        ENDIF
 
117
      ENDDO
 
118
      RETURN
 
119
      END
 
120