Subroutine Get_VTX(Vtx,Ierr) C----------------------------------------------------------------------- C C Get and store event vertex from VCTVTX table in NTuple. C Vertex is set to (x,y,z)=(0.,0.,0.) if no vertex is found. C C Input: NONE C Output: Ierr = 0 ok C = 1 no vertex found, set to 0/0/0 C----------------------------------------------------------------------- Implicit NONE #include "partap.inc" #include "vctvtx.inc" #include "vcevtctd.inc" Integer Iprint,Ierr,IerrVCget Real Vtx Dimension Vtx(3) Data Iprint/0/ Save Iprint C... Check whether there was a vertex at all. If yes, store vertex: Ierr = 0 if (COUTAB(VCEVTCTD).ne.0) then call VCGetctd(IerrVCget) endif If ( Coutab(VCTVTX) .eq. 0 ) Then Ierr = 1 Vtx(1) = 0. Vtx(2) = 0. Vtx(3) = 0. Return EndIf Call Fettab(VCTVTX,ID,1) Vtx(1) = VCTVTX_V(1) Vtx(2) = VCTVTX_V(2) Vtx(3) = VCTVTX_V(3) 999 Continue Return End