← Back to branch summary

~maddevelopers/mg5amcnlo/FKS_EW_granny-run_py8_at_evtgen

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%parton_lum_1.f

  • Committer: Marco Zaro
  • Date: 2017-03-08 10:31:32 UTC
  • Revision ID: marco.zaro@gmail.com-20170308103132-vsnzjovp9f90gc4y
fixing tests

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/IOExportFKSTest/test_pptt_fksreal/%SubProcesses%P0_uux_ttx%parton_lum_1.f
7
7
C     RETURNS PARTON LUMINOSITIES FOR MADFKS                          
8
8
C        
9
9
C     
10
 
C     Process: u u~ > t t~ g WEIGHTED<=3 [ real = QCD ]
11
 
C     Process: c c~ > t t~ g WEIGHTED<=3 [ real = QCD ]
12
 
C     Process: d d~ > t t~ g WEIGHTED<=3 [ real = QCD ]
13
 
C     Process: s s~ > t t~ g WEIGHTED<=3 [ real = QCD ]
 
10
C     Process: u u~ > t t~ g [ real = QED QCD ] QCD^2=6 QED^2=0
 
11
C     Process: c c~ > t t~ g [ real = QED QCD ] QCD^2=6 QED^2=0
 
12
C     Process: d d~ > t t~ g [ real = QED QCD ] QCD^2=6 QED^2=0
 
13
C     Process: s s~ > t t~ g [ real = QED QCD ] QCD^2=6 QED^2=0
14
14
C     
15
15
C     ****************************************************            
16
16
C         
27
27
C     ARGUMENTS                                                       
28
28
C         
29
29
C     
30
 
      DOUBLE PRECISION PP(0:3,NEXTERNAL), LUM
 
30
      DOUBLE PRECISION LUM
31
31
C     
32
32
C     LOCAL VARIABLES                                                 
33
33
C         
34
34
C     
35
 
      INTEGER I, ICROSS,ITYPE,LP
36
 
      DOUBLE PRECISION P1(0:3,NEXTERNAL)
 
35
      INTEGER I, ICROSS,LP
37
36
      DOUBLE PRECISION D1,U1,S1,C1
38
37
      DOUBLE PRECISION CX2,SX2,UX2,DX2
39
 
      DOUBLE PRECISION XPQ(-7:7)
40
38
C     
41
39
C     EXTERNAL FUNCTIONS                                              
42
40
C         
43
41
C     
44
 
      DOUBLE PRECISION ALPHAS2,REWGT,PDG2PDF
 
42
      DOUBLE PRECISION PDG2PDF
45
43
C     
46
44
C     GLOBAL VARIABLES                                                
47
45
C         
67
65
C     ----------                                                      
68
66
C         
69
67
      LUM = 0D0
70
 
      IF (IMIRROR.EQ.2) THEN
71
 
        IF (ABS(LPP(2)) .GE. 1) THEN
72
 
          LP=SIGN(1,LPP(2))
73
 
          D1=PDG2PDF(ABS(LPP(2)),1*LP,XBK(2),DSQRT(Q2FACT(2)))
74
 
          U1=PDG2PDF(ABS(LPP(2)),2*LP,XBK(2),DSQRT(Q2FACT(2)))
75
 
          S1=PDG2PDF(ABS(LPP(2)),3*LP,XBK(2),DSQRT(Q2FACT(2)))
76
 
          C1=PDG2PDF(ABS(LPP(2)),4*LP,XBK(2),DSQRT(Q2FACT(2)))
77
 
        ENDIF
78
 
        IF (ABS(LPP(1)) .GE. 1) THEN
79
 
          LP=SIGN(1,LPP(1))
80
 
          CX2=PDG2PDF(ABS(LPP(1)),-4*LP,XBK(1),DSQRT(Q2FACT(1)))
81
 
          SX2=PDG2PDF(ABS(LPP(1)),-3*LP,XBK(1),DSQRT(Q2FACT(1)))
82
 
          UX2=PDG2PDF(ABS(LPP(1)),-2*LP,XBK(1),DSQRT(Q2FACT(1)))
83
 
          DX2=PDG2PDF(ABS(LPP(1)),-1*LP,XBK(1),DSQRT(Q2FACT(1)))
84
 
        ENDIF
85
 
        PD(0) = 0D0
86
 
        IPROC = 0
87
 
        IPROC=IPROC+1  ! u u~ > t t~ g
88
 
        PD(IPROC) = U1*UX2
89
 
        IPROC=IPROC+1  ! c c~ > t t~ g
