c   create geant ntuples - column wise ntuple as per hbook manual p 25-30
      subroutine GUNTUPLE
c
      CHARACTER*8 CHTAGS(17),dirntup,nttitle
	character*80 namefil
	integer ntupid,ISTAT,lunnum
        real xvert,yvert,zvert,eex,eey,eez,
     1 ee,ppx,ppy,ppz,pproton,xe7,ye7,ze7,xp8,yp8,zp8,ppxv,ppyv,ppzv,
     1 ppv,eexv,eeyv,eezv,eev,xrast0,yrast0
      common/ntinfo/dirntup,ntupid,lunnum,xvert,yvert,zvert,eex,eey,eez,
     1 ee,ppx,ppy,ppz,pproton,xe7,ye7,ze7,xp8,yp8,zp8,ppxv,ppyv,ppzv,
     1 ppv,eexv,eeyv,eezv,eev,xrast0,yrast0

c e-p coincidence starts at vertex point. Electron detected in sp07, proton detected in sp08
c eepv(1) = xvert             vertex x
c eepv(2) = yvert
c eepv(3) = zvert
c eepv(4) = eex = px/ptot  direction cosines for electron in sp07
c eepv(5) = eey = py/ptot
c eepv(6) = eez = pz/ptot
c eepv(7) = ee  = ptot     total 3 momentum of electron in sp07
c eepv(8) = ppx = px/ptot  direction cosines for proton in sp08
c eepv(9) = ppy = py/ptot  
c eepv(10) = ppz = pz/ptot 
c eepv(11) = ptot =pproton  total 3 momentum of proton in sp08
c eepv(12) = xe7            electron x coordinate in sp07
c eepv(13) = ye7
c eepv(14) = ze7
c eepv(15) = xp8            proton x coordinate in sp08
c eepv(16) = yp8            
c eepv(17) = zp8
c eepv(18) = ppxv  proton x momentum at the vertex in target
c eepv(19) = ppyv    "    y   "
c eepv(20) = ppzv    "    z   "
c eepv(21) = ppv     "   total "
c eepv(22) = eex   electron x momentum at the vertex in the target
c eepv(23) = eey       "    y     "
c eepv(24) = eez       "    z     "
c eepv(25) = ee        " total    "
c eepv(26) = xrast  from SP01
c eepv(27) = yrast from SP01
*
*     ------------------------------------------------------------------
*
*             Open a CW NTUPLE direct access file
*
      dirntup = 'NTUPEL'
      namefil = 'eepnt.hbook'
      nttitle = 'eep coincidence'
      ntupid = 20
      lunnum=76
      call hropen(lunnum,dirntup,namefil,'N',1024,ISTAT)
      call hbnt(ntupid,nttitle,' ')
      call hbname(ntupid,'eepcoin',xvert,'xvert,yvert,zvert,eex,eey,
     1 eez,ee,ppx,ppy,ppz,pproton,xe7,ye7,ze7,xp8,yp8,zp8,ppxv,ppyv,
     1 ppzv,ppv,eexv,eeyv,eezv,eev,xrast0,yrast0')
      return
      END