PDF/spdf/grv94ho.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:35  cernlib
   9 * Version 7.04
  10 *
  11 * Revision 1.1.1.1  1996/04/12 15:29:22  plothow
  12 * Version 7.01
  13 *
  14 *
  15 #include "pdf/pilot.h"
  16 *
  17 *...NLO PARAMETRIZATION (MS(BAR)) :
  18 *
  19        SUBROUTINE GRV94HO (ZX,ZQ,ZUV,ZDV,ZUB,ZDB,ZSB,ZCB,ZBB,ZTB,ZGL)
  20        IMPLICIT REAL (A - Y)
  21 C+SEQ, IMPDP.
  22 #if defined(CERNLIB_DOUBLE)
  23       DOUBLE PRECISION
  24 #endif
  25 #if defined(CERNLIB_SINGLE)
  26       REAL
  27 #endif
  28      +      ZX,ZQ,ZUV,ZDV,ZDEL,ZUDB,ZUB,ZDB,ZSB,ZCB,ZBB,ZTB,ZGL
  29        REAL X, Q
  30        X = ZX
  31        Q = ZQ
  32        MU2  = 0.34
  33        LAM2 = 0.248   * 0.248
  34        Q2 = Q*Q
  35        ARG1 = MAX(MU2,Q2)/LAM2
  36        ARG2 = MU2/LAM2
  37        S  =  LOG ( LOG( ARG1  ) /  LOG( ARG2   ))
  38        DS =  SQRT (S)
  39        S2 = S * S
  40        S3 = S2 * S
  41 *...UV :
  42        NU  =  1.304 + 0.863 * S
  43        AKU =  0.558 - 0.020 * S
  44        BKU =          0.183 * S
  45        AU  = -0.113 + 0.283 * S - 0.321 * S2
  46        BU  =  6.843 - 5.089 * S + 2.647 * S2 - 0.527 * S3
  47        CU  =  7.771 - 10.09 * S + 2.630 * S2
  48        DU  =  3.315 + 1.145 * S - 0.583 * S2 + 0.154 * S3
  49        UV  = GR94FV (X, NU, AKU, BKU, AU, BU, CU, DU)
  50        ZUV = UV
  51 *...DV :
  52        ND  =  0.102 - 0.017 * S + 0.005 * S2
  53        AKD =  0.270 - 0.019 * S
  54        BKD =  0.260
  55        AD  =  2.393 + 6.228 * S - 0.881 * S2
  56        BD  =  46.06 + 4.673 * S - 14.98 * S2 + 1.331 * S3
  57        CD  =  17.83 - 53.47 * S + 21.24 * S2
  58        DD  =  4.081 + 0.976 * S - 0.485 * S2 + 0.152 * S3
  59        DV  = GR94FV (X, ND, AKD, BKD, AD, BD, CD, DD)
  60        ZDV = DV
  61 *...DEL :
  62        NE  =  0.070 + 0.042 * S - 0.011 * S2 + 0.004 * S3
  63        AKE =  0.409 - 0.007 * S
  64        BKE =  0.782 + 0.082 * S
  65        AE  = -29.65 + 26.49 * S + 5.429 * S2
  66        BE  =  90.20 - 74.97 * S + 4.526 * S2
  67        CE  =  0.0
  68        DE  =  8.122 + 2.120 * S - 1.088 * S2 + 0.231 * S3
  69        DEL = GR94FV (X, NE, AKE, BKE, AE, BE, CE, DE)
  70        ZDEL= DEL
  71 *...UDB :
  72        ALX =  0.877
  73        BEX =  0.561
  74        AKX =  0.275
  75        BKX =  0.0
  76        AGX =  0.997
  77        BGX =  3.210 - 1.866 * S
  78        CX  =  7.300
  79        DX  =  9.010 + 0.896 * DS + 0.222 * S2
  80        EX  =  3.077 + 1.446 * S
  81        ESX =  3.173 - 2.445 * DS + 2.207 * S
  82        UDB = GR94FW(X, S, ALX, BEX, AKX, BKX, AGX, BGX, CX, DX, EX, ESX)
  83        ZUDB= UDB
  84 *...UB
  85        ZUB = (ZUDB - ZDEL)/2.D0
  86 *...DB
  87        ZDB = (ZUDB + ZDEL)/2.D0
  88 *...SB :
  89        ALS =  0.756
  90        BES =  0.216
  91        AKS =  1.690 + 0.650 * DS - 0.922 * S
  92        AS  = -4.329 + 1.131 * S
  93        BS  =  9.568 - 1.744 * S
  94        DST =  9.377 + 1.088 * DS - 1.320 * S + 0.130 * S2
  95        EST =  3.031 + 1.639 * S
  96        ESS =  5.837 + 0.815 * S
  97        SB  = GR94FWS (X, S, ALS, BES, AKS, AS, BS, DST, EST, ESS)
  98        ZSB = SB
  99 *...GL :
 100        ALG =  1.014
 101        BEG =  1.738
 102        AKG =  1.724 + 0.157 * S
 103        BKG =  0.800 + 1.016 * S
 104        AG  =  7.517 - 2.547 * S
 105        BG  =  34.09 - 52.21 * DS + 17.47 * S
 106        CG  =  4.039 + 1.491 * S
 107        DG  =  3.404 + 0.830 * S
 108        EG  = -1.112 + 3.438 * S  - 0.302 * S2
 109        ESG =  3.256 - 0.436 * S
 110        GL  = GR94FW (X, S, ALG, BEG, AKG, BKG, AG, BG, CG, DG, EG, ESG)
 111        ZGL = GL
 112 *...
 113        ZCB = 0.D0
 114        ZBB = 0.D0
 115        ZTB = 0.D0
 116 *...
 117        RETURN
 118        END