]>
Commit | Line | Data |
---|---|---|
82cedcfe | 1 | #include "pdf/pilot.h" |
2 | C*************************************************************************** | |
3 | C | |
4 | C sfeks98.f | |
5 | C | |
6 | C An interface for calculating the SCALE DEPENDENT NUCLEAR PARTON | |
7 | C DISTRIBUTIONS | |
8 | C f_A(x,Q) | |
9 | C where f_A is the distribution of parton flavour f in a proton of a | |
10 | C nucleus A, and f_p is the corresponding parton distribution in the | |
11 | C free proton. | |
12 | C | |
13 | C Uses eks98 nuclear parton ratios and PDFLIB parton distributions. | |
14 | C | |
15 | C | |
16 | C------------------------------------------------------------------- | |
17 | subroutine sfeks98 | |
18 | + (x,Q,A,uvA,dvA,usA,dsA,stA,chA,botA,topA,glA) | |
19 | #include "pdf/impdp.inc" | |
20 | data zerod/0.d0/ | |
21 | C external structm,eks98 | |
22 | dx = x | |
23 | dq = q | |
24 | call structm(dx,dq,uv,dv,us,ds,st,ch,bot,top,gl) | |
25 | dx = x | |
26 | dq = q | |
27 | dA = A | |
28 | call eks98(dx,dq,dA,Ruv,Rdv,Rub,Rdb,Rs,Rc,Rb,Rt,Rg) | |
29 | ||
30 | uvA = Ruv*uv | |
31 | dvA = Rdv*dv | |
32 | usA = Rub*us | |
33 | dsA = Rdb*ds | |
34 | stA = Rs*st | |
35 | chA = Rc*ch | |
36 | botA= Rb*bot | |
37 | topA= Rt*top | |
38 | glA = Rg*gl | |
39 | ||
40 | return | |
41 | end | |
42 | ||
43 | ||
44 | ||
45 |