;-----
; %***% Obsid 457, Mkn 421
;
; COMMENTS:
;   The CXC evt1a.fits file was processed:
;     - fselect to eliminate grades 1,5,7 and keep
;       energy in 100eV to 10000eV range; "g02346Elt1"
;     - run Houck's destreak_s4.pro on result; "S4"
;     - run Glenn's adjust_pi_extrap.pro on result; "G110"
; adjust_pi_extrap, '/nfs/atum/d2/hetgs/obsid_457/acisf457_000N001_g02346Elt10S4_evt1a.fits', $
; '/nfs/atum/d2/hetgs/obsid_457/acisf457_000N001_g02346Elt10S4G110_evt1a.fits', $
; '/nfs/wiwaxia/d4/ASC/data/caldb/data/axaf/acis/bcf/gain/acisD1999-09-16gainN0002_110.fits'

;

; Process using TDET coordinates... to show in det coord.s
obs_anal, '/nfs/atum/d2/hetgs/obsid_457/'+ $
	'acisf457_000N001_g02346Elt10S4G110_evt1a.fits', $
	OUTPUT_PREFIX = 'obsid_457', OUT_DIR='Obs_out', $
	COORDS_NAME='TDET', $
	/OVERRIDE_ZERO, /ASPECT, /CENTER_FIT, $
	CD_WIDTH=0.5, $
	GRATING='HETG', DETECTOR='ACIS-S', $
	ORDER_SEL_ACCURACY = 0.5
	;;GAIN_FAC=[1.0,1.0,1.0,1.0, $   
	;;	1.54,1.11,1.18,1.00,1.11,1.18]

; Get roll angle for below from:
;; sol_plot, '/nfs/atum/d2/hetgs/obsid_457/pcadf457_000N001_asol1.fits'

; Process using Sky coordinates rotated to detector axes
obs_anal, '/nfs/atum/d2/hetgs/obsid_457/'+ $
	'acisf457_000N001_g02346Elt10S4G110_evt1a.fits', $
	OUTPUT_PREFIX = 'obsid_457', OUT_DIR='Obs_out', $
	COORDS_NAME='Sky', /FLIPY, ROLLAXAY=-54.204432, $
	/OVERRIDE_ZERO, /CENTER_FIT, $
	CD_WIDTH=0.25, $
	GRATING='HETG', DETECTOR='ACIS-S', $
	ORDER_SEL_ACCURACY = 0.3, $
	/EXPORT
	;;GAIN_FAC=[1.0,1.0,1.0,1.0, $
	;;	1.1,1.1,1.1,1.00,1.1,1.1], $


; Make image!
sel = where(oa_Etouse GT 0.3 AND oa_Etouse LT 10.0 AND $
	ABS(oa_AY - oa_aveAY) LT 300.0 AND $
	ABS(oa_AX - oa_aveAX) GT 2.0)   ; get rid of bright 
                                        ; (not piled-up) zero order core
Xtouse = oa_AX(sel)
Ytouse = oa_AY(sel)
Etouse = oa_Etouse(sel)

set_plot,'X'
xye_image, xtouse, ytouse, etouse, $
	xbin=6.0, ybin=6.0, gif_file='Obs_out/Obs457_image.gif', ONE=0.4