Change needed for G4
[u/mrichter/AliRoot.git] / LHAPDF / lhapdf5.3.1 / wrapdog.f
CommitLineData
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
29ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
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
49ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
50 entry DOGread(nset)
51 read(1,*)nmem(nset),ndef(nset)
52 return
53c
54ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
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
62c
63ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
64 entry DOGinit(Eorder,Q2fit)
65 return
66c
67ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
68 entry DOGpdf(mem)
69 imem = mem
70 return
71c
72 1000 format(5e13.5)
73 end
74ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
75 SUBROUTINE DOPHO1(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL)
76C********************************************************************
77C* *
78C* Parametrization of parton distribution functions *
79C* in the photon (LO analysis) - asymptotic solution of AP eq.! *
80C* *
81C* authors: D.Duke and H.Owens (DO) *
82C* /Phys.Rev. D26 (1982) 1600/ *
83C* *
84C* Prepared by: *
85C* Krzysztof Charchula, DESY *
86C* bitnet: F1PCHA@DHHDESY3 *
87C* decnet: 13313::CHARCHULA *
88C* *
89C* Modified by: *
90C* H. Plothow-Besch/CERN-PPE *
91C* *
92C********************************************************************
93C
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/
101C
102 Q2 = DQ*DQ
103 ALAM2=ALAM**2
104 FQ=ALPEM/(2.*PI)*LOG(Q2/ALAM2)
105C
106C...gluons
107 POMG=0.194*(1.-DX)**1.03/(DX**0.97)
108 DGL=POMG*FQ
109C
110C...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
123ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
124 SUBROUTINE DOPHO2(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL)
125C********************************************************************
126C* *
127C* Parametrization of parton distribution functions *
128C* in the photon (LO analysis) - asymptotic solution of AP eq.! *
129C* *
130C* authors: D.Duke and H.Owens (DO) *
131C* /Phys.Rev. D26 (1982) 1600/ *
132C* *
133C* Prepared by: *
134C* Krzysztof Charchula, DESY *
135C* bitnet: F1PCHA@DHHDESY3 *
136C* decnet: 13313::CHARCHULA *
137C* *
138C* Modified by: *
139C* H. Plothow-Besch/CERN-PPE *
140C* *
141C********************************************************************
142C
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/
150C
151 Q2 = DQ*DQ
152 ALAM2=ALAM**2
153 FQ=ALPEM/(2.*PI)*LOG(Q2/ALAM2)
154C
155C...gluons
156 POMG=0.194*(1.-DX)**1.03/(DX**0.97)
157 DGL=POMG*FQ
158C
159C...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
172ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc