***********************************************************************
*      RICHTUPS.INC                                                   *
*                                                                     *
*     last modified on 05/12/2008 by I. Koenig                        *
***********************************************************************

#include "hgeantmaxtrk.h"

      INTEGER MAXCKOV, MAXPART, MAXMIRR, NSECT

      PARAMETER (MAXCKOV=MAXCKOVRICH)
      PARAMETER (MAXPART=MAXPARTRICH)
      PARAMETER (MAXMIRR=MAXMIRRRICH)
      PARAMETER (NSECT=600)

      INTEGER NMIRCOUNT                 ! auxiliary
      INTEGER NUMPHOTALL(MAXMIRR)       ! number of reflected photons
      INTEGER LEPTRACKALL(MAXMIRR)      ! parent track number
      REAL XR(MAXMIRR)
      REAL YR(MAXMIRR)

      INTEGER NMIR                      ! number of leptons hits on mirror
      INTEGER NUMPHOT(MAXMIRR)          ! number of reflected photons
      INTEGER LEPTRACK(MAXMIRR)         ! parent track number
                                        ! for the reflected photons
      REAL XRING(MAXMIRR)               ! expected ring X-position
      REAL YRING(MAXMIRR)               ! expected ring Y-position
      REAL XLEP(MAXMIRR)
      REAL YLEP(MAXMIRR)
      REAL ZLEP(MAXMIRR)
      INTEGER LEPID(MAXMIRR)


      INTEGER NCER			! number of photon hits
      REAL XCER(MAXCKOV)		! x-coord of photon
      REAL YCER(MAXCKOV)		! y-coord of photon
      REAL ECER(MAXCKOV)		! photon momentum
      INTEGER PARENTCER(MAXCKOV)	! parent track

      INTEGER NHIT			! number of particle hits
      REAL XHIT(MAXPART)		! x-coord of particle
      REAL YHIT(MAXPART)		! y-coord of particle
      REAL ZHIT(MAXPART)		! z-coord of particle
      INTEGER PARTHIT(MAXPART)		! GEANT particle id
      REAL MOMHIT(MAXPART)		! momentum
      REAL THETAHIT(MAXPART)		! theta incident angle
      REAL PHIHIT(MAXPART)		! phi incident angle
      REAL ELOSHIT(MAXPART)		! energy lost in gas of detector
      REAL TLENHIT(MAXPART)		! track length in detector
      INTEGER TRKHIT(MAXPART)		! track number

      REAL ENLOS
      SAVE ENLOS
      REAL EVTNOW
      SAVE EVTNOW

      INTEGER MED_C4F10 	! C4F10 medium number
      INTEGER MED_C4H10 	! C4H10 medium number
      INTEGER MED_GLASS_PMT 	! Glass PMT medium number
      INTEGER MED_CAF2		! CaF2 medium number
      INTEGER MED_CH4		! methan medium number
      INTEGER MED_N2_DIS       	! richh700 (instead c4h10)
      INTEGER MED_DET		! methan (volume of photon detector)
      INTEGER MED_MIRROR	! carbon mirror medium number
      INTEGER MED_GLASSMIRROR	! glass mirror medium number
      INTEGER MED_PADS		! copper pads medium number
      INTEGER MED_SIO2A         ! SIO2A medium number
      INTEGER MED_SIO2B         ! SIO2B medium number
      INTEGER MED_MGF2A         ! MgF2A medium number
      INTEGER MED_MGF2B         ! MgF2B medium number
      INTEGER MED_N2            ! N2 medium number

      INTEGER ISRICHPMT         ! 1 if PMT version of RICH

      COMMON /RICHTUPS/ NCER, XCER, YCER, ECER, PARENTCER,
     &                  NHIT, XHIT, YHIT, ZHIT, PARTHIT, MOMHIT,
     &                  THETAHIT, PHIHIT, ELOSHIT, TLENHIT, TRKHIT,
     &	                NMIR, NUMPHOT, LEPTRACK,
     &			XRING, YRING,
     &                  XLEP, YLEP, ZLEP, LEPID
      COMMON /RICHDATA/ NMIRCOUNT,NUMPHOTALL,LEPTRACKALL,XR,YR
      COMMON /RICHMED/ MED_C4F10, MED_C4H10, MED_CAF2, MED_CH4,
     &                 MED_SIO2A, MED_SIO2B, MED_GLASS_PMT,
     &                 MED_MGF2A, MED_MGF2B, MED_N2,
     &                 MED_DET, MED_MIRROR,
     &                 MED_GLASSMIRROR, MED_PADS,
     &                 MED_N2_DIS
      COMMON /RICHFLAG/ ISRICHPMT