]>
Commit | Line | Data |
---|---|---|
e74335a4 | 1 | * $Id$ |
2 | C | |
3 | C | |
4 | C | |
5 | SUBROUTINE PARTON(F,X1,X2,QQ) | |
6 | IMPLICIT REAL*8(A-H,O-Z) | |
bc676b8e | 7 | #define BLANKET_SAVE |
e74335a4 | 8 | #include "hiparnt.inc" |
9 | #include "njet_hijing.inc" | |
10 | SAVE | |
11 | DIMENSION F(2,7) | |
12 | DLAM=HIPR1(15) | |
13 | Q0=HIPR1(16) | |
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 | |
29 | GMS=14.7*S+9.80*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 | |
43 | GMG=15.3*S-10.1*S**2 | |
44 | GO TO 300 | |
45 | C******************************************************** | |
46 | 200 AT1=0.374+0.014*S | |
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 | |
56 | APHS=15.7*S-2.83*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 | |
68 | AG=-1.16*S+0.476*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******************************************************** | |
79 | C FUD=X*(U+D) | |
80 | C FS=X*2(UBAR+DBAR+SBAR) AND UBAR=DBAR=SBAR | |
81 | C******************************************************* | |
45ad89de | 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 | |
e74335a4 | 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 | |
90 | F(1,2)=FS1/6.D0 | |
91 | F(1,4)=FS1/6.D0 | |
92 | F(1,5)=FS1/6.D0 | |
93 | F(1,6)=FS1/6.D0 | |
94 | F(1,7)=CAG*X1**AG*(1.D0-X1)**BG*(1.D0+APHG*X1 | |
95 | & +BTAG*X1**2+GMG*X1**3) | |
96 | C | |
45ad89de | 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 | |
e74335a4 | 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 | |
105 | F(2,2)=FS2/6.D0 | |
106 | F(2,4)=FS2/6.D0 | |
107 | F(2,5)=FS2/6.D0 | |
108 | F(2,6)=FS2/6.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**************** | |
112 | C | |
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) | |
119 | DO 400 I=1,7 | |
120 | F(1,I)=RRX*F(1,I) | |
121 | 400 CONTINUE | |
122 | ENDIF | |
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) | |
129 | DO 500 I=1,7 | |
130 | F(2,I)=RRX*F(2,I) | |
131 | 500 CONTINUE | |
132 | ENDIF | |
133 | c | |
134 | RETURN | |
135 | END |