CHARACTER BEAM*16,TARG*16
DIMENSION XSEC0(8,0:200),COEF0(8,200,20),INI(8),
& MINT44(8),MINT45(8)
+#define BLANKET_SAVE
#include "hijcrdn.inc"
#include "hiparnt.inc"
#include "histrng.inc"
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
ENDIF
ENDIF
c
+
IF(I_TRIG.NE.0) GO TO 160
IF(I_TRIG.EQ.I_LAST) GO TO 150
MSTP(2)=2
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
MSUB(28)=1
MSUB(53)=1
MSUB(68)=1
- MSUB(81)=1
- MSUB(82)=1
+ MSUB(81)=1
+ MSUB(82)=1
DO 110 J=1,MIN(8,MDCY(21,3))
110 MDME(MDCY(21,2)+J-1,1)=0
ISEL=4
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)
+C ALICE ->
+ IF (IHPR2(49) .EQ. 1) THEN
+C SWITCH OFF MASSIVE FINAL QUARKS
+ MSEL = 0
+C QCD on
+ MSUB(11) = 1
+ MSUB(12) = 1
+ MSUB(13) = 1
+ MSUB(28) = 1
+ MSUB(53) = 1
+ MSUB(68) = 1
+C HF OFF
+ MSUB(81) = 0
+ MSUB(82) = 0
+ MDME(127, 1) = 1
+ MDME(128, 1) = 1
+ MDME(129, 1) = 1
+C g->QQbar off
+ MDME(130, 1) = 0
+ MDME(131, 1) = 0
+C highest flavor in final state shower
+ MSTJ(45) = 3
+C highest flavor in initial state shower and pdf
+ MSTP(54) = 3
+ ENDIF
+C <- ALICE
CALL PYINIT_HIJING('CMS',BEAM,TARG,HINT1(1))
MINT4=MINT(44)
MINT5=MINT(45)
MINT(12)=NFT(JT,4)
RETURN
END
+
+
+
+
+
+
+
+
+