GEANT321/ggeom/gfstat.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:20:49  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.28  by  S.Giani
  11 *-- Author :
  12       SUBROUTINE GFSTAT(ICASE)
  13 C.
  14 C.    ******************************************************************
  15 C.    *                                                                *
  16 C.    *       FILL banks for volume statistics                         *
  17 C.    *                                                                *
  18 C.    *    ==>Called by : GINVOL, GMEDIA, GTNEXT, GTMEDI, GTRACK       *
  19 C.    *                                                                *
  20 C.    *         Author  R.Brun  *********                              *
  21 C.    *                                                                *
  22 C.    ******************************************************************
  23 C.
  24 #include "geant321/gcbank.inc"
  25 #include "geant321/gcvolu.inc"
  26 C.
  27 C.    ------------------------------------------------------------------
  28 C.
  29       IF(JGSTAT.LE.0)GO TO 999
  30       IF(NLEVEL.LE.0)GO TO 999
  31 C
  32       IF(ICASE.LE.1) THEN
  33 C
  34 C                 GINVOL   ICASE=0,1
  35 C
  36          IVO=LVOLUM(NLEVEL)
  37          JG=LQ(JGSTAT-IVO)
  38          IF(JG.GT.0) THEN
  39             IQ(JG+1+ICASE)=IQ(JG+1+ICASE)+1
  40          ENDIF
  41       ELSEIF(ICASE.LE.10) THEN
  42 C
  43 C         GMEDIA  (ICASE=2), GTNEXT (ICASE=3) or GTMEDI (ICASE=4)
  44 C
  45          IC = ICASE*2-1
  46          DO 10 ILEV=1,NLEVEL
  47             IVO=LVOLUM(ILEV)
  48             JG=LQ(JGSTAT-IVO)
  49             IF(JG.GT.0) THEN
  50                IQ(JG+IC)=IQ(JG+IC)+1
  51                IF(ILEV.EQ.NLEVEL) THEN
  52                   IQ(JG+IC+1)=IQ(JG+IC+1)+1
  53                ENDIF
  54             ENDIF
  55    10    CONTINUE
  56 C
  57 C               ICASE=10+ITRTYP  (steps)
  58 C
  59       ELSE
  60          IVO=LVOLUM(NLEVEL)
  61          JG=LQ(JGSTAT-IVO)
  62          IF(JG.GT.0) THEN
  63             IQ(JG+9)=IQ(JG+9)+1
  64          ENDIF
  65       ENDIF
  66 C
  67   999 END