← Back to branch summary

~madteam/mg5amcnlo/3.x

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_wprod_fksew/%SubProcesses%P0_dxu_veep%parton_lum_1.f

  • Committer: olivier-mattelaer
  • Date: 2021-08-23 07:11:46 UTC
  • mfrom: (965.4.17 3.2.0)
  • Revision ID: olivier-mattelaer-20210823071146-a9qwmf4etgzorjdl
pass to 3.2.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/IOExportFKSTest/test_wprod_fksew/%SubProcesses%P0_dxu_veep%parton_lum_1.f
50
50
      INTEGER IMIRROR
51
51
      COMMON/CMIRROR/IMIRROR
52
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 SX1_COMPONENTS(N_EE),DX1_COMPONENTS(N_EE)
 
59
      DOUBLE PRECISION U2_COMPONENTS(N_EE),C2_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     
53
70
C     DATA                                                            
54
71
C         
55
72
C     
64
81
C         
65
82
      LUM = 0D0
66
83
      IF (ABS(LPP(1)) .GE. 1) THEN
67
 
        LP=SIGN(1,LPP(1))
68
 
        SX1=PDG2PDF(ABS(LPP(1)),-3*LP,XBK(1),DSQRT(Q2FACT(1)))
69
 
        DX1=PDG2PDF(ABS(LPP(1)),-1*LP,XBK(1),DSQRT(Q2FACT(1)))
 
84
        SX1=PDG2PDF(LPP(1),-3,1,XBK(1),DSQRT(Q2FACT(1)))
 
85
        IF ((ABS(LPP(1)).EQ.4.OR.ABS(LPP(1)).EQ.3)
 
86
     $   .AND.PDLABEL.NE.'none') SX1_COMPONENTS(1:N_EE) =
 
87
     $    EE_COMPONENTS(1:N_EE)
 
88
        DX1=PDG2PDF(LPP(1),-1,1,XBK(1),DSQRT(Q2FACT(1)))
 
89
        IF ((ABS(LPP(1)).EQ.4.OR.ABS(LPP(1)).EQ.3)
 
90
     $   .AND.PDLABEL.NE.'none') DX1_COMPONENTS(1:N_EE) =
 
91
     $    EE_COMPONENTS(1:N_EE)
70
92
      ENDIF
71
93
      IF (ABS(LPP(2)) .GE. 1) THEN
72
 
        LP=SIGN(1,LPP(2))
73
 
        U2=PDG2PDF(ABS(LPP(2)),2*LP,XBK(2),DSQRT(Q2FACT(2)))
74
 
        C2=PDG2PDF(ABS(LPP(2)),4*LP,XBK(2),DSQRT(Q2FACT(2)))
 
94
        U2=PDG2PDF(LPP(2),2,2,XBK(2),DSQRT(Q2FACT(2)))
 
95
        IF ((ABS(LPP(2)).EQ.4.OR.ABS(LPP(2)).EQ.3)
 
96
     $   .AND.PDLABEL.NE.'none') U2_COMPONENTS(1:N_EE) =
 
97
     $    EE_COMPONENTS(1:N_EE)
 
98
        C2=PDG2PDF(LPP(2),4,2,XBK(2),DSQRT(Q2FACT(2)))
 
99
        IF ((ABS(LPP(2)).EQ.4.OR.ABS(LPP(2)).EQ.3)
 
100
     $   .AND.PDLABEL.NE.'none') C2_COMPONENTS(1:N_EE) =
 
101
     $    EE_COMPONENTS(1:N_EE)
75
102
      ENDIF
76
103
      PD(0) = 0D0
77
104
      IPROC = 0
78
105
      IPROC=IPROC+1  ! d~ u > ve e+ a
79
106
      PD(IPROC) = DX1*U2
 
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(DX1_COMPONENTS,U2_COMPONENTS)
80
110
      IPROC=IPROC+1  ! s~ c > ve e+ a
81
111
      PD(IPROC) = SX1*C2
 
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(SX1_COMPONENTS,C2_COMPONENTS)
82
115
      DO I=1,IPROC
83
116
        IF (NINCOMING.EQ.2) THEN
84
117
          LUM = LUM + PD(I) * CONV