c ** fill the eep-coin ntuple - a column wise ntuple - K. Aniol Feb. 12, 2005
c
      subroutine GUNTFIL
c
      character*8 dirntup
      real eepv,snth2,csth2,maprt22,mpart12
      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
      integer ntupid,lunnum
	common/eepvert/eepv(27),ipart1,ipart2,snth2,csth2,
     1  mpart22, mpart12
      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  --------------------------------------------------------------------
c eepv(1) = xvert             vertex x in target      
c eepv(2) = yvert
c eepv(3) = zvert
c eepv(4) = eex = px  electron in sp07; sp07 is immediately behind the electron spectrometer collimator
c eepv(5) = eey = py
c eepv(6) = eez = pz
c eepv(7) = ee  = ptot     for electron in sp07
c eepv(8) = ppx = px  proton in sp08; sp08 is immediately behind the proton spectrometer collimator
c eepv(9) = ppy = py
c eepv(10) = ppz = pz
c eepv(11) = ptot=pproton          for 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 ; SP01 is a vacuum region before the target in order to obtain the x,y raster values
c eepv(27) = yrast   from SP01
c mpart22 = (mass of ipart2)**2
c mpart12 = (mass of ipart1)**2


c fill variables for ntuple
        xvert=eepv(1)
        yvert=eepv(2)
        zvert=eepv(3)
        eex=eepv(4)
        eey=eepv(5)
        eez=eepv(6)
        ee =eepv(7)
        ppx=eepv(8)
        ppy=eepv(9)
        ppz=eepv(10)
        pproton=eepv(11)
        xe7=eepv(12)
        ye7=eepv(13)
        ze7=eepv(14)
        xp8=eepv(15)
        yp8=eepv(16)
        zp8=eepv(17)
        ppxv=eepv(18)
        ppyv=eepv(19)
        ppzv=eepv(20)
        ppv=eepv(21)
        eexv=eepv(22)
        eeyv=eepv(23)
        eezv=eepv(24)
        eev=eepv(25)
	xrast0 = eepv(26)
	yrast0 = eepv(27)
        call hfnt(ntupid)               !fill ntuple
        return
        end