← 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_ddx_ttx%parton_lum_3.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_ddx_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: d d~ > t t~ g [ real = QED QCD ] QCD^2<=4 QED^2<=2
 
11
C     Process: s s~ > t t~ g [ 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 D1,S1
 
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 D1_COMPONENTS(N_EE),S1_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 D1,S1/2*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
        D1=PDG2PDF(LPP(1),1,1,XBK(1),DSQRT(Q2FACT(1)))
 
85
        IF ((ABS(LPP(1)).EQ.4.OR.ABS(LPP(1)).EQ.3)
 
86
     $   .AND.PDLABEL.NE.'none') D1_COMPONENTS(1:N_EE) =
 
87
     $    EE_COMPONENTS(1:N_EE)
 
88
        S1=PDG2PDF(LPP(1),3,1,XBK(1),DSQRT(Q2FACT(1)))
 
89
        IF ((ABS(LPP(1)).EQ.4.OR.ABS(LPP(1)).EQ.3)
 
90
     $   .AND.PDLABEL.NE.'none') S1_COMPONENTS(1:N_EE) =
 
91
     $    EE_COMPONENTS(1:N_EE)
 
92
      ENDIF
 
93
      IF (ABS(LPP(2)) .GE. 1) THEN
 
94
        SX2=PDG2PDF(LPP(2),-3,2,XBK(2),DSQRT(Q2FACT(2)))
 
95
        IF ((ABS(LPP(2)).EQ.4.OR.ABS(LPP(2)).EQ.3)
 
96
     $   .AND.PDLABEL.NE.'none') SX2_COMPONENTS(1:N_EE) =
 
97
     $    EE_COMPONENTS(1:N_EE)
 
98
        DX2=PDG2PDF(LPP(2),-1,2,XBK(2),DSQRT(Q2FACT(2)))
 
99
        IF ((ABS(LPP(2)).EQ.4.OR.ABS(LPP(2)).EQ.3)
 
100
     $   .AND.PDLABEL.NE.'none') DX2_COMPONENTS(1:N_EE) =
 
101
     $    EE_COMPONENTS(1:N_EE)
 
102
      ENDIF
 
103
      PD(0) = 0D0
 
104
      IPROC = 0
 
105
      IPROC=IPROC+1  ! d d~ > t t~ g
 
106
      PD(IPROC) = D1*DX2
 
107
      IF (ABS(LPP(1)).EQ.ABS(LPP(2)).AND. (ABS(LPP(1))
 
108
     $ .EQ.3.OR.ABS(LPP(1)).EQ.4).AND.PDLABEL.NE.'none')PD(IPROC)
 
109
     $ =EE_COMP_PROD(D1_COMPONENTS,DX2_COMPONENTS)
 
110
      IPROC=IPROC+1  ! s s~ > t t~ g
 
111
      PD(IPROC) = S1*SX2
 
112
      IF (ABS(LPP(1)).EQ.ABS(LPP(2)).AND. (ABS(LPP(1))
 
113
     $ .EQ.3.OR.ABS(LPP(1)).EQ.4).AND.PDLABEL.NE.'none')PD(IPROC)
 
114
     $ =EE_COMP_PROD(S1_COMPONENTS,SX2_COMPONENTS)
 
115
      DO I=1,IPROC
 
116
        IF (NINCOMING.EQ.2) THEN
 
117
          LUM = LUM + PD(I) * CONV
 
118
        ELSE
 
119
          LUM = LUM + PD(I)
 
120
        ENDIF
 
121
      ENDDO
 
122
      RETURN
 
123
      END
 
124