PDF/spdf/grvgaho.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:28:40  cernlib
   6 * Version 7.04
   7 *
   8 * Revision 1.1.1.1  1996/04/12 15:29:23  plothow
   9 * Version 7.01
  10 *
  11 *
  12 #include "pdf/pilot.h"
  13        SUBROUTINE GRVGAHO (ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZGL)
  14 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
  15 *                                                                 *
  16 *      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      *
  17 *                                                                 *
  18 *                 FOR A DETAILED EXPLANATION SEE :                *
  19 *              M. GLUECK, E.REYA, A.VOGT: DO-TH 91/31             *
  20 *                                                                 *
  21 *    THE OUTPUT IS ALWAYS   1./ ALPHA(EM) * X * PARTON DENSITY    *
  22 *    output modified by HPB to be always    X * PARTON DENSITY    *
  23 *                                                                 *
  24 *   THE PARAMETRIZATIONS ARE FITTED TO THE PARTON DISTRIBUTIONS   *
  25 *   FOR Q ** 2 BETWEEN MU ** 2 (=  0.25 / 0.30  GEV ** 2  IN LO   *
  26 *   / HO) AND  1.E6 GEV ** 2  AND FOR X BETWEEN  1.E-5  AND  1.   *
  27 *                                                                 *
  28 *              HEAVY QUARK THRESHOLDS  Q(H) = M(H) :              *
  29 *         M(C)  =  1.5,  M(B)  =  4.5,  M(T)  =  100  GEV         *
  30 *                                                                 *
  31 *      CORRESPONDING LAMBDA(F) VALUES FOR F ACTIVE FLAVOURS :     *
  32 *      LO :   LAMBDA(3)  =  0.232,   LAMBDA(4)  =  0.200,         *
  33 *             LAMBDA(5)  =  0.153,   LAMBDA(6)  =  0.082  GEV     *
  34 *      HO :   LAMBDA(3)  =  0.248,   LAMBDA(4)  =  0.200,         *
  35 *             LAMBDA(5)  =  0.131,   LAMBDA(6)  =  0.053  GEV     *
  36 *                                                                 *
  37 *      HO DISTRIBUTIONS REFER TO THE DIS(GAMMA) SCHEME, SEE :     *
  38 *              M. GLUECK, E.REYA, A.VOGT: DO-TH 91/26             *
  39 *                                                                 *
  40 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
  41 C
  42        IMPLICIT REAL (A - Y)
  43 #include "pdf/expdp.inc"
  44      +        ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZGL
  45        DATA ALPHEM/7.29927D-3/
  46        REAL  X, Q
  47        X = ZX
  48        Q = ZQ
  49        MU2  = 0.3
  50        LAM2 = 0.248 * 0.248
  51        Q2 = Q*Q
  52        S  = ALOG (ALOG(Q2/LAM2) / ALOG(MU2/LAM2))
  53        SS = SQRT (S)
  54        S2 = S * S
  55 C...X * U = X * UBAR :
  56        AL =  0.583
  57        BE =  0.688
  58        AK =  0.449 - 0.025 * S  - 0.071 * S2
  59        BK =  5.060 - 1.116 * SS
  60        AG =  0.103
  61        BG =  0.319 + 0.422 * S
  62        C  =  1.508 + 4.792 * S  - 1.963 * S2
  63        D  =  1.075 + 0.222 * SS - 0.193 * S2
  64        E  =  4.147 + 1.131 * S
  65        ES =  1.661 + 0.874 * S
  66        UH =  GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
  67        ZUV = UH * ALPHEM
  68        ZUB = ZUV
  69 C...X * D = X * DBAR :
  70        AL =  0.591
  71        BE =  0.698
  72        AK =  0.442 - 0.132 * S  - 0.058 * S2
  73        BK =  5.437 - 1.916 * SS
  74        AG =  0.099
  75        BG =  0.311 - 0.059 * S
  76        C  =  0.800 + 0.078 * S  - 0.100 * S2
  77        D  =  0.862 + 0.294 * SS - 0.184 * S2
  78        E  =  4.202 + 1.352 * S
  79        ES =  1.841 + 0.990 * S
  80        DH  =  GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
  81        ZDV = DH * ALPHEM
  82        ZDB = ZDV
  83 C...X * G :
  84        AL =  1.161
  85        BE =  1.591
  86        AK =  0.530 - 0.742 * SS + 0.025 * S2
  87        BK =  5.662
  88        AG =  0.533 - 0.281 * SS + 0.218 * S2
  89        BG =  0.025 - 0.518 * S  + 0.156 * S2
  90        C  = -0.282              + 0.209 * S2
  91        D  =  0.107 + 1.058 * S  - 0.218 * S2
  92        E  =   0.0  + 2.704 * S
  93        ES =  3.071 - 0.378 * S
  94        GH =  GRVGF (X, S, AL, BE, AK, BK, AG, BG, C, D, E, ES)
  95        ZGL = GH * ALPHEM
  96 C...X * S = X * SBAR :
  97        SF =   0.0
  98        AL =  0.635
  99        BE =  0.456
 100        AK =  1.770 - 0.735 * SS - 0.079 * S2
 101        BK =  3.832
 102        AG =  0.084 - 0.023 * S
 103        BG =  0.136
 104        C  =  2.119 - 0.942 * S  + 0.063 * S2
 105        D  =  1.271 + 0.076 * S  - 0.190 * S2
 106        E  =  4.604 + 0.737 * S
 107        ES =  1.641 + 0.976 * S
 108        SH =  GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
 109        ZSB = SH * ALPHEM
 110 C...X * C = X * CBAR :
 111        SF =  0.820
 112        AL =  0.926
 113        BE =  0.152
 114        AK =  1.142 - 0.175 * S
 115        BK =  3.276
 116        AG =  0.504 + 0.317 * S
 117        BG = -0.433
 118        C  =  3.334
 119        D  =  0.398 + 0.326 * S  - 0.107 * S2
 120        E  =  5.493 + 0.408 * S
 121        ES =  2.426 + 1.277 * S
 122        CH =  GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
 123        ZCB = CH * ALPHEM
 124 C...X * B = X * BBAR :
 125        SF =  1.297
 126        AL =  0.969
 127        BE =  0.266
 128        AK =  1.953 - 0.391 * S
 129        BK =  1.657 - 0.161 * S
 130        AG =  1.076 + 0.034 * S
 131        BG = -2.015
 132        C  =  1.662
 133        D  =  0.353 + 0.016 * S
 134        E  =  5.713 + 0.249 * S
 135        ES =  3.456 + 0.673 * S
 136        BH =  GRVGFS (X, S, SF, AL, BE, AK, BK, AG, BG, C, D, E, ES)
 137        ZBB = BH * ALPHEM
 138 c
 139        RETURN
 140        END