Reading QA thresholds from external file II
[u/mrichter/AliRoot.git] / LHAPDF / lhapdf5.3.1 / wrapgrvpi.f
CommitLineData
4e9e3152 1 subroutine GRVP1evolve(xin,qin,pdf)
2 include 'parmsetup.inc'
3 real*8 xin,qin,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 nset
9
10 save
11
12 call grvpiho(xin,Qin,upv,dnv,usea,str,chm,bot,top,glu)
13
14 pdf(-6)= top
15 pdf(6)= top
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+usea
24 pdf(-1)= usea
25 pdf(1 )= dnv+usea
26 pdf(0 )= glu
27
28 return
29c
30 entry GRVP0evolve(xin,qin,pdf)
31
32 call grvpilo(xin,Qin,upv,dnv,usea,str,chm,bot,top,glu)
33
34 pdf(-6)= top
35 pdf(6)= top
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+usea
44 pdf(-1)= usea
45 pdf(1 )= dnv+usea
46 pdf(0 )= glu
47
48 return
49c
50ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
51 entry GRVPread(nset)
52 read(1,*)nmem(nset),ndef(nset)
53 return
54c
55ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
56 entry GRVPalfa(alfas,qalfa)
57 call getnset(iset)
58 call getnmem(iset,imem)
59 call GetOrderAsM(iset,iord)
60 call Getlam4M(iset,imem,qcdl4)
61 call Getlam5M(iset,imem,qcdl5)
62 call aspdflib(alfas,Qalfa,iord,qcdl5)
63 return
64c
65ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
66 entry GRVPinit(Eorder,Q2fit)
67 return
68c
69ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
70 entry GRVPpdf(mem)
71c imem = mem
72 call getnset(iset)
73 call setnmem(iset,mem)
74 return
75c
76 1000 format(5e13.5)
77 end
78c
79ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
80*
81* $Id: wrapgrvpi.f 70 2006-02-15 15:45:29Z whalley $
82*
83* $Log$
84* Revision 1.2 2005/10/07 15:15:05 whalley
85* Changes to most files for V5 - multiset initializations
86*
87* Revision 1.1.1.1 2005/05/06 14:54:43 whalley
88* Initial CVS import of the LHAPDF code and data sets
89*
90* Revision 1.1.1.2 1996/10/30 08:28:44 cernlib
91* Version 7.04
92*
93* Revision 1.1.1.1 1996/04/12 15:29:24 plothow
94* Version 7.01
95*
96*
97 SUBROUTINE GRVPIHO (ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL)
98* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
99* *
100* G R V - P I O N - P A R A M E T R I Z A T I O N S *
101* *
102* FOR A DETAILED EXPLANATION SEE : *
103* M. GLUECK, E.REYA, A.VOGT: DO-TH 91/16 *
104* *
105* THE PARAMETRIZATIONS ARE FITTED TO THE PARTON DISTRIBUTIONS *
106* FOR Q ** 2 BETWEEN MU ** 2 (= 0.25 / 0.30 GEV ** 2 IN LO *
107* / HO) AND 1.E8 GEV ** 2 AND FOR X BETWEEN 1.E-5 AND 1. *
108* REGIONS, WHERE THE DISTRIBUTION UNDER CONSIDERATION IS NEG- *
109* LIGIBLE, I.E. BELOW ABOUT 1.E-4, WERE EXCLUDED FROM THE FIT. *
110* *
111* HEAVY QUARK THRESHOLDS Q(H) = M(H) : *
112* M(C) = 1.5, M(B) = 4.5, M(T) = 100 GEV *
113* *
114* CORRESPONDING LAMBDA(F) VALUES FOR F ACTIVE FLAVOURS : *
115* LO : LAMBDA(3) = 0.232, LAMBDA(4) = 0.200, *
116* LAMBDA(5) = 0.153, LAMBDA(6) = 0.082 GEV *
117* HO : LAMBDA(3) = 0.248, LAMBDA(4) = 0.200, *
118* LAMBDA(5) = 0.131, LAMBDA(6) = 0.053 GEV *
119* *
120* HO DISTRIBUTION REFER TO THE MS-BAR SCHEME OF BARDEEN ET AL. *
121* *
122* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
123C
124 IMPLICIT REAL (A - Y)
125 double precision
126 + ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL
127 REAL X, Q
128 X = ZX
129 Q = ZQ
130 MU2 = 0.3
131 LAM2 = 0.248 * 0.248
132 Q2 = Q*Q
133 S = ALOG (ALOG(Q2/LAM2) / ALOG(MU2/LAM2))
134 DS = SQRT (S)
135 S2 = S * S
136C...X * VALENCE :
137 NV = 0.456 + 0.150 * DS + 0.112 * S - 0.019 * S2
138 AKV = 0.505 - 0.033 * S
139 AGV = 0.748 - 0.669 * DS - 0.133 * S
140 DV = 0.365 + 0.197 * DS + 0.394 * S
141 VAP = GRVFVP (X, NV, AKV, AGV, DV)
142 ZUV = VAP
143 ZDV = ZUV
144C...X * GLUON :
145 ALG = 1.096
146 BEG = 1.371
147 AKG = 0.437 - 0.689 * DS
148 BKG = -0.631
149 AGG = 1.324 - 0.441 * DS - 0.130 * S
150 BGG = -0.955 + 0.259 * S
151 CG = 1.075 - 0.302 * S
152 DG = 1.158 + 1.229 * S
153 EG = 0.0 + 2.510 * S
154 ESG = 2.604 + 0.165 * S
155 GLP = GRVFGP(X,S, ALG, BEG, AKG, BKG, AGG, BGG, CG, DG, EG, ESG)
156 ZGL = GLP
157C...X * QBAR (SU(3)-SYMMETRIC SEA) :
158 SL = 0.0
159 ALS = 0.85
160 BES = 0.96
161 AKS = -0.350 + 0.806 * S
162 AGS = -1.663
163 BS = 3.148
164 DS = 2.273 + 1.438 * S
165 EST = 3.214 + 1.545 * S
166 ESS = 1.341 + 1.938 * S
167 QBP = GRVFQBP (X, S, SL, ALS, BES, AKS, AGS, BS, DS, EST, ESS)
168 ZUDB = QBP
169 ZSB = ZUDB
170C...X * CBAR = X * C :
171 SC = 0.820
172 ALC = 0.98
173 BEC = 0.0
174 AKC = 0.0 - 0.457 * S
175 AGC = 0.0
176 BC = -1.00 + 1.40 * S
177 DC = 1.318 + 0.584 * S
178 EC = 4.45 + 1.235 * S
179 ESC = 1.496 + 1.010 * S
180 CBP = GRVFQBP (X, S, SC, ALC, BEC, AKC, AGC, BC, DC, EC, ESC)
181 ZCB = CBP
182C...X * BBAR = X * B :
183 SBO = 1.297
184 ALB = 0.99
185 BEB = 0.0
186 AKB = 0.0 - 0.172 * S
187 AGB = 0.0
188 BBO = 0.0
189 DB = 1.447 + 0.485 * S
190 EB = 4.79 + 1.164 * S
191 ESB = 1.724 + 2.121 * S
192 BBP = GRVFQBP (X, S, SBO, ALB, BEB, AKB, AGB, BBO, DB, EB, ESB)
193 ZBB = BBP
194C...X * TBAR = X * T :
195 TBP = 0.
196 ZTB = TBP
197 RETURN
198 END
199c=================================================================
200*
201* $Id: wrapgrvpi.f 70 2006-02-15 15:45:29Z whalley $
202*
203* $Log$
204* Revision 1.2 2005/10/07 15:15:05 whalley
205* Changes to most files for V5 - multiset initializations
206*
207* Revision 1.1.1.1 2005/05/06 14:54:43 whalley
208* Initial CVS import of the LHAPDF code and data sets
209*
210* Revision 1.1.1.2 1996/10/30 08:28:44 cernlib
211* Version 7.04
212*
213* Revision 1.1.1.1 1996/04/12 15:29:24 plothow
214* Version 7.01
215*
216*
217 SUBROUTINE GRVPILO (ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL)
218* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
219* *
220* G R V - P I O N - P A R A M E T R I Z A T I O N S *
221* *
222* FOR A DETAILED EXPLANATION SEE : *
223* M. GLUECK, E.REYA, A.VOGT: DO-TH 91/16 *
224* *
225* THE PARAMETRIZATIONS ARE FITTED TO THE PARTON DISTRIBUTIONS *
226* FOR Q ** 2 BETWEEN MU ** 2 (= 0.25 / 0.30 GEV ** 2 IN LO *
227* / HO) AND 1.E8 GEV ** 2 AND FOR X BETWEEN 1.E-5 AND 1. *
228* REGIONS, WHERE THE DISTRIBUTION UNDER CONSIDERATION IS NEG- *
229* LIGIBLE, I.E. BELOW ABOUT 1.E-4, WERE EXCLUDED FROM THE FIT. *
230* *
231* HEAVY QUARK THRESHOLDS Q(H) = M(H) : *
232* M(C) = 1.5, M(B) = 4.5, M(T) = 100 GEV *
233* *
234* CORRESPONDING LAMBDA(F) VALUES FOR F ACTIVE FLAVOURS : *
235* LO : LAMBDA(3) = 0.232, LAMBDA(4) = 0.200, *
236* LAMBDA(5) = 0.153, LAMBDA(6) = 0.082 GEV *
237* HO : LAMBDA(3) = 0.248, LAMBDA(4) = 0.200, *
238* LAMBDA(5) = 0.131, LAMBDA(6) = 0.053 GEV *
239* *
240* HO DISTRIBUTION REFER TO THE MS-BAR SCHEME OF BARDEEN ET AL. *
241* *
242* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
243C
244 IMPLICIT REAL (A - Y)
245 double precision
246 + ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL
247 REAL X, Q
248 X = ZX
249 Q = ZQ
250 MU2 = 0.25
251 LAM2 = 0.232 * 0.232
252 Q2 = Q*Q
253 S = ALOG (ALOG(Q2/LAM2) / ALOG(MU2/LAM2))
254 DS = SQRT (S)
255 S2 = S * S
256C...X * VALENCE :
257 NV = 0.519 + 0.180 * S - 0.011 * S2
258 AKV = 0.499 - 0.027 * S
259 AGV = 0.381 - 0.419 * S
260 DV = 0.367 + 0.563 * S
261 VAP = GRVFVP (X, NV, AKV, AGV, DV)
262 ZUV = VAP
263 ZDV = ZUV
264C...X * GLUON :
265 ALG = 0.599
266 BEG = 1.263
267 AKG = 0.482 + 0.341 * DS
268 BKG = 0.0
269 AGG = 0.678 + 0.877 * S - 0.175 * S2
270 BGG = 0.338 - 1.597 * S
271 CG = 0.0 - 0.233 * S + 0.406 * S2
272 DG = 0.390 + 1.053 * S
273 EG = 0.618 + 2.070 * S
274 ESG = 3.676
275 GLP = GRVFGP(X,S, ALG, BEG, AKG, BKG, AGG, BGG, CG, DG, EG, ESG)
276 ZGL = GLP
277C...X * QBAR (SU(3)-SYMMETRIC SEA) :
278 SL = 0.0
279 ALS = 0.55
280 BES = 0.56
281 AKS = 2.538 - 0.763 * S
282 AGS = -0.748
283 BS = 0.313 + 0.935 * S
284 DS = 3.359
285 EST = 4.433 + 1.301 * S
286 ESS = 9.30 - 0.887 * S
287 QBP = GRVFQBP (X, S, SL, ALS, BES, AKS, AGS, BS, DS, EST, ESS)
288 ZUDB = QBP
289 ZSB = ZUDB
290C...X * CBAR = X * C :
291 SC = 0.888
292 ALC = 1.02
293 BEC = 0.39
294 AKC = 0.0
295 AGC = 0.0
296 BC = 1.008
297 DC = 1.208 + 0.771 * S
298 EC = 4.40 + 1.493 * S
299 ESC = 2.032 + 1.901 * S
300 CBP = GRVFQBP (X, S, SC, ALC, BEC, AKC, AGC, BC, DC, EC, ESC)
301 ZCB = CBP
302C...X * BBAR = X * B :
303 SBO = 1.351
304 ALB = 1.03
305 BEB = 0.39
306 AKB = 0.0
307 AGB = 0.0
308 BBO = 0.0
309 DB = 0.697 + 0.855 * S
310 EB = 4.51 + 1.490 * S
311 ESB = 3.056 + 1.694 * S
312 BBP = GRVFQBP (X, S, SBO, ALB, BEB, AKB, AGB, BBO, DB, EB, ESB)
313 ZBB = BBP
314C...X * TBAR = X * T :
315 TBP = 0.
316 ZTB = TBP
317 RETURN
318 END
319c====================================================================
320*
321* $Id: wrapgrvpi.f 70 2006-02-15 15:45:29Z whalley $
322*
323* $Log$
324* Revision 1.2 2005/10/07 15:15:05 whalley
325* Changes to most files for V5 - multiset initializations
326*
327* Revision 1.1.1.1 2005/05/06 14:54:43 whalley
328* Initial CVS import of the LHAPDF code and data sets
329*
330* Revision 1.1.1.2 1996/10/30 08:28:38 cernlib
331* Version 7.04
332*
333* Revision 1.1.1.1 1996/04/12 15:29:23 plothow
334* Version 7.01
335*
336*
337 FUNCTION GRVFVP (X, N, AK, AG, D)
338 IMPLICIT REAL (A - Z)
339 DX = SQRT (X)
340 GRVFVP = N * X**AK * (1.+ AG*DX) * (1.- X)**D
341 RETURN
342 END
343c====================================================================
344*
345* $Id: wrapgrvpi.f 70 2006-02-15 15:45:29Z whalley $
346*
347* $Log$
348* Revision 1.2 2005/10/07 15:15:05 whalley
349* Changes to most files for V5 - multiset initializations
350*
351* Revision 1.1.1.1 2005/05/06 14:54:43 whalley
352* Initial CVS import of the LHAPDF code and data sets
353*
354* Revision 1.1.1.2 1996/10/30 08:28:37 cernlib
355* Version 7.04
356*
357* Revision 1.1.1.1 1996/04/12 15:29:23 plothow
358* Version 7.01
359*
360*
361 FUNCTION GRVFQBP (X, S, ST, AL, BE, AK, AG, B, D, E, ES)
362 IMPLICIT REAL (A - Z)
363 DX = SQRT (X)
364 LX = ALOG (1./X)
365 IF (S .LE. ST) THEN
366 GRVFQBP = 0.0
367 ELSE
368 GRVFQBP = (S-ST)**AL / LX**AK * (1.+ AG*DX + B*X) * (1.- X)**D
369 1 * EXP (-E + SQRT (ES * S**BE * LX))
370 END IF
371 RETURN
372 END
373c====================================================================
374*
375* $Id: wrapgrvpi.f 70 2006-02-15 15:45:29Z whalley $
376*
377* $Log$
378* Revision 1.2 2005/10/07 15:15:05 whalley
379* Changes to most files for V5 - multiset initializations
380*
381* Revision 1.1.1.1 2005/05/06 14:54:43 whalley
382* Initial CVS import of the LHAPDF code and data sets
383*
384* Revision 1.1.1.2 1996/10/30 08:28:36 cernlib
385* Version 7.04
386*
387* Revision 1.1.1.1 1996/04/12 15:29:23 plothow
388* Version 7.01
389*
390*
391 FUNCTION GRVFGP (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
392 IMPLICIT REAL (A - Z)
393 DX = SQRT (X)
394 LX = ALOG (1./X)
395 GRVFGP = (X**AK * (AG + BG*DX + C*X) * LX**BK + S**AL
396 1 * EXP (-E + SQRT (ES * S**BE * LX))) * (1.- X)**D
397 RETURN
398 END
399c====================================================================