#-------------------------------------------------------
# Note that the current response and background file and
# its systematic error is ONLY FOR AO-1 SIMULATION. 
# These are on-axis response, 
# and no off-axis response is available yet.
#-------------------------------------------------------

#----------------------------------------------
# Simulate 100 ks observation of 100mCrab source 
# (tested with Xspec 11.3.0t)
#----------------------------------------------

#xspec
#cpd /xw

# step-A: simulate your model spectra

model const powerlaw
   1e-1 -1
   2.1000      1.00000E-02  -3.0000      -2.0000       9.0000       10.000
   9.7000      1.00000E-02  0.00000E+00  0.00000E+00  1.00000E+24  1.00000E+24
data 1 hxd_pin_bgd_ao1.pha
back 1 hxd_pin_bgd_ao1.pha
resp 1 hxd_pin_ao1.rsp
data 2 hxd_gso_bgd_ao1.pha
back 2 hxd_gso_bgd_ao1.pha
resp 2 hxd_gso_ao1.rsp

fake
y
10mCrab
100mcrab_test_pin_100ks.fak
1e5,1
100mcrab_test_gso_100ks.fak
1e5,1


# default energy band of HXD
ignore 1:60.0-**
ignore 2:**-30.0 600.0-**


# step-B: have a glance at the HXD background level and then delete
data 3 hxd_pin_bgd_ao1.pha
resp 3 hxd_pin_ao1.rsp
data 4 hxd_gso_bgd_ao1.pha
resp 4 hxd_gso_ao1.rsp
ignore 3:60.0-**
ignore 4:**-30.0 600.0-**

setp e
pl ld 

data 4 none
data 3 none


# step-C: include correction file and see 
# 	at what energy band your source is detectable 
corf 1 hxd_pin_bgd_ao1.pha
corf 2 hxd_gso_bgd_ao1.pha
cornorm 1 0.1 ! BGD overestimate
cornorm 2 0.1 ! BGD overestimate
setp e
pl ld

# check spectra and ignore undetectable band
ig 2:250.0-**


# step-E: check the effect of background systematic error

cornorm 1 0.0
cornorm 2 0.0
fit
pl
cornorm 1 -0.1  ! BGD underestimate
cornorm 2 -0.1  ! BGD underestimate
fit
pl 
cornorm 1 0.1 ! BGD overestimate
cornorm 2 0.1 ! BGD overestimate
fit
pl 

exit
y