[u/mrichter/AliRoot.git] / HIJING / hijing1_36 / hijfrg.F

* $Id$

C
C
C
C
	SUBROUTINE HIJFRG(JTP,NTP,IERROR)
C	NTP=1, fragment proj string, NTP=2, targ string, 
C       NTP=3, independent 
C	strings from jets.  JTP is the line number of the string
C*******Fragment all leading strings of proj and targ**************
C	IHNT2(1)=atomic #, IHNT2(2)=proton #(=-1 if anti-proton)  *
C******************************************************************
#define BLANKET_SAVE
#include "hiparnt.inc"
#include "hijdat.inc"
#include "histrng.inc"
#include "hijjet1.inc"
#include "hijjet2.inc"
C
#include "lujets_hijing.inc"
#include "ludat1_hijing.inc"
	SAVE
	
	IERROR=0
	CALL LUEDIT_HIJING(0)
	N=0
C			********initialize the document lines
	IF(NTP.EQ.3) THEN
		ISG=JTP
		N=NJSG(ISG)
		DO 100 I=1,NJSG(ISG)
			K(I,1)=K1SG(ISG,I)
			K(I,2)=K2SG(ISG,I)
			P(I,1)=PXSG(ISG,I)
			P(I,2)=PYSG(ISG,I)
			P(I,3)=PZSG(ISG,I)
			P(I,4)=PESG(ISG,I)
			P(I,5)=PMSG(ISG,I)
100		CONTINUE
C		IF(IHPR2(1).GT.0) CALL ATTRAD(IERROR)
c		IF(IERROR.NE.0) RETURN
C		CALL LULIST_HIJING(1)
		CALL LUEXEC_HIJING
		RETURN
	ENDIF
C
	IF(NTP.EQ.2) GO TO 200
	IF(JTP.GT.IHNT2(1))   RETURN
	IF(NFP(JTP,5).NE.3.AND.NFP(JTP,3).NE.0
     &	    .AND.NPJ(JTP).EQ.0.AND.NFP(JTP,10).EQ.0) GO TO 1000
	IF(NFP(JTP,15).EQ.-1) THEN
		KF1=NFP(JTP,2)
		KF2=NFP(JTP,1)
		PQ21=PP(JTP,6)
		PQ22=PP(JTP,7)
		PQ11=PP(JTP,8)
		PQ12=PP(JTP,9)
		AM1=PP(JTP,15)
		AM2=PP(JTP,14)
	ELSE
		KF1=NFP(JTP,1)
		KF2=NFP(JTP,2)
		PQ21=PP(JTP,8)
		PQ22=PP(JTP,9)
		PQ11=PP(JTP,6)
		PQ12=PP(JTP,7)
		AM1=PP(JTP,14)
		AM2=PP(JTP,15)	
	ENDIF
C	********for NFP(JTP,15)=-1 NFP(JTP,1) IS IN -Z DIRECTION
	PB1=PQ11+PQ21
	PB2=PQ12+PQ22
	PB3=PP(JTP,3)
	PECM=PP(JTP,5)
	BTZ=PB3/PP(JTP,4)
	IF((ABS(PB1-PP(JTP,1)).GT.0.01.OR.
     &    ABS(PB2-PP(JTP,2)).GT.0.01).AND.IHPR2(10).NE.0)
     &	  WRITE(6,*) '  Pt of Q and QQ do not sum to the total'

	GO TO 300

200	IF(JTP.GT.IHNT2(3))  RETURN
	IF(NFT(JTP,5).NE.3.AND.NFT(JTP,3).NE.0
     &	   .AND.NTJ(JTP).EQ.0.AND.NFT(JTP,10).EQ.0) GO TO 1200
	IF(NFT(JTP,15).EQ.1) THEN
		KF1=NFT(JTP,1)
		KF2=NFT(JTP,2)
		PQ11=PT(JTP,6)
		PQ12=PT(JTP,7)
		PQ21=PT(JTP,8)
		PQ22=PT(JTP,9)
		AM1=PT(JTP,14)
		AM2=PT(JTP,15)
	ELSE
		KF1=NFT(JTP,2)
		KF2=NFT(JTP,1)
		PQ11=PT(JTP,8)
		PQ12=PT(JTP,9)
		PQ21=PT(JTP,6)
		PQ22=PT(JTP,7)
		AM1=PT(JTP,15)
		AM2=PT(JTP,14)
	ENDIF	
