EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
grv94ho.F
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file grv94ho.F
1 *
2 *...NLO PARAMETRIZATION (MS(BAR)) :
3 *
4 *72*********************************************************************
5  SUBROUTINE grv94ho (X, Q2, UV, DV, DEL, UDB, SB, GL)
6 
7  IMPLICIT DOUBLE PRECISION (a - z)
8  mu2 = 0.34
9  lam2 = 0.248 * 0.248
10  s = dlog(dlog(q2/lam2) / dlog(mu2/lam2))
11  ds = dsqrt(s)
12  s2 = s * s
13  s3 = s2 * s
14 *...UV :
15  nu = 1.304 + 0.863 * s
16  aku = 0.558 - 0.020 * s
17  bku = 0.183 * s
18  au = -0.113 + 0.283 * s - 0.321 * s2
19  bu = 6.843 - 5.089 * s + 2.647 * s2 - 0.527 * s3
20  cu = 7.771 - 10.09 * s + 2.630 * s2
21  du = 3.315 + 1.145 * s - 0.583 * s2 + 0.154 * s3
22  uv = fv(x, nu, aku, bku, au, bu, cu, du)
23 *...DV :
24  nd = 0.102 - 0.017 * s + 0.005 * s2
25  akd = 0.270 - 0.019 * s
26  bkd = 0.260
27  ad = 2.393 + 6.228 * s - 0.881 * s2
28  bd = 46.06 + 4.673 * s - 14.98 * s2 + 1.331 * s3
29  cd = 17.83 - 53.47 * s + 21.24 * s2
30  dd = 4.081 + 0.976 * s - 0.485 * s2 + 0.152 * s3
31  dv = fv(x, nd, akd, bkd, ad, bd, cd, dd)
32 *...DEL :
33  ne = 0.070 + 0.042 * s - 0.011 * s2 + 0.004 * s3
34  ake = 0.409 - 0.007 * s
35  bke = 0.782 + 0.082 * s
36  ae = -29.65 + 26.49 * s + 5.429 * s2
37  be = 90.20 - 74.97 * s + 4.526 * s2
38  ce = 0.0
39  de = 8.122 + 2.120 * s - 1.088 * s2 + 0.231 * s3
40  del = fv(x, ne, ake, bke, ae, be, ce, de)
41 *...UDB :
42  alx = 0.877
43  bex = 0.561
44  akx = 0.275
45  bkx = 0.0
46  agx = 0.997
47  bgx = 3.210 - 1.866 * s
48  cx = 7.300
49  dx = 9.010 + 0.896 * ds + 0.222 * s2
50  ex = 3.077 + 1.446 * s
51  esx = 3.173 - 2.445 * ds + 2.207 * s
52  udb = fw(x, s, alx, bex, akx, bkx, agx, bgx, cx, dx, ex, esx)
53 *...SB :
54  als = 0.756
55  bes = 0.216
56  aks = 1.690 + 0.650 * ds - 0.922 * s
57  as = -4.329 + 1.131 * s
58  bs = 9.568 - 1.744 * s
59  dst = 9.377 + 1.088 * ds - 1.320 * s + 0.130 * s2
60  est = 3.031 + 1.639 * s
61  ess = 5.837 + 0.815 * s
62  sb = fws(x, s, als, bes, aks, as, bs, dst, est, ess)
63 *...GL :
64  alg = 1.014
65  beg = 1.738
66  akg = 1.724 + 0.157 * s
67  bkg = 0.800 + 1.016 * s
68  ag = 7.517 - 2.547 * s
69  bg = 34.09 - 52.21 * ds + 17.47 * s
70  cg = 4.039 + 1.491 * s
71  dg = 3.404 + 0.830 * s
72  eg = -1.112 + 3.438 * s - 0.302 * s2
73  esg = 3.256 - 0.436 * s
74  gl = fw(x, s, alg, beg, akg, bkg, ag, bg, cg, dg, eg, esg)
75  RETURN
76  END