GEANT321/fluka/ankeka.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 ANKEKA.FOR
  13 *COPY ANKEKA
  14 *                                                                      *
  15 *=== ankeka ===========================================================*
  16 *                                                                      *
  17       FUNCTION ANKEKA(IT,TO,AMSS,SQAMSS)
  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     TO BE CALLED FROM THE HIGH ENERGY PRODUCTION
  30 C
  31 C     TO GIVE THE PROTON AND NEUTRON MULTIPLICITIES OF THE HIGH
  32 C     ENERGY PART OF THE INTRANUCLEAR CASCADE  (N-2P,N-2N).
  33 C     NOTE THAT IT IS SUPPOSED IN EKEKA AND RAKEKA THAT
  34 C     N-1P=9*N-2P, N-1N=9*N-2N.
  35 C
  36 C     NOTE!!!!!!! NON-STANDARD PARTICLE NUMBERING
  37 C
  38 C     INPUT VARIABLES:
  39 C        IT     = TYPE OF THE PARTICLE
  40 C                 1=INTRANUCLEAR PROTON
  41 C                 2=INTRANUCLEAR NEUTRON
  42 C        TO     = ENERGY OF THE COLLIDING PARTICLE IN GEV
  43 C        AMSS   = ATOMIC WEIGHT OF THE MEDIUM
  44 C        SQAMSS = SQRT(AMSS)
  45 C
  46 C     SEE RANFT/ROUTTI PARTICLE ACC. VOL 4 P 105
  47 C********************************************************************
  48 C
  49       DIMENSION A(2),B(2),C(2),D(2)
  50       SAVE A,B,C,D
  51 C
  52 C
  53 C     THESE ARE EFFECTIVE FLUKA81 MULTIPLICITIES
  54 C
  55       DATA A/1.D0,1.3D0/
  56       DATA B/0.06D0,0.09D0/
  57       DATA C/1.104D0,1.90D0/
  58 C     DATA D/0.14D0,0.19D0/
  59 C  This ones for continuity !!
  60       DATA D/0.245333D0,0.33333D0/
  61       IF(TO.GT.0.1D0)GO TO 1
  62       ANKEKA=SQAMSS*B(IT)*0.1D0
  63       RETURN
  64  1    CONTINUE
  65       IF(TO.GE.10.D0)GO TO 2
  66       ANKEKA=0.1D0*SQAMSS*(0.5D0+A(IT)*(1.D0+LOG10(TO))**2)*D(IT)
  67       RETURN
  68  2    CONTINUE
  69       ANKEKA=SQAMSS*0.1D0*C(IT)
  70       RETURN
  71       END