GEANT321/fluka/hinhev.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:19:57  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.43  by  S.Giani
  11 *-- Author :
  12 *$ CREATE HINHEV.FOR
  13 *COPY HINHEV
  14 *
  15 *=== hinhev ===========================================================*
  16 *
  17       SUBROUTINE HINHEV ( N, PLAB, ELAB, ITTA )
  18
  19 #include "geant321/dblprc.inc"
  20 #include "geant321/dimpar.inc"
  21 #include "geant321/iounit.inc"
  22 *
  23 *----------------------------------------------------------------------*
  24 *                                                                      *
  25 *     Created on 01 october 1991   by    Alfredo Ferrari & Paola Sala  *
  26 *                                                   Infn - Milan       *
  27 *                                                                      *
  28 *     Last change on 01-oct-91     by    Alfredo Ferrari               *
  29 *                                                                      *
  30 *                                                                      *
  31 *----------------------------------------------------------------------*
  32 *
  33 #include "geant321/finlsp.inc"
  34 #include "geant321/hadpar.inc"
  35 #include "geant321/part.inc"
  36 *
  37       EPROJ = ELAB
  38       PPROJ = PLAB
  39       KPROJ = N
  40       KTARG = ITTA
  41       AMTAR = AM (KTARG)
  42       UMO   = SQRT ( ( EPROJ + AMTAR - PPROJ ) * ( EPROJ + AMTAR
  43      &      + PPROJ ) )
  44       NHAD  = 0
  45       CALL HADEVV ( NHAD, KPROJ, KTARG, PPROJ, EPROJ, UMO )
  46       IR = NHAD
  47       DO 1000 I = 1, NHAD
  48          ELR (I) = HEPH (I)
  49          PLR (I) = SQRT ( ( HEPH (I) - AMH (I) ) * ( HEPH (I) + AMH (I)
  50      &           ) )
  51          CXR (I) = PXH (I) / PLR (I)
  52          CYR (I) = PYH (I) / PLR (I)
  53          CZR (I) = PZH (I) / PLR (I)
  54          ITR (I) = NREH (I)
  55  1000 CONTINUE
  56 *=== End of subroutine hinhev =========================================*
  57       RETURN
  58       END