formation time not added to vertex position
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 17 Nov 2010 00:31:17 +0000 (00:31 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 17 Nov 2010 00:31:17 +0000 (00:31 +0000)
HIJING/hijing1_36/hijing.F

index 9900358..7da5a65 100644 (file)
@@ -601,8 +601,6 @@ C********N_STR+1 is the number of strings in fragmentation
 C********the number of strings before a line is stored in K(I,4)
 C********IDSTR is id number of the string system (91,92 or 93)
 C
-        GAMMA0=HINT1(1)/2.0/0.93827
-        BETA0=sqrt(GAMMA0**2-1.0)/GAMMA0
 c
         IF(IHPR2(20).NE.0) THEN
           DO 360 ISG=1,NSG
@@ -659,31 +657,12 @@ C       ****** identify the mother particle
                   PATT(NATT,2)=P(I,2)
                   PATT(NATT,3)=P(I,3)
                   PATT(NATT,4)=P(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+P(I,4)
-                   VATT01=0.5*(YP(1,IASG(ISG,1))+YT(1,IASG(ISG,2)))
-                   VATT02=0.5*(YP(2,IASG(ISG,1))+YT(2,IASG(ISG,2)))
-                   VATT03=0.5*(YP(3,IASG(ISG,1))+YT(3,IASG(ISG,2)))
-     &                  /GAMMA0
-                   VATT04=-0.5*(YP(3,IASG(ISG,1))
-     &                  -YT(3,IASG(ISG,2)))/GAMMA0/BETA0
-                   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
 360       CONTINUE
 C              ********Fragment the q-qbar jets systems *****
 C