;-----
; %***% Capella Obsid 1318 SIMULATION
;
; COMMENTS:
;   Testing MARX 3.0
;   Simulated along with background!
;   3mm Z offset, -0.33 Y offset

; ---
; Process using TDET coordinates... to show in det coord.s
sim_prefix = 'Capella_wbkg'
obs_anal, '/nfs/spectra/d0/MARX/test2.4/'+ $
	sim_prefix+'.fits', $
	OUTPUT_PREFIX = sim_prefix, OUT_DIR='MARX_out', $
	COORDS_NAME='TDET', /ASPECT, $
	/ZB, OVERR_WIDTH=200.0, $
	/CENTER_STREAK, $
	CD_WIDTH = 0.5, /LINE, $
	GRATING='HETG', DETECTOR='ACIS-S', $
	ORDER_SEL_ACCURACY = 0.3


;-----
; Analyze using sky X,Y coord.s
; Using MARX 3.0 simulation

; Get roll angle (ROLLAXAY in call below):
sol_plot, '/spectra/d0/MARX/test2.4/acis_asol.fits'

; Process using Sky coordinates rotated to detector axes
sim_prefix = 'Capella_wbkg'  
obs_anal, '/nfs/spectra/d0/MARX/test2.4/'+ $
	sim_prefix+'.fits', $
	OUTPUT_PREFIX = sim_prefix, OUT_DIR='MARX_out', $
	COORDS_NAME='Sky', /FLIPY, ROLLAXAY=-98.6, $
	/CENTER_STREAK, /ZB, OVERR_WIDTH=200.0, $
	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 z-o

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='MARX_out/'+sim_prefix+'_image.gif', ONE=0.4