[u/mrichter/AliRoot.git] / GEANT321 / neutron / angcdf.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1995/10/24 10:21:54  cernlib
* Geant
*
*
#include "geant321/pilot.h"
*CMZ :  3.21/02 29/03/94  15.41.48  by  S.Giani
*-- Author :
      SUBROUTINE ANGCDF(D,LD,LZ)
C       THIS ROUTINE READS THE INPUT ANGULAR DISTRIBUTION FILES
C       AND CONVERTS THEM TO A NORMALIZED CDF
      DIMENSION D(*),LD(*)
      IPP=1
      NR=LD(IPP)
      NE=LD(IPP+1)
      NR2=2*NR
      II=2+NR2
   10 CONTINUE
      E=D(II+1)
      NP=LD(II+2)
      A1=-1.0
      PL=D(II+4)
      D(II+4)=0.0
      PROB=0.0
      DO 20 I=2,NP
         N=II+2*I+2
         A2=D(N-1)
         PH=D(N)
         PROB=PROB+(PH+PL)*(A2-A1)/2.0
         PL=PH
         D(N)=PROB
         A1=A2
   20 CONTINUE
      DO 30 I=1,NP
         N=II+2*I+2
         D(N)=D(N)/PROB
   30 CONTINUE
      II=II+2*NP+2
      IF(II.GE.LZ)GO TO 40
      GO TO 10
   40 RETURN
      END
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.1 1995/10/24 10:21:54 cernlib
	6	* Geant
	7	*
	8	*
	9	#include "geant321/pilot.h"
	10	*CMZ : 3.21/02 29/03/94 15.41.48 by S.Giani
	11	*-- Author :
	12	SUBROUTINE ANGCDF(D,LD,LZ)
	13	C THIS ROUTINE READS THE INPUT ANGULAR DISTRIBUTION FILES
	14	C AND CONVERTS THEM TO A NORMALIZED CDF
	15	DIMENSION D(),LD()
	16	IPP=1
	17	NR=LD(IPP)
	18	NE=LD(IPP+1)
	19	NR2=2*NR
	20	II=2+NR2
	21	10 CONTINUE
	22	E=D(II+1)
	23	NP=LD(II+2)
	24	A1=-1.0
	25	PL=D(II+4)
	26	D(II+4)=0.0
	27	PROB=0.0
	28	DO 20 I=2,NP
	29	N=II+2*I+2
	30	A2=D(N-1)
	31	PH=D(N)
	32	PROB=PROB+(PH+PL)*(A2-A1)/2.0
	33	PL=PH
	34	D(N)=PROB
	35	A1=A2
	36	20 CONTINUE
	37	DO 30 I=1,NP
	38	N=II+2*I+2
	39	D(N)=D(N)/PROB
	40	30 CONTINUE
	41	II=II+2*NP+2
	42	IF(II.GE.LZ)GO TO 40
	43	GO TO 10
	44	40 RETURN
	45	END