[u/mrichter/AliRoot.git] / GEANT321 / gheisha / cinema.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1995/10/24 10:21:03  cernlib
* Geant
*
*
#include "geant321/pilot.h"
*CMZ :  3.21/02 29/03/94  15.41.39  by  S.Giani
*-- Author :
      FUNCTION CINEMA(EK1)
C
C *** INELASTICITY IN NUCLEAR INTERACTIONS AS A FUNCTION ***
C *** OF ATOMIC NUMBER ATNO2 AND KINETIC ENERGY EK1 ***
C *** NVE 12-JUL-1988 CERN GENEVA ***
C
C ORIGIN : H.FESEFELDT (14-OCT-1987)
C
C THE FUNCTIONAL DEPENDENCE AND THE PARAMETERS HAVE BEEN OBTAINED
C BY STUDY OF VARIOUS NUCLEAR STRUCTURE MODELS.
C BUT: IT IS OF COURSE AN INTERPOLATION AS FUNCTION OF ATOMIC
C      NUMBER, FOR CERTAIN NUCLEI A DIFFERENT DESCRIPTION MAY BE
C      MORE ADEQUATE. DETAILED TESTS HAVE BEEN PERFORMED FOR
C      FE, CU, PB ,U AND SOME MIXTURES LIKE NAI, BGO, CONCRETE.
C
#include "geant321/mxgkgh.inc"
#include "geant321/s_curpar.inc"
#include "geant321/s_result.inc"
#include "geant321/s_mat.inc"
#include "geant321/limits.inc"
C
      CINEMA=0.0
      ND=IND+1
      ALA=LOG(ATNO2)
      ALEK1=LOG(EK1)
      SIG1=0.50
      SIG2=0.50
      EM=0.2390+0.0408*ALA**2
      IF (EM. GT. 1.0) EM=1.0
      CINEM=0.0019*ALA**3
      IF(CINEM.GT.0.15) CINEM=0.15
      IF (PARMAT(ND,10) .GE. 0.01) CINEM=CINEM*PARMAT(ND,10)
C
      IF (ALEK1 .GT. EM) GO TO 1
C
      CORR=-(ALEK1-EM)**2/(2.0*SIG1**2)
      IF (CORR .LT. EXPXL) CORR=EXPXL
      IF (CORR .GT. EXPXU) CORR=EXPXU
      DUM1=-EK1*CINEM
      DUM2=ABS(DUM1)
      DUM3=EXP(CORR)
      CINEMA=0.0
      IF (DUM2 .GE. 1.0) CINEMA=DUM1*DUM3
      IF ((DUM2 .LT. 1.0) .AND. (DUM3 .GT. 1.0E-10)) CINEMA=DUM1*DUM3
      GO TO 2
C
 1    CONTINUE
      CORR=-(ALEK1-EM)**2/(2.0*SIG2**2)
      IF (CORR .LT. EXPXL) CORR=EXPXL
      IF (CORR .GT. EXPXU) CORR=EXPXU
      DUM1=-EK1*CINEM
      DUM2=ABS(DUM1)
      DUM3=EXP(CORR)
      CINEMA=0.0
      IF (DUM2 .GE. 1.0) CINEMA=DUM1*DUM3
      IF ((DUM2 .LT. 1.0) .AND. (DUM3 .GT. 1.0E-10)) CINEMA=DUM1*DUM3
C
 2    CONTINUE
      IF (CINEMA .LT. -EK1) CINEMA=-EK1
C
      END
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.1 1995/10/24 10:21:03 cernlib
	6	* Geant
	7	*
	8	*
	9	#include "geant321/pilot.h"
	10	*CMZ : 3.21/02 29/03/94 15.41.39 by S.Giani
	11	*-- Author :
	12	FUNCTION CINEMA(EK1)
	13	C
	14	C * INELASTICITY IN NUCLEAR INTERACTIONS AS A FUNCTION *
	15	C * OF ATOMIC NUMBER ATNO2 AND KINETIC ENERGY EK1 *
	16	C * NVE 12-JUL-1988 CERN GENEVA *
	17	C
	18	C ORIGIN : H.FESEFELDT (14-OCT-1987)
	19	C
	20	C THE FUNCTIONAL DEPENDENCE AND THE PARAMETERS HAVE BEEN OBTAINED
	21	C BY STUDY OF VARIOUS NUCLEAR STRUCTURE MODELS.
	22	C BUT: IT IS OF COURSE AN INTERPOLATION AS FUNCTION OF ATOMIC
	23	C NUMBER, FOR CERTAIN NUCLEI A DIFFERENT DESCRIPTION MAY BE
	24	C MORE ADEQUATE. DETAILED TESTS HAVE BEEN PERFORMED FOR
	25	C FE, CU, PB ,U AND SOME MIXTURES LIKE NAI, BGO, CONCRETE.
	26	C
	27	#include "geant321/mxgkgh.inc"
	28	#include "geant321/s_curpar.inc"
	29	#include "geant321/s_result.inc"
	30	#include "geant321/s_mat.inc"
	31	#include "geant321/limits.inc"
	32	C
	33	CINEMA=0.0
	34	ND=IND+1
	35	ALA=LOG(ATNO2)
	36	ALEK1=LOG(EK1)
	37	SIG1=0.50
	38	SIG2=0.50
	39	EM=0.2390+0.0408ALA*2
	40	IF (EM. GT. 1.0) EM=1.0
	41	CINEM=0.0019ALA*3
	42	IF(CINEM.GT.0.15) CINEM=0.15
	43	IF (PARMAT(ND,10) .GE. 0.01) CINEM=CINEM*PARMAT(ND,10)
	44	C
	45	IF (ALEK1 .GT. EM) GO TO 1
	46	C
	47	CORR=-(ALEK1-EM)*2/(2.0SIG1**2)
	48	IF (CORR .LT. EXPXL) CORR=EXPXL
	49	IF (CORR .GT. EXPXU) CORR=EXPXU
	50	DUM1=-EK1*CINEM
	51	DUM2=ABS(DUM1)
	52	DUM3=EXP(CORR)
	53	CINEMA=0.0
	54	IF (DUM2 .GE. 1.0) CINEMA=DUM1*DUM3
	55	IF ((DUM2 .LT. 1.0) .AND. (DUM3 .GT. 1.0E-10)) CINEMA=DUM1*DUM3
	56	GO TO 2
	57	C
	58	1 CONTINUE
	59	CORR=-(ALEK1-EM)*2/(2.0SIG2**2)
	60	IF (CORR .LT. EXPXL) CORR=EXPXL
	61	IF (CORR .GT. EXPXU) CORR=EXPXU
	62	DUM1=-EK1*CINEM
	63	DUM2=ABS(DUM1)
	64	DUM3=EXP(CORR)
65	CINEMA=0.0
66	IF (DUM2 .GE. 1.0) CINEMA=DUM1*DUM3
67	IF ((DUM2 .LT. 1.0) .AND. (DUM3 .GT. 1.0E-10)) CINEMA=DUM1*DUM3
68	C
69	2 CONTINUE
70	IF (CINEMA .LT. -EK1) CINEMA=-EK1
71	C
72	END