GEANT321/giface/gmunui.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:21:15  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.38  by  S.Giani
  11 *-- Author :
  12       SUBROUTINE GMUNUI
  13 C
  14 C *** CALCULATION OF MUON-NUCLEUS CROSS SECTIONS IN CURRENT MATERIAL ***
  15 C *** NVE 16-MAR-1988 CERN GENEVA ***
  16 C
  17 C CALLED BY : GPHYSI
  18 C ORIGIN : F.CARMINATI, H.FESEFELDT
  19 C
  20 #include "geant321/gcbank.inc"
  21 #include "geant321/gcjloc.inc"
  22 #include "geant321/gcmulo.inc"
  23 #include "geant321/gcmate.inc"
  24 #include "geant321/gconsp.inc"
  25 #include "geant321/gctrak.inc"
  26 C
  27 C
  28 C --- GET MUON MASS FROM GEANT PARTICLE DATA BANK ---
  29       LNVE=LQ(JPART-5)
  30       HMAMU=Q(LNVE+7)
  31 C
  32 C --- Infinite mean free patch for EKIN < 1 GeV (HJT/NVE 27-sep-1990) ---
  33       SIG=0.
  34       IF(JMUNU.EQ.0)GO TO 9999
  35       IF (ELOW(IEKBIN) .LT. 1.) GO TO 90
  36 C
  37       IF(Z.LT.1.)GO TO 90
  38       EN     = ELOW(IEKBIN)+HMAMU
  39       SIG    = 0.0003
  40       IF(EN.GT.30.) SIG = SIG*(EN/30.)**0.25
  41       SIG    = SIG*DENS*AVO*1E-3
  42 C
  43   90  IF(SIG.GT.0.)THEN
  44          Q(JMUNU+IEKBIN)=1./SIG
  45       ELSE
  46          Q(JMUNU+IEKBIN)=BIG
  47       ENDIF
  48 C
  49  9999 CONTINUE
  50       END