]>
Commit | Line | Data |
---|---|---|
21886bb6 | 1 | C----------------------------------------------------------------------- |
2 | #include "pdf/pilot.h" | |
3 | SUBROUTINE SRCTQ47(DX,DSCAL, | |
4 | + DUPV,DDNV,DSEA,DSEAD,DSTR,DCHM,DBOT,DTOP,DGL) | |
5 | C | |
6 | C ******************************************************************** | |
7 | C * * | |
8 | C * Interface to CTEQset 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 CTEQ collaboration version 4 | |
15 | C :::::::::::: Lambda = 0.385 GeV (NF=3) (CTEQ4F3) | |
16 | C | |
17 | #include "pdf/impdp.inc" | |
18 | PARAMETER (ISET=13) | |
19 | #include "pdf/w50511.inc" | |
20 | C | |
21 | X = DX | |
22 | Q = DSCAL | |
23 | C | |
24 | C generate the individual structure fcn calls | |
25 | C | |
26 | UPV = Ctq4Fn(ISET, 1,X,Q) | |
27 | DNV = Ctq4Fn(ISET, 2,X,Q) | |
28 | SEAU = Ctq4Fn(ISET,-1,X,Q) | |
29 | DSEA = SEAU*X | |
30 | DUPV = (UPV-SEAU)*X | |
31 | SEAD = Ctq4Fn(ISET,-2,X,Q) | |
32 | DSEAD = SEAD*X | |
33 | DDNV = (DNV-SEAD)*X | |
34 | STR = Ctq4Fn(ISET,-3,X,Q) | |
35 | DSTR = STR*X | |
36 | CHM = Ctq4Fn(ISET,-4,X,Q) | |
37 | DCHM = CHM*X | |
38 | BOT = Ctq4Fn(ISET,-5,X,Q) | |
39 | DBOT = BOT*X | |
40 | TOP = 0. | |
41 | C IF (DSCAL.GT.TMAS) TOP = Ctq4Fn(ISET,-6,X,Q) | |
42 | DTOP = TOP*X | |
43 | GL = Ctq4Fn(ISET, 0,X,Q) | |
44 | DGL = GL*X | |
45 | C | |
46 | RETURN | |
47 | END |