← 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_ag_ttx%parton_lum_1.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_ag_ttx%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: a g > t t~ g [ real = QED QCD ] QCD^2<=4 QED^2<=2
 
11
C     
 
12
C     ****************************************************            
 
13
C         
 
14
      IMPLICIT NONE
 
15
C     
 
16
C     CONSTANTS                                                       
 
17
C         
 
18
C     
 
19
      INCLUDE 'genps.inc'
 
20
      INCLUDE 'nexternal.inc'
 
21
      DOUBLE PRECISION       CONV
 
22
      PARAMETER (CONV=389379660D0)  !CONV TO PICOBARNS             
 
23
C     
 
24
C     ARGUMENTS                                                       
 
25
C         
 
26
C     
 
27
      DOUBLE PRECISION LUM
 
28
C     
 
29
C     LOCAL VARIABLES                                                 
 
30
C         
 
31
C     
 
32
      INTEGER I, ICROSS,LP
 
33
      DOUBLE PRECISION A1
 
34
      DOUBLE PRECISION G2
 
35
C     
 
36
C     EXTERNAL FUNCTIONS                                              
 
37
C         
 
38
C     
 
39
      DOUBLE PRECISION PDG2PDF
 
40
C     
 
41
C     GLOBAL VARIABLES                                                
 
42
C         
 
43
C     
 
44
      INTEGER              IPROC
 
45
      DOUBLE PRECISION PD(0:MAXPROC)
 
46
      COMMON /SUBPROC/ PD, IPROC
 
47
      INCLUDE 'coupl.inc'
 
48
      INCLUDE 'run.inc'
 
49
      INTEGER IMIRROR
 
50
      COMMON/CMIRROR/IMIRROR
 
51
C     
 
52
C     STUFF FOR DRESSED EE COLLISIONS
 
53
C     
 
54
      INCLUDE 'eepdf.inc'
 
55
      DOUBLE PRECISION EE_COMP_PROD
 
56
      DOUBLE PRECISION DUMMY_COMPONENTS(N_EE)
 
57
      DOUBLE PRECISION A1_COMPONENTS(N_EE)
 
58
      DOUBLE PRECISION G2_COMPONENTS(N_EE)
 
59
 
 
60
      INTEGER I_EE
 
61
C     
 
62
C     
 
63
C     
 
64
C     Common blocks
 
65
      CHARACTER*7         PDLABEL,EPA_LABEL
 
66
      INTEGER       LHAID
 
67
      COMMON/TO_PDF/LHAID,PDLABEL,EPA_LABEL
 
68
C     
 
69
C     DATA                                                            
 
70
C         
 
71
C     
 
72
      DATA A1/1*1D0/
 
73
      DATA G2/1*1D0/
 
74
      DATA ICROSS/1/
 
75
C     ----------                                                      
 
76
C         
 
77
C     BEGIN CODE                                                      
 
78
C         
 
79
C     ----------                                                      
 
80
C         
 
81
      LUM = 0D0
 
82
      IF (ABS(LPP(1)) .GE. 1) THEN
 
83
        A1=PDG2PDF(LPP(1),7,1,XBK(1),DSQRT(Q2FACT(1)))
 
84
        IF ((ABS(LPP(1)).EQ.4.OR.ABS(LPP(1)).EQ.3)
 
85
     $   .AND.PDLABEL.NE.'none') A1_COMPONENTS(1:N_EE) =
 
86
     $    EE_COMPONENTS(1:N_EE)
 
87
      ENDIF
 
88
      IF (ABS(LPP(2)) .GE. 1) THEN
 
89
        G2=PDG2PDF(LPP(2),0,2,XBK(2),DSQRT(Q2FACT(2)))
 
90
        IF ((ABS(LPP(2)).EQ.4.OR.ABS(LPP(2)).EQ.3)
 
91
     $   .AND.PDLABEL.NE.'none') G2_COMPONENTS(1:N_EE) =
 
92
     $    EE_COMPONENTS(1:N_EE)
 
93
      ENDIF
 
94
      PD(0) = 0D0
 
95
      IPROC = 0
 
96
      IPROC=IPROC+1  ! a g > t t~ g
 
97
      PD(IPROC) = A1*G2
 
98
      IF (ABS(LPP(1)).EQ.ABS(LPP(2)).AND. (ABS(LPP(1))
 
99
     $ .EQ.3.OR.ABS(LPP(1)).EQ.4).AND.PDLABEL.NE.'none')PD(IPROC)
 
100
     $ =EE_COMP_PROD(A1_COMPONENTS,G2_COMPONENTS)
 
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