[u/mrichter/AliRoot.git] / GEANT321 / gphys / gfrdt.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1995/10/24 10:21:25  cernlib
* Geant
*
*
#include "geant321/pilot.h"
*CMZ :  3.21/02 29/03/94  15.41.21  by  S.Giani
*-- Author :
      SUBROUTINE GFRDT(ZZ,ESHL,NSHL,NWORD,NRAD,PRB,ENER)
C.
C.    ******************************************************************
C.    *                                                                *
C     *       Fetch Radiative Decay probability Tables and the modes   *
C.    *       the radiative shells decay                               *
C.    *                                                                *
C.    *    ==>Called by : GFSHDC                                       *
C.    *       Author    J. Chwastowski                                 *
C.    *                                                                *
C.    ******************************************************************
C.
#include "geant321/gccuts.inc"
#include "geant321/gcpmxz.inc"
#include "geant321/gcphrd.inc"
      REAL ONEEV
      PARAMETER (ONEEV = 1.E-9)
      DIMENSION ESHL(*),PRB(*),ENER(*),NRAD(*)
      IZ = ZZ
      NWORD = 0
      DO 30 J = 1,NSHL
         NRAD(J) = 0
         ESHELL = ESHL(J)
         IF(ESHELL.GT.0.0) THEN
            IF(J.EQ.1) THEN
               NR = KSHLS
               NS = 0
            ELSEIF(J.EQ.2) THEN
               NR = L1SHLS
               NS = KSHLS
            ELSEIF(J.EQ.3) THEN
               NR = L2SHLS
               NS = KSHLS+L1SHLS
            ELSEIF(J.EQ.4) THEN
               NR = L3SHLS
               NS = KSHLS+L1SHLS+L2SHLS
            ENDIF
            NWS = NWORD+1
            QS = 0.0
            DO 10 I = 1,NR
               QS = QS+GPHRAT(NS+I,IZ)
               IF(GPHRAT(NS+I,IZ).GT.0.0) THEN
                  NWORD = NWORD+1
                  NRAD(J) = NRAD(J)+1
                  PRB(NWORD) = QS
                  ENER(NWORD) = ESHELL-ESHL(ISHLTR(NS+I))
                  ENER(NWORD) = ENER(NWORD)*ONEEV
                  IF(ENER(NWORD).LE.CUTGAM) ENER(NWORD) = -ENER(NWORD)
               ENDIF
   10       CONTINUE
            IF(QS.GT.0.0) THEN
               QS1 = 1./QS
               DO 20 I = NWS,NWORD
                  PRB(I) = PRB(I)*QS1
   20          CONTINUE
            ENDIF
         ENDIF
   30 CONTINUE
      END
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.1 1995/10/24 10:21:25 cernlib
	6	* Geant
	7	*
	8	*
	9	#include "geant321/pilot.h"
	10	*CMZ : 3.21/02 29/03/94 15.41.21 by S.Giani
	11	*-- Author :
	12	SUBROUTINE GFRDT(ZZ,ESHL,NSHL,NWORD,NRAD,PRB,ENER)
	13	C.
	14	C. ******************************************************************
	15	C. * *
	16	C * Fetch Radiative Decay probability Tables and the modes *
	17	C. * the radiative shells decay *
	18	C. * *
	19	C. * ==>Called by : GFSHDC *
	20	C. * Author J. Chwastowski *
	21	C. * *
	22	C. ******************************************************************
	23	C.
	24	#include "geant321/gccuts.inc"
	25	#include "geant321/gcpmxz.inc"
	26	#include "geant321/gcphrd.inc"
	27	REAL ONEEV
	28	PARAMETER (ONEEV = 1.E-9)
	29	DIMENSION ESHL(),PRB(),ENER(),NRAD()
	30	IZ = ZZ
	31	NWORD = 0
	32	DO 30 J = 1,NSHL
	33	NRAD(J) = 0
	34	ESHELL = ESHL(J)
	35	IF(ESHELL.GT.0.0) THEN
	36	IF(J.EQ.1) THEN
	37	NR = KSHLS
	38	NS = 0
	39	ELSEIF(J.EQ.2) THEN
	40	NR = L1SHLS
	41	NS = KSHLS
	42	ELSEIF(J.EQ.3) THEN
	43	NR = L2SHLS
	44	NS = KSHLS+L1SHLS
	45	ELSEIF(J.EQ.4) THEN
	46	NR = L3SHLS
	47	NS = KSHLS+L1SHLS+L2SHLS
	48	ENDIF
	49	NWS = NWORD+1
	50	QS = 0.0
	51	DO 10 I = 1,NR
	52	QS = QS+GPHRAT(NS+I,IZ)
	53	IF(GPHRAT(NS+I,IZ).GT.0.0) THEN
	54	NWORD = NWORD+1
	55	NRAD(J) = NRAD(J)+1
	56	PRB(NWORD) = QS
	57	ENER(NWORD) = ESHELL-ESHL(ISHLTR(NS+I))
	58	ENER(NWORD) = ENER(NWORD)*ONEEV
	59	IF(ENER(NWORD).LE.CUTGAM) ENER(NWORD) = -ENER(NWORD)
	60	ENDIF
	61	10 CONTINUE
	62	IF(QS.GT.0.0) THEN
	63	QS1 = 1./QS
	64	DO 20 I = NWS,NWORD
65	PRB(I) = PRB(I)*QS1
66	20 CONTINUE
67	ENDIF
68	ENDIF
69	30 CONTINUE
70	END