← 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_dxd_ttx%parton_lum_6.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_dxd_ttx%parton_lum_6.f
7
7
C     RETURNS PARTON LUMINOSITIES FOR MADFKS                          
8
8
C        
9
9
C     
10
 
C     Process: d~ g > t t~ d~ [ real = QED QCD ] QCD^2<=4 QED^2<=2
11
 
C     Process: s~ g > t t~ s~ [ real = QED QCD ] QCD^2<=4 QED^2<=2
 
10
C     Process: d~ a > t t~ d~ [ real = QCD QED ] QCD^2<=4 QED^2<=2
 
11
C     Process: s~ a > t t~ s~ [ real = QCD QED ] QCD^2<=4 QED^2<=2
12
12
C     
13
13
C     ****************************************************            
14
14
C         
32
32
C     
33
33
      INTEGER I, ICROSS,LP
34
34
      DOUBLE PRECISION SX1,DX1
35
 
      DOUBLE PRECISION G2
 
35
      DOUBLE PRECISION A2
36
36
C     
37
37
C     EXTERNAL FUNCTIONS                                              
38
38
C         
56
56
      DOUBLE PRECISION EE_COMP_PROD
57
57
      DOUBLE PRECISION DUMMY_COMPONENTS(N_EE)
58
58
      DOUBLE PRECISION SX1_COMPONENTS(N_EE),DX1_COMPONENTS(N_EE)
59
 
      DOUBLE PRECISION G2_COMPONENTS(N_EE)
 
59
      DOUBLE PRECISION A2_COMPONENTS(N_EE)
60
60
 
61
61
      INTEGER I_EE
62
62
C     
71
71
C         
72
72
C     
73
73
      DATA SX1,DX1/2*1D0/
74
 
      DATA G2/1*1D0/
 
74
      DATA A2/1*1D0/
75
75
      DATA ICROSS/1/
76
76
C     ----------                                                      
77
77
C         
91
91
     $    EE_COMPONENTS(1:N_EE)
92
92
      ENDIF
93
93
      IF (ABS(LPP(2)) .GE. 1) THEN
94
 
        G2=PDG2PDF(LPP(2),0,2,XBK(2),DSQRT(Q2FACT(2)))
 
94
        A2=PDG2PDF(LPP(2),7,2,XBK(2),DSQRT(Q2FACT(2)))
95
95
        IF ((ABS(LPP(2)).EQ.4.OR.ABS(LPP(2)).EQ.3)
96
 
     $   .AND.PDLABEL.NE.'none') G2_COMPONENTS(1:N_EE) =
 
96
     $   .AND.PDLABEL.NE.'none') A2_COMPONENTS(1:N_EE) =
97
97
     $    EE_COMPONENTS(1:N_EE)
98
98
      ENDIF
99
99
      PD(0) = 0D0
100
100
      IPROC = 0
101
 
      IPROC=IPROC+1  ! d~ g > t t~ d~
102
 
      PD(IPROC) = DX1*G2
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(DX1_COMPONENTS,G2_COMPONENTS)
106
 
      IPROC=IPROC+1  ! s~ g > t t~ s~
107
 
      PD(IPROC) = SX1*G2
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(SX1_COMPONENTS,G2_COMPONENTS)
 
101
      IPROC=IPROC+1  ! d~ a > t t~ d~
 
102
      PD(IPROC) = DX1*A2
 
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(DX1_COMPONENTS,A2_COMPONENTS)
 
106
      IPROC=IPROC+1  ! s~ a > t t~ s~
 
107
      PD(IPROC) = SX1*A2
 
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(SX1_COMPONENTS,A2_COMPONENTS)
111
111
      DO I=1,IPROC
112
112
        IF (NINCOMING.EQ.2) THEN
113
113
          LUM = LUM + PD(I) * CONV