paramprompt off
sherpa.dataplot.x_log = 1 
sherpa.dataplot.y_log = 1 

data 1 ae_hxd_pinbkg_20051105.pha
back 1 ae_hxd_pinbkg_20051105.pha
data 2 ae_hxd_gsobkg_20051105.pha
back 2 ae_hxd_gsobkg_20051105.pha

rsp[pin]
rsp[gso]
pin.rmf = ae_hxd_pinhxnom_20051104.rsp
gso.rmf = ae_hxd_gso_20051019.rsp
instrument 1 = pin
instrument 2 = gso

source 1,2 = const1d[c1]*powlaw1d[p1]
c1 integrate off
c1.c0 = 0.1
p1.gamma = 2.1
p1.ref = 1.0
p1.ampl.min = 0.0
p1.ampl.max = 100.0
p1.ampl = 9.7
freeze p1.ampl

fakeit 1 time = 1e5
fakeit 1
fakeit 2 time = 1e5
fakeit 2

ignore 1 energy 60.0:
ignore 2 energy :30.0,600.0:

subtract 1,2
oplot data 1 data 2 back 1 back 2
limits x 10 600
limits y 0.001 1   
c 2 blue
c 3 green
c 4 red 
redraw

ignore 2 energy 250.0:

unsubtract 1,2 
setdata 1 backscale = 1.0 
setdata 2 backscale = 1.0 
setback 1 backscale = 1.0 
setback 2 backscale = 1.0 
subtract 1,2
fit 
projection
oplot data 1 data 2 back 1 back 2
limits x 10 600
limits y 0.001 1   
c 2 blue
c 3 green
c 4 red
redraw

unsubtract 1,2
setback 1 backscale = 0.9  
setback 2 backscale = 0.9
subtract 1,2
fit
projection
oplot data 1 data 2 back 1 back 2
limits x 10 600
limits y 0.001 1   
c 2 blue
c 3 green
c 4 red
redraw

unsubtract 1,2 
setback 1 backscale = 1.1
setback 2 backscale = 1.1
subtract 1,2
fit
projection
oplot data 1 data 2 back 1 back 2
limits x 10 600
limits y 0.001 1   
c 2 blue
c 3 green
c 4 red
redraw

write data fake1.pha PHA
write data 2 fake2.pha PHA
save all suzaku.shp
write mdl "suzaku_mdl.fits"