8 parameter(nmxhep=10000)
9 common/hepevt/nevhep,nhep,isthep(nmxhep),idhep(nmxhep),
10 &jmohep(2,nmxhep),jdahep(2,nmxhep),phep(5,nmxhep),vhep(4,nmxhep)
11 DOUBLE PRECISION phep,vhep
13 common/lujets/
n,
k(9000,5),
p(9000,5),
v(9000,5)
15 common/ludat1/mstu(200),paru(200),mstj(200),parj(200)
17 common/ludat2/kchg(500,3),pmas(500,4),parf(2000),vckm(4,4)
24 &
'(LUHEPC:) no more space in /HEPEVT/')
28 IF(
k(
i,1).GE.1.AND.
k(
i,1).LE.10) isthep(
i)=1
29 IF(
k(
i,1).GE.11.AND.
k(
i,1).LE.20) isthep(
i)=2
30 IF(
k(
i,1).GE.21.AND.
k(
i,1).LE.30) isthep(
i)=3
31 IF(
k(
i,1).GE.31.AND.
k(
i,1).LE.100) isthep(
i)=
k(
i,1)
35 IF(
k(
i,1).NE.3.AND.
k(
i,1).NE.13.AND.
k(
i,1).NE.14)
THEN
48 IF(
k(
i,2).GE.91.AND.
k(
i,2).LE.93)
THEN
52 &
'(LUHEPC:) translation of inconsistent event history')
53 IF(
i1.LT.
i.AND.
k(
i1,1).NE.1.AND.
k(
i1,1).NE.11) goto 120
55 IF(
i1.LT.
i.AND.kc.EQ.0) goto 120
56 IF(
i1.LT.
i.AND.kchg(kc,2).EQ.0) goto 120
58 ELSEIF(
k(
i,2).EQ.94)
THEN
60 IF(nhep.GE.
i+3.AND.
k(
i+3,3).LE.
i) njet=3
61 IF(nhep.GE.
i+4.AND.
k(
i+4,3).LE.
i) njet=4
62 jmohep(2,
i)=mod(
k(
i+njet,4)/mstu(5),mstu(5))
63 IF(jmohep(2,
i).EQ.jmohep(1,
i)) jmohep(2,
i)=
64 & mod(
k(
i+1,4)/mstu(5),mstu(5))
68 IF(
k(
i,2).EQ.94.AND.mstu(16).NE.2)
THEN
69 DO 130
i1=jdahep(1,
i),jdahep(2,
i)
70 i2=mod(
k(
i1,4)/mstu(5),mstu(5))
73 IF(
k(
i,2).GE.91.AND.
k(
i,2).LE.94) goto 140
75 IF(
i1.LE.0.OR.
i1.GT.nhep) goto 140
76 IF(
k(
i1,1).NE.13.AND.
k(
i1,1).NE.14) goto 140
77 IF(jdahep(1,
i1).EQ.0)
THEN
84 IF(
k(
i,1).NE.13.AND.
k(
i,1).NE.14) goto 150
85 IF(jdahep(2,
i).EQ.0) jdahep(2,
i)=jdahep(1,
i)
90 IF(nhep.GT.mstu(4)) CALL
luerrm(8,
91 &
'(LUHEPC:) no more space in /LUJETS/')
97 IF(isthep(
i).EQ.1)
k(
i,1)=1
98 IF(isthep(
i).EQ.2)
k(
i,1)=11
99 IF(isthep(
i).EQ.3)
k(
i,1)=21
109 IF(isthep(
i).EQ.2.AND.phep(4,
i).GT.phep(5,
i))
THEN
111 IF(
i1.GT.0.AND.
i1.LE.nhep)
v(
i,5)=(vhep(4,
i1)-vhep(4,
i))*
112 & phep(5,
i)/phep(4,
i)
116 IF(isthep(
i).EQ.1)
THEN
119 IF(kc.NE.0) kq=kchg(kc,2)*isign(1,
k(
i,2))
120 IF(kq.NE.0) nkq=nkq+1
121 IF(kq.NE.2) kqsum=kqsum+kq
122 IF(kq.NE.0.AND.kqsum.NE.0)
THEN
124 ELSEIF(kq.EQ.2.AND.
i.LT.
n)
THEN
125 IF(
k(
i+1,2).EQ.21)
k(
i,1)=2
129 IF(nkq.EQ.1.OR.kqsum.NE.0) CALL
luerrm(8,
130 &
'(LUHEPC:) input parton configuration not colour singlet')