← Back to branch summary

~maddevelopers/mg5amcnlo/3.3.0

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksrealew/%SubProcesses%P0_uxu_ttx%parton_lum_4.f

  • Committer: olivier-mattelaer
  • Date: 2021-11-08 09:11:22 UTC
  • Revision ID: olivier-mattelaer-20211108091122-a88ei0x7xo56j55o
update IOTest (now should be reproducible version to version)

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksrealew/%SubProcesses%P0_uxu_ttx%parton_lum_4.f
7
7
C     RETURNS PARTON LUMINOSITIES FOR MADFKS                          
8
8
C        
9
9
C     
10
 
C     Process: g u > t t~ u [ real = QED QCD ] QCD^2<=4 QED^2<=2
11
 
C     Process: g c > t t~ c [ real = QED QCD ] QCD^2<=4 QED^2<=2
 
10
C     Process: a u > t t~ u [ real = QCD QED ] QCD^2<=4 QED^2<=2
 
11
C     Process: a c > t t~ c [ real = QCD QED ] QCD^2<=4 QED^2<=2
12
12
C     
13
13
C     ****************************************************            
14
14
C         
31
31
C         
32
32
C     
33
33
      INTEGER I, ICROSS,LP
34
 
      DOUBLE PRECISION G1
 
34
      DOUBLE PRECISION A1
35
35
      DOUBLE PRECISION U2,C2
36
36
C     
37
37
C     EXTERNAL FUNCTIONS                                              
55
55
      INCLUDE 'eepdf.inc'
56
56
      DOUBLE PRECISION EE_COMP_PROD
57
57
      DOUBLE PRECISION DUMMY_COMPONENTS(N_EE)
58
 
      DOUBLE PRECISION G1_COMPONENTS(N_EE)
 
58
      DOUBLE PRECISION A1_COMPONENTS(N_EE)
59
59
      DOUBLE PRECISION U2_COMPONENTS(N_EE),C2_COMPONENTS(N_EE)
60
60
 
61
61
      INTEGER I_EE
70
70
C     DATA                                                            
71
71
C         
72
72
C     
73
 
      DATA G1/1*1D0/
 
73
      DATA A1/1*1D0/
74
74
      DATA U2,C2/2*1D0/
75
75
      DATA ICROSS/1/
76
76
C     ----------                                                      
81
81
C         
82
82
      LUM = 0D0
83
83
      IF (ABS(LPP(1)) .GE. 1) THEN
84
 
        G1=PDG2PDF(LPP(1),0,1,XBK(1),DSQRT(Q2FACT(1)))
 
84
        A1=PDG2PDF(LPP(1),7,1,XBK(1),DSQRT(Q2FACT(1)))
85
85
        IF ((ABS(LPP(1)).EQ.4.OR.ABS(LPP(1)).EQ.3)
86
 
     $   .AND.PDLABEL.NE.'none') G1_COMPONENTS(1:N_EE) =
 
86
     $   .AND.PDLABEL.NE.'none') A1_COMPONENTS(1:N_EE) =
87
87
     $    EE_COMPONENTS(1:N_EE)
88
88
      ENDIF
89
89
      IF (ABS(LPP(2)) .GE. 1) THEN
98
98
      ENDIF
99
99
      PD(0) = 0D0
100
100
      IPROC = 0
101
 
      IPROC=IPROC+1  ! g u > t t~ u
102
 
      PD(IPROC) = G1*U2
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,U2_COMPONENTS)
106
 
      IPROC=IPROC+1  ! g c > t t~ c
107
 
      PD(IPROC) = G1*C2
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,C2_COMPONENTS)
 
101
      IPROC=IPROC+1  ! a u > t t~ u
 
102
      PD(IPROC) = A1*U2
 
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(A1_COMPONENTS,U2_COMPONENTS)
 
106
      IPROC=IPROC+1  ! a c > t t~ c
 
107
      PD(IPROC) = A1*C2
 
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(A1_COMPONENTS,C2_COMPONENTS)
111
111
      DO I=1,IPROC
112
112
        IF (NINCOMING.EQ.2) THEN
113
113
          LUM = LUM + PD(I) * CONV