]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | C----------------------------------------------------------------------- |
2 | #include "pdf/pilot.h" | |
3 | SUBROUTINE SFSAS23(DX,DQ2,DP2,IP2, | |
4 | + DUPV,DDNV,DSEA,DSEAD,DSTR,DCHM,DBOT,DTOP,DGL) | |
5 | C | |
6 | C ******************************************************************** | |
7 | C * * | |
8 | C * Interface to SASset of structure functions * | |
9 | C * * | |
10 | C * Author: H. Plothow-Besch (CERN-PPE) * | |
11 | C * * | |
12 | C ******************************************************************** | |
13 | C | |
14 | C :::::::::::: Structure functions from the SAS group version 2 | |
15 | C :::::::::::: Lambda = 0.200 GeV, Q**2 = 4 GeV**2 (DIS) | |
16 | C | |
17 | #include "pdf/expdp.inc" | |
18 | + DX,DQ2,DP2, | |
19 | + DUPV,DDNV,DSEA,DSEAD,DSTR,DCHM,DBOT,DTOP,DGL | |
20 | DIMENSION XPDFGM(-6:6) | |
21 | REAL X, Q, Q2, P2, F2GAM, XPDFGM | |
22 | PARAMETER (ISET=3) | |
23 | #include "pdf/w50511.inc" | |
24 | C | |
25 | X = DX | |
26 | Q = SQRT(DQ2) | |
27 | Q2 = DQ2 | |
28 | P2 = DP2 | |
29 | IP = IP2 | |
30 | IF(IP.LT.0 .OR. IP.GT.7) IP=0 | |
31 | C | |
32 | C generate the individual structure fcn calls | |
33 | C | |
34 | CALL SASGAM2(ISET,X,Q2,P2,IP,F2GAM,XPDFGM) | |
35 | UPV = XPDFGM(2) | |
36 | DUPV = UPV | |
37 | DNV = XPDFGM(1) | |
38 | DDNV = DNV | |
39 | SEAU = XPDFGM(-2) | |
40 | DSEA = SEAU | |
41 | SEAD = XPDFGM(-1) | |
42 | DSEAD = SEAD | |
43 | STR = XPDFGM(-3) | |
44 | DSTR = STR | |
45 | CHM = XPDFGM(-4) | |
46 | DCHM = CHM | |
47 | BOT = XPDFGM(-5) | |
48 | DBOT = BOT | |
49 | TOP = 0. | |
50 | C IF (DSCAL.GT.TMAS) TOP = XPDFGM(6) | |
51 | DTOP = TOP | |
52 | GL = XPDFGM(0) | |
53 | DGL = GL | |
54 | C | |
55 | RETURN | |
56 | END |