; Play with the data in H-HAS-PI-18.001, AL-K-HETG-ACIS-S

;; restore, '/nfs/spectra/d6/ACIS_anal/H-HAS-PI-18.001_l1.idlsav'

;;sel = where(l1.grade EQ 0 OR l1.grade EQ 2 OR l1.grade EQ 3 OR $
;;		l1.grade EQ 4 OR l1.grade EQ 6)

; see the events...
;;plot, l1(sel).detx, l1(sel).dety, PSYM=3

; Make a crude histogram to count events in each order(s)
;;lin_hist, l1(sel).detx, 2.0, bins, counts
;;plot, bins, counts, PSYM=10

; Plot the MEG 3rd orders
!p.multi=[0,1,2]
;;pre_print_portrait
;;plot, l1(sel).detx, l1(sel).dety, PSYM=3, XRANGE=[54.,55.], YRANGE=[4.,5.]
;;plot, l1(sel).detx, l1(sel).dety, PSYM=3, XRANGE=[-55.,-54.], YRANGE=[-5.,-4.]
;;device, /close
;;$cp idl.ps h-has-pi-18.001.meg3.ps
;;set_plot, 'X'

; Plot the MEG 1st orders - pileup-deleted cores!!!
;;!p.multi=[0,1,2]
;;pre_print_portrait
;;plot, l1(sel).detx, l1(sel).dety, PSYM=3, $
;;	XRANGE=[52.,57.]*0.333, YRANGE=[2.,7.]*0.333
;;plot, l1(sel).detx, l1(sel).dety, PSYM=3, $
;;	XRANGE=[-57.,-52.]*0.333, YRANGE=[-7.,-2.]*0.333
;;device, /close
;;$cp idl.ps h-has-pi-18.001.meg1.ps
;;set_plot, 'X'

; Plot the HEG 1st orders - pileup-deleted cores!!!
;;!p.multi=[0,1,2]
;;pre_print_portrait
;;plot, l1(sel).detx, l1(sel).dety, PSYM=3, $
;;	XRANGE=[35.,38.], YRANGE=[-5.,-2.]
;;plot, l1(sel).detx, l1(sel).dety, PSYM=3, $
;;	XRANGE=[-38.,-35.], YRANGE=[2.,5.]
;;device, /close
;;$cp idl.ps h-has-pi-18.001.heg1.ps
;;set_plot, 'X'

; Plot the HEG 2nd order
!p.multi=[0,1,2]
pre_print_portrait
plot, l1(sel).detx, l1(sel).dety, PSYM=3, $
	XRANGE=[72.,73.5], YRANGE=[-7.5,-6.]
plot, l1(sel).detx, l1(sel).dety, PSYM=3, $
	XRANGE=[-73.5,-72.], YRANGE=[6.,7.5]
device, /close
$cp idl.ps h-has-pi-18.001.heg2.ps
set_plot, 'X'

!p.multi=0