program tracker_acceptance
  use precision_def
  implicit none
  type tracker_voxel
     real(rp) pVtx
     integer nGen, nVtx
  end type tracker_voxel
  type (tracker_voxel)  tvox(0:180,0:18,0:18)
  integer ix,iy,iz

  call get_tracker_acceptance(tvox)

!  do iz=0,180
!     do ix=0,18
!        do iy=0,18
!           write(11, '(3es12.4,2i12,es12.4)') iz*10+1700., ix*5-42.5,iy*5-42.5, tvox(iz,ix,iy)%nGen,tvox(iz,ix,iy)%nVtx,tvox(iz,ix,iy)%pVtx
!        end do
!     end do
!  end do

  ! iz = (z-1700)/10
  ! iy = (y+42.5)/5
  ! ixx = (x+42.5)/5

 call get_data_from_all_EndOfTracking_phase_space(tvox)


end program tracker_acceptance