1 subroutine DGGevolvep(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 getnmem(iset,imem)
14 if(imem.eq.1.or.imem.eq.0) then
15 call DGPHO1(xin,qin,upv,dnv,usea,dsea,str,chm,bot,glu)
17 elseif(imem.eq.2) then
18 call DGPHO2(xin,qin,upv,dnv,usea,dsea,str,chm,bot,glu)
20 elseif(imem.eq.3) then
21 call DGPHO3(xin,qin,upv,dnv,usea,dsea,str,chm,bot,glu)
23 elseif(imem.eq.4) then
24 call DGPHO4(xin,qin,upv,dnv,usea,dsea,str,chm,bot,glu)
45 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
47 read(1,*)nmem(nset),ndef(nset)
50 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
51 entry DGGalfa(alfas,qalfa)
53 call getnmem(iset,imem)
54 call GetOrderAsM(iset,iord)
55 call Getlam4M(iset,imem,qcdl4)
56 call Getlam5M(iset,imem,qcdl5)
57 call aspdflib(alfas,Qalfa,iord,qcdl5)
60 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
61 entry DGGinit(Eorder,Q2fit)
64 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
67 call setnmem(iset,mem)
73 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
74 SUBROUTINE DGPHO1(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL)
75 C********************************************************************
77 C* Parametrization of parton distribution functions *
78 C* in the photon (LO analysis) - full solution of AP eq.! *
80 C* authors: M.Drees and K.Grassie (DG) *
81 C* /Z. Phys. C28 (1985) 451/ *
84 C* Krzysztof Charchula, DESY *
85 C* bitnet: F1PCHA@DHHDESY3 *
86 C* decnet: 13313::CHARCHULA *
89 C* H. Plothow-Besch/CERN-PPE *
91 C********************************************************************
93 implicit real*8 (a-h,o-z)
96 + B(5,4,2,3),BT(5,2),XQPOM(2),E(2),
97 + DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL
98 PARAMETER (ALPEM=7.29927D-3, PI=3.141592D0)
99 PARAMETER (ALAM=0.4D0)
101 C...--------------------------------------------------
102 C... nf=3 for 1< Q2 <32 GeV2
103 C... nf=4 for 32< Q2 <200 GeV2
104 C... nf=5 for 200< Q2 <1D4 GeV2
105 C...--------------------------------------------------
107 C...initialization of gluon parameters array for DG
108 DATA (((A(I,J,K),I=1,3),J=1,4),K=1,3)/
109 > -0.20700, -0.19870, 5.1190,
110 > 0.61580, 0.62570, -0.2752,
111 > 1.07400, 8.35200, -6.9930,
112 > 0.00000, 5.02400, 2.2980,
113 > 0.8926D-2,0.0509, -0.2313,
114 > 0.65940, 0.27740, 0.1382,
115 > 0.47660, -0.39060, 6.5420,
116 > 0.01975, -0.32120, 0.5162,
117 > 0.03197, -0.618D-2,-0.1216,
118 > 1.01800, 0.94760, 0.9047,
119 > 0.24610, -0.60940, 2.6530,
120 > 0.02707, -0.01067, 0.2003D-2/
122 C...initialization of quark parameters array for DG
123 DATA (((B(I,J,K,1),I=1,5),J=1,4),K=1,2)/
124 > 2.2850, 6.0730, -0.4202, -0.0808, 0.0553,
125 > -0.0153, -0.8132, 0.0178, 0.6346, 1.1360,
126 > 1.33D3, -41.310, 0.9216, 1.2080, 0.9512,
127 > 4.2190, 3.1650, 0.1800, 0.2030, 0.0116,
128 > 16.690, 0.1760, -0.0208, -0.0168, -0.1986,
129 > -0.7916, 0.0479, 0.3386D-2, 1.3530, 1.1000,
130 > 1.0990D3, 1.0470, 4.8530, 1.4260, 1.1360,
131 > 4.4280, 0.0250, 0.8404, 1.2390, -0.2779/
132 DATA (((B(I,J,K,2),I=1,5),J=1,4),K=1,2)/
133 > -0.3711, -0.1717, 0.08766, -0.8915, -0.1816,
134 > 1.0610, 0.7815, 0.02197, 0.2857, 0.5866,
135 > 4.7580, 1.5350, 0.10960, 2.9730, 2.4210,
136 > -0.0150, 0.7067D-2,0.20400, 0.1185, 0.4059,
137 > -0.1207, 25.000, -0.01230, -0.0919, 0.02015,
138 > 1.0710, -1.6480, 1.16200, 0.7912, 0.9869,
139 > 1.9770, -0.01563, 0.48240, 0.6397, -0.07036,
140 > -0.8625D-2,6.4380, -0.01100, 2.3270, 0.01694/
141 DATA (((B(I,J,K,3),I=1,5),J=1,4),K=1,2)/
142 > 15.8, 2.742, 0.02917, -0.0342, -0.02302,
143 > -0.9464, -0.7332, 0.04657, 0.7196, 0.9229,
144 > -0.5, 0.7148, 0.1785, 0.7338, 0.5873,
145 > -0.2118, 3.287, 0.04811, 0.08139,-0.79D-4,
146 > 6.734, 59.88, -0.3226D-2,-0.03321, 0.1059,
147 > -1.008, -2.983, 0.8432, 0.9475, 0.6954,
148 > -0.08594, 4.48, 0.3616, -0.3198, -0.6663,
149 > 0.07625, 0.9686, 0.1383D-2, 0.02132, 0.3683/
151 C...specification of sets
153 IF (Q2.LT.32.0D0) NFL=3
154 IF((Q2.GE.32.0D0).AND.(Q2.LT.200.0D0)) NFL=4
155 IF (Q2.GE.200.0D0) NFL=5
164 AT(I)=A(I,1,LF)*T**A(I,2,LF)+A(I,3,LF)*T**(-A(I,4,LF))
166 POMG=AT(1)*DX**AT(2)*(1.D0-DX)**AT(3)
173 ELSEIF(NFL.EQ.4) THEN
175 ELSEIF(NFL.EQ.5) THEN
180 BTP=B(I,1,J,LF)*T**B(I,2,J,LF)
181 BT(I,J)=BTP+B(I,3,J,LF)*T**(-B(I,4,J,LF))
185 C...singlet & non-singlet combinations
187 POM1=DX*(DX*DX+(1.D0-DX)**2)/(BT(1,J)-BT(2,J)*LOG(1.D0-DX))
188 POM2=BT(3,J)*DX**BT(4,J)*(1.D0-DX)**BT(5,J)
189 XQPOM(J)=E(J)*POM1+POM2
194 DUB=ALPEM*1.D0/6.D0*(XQPOM(2)+9.D0*XQPOM(1))
195 DDB=ALPEM*1.D0/6.D0*(XQPOM(2)-9.D0/2.D0*XQPOM(1))
199 ELSEIF (NFL.EQ.4) THEN
200 DUB=ALPEM*1.D0/8.D0*(XQPOM(2)+6.D0*XQPOM(1))
202 DDB=ALPEM*1.D0/8.D0*(XQPOM(2)-6.D0*XQPOM(1))
205 ELSEIF (NFL.EQ.5) THEN
206 DUB=ALPEM*1.D0/10.D0*(XQPOM(2)+15.D0/2.D0*XQPOM(1))
208 DDB=ALPEM*1.D0/10.D0*(XQPOM(2)-5.D0*XQPOM(1))
217 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
218 SUBROUTINE DGPHO2(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL)
219 C********************************************************************
221 C* Parametrization of parton distribution functions *
222 C* in the photon (LO analysis) - full solution of AP eq.! *
224 C* authors: M.Drees and K.Grassie (DG) *
225 C* /Z. Phys. C28 (1985) 451/ *
228 C* Krzysztof Charchula, DESY *
229 C* bitnet: F1PCHA@DHHDESY3 *
230 C* decnet: 13313::CHARCHULA *
233 C* H. Plothow-Besch/CERN-PPE *
235 C********************************************************************
237 implicit real*8 (a-h,o-z)
240 + B(5,4,2,3),BT(5,2),XQPOM(2),E(2),
241 + DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL
242 PARAMETER (ALPEM=7.29927D-3, PI=3.141592D0)
243 PARAMETER (ALAM=0.4D0)
245 C...--------------------------------------------------
246 C... with nf=3 (valid for 1< Q2 <50 GeV2)
247 C...--------------------------------------------------
249 C...initialization of gluon parameters array for DG
250 DATA (((A(I,J,K),I=1,3),J=1,4),K=1,3)/
251 > -0.20700, -0.19870, 5.1190,
252 > 0.61580, 0.62570, -0.2752,
253 > 1.07400, 8.35200, -6.9930,
254 > 0.00000, 5.02400, 2.2980,
255 > 0.8926D-2,0.0509, -0.2313,
256 > 0.65940, 0.27740, 0.1382,
257 > 0.47660, -0.39060, 6.5420,
258 > 0.01975, -0.32120, 0.5162,
259 > 0.03197, -0.618D-2,-0.1216,
260 > 1.01800, 0.94760, 0.9047,
261 > 0.24610, -0.60940, 2.6530,
262 > 0.02707, -0.01067, 0.2003D-2/
264 C...initialization of quark parameters array for DG
265 DATA (((B(I,J,K,1),I=1,5),J=1,4),K=1,2)/
266 > 2.2850, 6.0730, -0.4202, -0.0808, 0.0553,
267 > -0.0153, -0.8132, 0.0178, 0.6346, 1.1360,
268 > 1.33D3, -41.310, 0.9216, 1.2080, 0.9512,
269 > 4.2190, 3.1650, 0.1800, 0.2030, 0.0116,
270 > 16.690, 0.1760, -0.0208, -0.0168, -0.1986,
271 > -0.7916, 0.0479, 0.3386D-2, 1.3530, 1.1000,
272 > 1.0990D3, 1.0470, 4.8530, 1.4260, 1.1360,
273 > 4.4280, 0.0250, 0.8404, 1.2390, -0.2779/
274 DATA (((B(I,J,K,2),I=1,5),J=1,4),K=1,2)/
275 > -0.3711, -0.1717, 0.08766, -0.8915, -0.1816,
276 > 1.0610, 0.7815, 0.02197, 0.2857, 0.5866,
277 > 4.7580, 1.5350, 0.10960, 2.9730, 2.4210,
278 > -0.0150, 0.7067D-2,0.20400, 0.1185, 0.4059,
279 > -0.1207, 25.000, -0.01230, -0.0919, 0.02015,
280 > 1.0710, -1.6480, 1.16200, 0.7912, 0.9869,
281 > 1.9770, -0.01563, 0.48240, 0.6397, -0.07036,
282 > -0.8625D-2,6.4380, -0.01100, 2.3270, 0.01694/
283 DATA (((B(I,J,K,3),I=1,5),J=1,4),K=1,2)/
284 > 15.8, 2.742, 0.02917, -0.0342, -0.02302,
285 > -0.9464, -0.7332, 0.04657, 0.7196, 0.9229,
286 > -0.5, 0.7148, 0.1785, 0.7338, 0.5873,
287 > -0.2118, 3.287, 0.04811, 0.08139,-0.79D-4,
288 > 6.734, 59.88, -0.3226D-2,-0.03321, 0.1059,
289 > -1.008, -2.983, 0.8432, 0.9475, 0.6954,
290 > -0.08594, 4.48, 0.3616, -0.3198, -0.6663,
291 > 0.07625, 0.9686, 0.1383D-2, 0.02132, 0.3683/
293 C...specification of sets
304 AT(I)=A(I,1,LF)*T**A(I,2,LF)+A(I,3,LF)*T**(-A(I,4,LF))
306 POMG=AT(1)*DX**AT(2)*(1.D0-DX)**AT(3)
314 BTP=B(I,1,J,LF)*T**B(I,2,J,LF)
315 BT(I,J)=BTP+B(I,3,J,LF)*T**(-B(I,4,J,LF))
319 C...singlet & non-singlet combinations
321 POM1=DX*(DX*DX+(1.D0-DX)**2)/(BT(1,J)-BT(2,J)*LOG(1.D0-DX))
322 POM2=BT(3,J)*DX**BT(4,J)*(1.D0-DX)**BT(5,J)
323 XQPOM(J)=E(J)*POM1+POM2
327 DUB=ALPEM*1.D0/6.D0*(XQPOM(2)+9.D0*XQPOM(1))
329 DDB=ALPEM*1.D0/6.D0*(XQPOM(2)-9.D0/2.D0*XQPOM(1))
337 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
338 SUBROUTINE DGPHO3(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL)
339 C********************************************************************
341 C* Parametrization of parton distribution functions *
342 C* in the photon (LO analysis) - full solution of AP eq.! *
344 C* authors: M.Drees and K.Grassie (DG) *
345 C* /Z. Phys. C28 (1985) 451/ *
348 C* Krzysztof Charchula, DESY *
349 C* bitnet: F1PCHA@DHHDESY3 *
350 C* decnet: 13313::CHARCHULA *
353 C* H. Plothow-Besch/CERN-PPE *
355 C********************************************************************
357 implicit real*8 (a-h,o-z)
360 + B(5,4,2,3),BT(5,2),XQPOM(2),E(2),
361 + DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL
362 PARAMETER (ALPEM=7.29927D-3, PI=3.141592D0)
363 PARAMETER (ALAM=0.4D0)
365 C...--------------------------------------------------
366 C... with nf=4 (valid for 20< Q2 <500 GeV2)
367 C...--------------------------------------------------
369 C...initialization of gluon parameters array for DG
370 DATA (((A(I,J,K),I=1,3),J=1,4),K=1,3)/
371 > -0.20700, -0.19870, 5.1190,
372 > 0.61580, 0.62570, -0.2752,
373 > 1.07400, 8.35200, -6.9930,
374 > 0.00000, 5.02400, 2.2980,
375 > 0.8926D-2,0.0509, -0.2313,
376 > 0.65940, 0.27740, 0.1382,
377 > 0.47660, -0.39060, 6.5420,
378 > 0.01975, -0.32120, 0.5162,
379 > 0.03197, -0.618D-2,-0.1216,
380 > 1.01800, 0.94760, 0.9047,
381 > 0.24610, -0.60940, 2.6530,
382 > 0.02707, -0.01067, 0.2003D-2/
384 C...initialization of quark parameters array for DG
385 DATA (((B(I,J,K,1),I=1,5),J=1,4),K=1,2)/
386 > 2.2850, 6.0730, -0.4202, -0.0808, 0.0553,
387 > -0.0153, -0.8132, 0.0178, 0.6346, 1.1360,
388 > 1.33D3, -41.310, 0.9216, 1.2080, 0.9512,
389 > 4.2190, 3.1650, 0.1800, 0.2030, 0.0116,
390 > 16.690, 0.1760, -0.0208, -0.0168, -0.1986,
391 > -0.7916, 0.0479, 0.3386D-2, 1.3530, 1.1000,
392 > 1.0990D3, 1.0470, 4.8530, 1.4260, 1.1360,
393 > 4.4280, 0.0250, 0.8404, 1.2390, -0.2779/
394 DATA (((B(I,J,K,2),I=1,5),J=1,4),K=1,2)/
395 > -0.3711, -0.1717, 0.08766, -0.8915, -0.1816,
396 > 1.0610, 0.7815, 0.02197, 0.2857, 0.5866,
397 > 4.7580, 1.5350, 0.10960, 2.9730, 2.4210,
398 > -0.0150, 0.7067D-2,0.20400, 0.1185, 0.4059,
399 > -0.1207, 25.000, -0.01230, -0.0919, 0.02015,
400 > 1.0710, -1.6480, 1.16200, 0.7912, 0.9869,
401 > 1.9770, -0.01563, 0.48240, 0.6397, -0.07036,
402 > -0.8625D-2,6.4380, -0.01100, 2.3270, 0.01694/
403 DATA (((B(I,J,K,3),I=1,5),J=1,4),K=1,2)/
404 > 15.8, 2.742, 0.02917, -0.0342, -0.02302,
405 > -0.9464, -0.7332, 0.04657, 0.7196, 0.9229,
406 > -0.5, 0.7148, 0.1785, 0.7338, 0.5873,
407 > -0.2118, 3.287, 0.04811, 0.08139,-0.79D-4,
408 > 6.734, 59.88, -0.3226D-2,-0.03321, 0.1059,
409 > -1.008, -2.983, 0.8432, 0.9475, 0.6954,
410 > -0.08594, 4.48, 0.3616, -0.3198, -0.6663,
411 > 0.07625, 0.9686, 0.1383D-2, 0.02132, 0.3683/
413 C...specification of sets
424 AT(I)=A(I,1,LF)*T**A(I,2,LF)+A(I,3,LF)*T**(-A(I,4,LF))
426 POMG=AT(1)*DX**AT(2)*(1.D0-DX)**AT(3)
434 BTP=B(I,1,J,LF)*T**B(I,2,J,LF)
435 BT(I,J)=BTP+B(I,3,J,LF)*T**(-B(I,4,J,LF))
439 C...singlet & non-singlet combinations
441 POM1=DX*(DX*DX+(1.D0-DX)**2)/(BT(1,J)-BT(2,J)*LOG(1.D0-DX))
442 POM2=BT(3,J)*DX**BT(4,J)*(1.D0-DX)**BT(5,J)
443 XQPOM(J)=E(J)*POM1+POM2
447 DUB=ALPEM*1.D0/8.D0*(XQPOM(2)+6.D0*XQPOM(1))
449 DDB=ALPEM*1.D0/8.D0*(XQPOM(2)-6.D0*XQPOM(1))
457 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
458 SUBROUTINE DGPHO4(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL)
459 C********************************************************************
461 C* Parametrization of parton distribution functions *
462 C* in the photon (LO analysis) - full solution of AP eq.! *
464 C* authors: M.Drees and K.Grassie (DG) *
465 C* /Z. Phys. C28 (1985) 451/ *
468 C* Krzysztof Charchula, DESY *
469 C* bitnet: F1PCHA@DHHDESY3 *
470 C* decnet: 13313::CHARCHULA *
473 C* H. Plothow-Besch/CERN-PPE *
475 C********************************************************************
477 implicit real*8 (a-h,o-z)
480 + B(5,4,2,3),BT(5,2),XQPOM(2),E(2),
481 + DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL
482 PARAMETER (ALPEM=7.29927D-3, PI=3.141592D0)
483 PARAMETER (ALAM=0.4D0)
485 C...--------------------------------------------------
486 C... with nf=5 (valid for 200< Q2 <1D4 GeV2)
487 C...--------------------------------------------------
489 C...initialization of gluon parameters array for DG
490 DATA (((A(I,J,K),I=1,3),J=1,4),K=1,3)/
491 > -0.20700, -0.19870, 5.1190,
492 > 0.61580, 0.62570, -0.2752,
493 > 1.07400, 8.35200, -6.9930,
494 > 0.00000, 5.02400, 2.2980,
495 > 0.8926D-2,0.0509, -0.2313,
496 > 0.65940, 0.27740, 0.1382,
497 > 0.47660, -0.39060, 6.5420,
498 > 0.01975, -0.32120, 0.5162,
499 > 0.03197, -0.618D-2,-0.1216,
500 > 1.01800, 0.94760, 0.9047,
501 > 0.24610, -0.60940, 2.6530,
502 > 0.02707, -0.01067, 0.2003D-2/
504 C...initialization of quark parameters array for DG
505 DATA (((B(I,J,K,1),I=1,5),J=1,4),K=1,2)/
506 > 2.2850, 6.0730, -0.4202, -0.0808, 0.0553,
507 > -0.0153, -0.8132, 0.0178, 0.6346, 1.1360,
508 > 1.33D3, -41.310, 0.9216, 1.2080, 0.9512,
509 > 4.2190, 3.1650, 0.1800, 0.2030, 0.0116,
510 > 16.690, 0.1760, -0.0208, -0.0168, -0.1986,
511 > -0.7916, 0.0479, 0.3386D-2, 1.3530, 1.1000,
512 > 1.0990D3, 1.0470, 4.8530, 1.4260, 1.1360,
513 > 4.4280, 0.0250, 0.8404, 1.2390, -0.2779/
514 DATA (((B(I,J,K,2),I=1,5),J=1,4),K=1,2)/
515 > -0.3711, -0.1717, 0.08766, -0.8915, -0.1816,
516 > 1.0610, 0.7815, 0.02197, 0.2857, 0.5866,
517 > 4.7580, 1.5350, 0.10960, 2.9730, 2.4210,
518 > -0.0150, 0.7067D-2,0.20400, 0.1185, 0.4059,
519 > -0.1207, 25.000, -0.01230, -0.0919, 0.02015,
520 > 1.0710, -1.6480, 1.16200, 0.7912, 0.9869,
521 > 1.9770, -0.01563, 0.48240, 0.6397, -0.07036,
522 > -0.8625D-2,6.4380, -0.01100, 2.3270, 0.01694/
523 DATA (((B(I,J,K,3),I=1,5),J=1,4),K=1,2)/
524 > 15.8, 2.742, 0.02917, -0.0342, -0.02302,
525 > -0.9464, -0.7332, 0.04657, 0.7196, 0.9229,
526 > -0.5, 0.7148, 0.1785, 0.7338, 0.5873,
527 > -0.2118, 3.287, 0.04811, 0.08139,-0.79D-4,
528 > 6.734, 59.88, -0.3226D-2,-0.03321, 0.1059,
529 > -1.008, -2.983, 0.8432, 0.9475, 0.6954,
530 > -0.08594, 4.48, 0.3616, -0.3198, -0.6663,
531 > 0.07625, 0.9686, 0.1383D-2, 0.02132, 0.3683/
533 C...specification of sets
544 AT(I)=A(I,1,LF)*T**A(I,2,LF)+A(I,3,LF)*T**(-A(I,4,LF))
546 POMG=AT(1)*DX**AT(2)*(1.D0-DX)**AT(3)
554 BTP=B(I,1,J,LF)*T**B(I,2,J,LF)
555 BT(I,J)=BTP+B(I,3,J,LF)*T**(-B(I,4,J,LF))
559 C...singlet & non-singlet combinations
561 POM1=DX*(DX*DX+(1.D0-DX)**2)/(BT(1,J)-BT(2,J)*LOG(1.D0-DX))
562 POM2=BT(3,J)*DX**BT(4,J)*(1.D0-DX)**BT(5,J)
563 XQPOM(J)=E(J)*POM1+POM2
567 DUB=ALPEM*1.D0/10.D0*(XQPOM(2)+15.D0/2.D0*XQPOM(1))
570 DDB=ALPEM*1.D0/10.D0*(XQPOM(2)-5.D0*XQPOM(1))
578 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc