This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / 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