[u/mrichter/AliRoot.git] / GEANT321 / fluka / sitsao.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1995/10/24 10:20:05  cernlib
* Geant
*
*
#include "geant321/pilot.h"
*CMZ :  3.21/02 29/03/94  15.41.45  by  S.Giani
*-- Author :
*$ CREATE SITSAO.FOR
*COPY SITSAO
*
*=== sitsao ===========================================================*
*
      FUNCTION SITSAO ( E, IZ, A )
 
#include "geant321/dblprc.inc"
#include "geant321/dimpar.inc"
#include "geant321/iounit.inc"
*
*----------------------------------------------------------------------*
*                                                                      *
*     Original from A. Fasso`, CERN-TIS                                *
*                                                                      *
*     Created on 20 september 1991 by    Alfredo Ferrari & Paola Sala  *
*                                                   Infn - Milan       *
*                                                                      *
*     Last change on 20-sep-91     by    Alfredo Ferrari               *
*                                                                      *
*                                                                      *
*----------------------------------------------------------------------*
*
*----------------------------------------------------------------------*
*   Calculates proton-nucleus inelastic cross-section by the formula   *
*   of Letaw, Silberberg and Tsao, Astrophys. J. Suppl. 51, 271 (1983) *
*----------------------------------------------------------------------*
*   E = Energy in MeV, IZ = Atomic number, A = Mass number (At. weight)*
*----------------------------------------------------------------------*
*                                                formula (2)
      SITSAO = 45.D+00 * A**0.7D+00 * ( 1.D+00 + 0.016D+00 * SIN
     &       ( 5.3D+00 - 2.63D+00 * LOG(A)) )
*                                                formula (5)
      IF(E .LT. 2000.D+00) SITSAO = SITSAO * ( 1.D+00 - 0.62D+00
     &                            / EXP ( E / 200.D+00 )
     &                            * SIN ( 10.9D+00 / E**0.28D+00 ) )
      IF (IZ .EQ. 2) THEN
         SITSAO=SITSAO*0.8D+00
      ELSE IF (IZ .EQ. 4) THEN
         SITSAO=SITSAO*(1.D+00 + 0.75D+00/EXP(E/75.D+00))
      END IF
      RETURN
      END
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.1 1995/10/24 10:20:05 cernlib
	6	* Geant
	7	*
	8	*
	9	#include "geant321/pilot.h"
	10	*CMZ : 3.21/02 29/03/94 15.41.45 by S.Giani
	11	*-- Author :
	12	*$ CREATE SITSAO.FOR
	13	*COPY SITSAO
	14	*
	15	=== sitsao ===========================================================
	16	*
	17	FUNCTION SITSAO ( E, IZ, A )
	18
	19	#include "geant321/dblprc.inc"
	20	#include "geant321/dimpar.inc"
	21	#include "geant321/iounit.inc"
	22	*
	23	----------------------------------------------------------------------
	24	* *
	25	* Original from A. Fasso`, CERN-TIS *
	26	* *
	27	* Created on 20 september 1991 by Alfredo Ferrari & Paola Sala *
	28	* Infn - Milan *
	29	* *
	30	* Last change on 20-sep-91 by Alfredo Ferrari *
	31	* *
	32	* *
	33	----------------------------------------------------------------------
	34	*
	35	----------------------------------------------------------------------
	36	* Calculates proton-nucleus inelastic cross-section by the formula *
	37	* of Letaw, Silberberg and Tsao, Astrophys. J. Suppl. 51, 271 (1983) *
	38	----------------------------------------------------------------------
	39	* E = Energy in MeV, IZ = Atomic number, A = Mass number (At. weight)*
	40	----------------------------------------------------------------------
	41	* formula (2)
	42	SITSAO = 45.D+00 * A*0.7D+00 ( 1.D+00 + 0.016D+00 * SIN
	43	& ( 5.3D+00 - 2.63D+00 * LOG(A)) )
	44	* formula (5)
	45	IF(E .LT. 2000.D+00) SITSAO = SITSAO * ( 1.D+00 - 0.62D+00
	46	& / EXP ( E / 200.D+00 )
	47	& * SIN ( 10.9D+00 / E**0.28D+00 ) )
	48	IF (IZ .EQ. 2) THEN
	49	SITSAO=SITSAO*0.8D+00
	50	ELSE IF (IZ .EQ. 4) THEN
	51	SITSAO=SITSAO*(1.D+00 + 0.75D+00/EXP(E/75.D+00))
	52	END IF
	53	RETURN
	54	END