5 * Revision 1.1.1.1 1995/10/24 10:20:58 cernlib
9 #include "geant321/pilot.h"
10 *CMZ : 3.21/02 29/03/94 15.41.38 by S.Giani
12 SUBROUTINE KMABS(NOPT)
14 C *** CHARGED KAON ABSORPTION BY A NUCLEUS ***
15 C *** NVE 04-MAR-1988 CERN GENEVA ***
17 C ORIGIN : H.FESEFELDT (09-JULY-1987)
19 C PRODUCTION OF A HYPERFRAGMENT WITH SUBSEQUENT DECAY
20 C PANOFSKY RATIO (K- P --> LAMBDA PI0/K- P --> LAMBDA GAMMA) = 3/2
22 #include "geant321/s_defcom.inc"
37 IF(ATNO2.GT.1.5) GOTO 30
56 PCM=RMASS(13)+RMASS(14)-RMASS(18)
58 SINT=SQRT(ABS(1.-COST*COST))
61 PV(1,2)=PCM*COST*SIN(PHI)
62 PV(2,2)=PCM*COST*COS(PHI)
72 1 PCM=PCM*PCM-RMASS(8)*RMASS(8)
75 PV(1,2)=PCM*COST*SIN(PHI)
76 PV(2,2)=PCM*COST*COS(PHI)
78 PV(4,2)=SQRT(PCM*PCM+RMASS(8)*RMASS(8))
90 *WRITE(NEWBCD,1002) XEND,YEND,ZEND,P,ISW
91 1002 FORMAT(1H0,'KAON ABSORBTION POSITION',3(2X,F8.2),2X,
92 * 'K MOMENTUM',2X,F8.4,2X,'INT CODE',2X,I2)
95 C** STAR PRODUCTION FOR PION ABSORPTION IN HEAVY ELEMENTS
102 CALL POISSO(BLACK,NBL)
104 *WRITE(NEWBCD,3003) NBL,TEX
105 IF(NT+NBL.GT.MXGKPV-2) NBL=MXGKPV-2-NT
110 IF(NT.EQ.MXGKPV-2) GOTO 31
113 EKIN1=-EKIN*LOG(RAN2)
117 IF(RNDM(2).GT.PNRAT) IPA1=14
120 SINT=SQRT(ABS(1.-COST*COST))
123 PV(5,NT)=ABS(RMASS(IPA1))
124 PV(6,NT)=RCHARG(IPA1)
126 PV(4,NT)=EKIN1+PV(5,NT)
127 PP=SQRT(ABS(PV(4,NT)**2-PV(5,NT)**2))
128 PV(1,NT)=PP*SINT*SIN(PHI)
129 PV(2,NT)=PP*SINT*COS(PHI)
131 IF(EKIN2.GT.TEX) GOTO 33
134 BLACK=0.50*LOG(ATNO2)
135 CALL POISSO(BLACK,NBL)
136 IF(NT+NBL.GT.MXGKPV-2) NBL=MXGKPV-2-NT
141 *WRITE(NEWBCD,3004) NBL,TEX
143 IF(NT.EQ.MXGKPV-2) GOTO 32
146 EKIN1=-EKIN*LOG(RAN2)
150 SINT=SQRT(ABS(1.-COST*COST))
154 IF(RAN.GT.0.60) IPA(NT)=-31
155 IF(RAN.GT.0.90) IPA(NT)=-32
156 PV(5,NT)=(ABS(IPA(NT))-28)*RMASS(14)
158 IF(IPA(NT).EQ.-32) PV(6,NT)=2.
160 PV(4,NT)=PV(5,NT)+EKIN1
161 PP=SQRT(ABS(PV(4,NT)**2-PV(5,NT)**2))
162 PV(1,NT)=PP*SINT*SIN(PHI)
163 PV(2,NT)=PP*SINT*COS(PHI)
165 IF(EKIN2.GT.TEX) GOTO 40
168 C** STORE ON EVENT COMMON
170 40 CALL GRNDM(RNDM,1)
175 IF(PV(7,I).LT.0.) PV(5,I)=-PV(5,I)
183 IF(NT.EQ.2) GO TO 9999
185 IF(NTOT.LT.NSIZE/12) GOTO 43
190 3003 FORMAT(1H ,I3,' BLACK TRACK PARTICLES PRODUCED WITH TOTAL KINETIC
191 * ENERGY OF ',F8.3,' GEV')
192 3004 FORMAT(1H ,I5,' HEAVY FRAGMENTS WITH TOTAL ENERGY OF',F8.4,' GEV')