DATA INI/8*0/I_LAST/-1/
SAVE
C
-
IHNT2(11)=JP
IHNT2(12)=JT
IF(IHNT2(5).NE.0 .AND. IHNT2(6).NE.0) THEN
C ********FINAL STATE RADIATION
IF(IHPR2(2).EQ.0.OR.IHPR2(2).EQ.2) MSTP(61)=0
IF(IHPR2(2).EQ.0.OR.IHPR2(2).EQ.1) MSTP(71)=0
+ MSTP(91) = 1
c
MSTP(81)=0
C ******** NO MULTIPLE INTERACTION
CKIN(3)=ABS(HIPR1(10))-0.25
CKIN(4)=ABS(HIPR1(10))+0.25
IF(HIPR1(10).LT.HIPR1(8)) CKIN(4)=-1.0
+C WRITE(6,*) "TRIGGERED", CKIN(3), CKIN(4), CKIN(5), IHPR2(3)
+C
+C ALICE
+C
+ IF (HIPR1(50) .EQ. 1) THEN
+ MSTP(61) = 0
+ MSTP(71) = 0
+ MSTP(91) = 0
+ ENDIF
+
c
MSEL=0
DO 101 ISUB=1,200
MDME(MDCY(21,2)+ISEL-1,1)=1
C ********QCD subprocesses
ELSE IF(IHPR2(3).EQ.2) THEN
+C WRITE(6,*) "DIRECT PHOTON PRODUCTION"
MSUB(14)=1
MSUB(18)=1
MSUB(29)=1
C ******************indicate for initialization use when
C structure functions are called in PYTHIA
C
+C WRITE(6,*) "CALL PYINIT WITH", CKIN(3), CKIN(4), MSTP(61), MSTP(71), MSTP(91)
CALL PYINIT_HIJING('CMS',BEAM,TARG,HINT1(1))
MINT4=MINT(44)
MINT5=MINT(45)
MINT(12)=NFT(JT,4)
RETURN
END
+
+
+
+
+
+
+
+
+