[u/mrichter/AliRoot.git] / GEANT321 / gheisha / selfab.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1995/10/24 10:20:59  cernlib
* Geant
*
*
#include "geant321/pilot.h"
*CMZ :  3.21/02 29/03/94  15.41.38  by  S.Giani
*-- Author :
      SUBROUTINE SELFAB(SPROB)
C
C *** SELF-ABSORBTION IN HEAVY MOLECULES ***
C *** NVE 16-MAR-1988 CERN GENEVA ***
C
C ORIGIN : H.FESEFELDT (11-OCT-1987)
C
#include "geant321/mxgkgh.inc"
#include "geant321/s_curpar.inc"
#include "geant321/s_blank.inc"
C
      SPROB=0.
      EKW=PV(4,MXGKPV)-ABS(PV(5,MXGKPV))
      IF(EKW.LT.5.) RETURN
      ALEKW=LOG(EKW-4.)
      SPROB=0.6*ALEKW
      IF(SPROB.GT.1.) SPROB=1.
      RETURN
      END
Commit	Line	Data
fe4da5cc	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