7 SUBROUTINE HIJFRG(JTP,NTP,IERROR)
8 C NTP=1, fragment proj string, NTP=2, targ string,
10 C strings from jets. JTP is the line number of the string
11 C*******Fragment all leading strings of proj and targ**************
12 C IHNT2(1)=atomic #, IHNT2(2)=proton #(=-1 if anti-proton) *
13 C******************************************************************
15 #include "hiparnt.inc"
17 #include "histrng.inc"
18 #include "hijjet1.inc"
19 #include "hijjet2.inc"
21 #include "lujets_hijing.inc"
22 #include "ludat1_hijing.inc"
28 C ********initialize the document lines
41 C IF(IHPR2(1).GT.0) CALL ATTRAD(IERROR)
42 c IF(IERROR.NE.0) RETURN
43 C CALL LULIST_HIJING(1)
48 IF(NTP.EQ.2) GO TO 200
49 IF(JTP.GT.IHNT2(1)) RETURN
50 IF(NFP(JTP,5).NE.3.AND.NFP(JTP,3).NE.0
51 & .AND.NPJ(JTP).EQ.0.AND.NFP(JTP,10).EQ.0) GO TO 1000
52 IF(NFP(JTP,15).EQ.-1) THEN
71 C ********for NFP(JTP,15)=-1 NFP(JTP,1) IS IN -Z DIRECTION
77 IF((ABS(PB1-PP(JTP,1)).GT.0.01.OR.
78 & ABS(PB2-PP(JTP,2)).GT.0.01).AND.IHPR2(10).NE.0)
79 & WRITE(6,*) ' Pt of Q and QQ do not sum to the total'
83 200 IF(JTP.GT.IHNT2(3)) RETURN
84 IF(NFT(JTP,5).NE.3.AND.NFT(JTP,3).NE.0
85 & .AND.NTJ(JTP).EQ.0.AND.NFT(JTP,10).EQ.0) GO TO 1200
86 IF(NFT(JTP,15).EQ.1) THEN
105 C ********for NFT(JTP,15)=1 NFT(JTP,1) IS IN +Z DIRECTION
112 IF((ABS(PB1-PT(JTP,1)).GT.0.01.OR.
113 & ABS(PB2-PT(JTP,2)).GT.0.01).AND.IHPR2(10).NE.0)
114 & WRITE(6,*) ' Pt of Q and QQ do not sum to the total'
116 300 IF(PECM.LT.HIPR1(1)) THEN
118 IF(IHPR2(10).EQ.0) RETURN
119 WRITE(6,*) ' ECM=',PECM,' energy of the string is too small'
122 AMT=PECM**2+PB1**2+PB2**2
123 AMT1=AM1**2+PQ11**2+PQ12**2
124 AMT2=AM2**2+PQ21**2+PQ22**2
125 PZCM=SQRT(ABS(AMT**2+AMT1**2+AMT2**2-2.0*AMT*AMT1
126 & -2.0*AMT*AMT2-2.0*AMT1*AMT2))/2.0/SQRT(AMT)
127 C *******PZ of end-partons in c.m. frame of the string
133 P(1,4)=SQRT(AMT1+PZCM**2)
140 P(2,4)=SQRT(AMT2+PZCM**2)
144 CALL HIROBO(0.0,0.0,0.0,0.0,BTZ)
146 IF((PQ21**2+PQ22**2).GT.(PQ11**2+PQ12**2)) THEN
159 ELSE IF(NTP.EQ.2) THEN
164 C*******************attach produced jets to the leading partons****
165 IF(NTP.EQ.1.AND.NPJ(JTP).NE.0) THEN
167 C IF(NPJ(JTP).GE.2) CALL HIJSRT(JTP,1)
168 C ********sort jets in order of y
170 IF((ABS(KF1).GT.1000.AND.KF1.LT.0)
171 & .OR.(ABS(KF1).LT.1000.AND.KF1.GT.0)) IEX=1
185 IF(IEX.EQ.1) I0=NPJ(JTP)-I+1
186 C ********reverse the order of jets
190 IF(KK1.NE.21 .AND. KK1.NE.0) K(I+1,1)=
191 & 1+(ABS(KK1)+(2*IEX-1)*KK1)/2/ABS(KK1)
192 P(I+1,1)=PJPX(JTP,I0)
193 P(I+1,2)=PJPY(JTP,I0)
194 P(I+1,3)=PJPZ(JTP,I0)
195 P(I+1,4)=PJPE(JTP,I0)
196 P(I+1,5)=PJPM(JTP,I0)
199 ELSE IF(NTP.EQ.2.AND.NTJ(JTP).NE.0) THEN
201 c IF(NTJ(JTP).GE.2) CALL HIJSRT(JTP,2)
202 C ********sort jets in order of y
204 IF((ABS(KF2).GT.1000.AND.KF2.LT.0)
205 & .OR.(ABS(KF2).LT.1000.AND.KF2.GT.0)) IEX=0
219 IF(IEX.EQ.1) I0=NTJ(JTP)-I+1
220 C ********reverse the order of jets
224 IF(KK1.NE.21 .AND. KK1.NE.0) K(I+1,1)=
225 & 1+(ABS(KK1)+(2*IEX-1)*KK1)/2/ABS(KK1)
226 P(I+1,1)=PJTX(JTP,I0)
227 P(I+1,2)=PJTY(JTP,I0)
228 P(I+1,3)=PJTZ(JTP,I0)
229 P(I+1,4)=PJTE(JTP,I0)
230 P(I+1,5)=PJTM(JTP,I0)
234 IF(IHPR2(1).GT.0.AND.RLU_HIJING(0).LE.HIDAT(3)) THEN
237 IF(IHPR2(8).EQ.0.AND.IHPR2(3).EQ.0.AND.IHPR2(9).EQ.0)
239 IF(HINT1(1).GE.1000.0.AND.JETOT.EQ.0)THEN
246 ELSE IF(JETOT.EQ.0.AND.IHPR2(1).GT.0.AND.
247 & HINT1(1).GE.1000.0.AND.
248 & RLU_HIJING(0).LE.0.8) THEN
253 IF(IHPR2(8).EQ.0.AND.IHPR2(3).EQ.0.AND.IHPR2(9).EQ.0)
259 IF(IERROR.NE.0) RETURN
260 C ******** conduct soft radiations
261 C****************************
264 C CALL LULIST_HIJING(1)
274 C ********proj remain as a nucleon or delta
276 C call LULIST_HIJING(1)
285 C ********targ remain as a nucleon or delta
287 C call LULIST_HIJING(1)