GEANT321/neutron/trebod.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:21:59  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 TREBOD(D,LD,KZ1,KZ2,KZ3,A1,A2,A3,Z1,Z2,Z3,
  13      +                  ATAR,Q,MT)
  14 C CZ July 30,1992 Simple aproach to get (N,PA), (N,T2A),(N,D2A)
  15 C processes. This is TWOBOD extended to a third particle
  16 C       THIS ROUTINE CALCULATES THE EXIT ENERGIES AND DIRECTIONAL
  17 C       COSINES FOR THE CHARGED PARTICLE AND RECOIL NUCLEUS FOR
  18 C       A THREE-BODY REACTION USING AN EVAPORATION SPECTRUM AND
  19 C       MOMEMTUM BALANCE.  IT ALSO SETS ALL EXIT PARAMETERS FOR
  20 C       THE COLLISION PRODUCTS AND STORES THEM IN THE RECOIL BANK.
  21 #include "geant321/minput.inc"
  22 #include "geant321/mconst.inc"
  23 #include "geant321/mnutrn.inc"
  24 #include "geant321/mrecoi.inc"
  25 #include "geant321/mapoll.inc"
  26 #include "geant321/mmass.inc"
  27 #include "geant321/mpstor.inc"
  28       DIMENSION D(*),LD(*),ER1(3)
  29       SAVE
  30 C loop over no. of emmitted particles CZ July 30,1992
  31       NPN = 1
  32       IF(MT.EQ.112) NPN = 2
  33       IF(MT.EQ.113) NPN = 3
  34       IF(MT.EQ.114) NPN = 3
  35       PRXO = 0.0
  36       PRYO = 0.0
  37       PRZO = 0.0
  38       DO 10  NP=1,NPN
  39 C       CALCULATE THE COULOMB BARRIER (CB)
  40          CALL BARIER(KZ1,KZ2,A1,A3,CB)
  41 C       CALCULATE THE CHARGED PARTICLE EXIT ENERGY (EX)
  42          CALL CEVAP(EOLD,Q,ATAR,CB,EX)
  43          E1=EX+CB
  44          ZMSS = Z2
  45          AMSS = A2
  46          KZZ = KZ2
  47          IF(NP.EQ.1) THEN
  48             ZMSS = Z1
  49             AMSS = A1
  50             KZZ = KZ1
  51          ENDIF
  52 C       ASSUME ISOTROPIC CHARGED PARTICLE EMISSION IN THE LABORATORY
  53          CALL GTISO(U1,V1,W1)
  54          PPN = SQRT(2.0*ZMSS*E1)
  55          PRXO = PRXO + U1*PPN
  56          PRYO = PRYO + V1*PPN
  57          PRZO = PRZO + W1*PPN
  58 C       CALCULATE AND SET THE CHARGED PARTICLE EXIT PARAMETERS
  59          XR=X
  60          YR=Y
  61          ZR=Z
  62          WATER=WTBC
  63          NZR=KZZ
  64          AGER=AGE
  65          NCOLR=NCOL
  66          MTNR=MT
  67          AR=AMSS
  68          ENIR=EOLD
  69          UNIR=UOLD
  70          VNIR=VOLD
  71          WNIR=WOLD
  72          ENOR=0.0
  73          UNOR=0.0
  74          VNOR=0.0
  75          WNOR=0.0
  76          WTNR=0.0
  77          QR=Q
  78          UR=U1
  79          VR=V1
  80          WR=W1
  81          ER=E1
  82 C       STORE THE CHARGED PARTICLE IN THE RECOIL BANK
  83          EP = ER
  84          UP = UR
  85          VP = VR
  86          WP = WR
  87          AMP = AR
  88          ZMP = FLOAT(NZR)
  89          AGEP = AGE
  90          MTP = MT
  91          CALL STOPAR(IDHEVY,NHEVY)
  92          A3 = A3 - A2
  93          Z3 = Z3 - Z2
  94          KZ3 = KZ3 - KZ2
  95    10 CONTINUE
  96       A3 = A3 + A2
  97       Z3 = Z3 + Z2
  98       KZ3 = KZ3 + KZ2
  99 C       CALCULATE THE TOTAL MOMENTUM BEFORE THE COLLISION
 100 C       NEUTRON MOMENTUM BEFORE COLLISION (PI) EQUALS TOTAL MOMENTUM
 101       PI=SQRT(2.0*ZN*EOLD)
 102 C       CALCULATE THE DIRECTIONAL MOMENTUM OF THE RECOIL NUCLEUS
 103       PRX=PI*UOLD - PRXO
 104       PRY=PI*VOLD - PRYO
 105       PRZ=PI*WOLD - PRZO
 106 C       CALCULATE THE TOTAL MOMENTUM OF THE RECOIL NUCLEUS
 107       PR=SQRT(PRX**2+PRY**2+PRZ**2)
 108 C       CALCULATE THE RECOIL NUCLEUS DIRECTIONAL COSINES
 109       U2=PRX/PR
 110       V2=PRY/PR
 111       W2=PRZ/PR
 112 C       CALCULATE THE RECOIL NUCLEUS EXIT ENERGY
 113       XM = A2 * 931.075E6
 114       E2 = SQRT(PR**2+XM**2) - XM
 115 C       CALCULATE AND SET THE CHARGED PARTICLE EXIT PARAMETERS
 116       XR=X
 117       YR=Y
 118       ZR=Z
 119       WATER=WTBC
 120       NZR=KZ3
 121       AGER=AGE
 122       NCOLR=NCOL
 123       MTNR=MT
 124       AR=A3
 125       ENIR=EOLD
 126       UNIR=UOLD
 127       VNIR=VOLD
 128       WNIR=WOLD
 129       ENOR=0.0
 130       UNOR=0.0
 131       VNOR=0.0
 132       WNOR=0.0
 133       WTNR=0.0
 134       QR=Q
 135       UR=U2
 136       VR=V2
 137       WR=W2
 138       ER=E2
 139 C       STORE THE RECOIL HEAVY ION IN THE RECOIL BANK
 140       EP = ER
 141       UP = UR
 142       VP = VR
 143       WP = WR
 144       AMP = AR
 145       ZMP = FLOAT(NZR)
 146       AGEP = AGE
 147       MTP = MT
 148       CALL STOPAR(IDHEVY,NHEVY)
 149       RETURN
 150       END