subroutine lumisel(elume,xlume,ylume,elumg,xlumg,ylumg)
      implicit none
#include "partap.inc"
#include "lmresu.inc"
      real elume,xlume,ylume,elumg,xlumg,ylumg
      real xe,ye
      integer fle,flg

C get energy from LMRESU table
 
      if(coutab(LMRESU) .ne. 0 ) then
        call fettab(LMRESU,id,1)
          elumg = LMRESU_enrg
          elume = LMRESU_enre
      else
          elumg = 0.0
          elume = 0.0
      endif

C call PHANTOM routines lmpose,lmposg to get hit position 

      call lmpose(xlume,xe,ylume,ye,fle)
      call lmposg(xlumg,xe,ylumg,ye,flg)

C put lumi position to -100. if it is not measured.

      if(fle.eq.0) then
         xlume = -100.
         ylume = -100.
      endif

      if(flg.eq.0) then
         xlumg = -100. 
         ylumg = -100. 
      endif

      return
      end