This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / GEANT321 / gphys / gmolio.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 GMOLIO(AC,ZC,WMAT,NLM,DENS,BETA2,CHARG2,OMC)
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.    *                                                                *
23 C.    *    ==>Called by : GMULOF                                       *
24 C.    *         Author M.S. Dixit NRCC Ottawa    *********             *
25 C.    *                                                                *
26 C.    ******************************************************************
27 C.
28 C     ZETA=factor for scattering by atomic electrons
29 C
30       DIMENSION WMAT(*),AC(*),ZC(*)
31       PARAMETER (ZETA=1.)
32 *
33 * *** This is equal to 3.34/137.036**2
34 *
35       PARAMETER (CONST1 = 1.778595E-4)
36 *
37 *     ------------------------------------------------------------------
38 *
39       ZS     = 0.
40       ZE     = 0.
41       ZX     = 0.
42       CONST  = CONST1*CHARG2/BETA2
43       DO 10 I=1,NLM
44          TMP   = (WMAT(I)/AC(I))*ZC(I)*(ZC(I)+ZETA)
45          ZS    = ZS+TMP
46          ZE    = ZE-0.66666*TMP*LOG(ZC(I))
47          ZX    = ZX+TMP*LOG(1.+CONST*ZC(I)**2)
48   10  CONTINUE
49       OMC    = 6702.33*ZS*DENS*EXP((ZE-ZX)/ZS) * CHARG2
50 *
51       END
52