EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
grv94di.F
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file grv94di.F
1 *
2 *...NLO PARAMETRIZATION (DIS) :
3 *
4 *72*********************************************************************
5  SUBROUTINE grv94di (X, Q2, UV, DV, DEL, UDB, SB, GL)
6  IMPLICIT DOUBLE PRECISION (a - z)
7  mu2 = 0.34
8  lam2 = 0.248 * 0.248
9  s = dlog(dlog(q2/lam2) / dlog(mu2/lam2))
10  ds = dsqrt(s)
11  s2 = s * s
12  s3 = s2 * s
13 *...UV :
14  nu = 2.484 + 0.116 * s + 0.093 * s2
15  aku = 0.563 - 0.025 * s
16  bku = 0.054 + 0.154 * s
17  au = -0.326 - 0.058 * s - 0.135 * s2
18  bu = -3.322 + 8.259 * s - 3.119 * s2 + 0.291 * s3
19  cu = 11.52 - 12.99 * s + 3.161 * s2
20  du = 2.808 + 1.400 * s - 0.557 * s2 + 0.119 * s3
21  uv = fv(x, nu, aku, bku, au, bu, cu, du)
22 *...DV :
23  nd = 0.156 - 0.017 * s
24  akd = 0.299 - 0.022 * s
25  bkd = 0.259 - 0.015 * s
26  ad = 3.445 + 1.278 * s + 0.326 * s2
27  bd = -6.934 + 37.45 * s - 18.95 * s2 + 1.463 * s3
28  cd = 55.45 - 69.92 * s + 20.78 * s2
29  dd = 3.577 + 1.441 * s - 0.683 * s2 + 0.179 * s3
30  dv = fv(x, nd, akd, bkd, ad, bd, cd, dd)
31 *...DEL :
32  ne = 0.099 + 0.019 * s + 0.002 * s2
33  ake = 0.419 - 0.013 * s
34  bke = 1.064 - 0.038 * s
35  ae = -44.00 + 98.70 * s - 14.79 * s2
36  be = 28.59 - 40.94 * s - 13.66 * s2 + 2.523 * s3
37  ce = 84.57 - 108.8 * s + 31.52 * s2
38  de = 7.469 + 2.480 * s - 0.866 * s2
39  del = fv(x, ne, ake, bke, ae, be, ce, de)
40 *...UDB :
41  alx = 1.215
42  bex = 0.466
43  akx = 0.326 + 0.150 * s
44  bkx = 0.956 + 0.405 * s
45  agx = 0.272
46  bgx = 3.794 - 2.359 * ds
47  cx = 2.014
48  dx = 7.941 + 0.534 * ds - 0.940 * s + 0.410 * s2
49  ex = 3.049 + 1.597 * s
50  esx = 4.396 - 4.594 * ds + 3.268 * s
51  udb = fw(x, s, alx, bex, akx, bkx, agx, bgx, cx, dx, ex, esx)
52 *...SB :
53  als = 0.175
54  bes = 0.344
55  aks = 1.415 - 0.641 * ds
56  as = 0.580 - 9.763 * ds + 6.795 * s - 0.558 * s2
57  bs = 5.617 + 5.709 * ds - 3.972 * s
58  dst = 13.78 - 9.581 * s + 5.370 * s2 - 0.996 * s3
59  est = 4.546 + 0.372 * s2
60  ess = 5.053 - 1.070 * s + 0.805 * s2
61  sb = fws(x, s, als, bes, aks, as, bs, dst, est, ess)
62 *...GL :
63  alg = 1.258
64  beg = 1.846
65  akg = 2.423
66  bkg = 2.427 + 1.311 * s - 0.153 * s2
67  ag = 25.09 - 7.935 * s
68  bg = -14.84 - 124.3 * ds + 72.18 * s
69  cg = 590.3 - 173.8 * s
70  dg = 5.196 + 1.857 * s
71  eg = -1.648 + 3.988 * s - 0.432 * s2
72  esg = 3.232 - 0.542 * s
73  gl = fw(x, s, alg, beg, akg, bkg, ag, bg, cg, dg, eg, esg)
74  RETURN
75  END