; This code computes and plots the response the firdem code to a log-normal distribution at ; several temperatures, using AIA data. .compile ../estimate_quantile.pro .compile ../firdem.pro .compile aia_firdem_wrapper.pro .compile ../firdem_get_pixel.pro .compile aia_firdem_test_1.pro set_plot,'ps' !p.font=0 device,file='plots/aia_dem_iterative_test_plots.eps',/encapsulated,/times device,/inches,xsize=10.0,ysize=6 !p.multi = [0,2,2] tempsin = [5e5,1e6,5e6,1e7] tempin=5e5 aia_firdem_test_1, tempsin[0], temps_in, dem_in, temps_out, dem_out, niter=ni, tr_struct=tr_struct,its=its1,chi2s=chi2s1 tempin=10e5 aia_firdem_test_1, tempsin[1], temps_in, dem_in, temps_out, dem_out, niter=ni, tr_struct=tr_struct,its=its2,chi2s=chi2s2 tempin=50e5 aia_firdem_test_1, tempsin[2], temps_in, dem_in, temps_out, dem_out, niter=ni, tr_struct=tr_struct,its=its3,chi2s=chi2s3 tempin=100e5 aia_firdem_test_1, tempsin[3], temps_in, dem_in, temps_out, dem_out, niter=ni, tr_struct=tr_struct,its=its4,chi2s=chi2s4 device,/close set_plot,'x' ;window,1 ;device,file='plots/aia_dem_iterative_test_chi2_plots.eps',/encapsulated !p.multi = [0,2,2] hist = histogram(chi2s1,nbins=40,locations=locs,min=0,max=10) plot,locs,hist,psym=10,title=string("Chi Squared histogram for T=",tempsin[0]),xtitle="Reduced Chi Squared",ytitle = "Counts" hist = histogram(chi2s2,nbins=40,locations=locs,min=0,max=10) plot,locs,hist,psym=10,title=string("Chi Squared histogram for T=",tempsin[1]),xtitle="Reduced Chi Squared",ytitle = "Counts" hist = histogram(chi2s3,nbins=40,locations=locs,min=0,max=10) plot,locs,hist,psym=10,title=string("Chi Squared histogram for T=",tempsin[2]),xtitle="Reduced Chi Squared",ytitle = "Counts" hist = histogram(chi2s4,nbins=40,locations=locs,min=0,max=10) plot,locs,hist,psym=10,title=string("Chi Squared histogram for T=",tempsin[3]),xtitle="Reduced Chi Squared",ytitle = "Counts" ;device,/close set_plot,'x' !p.multi=0 end