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

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1995/10/24 10:20:04  cernlib
* Geant
*
*
#include "geant321/pilot.h"
*CMZ :  3.21/02 29/03/94  15.41.44  by  S.Giani
*-- Author :
*$ CREATE NIZLNW.FOR
*COPY NIZLNW
*                                                                      *
*=== nizlnw ===========================================================*
*                                                                      *
      SUBROUTINE NIZLNW ( IT, ZZZ, AAA, EKE, PO, SI, ZL )
 
#include "geant321/dblprc.inc"
#include "geant321/dimpar.inc"
#include "geant321/iounit.inc"
C
#include "geant321/paprop.inc"
C
      PARAMETER ( AVOGMB = 1.0D+27 / AVOGAD )
C
      DIMENSION ITT (39)
      SAVE ITT
      DATA ITT/   1,   7,   0,   0,   0,   0,   0,   2,   8,   0,   0,
*              Klong pi+  pi-   K+   K- Lam  Alam Kshrt Sig- Sig+ Sig0
     &           10,   3,   4,   6,   5,   2,   8,  10,   2,   1,   2,
*               pi0   K0  AK0  pi0  res. res. res. res. Asi- Asi0 Asi+
     &            3,  10,   9,   3,   0,   0,   0,   0,   8,   8,   7,
*                X0  Ax0   X-  AX-  Om- Aom+
     &            2,   8,   2,   8,   2,   8 /
*
      SI = AZRZRZ
      ZL = AINFNT
C
C     CALCULATE THE NEW PARTICLE NUMBER IIT:   1=P,2=N,3=PI+,4=PI-,
C     5=K-,6=K+,7=P BAR,8=N BAR,9=K ZERO BAR,10=K ZERO
C
      IIT = ITT (IT)
      IF ( IIT .LE. 0 ) RETURN
*  Check the kinetic energy: no interaction below 50 MeV at present
      IF ( IIT .LE. 2 .AND. EKE .LT. 0.400D+00 .AND. AAA .GT. 1.5D+00 )
     &   THEN
         IF ( IIT .EQ. 1 ) THEN
            SI = XSEPRO ( EKE, ZZZ, AAA )
         ELSE
            SI = XSENEU ( EKE, ZZZ, AAA )
         END IF
         IF ( EKE .GT. 0.200D+00 ) THEN
            CALL NIZL ( IT, AAA, EKE, PO, SI2, ZL2 )
            WEIGH1 = ( 0.400D+00 - EKE ) / 0.200D+00
            SI = WEIGH1 * SI + ( 1.D+00 - WEIGH1 ) * SI2
         END IF
         IF ( SI .GT. 0.D+00 ) THEN
            ZL = AVOGMB * AAA / SI
         ELSE
            SI = 0.D+00
            ZL = AINFNT
         END IF
      ELSE
         CALL NIZL ( IT, AAA, EKE, PO, SI, ZL )
      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:04 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 NIZLNW.FOR
	13	*COPY NIZLNW
	14	* *
	15	=== nizlnw ===========================================================
	16	* *
	17	SUBROUTINE NIZLNW ( IT, ZZZ, AAA, EKE, PO, SI, ZL )
	18
	19	#include "geant321/dblprc.inc"
	20	#include "geant321/dimpar.inc"
	21	#include "geant321/iounit.inc"
	22	C
	23	#include "geant321/paprop.inc"
	24	C
	25	PARAMETER ( AVOGMB = 1.0D+27 / AVOGAD )
	26	C
	27	DIMENSION ITT (39)
	28	SAVE ITT
	29	DATA ITT/ 1, 7, 0, 0, 0, 0, 0, 2, 8, 0, 0,
	30	* Klong pi+ pi- K+ K- Lam Alam Kshrt Sig- Sig+ Sig0
	31	& 10, 3, 4, 6, 5, 2, 8, 10, 2, 1, 2,
	32	* pi0 K0 AK0 pi0 res. res. res. res. Asi- Asi0 Asi+
	33	& 3, 10, 9, 3, 0, 0, 0, 0, 8, 8, 7,
	34	* X0 Ax0 X- AX- Om- Aom+
	35	& 2, 8, 2, 8, 2, 8 /
	36	*
	37	SI = AZRZRZ
	38	ZL = AINFNT
	39	C
	40	C CALCULATE THE NEW PARTICLE NUMBER IIT: 1=P,2=N,3=PI+,4=PI-,
	41	C 5=K-,6=K+,7=P BAR,8=N BAR,9=K ZERO BAR,10=K ZERO
	42	C
	43	IIT = ITT (IT)
	44	IF ( IIT .LE. 0 ) RETURN
	45	* Check the kinetic energy: no interaction below 50 MeV at present
	46	IF ( IIT .LE. 2 .AND. EKE .LT. 0.400D+00 .AND. AAA .GT. 1.5D+00 )
	47	& THEN
	48	IF ( IIT .EQ. 1 ) THEN
	49	SI = XSEPRO ( EKE, ZZZ, AAA )
	50	ELSE
	51	SI = XSENEU ( EKE, ZZZ, AAA )
	52	END IF
	53	IF ( EKE .GT. 0.200D+00 ) THEN
	54	CALL NIZL ( IT, AAA, EKE, PO, SI2, ZL2 )
	55	WEIGH1 = ( 0.400D+00 - EKE ) / 0.200D+00
	56	SI = WEIGH1 * SI + ( 1.D+00 - WEIGH1 ) * SI2
	57	END IF
	58	IF ( SI .GT. 0.D+00 ) THEN
	59	ZL = AVOGMB * AAA / SI
	60	ELSE
	61	SI = 0.D+00
	62	ZL = AINFNT
	63	END IF
	64	ELSE
65	CALL NIZL ( IT, AAA, EKE, PO, SI, ZL )
66	END IF
67	RETURN
68	END