90
 
        PD(IPROC) = C1*CX2
91
 
        IPROC=IPROC+1  ! d d~ > t t~ g
92
 
        PD(IPROC) = D1*DX2
93
 
        IPROC=IPROC+1  ! s s~ > t t~ g
94
 
        PD(IPROC) = S1*SX2
95
 
      ELSE
96
 
        IF (ABS(LPP(1)) .GE. 1) THEN
97
 
          LP=SIGN(1,LPP(1))
98
 
          D1=PDG2PDF(ABS(LPP(1)),1*LP,XBK(1),DSQRT(Q2FACT(1)))
99
 
          U1=PDG2PDF(ABS(LPP(1)),2*LP,XBK(1),DSQRT(Q2FACT(1)))
100
 
          S1=PDG2PDF(ABS(LPP(1)),3*LP,XBK(1),DSQRT(Q2FACT(1)))
101
 
          C1=PDG2PDF(ABS(LPP(1)),4*LP,XBK(1),DSQRT(Q2FACT(1)))
102
 
        ENDIF
103
 
        IF (ABS(LPP(2)) .GE. 1) THEN
104
 
          LP=SIGN(1,LPP(2))
105
 
          CX2=PDG2PDF(ABS(LPP(2)),-4*LP,XBK(2),DSQRT(Q2FACT(2)))
106
 
          SX2=PDG2PDF(ABS(LPP(2)),-3*LP,XBK(2),DSQRT(Q2FACT(2)))
107
 
          UX2=PDG2PDF(ABS(LPP(2)),-2*LP,XBK(2),DSQRT(Q2FACT(2)))
108
 
          DX2=PDG2PDF(ABS(LPP(2)),-1*LP,XBK(2),DSQRT(Q2FACT(2)))
109
 
        ENDIF
110
 
        PD(0) = 0D0
111
 
        IPROC = 0
112
 
        IPROC=IPROC+1  ! u u~ > t t~ g
113
 
        PD(IPROC) = U1*UX2
114
 
        IPROC=IPROC+1  ! c c~ > t t~ g
115
 
        PD(IPROC) = C1*CX2
116
 
        IPROC=IPROC+1  ! d d~ > t t~ g
117
 
        PD(IPROC) = D1*DX2
118
 
        IPROC=IPROC+1  ! s s~ > t t~ g
119
 
        PD(IPROC) = S1*SX2
120
 
      ENDIF
 
68
      IF (ABS(LPP(1)) .GE. 1) THEN
 
69
        LP=SIGN(1,LPP(1))
 
70
        D1=PDG2PDF(ABS(LPP(1)),1*LP,XBK(1),DSQRT(Q2FACT(1)))
 
71
        U1=PDG2PDF(ABS(LPP(1)),2*LP,XBK(1),DSQRT(Q2FACT(1)))
 
72
        S1=PDG2PDF(ABS(LPP(1)),3*LP,XBK(1),DSQRT(Q2FACT(1)))
 
73
        C1=PDG2PDF(ABS(LPP(1)),4*LP,XBK(1),DSQRT(Q2FACT(1)))
 
74
      ENDIF
 
75
      IF (ABS(LPP(2)) .GE. 1) THEN
 
76
        LP=SIGN(1,LPP(2))
 
77
        CX2=PDG2PDF(ABS(LPP(2)),-4*LP,XBK(2),DSQRT(Q2FACT(2)))
 
78
        SX2=PDG2PDF(ABS(LPP(2)),-3*LP,XBK(2),DSQRT(Q2FACT(2)))
 
79
        UX2=PDG2PDF(ABS(LPP(2)),-2*LP,XBK(2),DSQRT(Q2FACT(2)))
 
80
        DX2=PDG2PDF(ABS(LPP(2)),-1*LP,XBK(2),DSQRT(Q2FACT(2)))
 
81
      ENDIF
 
82
      PD(0) = 0D0
 
83
      IPROC = 0
 
84
      IPROC=IPROC+1  ! u u~ > t t~ g
 
85
      PD(IPROC) = U1*UX2
 
86
      IPROC=IPROC+1  ! c c~ > t t~ g
 
87
      PD(IPROC) = C1*CX2
 
88
      IPROC=IPROC+1  ! d d~ > t t~ g
 
89
      PD(IPROC) = D1*DX2
 
90
      IPROC=IPROC+1  ! s s~ > t t~ g
 
91
      PD(IPROC) = S1*SX2
121
92
      DO I=1,IPROC
122
93
        IF (NINCOMING.EQ.2) THEN
123
94
          LUM = LUM + PD(I) * CONV