]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 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 |