5 * Revision 1.1.1.1 1995/10/24 10:20:59 cernlib
9 #include "geant321/pilot.h"
10 *CMZ : 3.21/02 29/03/94 15.41.38 by S.Giani
12 SUBROUTINE STEEQ(XXH,IPV)
14 C *** CORRECTIONS FOR SINGLE PARTICLE SPECTRA (SHOWER PARTICLES) ***
15 C *** NVE 16-MAR-1988 CERN GENEVA
17 C ORIGIN : H.FESEFELDT (06-SEP-1985)
19 #include "geant321/mxgkgh.inc"
20 #include "geant321/s_curpar.inc"
21 #include "geant321/s_blank.inc"
22 #include "geant321/s_result.inc"
24 DIMENSION ALEM(7),VAL0(7)
26 C** DATA EM/ 4.0 , 10. , 15. , 20. , 30. , 100. , 1000./
27 DATA ALEM/ 1.40, 2.30 , 2.70 , 3.00 , 3.40 , 4.60 , 7.00 /
28 DATA VAL0/ 0.00, 0.40 , 0.48 , 0.51 , 0.54 , 0.60 , 0.65 /
32 IF ((IPART .NE. 7) .AND. (IPART .NE. 9)) GO TO 9999
33 IF (ABS(IPA(IPV)) .NE. 8) GO TO 9999
35 IF (RNDM(1) .GT. LOG(ATNO2)) GO TO 9999
36 EKW=PV(4,MXGKPV)-ABS(PV(5,MXGKPV))
38 IF (ALEKW .LE. ALEM(1)) GO TO 9999
40 C --- GET ENERGY BIN ---
42 IF (ALEKW .LT. ALEM(I)) GO TO 2
47 C *** USE LINEAR INTERPOLATION OR EXTRAPOLATION BY Y=RC*X+B ***
51 DXNVE=ALEM(I2)-ALEM(I1)
52 DYNVE=VAL0(I2)-VAL0(I1)
54 BNVE=VAL0(I1)-RCNVE*ALEM(I1)