C	********for NFT(JTP,15)=1 NFT(JTP,1) IS IN +Z DIRECTION
	PB1=PQ11+PQ21
	PB2=PQ12+PQ22
	PB3=PT(JTP,3)
	PECM=PT(JTP,5)
	BTZ=PB3/PT(JTP,4)

	IF((ABS(PB1-PT(JTP,1)).GT.0.01.OR.
     &     ABS(PB2-PT(JTP,2)).GT.0.01).AND.IHPR2(10).NE.0)
     &     WRITE(6,*) '  Pt of Q and QQ do not sum to the total'

300	IF(PECM.LT.HIPR1(1)) THEN
	   IERROR=1
	   IF(IHPR2(10).EQ.0) RETURN
	   WRITE(6,*) ' ECM=',PECM,' energy of the string is too small'
	   RETURN
	ENDIF
	AMT=PECM**2+PB1**2+PB2**2
	AMT1=AM1**2+PQ11**2+PQ12**2
	AMT2=AM2**2+PQ21**2+PQ22**2
	PZCM=SQRT(ABS(AMT**2+AMT1**2+AMT2**2-2.0*AMT*AMT1
     &       -2.0*AMT*AMT2-2.0*AMT1*AMT2))/2.0/SQRT(AMT)
C		*******PZ of end-partons in c.m. frame of the string
	K(1,1)=2
	K(1,2)=KF1
	P(1,1)=PQ11
	P(1,2)=PQ12
	P(1,3)=PZCM
	P(1,4)=SQRT(AMT1+PZCM**2)
	P(1,5)=AM1
	K(2,1)=1
	K(2,2)=KF2
	P(2,1)=PQ21
	P(2,2)=PQ22
	P(2,3)=-PZCM
	P(2,4)=SQRT(AMT2+PZCM**2)
	P(2,5)=AM2
	N=2
C*****
	CALL HIROBO(0.0,0.0,0.0,0.0,BTZ)
	JETOT=0
	IF((PQ21**2+PQ22**2).GT.(PQ11**2+PQ12**2)) THEN
		PMAX1=P(2,1)
		PMAX2=P(2,2)
		PMAX3=P(2,3)
	ELSE
		PMAX1=P(1,1)
		PMAX2=P(1,2)
		PMAX3=P(1,3)
	ENDIF
	IF(NTP.EQ.1) THEN
		PP(JTP,10)=PMAX1
		PP(JTP,11)=PMAX2
		PP(JTP,12)=PMAX3
	ELSE IF(NTP.EQ.2) THEN
		PT(JTP,10)=PMAX1
		PT(JTP,11)=PMAX2
		PT(JTP,12)=PMAX3
	ENDIF
C*******************attach produced jets to the leading partons****
	IF(NTP.EQ.1.AND.NPJ(JTP).NE.0) THEN
		JETOT=NPJ(JTP)
C		IF(NPJ(JTP).GE.2) CALL HIJSRT(JTP,1)
C			********sort jets in order of y
		IEX=0
		IF((ABS(KF1).GT.1000.AND.KF1.LT.0)
     &			.OR.(ABS(KF1).LT.1000.AND.KF1.GT.0)) IEX=1
		DO 520 I=N,2,-1
		DO 520 J=1,5
			II=NPJ(JTP)+I
			K(II,J)=K(I,J)
			P(II,J)=P(I,J)
			V(II,J)=V(I,J)
520		CONTINUE
		DO 540 I=1,NPJ(JTP)
			DO 542 J=1,5
				K(I+1,J)=0
				V(I+1,J)=0
