GEANT321/gheisha/selfab.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:20:59  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.38  by  S.Giani
  11 *-- Author :
  12       SUBROUTINE SELFAB(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 #include "geant321/mxgkgh.inc"
  20 #include "geant321/s_curpar.inc"
  21 #include "geant321/s_blank.inc"
  22 C
  23       SPROB=0.
  24       EKW=PV(4,MXGKPV)-ABS(PV(5,MXGKPV))
  25       IF(EKW.LT.5.) RETURN
  26       ALEKW=LOG(EKW-4.)
  27       SPROB=0.6*ALEKW
  28       IF(SPROB.GT.1.) SPROB=1.
  29       RETURN
  30       END