;-----
; %***% Obsid 459, 3C 273
;
; COMMENTS:
;  - started with grating level 2 evt file
;  - fselect grade and energy select (g02346Elt10)
;  - destreak_s4.pro applied
;  - adjust_pi_dd.pro applied using -110 gains.

; Process using TDET coordinates... to show in det coord.s
obs_anal, '/nfs/spectra/d2/3c273/'+ $
	'acisf459_000N001_g02346Elt10s4G110_evt2.fits', $
	OUTPUT_PREFIX = 'obsid_459', OUT_DIR='Obs_out', $
	COORDS_NAME='TDET', $
	/OVERRIDE_ZERO, /ASPECT, /CENTER_STREAK, $
	CD_WIDTH=1.2,  $
	GRATING='HETG', DETECTOR='ACIS-S', $
	ORDER_SEL_ACCURACY = 0.4

; Get roll angle for below from:
;; sol_plot, '/nfs/spectra/d2/3c273/pcadf459aN000_asol1.fits'
;; sol_plot, '/nfs/spectra/d2/3c273/pcadf459bN000_asol1.fits'
; a) 27.3 ks : -68.000341
; b) 12.2 ks : -68.000318
; DOUBLE(27.3*(-68.000341) + 12.2*(-68.000318) )/ (27.3+12.2)
; = -68.000334
; Process using Sky coordinates rotated to detector axes
obs_anal, '/nfs/spectra/d2/3c273/'+ $
	'acisf459_000N001_g02346Elt10s4G110_evt2.fits', $
	OUTPUT_PREFIX = 'obsid_459', OUT_DIR='Obs_out', $
	COORDS_NAME='Sky', /FLIPY, ROLLAXAY=-68.000334, $
	/OVERRIDE_ZERO, /CENTER_STREAK, $
	CD_WIDTH=1.2, $
	GRATING='HETG', DETECTOR='ACIS-S', $
	ORDER_SEL_ACCURACY = 0.4, $
	/EXPORT

; Make standard image!
sel = where(oa_Etouse GT 0.3 AND oa_Etouse LT 10.0 AND $
	ABS(oa_AY - oa_aveAY) LT 350.0 )
;	ABS(oa_AX - oa_aveAX) GT 2.0)   ; if needed: 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, /LOG, $
	xbin=6.0, ybin=6.0, gif_file='Obs_out/Obs459_image.gif', ONE=0.4

; Make zero-order 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) LT 300.0 )
;	ABS(oa_AX - oa_aveAX) GT 2.0)   ; if needed: 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, /LOG, $
	xbin=1.5, ybin=1.5, gif_file='Obs_out/Obs459_zo_image.gif', ONE=0.4

; Super close-up for JET
sel = where(oa_Etouse GT 0.3 AND oa_Etouse LT 10.0 AND $
	ABS(oa_AY - oa_aveAY) LT 150.0 AND $
	ABS(oa_AX - oa_aveAX) LT 150.0 )
;	ABS(oa_AX - oa_aveAX) GT 2.0)   ; if needed: 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, /LOG, $
	xbin=0.5, ybin=0.5, gif_file='Obs_out/Obs459_jet_image.gif', ONE=0.4