This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / GEANT321 / neutron / nghevy.F
1 *
2 * $Id$
3 *
4 * $Log$
5 * Revision 1.1.1.1  1995/10/24 10:21:58  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 NGHEVY(D,LD,KZ,AWR,Q,MT)
13 C       THIS ROUTINE CALCULATES THE EXIT ENERGY AND DIRECTIONAL
14 C       COSINES FOR THE RECOIL NUCLEUS RESULTING FROM THE (N,G)
15 C       REACTION MT-102, AND STORES THE RECOIL NUCLEUS IN THE
16 C       HEAVY ION BANK.  THE ENERGY AND DIRECTIONAL COSINES ARE
17 C       DETERMINED BY A MOMENTUM BALANCE IN THE LABORATORY SYSTEM
18 C       WITH THE PHOTONS MOMENTUM EQUAL TO ITS ENERGY.
19 #include "geant321/minput.inc"
20 #include "geant321/mconst.inc"
21 #include "geant321/mnutrn.inc"
22 #include "geant321/mrecoi.inc"
23 #include "geant321/mapoll.inc"
24 #include "geant321/mmass.inc"
25 #include "geant321/mpstor.inc"
26 #include "geant321/mgamma.inc"
27       DIMENSION D(*),LD(*)
28       SAVE
29       AR=AWR*AN+AN
30 C       CALCULATE THE TOTAL MOMENTUM BEFORE THE COLLISION
31 C       NEUTRON MOMENTUM BEFORE COLLISION (PI) EQUALS TOTAL MOMENTUM
32       PI=SQRT(2.0*ZN*EOLD)
33 C       CALCULATE THE TOTAL MOMEMTUM OF THE EXIT PHOTON
34       PO=EG*1.00E+06
35 C       CALCULATE THE DIRECTIONAL MOMENTUM OF THE RECOIL NUCLEUS
36       PRX=PI*UOLD-PO*UG
37       PRY=PI*VOLD-PO*VG
38       PRZ=PI*WOLD-PO*WG
39 C       CALCULATE THE TOTAL MOMENTUM OF THE RECOIL NUCLEUS
40       PR=SQRT(PRX**2+PRY**2+PRZ**2)
41 C       CALCULATE THE RECOIL NUCLEUS DIRECTIONAL COSINES
42       UR=PRX/PR
43       VR=PRY/PR
44       WR=PRZ/PR
45 C       CALCULATE THE RECOIL NUCLEUS EXIT ENERGY
46       ER=PR**2/(2*AR*9.31075E+08)
47 C       CALCULATE AND SET THE CHARGED PARTICLE EXIT PARAMETERS
48       XR=X
49       YR=Y
50       ZR=Z
51       WATER=WTBC
52       NZR=KZ
53       AGER=AGE
54       NCOLR=NCOL
55       MTNR=MT
56       ENIR=EOLD
57       UNIR=UOLD
58       VNIR=VOLD
59       WNIR=WOLD
60       ENOR=0.0
61       UNOR=0.0
62       VNOR=0.0
63       WNOR=0.0
64       WTNR=0.0
65       QR=Q
66 C       STORE THE RECOIL HEAVY ION IN THE RECOIL BANK
67       EP = ER
68       UP = UR
69       VP = VR
70       WP = WR
71       AMP = AR
72       ZMP = FLOAT(NZR)
73       AGEP = AGE
74       MTP = MT
75       CALL STOPAR(IDHEVY,NHEVY)
76       RETURN
77       END