PATT(NATT,3)=P(I,3)
PATT(NATT,4)=P(I,4)
EATT=EATT+P(I,4)
-C ******* the following calculate the formation time and position
-C of the produced hadrons
- IF(NTP.EQ.1) THEN
- VATT01=YP(1,J_JTP)
- VATT02=YP(2,J_JTP)
- VATT03=YP(3,J_JTP)/GAMMA0
- ELSE IF(NTP.EQ.2) THEN
- VATT01=YT(1,J_JTP)
- VATT02=YT(2,J_JTP)
- VATT03=YT(3,J_JTP)/GAMMA0
- ENDIF
- RARB=1.12*(IHNT2(1)**0.33333+IHNT2(3)**0.33333)
- V3MIN1=RARB/GAMMA0
- V3MIN2=1.0/MAX(1.0,5.08*ABS(PATT(I,3)))
- VATT_MIN=MAX(V3MIN1,V3MIN2)
- VATT03=VATT03+(0.5-RLU_HIJING(0))*VATT_MIN
- amt2=P(I,1)**2+P(I,2)**2+P(I,5)**2
- IF(amt2.GT.0.0) THEN
- tauf=0.2*2.0*P(I,3)/amt2
- VATT(NATT,1)=VATT01
- VATT(NATT,2)=VATT02
- VATT(NATT,3)=VATT03+tauf
- VATT(NATT,4)=abs(VATT(NATT,3))
- ELSE
- VATT(NATT,4)=abs(VATT03)
- VATT(NATT,1)=VATT01
- VATT(NATT,2)=VATT02
- VATT(NATT,3)=VATT03
- ENDIF
+ VATT(NATT,1)=V(I,1)
+ VATT(NATT,2)=V(I,2)
+ VATT(NATT,3)=V(I,3)
+ VATT(NATT,4)=V(I,4)
390 CONTINUE
400 CONTINUE
C ********Fragment the q-qq related string systems
PATT(NATT,2)=PDR(I,2)
PATT(NATT,3)=PDR(I,3)
PATT(NATT,4)=PDR(I,4)
+ VATT(NATT,1)=V(I,1)
+ VATT(NATT,2)=V(I,2)
+ VATT(NATT,3)=V(I,3)
+ VATT(NATT,4)=V(I,4)
EATT=EATT+PDR(I,4)
- amt2=P(I,1)**2+P(I,2)**2+P(I,5)**2
- IF(amt2.GT.0.0) THEN
- tauf=0.2*2.0*P(I,3)/amt2
- VATT(NATT,1)=0.0
- VATT(NATT,2)=0.0
- VATT(NATT,3)=tauf
- VATT(NATT,4)=abs(VATT(NATT,3))
- ELSE
- VATT(NATT,4)=0.0
- VATT(NATT,1)=0.0
- VATT(NATT,2)=0.0
- VATT(NATT,3)=0.0
- ENDIF
450 CONTINUE
C ********store the direct-produced particles
C