9 dimension dndpt(50), dndy(50)
50 write(6,*)
'===================================================='
51 write(6,*)
' Calculation of transverse momentum and rapidity '
52 write(6,*)
' distributions of charged particles in '
53 write(6,*)
'hh-, hA- and AA-collisions at fixed impact parameter'
54 write(6,*)
' Calculation will be performed in CM System '
55 write(6,*)
'===================================================='
58 write(6,*)
'Enter the energy per NN-collision (GeV)'
62 write(6,*)
'Enter a type of the "projectile" particle'
64 write(6,*)
' P proton, PBAR anti-proton,'
65 write(6,*)
' N neutron, NBAR anti-neutron,'
66 write(6,*)
' PI+ - positive pion, PI- negative pion,'
67 write(6,*)
' K+ positive kaon, K- negative kaon'
69 write(6,*)
' A - nucleus --------------------------'
86 write(6,*)
'Enter mass number and charge of the proj. nucleus'
91 write(6,*)
'Enter a type of the "target" particle (same notations)'
106 write(6,*)
'Enter mass number and charge of the target nucleus'
111 write(6,*)
'Enter number of events'
117 write(6,*)
' Simulation of interactions with'
119 write(6,*)
' Proj = ',
proj,
' and Targ = ',
targ
120 write(6,*)
' IAP =',
iap ,
' IAT =',
iat
121 write(6,*)
' IZP =',
izp ,
' IZT =',
izt
123 write(6,*)
' Reference frame - ',frame
124 write(6,*)
' ENERGY ',
efrm,
' (GeV)'
125 write(6,*)
' Number of generated events -',n_event
129 write(6,*)
'Enter Min. and Max. values of impact parameter (fm)'
138 write(6,*)
' Event # ',j,
' ------------------------------'
140 CALL
hijing(frame,bmin,bmax)
149 if(
katt(i,2).EQ.0 .OR.
katt(i,2).EQ.10) go to 1000
156 IF (ptr .GE. 10.0) go to 100
159 dndpt(ipt)=dndpt(ipt)+1.0/float(n_event)/0.2/2.0/ptr
162 IF(abs(
y) .GE. 10.0) go to 1000
165 dndy(iy)=dndy(iy)+1.0/float(n_event)/0.2/2.0
172 WRITE(6,2)0.2*(
k-1),dndpt(
k),dndy(
k)
173 2
format(1
x,f5.1,2(2
x,e11.4))