| 3c5d1739 |
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==================================================================== |
git://git.uio.no / u / mrichter / AliRoot.git / blame