GEANT321/gphys/gmoli.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 *CMZ :  3.21/02 29/03/94  15.41.22  by  S.Giani
  11 *-- Author :
  12       SUBROUTINE GMOLI (AC,ZC,WMAT,NLM,DENS,OMC,CHC)
  13 C.
  14 C.    ******************************************************************
  15 C.    *                                                                *
  16 C.    *     Initializes material constants OMC,CHC for Moliere         *
  17 C.    *     scattering.                                                *
  18 C.    *                                                                *
  19 C.    * WMAT=proportion by weight of NLM  individual constituents of   *
  20 C.    *      charge ZC and atomic weight AC for the mixture/compound   *
  21 C.    * DENS=density in gm/cm**3                                       *
  22 C.    * Warning : the beta and charge dependance of ZX is omitted      *
  23 C.    *                                                                *
  24 C.    *    ==>Called by : GPROBI                                       *
  25 C.    *         Author M.S. Dixit NRCC Ottawa    *********             *
  26 C.    *                                                                *
  27 C.    ******************************************************************
  28 C.
  29 C     ZETA=factor for scattering by atomic electrons
  30 C
  31       DIMENSION WMAT(*),AC(*),ZC(*)
  32       PARAMETER (ZETA=1)
  33 *
  34 *     ------------------------------------------------------------------
  35 *
  36       ZS=0.
  37       ZE=0.
  38       ZX=0.
  39       DO 15 I=1,NLM
  40          TMP=(WMAT(I)/AC(I))*ZC(I)*(ZC(I)+ZETA)
  41          ZS=ZS+TMP
  42          ZE=ZE-0.66666*TMP*LOG(ZC(I))
  43          ZX=ZX+TMP*LOG(1.+3.34*(ZC(I)/137.036)**2)
  44 15    CONTINUE
  45       CHC=0.39612E-3*SQRT(DENS*ZS)
  46       OMC=6702.33*ZS*DENS*EXP((ZE-ZX)/ZS)
  47 *
  48       END