GEANT321/gphys/gmoli1.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:21:27  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 #if defined(CERNLIB_VER314)
  11 *CMZ :  3.21/02 29/03/94  15.41.22  by  S.Giani
  12 *-- Author :
  13       SUBROUTINE GMOLI1(AC,ZC,DENS,BETA2,OMC)
  14 C.
  15 C.    ******************************************************************
  16 C.    *                                                                *
  17 C.    *     Initializes material constants OMC,CHC for Moliere         *
  18 C.    *     scattering (case of a simple material)                     *
  19 C.    *                                                                *
  20 C.    *    ==>Called by : GMULOF                                       *
  21 C.    *         Author M.S. Dixit NRCC Ottawa    *********             *
  22 C.    *                                                                *
  23 C.    ******************************************************************
  24 C.
  25 C     ZETA=factor for scattering by atomic electrons
  26 C
  27       DATA ZETA / 1. /
  28 *
  29 *     ------------------------------------------------------------------
  30 *
  31       ZS=(1./AC)*ZC*(ZC+ZETA)
  32       ZE=-0.66666*ZS*LOG(ZC)
  33       ZX=ZS*LOG(1.+3.34*(ZC/(BETA2*137.036))**2)
  34       OMC=6702.33*ZS*DENS*EXP((ZE-ZX)/ZS)
  35 *
  36       END
  37 #endif