GEANT321/neutron/barier.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 BARIER(KZ1,KZ2,A1,A2,CB)
  13 C       THIS ROUTINE CALCULATES THE COULOMB BARRIER FOR A
  14 C       COLLISION INVOLVING CHARGED PARTICLE EMISSION
  15       IFLG=0
  16 C       CALCULATE THE RADIUS OF THE NUCLEUS AND CHARGED PARTICLE
  17       A=A1
  18    10 IF(A.LT.5.5)R=1.20E-13
  19       IF((A.GE.5.5).AND.(A.LT.6.5))R=2.02E-13
  20       IF((A.GE.6.5).AND.(A.LT.7.5))R=2.43E-13
  21       IF((A.GE.7.5).AND.(A.LT.8.5))R=2.84E-13
  22       IF((A.GE.8.5).AND.(A.LT.9.5))R=3.25E-13
  23       IF(A.GE.9.5)R=(A**(1.0/3.0))*1.70E-13
  24       IF(IFLG.EQ.0)R1=R
  25       IF(IFLG.EQ.1)GO TO 20
  26       IFLG=1
  27       A=A2
  28       GO TO 10
  29    20 R2=R
  30 C       CALCULATE THE COULOMB BARRIER (UNITS=MEV)
  31 C       THE FACTOR 0.75 IS ARBITRARYLY SET TO ACCOUNT FOR CHARGED
  32 C       PARTICLE EMISSION BELOW THE COULOMB BARRIER
  33       CB=((KZ1*KZ2*1.44E-13)/(R1+R2))*0.75
  34       IF(CB.LT.0.0) CB = 0.0
  35       RETURN
  36       END