[u/mrichter/AliRoot.git] / PYTHIA6 / pdf_alice.F

C...  ALICE interface to PDFLIB with possibility to select nuclear structure 
C...  functions. 
C...  
C...  The MSTP array in the PYPARS common block is used to enable and 
C...  select the nuclear structure functions. 
C...  MSTP(52)  : (D=1) choice of proton and nuclear structure-function library
C...          =1: internal PYTHIA acording to MSTP(51) 
C...          =2: PDFLIB proton  s.f., with MSTP(51)  = 1000xNGROUP+NSET
C...              MSTP( 51)  = 1000xNPGROUP+NPSET
C...              MSTP(151)  = 1000xNAGROUP+NASET
C...  MSTP(192) : Mass number of nucleus side 1
C...  MSTP(193) : Mass number of nucleus side 2
C...
C...
C...  MINT(124) : side (1 or 2)


      SUBROUTINE PDFSET_ALICE(PARM, VALUE)
C...
      IMPLICIT DOUBLE PRECISION(A-H, O-Z)
      IMPLICIT INTEGER(I-N)
      COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)  
      DOUBLE PRECISION VALUE(20)
      CHARACTER*20 PARM(20)
      write(6,*) MSTP(52)
      write(6,*) PARM
      write(6,*) VALUE

      IF (MSTP(192) .GT. 0 .AND. MSTP(193) .GT. 0) THEN
         PARM(5)='NATYPE'
         VALUE(5)=4
         PARM(6)='NAGROUP'
         VALUE(6)=MSTP(191)/1000
         PARM(7)='NASET'
         VALUE(7)=MOD(MSTP(191),1000)
         CALL PDFSET(PARM,VALUE)
      ELSE 
         write(6,*) "-> pdfset"
         CALL PDFSET(PARM,VALUE)
      ENDIF
      write(6,*) "done"
      END


      SUBROUTINE STRUCTM_ALICE
     +     (XX,QQ,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)
C...
      IMPLICIT DOUBLE PRECISION(A-H, O-Z)
      IMPLICIT INTEGER(I-N)
      COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
      COMMON/PYINT1/MINT(400),VINT(400)
C      write(6,*) "structm_alice->"
      IF (MSTP(192) .GT. 0 .AND. MSTP(193) .GT. 0) THEN
         A=MSTP(191+MINT(124))
C         write(6,*) mint(124), "-> structa ", A
          CALL STRUCTA(XX,QQ,A,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)
      ELSE
C         write(6,*) mint(124), "-> structm "
         CALL STRUCTM(XX,QQ,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)
      ENDIF
      END
Commit	Line	Data
2eef49c0	1	C... ALICE interface to PDFLIB with possibility to select nuclear structure
	2	C... functions.
	3	C...
	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
	13	C...
	14	C...
	15	C... MINT(124) : side (1 or 2)
	16
	17
	18	SUBROUTINE PDFSET_ALICE(PARM, VALUE)
	19	C...
	20	IMPLICIT DOUBLE PRECISION(A-H, O-Z)
	21	IMPLICIT INTEGER(I-N)
	22	COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
	23	DOUBLE PRECISION VALUE(20)
	24	CHARACTER*20 PARM(20)
	25	write(6,*) MSTP(52)
	26	write(6,*) PARM
	27	write(6,*) VALUE
	28
	29	IF (MSTP(192) .GT. 0 .AND. MSTP(193) .GT. 0) THEN
	30	PARM(5)='NATYPE'
	31	VALUE(5)=4
	32	PARM(6)='NAGROUP'
	33	VALUE(6)=MSTP(191)/1000
	34	PARM(7)='NASET'
	35	VALUE(7)=MOD(MSTP(191),1000)
	36	CALL PDFSET(PARM,VALUE)
	37	ELSE
	38	write(6,*) "-> pdfset"
	39	CALL PDFSET(PARM,VALUE)
	40	ENDIF
	41	write(6,*) "done"
	42	END
	43
	44
	45
	46	SUBROUTINE STRUCTM_ALICE
	47	+ (XX,QQ,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)
	48	C...
	49	IMPLICIT DOUBLE PRECISION(A-H, O-Z)
	50	IMPLICIT INTEGER(I-N)
	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
	55	A=MSTP(191+MINT(124))
	56	C write(6,*) mint(124), "-> structa ", A
	57	CALL STRUCTA(XX,QQ,A,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)
	58	ELSE
	59	C write(6,*) mint(124), "-> structm "
	60	CALL STRUCTM(XX,QQ,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)
	61	ENDIF
	62	END
	63
	64