]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PDF/spdf/spctq51.F
Default shish-kebab geometry for EMCAL
[u/mrichter/AliRoot.git] / PDF / spdf / spctq51.F
CommitLineData
21886bb6 1#include "pdf/pilot.h"
2C-----------------------------------------------------------------------
3 SUBROUTINE SPCTQ51(DX,DSCAL,
4 + DUPV,DDNV,DSEA,DSEAD,DSTR,DCHM,DBOT,DTOP,DGL)
5C
6C ********************************************************************
7C * *
8C * Interface to CTEQset of structure functions *
9C * *
10C * Author: H. Plothow-Besch (CERN-PPE) *
11C * *
12C ********************************************************************
13C
14C :::::::::::: Structure functions from CETQ5M1 Parametrization
15C :::::::::::: Lambda = 0.326 GeV (MSbar) (CTEQ5M1)
16C
17#include "pdf/impdp.inc"
18C
19 PARAMETER (ISET=1)
20#include "pdf/w50511.inc"
21C
22 X = DX
23 Q = DSCAL
24C
25C generate the individual structure fcn calls
26C
27 UPV = Ctq5Mi(1,X,Q)
28 DNV = Ctq5Mi(2,X,Q)
29 SEAU = Ctq5Mi(-1,X,Q)
30 DSEA = SEAU * X
31 DUPV = (UPV-SEAU) * X
32 SEAD = Ctq5Mi(-2,X,Q)
33 DSEAD = SEAD * X
34 DDNV = (DNV-SEAD) * X
35 STR = Ctq5Mi(-3,X,Q)
36 DSTR = STR * X
37 CHM = Ctq5Mi(-4,X,Q)
38 DCHM = CHM * X
39 BOT = Ctq5Mi(-5,X,Q)
40 DBOT = BOT * X
41 TOP = 0.
42 IF (DSCAL.GT.TMAS) TOP = Ctq5Mi(-6,X,QL)
43 DTOP = TOP * X
44 GL = Ctq5Mi(0,X,Q)
45 DGL = GL * X
46C
47 RETURN
48 END