PDF/spdf/grvpiho.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.4  2000/09/18 10:02:35  hristov
   6 * Makefile added to PDF8
   7 *
   8 * Revision 1.1.1.2  1996/10/30 08:28:44  cernlib
   9 * Version 7.04
  10 *
  11 * Revision 1.1.1.1  1996/04/12 15:29:24  plothow
  12 * Version 7.01
  13 *
  14 *
  15 #include "pdf/pilot.h"
  16        SUBROUTINE GRVPIHO (ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL)
  17 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
  18 *                                                                 *
  19 *         G R V - P I 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/16             *
  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.E8 GEV ** 2  AND FOR X BETWEEN  1.E-5  AND  1.   *
  27 *   REGIONS, WHERE THE DISTRIBUTION UNDER CONSIDERATION IS NEG-   *
  28 *   LIGIBLE, I.E. BELOW ABOUT 1.E-4, WERE EXCLUDED FROM THE FIT.  *
  29 *                                                                 *
  30 *              HEAVY QUARK THRESHOLDS  Q(H) = M(H) :              *
  31 *         M(C)  =  1.5,  M(B)  =  4.5,  M(T)  =  100  GEV         *
  32 *                                                                 *
  33 *      CORRESPONDING LAMBDA(F) VALUES FOR F ACTIVE FLAVOURS :     *
  34 *      LO :   LAMBDA(3)  =  0.232,   LAMBDA(4)  =  0.200,         *
  35 *             LAMBDA(5)  =  0.153,   LAMBDA(6)  =  0.082  GEV     *
  36 *      HO :   LAMBDA(3)  =  0.248,   LAMBDA(4)  =  0.200,         *
  37 *             LAMBDA(5)  =  0.131,   LAMBDA(6)  =  0.053  GEV     *
  38 *                                                                 *
  39 *   HO DISTRIBUTION REFER TO THE MS-BAR SCHEME OF BARDEEN ET AL.  *
  40 *                                                                 *
  41 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
  42 C
  43        IMPLICIT REAL (A - Y)
  44 #if defined(CERNLIB_DOUBLE)
  45       DOUBLE PRECISION
  46 #endif
  47 #if defined(CERNLIB_SINGLE)
  48       REAL
  49 #endif
  50      +        ZX,ZQ,ZUV,ZDV,ZUDB,ZSB,ZCB,ZBB,ZTB,ZGL
  51        REAL  X, Q
  52        X = ZX
  53        Q = ZQ
  54        MU2  = 0.3
  55        LAM2 = 0.248 * 0.248
  56        Q2 = Q*Q
  57        S  = ALOG (ALOG(Q2/LAM2) / ALOG(MU2/LAM2))
  58        DS = SQRT (S)
  59        S2 = S * S
  60 C...X * VALENCE :
  61        NV  =  0.456 + 0.150 * DS + 0.112 * S - 0.019 * S2
  62        AKV =  0.505 - 0.033 * S
  63        AGV =  0.748 - 0.669 * DS - 0.133 * S
  64        DV  =  0.365 + 0.197 * DS + 0.394 * S
  65        VAP =  GRVFVP (X, NV, AKV, AGV, DV)
  66        ZUV = VAP
  67        ZDV = ZUV
  68 C...X * GLUON :
  69        ALG =  1.096
  70        BEG =  1.371
  71        AKG =  0.437 - 0.689 * DS
  72        BKG = -0.631
  73        AGG =  1.324 - 0.441 * DS - 0.130 * S
  74        BGG = -0.955 + 0.259 * S
  75        CG  =  1.075 - 0.302 * S
  76        DG  =  1.158 + 1.229 * S
  77        EG  =   0.0  + 2.510 * S
  78        ESG =  2.604 + 0.165 * S
  79        GLP =  GRVFGP(X,S, ALG, BEG, AKG, BKG, AGG, BGG, CG, DG, EG, ESG)
  80        ZGL = GLP
  81 C...X * QBAR (SU(3)-SYMMETRIC SEA) :
  82        SL  =   0.0
  83        ALS =   0.85
  84        BES =   0.96
  85        AKS = -0.350 + 0.806 * S
  86        AGS = -1.663
  87        BS  =  3.148
  88        DS  =  2.273 + 1.438 * S
  89        EST =  3.214 + 1.545 * S
  90        ESS =  1.341 + 1.938 * S
  91        QBP =  GRVFQBP (X, S, SL, ALS, BES, AKS, AGS, BS, DS, EST, ESS)
  92        ZUDB = QBP
  93        ZSB = ZUDB
  94 C...X * CBAR = X * C :
  95        SC  =  0.820
  96        ALC =   0.98
  97        BEC =   0.0
  98        AKC =   0.0  - 0.457 * S
  99        AGC =   0.0
 100        BC  =  -1.00 +  1.40 * S
 101        DC  =  1.318 + 0.584 * S
 102        EC  =   4.45 + 1.235 * S
 103        ESC =  1.496 + 1.010 * S
 104        CBP =  GRVFQBP (X, S, SC, ALC, BEC, AKC, AGC, BC, DC, EC, ESC)
 105        ZCB = CBP
 106 C...X * BBAR = X * B :
 107        SBO =  1.297
 108        ALB =   0.99
 109        BEB =   0.0
 110        AKB =   0.0  - 0.172 * S
 111        AGB =   0.0
 112        BBO =   0.0
 113        DB  =  1.447 + 0.485 * S
 114        EB  =   4.79 + 1.164 * S
 115        ESB =  1.724 + 2.121 * S
 116        BBP =  GRVFQBP (X, S, SBO, ALB, BEB, AKB, AGB, BBO, DB, EB, ESB)
 117        ZBB = BBP
 118 C...X * TBAR = X * T :
 119        TBP = 0.
 120        ZTB = TBP
 121        RETURN
 122        END