PDF/spdf/grvgah0.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:32:35  cernlib
   6 * Version 7.04
   7 *
   8 * Revision 1.2  1996/10/28 15:02:12  pdflib
   9 * PDFLIB new version 7.04
  10 *
  11 * Revision 1.1.1.1  1996/04/12 15:29:23  plothow
  12 * Version 7.01
  13 *
  14 *
  15 #include "pdf/pilot.h"
  16        SUBROUTINE GRVGAH0 (ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZGL)
  17 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
  18 *                                                                 *
  19 *      G R V - P H O T O N - P A R A M E T R I Z A T I O N S      *
  20 *                                                                 *
  21 *                 FOR A DETAILED EXPLANATION SEE :                *
  22 *              M. GLUECK, E.REYA, A.VOGT: DO-TH 91/31             *
  23 *                                                                 *
  24 *    THE OUTPUT IS ALWAYS   1./ ALPHA(EM) * X * PARTON DENSITY    *
  25 *    output modified by HPB to be always    X * PARTON DENSITY    *
  26 *                                                                 *
  27 *   THE PARAMETRIZATIONS ARE FITTED TO THE PARTON DISTRIBUTIONS   *
  28 *   FOR Q ** 2 BETWEEN MU ** 2 (=  0.25 / 0.30  GEV ** 2  IN LO   *
  29 *   / HO) AND  1.E6 GEV ** 2  AND FOR X BETWEEN  1.E-5  AND  1.   *
  30 *                                                                 *
  31 *              HEAVY QUARK THRESHOLDS  Q(H) = M(H) :              *
  32 *         M(C)  =  1.5,  M(B)  =  4.5,  M(T)  =  100  GEV         *
  33 *                                                                 *
  34 *      CORRESPONDING LAMBDA(F) VALUES FOR F ACTIVE FLAVOURS :     *
  35 *      LO :   LAMBDA(3)  =  0.232,   LAMBDA(4)  =  0.200,         *
  36 *             LAMBDA(5)  =  0.153,   LAMBDA(6)  =  0.082  GEV     *
  37 *      HO :   LAMBDA(3)  =  0.248,   LAMBDA(4)  =  0.200,         *
  38 *             LAMBDA(5)  =  0.131,   LAMBDA(6)  =  0.053  GEV     *
  39 *                                                                 *
  40 *      HO DISTRIBUTIONS REFER TO THE DIS(GAMMA) SCHEME, SEE :     *
  41 *              M. GLUECK, E.REYA, A.VOGT: DO-TH 91/26             *
  42 *                                                                 *
  43 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
  44 C
  45        IMPLICIT REAL (A - Y)
  46 #include "pdf/expdp.inc"
  47      +        ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZGL
  48        REAL  X, Q
  49        DATA ALPHEM/7.29927D-3/
  50        X = ZX
  51        Q = ZQ
  52        MU2  = 0.3
  53        LAM2 = 0.248 * 0.248
  54        Q2 = Q*Q
  55        S  = ALOG (ALOG(Q2/LAM2) / ALOG(MU2/LAM2))
  56        SS = SQRT (S)
  57        S2 = S * S
  58 C...X * U = X * UBAR :
  59        AL =  1.447
  60        BE =  0.848
  61        AK =  0.527 + 0.200 * S  - 0.107 * S2
  62        BK =  7.106 - 0.310 * SS - 0.786 * S2
  63        AG =  0.197 + 0.533 * S
  64        BG =  0.062 - 0.398 * S  + 0.109 * S2
  65        C  =          0.755 * S  - 0.112 * S2
  66        D  =  0.318 - 0.059 * S
  67        E  =  4.225 + 1.708 * S
  68        ES =  1.752 + 0.866 * S
  69        U0 =  GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
  70        ZUV = U0 * ALPHEM
  71        ZUB = ZUV
  72 C...X * D = X * DBAR :
  73        AL =  1.424
  74        BE =  0.770
  75        AK =  0.500 + 0.067 * SS - 0.055 * S2
  76        BK =  0.376 - 0.453 * SS + 0.405 * S2
  77        AG =  0.156 + 0.184 * S
  78        BG =   0.0  - 0.528 * S  + 0.146 * S2
  79        C  =  0.121 + 0.092 * S
  80        D  =  0.379 - 0.301 * S  + 0.081 * S2
  81        E  =  4.346 + 1.638 * S
  82        ES =  1.645 + 1.016 * S
  83        D0  =  GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
  84        ZDV = D0 * ALPHEM
  85        ZDB = ZDV
  86 C...X * G :
  87        AL =  0.661
  88        BE =  0.793
  89        AK =  0.537 - 0.600 * SS
  90        BK =  6.389              - 0.953 * S2
  91        AG =  0.558 - 0.383 * SS + 0.261 * S2
  92        BG =   0.0  - 0.305 * S
  93        C  = -0.222              + 0.078 * S2
  94        D  =  0.153 + 0.978 * S  - 0.209 * S2
  95        E  =  1.429 + 1.772 * S
  96        ES =  3.331 + 0.806 * S
  97        G0 =  GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
  98        ZGL = G0 * ALPHEM
  99 C...X * S = X * SBAR :
 100        SF =   0.0
 101        AL =  1.578
 102        BE =  0.863
 103        AK =  0.622 + 0.332 * S  - 0.300 * S2
 104        BK =  2.469
 105        AG =  0.211 - 0.064 * SS - 0.018 * S2
 106        BG = -0.215 + 0.122 * S
 107        C  =  0.153
 108        D  =   0.0  + 0.253 * S  - 0.081 * S2
 109        E  =  3.990 + 2.014 * S
 110        ES =  1.720 + 0.986 * S
 111        S0 =  GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
 112        ZSB = S0 * ALPHEM
 113 C...X * C = X * CBAR :
 114        SF =  0.820
 115        AL =  0.929
 116        BE =  0.381
 117        AK =  1.228 - 0.231 * S
 118        BK =  3.806             - 0.337 * S2
 119        AG =  0.932 + 0.150 * S
 120        BG = -0.906
 121        C  =  1.133
 122        D  =   0.0  + 0.138 * S  - 0.028 * S2
 123        E  =  5.588 + 0.628 * S
 124        ES =  2.665 + 1.054 * S
 125        C0 =  GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
 126        ZCB = C0 * ALPHEM
 127 C...X * B = X * BBAR :
 128        SF =  1.297
 129        AL =  0.970
 130        BE =  0.207
 131        AK =  1.719 - 0.292 * S
 132        BK =  0.928 + 0.096 * S
 133        AG =  0.845 + 0.178 * S
 134        BG = -2.310
 135        C  =  1.558
 136        D  = -0.191 + 0.151 * S
 137        E  =  6.089 + 0.282 * S
 138        ES =  3.379 + 1.062 * S
 139        B0 =  GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
 140        ZBB = B0 * ALPHEM
 141 C
 142        RETURN
 143        END