offset=0.0
CIC=1
e=fltarr(9,25)
s=fltarr(9,25)
                                                                                  
                                                                                  
set_plot,'ps'
device,filename='/afs/slac/u/ki/lances/idl/lsst2.ps'
!p.multi=[0,2,2]
!p.multi=0
loadct,39
device,/color
device,ysize=8.5,/inches
device,xsize=7.0,/inches
device,yoffset=1.0,/INCHES
                                                                                  
                                                                                  
;for q=0,8 do begin
                                                                                  
tt=0
CIC=1
cc=fltarr(15)
cc(1)=60
cc(2)=240
cc(3)=160
cc(4)=240
cc(5)=60
cc(6)=240
cc(7)=60
cc(8)=240
cc(9)=60
cc(10)=240
cc(11)=60
cc(12)=240
cc(13)=200
cc(14)=200
                                                                                  
                                                                                  
;Run rdfloat.pro in my directory
;.com /afs/slac/u/ki/lances/idl/rdfloat.pro       
;@/afs/slac/u/ki/lances/idl/batch.pro
                                                                           
;if q lt 10 then filename='output_'+string(q,format='(I1)')
;if q ge 10 then filename='output_'+string(q,format='(I2)')
filename='/afs/slac.stanford.edu/u/ki/lances/lsst_sim/output_0'
rdfloat,filename,x,y,z,vx,vy,vz,xo,yo,zo,vxo,vyo,vzo,n
;if q ne 0 then goto,jump2
jump1:
                                                                                  
;r=sqrt(x*x+y*y)*abs(y)/y
;ro=sqrt(xo*xo+yo*yo)*abs(yo)/yo

if tt eq 0 then begin
plot,[-11000*7.0/8.5/2.,11000*7.0/8.5/2.],[-2000,9000],/NODATA,/xstyle,/ystyle
endif else begin
plot,[-2000*7.0/8.5/2.,2000*7.0/8.5/2.],[3400,5400],/NODATA,/xstyle,/ystyle
endelse
                                                                                  
for i=0L,N_elements(x)-1 do begin
                                                                                  
    uu=randomu(seed,1)
                                                                                  
    if (uu(0) lt (1000./float(N_elements(x))) ) then begin
                                                                                  
    oplot,[x(i),xo(i)],[z(i),zo(i)],color=cc(n(i))*CIC
                                                                                  
endif
                                                                                  
                                                                                  
endfor

aaa=[4180,1625,2533,800,800,550,530,385,385,365,365,320]
bbb=[2613,800, 0 ,   0,   0, 0,   0,  0,  0,  0,  0,  0]
ccc=[19200,6032,8577,2739,3803,5198,2058,5630,5630,3625,17192,0]
ddd=[0,6.034,6.034-1.822-4.436,6.035-1.822-4.436+4.038,6.035-1.822-4.436+4.038+0.068,6.035-1.822-4.436+4.038+0.068+0.508,6.035-1.822-4.436+4.038+0.068+0.508+0.03,6.035-1.822-4.436+4.038+0.068+0.508+0.03+0.368,6.035-1.822-4.436+4.038+0.068+0.508+0.03+0.368+0.016,6.035-1.822-4.436+4.038+0.068+0.508+0.03+0.368+0.016+0.045,6.035-1.822-4.436+4.038+0.068+0.508+0.03+0.368+0.016+0.045+0.06,6.035-1.822-4.436+4.038+0.068+0.508+0.03+0.368+0.016+0.045+0.06+0.025]*1e3
                                                                           
for iii=0,N_elements(aaa)-1 do begin
                                                                           
if ccc(iii) ne 0 then begin
xq=findgen(1000)/1000*(aaa(iii)-bbb(iii))+bbb(iii)
yq=-xq
zq=ddd(iii)+xq*xq/2/(ccc(iii))
oplot,xq,zq,thick=3
oplot,yq,zq,thick=3
xq=findgen(1000)/1000*aaa(iii)*2-aaa(iii)
oplot,xq,fltarr(N_elements(xq))+zq(N_elements(zq)-1),linestyle=2
xq=findgen(1000)/1000*bbb(iii)*2-bbb(iii)
oplot,xq,fltarr(N_elements(xq))+zq(0),linestyle=2
endif else begin
                                                                           
xq=findgen(1000)/1000*(aaa(iii)-bbb(iii))+bbb(iii)
yq=-xq
zq=fltarr(N_elements(xq))+ddd(iii)
oplot,xq,zq,thick=3
oplot,yq,zq,thick=3
endelse
                                                                           
endfor
                                                                           


END