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 PBANH(NOPT)
13 C *** ANTI PROTON ANNIHILATION AT REST ***
14 C *** NVE 04-MAR-1988 CERN GENEVA ***
16 C ORIGIN : H.FESEFELDT (09-JULY-1987)
18 C NOPT=0 NO ANNIHILATION
19 C NOPT=1 ANNIH.IN PI+ PI-
20 C NOPT=2 ANNIH.IN PI0 PI0
21 C NOPT=3 ANNIH.IN PI- PI0
22 C NOPT=4 ANNIH.IN GAMMA GAMMA
24 #include "geant321/s_defcom.inc"
28 DATA BRR/0.125,0.25,0.5/
33 PV(4,1)=ABS(RMASS(15))
44 IF(RAN.GT.BRR(1)) ISW=2
45 IF(RAN.GT.BRR(2)) ISW=3
46 IF(RAN.GT.BRR(3)) ISW=4
55 ELSEIF (ISW .EQ. 2) THEN
58 ELSEIF (ISW .EQ. 3) THEN
61 ELSEIF (ISW .EQ. 4) THEN
65 EK=RMASS(14)+ABS(RMASS(15))-RMNVE1-RMNVE2
68 IF(EK.LT.0.0001) EK=0.0001
70 EN=EK+0.5*(RMNVE1+RMNVE2)
71 S=AMAS*AMAS+RMASS(14)**2+2.0*RMASS(14)*EN
73 PCM=SQRT(ABS(EN*EN-RMNVE1*RMNVE2))
77 SINT=SQRT(ABS((1.-COST)*(1.+COST)))
78 PV(1,2)=PCM*SINT*COS(PHI)
79 PV(2,2)=PCM*SINT*SIN(PHI)
88 2 PV(4,2)=SQRT(PV(5,2)**2+PCM**2)
89 PV(4,3)=SQRT(PV(5,3)**2+PCM**2)
96 GOTO (21,22,23,24),ISW
117 IF(ATNO2.LT.1.5) GOTO 40
118 AFC=0.312+0.200*LOG(LOG(S))+S**1.5/6000.
119 C TARG=AFC*(ATNO2**0.33 -1.0)
120 CFA=0.025*((ATNO2-1.)/120.)*EXP(-(ATNO2-1.)/120.)
123 IF(TEX.LT.0.001) GOTO 445
124 BLACK=(1.5+1.25*TARG)*ENP(1)/(ENP(1)+ENP(3))
125 CALL POISSO(BLACK,NBL)
126 IF(IFIX(TARG)+NBL.GT.ATNO2) NBL=ATNO2-TARG
127 IF(NT+NBL.GT.MXGKPV-2) NBL=MXGKPV-2-NT
128 IF(NBL.LE.0) GOTO 445
133 IF(NT.EQ.MXGKPV-2) GOTO 441
134 IF(EKIN2.GT.TEX) GOTO 443
138 EKIN1=-EKIN*LOG(RAN1)-CFA*(1.+0.5*RAN2)
139 IF(EKIN1.LT.0.0) EKIN1=-0.010*LOG(RAN1)
142 IF(EKIN2.GT.TEX) EKIN1=TEX-(EKIN2-EKIN1)
143 IF(EKIN1.LT.0.) EKIN1=0.001
147 IF(RNDM(1).GT.PNRAT) IPA1=14
150 SINT=SQRT(ABS(1.-COST*COST))
153 PV(5,NT)=ABS(RMASS(IPA1))
154 PV(6,NT)=RCHARG(IPA1)
159 PV(4,NT)=EKIN1+PV(5,NT)
160 PP=SQRT(ABS(PV(4,NT)**2-PV(5,NT)**2))
161 PV(1,NT)=PP*SINT*SIN(PHI)
162 PV(2,NT)=PP*SINT*COS(PHI)
165 443 IF(ATNO2.LT.230.) GOTO 445
166 IF(EK.GT.2.0) GOTO 445
170 IF(EKA.GT.1.) EKA=EKA*EKA
171 IF(EKA.LT.0.1) EKA=0.1
175 IF(IPA(II).NE.-14) GOTO 444
178 PV(5,II)=ABS(RMASS(IPA1))
179 PV(6,II)=RCHARG(IPA1)
182 IF(KK.GT.IKA) GOTO 445
185 C** THEN ALSO DEUTERONS, TRITONS AND ALPHAS
188 IF(TEX.LT.0.001) GOTO 40
189 BLACK=(1.5+1.25*TARG)*ENP(3)/(ENP(1)+ENP(3))
190 CALL POISSO(BLACK,NBL)
191 IF(NT+NBL.GT.MXGKPV-2) NBL=MXGKPV-2-NT
197 IF(NT.EQ.MXGKPV-2) GOTO 442
198 IF(EKIN2.GT.TEX) GOTO 40
202 EKIN1=-EKIN*LOG(RAN1)-CFA*(1.+0.5*RAN2)
203 IF(EKIN1.LT.0.0) EKIN1=-0.010*LOG(RAN1)
206 IF(EKIN2.GT.TEX) EKIN1=TEX-(EKIN2-EKIN1)
207 IF(EKIN1.LT.0.) EKIN1=0.001
210 SINT=SQRT(ABS(1.-COST*COST))
214 IF(RAN.GT.0.60) IPA(NT+1)=-31
215 IF(RAN.GT.0.90) IPA(NT+1)=-32
217 PV(5,NT+1)=RMASS(INVE)
219 PV(6,NT)=RCHARG(INVE)
221 PV(8,NT)=ABS(IPA(NT))
224 PV(4,NT)=PV(5,NT)+EKIN1
225 PP=SQRT(ABS(PV(4,NT)**2-PV(5,NT)**2))
226 PV(1,NT)=PP*SINT*SIN(PHI)
227 PV(2,NT)=PP*SINT*COS(PHI)
233 IF(NT.EQ.2) GO TO 9999
235 IF(NTOT.LT.NSIZE/12) GOTO 43
241 *WRITE(NEWBCD,1001) XEND,YEND,ZEND,P,NCH,PP,PV(6,3)
242 1001 FORMAT(1H0,'PB ANNIHILATION AT REST POSITION',3(1X,F8.2),1X,
243 * 'PI MOMENTA,CHARGES',2(1X,F8.4,1X,F4.1))