]>
Commit | Line | Data |
---|---|---|
4e9e3152 | 1 | subroutine DOGevolvep0(xin,qin,p2in,ip2in,pdf) |
2 | include 'parmsetup.inc' | |
3 | real*8 xin,qin,q2in,p2in,pdf(-6:6),xval(45),qcdl4,qcdl5 | |
4 | real*8 upv,dnv,usea,dsea,str,chm,bot,top,glu | |
5 | character*16 name(nmxset) | |
6 | integer nmem(nmxset),ndef(nmxset),mmem | |
7 | common/NAME/name,nmem,ndef,mmem | |
8 | integer ns | |
9 | ||
10 | save | |
11 | ||
12 | call DOPHO1(xin,qin,upv,dnv,usea,dsea,str,chm,bot,glu) | |
13 | ||
14 | pdf(-6)= 0.0d0 | |
15 | pdf(6)= 0.0d0 | |
16 | pdf(-5)= bot | |
17 | pdf(5 )= bot | |
18 | pdf(-4)= chm | |
19 | pdf(4 )= chm | |
20 | pdf(-3)= str | |
21 | pdf(3 )= str | |
22 | pdf(-2)= usea | |
23 | pdf(2 )= upv | |
24 | pdf(-1)= dsea | |
25 | pdf(1 )= dnv | |
26 | pdf(0 )= glu | |
27 | ||
28 | return | |
29 | ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc | |
30 | entry DOGevolvep1(xin,qin,p2in,ip2in,pdf) | |
31 | ||
32 | call DOPHO2(xin,qin,upv,dnv,usea,dsea,str,chm,bot,glu) | |
33 | ||
34 | pdf(-6)= 0.0d0 | |
35 | pdf(6)= 0.0d0 | |
36 | pdf(-5)= bot | |
37 | pdf(5 )= bot | |
38 | pdf(-4)= chm | |
39 | pdf(4 )= chm | |
40 | pdf(-3)= str | |
41 | pdf(3 )= str | |
42 | pdf(-2)= usea | |
43 | pdf(2 )= upv | |
44 | pdf(-1)= dsea | |
45 | pdf(1 )= dnv | |
46 | pdf(0 )= glu | |
47 | ||
48 | return | |
49 | ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc | |
50 | entry DOGread(nset) | |
51 | read(1,*)nmem(nset),ndef(nset) | |
52 | return | |
53 | c | |
54 | ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc | |
55 | entry DOGalfa(alfas,qalfa) | |
56 | call getnset(iset) | |
57 | call GetOrderAsM(iset,iord) | |
58 | call Getlam4M(iset,imem,qcdl4) | |
59 | call Getlam5M(iset,imem,qcdl5) | |
60 | call aspdflib(alfas,Qalfa,iord,qcdl5) | |
61 | return | |
62 | c | |
63 | ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc | |
64 | entry DOGinit(Eorder,Q2fit) | |
65 | return | |
66 | c | |
67 | ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc | |
68 | entry DOGpdf(mem) | |
69 | imem = mem | |
70 | return | |
71 | c | |
72 | 1000 format(5e13.5) | |
73 | end | |
74 | ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc | |
75 | SUBROUTINE DOPHO1(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL) | |
76 | C******************************************************************** | |
77 | C* * | |
78 | C* Parametrization of parton distribution functions * | |
79 | C* in the photon (LO analysis) - asymptotic solution of AP eq.! * | |
80 | C* * | |
81 | C* authors: D.Duke and H.Owens (DO) * | |
82 | C* /Phys.Rev. D26 (1982) 1600/ * | |
83 | C* * | |
84 | C* Prepared by: * | |
85 | C* Krzysztof Charchula, DESY * | |
86 | C* bitnet: F1PCHA@DHHDESY3 * | |
87 | C* decnet: 13313::CHARCHULA * | |
88 | C* * | |
89 | C* Modified by: * | |
90 | C* H. Plothow-Besch/CERN-PPE * | |
91 | C* * | |
92 | C******************************************************************** | |
93 | C | |
94 | implicit real*8 (a-h,o-z) | |
95 | double precision | |
96 | + CQ(5), | |
97 | + DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL | |
98 | PARAMETER (ALPEM=7.29927D-3, PI=3.141592D0) | |
99 | PARAMETER (ALAM=0.2D0) | |
100 | DATA CQ/0.33333D0,0.66666D0,0.33333D0,0.66666D0,0.33333D0/ | |
101 | C | |
102 | Q2 = DQ*DQ | |
103 | ALAM2=ALAM**2 | |
104 | FQ=ALPEM/(2.*PI)*LOG(Q2/ALAM2) | |
105 | C | |
106 | C...gluons | |
107 | POMG=0.194*(1.-DX)**1.03/(DX**0.97) | |
108 | DGL=POMG*FQ | |
109 | C | |
110 | C...quarks | |
111 | POM1=(1.81-1.67*DX+2.16*DX**2) | |
112 | POM2=DX**0.7/(1.-0.4*LOG(1.-DX)) | |
113 | POM3=38.D-4*(1.-DX)**1.82/(DX**1.18) | |
114 | DDB=(CQ(1)**2*POM1*POM2+POM3)*FQ | |
115 | DDV=DDB | |
116 | DUB=(CQ(2)**2*POM1*POM2+POM3)*FQ | |
117 | DUV=DUB | |
118 | DSB=(CQ(3)**2*POM1*POM2+POM3)*FQ | |
119 | DCB=(CQ(4)**2*POM1*POM2+POM3)*FQ | |
120 | DBB=(CQ(5)**2*POM1*POM2+POM3)*FQ | |
121 | RETURN | |
122 | END | |
123 | ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc | |
124 | SUBROUTINE DOPHO2(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL) | |
125 | C******************************************************************** | |
126 | C* * | |
127 | C* Parametrization of parton distribution functions * | |
128 | C* in the photon (LO analysis) - asymptotic solution of AP eq.! * | |
129 | C* * | |
130 | C* authors: D.Duke and H.Owens (DO) * | |
131 | C* /Phys.Rev. D26 (1982) 1600/ * | |
132 | C* * | |
133 | C* Prepared by: * | |
134 | C* Krzysztof Charchula, DESY * | |
135 | C* bitnet: F1PCHA@DHHDESY3 * | |
136 | C* decnet: 13313::CHARCHULA * | |
137 | C* * | |
138 | C* Modified by: * | |
139 | C* H. Plothow-Besch/CERN-PPE * | |
140 | C* * | |
141 | C******************************************************************** | |
142 | C | |
143 | implicit real*8 (a-h,o-z) | |
144 | double precision | |
145 | + CQ(5), | |
146 | + DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL | |
147 | PARAMETER (ALPEM=7.29927D-3,PI=3.141592D0) | |
148 | PARAMETER (ALAM=0.4D0) | |
149 | DATA CQ/0.33333D0,0.66666D0,0.33333D0,0.66666D0,0.33333D0/ | |
150 | C | |
151 | Q2 = DQ*DQ | |
152 | ALAM2=ALAM**2 | |
153 | FQ=ALPEM/(2.*PI)*LOG(Q2/ALAM2) | |
154 | C | |
155 | C...gluons | |
156 | POMG=0.194*(1.-DX)**1.03/(DX**0.97) | |
157 | DGL=POMG*FQ | |
158 | C | |
159 | C...quarks | |
160 | POM1=(1.81-1.67*DX+2.16*DX**2) | |
161 | POM2=DX**0.7/(1.-0.4*LOG(1.-DX)) | |
162 | POM3=38.D-4*(1.-DX)**1.82/(DX**1.18) | |
163 | DDB=(CQ(1)**2*POM1*POM2+POM3)*FQ | |
164 | DDV=DDB | |
165 | DUB=(CQ(2)**2*POM1*POM2+POM3)*FQ | |
166 | DUV=DUB | |
167 | DSB=(CQ(3)**2*POM1*POM2+POM3)*FQ | |
168 | DCB=(CQ(4)**2*POM1*POM2+POM3)*FQ | |
169 | DBB=(CQ(5)**2*POM1*POM2+POM3)*FQ | |
170 | RETURN | |
171 | END | |
172 | ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc |