LHAPDF 5.2.2 source code.
[u/mrichter/AliRoot.git] / LHAPDF / lhapdf5.2.2 / wrapgrvpi.f
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
29 c      
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
49 c      
50 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
51       entry GRVPread(nset)
52       read(1,*)nmem(nset),ndef(nset)
53       return
54 c
55 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
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
64 c
65 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
66       entry GRVPinit(Eorder,Q2fit)
67       return
68 c
69 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
70       entry GRVPpdf(mem)
71 c      imem = mem
72       call getnset(iset)
73       call setnmem(iset,mem)
74       return
75 c
76  1000 format(5e13.5)
77       end
78 c
79 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
80 *
81 * $Id$
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 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
123 C
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
136 C...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
144 C...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
157 C...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
170 C...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
182 C...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
194 C...X * TBAR = X * T :
195        TBP = 0.
196        ZTB = TBP
197        RETURN
198        END
199 c=================================================================
200 *
201 * $Id$
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 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
243 C
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
256 C...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
264 C...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
277 C...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
290 C...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
302 C...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
314 C...X * TBAR = X * T :
315        TBP = 0.
316        ZTB = TBP
317        RETURN
318        END
319 c====================================================================
320 *
321 * $Id$
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
343 c====================================================================
344 *
345 * $Id$
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
373 c====================================================================
374 *
375 * $Id$
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
399 c====================================================================