GEANT321/fluka/ainel.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:20:03  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.44  by  S.Giani
  11 *-- Author :
  12 *$ CREATE AINEL.FOR
  13 *COPY AINEL
  14 *                                                                      *
  15 *=== ainel ============================================================*
  16 *                                                                      *
  17       FUNCTION AINEL(IJ,ILO,E,A,SQA)
  18
  19 #include "geant321/dblprc.inc"
  20 #include "geant321/dimpar.inc"
  21 #include "geant321/iounit.inc"
  22 C********************************************************************
  23 C     VERSION BY                     J. RANFT
  24 C                                    LEIPZIG
  25 C     LAST CHANGE 16. JULY 81 BY     PERTTI AARNIO
  26 C                                    HELSINKI UNIVERSITY OF
  27 C                                    TECHNOLOGY, FINLAND
  28 C
  29 C
  30 C     FUNCTION OF FLUKA TO GIVE INELASTICITIES FOR PRODUCTION OF
  31 C     PARTICLES ILO BY PARTICLE IJ
  32 C
  33 C     NOTE!!!!!!! NON-STANDARD PARTICLE NUMBERING
  34 C
  35 C     INPUT VARIABLES:
  36 C        IJ    = TYPE OF THE PRIMARY
  37 C                1=PROTON
  38 C                2=NEUTRON
  39 C                3=CHARGED PION
  40 C
  41 C        ILO   = TYPE OF THE SECONDARY
  42 C                1=PROTON
  43 C                2=NEUTRON
  44 C                3=CHARGED PION
  45 C                4=NEUTRAL PION
  46 C                5=NUCLEAR EXCITATION EXCLUDING LOW ENERGY
  47 C                  SECONDARIES FROM INTRANUCLEAR CASCADE
  48 C                6=INTRANUCLEAR PROTON
  49 C                7=INTRANUCLEAR NEUTRON
  50 C
  51 C        E     = KINETIC ENERGY OF THE PRIMARY IN GEV
  52 C        A     = ATOMIC WEIGHT OF THE MEDIUM
  53 C        SQA   = SQRT(A)
  54 C
  55 C     OTHER VARIABLES:
  56 C        AK    = LOOK-UP TABLE FOR INELASTCITIES
  57 C
  58 C********************************************************************
  59 C
  60       DIMENSION AK(4,3)
  61       SAVE AK
  62       DATA AK/
  63      10.35D0,0.3D0,0.25D0,0.1D0,
  64      10.3D0,0.35D0,0.25D0,0.1D0,
  65      10.15D0,0.15D0,0.5D0,0.2D0/
  66 C
  67 C
  68       GO TO (1234,1234,1234,1234,5,6,7),ILO
  69 C
  70 C
  71  1234 CONTINUE
  72       IF(E.LT.0.125D0) GO TO 10
  73       AINEL=AK(ILO,IJ)*(1.D0-EEXI(IJ,E,A)/E)
  74       RETURN
  75 C
  76 C
  77  5    CONTINUE
  78       IF(E.LT.0.125D0) GO TO 11
  79       ANEL=EEXI(IJ,E,A)/E
  80       AKEK=EKEKA(2,E,A,SQA)/E+EKEKA(3,E,A,SQA)/E
  81       AINEL=ANEL-AKEK
  82       IF (AKEK.GE.ANEL) AINEL=EKEKA(1,E,A,SQA)/E
  83       RETURN
  84 C
  85 C
  86     6 CONTINUE
  87       IF(E.LT.0.125D0) GO TO 10
  88       AINEL=EKEKA(2,E,A,SQA)/E
  89       RETURN
  90 C
  91 C
  92     7 CONTINUE
  93       IF(E.LT.0.125D0) GO TO 10
  94       AINEL=EKEKA(3,E,A,SQA)/E
  95       RETURN
  96 C
  97 C
  98    10 CONTINUE
  99       AINEL=0.D0
 100       RETURN
 101 C
 102 C
 103  11   CONTINUE
 104       AINEL=1.D0
 105       RETURN
 106       END