9 dimension kf(100),
px(100),
py(100),
pz(100),pe(100),pm(100)
10 dimension
y(100),ip(100,2)
17 IF(npt.EQ.2) go to 500
27 y(i-iq)=0.5*alog((abs(pe(i)+
pz(i))+1.
e-5)
28 & /(abs(pe(i)-
pz(i))+1.
e-5))
43 IF(i.LE.
npj(jp)) go to 100
46 DO 200 j=i+1,
npj(jp)-iq
47 IF(
y(i).GT.
y(j)) go to 200
58 300
kfpj(jp,i)=kf(ip(i-iqq,1))
62 pjpe(jp,i)=pe(ip(i-iqq,1))
63 pjpm(jp,i)=pm(ip(i-iqq,1))
64 IF(ip(i-iqq,2).EQ.1)
THEN
65 kfpj(jp,i+1)=kf(ip(i-iqq,1)+1)
66 pjpx(jp,i+1)=
px(ip(i-iqq,1)+1)
67 pjpy(jp,i+1)=
py(ip(i-iqq,1)+1)
68 pjpz(jp,i+1)=
pz(ip(i-iqq,1)+1)
69 pjpe(jp,i+1)=pe(ip(i-iqq,1)+1)
70 pjpm(jp,i+1)=pm(ip(i-iqq,1)+1)
75 IF(i.LE.
npj(jp)) go to 300
88 y(i-iq)=0.5*alog((abs(pe(i)+
pz(i))+1.
e-5)
89 & /(abs(pe(i)-
pz(i))+1.
e-5))
104 IF(i.LE.
ntj(jt)) go to 600
106 DO 700 i=1,
ntj(jt)-iq
107 DO 700 j=i+1,
ntj(jt)-iq
108 IF(
y(i).LT.
y(j)) go to 700
119 800
kftj(jt,i)=kf(ip(i-iqq,1))
120 pjtx(jt,i)=
px(ip(i-iqq,1))
121 pjty(jt,i)=
py(ip(i-iqq,1))
122 pjtz(jt,i)=
pz(ip(i-iqq,1))
123 pjte(jt,i)=pe(ip(i-iqq,1))
124 pjtm(jt,i)=pm(ip(i-iqq,1))
125 IF(ip(i-iqq,2).EQ.1)
THEN
126 kftj(jt,i+1)=kf(ip(i-iqq,1)+1)
127 pjtx(jt,i+1)=
px(ip(i-iqq,1)+1)
128 pjty(jt,i+1)=
py(ip(i-iqq,1)+1)
129 pjtz(jt,i+1)=
pz(ip(i-iqq,1)+1)
130 pjte(jt,i+1)=pe(ip(i-iqq,1)+1)
131 pjtm(jt,i+1)=pm(ip(i-iqq,1)+1)
136 IF(i.LE.
ntj(jt)) go to 800