542			CONTINUE				
			I0=I
			IF(IEX.EQ.1) I0=NPJ(JTP)-I+1
C				********reverse the order of jets
			KK1=KFPJ(JTP,I0)
			K(I+1,1)=2
			K(I+1,2)=KK1
			IF(KK1.NE.21 .AND. KK1.NE.0)  K(I+1,1)=
     &			  1+(ABS(KK1)+(2*IEX-1)*KK1)/2/ABS(KK1)
			P(I+1,1)=PJPX(JTP,I0)
			P(I+1,2)=PJPY(JTP,I0)
			P(I+1,3)=PJPZ(JTP,I0)
			P(I+1,4)=PJPE(JTP,I0)
			P(I+1,5)=PJPM(JTP,I0)
540		CONTINUE
		N=N+NPJ(JTP)
	ELSE IF(NTP.EQ.2.AND.NTJ(JTP).NE.0) THEN
		JETOT=NTJ(JTP)
c		IF(NTJ(JTP).GE.2)  CALL HIJSRT(JTP,2)
C			********sort jets in order of y
		IEX=1
		IF((ABS(KF2).GT.1000.AND.KF2.LT.0)
     &			.OR.(ABS(KF2).LT.1000.AND.KF2.GT.0)) IEX=0
		DO 560 I=N,2,-1
		DO 560 J=1,5
			II=NTJ(JTP)+I
			K(II,J)=K(I,J)
			P(II,J)=P(I,J)
			V(II,J)=V(I,J)
560		CONTINUE
		DO 580 I=1,NTJ(JTP)
			DO 582 J=1,5
				K(I+1,J)=0
				V(I+1,J)=0
582			CONTINUE				
			I0=I
			IF(IEX.EQ.1) I0=NTJ(JTP)-I+1
C				********reverse the order of jets
			KK1=KFTJ(JTP,I0)
			K(I+1,1)=2
			K(I+1,2)=KK1
			IF(KK1.NE.21 .AND. KK1.NE.0) K(I+1,1)=
     &		           1+(ABS(KK1)+(2*IEX-1)*KK1)/2/ABS(KK1)
			P(I+1,1)=PJTX(JTP,I0)
			P(I+1,2)=PJTY(JTP,I0)
			P(I+1,3)=PJTZ(JTP,I0)
			P(I+1,4)=PJTE(JTP,I0)
			P(I+1,5)=PJTM(JTP,I0)
580		CONTINUE
		N=N+NTJ(JTP)
	ENDIF
	IF(IHPR2(1).GT.0.AND.RLU_HIJING(0).LE.HIDAT(3)) THEN
	     HIDAT20=HIDAT(2)
	     HIPR150=HIPR1(5)
	     IF(IHPR2(8).EQ.0.AND.IHPR2(3).EQ.0.AND.IHPR2(9).EQ.0)
     &			HIDAT(2)=2.0
	     IF(HINT1(1).GE.1000.0.AND.JETOT.EQ.0)THEN
		HIDAT(2)=3.0
		HIPR1(5)=5.0
	     ENDIF
	     CALL ATTRAD(IERROR)
	     HIDAT(2)=HIDAT20
	     HIPR1(5)=HIPR150
	ELSE IF(JETOT.EQ.0.AND.IHPR2(1).GT.0.AND.
     &                       HINT1(1).GE.1000.0.AND.
     &		RLU_HIJING(0).LE.0.8) THEN
		HIDAT20=HIDAT(2)
		HIPR150=HIPR1(5)
		HIDAT(2)=3.0
		HIPR1(5)=5.0
	     IF(IHPR2(8).EQ.0.AND.IHPR2(3).EQ.0.AND.IHPR2(9).EQ.0)
     &			HIDAT(2)=2.0
		CALL ATTRAD(IERROR)
		HIDAT(2)=HIDAT20
		HIPR1(5)=HIPR150
	ENDIF
	IF(IERROR.NE.0) RETURN
