5 C********************************************************************
6 C Sort the jets associated with a nucleon in order of their
8 C********************************************************************
9 SUBROUTINE HIJSRT(JPJT,NPT)
10 DIMENSION KF(100),PX(100),PY(100),PZ(100),PE(100),PM(100)
11 DIMENSION Y(100),IP(100,2)
13 #include "hijjet1.inc"
15 IF(NPT.EQ.2) GO TO 500
25 Y(I-IQ)=0.5*ALOG((ABS(PE(I)+PZ(I))+1.E-5)
26 & /(ABS(PE(I)-PZ(I))+1.E-5))
41 IF(I.LE.NPJ(JP)) GO TO 100
44 DO 200 J=I+1,NPJ(JP)-IQ
45 IF(Y(I).GT.Y(J)) GO TO 200
53 C ********sort in decending y
56 300 KFPJ(JP,I)=KF(IP(I-IQQ,1))
57 PJPX(JP,I)=PX(IP(I-IQQ,1))
58 PJPY(JP,I)=PY(IP(I-IQQ,1))
59 PJPZ(JP,I)=PZ(IP(I-IQQ,1))
60 PJPE(JP,I)=PE(IP(I-IQQ,1))
61 PJPM(JP,I)=PM(IP(I-IQQ,1))
62 IF(IP(I-IQQ,2).EQ.1) THEN
63 KFPJ(JP,I+1)=KF(IP(I-IQQ,1)+1)
64 PJPX(JP,I+1)=PX(IP(I-IQQ,1)+1)
65 PJPY(JP,I+1)=PY(IP(I-IQQ,1)+1)
66 PJPZ(JP,I+1)=PZ(IP(I-IQQ,1)+1)
67 PJPE(JP,I+1)=PE(IP(I-IQQ,1)+1)
68 PJPM(JP,I+1)=PM(IP(I-IQQ,1)+1)
73 IF(I.LE.NPJ(JP)) GO TO 300
86 Y(I-IQ)=0.5*ALOG((ABS(PE(I)+PZ(I))+1.E-5)
87 & /(ABS(PE(I)-PZ(I))+1.E-5))
102 IF(I.LE.NTJ(JT)) GO TO 600
104 DO 700 I=1,NTJ(JT)-IQ
105 DO 700 J=I+1,NTJ(JT)-IQ
106 IF(Y(I).LT.Y(J)) GO TO 700
114 C ********sort in acending y
117 800 KFTJ(JT,I)=KF(IP(I-IQQ,1))
118 PJTX(JT,I)=PX(IP(I-IQQ,1))
119 PJTY(JT,I)=PY(IP(I-IQQ,1))
120 PJTZ(JT,I)=PZ(IP(I-IQQ,1))
121 PJTE(JT,I)=PE(IP(I-IQQ,1))
122 PJTM(JT,I)=PM(IP(I-IQQ,1))
123 IF(IP(I-IQQ,2).EQ.1) THEN
124 KFTJ(JT,I+1)=KF(IP(I-IQQ,1)+1)
125 PJTX(JT,I+1)=PX(IP(I-IQQ,1)+1)
126 PJTY(JT,I+1)=PY(IP(I-IQQ,1)+1)
127 PJTZ(JT,I+1)=PZ(IP(I-IQQ,1)+1)
128 PJTE(JT,I+1)=PE(IP(I-IQQ,1)+1)
129 PJTM(JT,I+1)=PM(IP(I-IQQ,1)+1)
134 IF(I.LE.NTJ(JT)) GO TO 800