5 SUBROUTINE PARTON(F,X1,X2,QQ)
6 IMPLICIT REAL*8(A-H,O-Z)
9 #include "njet_hijing.inc"
14 S=DLOG(DLOG(QQ/DLAM**2)/DLOG(Q0**2/DLAM**2))
15 IF(IHPR2(7).EQ.2) GO TO 200
16 C*******************************************************
17 AT1=0.419+0.004*S-0.007*S**2
18 AT2=3.460+0.724*S-0.066*S**2
19 GMUD=4.40-4.86*S+1.33*S**2
20 AT3=0.763-0.237*S+0.026*S**2
21 AT4=4.00+0.627*S-0.019*S**2
22 GMD=-0.421*S+0.033*S**2
23 C*******************************************************
24 CAS=1.265-1.132*S+0.293*S**2
25 AS=-0.372*S-0.029*S**2
26 BS=8.05+1.59*S-0.153*S**2
27 APHS=6.31*S-0.273*S**2
28 BTAS=-10.5*S-3.17*S**2
30 C********************************************************
31 C CAC=0.135*S-0.075*S**2
32 C AC=-0.036-0.222*S-0.058*S**2
33 C BC=6.35+3.26*S-0.909*S**2
34 C APHC=-3.03*S+1.50*S**2
35 C BTAC=17.4*S-11.3*S**2
36 C GMC=-17.9*S+15.6*S**2
37 C***********************************************************
38 CAG=1.56-1.71*S+0.638*S**2
39 AG=-0.949*S+0.325*S**2
40 BG=6.0+1.44*S-1.05*S**2
41 APHG=9.0-7.19*S+0.255*S**2
42 BTAG=-16.5*S+10.9*S**2
45 C********************************************************
47 AT2=3.33+0.753*S-0.076*S**2
48 GMUD=6.03-6.22*S+1.56*S**2
49 AT3=0.761-0.232*S+0.023*S**2
50 AT4=3.83+0.627*S-0.019*S**2
51 GMD=-0.418*S+0.036*S**2
52 C************************************
53 CAS=1.67-1.92*S+0.582*S**2
54 AS=-0.273*S-0.164*S**2
55 BS=9.15+0.530*S-0.763*S**2
57 BTAS=-101.0*S+44.7*S**2
58 GMS=223.0*S-117.0*S**2
59 C*********************************
60 C CAC=0.067*S-0.031*S**2
61 C AC=-0.120-0.233*S-0.023*S**2
62 C BC=3.51+3.66*S-0.453*S**2
63 C APHC=-0.474*S+0.358*S**2
64 C BTAC=9.50*S-5.43*S**2
65 C GMC=-16.6*S+15.5*S**2
66 C**********************************
67 CAG=0.879-0.971*S+0.434*S**2
69 BG=4.0+1.23*S-0.254*S**2
70 APHG=9.0-5.64*S-0.817*S**2
71 BTAG=-7.54*S+5.50*S**2
72 GMG=-0.596*S+1.26*S**2
73 C*********************************
74 300 B12=DEXP(GMRE(AT1)+GMRE(AT2+1.D0)-GMRE(AT1+AT2+1.D0))
75 B34=DEXP(GMRE(AT3)+GMRE(AT4+1.D0)-GMRE(AT3+AT4+1.D0))
76 CNUD=3.D0/B12/(1.D0+GMUD*AT1/(AT1+AT2+1.D0))
77 CND=1.D0/B34/(1.D0+GMD*AT3/(AT3+AT4+1.D0))
78 C********************************************************
80 C FS=X*2(UBAR+DBAR+SBAR) AND UBAR=DBAR=SBAR
81 C*******************************************************
82 C PH Make sure X1 belongs to [0,1]
83 IF (X1.LT.0.D0) X1=0.D0
84 IF (X1.GT.1.D0) X1=1.D0
85 FUD1=CNUD*X1**AT1*(1.D0-X1)**AT2*(1.D0+GMUD*X1)
86 FS1=CAS*X1**AS*(1.D0-X1)**BS*(1.D0+APHS*X1
87 & +BTAS*X1**2+GMS*X1**3)
88 F(1,3)=CND*X1**AT3*(1.D0-X1)**AT4*(1.D0+GMD*X1)+FS1/6.D0
89 F(1,1)=FUD1-F(1,3)+FS1/3.D0
94 F(1,7)=CAG*X1**AG*(1.D0-X1)**BG*(1.D0+APHG*X1
95 & +BTAG*X1**2+GMG*X1**3)
97 C PH Make sure X2 belongs to [0,1]
98 IF (X2.LT.0.D0) X2=0.D0
99 IF (X2.GT.1.D0) X2=1.D0
100 FUD2=CNUD*X2**AT1*(1.D0-X2)**AT2*(1.D0+GMUD*X2)
101 FS2=CAS*X2**AS*(1.D0-X2)**BS*(1.D0+APHS*X2
102 & +BTAS*X2**2+GMS*X2**3)
103 F(2,3)=CND*X2**AT3*(1.D0-X2)**AT4*(1.D0+GMD*X2)+FS2/6.D0
104 F(2,1)=FUD2-F(2,3)+FS2/3.D0
109 F(2,7)=CAG*X2**AG*(1.D0-X2)**BG*(1.D0+APHG*X2
110 & +BTAG*X2**2+GMG*X2**3)
111 C***********Nuclear effect on the structure function****************
113 IF(IHPR2(6).EQ.1 .AND. IHNT2(1).GT.1) THEN
114 AAX=1.193*ALOG(FLOAT(IHNT2(1)))**0.16666666
115 RRX=AAX*(X1**3-1.2*X1**2+0.21*X1)+1.0
116 & +1.079*(FLOAT(IHNT2(1))**0.33333333-1.0)
117 & /DLOG(IHNT2(1)+1.0D0)*DSQRT(X1)*DEXP(-X1**2/0.01)
118 IF(IP_CRS.EQ.1 .OR.IP_CRS.EQ.3) RRX=DEXP(-X1**2/0.01)
123 IF(IHPR2(6).EQ.1 .AND. IHNT2(3).GT.1) THEN
124 AAX=1.193*ALOG(FLOAT(IHNT2(3)))**0.16666666
125 RRX=AAX*(X2**3-1.2*X2**2+0.21*X2)+1.0
126 & +1.079*(FLOAT(IHNT2(3))**0.33333-1.0)
127 & /DLOG(IHNT2(3)+1.0D0)*DSQRT(X2)*DEXP(-X2**2/0.01)
128 IF(IP_CRS.EQ.2 .OR. IP_CRS.EQ.3) RRX=DEXP(-X2**2/0.01)