C		******** conduct soft radiations
C****************************
C
C
C	CALL LULIST_HIJING(1)
	CALL LUEXEC_HIJING
	RETURN

1000	N=1
	K(1,1)=1
       	K(1,2)=NFP(JTP,3)
	DO 1100 JJ=1,5
       		P(1,JJ)=PP(JTP,JJ)
1100		CONTINUE
C			********proj remain as a nucleon or delta
	CALL LUEXEC_HIJING
C	call LULIST_HIJING(1)
	RETURN
C
1200	N=1
	K(1,1)=1
	K(1,2)=NFT(JTP,3)
	DO 1300 JJ=1,5
		P(1,JJ)=PT(JTP,JJ)
1300	CONTINUE
C			********targ remain as a nucleon or delta
	CALL LUEXEC_HIJING
C	call LULIST_HIJING(1)
	RETURN
	END
Commit	Line	Data
e74335a4	1	* $Id$
	2
	3	C
	4	C
	5	C
	6	C
	7	SUBROUTINE HIJFRG(JTP,NTP,IERROR)
	8	C NTP=1, fragment proj string, NTP=2, targ string,
	9	C NTP=3, independent
	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******************************************************************
bc676b8e	14	#define BLANKET_SAVE
e74335a4	15	#include "hiparnt.inc"
	16	#include "hijdat.inc"
	17	#include "histrng.inc"
	18	#include "hijjet1.inc"
	19	#include "hijjet2.inc"
	20	C
	21	#include "lujets_hijing.inc"
	22	#include "ludat1_hijing.inc"
	23	SAVE
	24
	25	IERROR=0
	26	CALL LUEDIT_HIJING(0)
	27	N=0
	28	C ********initialize the document lines
	29	IF(NTP.EQ.3) THEN
	30	ISG=JTP
	31	N=NJSG(ISG)
	32	DO 100 I=1,NJSG(ISG)
	33	K(I,1)=K1SG(ISG,I)
	34	K(I,2)=K2SG(ISG,I)
	35	P(I,1)=PXSG(ISG,I)
	36	P(I,2)=PYSG(ISG,I)
	37	P(I,3)=PZSG(ISG,I)
	38	P(I,4)=PESG(ISG,I)
	39	P(I,5)=PMSG(ISG,I)
	40	100 CONTINUE
	41	C IF(IHPR2(1).GT.0) CALL ATTRAD(IERROR)
	42	c IF(IERROR.NE.0) RETURN
	43	C CALL LULIST_HIJING(1)
	44	CALL LUEXEC_HIJING
	45	RETURN
	46	ENDIF
	47	C
	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
	53	KF1=NFP(JTP,2)
	54	KF2=NFP(JTP,1)
	55	PQ21=PP(JTP,6)
	56	PQ22=PP(JTP,7)
	57	PQ11=PP(JTP,8)
	58	PQ12=PP(JTP,9)
	59	AM1=PP(JTP,15)
	60	AM2=PP(JTP,14)
	61	ELSE
	62	KF1=NFP(JTP,1)
	63	KF2=NFP(JTP,2)
	64	PQ21=PP(JTP,8)
	65	PQ22=PP(JTP,9)
	66	PQ11=PP(JTP,6)
	67	PQ12=PP(JTP,7)
	68	AM1=PP(JTP,14)
	69	AM2=PP(JTP,15)
	70	ENDIF
	71	C ********for NFP(JTP,15)=-1 NFP(JTP,1) IS IN -Z DIRECTION
	72	PB1=PQ11+PQ21
	73	PB2=PQ12+PQ22
	74	PB3=PP(JTP,3)
	75	PECM=PP(JTP,5)
	76	BTZ=PB3/PP(JTP,4)
	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'
