pro puffy_events_radius_histogram ;Pretty much copied directly from Ernie's makehistogram_crosstalk.pro readcol,'C:\H2RG_15_EVENTS\100star.txt',radius,x_pix,y_pix,file_name,format='i,i,i,a' psfile='c:\H2RG_15_EVENTS\radius.ps' white=255 black=0 fg=0 _min=4 _max=10 _nbins=6 _binsize=(_max-_min)/(_nbins-1) theHist=histogram(radius,nbins=_nbins,max=_max,min=_min,locations=theVals) !x.minor = 3 !y.minor = 3 set_plot, 'ps' !p.font=-1 device, filename= psfile fsub = psfile ;filename to be written at bottom of plot ;set line and text characteristics for publication quality !P.CHARSIZE=1.2 !P.CHARTHICK=4. !X.THICK=8. !Y.THICK=8. !P.THICK=8. !P.MULTI=0 !P.TICKLEN=0.03 symsize=1.5 filenamecharsize=1.1 erase plot,theVals,theHist,psym=10,title="Cosmic Rays in H2RG-015",xtitle="Radius (pixels)",charsize=1.10*!P.CHARSIZE,ytitle='Number of Events', $ color=black,background=white,position=[0.15,0.15,0.8,0.9],ystyle=1,$ xrange=[0, _max],xstyle=1 ; xyouts,'hello','Total Exposure time',color=0,$ ; align=1.0,charsize=!P.CHARSIZE/1.25 device, /close set_plot, 'z' device, set_resolution=[8000,6000] device, set_font='Courier' ;fsub = jpgfile ;filename to be written at bottom of plot ; set line and text characteristics for publication quality !P.CHARSIZE=10 !P.CHARTHICK=15. !X.THICK=20. !Y.THICK=20. !P.THICK=20. symsize=10. filenamecharsize=8. erase plot,theVals,theHist,psym=10,title="Cosmic Rays in H2RG-015",xtitle="Radius (pixels)",charsize=1.10*!P.CHARSIZE,ytitle='Number of Events', $ color=black,background=white,position=[0.15,0.15,0.8,0.9],ystyle=1,$ xrange=[4, _max],xstyle=1 jpgimg = tvrd() write_jpeg, 'C:\H2RG_15_EVENTS\radius_histogram1.jpg', congrid(jpgimg, 1600, 1200, /interp, /center), quality=100 device,/close end