vector i1 i2 i3 i4 i5 tp x y z px py pz bx by bz wt ex ey ez sarc 500000
read\fast extreg9.dat\13 i1 i2 i3 i4 i5 tp x y z px py pz bx by bz wt ex ey ez sarc
set 
   linthk 2.5
   lintyp 2
   pchar 0
   histyp 1

clear
window\tile 1 3
color white
p=sqrt(px*px+py*py+pz*pz)
display `Long. phase space'
label\yaxis `No. of <mu>s'
label\xaxis `t (ns)'
display `Now t ns '
tmin=min(tp)
tmax=max(tp)
copy tp tp2 iff (i3=2)
period=1e3/201.25 ! nano sec
tp2=tp2*1e9 ! nano sec
p=sqrt(px*px+py*py+pz*pz)
copy p p2 iff (i3=2)
copy wt wt2 iff (i3=2)
bin\nbins\weight wt2 tp2 qp ep 100
qmin=min(qp)
qp=qp-qmin
statistics\weights wt2 tp2
statistics\weights wt2 p
!show qp
!show ep
show tp2
show p2
copy wt2 r iff (p2<=0.35)
copy tp2 q iff (p2<=0.35)
copy p2  w iff (p2<=0.35)
qnew=q+1.5
qq=mod(qnew,period) ! nano second
display ` STATISTICS time length; 0.1 ns per division'
statistics\weights r qq
display ` STATISTICS total momentum; 10 MeV per division'
statistics\weights r w
!qmin=min(qq)
!qp=qp-qmin
window 1
bin\nbins\weight r qq qp ep 50
graph\axes qp ep
graph\-axes qp ep
set
   histyp 0
   pchar -16

window 2
label\yaxis `p (MeV/c)'
label\xaxis `<tau> (ns)'
graph\axes qq w
graph\-axes qq w
window 3
bin2d\weights r qq w xi yi mc 50 35
!scale 0 5 0 350
label\yaxis `p (A.U.)'
label\xaxis `<tau> (A.U.)'
density\profile mc