author morsch Fri, 10 Sep 2010 07:02:08 +0000 (07:02 +0000) committer morsch Fri, 10 Sep 2010 07:02:08 +0000 (07:02 +0000)

index 47c2474..9900358 100644 (file)
@@ -745,35 +745,10 @@ C       ****** identify the mother particle
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
-                   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
@@ -788,20 +763,11 @@ 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