GEANT321/fluka/akeka.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 AKEKA.FOR
  13 *COPY AKEKA
  14 *                                                                      *
  15 *=== akeka ============================================================*
  16 *                                                                      *
  17       FUNCTION AKEKA(IT,TO,AMSS)
  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 05. DECEMBER 83    PERTTI AARNIO
  26 C                                    HELSINKI UNIVERSITY OF
  27 C                                    TECHNOLOGY, FINLAND
  28 C
  29 C
  30 C     TO BE CALLED FROM HIGH ENERGY PRODUCTION
  31 C
  32 C
  33 C     SLOPE OF THE ENERGY DISTRIBUTION OF THE HIGHER ENERGY PARTICLES
  34 C     OF THE INTRANUCLEAR CASCADE PROTONS AND NEUTRONS (ALFA-2P
  35 C     AND ALFA-2N). NOTE THAT IT IS ASSUMED IN EKEKA AND RAKEKA
  36 C     THAT ALFA-2N=6*ALFA-1N AND ALFA-2P=6*ALFA-1P.
  37 C
  38 C     NOTE!!!!!!! NON-STANDARD PARTICLE NUMBERING
  39 C
  40 C     INPUT VARIABLES:
  41 C        IT   = TYPE OF THE SECONDARY; 1=PROTON, 2=NEUTRON
  42 C        TO   = KINETIC ENERGY OF THE PRIMARY HADRON IN GEV
  43 C        AMSS = ATOMIC WEIGHT OF THE NUCLEUS
  44 C
  45 C     SEE RANFT/ROUTTI PARTICLE ACC. VOL 4 P 105
  46 C********************************************************************
  47 C
  48       DIMENSION A(2),B(2)
  49       SAVE A,B
  50       DATA  A/0.11D0,0.1D0/
  51       DATA B/0.21D0,0.20D0/
  52       IF (TO .LT. 10.D0) THEN
  53          AKEKA=(1.D0-0.001D0*AMSS)*(A(IT)+0.01D0*TO)
  54       ELSE
  55          AKEKA=B(IT)*(1.D0-0.001D0*AMSS)
  56       END IF
  57       RETURN
  58       END