This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / GEANT321 / ghrout / highab.F
1 *
2 * $Id$
3 *
4 * $Log$
5 * Revision 1.1.1.1  1995/10/24 10:21:13  cernlib
6 * Geant
7 *
8 *
9 #include "geant321/pilot.h"
10 *CMZ :  3.21/02 29/03/94  15.41.40  by  S.Giani
11 *-- Author :
12       SUBROUTINE HIGHAB(SPROB)
13 C
14 C *** SELF-ABSORBTION IN HEAVY MOLECULES ***
15 C *** NVE 16-MAR-1988 CERN GENEVA ***
16 C
17 C ORIGIN : H.FESEFELDT (11-OCT-1987)
18 C
19 C
20 #include "geant321/mxgkgh.inc"
21 #include "geant321/s_curpar.inc"
22 #include "geant321/s_blank.inc"
23 C
24       SPROB=0.
25       EKW=PV(4,MXGKPV)-ABS(PV(5,MXGKPV))
26       IF(EKW.LT.5.) RETURN
27       ALEKW=LOG(EKW-4.)
28       BLEKW=0.25-0.02*LOG(ATNO2)
29       SPROB=BLEKW*ALEKW
30       IF(SPROB.GT.1.) SPROB=1.
31       RETURN
32       END