#include "mt_rayor.i"
#include "image_profile.i"

mt_setup_system,"system_f15m.scm";

// First simple case to see system works with no scatter

mt_log,logflag=3;
mt_run,20.,0,0,dphot=0.1,no_scatter=1,no_mdeform=1,chat=2;
mt_analysis;
mt_detector;
im = mt_det_image(outfile="im15.fits");
disp,im,pane=0;

d = mt_sel_scatter(20,0.001);
plot,Anglesarr*206265,d; // convert radians to arcsec

// From now on the mirror deformations will not be considered
// More photons and now with scattering

mt_run,20.,0,0,dphot=0.5,no_scatter=0,chat=2;
mt_analysis;
mt_detector;
jm = mt_det_image(outfile="jm15.fits");
disp,jm,pane=1;