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 COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
23 DOUBLE PRECISION VALUE(20)
29 IF (MSTP(192) .GT. 0 .AND. MSTP(193) .GT. 0) THEN
33 VALUE(6)=MSTP(191)/1000
35 VALUE(7)=MOD(MSTP(191),1000)
36 CALL PDFSET(PARM,VALUE)
38 write(6,*) "-> pdfset"
39 CALL PDFSET(PARM,VALUE)
46 SUBROUTINE STRUCTM_ALICE
47 + (XX,QQ,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)
49 IMPLICIT DOUBLE PRECISION(A-H, O-Z)
51 COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
52 COMMON/PYINT1/MINT(400),VINT(400)
53 C write(6,*) "structm_alice->"
54 IF (MSTP(192) .GT. 0 .AND. MSTP(193) .GT. 0) THEN
56 C write(6,*) mint(124), "-> structa ", A
57 CALL STRUCTA(XX,QQ,A,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)
59 C write(6,*) mint(124), "-> structm "
60 CALL STRUCTM(XX,QQ,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)