18 dimension gb(101),
xb(101), dndp(50)
31 write(6,*)
'===================================================='
32 write(6,*)
' Calculation of transverse momentum distribution '
33 write(6,*)
' of charged pions in minimum bias '
34 write(6,*)
' hh-, hA- and AA-collisions '
35 write(6,*)
' Calculation will be performed in CM System '
36 write(6,*)
'===================================================='
39 write(6,*)
'Enter the energy per NN-collision (GeV)'
43 write(6,*)
'Enter a type of the "projectile" particle'
45 write(6,*)
' P proton, PBAR anti-proton,'
46 write(6,*)
' N neutron, NBAR anti-neutron,'
47 write(6,*)
' PI+ - positive pion, PI- negative pion,'
48 write(6,*)
' K+ positive kaon, K- negative kaon'
50 write(6,*)
' A - nucleus --------------------------'
67 write(6,*)
'Enter mass number and charge of the proj. nucleus'
72 write(6,*)
'Enter a type of the "target" particle (same notations)'
87 write(6,*)
'Enter mass number and charge of the target nucleus'
92 write(6,*)
'Enter number of events per each value of impact',
93 ,
' parameter (e.a. 10)'
100 write(6,*)
' Simulation of interactions with'
102 write(6,*)
' Proj = ',
proj,
' and Targ = ',
targ
103 write(6,*)
' IAP =',
iap ,
' IAT =',
iat
104 write(6,*)
' IZP =',
izp ,
' IZT =',
izt
106 write(6,*)
' Reference frame - ',frame
107 write(6,*)
' ENERGY ',
efrm,
' GeV'
108 write(6,*)
' Number of generated events per B interval -',n_event
118 dip=(bmax-bmin)/100.0
126 write(6,*)
'Inelastic X-section (mb) ',gbtot*
dip*10.*6.28
140 w_gb=(gb(ib)+gb(ib+1))/2.0/gbtot
148 IF(
np+
nt .EQ. 0)
THEN
150 gnorm=gnorm-gb(ib)/float(n_event)
158 IF(abs(
katt(
k,1)) .NE. 211) go to 150
166 IF(ptr .GE. 10.0) go to 150
168 dndp(ipt)=dndp(ipt)+w_gb/float(n_event)/0.2
178 dndp(i)=dndp(i)*gbtot/gnorm
179 if(dndp(i).ne.0.)
write(6,350)0.2*(i-1),dndp(i)
180 350
format(1
x,f5.1,2
x,e11.4)