# plot the absorption data with various normalizations
sigmaAbs = getresult("abs","sigma");
f=sigmaAbs.f;
sigma=sigmaAbs.sigma;
SP=sourcepower(f);
SI=sourceintensity(f);

plot(c/f*1e9,-sigma*SI,"wavelength (nm)","Absorbed Power (Watts)","","linewidth=2");
legend("");
plot(c/f*1e9,-sigma*SI/SP,"wavelength (nm)","Absorbed Power (a.u.)","Normalized to sourcepower","linewidth=2");
legend("");
plot(c/f*1e9,-sigma*1e18,sigma*0+pi*50^2,sigma*0+140^2,"wavelength (nm)","Absorption Corss-section (nm^2)","","linewidth=2");
legend("Cross section","Sphere cross sectional area","TFSF injection area" );