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