load_pha("source_bin10.pi") 
show_all()

subtract()
plot_data()

ignore()
notice(1,7) 
plot_data()


set_model(xswabs.abs1*xspowerlaw.power1)
show_model()

guess(get_model())
guess(power1)
show_model()

thaw(abs1.nh)
show_model()

set_pileup_model(jdpileup.jdp)
print(get_pileup_model())
jdp.f.min=0.85
!dmkeypar acisf01618_evt2.fits EXPTIME echo+
jdp.ftime=3.2
jdp.fracexp=1   
jdp.nterms=5 
print(get_pileup_model())

set_method("neldermead")
set_stat("chi2datavar")
fit()

set_method("moncar")
fit()

set_method("neldermead")
set_proj_opt("fast", False)
proj()

plot_fit_delchi()

print(get_pileup_model())

save("pileup_fit.save")