GEANT321/neutron/nn2bod.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 NN2BOD(D,LD,KZ1,KZ2,A1,A2,Z1,Z2,ATAR,Q,MT)
  13 C       THIS ROUTINE CALCULATES THE EXIT ENERGIES AND DIRECTIONAL
  14 C       COSINES FOR THE CHARGED PARTICLE AND RECOIL NUCLEUS FOR
  15 C       A TWO-BODY REACTION USING AN EVAPORATION SPECTRUM AND
  16 C       MOMEMTUM BALANCE.  IT ALSO SETS ALL EXIT PARAMETERS FOR
  17 C       THE COLLISION PRODUCTS AND STORES THEM IN THE RECOIL BANK.
  18 C       THE TWO BODY REACTION RESULTS FROM THE BREAK-UP OF A NUCLEUS
  19 C       LEFT IN AN EXCITED STATE BY AN INELASTIC COLLISION OR A
  20 C       N,2N REACTION (I.E. MT-24).
  21 #include "geant321/minput.inc"
  22 #include "geant321/mconst.inc"
  23 #include "geant321/mrecoi.inc"
  24 #include "geant321/mapoll.inc"
  25 #include "geant321/mmass.inc"
  26 #include "geant321/mpstor.inc"
  27 #include "geant321/mnutrn.inc"
  28       DIMENSION D(*),LD(*)
  29       SAVE
  30 C       TRANSFER THE RECOILING COMPOUND NUCLEUS PARAMETERS OUT OF
  31 C       COMMON RECOIL FOR USE IN THE MOMENTUM BALANCE EQUATIONS
  32       ERCN=ER
  33       URCN=UR
  34       VRCN=VR
  35       WRCN=WR
  36       ARCN=AR
  37       NZRCN=NZR
  38       ZARCN=ARCN*9.31075E+08
  39 C       CALCULATE THE COULOMB BARRIER (CB)
  40       CALL BARIER(KZ1,KZ2,A1,A2,CB)
  41 C       CALCULATE THE CHARGED PARTICLE EXIT ENERGY (EX)
  42       CALL CEVAP1(EOLD,E,Q,ATAR,CB,EX)
  43       E1=EX+CB
  44 C       ASSUME ISOTROPIC CHARGED PARTICLE EMISSION IN THE LABORATORY
  45       CALL GTISO(U1,V1,W1)
  46 C       CALCULATE AND SET THE CHARGED PARTICLE EXIT PARAMETERS
  47       XR=X
  48       YR=Y
  49       ZR=Z
  50       WATER=WTBC
  51       NZR=KZ1
  52       AGER=AGE
  53       NCOLR=NCOL
  54       MTNR=MT
  55       AR=A1
  56       ENIR=EOLD
  57       UNIR=UOLD
  58       VNIR=VOLD
  59       WNIR=WOLD
  60       ENOR=E
  61       UNOR=U
  62       VNOR=V
  63       WNOR=W
  64       WTNR=WATE
  65       QR=Q
  66       UR=U1
  67       VR=V1
  68       WR=W1
  69       ER=E1
  70 C       STORE THE CHARGED PARTICLE IN THE RECOIL BANK
  71       EP = ER
  72       UP = UR
  73       VP = VR
  74       WP = WR
  75       AGEP = AGE
  76       MTP = MT
  77       AMP = AR
  78       ZMP = FLOAT(NZR)
  79       CALL STOPAR(IDHEVY,NHEVY)
  80 C       CALCULATE THE TOTAL MOMENTUM BEFORE THE COLLISION
  81 C       COMPOUND NUCLEUS MOMENTUM BEFORE THE COLLISION (PI) EQUALS
  82 C       THE TOTAL MOMENTUM
  83       PI=SQRT(2.0*ZARCN*ERCN)
  84 C       CALCULATE THE TOTAL MOMEMTUM OF THE EXIT CHARGED PARTICLE
  85       PO=SQRT(2.0*Z1*E1)
  86 C       CALCULATE THE DIRECTIONAL MOMENTUM OF THE RECOIL NUCLEUS
  87       PRX=PI*URCN-PO*U1
  88       PRY=PI*VRCN-PO*V1
  89       PRZ=PI*WRCN-PO*W1
  90 C       CALCULATE THE TOTAL MOMENTUM OF THE RECOIL NUCLEUS
  91       PR=SQRT(PRX**2+PRY**2+PRZ**2)
  92 C       CALCULATE THE RECOIL NUCLEUS DIRECTIONAL COSINES
  93       U2=PRX/PR
  94       V2=PRY/PR
  95       W2=PRZ/PR
  96 C       CALCULATE THE RECOIL NUCLEUS EXIT ENERGY
  97       XM = A2 * 931.075E6
  98       E2 = SQRT(PR**2 + XM**2) - XM
  99 C       CALCULATE AND SET THE CHARGED PARTICLE EXIT PARAMETERS
 100       XR=X
 101       YR=Y
 102       ZR=Z
 103       WATER=WTBC
 104       NZR=KZ2
 105       AGER=AGE
 106       NCOLR=NCOL
 107       MTNR=MT
 108       AR=A2
 109       ENIR=EOLD
 110       UNIR=UOLD
 111       VNIR=VOLD
 112       WNIR=WOLD
 113       ENOR=E
 114       UNOR=U
 115       VNOR=V
 116       WNOR=W
 117       WTNR=WATE
 118       QR=Q
 119       UR=U2
 120       VR=V2
 121       WR=W2
 122       ER=E2
 123 C       STORE THE RECOIL HEAVY ION IN THE RECOIL BANK
 124       EP = ER
 125       UP = UR
 126       VP = VR
 127       WP = WR
 128       AGEP = AGE
 129       MTP = MT
 130       AMP = AR
 131       ZMP = FLOAT(NZR)
 132       CALL STOPAR(IDHEVY,NHEVY)
 133       RETURN
 134       END