]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PYTHIA6/pdf_alice.F
Correct function Compare() for "pixels" from MLEM cluster finder.
[u/mrichter/AliRoot.git] / PYTHIA6 / pdf_alice.F
CommitLineData
2eef49c0 1C... ALICE interface to PDFLIB with possibility to select nuclear structure
2C... functions.
3C...
4C... The MSTP array in the PYPARS common block is used to enable and
5C... select the nuclear structure functions.
6C... MSTP(52) : (D=1) choice of proton and nuclear structure-function library
7C... =1: internal PYTHIA acording to MSTP(51)
8C... =2: PDFLIB proton s.f., with MSTP(51) = 1000xNGROUP+NSET
9C... MSTP( 51) = 1000xNPGROUP+NPSET
10C... MSTP(151) = 1000xNAGROUP+NASET
11C... MSTP(192) : Mass number of nucleus side 1
12C... MSTP(193) : Mass number of nucleus side 2
13C...
14C...
15C... MINT(124) : side (1 or 2)
16
17
18 SUBROUTINE PDFSET_ALICE(PARM, VALUE)
19C...
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)
48C...
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)
53C write(6,*) "structm_alice->"
54 IF (MSTP(192) .GT. 0 .AND. MSTP(193) .GT. 0) THEN
55 A=MSTP(191+MINT(124))
56C write(6,*) mint(124), "-> structa ", A
57 CALL STRUCTA(XX,QQ,A,UPV,DNV,USEA,DSEA,STR,CHM,BOT,TOP,GL)
58 ELSE
59C 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