80
81	GO TO 300
82
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
87	KF1=NFT(JTP,1)
88	KF2=NFT(JTP,2)
89	PQ11=PT(JTP,6)
90	PQ12=PT(JTP,7)
91	PQ21=PT(JTP,8)
92	PQ22=PT(JTP,9)
93	AM1=PT(JTP,14)
94	AM2=PT(JTP,15)
95	ELSE
96	KF1=NFT(JTP,2)
97	KF2=NFT(JTP,1)
98	PQ11=PT(JTP,8)
99	PQ12=PT(JTP,9)
100	PQ21=PT(JTP,6)
101	PQ22=PT(JTP,7)
102	AM1=PT(JTP,15)
103	AM2=PT(JTP,14)
104	ENDIF
105	C ********for NFT(JTP,15)=1 NFT(JTP,1) IS IN +Z DIRECTION
106	PB1=PQ11+PQ21
107	PB2=PQ12+PQ22
108	PB3=PT(JTP,3)
109	PECM=PT(JTP,5)
110	BTZ=PB3/PT(JTP,4)
111
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'
115
116	300 IF(PECM.LT.HIPR1(1)) THEN
117	IERROR=1
118	IF(IHPR2(10).EQ.0) RETURN
119	WRITE(6,*) ' ECM=',PECM,' energy of the string is too small'
120	RETURN
121	ENDIF
122	AMT=PECM2+PB12+PB2**2
123	AMT1=AM12+PQ112+PQ12**2
124	AMT2=AM22+PQ212+PQ22**2
125	PZCM=SQRT(ABS(AMT2+AMT12+AMT2*2-2.0AMT*AMT1
126	& -2.0AMTAMT2-2.0AMT1AMT2))/2.0/SQRT(AMT)
127	C *******PZ of end-partons in c.m. frame of the string
128	K(1,1)=2
129	K(1,2)=KF1
130	P(1,1)=PQ11
131	P(1,2)=PQ12
132	P(1,3)=PZCM
133	P(1,4)=SQRT(AMT1+PZCM**2)
134	P(1,5)=AM1
135	K(2,1)=1
136	K(2,2)=KF2
137	P(2,1)=PQ21
138	P(2,2)=PQ22
139	P(2,3)=-PZCM
140	P(2,4)=SQRT(AMT2+PZCM**2)
141	P(2,5)=AM2
142	N=2
143	C*****
144	CALL HIROBO(0.0,0.0,0.0,0.0,BTZ)
145	JETOT=0
146	IF((PQ212+PQ222).GT.(PQ112+PQ122)) THEN
147	PMAX1=P(2,1)
148	PMAX2=P(2,2)
149	PMAX3=P(2,3)
150	ELSE
151	PMAX1=P(1,1)
152	PMAX2=P(1,2)
153	PMAX3=P(1,3)
154	ENDIF
155	IF(NTP.EQ.1) THEN
156	PP(JTP,10)=PMAX1
157	PP(JTP,11)=PMAX2
158	PP(JTP,12)=PMAX3
159	ELSE IF(NTP.EQ.2) THEN
160	PT(JTP,10)=PMAX1
161	PT(JTP,11)=PMAX2
162	PT(JTP,12)=PMAX3
163	ENDIF
164	C*****************attach produced jets to the leading partons**
165	IF(NTP.EQ.1.AND.NPJ(JTP).NE.0) THEN
166	JETOT=NPJ(JTP)
167	C IF(NPJ(JTP).GE.2) CALL HIJSRT(JTP,1)
168	C ********sort jets in order of y
169	IEX=0
170	IF((ABS(KF1).GT.1000.AND.KF1.LT.0)
171	& .OR.(ABS(KF1).LT.1000.AND.KF1.GT.0)) IEX=1
172	DO 520 I=N,2,-1
173	DO 520 J=1,5
174	II=NPJ(JTP)+I
175	K(II,J)=K(I,J)
176	P(II,J)=P(I,J)
177	V(II,J)=V(I,J)
178	520 CONTINUE
179	DO 540 I=1,NPJ(JTP)
180	DO 542 J=1,5
181	K(I+1,J)=0
182	V(I+1,J)=0
183	542 CONTINUE
184	I0=I
185	IF(IEX.EQ.1) I0=NPJ(JTP)-I+1
186	C ********reverse the order of jets
187	KK1=KFPJ(JTP,I0)
188	K(I+1,1)=2
189	K(I+1,2)=KK1
190	IF(KK1.NE.21 .AND. KK1.NE.0) K(I+1,1)=
191	& 1+(ABS(KK1)+(2IEX-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)
197	540 CONTINUE
198	N=N+NPJ(JTP)
199	ELSE IF(NTP.EQ.2.AND.NTJ(JTP).NE.0) THEN
200	JETOT=NTJ(JTP)
201	c IF(NTJ(JTP).GE.2) CALL HIJSRT(JTP,2)
202	C ********sort jets in order of y
203	IEX=1
204	IF((ABS(KF2).GT.1000.AND.KF2.LT.0)
205	& .OR.(ABS(KF2).LT.1000.AND.KF2.GT.0)) IEX=0
206	DO 560 I=N,2,-1
207	DO 560 J=1,5
208	II=NTJ(JTP)+I
209	K(II,J)=K(I,J)
210	P(II,J)=P(I,J)
211	V(II,J)=V(I,J)
212	560 CONTINUE
213	DO 580 I=1,NTJ(JTP)
214	DO 582 J=1,5
215	K(I+1,J)=0
216	V(I+1,J)=0
217	582 CONTINUE
218	I0=I
219	IF(IEX.EQ.1) I0=NTJ(JTP)-I+1
220	C ********reverse the order of jets
221	KK1=KFTJ(JTP,I0)
222	K(I+1,1)=2
223	K(I+1,2)=KK1
224	IF(KK1.NE.21 .AND. KK1.NE.0) K(I+1,1)=
225	& 1+(ABS(KK1)+(2IEX-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)
231	580 CONTINUE
232	N=N+NTJ(JTP)
233	ENDIF
234	IF(IHPR2(1).GT.0.AND.RLU_HIJING(0).LE.HIDAT(3)) THEN
235	HIDAT20=HIDAT(2)
236	HIPR150=HIPR1(5)
237	IF(IHPR2(8).EQ.0.AND.IHPR2(3).EQ.0.AND.IHPR2(9).EQ.0)
238	& HIDAT(2)=2.0
239	IF(HINT1(1).GE.1000.0.AND.JETOT.EQ.0)THEN
240	HIDAT(2)=3.0
241	HIPR1(5)=5.0
242	ENDIF
243	CALL ATTRAD(IERROR)
244	HIDAT(2)=HIDAT20
245	HIPR1(5)=HIPR150
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
249	HIDAT20=HIDAT(2)
250	HIPR150=HIPR1(5)
251	HIDAT(2)=3.0
252	HIPR1(5)=5.0
253	IF(IHPR2(8).EQ.0.AND.IHPR2(3).EQ.0.AND.IHPR2(9).EQ.0)
254	& HIDAT(2)=2.0
255	CALL ATTRAD(IERROR)
256	HIDAT(2)=HIDAT20
257	HIPR1(5)=HIPR150
258	ENDIF
259	IF(IERROR.NE.0) RETURN
260	C ******** conduct soft radiations
261	C****************************
262	C
263	C
264	C CALL LULIST_HIJING(1)
265	CALL LUEXEC_HIJING
266	RETURN
267
268	1000 N=1
269	K(1,1)=1
270	K(1,2)=NFP(JTP,3)
271	DO 1100 JJ=1,5
272	P(1,JJ)=PP(JTP,JJ)
273	1100 CONTINUE
274	C ********proj remain as a nucleon or delta
275	CALL LUEXEC_HIJING
276	C call LULIST_HIJING(1)
277	RETURN
278	C
279	1200 N=1
280	K(1,1)=1
281	K(1,2)=NFT(JTP,3)
282	DO 1300 JJ=1,5
283	P(1,JJ)=PT(JTP,JJ)
284	1300 CONTINUE
285	C ********targ remain as a nucleon or delta
286	CALL LUEXEC_HIJING
287	C call LULIST_HIJING(1)
288	RETURN
289	END