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
12 call grvpiho(xin,Qin,upv,dnv,usea,str,chm,bot,top,glu)
30 entry GRVP0evolve(xin,qin,pdf)
32 call grvpilo(xin,Qin,upv,dnv,usea,str,chm,bot,top,glu)
50 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
52 read(1,*)nmem(nset),ndef(nset)
55 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
56 entry GRVPalfa(alfas,qalfa)
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)
65 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
66 entry GRVPinit(Eorder,Q2fit)
69 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
73 call setnmem(iset,mem)
79 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
81 * $Id: wrapgrvpi.f 70 2006-02-15 15:45:29Z whalley $
84 * Revision 1.2 2005/10/07 15:15:05 whalley
85 * Changes to most files for V5 - multiset initializations
87 * Revision 1.1.1.1 2005/05/06 14:54:43 whalley
88 * Initial CVS import of the LHAPDF code and data sets
90 * Revision 1.1.1.2 1996/10/30 08:28:44 cernlib
93 * Revision 1.1.1.1 1996/04/12 15:29:24 plothow
97 SUBROUTINE GRVPIHO (ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL)
98 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
100 * G R V - P I O N - P A R A M E T R I Z A T I O N S *
102 * FOR A DETAILED EXPLANATION SEE : *
103 * M. GLUECK, E.REYA, A.VOGT: DO-TH 91/16 *
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. *
111 * HEAVY QUARK THRESHOLDS Q(H) = M(H) : *
112 * M(C) = 1.5, M(B) = 4.5, M(T) = 100 GEV *
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 *
120 * HO DISTRIBUTION REFER TO THE MS-BAR SCHEME OF BARDEEN ET AL. *
122 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
124 IMPLICIT REAL (A - Y)
126 + ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL
133 S = ALOG (ALOG(Q2/LAM2) / ALOG(MU2/LAM2))
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)
147 AKG = 0.437 - 0.689 * DS
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
154 ESG = 2.604 + 0.165 * S
155 GLP = GRVFGP(X,S, ALG, BEG, AKG, BKG, AGG, BGG, CG, DG, EG, ESG)
157 C...X * QBAR (SU(3)-SYMMETRIC SEA) :
161 AKS = -0.350 + 0.806 * S
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)
170 C...X * CBAR = X * C :
174 AKC = 0.0 - 0.457 * S
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)
182 C...X * BBAR = X * B :
186 AKB = 0.0 - 0.172 * S
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)
194 C...X * TBAR = X * T :
199 c=================================================================
201 * $Id: wrapgrvpi.f 70 2006-02-15 15:45:29Z whalley $
204 * Revision 1.2 2005/10/07 15:15:05 whalley
205 * Changes to most files for V5 - multiset initializations
207 * Revision 1.1.1.1 2005/05/06 14:54:43 whalley
208 * Initial CVS import of the LHAPDF code and data sets
210 * Revision 1.1.1.2 1996/10/30 08:28:44 cernlib
213 * Revision 1.1.1.1 1996/04/12 15:29:24 plothow
217 SUBROUTINE GRVPILO (ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL)
218 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
220 * G R V - P I O N - P A R A M E T R I Z A T I O N S *
222 * FOR A DETAILED EXPLANATION SEE : *
223 * M. GLUECK, E.REYA, A.VOGT: DO-TH 91/16 *
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. *
231 * HEAVY QUARK THRESHOLDS Q(H) = M(H) : *
232 * M(C) = 1.5, M(B) = 4.5, M(T) = 100 GEV *
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 *
240 * HO DISTRIBUTION REFER TO THE MS-BAR SCHEME OF BARDEEN ET AL. *
242 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
244 IMPLICIT REAL (A - Y)
246 + ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL
253 S = ALOG (ALOG(Q2/LAM2) / ALOG(MU2/LAM2))
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)
267 AKG = 0.482 + 0.341 * DS
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
275 GLP = GRVFGP(X,S, ALG, BEG, AKG, BKG, AGG, BGG, CG, DG, EG, ESG)
277 C...X * QBAR (SU(3)-SYMMETRIC SEA) :
281 AKS = 2.538 - 0.763 * S
283 BS = 0.313 + 0.935 * S
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)
290 C...X * CBAR = X * C :
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)
302 C...X * BBAR = X * B :
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)
314 C...X * TBAR = X * T :
319 c====================================================================
321 * $Id: wrapgrvpi.f 70 2006-02-15 15:45:29Z whalley $
324 * Revision 1.2 2005/10/07 15:15:05 whalley
325 * Changes to most files for V5 - multiset initializations
327 * Revision 1.1.1.1 2005/05/06 14:54:43 whalley
328 * Initial CVS import of the LHAPDF code and data sets
330 * Revision 1.1.1.2 1996/10/30 08:28:38 cernlib
333 * Revision 1.1.1.1 1996/04/12 15:29:23 plothow
337 FUNCTION GRVFVP (X, N, AK, AG, D)
338 IMPLICIT REAL (A - Z)
340 GRVFVP = N * X**AK * (1.+ AG*DX) * (1.- X)**D
343 c====================================================================
345 * $Id: wrapgrvpi.f 70 2006-02-15 15:45:29Z whalley $
348 * Revision 1.2 2005/10/07 15:15:05 whalley
349 * Changes to most files for V5 - multiset initializations
351 * Revision 1.1.1.1 2005/05/06 14:54:43 whalley
352 * Initial CVS import of the LHAPDF code and data sets
354 * Revision 1.1.1.2 1996/10/30 08:28:37 cernlib
357 * Revision 1.1.1.1 1996/04/12 15:29:23 plothow
361 FUNCTION GRVFQBP (X, S, ST, AL, BE, AK, AG, B, D, E, ES)
362 IMPLICIT REAL (A - Z)
368 GRVFQBP = (S-ST)**AL / LX**AK * (1.+ AG*DX + B*X) * (1.- X)**D
369 1 * EXP (-E + SQRT (ES * S**BE * LX))
373 c====================================================================
375 * $Id: wrapgrvpi.f 70 2006-02-15 15:45:29Z whalley $
378 * Revision 1.2 2005/10/07 15:15:05 whalley
379 * Changes to most files for V5 - multiset initializations
381 * Revision 1.1.1.1 2005/05/06 14:54:43 whalley
382 * Initial CVS import of the LHAPDF code and data sets
384 * Revision 1.1.1.2 1996/10/30 08:28:36 cernlib
387 * Revision 1.1.1.1 1996/04/12 15:29:23 plothow
391 FUNCTION GRVFGP (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
392 IMPLICIT REAL (A - Z)
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
399 c====================================================================