7 CHARACTER frame*8,
proj*8,targ*8
15 common/himain2/
katt(130000,4),
patt(130000,4)
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)'
77 if(targ.eq.
'P' )
izt= 1
78 if(targ.eq.
'PBAR')
izt=-1
79 if(targ.eq.
'N' )
izt= 0
80 if(targ.eq.
'NBAR')
izt= 0
81 if(targ.eq.
'PI+' )
izt= 1
82 if(targ.eq.
'PI-' )
izt=-1
83 if(targ.eq.
'K+' )
izt= 1
84 if(targ.eq.
'K-' )
izt=-1
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)
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)