PDF/spdf/grvpiho.F

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