]> git.uio.no Git - u/mrichter/AliRoot.git/blob - PDF/spdf/sfsas11.F
Introduce factor for magnetic field
[u/mrichter/AliRoot.git] / PDF / spdf / sfsas11.F
1 C-----------------------------------------------------------------------
2 #include "pdf/pilot.h"
3       SUBROUTINE SFSAS11(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 (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=1)
23 #include "pdf/w50511.inc"
24 C
25       X  = DX
26       Q  = SQRT(DQ2)
27       Q2 = DQ2
28       P2 = DP2
29 C
30 C     generate the individual structure fcn calls
31 C
32       CALL SASGAM1(ISET,X,Q2,P2,F2GAM,XPDFGM)
33       UPV = XPDFGM(2)
34       DUPV = UPV
35       DNV = XPDFGM(1)
36       DDNV = DNV
37       SEAU = XPDFGM(-2)
38       DSEA = SEAU
39       SEAD = XPDFGM(-1)
40       DSEAD = SEAD
41       STR = XPDFGM(-3)
42       DSTR = STR
43       CHM = XPDFGM(-4)
44       DCHM = CHM
45       BOT = XPDFGM(-5)
46       DBOT = BOT
47       TOP = 0.
48 C      IF (DSCAL.GT.TMAS) TOP = XPDFGM(6)
49       DTOP = TOP
50       GL = XPDFGM(0)
51       DGL = GL
52 C
53       RETURN
54       END