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
12 call DOPHO1(xin,qin,upv,dnv,usea,dsea,str,chm,bot,glu)
29 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
30 entry DOGevolvep1(xin,qin,p2in,ip2in,pdf)
32 call DOPHO2(xin,qin,upv,dnv,usea,dsea,str,chm,bot,glu)
49 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
51 read(1,*)nmem(nset),ndef(nset)
54 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
55 entry DOGalfa(alfas,qalfa)
57 call GetOrderAsM(iset,iord)
58 call Getlam4M(iset,imem,qcdl4)
59 call Getlam5M(iset,imem,qcdl5)
60 call aspdflib(alfas,Qalfa,iord,qcdl5)
63 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
64 entry DOGinit(Eorder,Q2fit)
67 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
74 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
75 SUBROUTINE DOPHO1(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL)
76 C********************************************************************
78 C* Parametrization of parton distribution functions *
79 C* in the photon (LO analysis) - asymptotic solution of AP eq.! *
81 C* authors: D.Duke and H.Owens (DO) *
82 C* /Phys.Rev. D26 (1982) 1600/ *
85 C* Krzysztof Charchula, DESY *
86 C* bitnet: F1PCHA@DHHDESY3 *
87 C* decnet: 13313::CHARCHULA *
90 C* H. Plothow-Besch/CERN-PPE *
92 C********************************************************************
94 implicit real*8 (a-h,o-z)
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/
104 FQ=ALPEM/(2.*PI)*LOG(Q2/ALAM2)
107 POMG=0.194*(1.-DX)**1.03/(DX**0.97)
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
116 DUB=(CQ(2)**2*POM1*POM2+POM3)*FQ
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
123 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
124 SUBROUTINE DOPHO2(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL)
125 C********************************************************************
127 C* Parametrization of parton distribution functions *
128 C* in the photon (LO analysis) - asymptotic solution of AP eq.! *
130 C* authors: D.Duke and H.Owens (DO) *
131 C* /Phys.Rev. D26 (1982) 1600/ *
134 C* Krzysztof Charchula, DESY *
135 C* bitnet: F1PCHA@DHHDESY3 *
136 C* decnet: 13313::CHARCHULA *
139 C* H. Plothow-Besch/CERN-PPE *
141 C********************************************************************
143 implicit real*8 (a-h,o-z)
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/
153 FQ=ALPEM/(2.*PI)*LOG(Q2/ALAM2)
156 POMG=0.194*(1.-DX)**1.03/(DX**0.97)
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
165 DUB=(CQ(2)**2*POM1*POM2+POM3)*FQ
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
172 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc