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 C**
  28 C**  DO NOT BE CONFUSED, THIS HAS NOTHING TO DO WITH RELATIVISTIC
  29 C**  KINEMATIC
  30 C
  31       TARMAS=RMASS(14)
  32       IF (NFL .EQ. 2) TARMAS=RMASS(16)
  33       EKCOR=1.
  34       IF(EK.GT.1.) EKCOR=1./EK
  35       EK=2.*TARMAS+EK*(1.+EKCOR/ATNO2)
  36       EN=EK+ABS(AMAS)
  37       P =SQRT(ABS(EN*EN-AMAS*AMAS))
  38       S =AMAS*AMAS+TARMAS**2+2.0*TARMAS*EN
  39       RS=SQRT(S)
  40       ENP(5)=EK
  41       ENP(6)=EN
  42       ENP(7)=P
  43       ENP(8)=S
  44       ENP(9)=RS
  45 C**
  46 C**  EVAPORATION
  47 C**
  48       TKIN=EXNU(EK)
  49       ENP(5)=EK-TKIN
  50       IF(ENP(5).LT.0.0001) ENP(5)=0.0001
  51       ENP(6)=ENP(5)+ABS(AMAS)
  52       ENP(7)=ENP(6)*ENP(6)-AMASQ
  53       ENP(7)=SQRT(ABS(ENP(7)))
  54       ENP(8)=AMASQ+RMASS(14)**2+2.*RMASS(14)*ENP(6)
  55       ENP(9)=SQRT(ENP(8))
  56 C**  CHECK AVAILABLE ENERGY FOR FIRST INTERACTION
  57       IF(ENP(5).GT.CENG(3)) GO TO 9999
  58       ENP(5)=0.
  59       ENP(6)=ABS(AMAS)
  60       ENP(7)=0.
  61       ENP(8)=4.*RMASS(14)**2
  62       ENP(9)=2.*RMASS(14)
  63 C
  64  9999 CONTINUE
  65       END