3 C*********************************************************************
5 FUNCTION KLU_HIJING(I,J)
7 C...Purpose: to provide various integer-valued event related data.
8 #include "lujets_hijing.inc"
9 #include "ludat1_hijing.inc"
10 #include "ludat2_hijing.inc"
12 C...Default value. For I=0 number of entries, number of stable entries
13 C...or 3 times total charge.
15 IF(I.LT.0.OR.I.GT.MSTU(4).OR.J.LE.0) THEN
16 ELSEIF(I.EQ.0.AND.J.EQ.1) THEN
18 ELSEIF(I.EQ.0.AND.(J.EQ.2.OR.J.EQ.6)) THEN
20 IF(J.EQ.2.AND.K(I1,1).GE.1.AND.K(I1,1).LE.10) KLU_HIJING
22 IF(J.EQ.6.AND.K(I1,1).GE.1.AND.K(I1,1).LE.10) KLU_HIJING
23 $ =KLU_HIJING+LUCHGE_HIJING(K(I1,2))
27 C...For I > 0 direct readout of K matrix or charge.
31 KLU_HIJING=LUCHGE_HIJING(K(I,2))
33 C...Status (existing/fragmented/decayed), parton/hadron separation.
35 IF(K(I,1).GE.1.AND.K(I,1).LE.10) KLU_HIJING=1
36 IF(J.EQ.8) KLU_HIJING=KLU_HIJING*K(I,2)
41 IF(KC.NE.0) KQ=KCHG(KC,2)
42 IF(J.EQ.9.AND.KC.NE.0.AND.KQ.NE.0) KLU_HIJING=K(I,2)
43 IF(J.EQ.10.AND.KC.NE.0.AND.KQ.EQ.0) KLU_HIJING=K(I,2)
44 IF(J.EQ.11) KLU_HIJING=KC
45 IF(J.EQ.12) KLU_HIJING=KQ*ISIGN(1,K(I,2))
47 C...Heaviest flavour in hadron/diquark.
50 KLU_HIJING=MOD(KFA/100,10)*(-1)**MOD(KFA/100,10)
51 IF(KFA.LT.10) KLU_HIJING=KFA
52 IF(MOD(KFA/1000,10).NE.0) KLU_HIJING=MOD(KFA/1000,10)
53 KLU_HIJING=KLU_HIJING*ISIGN(1,K(I,2))
55 C...Particle history: generation, ancestor, rank.
59 110 KLU_HIJING=KLU_HIJING+1
63 IF(I1.GT.0.AND.K(I1,1).GT.0.AND.K(I1,1).LE.20) GOTO 110
64 IF(J.EQ.15) KLU_HIJING=I2
68 120 IF(K(I1,3).EQ.I2.AND.K(I1,1).GT.0.AND.K(I1,1).LE.20)
69 $ KLU_HIJING=KLU_HIJING+1
72 C...Particle coming from collapsing jet system or not.
75 130 KLU_HIJING=KLU_HIJING+1
79 KC=LUCOMP_HIJING(K(I0,2))
80 IF(I1.EQ.0.OR.K(I0,1).LE.0.OR.K(I0,1).GT.20.OR.KC.EQ.0) THEN
81 IF(KLU_HIJING.EQ.1) KLU_HIJING=-1
82 IF(KLU_HIJING.GT.1) KLU_HIJING=0
85 IF(KCHG(KC,2).EQ.0) GOTO 130
86 IF(K(I1,1).NE.12) KLU_HIJING=0
87 IF(K(I1,1).NE.12) RETURN
90 IF(I2.LT.N.AND.K(I2,1).NE.11) GOTO 140
92 IF(K3M.GE.I1.AND.K3M.LE.I2) KLU_HIJING=0
94 IF(I3.LT.N.AND.K3P.GE.I1.AND.K3P.LE.I2) KLU_HIJING=0
96 C...Number of decay products. Colour flow.
98 IF(K(I,1).EQ.11.OR.K(I,1).EQ.12) KLU_HIJING=MAX(0,K(I,5)-K(I,4)
100 IF(K(I,4).EQ.0.OR.K(I,5).EQ.0) KLU_HIJING=0
102 IF(K(I,1).NE.3.AND.K(I,1).NE.13.AND.K(I,1).NE.14) RETURN
103 IF(J.EQ.19) KLU_HIJING=MOD(K(I,4)/MSTU(5),MSTU(5))
104 IF(J.EQ.20) KLU_HIJING=MOD(K(I,5)/MSTU(5),MSTU(5))
105 IF(J.EQ.21) KLU_HIJING=MOD(K(I,4),MSTU(5))
106 IF(J.EQ.22) KLU_HIJING=MOD(K(I,5),MSTU(5))