]>
Commit | Line | Data |
---|---|---|
21886bb6 | 1 | C----------------------------------------------------------------------- |
2 | #include "pdf/pilot.h" | |
3 | SUBROUTINE SFSAS12(DX,DQ2,DP2, | |
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 = 0.36 GeV**2 (MS_bar) | |
16 | C | |
7ef50f50 | 17 | #if defined(CERNLIB_DOUBLE) |
18 | DOUBLE PRECISION | |
19 | #endif | |
20 | #if defined(CERNLIB_SINGLE) | |
21 | REAL | |
22 | #endif | |
21886bb6 | 23 | + DX,DQ2,DP2, |
24 | + DUPV,DDNV,DSEA,DSEAD,DSTR,DCHM,DBOT,DTOP,DGL | |
25 | DIMENSION XPDFGM(-6:6) | |
26 | REAL X, Q, Q2, P2, F2GAM, XPDFGM | |
27 | PARAMETER (ISET=2) | |
28 | #include "pdf/w50511.inc" | |
29 | C | |
30 | X = DX | |
31 | Q = SQRT(DQ2) | |
32 | Q2 = DQ2 | |
33 | P2 = DP2 | |
34 | C | |
35 | C generate the individual structure fcn calls | |
36 | C | |
37 | CALL SASGAM1(ISET,X,Q2,P2,F2GAM,XPDFGM) | |
38 | UPV = XPDFGM(2) | |
39 | DUPV = UPV | |
40 | DNV = XPDFGM(1) | |
41 | DDNV = DNV | |
42 | SEAU = XPDFGM(-2) | |
43 | DSEA = SEAU | |
44 | SEAD = XPDFGM(-1) | |
45 | DSEAD = SEAD | |
46 | STR = XPDFGM(-3) | |
47 | DSTR = STR | |
48 | CHM = XPDFGM(-4) | |
49 | DCHM = CHM | |
50 | BOT = XPDFGM(-5) | |
51 | DBOT = BOT | |
52 | TOP = 0. | |
53 | C IF (DSCAL.GT.TMAS) TOP = XPDFGM(6) | |
54 | DTOP = TOP | |
55 | GL = XPDFGM(0) | |
56 | DGL = GL | |
57 | C | |
58 | RETURN | |
59 | END |