LHAPDF 5.2.2 source code.
[u/mrichter/AliRoot.git] / LHAPDF / lhapdf5.2.2 / wrapdog.f
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