Bugfix in AliPoints2Memory
[u/mrichter/AliRoot.git] / GEANT321 / ghrout / highab.F
CommitLineData
fe4da5cc 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)
13C
14C *** SELF-ABSORBTION IN HEAVY MOLECULES ***
15C *** NVE 16-MAR-1988 CERN GENEVA ***
16C
17C ORIGIN : H.FESEFELDT (11-OCT-1987)
18C
19C
20#include "geant321/mxgkgh.inc"
21#include "geant321/s_curpar.inc"
22#include "geant321/s_blank.inc"
23C
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