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");