c**********************************************************************
c
      subroutine richbook
c
c    open 6 RICH n-tuples (2001-2006) with layout:
c
c     IDEVT     <--  GEANT event id
c     NCER      <--  nb. of Cerenkov photons in this sector
c     XCER      <--  x in padplane (cm)
c     YCER      <--  y in padplane (cm)
c     ECER      <--  energy in eV
c     PARENTCER <--  parent track nb.
c     NHIT      <--  nb. of direct hits into pad plane
c     XHIT      <--  x of direct hit in pad plane (cm)
c     YHIT      <--  y (cm)
c     ZHIT      <--  z (cm)
c     PARTHIT   <--  particle GEANT ID of hit
c     MOMHIT    <--  momentum (MeV)
c     THETAHIT  <--  theta of incidence  (0-180 deg)
c     PHIHIT    <--  azimuthal angle (0-360 deg)
c     ELOSHIT   <--  energy loss of hit in MeV
c     TLENHIT   <--  track length of hit in gas detector
c     TRKHIT    <--  track number of direct hit
c
c
c    designed by W. Przygoda (TUM)
c    last modified on: 07/03/2001  by R. Holzmann
c**********************************************************************
      implicit none

#include "richtups.inc"
#include "geant321/gcflag.inc"

      integer icycle, err, isect

      call hbset('BSIZE',1024,err)
      if (err.ne.0) write(6,*) 'RICH ntuple: Memory allocation error!'

      do isect=1,nsect
         call hbnt(2000+isect,'RICH'//char(48+isect),' ')
         if(isect.ne.0)
     &        call hbname(2000+isect,'IEVT',idevt,'IDEVT:U')  ! event number
         call hbname(2000+isect,'Cerenkov',ncer,      ! Cerenkov photons
     &             'NCER[0,4000]:I,'//
     &             'XCER(NCER):R,YCER(NCER):R,ECER(NCER):R,'//
     &             'PARENTCER(NCER)[0,4000]:U')
         call hbname(2000+isect,'Direct',nhit,        ! direct hits into RICH
     &             'NHIT[0,400]:I,'//
     &             'XHIT(NHIT):R,YHIT(NHIT):R,ZHIT(NHIT):R,'//
     &             'PARTHIT(NHIT)[0,4000]:U,MOMHIT(NHIT):R,'//
     &             'THETAHIT(NHIT):R,PHIHIT(NHIT):R,'//
     &             'ELOSHIT(NHIT):R,TLENHIT(NHIT):R,'//
     &             'TRKHIT(NHIT)[0,4000]:U')

      enddo
      return

**********************************************************

      entry richbooklast

**********************************************************

      do isect=1,nsect
         call hrout(2000+isect,icycle,' ')
c        call hprint(2000+isect)
      enddo

      return
      end