1 C... ALICE interface to PDFLIB with possibility to select nuclear structure
4 C... The MSTP array in the PYPARS common block is used to enable and
5 C... select the nuclear structure functions.
6 C... MSTP(52) : (D=1) choice of proton and nuclear structure-function library
7 C... =1: internal PYTHIA acording to MSTP(51)
8 C... =2: PDFLIB proton s.f., with MSTP(51) = 1000xNGROUP+NSET
9 C... MSTP( 51) = 1000xNPGROUP+NPSET
10 C... MSTP(151) = 1000xNAGROUP+NASET
11 C... MSTP(192) : Mass number of nucleus side 1
12 C... MSTP(193) : Mass number of nucleus side 2
15 C... MINT(124) : side (1 or 2)
18 SUBROUTINE PDFSET_ALICE(PARM, VALUE)
20 IMPLICIT DOUBLE PRECISION(A-H, O-Z)
22 C...Interface to PDFLIB.
23 COMMON/LW50512/QCDL4,QCDL5
25 DOUBLE PRECISION QCDL4,QCDL5
26 COMMON/LW50513/XMIN,XMAX,Q2MIN,Q2MAX
28 DOUBLE PRECISION XMIN,XMAX,Q2MIN,Q2MAX
30 COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
31 COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
32 DOUBLE PRECISION VALUE(20)
38 IF (MSTP(192) .GT. 0 .AND. MSTP(193) .GT. 0) THEN
42 VALUE(6)=MSTP(191)/1000
44 VALUE(7)=MOD(MSTP(191),1000)
45 CALL PDFSET(PARM,VALUE,
46 > MSTU(11),MSTP(51),MSTP(53),MSTP(55),
48 > XMIN,XMAX,Q2MIN,Q2MAX)
50 IF (MSTP(194) .EQ. 0) THEN
51 CALL SETLHAPARM("EKS98")
52 ELSE IF (MSTP(194) .EQ. 9) THEN
53 CALL SETLHAPARM("EPS09LO")
54 ELSE IF (MSTP(194) .EQ. 19) THEN
55 CALL SETLHAPARM("EPS09NLO")
56 ELSE IF (MSTP(194) .EQ. 8) THEN
57 CALL SETLHAPARM("EPS08")
59 CALL SETLHAPARM("EPS09LO")
62 write(6,*) "-> pdfset"
63 CALL PDFSET(PARM,VALUE,
64 > MSTU(11),MSTP(51),MSTP(53),MSTP(55),
66 > XMIN,XMAX,Q2MIN,Q2MAX)
73 SUBROUTINE STRUCTM_ALICE
74 + (XX,QQ,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)
76 IMPLICIT DOUBLE PRECISION(A-H, O-Z)
78 COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
79 COMMON/PYINT1/MINT(400),VINT(400)
80 C write(6,*) "structm_alice->"
81 IF (MSTP(192) .GT. 0 .AND. MSTP(193) .GT. 0) THEN
83 C write(6,*) mint(124), "-> structa ", A
84 CALL STRUCTA(XX,QQ,A,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)
86 C write(6,*) mint(124), "-> structm "
87 CALL STRUCTM(XX,QQ,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)