# use xspec to make an approximation to the XRCF carbon continuum source

# output selection
cpd /XTERM

# by eye this power law with some
# absorbtion looks like the predicted
# Kramers spectrum - however the absortion
# has edges that are not expect...
model wabs po
0.1
1.0
7.6

# dummy response
dummy 0.2 10.0

flux

plot model

# Make a fake data file
# make sure its not there 'cause we'll get a prompt other wise...
exec rm carbspec_heg.pha

# the .rmf and .arf files have been copied locally
# the long_fake_file_name forces a file name prompt...
# a simulated time of 135 seconds is used
fakeit none
heg_first_0518.rmf
ea_heg1_hybrid.arf
y
long_fake_file_name
carbspec_heg
135.0, 1.0, 0, 0

plot counts

# go to energy along the X axis
setplot energy
plot ldata delchi

# save the figure to a .ps file
cpd /PS
plot ldata delchi
cpd /XTERM
exec cp pgplot.ps carbspec_sim.ps
# leave a plot in the window
plot ldata delchi