GEANT321/gheisha/coranh.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:21:01  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.39  by  S.Giani
  11 *-- Author :
  12       SUBROUTINE CORANH(NIHIL,NFL)
  13 C
  14 C *** NUCLEAR INTERACTIONS FOR HEAVY FRAGMENTS ***
  15 C *** NVE 06-MAY-1988 CERN GENEVA ***
  16 C
  17 C ORIGIN : H.FESEFELDT (09-JULY-1987)
  18 C
  19 #include "geant321/s_defcom.inc"
  20 C
  21       NIHIL=0
  22       IF(AMAS.GT.0.)   GO TO 9999
  23       IF(IPART.LT.14)  GO TO 9999
  24       IF(IPA(1).GE.14) GO TO 9999
  25       IF(IPA(2).GE.14) GO TO 9999
  26       NIHIL=1
  27       RETURN
  28 C**
  29 C**  DO NOT BE CONFUSED, THIS HAS NOTHING TO DO WITH RELATIVISTIC
  30 C**  KINEMATIC
  31 C
  32       TARMAS=RMASS(14)
  33       IF (NFL .EQ. 2) TARMAS=RMASS(16)
  34       EKCOR=1.
  35       IF(EK.GT.1.) EKCOR=1./EK
  36       EK=2.*TARMAS+EK*(1.+EKCOR/ATNO2)
  37       EN=EK+ABS(AMAS)
  38       P =SQRT(ABS(EN*EN-AMAS*AMAS))
  39       S =AMAS*AMAS+TARMAS**2+2.0*TARMAS*EN
  40       RS=SQRT(S)
  41       ENP(5)=EK
  42       ENP(6)=EN
  43       ENP(7)=P
  44       ENP(8)=S
  45       ENP(9)=RS
  46 C**
  47 C**  EVAPORATION
  48 C**
  49       TKIN=EXNU(EK)
  50       ENP(5)=EK-TKIN
  51       IF(ENP(5).LT.0.0001) ENP(5)=0.0001
  52       ENP(6)=ENP(5)+ABS(AMAS)
  53       ENP(7)=ENP(6)*ENP(6)-AMASQ
  54       ENP(7)=SQRT(ABS(ENP(7)))
  55       ENP(8)=AMASQ+RMASS(14)**2+2.*RMASS(14)*ENP(6)
  56       ENP(9)=SQRT(ENP(8))
  57 C**  CHECK AVAILABLE ENERGY FOR FIRST INTERACTION
  58       IF(ENP(5).GT.CENG(3)) GO TO 9999
  59       ENP(5)=0.
  60       ENP(6)=ABS(AMAS)
  61       ENP(7)=0.
  62       ENP(8)=4.*RMASS(14)**2
  63       ENP(9)=2.*RMASS(14)
  64 C
  65  9999 CONTINUE
  66       END