LHAPDF/lhapdf5.3.1/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: 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 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
 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: 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 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
 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: 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
 343 c====================================================================
 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
 373 c====================================================================
 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
 399 c====================================================================