5 * Revision 1.1.1.1 1995/10/24 10:21:01 cernlib
9 #include "geant321/pilot.h"
10 *CMZ : 3.21/02 29/03/94 15.41.39 by S.Giani
12 SUBROUTINE CORANH(NIHIL,NFL)
14 C *** NUCLEAR INTERACTIONS FOR HEAVY FRAGMENTS ***
15 C *** NVE 06-MAY-1988 CERN GENEVA ***
17 C ORIGIN : H.FESEFELDT (09-JULY-1987)
19 #include "geant321/s_defcom.inc"
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
28 C** DO NOT BE CONFUSED, THIS HAS NOTHING TO DO WITH RELATIVISTIC
32 IF (NFL .EQ. 2) TARMAS=RMASS(16)
34 IF(EK.GT.1.) EKCOR=1./EK
35 EK=2.*TARMAS+EK*(1.+EKCOR/ATNO2)
37 P =SQRT(ABS(EN*EN-AMAS*AMAS))
38 S =AMAS*AMAS+TARMAS**2+2.0*TARMAS*EN
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)
56 C** CHECK AVAILABLE ENERGY FOR FIRST INTERACTION
57 IF(ENP(5).GT.CENG(3)) GO TO 9999
61 ENP(8)=4.*RMASS(14)**2