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