Change needed for G4
[u/mrichter/AliRoot.git] / LHAPDF / lhapdf5.3.1 / wrapgrvg.f
CommitLineData
4e9e3152 1 subroutine GRVGevolvep0(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,zbot,zchm
5 character*16 name(nmxset)
6 integer nmem(nmxset),ndef(nmxset),mmem
7 common/NAME/name,nmem,ndef,mmem
8 integer nset
9
10 save
11
12 call getnset(iset)
13 call getnmem(iset,imem)
14 if(imem.eq.1) then
15 call GRVGALO (xin,qin,upv,dnv,usea,dsea,str,chm,bot,glu)
16 elseif(imem.eq.2.or.imem.eq.0) then
17 q2in = qin*qin
18c calls GRVGALO for charm and bottom, rest from GRSGALO
19 call GRVGALO(xin,qin,upv,dnv,usea,dsea,str,chm,bot,glu)
20 call GRSGALO(xin,q2in,p2in,
21 + upv,dnv,usea,dsea,str,zchm,zbot,glu)
22 else
23 CONTINUE
24 endif
25 pdf(-6)= 0.0d0
26 pdf(6)= 0.0d0
27 pdf(-5)= bot
28 pdf(5 )= bot
29 pdf(-4)= chm
30 pdf(4 )= chm
31 pdf(-3)= str
32 pdf(3 )= str
33 pdf(-2)= usea
34 pdf(2 )= upv
35 pdf(-1)= dsea
36 pdf(1 )= dnv
37 pdf(0 )= glu
38
39 return
40ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
41 entry GRVGevolvep1(xin,qin,p2in,ip2in,pdf)
42
43 if(imem.eq.1) then
44 call GRVGAH0 (xin,qin,upv,dnv,usea,dsea,str,chm,bot,glu)
45 elseif(imem.eq.2 .or. imem.eq.0) then
46 call GRVGAHO (xin,qin,upv,dnv,usea,dsea,str,chm,bot,glu)
47 else
48 CONTINUE
49 endif
50
51 pdf(-6)= 0.0d0
52 pdf(6)= 0.0d0
53 pdf(-5)= bot
54 pdf(5 )= bot
55 pdf(-4)= chm
56 pdf(4 )= chm
57 pdf(-3)= str
58 pdf(3 )= str
59 pdf(-2)= usea
60 pdf(2 )= upv
61 pdf(-1)= dsea
62 pdf(1 )= dnv
63 pdf(0 )= glu
64
65 return
66ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
67 entry GRVGread(nset)
68 read(1,*)nmem(nset),ndef(nset)
69 return
70c
71ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
72 entry GRVGalfa(alfas,qalfa)
73 call getnset(iset)
74 call GetOrderAsM(iset,iord)
75 call Getlam4M(iset,imem,qcdl4)
76 call Getlam5M(iset,imem,qcdl5)
77 call aspdflib(alfas,Qalfa,iord,qcdl5)
78 return
79c
80ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
81 entry GRVGinit(Eorder,Q2fit)
82 return
83c
84ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
85 entry GRVGpdf(mem)
86 call getnset(iset)
87 call setnmem(iset,mem)
88
89c imem = mem
90 return
91c
92 1000 format(5e13.5)
93 end
94c
95 SUBROUTINE GRVGAH0 (ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZGL)
96* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
97* *
98* G R V - P H O T O N - P A R A M E T R I Z A T I O N S *
99* *
100* FOR A DETAILED EXPLANATION SEE : *
101* M. GLUECK, E.REYA, A.VOGT: DO-TH 91/31 *
102* *
103* THE OUTPUT IS ALWAYS 1./ ALPHA(EM) * X * PARTON DENSITY *
104* output modified by HPB to be always X * PARTON DENSITY *
105* *
106* THE PARAMETRIZATIONS ARE FITTED TO THE PARTON DISTRIBUTIONS *
107* FOR Q ** 2 BETWEEN MU ** 2 (= 0.25 / 0.30 GEV ** 2 IN LO *
108* / HO) AND 1.E6 GEV ** 2 AND FOR X BETWEEN 1.E-5 AND 1. *
109* *
110* HEAVY QUARK THRESHOLDS Q(H) = M(H) : *
111* M(C) = 1.5, M(B) = 4.5, M(T) = 100 GEV *
112* *
113* CORRESPONDING LAMBDA(F) VALUES FOR F ACTIVE FLAVOURS : *
114* LO : LAMBDA(3) = 0.232, LAMBDA(4) = 0.200, *
115* LAMBDA(5) = 0.153, LAMBDA(6) = 0.082 GEV *
116* HO : LAMBDA(3) = 0.248, LAMBDA(4) = 0.200, *
117* LAMBDA(5) = 0.131, LAMBDA(6) = 0.053 GEV *
118* *
119* HO DISTRIBUTIONS REFER TO THE DIS(GAMMA) SCHEME, SEE : *
120* M. GLUECK, E.REYA, A.VOGT: DO-TH 91/26 *
121* *
122* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
123C
124 IMPLICIT REAL (A - Y)
125 double precision
126 + ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZGL
127 REAL X, Q
128 DATA ALPHEM/7.29927D-3/
129 X = ZX
130 Q = ZQ
131 MU2 = 0.3
132 LAM2 = 0.248 * 0.248
133 Q2 = Q*Q
134 S = ALOG (ALOG(Q2/LAM2) / ALOG(MU2/LAM2))
135 SS = SQRT (S)
136 S2 = S * S
137C...X * U = X * UBAR :
138 AL = 1.447
139 BE = 0.848
140 AK = 0.527 + 0.200 * S - 0.107 * S2
141 BK = 7.106 - 0.310 * SS - 0.786 * S2
142 AG = 0.197 + 0.533 * S
143 BG = 0.062 - 0.398 * S + 0.109 * S2
144 C = 0.755 * S - 0.112 * S2
145 D = 0.318 - 0.059 * S
146 E = 4.225 + 1.708 * S
147 ES = 1.752 + 0.866 * S
148 U0 = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
149 ZUV = U0 * ALPHEM
150 ZUB = ZUV
151C...X * D = X * DBAR :
152 AL = 1.424
153 BE = 0.770
154 AK = 0.500 + 0.067 * SS - 0.055 * S2
155 BK = 0.376 - 0.453 * SS + 0.405 * S2
156 AG = 0.156 + 0.184 * S
157 BG = 0.0 - 0.528 * S + 0.146 * S2
158 C = 0.121 + 0.092 * S
159 D = 0.379 - 0.301 * S + 0.081 * S2
160 E = 4.346 + 1.638 * S
161 ES = 1.645 + 1.016 * S
162 D0 = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
163 ZDV = D0 * ALPHEM
164 ZDB = ZDV
165C...X * G :
166 AL = 0.661
167 BE = 0.793
168 AK = 0.537 - 0.600 * SS
169 BK = 6.389 - 0.953 * S2
170 AG = 0.558 - 0.383 * SS + 0.261 * S2
171 BG = 0.0 - 0.305 * S
172 C = -0.222 + 0.078 * S2
173 D = 0.153 + 0.978 * S - 0.209 * S2
174 E = 1.429 + 1.772 * S
175 ES = 3.331 + 0.806 * S
176 G0 = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
177 ZGL = G0 * ALPHEM
178C...X * S = X * SBAR :
179 SF = 0.0
180 AL = 1.578
181 BE = 0.863
182 AK = 0.622 + 0.332 * S - 0.300 * S2
183 BK = 2.469
184 AG = 0.211 - 0.064 * SS - 0.018 * S2
185 BG = -0.215 + 0.122 * S
186 C = 0.153
187 D = 0.0 + 0.253 * S - 0.081 * S2
188 E = 3.990 + 2.014 * S
189 ES = 1.720 + 0.986 * S
190 S0 = GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
191 ZSB = S0 * ALPHEM
192C...X * C = X * CBAR :
193 SF = 0.820
194 AL = 0.929
195 BE = 0.381
196 AK = 1.228 - 0.231 * S
197 BK = 3.806 - 0.337 * S2
198 AG = 0.932 + 0.150 * S
199 BG = -0.906
200 C = 1.133
201 D = 0.0 + 0.138 * S - 0.028 * S2
202 E = 5.588 + 0.628 * S
203 ES = 2.665 + 1.054 * S
204 C0 = GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
205 ZCB = C0 * ALPHEM
206C...X * B = X * BBAR :
207 SF = 1.297
208 AL = 0.970
209 BE = 0.207
210 AK = 1.719 - 0.292 * S
211 BK = 0.928 + 0.096 * S
212 AG = 0.845 + 0.178 * S
213 BG = -2.310
214 C = 1.558
215 D = -0.191 + 0.151 * S
216 E = 6.089 + 0.282 * S
217 ES = 3.379 + 1.062 * S
218 B0 = GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
219 ZBB = B0 * ALPHEM
220C
221 RETURN
222 END
223ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
224 SUBROUTINE GRVGAHO (ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZGL)
225* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
226* *
227* G R V - P H O T O N - P A R A M E T R I Z A T I O N S *
228* *
229* FOR A DETAILED EXPLANATION SEE : *
230* M. GLUECK, E.REYA, A.VOGT: DO-TH 91/31 *
231* *
232* THE OUTPUT IS ALWAYS 1./ ALPHA(EM) * X * PARTON DENSITY *
233* output modified by HPB to be always X * PARTON DENSITY *
234* *
235* THE PARAMETRIZATIONS ARE FITTED TO THE PARTON DISTRIBUTIONS *
236* FOR Q ** 2 BETWEEN MU ** 2 (= 0.25 / 0.30 GEV ** 2 IN LO *
237* / HO) AND 1.E6 GEV ** 2 AND FOR X BETWEEN 1.E-5 AND 1. *
238* *
239* HEAVY QUARK THRESHOLDS Q(H) = M(H) : *
240* M(C) = 1.5, M(B) = 4.5, M(T) = 100 GEV *
241* *
242* CORRESPONDING LAMBDA(F) VALUES FOR F ACTIVE FLAVOURS : *
243* LO : LAMBDA(3) = 0.232, LAMBDA(4) = 0.200, *
244* LAMBDA(5) = 0.153, LAMBDA(6) = 0.082 GEV *
245* HO : LAMBDA(3) = 0.248, LAMBDA(4) = 0.200, *
246* LAMBDA(5) = 0.131, LAMBDA(6) = 0.053 GEV *
247* *
248* HO DISTRIBUTIONS REFER TO THE DIS(GAMMA) SCHEME, SEE : *
249* M. GLUECK, E.REYA, A.VOGT: DO-TH 91/26 *
250* *
251* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
252C
253 IMPLICIT REAL (A - Y)
254 double precision
255 + ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZGL
256 DATA ALPHEM/7.29927D-3/
257 REAL X, Q
258 X = ZX
259 Q = ZQ
260 MU2 = 0.3
261 LAM2 = 0.248 * 0.248
262 Q2 = Q*Q
263 S = ALOG (ALOG(Q2/LAM2) / ALOG(MU2/LAM2))
264 SS = SQRT (S)
265 S2 = S * S
266C...X * U = X * UBAR :
267 AL = 0.583
268 BE = 0.688
269 AK = 0.449 - 0.025 * S - 0.071 * S2
270 BK = 5.060 - 1.116 * SS
271 AG = 0.103
272 BG = 0.319 + 0.422 * S
273 C = 1.508 + 4.792 * S - 1.963 * S2
274 D = 1.075 + 0.222 * SS - 0.193 * S2
275 E = 4.147 + 1.131 * S
276 ES = 1.661 + 0.874 * S
277 UH = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
278 ZUV = UH * ALPHEM
279 ZUB = ZUV
280C...X * D = X * DBAR :
281 AL = 0.591
282 BE = 0.698
283 AK = 0.442 - 0.132 * S - 0.058 * S2
284 BK = 5.437 - 1.916 * SS
285 AG = 0.099
286 BG = 0.311 - 0.059 * S
287 C = 0.800 + 0.078 * S - 0.100 * S2
288 D = 0.862 + 0.294 * SS - 0.184 * S2
289 E = 4.202 + 1.352 * S
290 ES = 1.841 + 0.990 * S
291 DH = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
292 ZDV = DH * ALPHEM
293 ZDB = ZDV
294C...X * G :
295 AL = 1.161
296 BE = 1.591
297 AK = 0.530 - 0.742 * SS + 0.025 * S2
298 BK = 5.662
299 AG = 0.533 - 0.281 * SS + 0.218 * S2
300 BG = 0.025 - 0.518 * S + 0.156 * S2
301 C = -0.282 + 0.209 * S2
302 D = 0.107 + 1.058 * S - 0.218 * S2
303 E = 0.0 + 2.704 * S
304 ES = 3.071 - 0.378 * S
305 GH = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
306 ZGL = GH * ALPHEM
307C...X * S = X * SBAR :
308 SF = 0.0
309 AL = 0.635
310 BE = 0.456
311 AK = 1.770 - 0.735 * SS - 0.079 * S2
312 BK = 3.832
313 AG = 0.084 - 0.023 * S
314 BG = 0.136
315 C = 2.119 - 0.942 * S + 0.063 * S2
316 D = 1.271 + 0.076 * S - 0.190 * S2
317 E = 4.604 + 0.737 * S
318 ES = 1.641 + 0.976 * S
319 SH = GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
320 ZSB = SH * ALPHEM
321C...X * C = X * CBAR :
322 SF = 0.820
323 AL = 0.926
324 BE = 0.152
325 AK = 1.142 - 0.175 * S
326 BK = 3.276
327 AG = 0.504 + 0.317 * S
328 BG = -0.433
329 C = 3.334
330 D = 0.398 + 0.326 * S - 0.107 * S2
331 E = 5.493 + 0.408 * S
332 ES = 2.426 + 1.277 * S
333 CH = GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
334 ZCB = CH * ALPHEM
335C...X * B = X * BBAR :
336 SF = 1.297
337 AL = 0.969
338 BE = 0.266
339 AK = 1.953 - 0.391 * S
340 BK = 1.657 - 0.161 * S
341 AG = 1.076 + 0.034 * S
342 BG = -2.015
343 C = 1.662
344 D = 0.353 + 0.016 * S
345 E = 5.713 + 0.249 * S
346 ES = 3.456 + 0.673 * S
347 BH = GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
348 ZBB = BH * ALPHEM
349c
350 RETURN
351 END
352cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
353 SUBROUTINE GRVGALO (ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZGL)
354* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
355* *
356* G R V - P H O T O N - P A R A M E T R I Z A T I O N S *
357* *
358* FOR A DETAILED EXPLANATION SEE : *
359* M. GLUECK, E.REYA, A.VOGT: DO-TH 91/31 *
360* *
361* THE OUTPUT IS ALWAYS 1./ ALPHA(EM) * X * PARTON DENSITY *
362* output modified by HPB to be always X * PARTON DENSITY *
363* *
364* THE PARAMETRIZATIONS ARE FITTED TO THE PARTON DISTRIBUTIONS *
365* FOR Q ** 2 BETWEEN MU ** 2 (= 0.25 / 0.30 GEV ** 2 IN LO *
366* / HO) AND 1.E6 GEV ** 2 AND FOR X BETWEEN 1.E-5 AND 1. *
367* *
368* HEAVY QUARK THRESHOLDS Q(H) = M(H) : *
369* M(C) = 1.5, M(B) = 4.5, M(T) = 100 GEV *
370* *
371* CORRESPONDING LAMBDA(F) VALUES FOR F ACTIVE FLAVOURS : *
372* LO : LAMBDA(3) = 0.232, LAMBDA(4) = 0.200, *
373* LAMBDA(5) = 0.153, LAMBDA(6) = 0.082 GEV *
374* HO : LAMBDA(3) = 0.248, LAMBDA(4) = 0.200, *
375* LAMBDA(5) = 0.131, LAMBDA(6) = 0.053 GEV *
376* *
377* HO DISTRIBUTIONS REFER TO THE DIS(GAMMA) SCHEME, SEE : *
378* M. GLUECK, E.REYA, A.VOGT: DO-TH 91/26 *
379* *
380* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
381C
382 IMPLICIT REAL (A - Y)
383 double precision
384 + ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZGL
385 REAL X, Q
386 DATA ALPHEM/7.29927D-3/
387 X = ZX
388 Q = ZQ
389 MU2 = 0.25
390 LAM2 = 0.232 * 0.232
391 Q2 = Q*Q
392 S = ALOG (ALOG(Q2/LAM2) / ALOG(MU2/LAM2))
393 SS = SQRT (S)
394 S2 = S * S
395C...X * U = X * UBAR :
396 AL = 1.717
397 BE = 0.641
398 AK = 0.500 - 0.176 * S
399 BK = 15.00 - 5.687 * SS - 0.552 * S2
400 AG = 0.235 + 0.046 * SS
401 BG = 0.082 - 0.051 * S + 0.168 * S2
402 C = 0.0 + 0.459 * S
403 D = 0.354 - 0.061 * S
404 E = 4.899 + 1.678 * S
405 ES = 2.046 + 1.389 * S
406 UL = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
407 ZUV = UL * ALPHEM
408 ZUB = ZUV
409C...X * D = X * DBAR :
410 AL = 1.549
411 BE = 0.782
412 AK = 0.496 + 0.026 * S
413 BK = 0.685 - 0.580 * SS + 0.608 * S2
414 AG = 0.233 + 0.302 * S
415 BG = 0.0 - 0.818 * S + 0.198 * S2
416 C = 0.114 + 0.154 * S
417 D = 0.405 - 0.195 * S + 0.046 * S2
418 E = 4.807 + 1.226 * S
419 ES = 2.166 + 0.664 * S
420 DL = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
421 ZDV = DL * ALPHEM
422 ZDB = ZDV
423C...X * G :
424 AL = 0.676
425 BE = 1.089
426 AK = 0.462 - 0.524 * SS
427 BK = 5.451 - 0.804 * S2
428 AG = 0.535 - 0.504 * SS + 0.288 * S2
429 BG = 0.364 - 0.520 * S
430 C = -0.323 + 0.115 * S2
431 D = 0.233 + 0.790 * S - 0.139 * S2
432 E = 0.893 + 1.968 * S
433 ES = 3.432 + 0.392 * S
434 GL = GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
435 ZGL = GL * ALPHEM
436C...X * S = X * SBAR :
437 SF = 0.0
438 AL = 1.609
439 BE = 0.962
440 AK = 0.470 - 0.099 * S2
441 BK = 3.246
442 AG = 0.121 - 0.068 * SS
443 BG = -0.090 + 0.074 * S
444 C = 0.062 + 0.034 * S
445 D = 0.0 + 0.226 * S - 0.060 * S2
446 E = 4.288 + 1.707 * S
447 ES = 2.122 + 0.656 * S
448 SL = GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
449 ZSB = SL * ALPHEM
450C...X * C = X * CBAR :
451 SF = 0.888
452 AL = 0.970
453 BE = 0.545
454 AK = 1.254 - 0.251 * S
455 BK = 3.932 - 0.327 * S2
456 AG = 0.658 + 0.202 * S
457 BG = -0.699
458 C = 0.965
459 D = 0.0 + 0.141 * S - 0.027 * S2
460 E = 4.911 + 0.969 * S
461 ES = 2.796 + 0.952 * S
462 CL = GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
463 ZCB = CL * ALPHEM
464C...X * B = X * BBAR :
465 SF = 1.351
466 AL = 1.016
467 BE = 0.338
468 AK = 1.961 - 0.370 * S
469 BK = 0.923 + 0.119 * S
470 AG = 0.815 + 0.207 * S
471 BG = -2.275
472 C = 1.480
473 D = -0.223 + 0.173 * S
474 E = 5.426 + 0.623 * S
475 ES = 3.819 + 0.901 * S
476 BL = GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
477 ZBB = BL * ALPHEM
478C
479 RETURN
480 END
481cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
482C-----------------------------------------------------------------------
483* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
484* *
485* G R S - LO - VIRTUAL PHOTON PARAMETRIZATIONS *
486* *
487* FOR A DETAILED EXPLANATION SEE *
488* M. GLUECK, E.REYA, M. STRATMANN : *
489* PHYS. REV. D51 (1995) 3220 *
490* *
491* THE PARAMETRIZATIONS ARE FITTED TO THE EVOLVED PARTONS FOR *
492* Q**2 / GEV**2 BETWEEN 0.6 AND 5.E4 *
493* AND (!) Q**2 > 5 P**2 *
494* P**2 / GEV**2 BETWEEN 0.0 AND 10. *
495* P**2 = 0 <=> REAL PHOTON *
496* X BETWEEN 1.E-4 AND 1. *
497* *
498* HEAVY QUARK THRESHOLDS Q(H) = M(H) IN THE BETA FUNCTION : *
499* M(C) = 1.5, M(B) = 4.5 *
500* CORRESPONDING LAMBDA(F) VALUES IN GEV FOR Q**2 > M(H)**2 : *
501* LO : LAMBDA(3) = 0.232, LAMBDA(4) = 0.200, *
502* LAMBDA(5) = 0.153, *
503* THE NUMBER OF ACTIVE QUARK FLAVOURS IS NF = 3 EVERYWHERE *
504* EXCEPT IN THE BETA FUNCTION, I.E. THE HEAVY QUARKS C,B,... *
505* ARE NOT PRESENT AS PARTONS IN THE Q2-EVOLUTION. *
506* *
507* PLEASE REPORT ANY STRANGE BEHAVIOUR TO : *
508* STRAT@HAL1.PHYSIK.UNI-DORTMUND.DE *
509* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
510*
511*...INPUT PARAMETERS :
512*
513* X = MOMENTUM FRACTION
514* Q2 = SCALE Q**2 IN GEV**2
515* P2 = VIRTUALITY OF THE PHOTON IN GEV**2
516*
517*...OUTPUT (ALWAYS X TIMES THE DISTRIBUTION DIVIDED BY ALPHA_EM) :
518*...OUTPUT (ALWAYS X TIMES THE DISTRIBUTION) : modified H.P.-B. 10.9.1996
519*
520********************************************************
521 SUBROUTINE GRSGALO(DX,DQ2,DP2,
522 + DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL)
523C subroutine grsgalo(x,q2,p2,ugam,dgam,sgam,ggam)
524 implicit real*8 (a-h,o-z)
525 double precision
526 + x, q2, p2, mu2, lam2,
527 + ugam, dgam, sgam, ggam,
528 + DUPV,DDNV,DUSEA,DDSEA,DSTR,DCHM,DBOT,DGL
529C
530 dimension u1(40),ds1(40),g1(40)
531 dimension ud2(20),s2(20),g2(20)
532 dimension up0(20),dsp0(20),gp0(20)
533 DATA ALPHEM/7.29927D-3/
534c
535 data u1/-0.139d0,0.783d0,0.132d0,0.087d0,0.003d0,-0.0134d0,
536 + 0.009d0,-0.017d0,0.092d0,-0.516d0,-0.085d0,0.439d0,
537 + 0.013d0,0.108d0,-0.019d0,-0.272d0,-0.167d0,0.138d0,
538 + 0.076d0,0.026d0,-0.013d0,0.27d0,0.107d0,-0.097d0,0.04d0,
539 + 0.064d0,0.011d0,0.002d0,0.057d0,-0.057d0,0.162d0,
540 + -0.172d0,0.124d0,-0.016d0,-0.065d0,0.044d0,-1.009d0,
541 + 0.622d0,0.227d0,-0.184d0/
542 data ds1/0.033d0,0.007d0,-0.0516d0,0.12d0,0.001d0,-0.013d0,
543 + 0.018d0,-0.028d0,0.102d0,-0.595d0,-0.114d0,0.669d0,
544 + 0.022d0,0.001d0,-0.003d0,-0.0583d0,-0.041d0,0.035d0,
545 + 0.009d0,0.009d0,0.004d0,0.054d0,0.025d0,-0.02d0,
546 + 0.007d0,0.021d0,0.01d0,0.004d0,-0.067d0,0.06d0,-0.148d0,
547 + 0.13d0,0.032d0,-0.009d0,-0.06d0,0.036d0,-0.39d0,0.033d0,
548 + 0.245d0,-0.171d0/
549 data g1/0.025d0,0.d0,-0.018d0,0.112d0,-0.025d0,0.177d0,
550 + -0.022d0,0.024d0,0.001d0,-0.0104d0,0.d0,0.d0,-1.082d0,
551 + -1.666d0,0.d0,0.086d0,0.d0,0.053d0,0.005d0,-0.058d0,
552 + 0.034d0,0.073d0,1.08d0,1.63d0,-0.0256d0,-0.088d0,0.d0,
553 + 0.d0,-0.004d0,0.016d0,0.007d0,-0.012d0,0.01d0,-0.673d0,
554 + 0.126d0,-0.167d0,0.032d0,-0.227d0,0.086d0,-0.159d0/
555 data ud2/0.756d0,0.187d0,0.109d0,-0.163d0,0.002d0,0.004d0,
556 + 0.054d0,-0.039d0,22.53d0,-21.02d0,5.608d0,0.332d0,
557 + -0.008d0,-0.021d0,0.381d0,0.572d0,4.774d0,1.436d0,
558 + -0.614d0,3.548d0/
559 data s2/0.902d0,0.182d0,0.271d0,-0.346d0,0.017d0,-0.01d0,
560 + -0.011d0,0.0065d0,17.1d0,-13.29d0,6.519d0,0.031d0,
561 + -0.0176d0,0.003d0,1.243d0,0.804d0,4.709d0,1.499d0,
562 + -0.48d0,3.401d0/
563 data g2/0.364d0,1.31d0,0.86d0,-0.254d0,0.611d0,0.008d0,
564 + -0.097d0,-2.412d0,-0.843d0,2.248d0,-0.201d0,1.33d0,
565 + 0.572d0,0.44d0,1.233d0,0.009d0,0.954d0,1.862d0,3.791d0,
566 + -0.079d0/
567 data up0/1.551d0,0.105d0,1.089d0,-0.172d0,3.822d0,-2.162d0,
568 + 0.533d0,-0.467d0,-0.412d0,0.2d0,0.377d0,0.299d0,0.487d0,
569 + 0.0766d0,0.119d0,0.063d0,7.605d0,0.234d0,-0.567d0,
570 + 2.294d0/
571 data dsp0/2.484d0,1.214d0,1.088d0,-0.1735d0,4.293d0,
572 + -2.802d0,0.5975d0,-0.1193d0,-0.0872d0,0.0418d0,0.128d0,
573 + 0.0337d0,0.127d0,0.0135d0,0.14d0,0.0423d0,6.946d0,
574 + 0.814d0,1.531d0,0.124d0/
575 data gp0/1.682d0,1.1d0,0.5888d0,-0.4714d0,0.5362d0,0.0127d0,
576 + -2.438d0,0.03399d0,0.07825d0,0.05842d0,0.08393d0,2.348d0,
577 + -0.07182d0,1.084d0,0.3098d0,-0.07514d0,3.327d0,1.1d0,
578 + 2.264d0,0.2675d0/
579c
580 save u1,ds1,g1,ud2,s2,g2,up0,dsp0,gp0
581c
582 x = DX
583 q = SQRT(DQ2)
584 q2 = DQ2
585 p2 = DP2
586 mu2=0.25d0
587 lam2=0.232d0*0.232d0
588c
589 if(p2.le.0.25d0) then
590 s=log(log(q2/lam2)/log(mu2/lam2))
591 lp1=0.d0
592 lp2=0.d0
593 else
594 if(q2.lt.p2) then
595 write(*,1000)
596 1000 format
597 + (' WARNING: GRSGALO has been called with Q2 < P2 !',/,
598 + ' GRSGALO is about to blow up, therefore',/,
599 + ' Q2 is set equal to P2')
600 q2=p2
601 endif
602 s=log(log(q2/lam2)/log(p2/lam2))
603 lp1=log(p2/mu2)*log(p2/mu2)
604 lp2=log(p2/mu2+log(p2/mu2))
605 endif
606c
607 alp=up0(1)+lp1*u1(1)+lp2*u1(2)
608 bet=up0(2)+lp1*u1(3)+lp2*u1(4)
609 a=up0(3)+lp1*u1(5)+lp2*u1(6)+
610 + (up0(4)+lp1*u1(7)+lp2*u1(8))*s
611 b=up0(5)+lp1*u1(9)+lp2*u1(10)+
612 + (up0(6)+lp1*u1(11)+lp2*u1(12))*s**0.5+
613 + (up0(7)+lp1*u1(13)+lp2*u1(14))*s**2
614 gb=up0(8)+lp1*u1(15)+lp2*u1(16)+
615 + (up0(9)+lp1*u1(17)+lp2*u1(18))*s+
616 + (up0(10)+lp1*u1(19)+lp2*u1(20))*s**2
617 ga=up0(11)+lp1*u1(21)+lp2*u1(22)+
618 + (up0(12)+lp1*u1(23)+lp2*u1(24))*s**0.5
619 gc=up0(13)+lp1*u1(25)+lp2*u1(33)+
620 + (up0(14)+lp1*u1(26)+lp2*u1(34))*s
621 gd=up0(15)+lp1*u1(27)+lp2*u1(35)+
622 + (up0(16)+lp1*u1(28)+lp2*u1(36))*s
623 ge=up0(17)+lp1*u1(29)+lp2*u1(37)+
624 + (up0(18)+lp1*u1(30)+lp2*u1(38))*s
625 gep=up0(19)+lp1*u1(31)+lp2*u1(39)+
626 + (up0(20)+lp1*u1(32)+lp2*u1(40))*s
627 upart1=grsf2(x,s,alp,bet,a,b,ga,gb,gc,gd,ge,gep)
628c
629 alp=dsp0(1)+lp1*ds1(1)+lp2*ds1(2)
630 bet=dsp0(2)+lp1*ds1(3)+lp2*ds1(4)
631 a=dsp0(3)+lp1*ds1(5)+lp2*ds1(6)+
632 + (dsp0(4)+lp1*ds1(7)+lp2*ds1(8))*s
633 b=dsp0(5)+lp1*ds1(9)+lp2*ds1(10)+
634 + (dsp0(6)+lp1*ds1(11)+lp2*ds1(12))*s**0.5+
635 + (dsp0(7)+lp1*ds1(13)+lp2*ds1(14))*s**2
636 gb=dsp0(8)+lp1*ds1(15)+lp2*ds1(16)+
637 + (dsp0(9)+lp1*ds1(17)+lp2*ds1(18))*s+
638 + (dsp0(10)+lp1*ds1(19)+lp2*ds1(20))*s**2
639 ga=dsp0(11)+lp1*ds1(21)+lp2*ds1(22)+
640 + (dsp0(12)+lp1*ds1(23)+lp2*ds1(24))*s
641 gc=dsp0(13)+lp1*ds1(25)+lp2*ds1(33)+
642 + (dsp0(14)+lp1*ds1(26)+lp2*ds1(34))*s
643 gd=dsp0(15)+lp1*ds1(27)+lp2*ds1(35)+
644 + (dsp0(16)+lp1*ds1(28)+lp2*ds1(36))*s
645 ge=dsp0(17)+lp1*ds1(29)+lp2*ds1(37)+
646 + (dsp0(18)+lp1*ds1(30)+lp2*ds1(38))*s
647 gep=dsp0(19)+lp1*ds1(31)+lp2*ds1(39)+
648 + (dsp0(20)+lp1*ds1(32)+lp2*ds1(40))*s
649 dspart1=grsf2(x,s,alp,bet,a,b,ga,gb,gc,gd,ge,gep)
650c
651 alp=gp0(1)+lp1*g1(1)+lp2*g1(2)
652 bet=gp0(2)+lp1*g1(3)+lp2*g1(4)
653 a=gp0(3)+lp1*g1(5)+lp2*g1(6)+
654 + (gp0(4)+lp1*g1(7)+lp2*g1(8))*s**0.5
655 b=gp0(5)+lp1*g1(9)+lp2*g1(10)+
656 + (gp0(6)+lp1*g1(11)+lp2*g1(12))*s**2
657 gb=gp0(7)+lp1*g1(13)+lp2*g1(14)+
658 + (gp0(8)+lp1*g1(15)+lp2*g1(16))*s
659 ga=gp0(9)+lp1*g1(17)+lp2*g1(18)+
660 + (gp0(10)+lp1*g1(19)+lp2*g1(20))*s**0.5+
661 + (gp0(11)+lp1*g1(21)+lp2*g1(22))*s**2
662 gc=gp0(12)+lp1*g1(23)+lp2*g1(24)+
663 + (gp0(13)+lp1*g1(25)+lp2*g1(26))*s**2
664 gd=gp0(14)+lp1*g1(27)+lp2*g1(28)+
665 + (gp0(15)+lp1*g1(29)+lp2*g1(30))*s+
666 + (gp0(16)+lp1*g1(31)+lp2*g1(32))*s**2
667 ge=gp0(17)+lp1*g1(33)+lp2*g1(34)+
668 + (gp0(18)+lp1*g1(35)+lp2*g1(36))*s
669 gep=gp0(19)+lp1*g1(37)+lp2*g1(38)+
670 + (gp0(20)+lp1*g1(39)+lp2*g1(40))*s
671 gpart1=grsf2(x,s,alp,bet,a,b,ga,gb,gc,gd,ge,gep)
672c
673 s=log(log(q2/lam2)/log(mu2/lam2))
674 suppr=1.d0/(1.d0+p2/0.59d0)**2
675c
676 alp=ud2(1)
677 bet=ud2(2)
678 a=ud2(3)+ud2(4)*s
679 ga=ud2(5)+ud2(6)*s**0.5
680 gc=ud2(7)+ud2(8)*s
681 b=ud2(9)+ud2(10)*s+ud2(11)*s**2
682 gb=ud2(12)+ud2(13)*s+ud2(14)*s**2
683 gd=ud2(15)+ud2(16)*s
684 ge=ud2(17)+ud2(18)*s
685 gep=ud2(19)+ud2(20)*s
686 udpart2=suppr*grsf1(x,s,alp,bet,a,b,ga,gb,gc,gd,ge,gep)
687c
688 alp=s2(1)
689 bet=s2(2)
690 a=s2(3)+s2(4)*s
691 ga=s2(5)+s2(6)*s**0.5
692 gc=s2(7)+s2(8)*s
693 b=s2(9)+s2(10)*s+s2(11)*s**2
694 gb=s2(12)+s2(13)*s+s2(14)*s**2
695 gd=s2(15)+s2(16)*s
696 ge=s2(17)+s2(18)*s
697 gep=s2(19)+s2(20)*s
698 spart2=suppr*grsf2(x,s,alp,bet,a,b,ga,gb,gc,gd,ge,gep)
699c
700 alp=g2(1)
701 bet=g2(2)
702 a=g2(3)+g2(4)*s**0.5
703 b=g2(5)+g2(6)*s**2
704 gb=g2(7)+g2(8)*s
705 ga=g2(9)+g2(10)*s**0.5+g2(11)*s**2
706 gc=g2(12)+g2(13)*s**2
707 gd=g2(14)+g2(15)*s+g2(16)*s**2
708 ge=g2(17)+g2(18)*s
709 gep=g2(19)+g2(20)*s
710 gpart2=suppr*grsf1(x,s,alp,bet,a,b,ga,gb,gc,gd,ge,gep)
711c
712 ugam=upart1+udpart2
713 DUPV = UGAM * ALPHEM
714 DUSEA = DUPV
715 dgam=dspart1+udpart2
716 DDNV = DGAM * ALPHEM
717 DDSEA = DDNV
718 sgam=dspart1+spart2
719 DSTR = SGAM * ALPHEM
720 ggam=gpart1+gpart2
721 DGL = GGAM * ALPHEM
722C
723 DCHM = 0.D0
724 DBOT = 0.D0
725c
726 return
727 end
728ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
729 FUNCTION GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
730 IMPLICIT REAL (A - Z)
731 SX = SQRT (X)
732 LX = ALOG (1./X)
733 GRVGF = (X**AK * (AG + BG * SX + C * X**BK) + S**AL
734 1 * EXP (-E + SQRT (ES * S**BE * LX))) * (1.- X)**D
735 RETURN
736 END
737ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
738 FUNCTION GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
739 IMPLICIT REAL (A - Z)
740 IF (S .LE. SF) THEN
741 GRVGFS = 0.0
742 ELSE
743 SX = SQRT (X)
744 LX = ALOG (1./X)
745 DS = S - SF
746 GRVGFS = (DS * X**AK * (AG + BG * SX + C * X**BK) + DS**AL
747 1 * EXP (-E + SQRT (ES * S**BE * LX))) * (1.- X)**D
748 END IF
749 RETURN
750 END
751ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
752 double precision function grsf1(x,s,alp,bet,a,b,ga,gb,gc,gd,
753 + ge,gep)
754 implicit real*8 (a-h,o-z)
755C
756 grsf1=(x**a*(ga+gb*sqrt(x)+gc*x**b)+
757 + s**alp*exp(-ge+sqrt(gep*s**bet*log(1.d0/x))))*
758 + (1.d0-x)**gd
759 return
760 end
761ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
762 double precision function grsf2(x,s,alp,bet,a,b,ga,gb,gc,gd,
763 + ge,gep)
764 implicit real*8 (a-h,o-z)
765C
766 grsf2=(s*x**a*(ga+gb*sqrt(x)+gc*x**b)+
767 + s**alp*exp(-ge+sqrt(gep*s**bet*log(1.d0/x))))*
768 + (1.d0-x)**gd
769 return
770 end
771ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
772