;-----
; %***% Obsid 1318, Capella
;
; COMMENTS:

; Process using TDET coordinates... to show in det coord.s
obs_anal, '/nfs/atum/d2/hetgs/obsid_1318/'+ $
	'acisf1318_000N001_g02346EsS4G110_evt1a.fits', $
	OUTPUT_PREFIX = 'obsid_1318', OUT_DIR='Obs_out', $
	COORDS_NAME='TDET', /ASPECT, $
	/OVERRIDE_ZERO, OVERR_WIDTH=200.0, $
	/CENTER_STREAK, $
	CD_WIDTH=0.5, /LINE, $
	GRATING='HETG', DETECTOR='ACIS-S', $
	ORDER_SEL_ACCURACY = 0.3


; Process using Sky coordinates rotated to detector axes

; Get roll angle for below from:
sol_plot, '/nfs/atum/d2/hetgs/obsid_1318/pcadf1318_000N001_osol1.fits'
;
; rot_angle =       -98.585223
;
obs_anal, '/nfs/atum/d2/hetgs/obsid_1318/'+ $
	'acisf1318_000N001_g02346EsS4G110_evt1a.fits', $
	OUTPUT_PREFIX = 'obsid_1318', OUT_DIR='Obs_out', $
	COORDS_NAME='Sky', /FLIPY, ROLLAXAY=-98.585223, $
	/OVERRIDE_ZERO, OVERR_WIDTH=200.0, $
	/CENTER_STREAK, $
	CD_WIDTH=0.5, /LINE, $
	GRATING='HETG', DETECTOR='ACIS-S', $
	ORDER_SEL_ACCURACY = 0.3, $
	/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/Obs1318_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/Obs1318_zo_image.gif', ONE=0.4

; Look at some lines:
; MEG 15A:
;;meas_line_angle, 2, 15.01
;;meas_line_angle, 2, 15.01, /VERB
; HEG 15A:
;;meas_line_angle, 1, 15.01
;;meas_line_angle, 1, 15.01, /VERB