opt stat
opt logy
zone 2 2
Ntuple/Cuts $1 theta<90
Ntuple/Cuts $2 theta>1
Ntuple/Cuts $3 iconef=1
1d 5000 ' ' 100 0 90
1d 5001 ' ' 100 0 90
1d 5002 ' ' 100 0 90
Ntuple/Plot //LUN1/9901.theta $1.AND.$2.AND.$3 300000 1 300001 ' ' 5000
nt/cuts $3 iconef=0
Ntuple/Plot //LUN1/9901.theta $1.AND.$2.AND.$3 300000 1 300001 ' ' 5001
nbins=55
1d 6000 ' ' [nbins] -40 400
1d 6001 ' ' [nbins] -40 400
nt/cuts $3 iconef=1
opt liny
Ntuple/Plot //LUN1/9901.phi $1.AND.$2.AND.$3 300000 1 300001 ' ' 6000
nt/cuts $3 iconef=0
Ntuple/Plot //LUN1/9901.phi $1.AND.$2.AND.$3 300000 1 300001 ' ' 6001
zone 2 2
opt liny
fort/file 66 conetest.ps
meta 66
zone 1 1
add 5000 5001 5002 1.0 1.0
hist/plot 5002
hist/plot 5001 s
opt liny
zone 1 2
hist/plot 6000
hist/plot 6001
1d 6002 ' ' [nbins] -40 400
add 6000 6001 6002 1.0 1.0
normfac=$hinfo(6000,'ENTRIES')
nfac=1.0/[normfac]*[nbins]*9.0/11.0
add 6000 6000 6000 [nfac] 0.0 
normfac=$hinfo(6001,'ENTRIES')
nfac=1.0/[normfac]*[nbins]*9.0/11.0
add 6001 6001 6001 [nfac] 0.0 
normfac=$hinfo(6002,'ENTRIES')
nfac=1.0/[normfac]*[nbins]*9.0/11.0
add 6002 6002 6002 [nfac] 0.0 
zone 1 3
hist/set/max 6000 1.2
hist/set/min 6000 0.8
hist/set/max 6001 1.2
hist/set/min 6001 0.8
hist/set/max 6002 1.2
hist/set/min 6002 0.8
hist/plot 6000
hist/plot 6001
hist/plot 6002
close 66
shell ps2pdf conetest.ps