Initialize and fill VATT array for vertex positions.

diff --git a/HIJING/hijing1_36/hidata.F b/HIJING/hijing1_36/hidata.F
index ad614bd5a1e0778959ad4cf8847c1aac43fedc7c..0bfd1d25034ee108e8b8015804adc40a2581615c 100644 (file)
--- a/HIJING/hijing1_36/hidata.F
+++ b/HIJING/hijing1_36/hidata.F
@@ -37,6 +37,7 @@ C...give all the switchs and parameters the default values
 C...initialize all the data common blocks
        DATA NATT/0/EATT/0.0/JATT/0/NT/0/NP/0/N0/0/N01/0/N10/0/N11/0/
        DATA KATT/800000*0/PATT/800000*0.0/
+       DATA VATT/800000*0/
 
        DATA NFP/4500*0/PP/4500*0.0/NFT/4500*0/PT/4500*0.0/
 

diff --git a/HIJING/hijing1_36/hijing.F b/HIJING/hijing1_36/hijing.F
index b904fdd206a49ce52c06348cfdb6087061830313..7274c64d3cb09b7d1b1627422e33705ba6ef8250 100644 (file)
--- a/HIJING/hijing1_36/hijing.F
+++ b/HIJING/hijing1_36/hijing.F
@@ -571,6 +571,11 @@ 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)
 360       CONTINUE
 C              ********Fragment the q-qbar jets systems *****
@@ -629,6 +634,11 @@ 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)
 390            CONTINUE 
 400       CONTINUE
@@ -644,6 +654,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)
 450    CONTINUE
 C                      ********store the direct-produced particles