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)
49 IF (MSTP(194) .EQ. 0) THEN
50 CALL SETLHAPARM("EKS98")
52 CALL SETLHAPARM("EPS08")
55 write(6,*) "-> pdfset"
56 CALL PDFSET(PARM,VALUE,
57 > MSTU(11),MSTP(51),MSTP(53),MSTP(55),
59 > XMIN,XMAX,Q2MIN,Q2MAX)
66 SUBROUTINE STRUCTM_ALICE
67 + (XX,QQ,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)
69 IMPLICIT DOUBLE PRECISION(A-H, O-Z)
71 COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
72 COMMON/PYINT1/MINT(400),VINT(400)
73 C write(6,*) "structm_alice->"
74 IF (MSTP(192) .GT. 0 .AND. MSTP(193) .GT. 0) THEN
76 C write(6,*) mint(124), "-> structa ", A
77 CALL STRUCTA(XX,QQ,A,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)
79 C write(6,*) mint(124), "-> structm "
80 CALL STRUCTM(XX,QQ,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)