GEANT321/neutron/angcdf.F

   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