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