# this script uses the "polarized transmission" group to 
# calculate the fraction of power that is polarized in the
# Ex direction and the fraction that is polarized in the Ey 
# direction. This is similar to measuring the transmission 
# after the passing the light through a polarizer.


runanalysis;
trans = getresult("polarized transmission","polarized_transmission");

plot(trans.f/1e12,trans.T_total,trans.T_Ex_polarized,trans.T_Ey_polarized,
     "frequency (THz)","transmission");
legend("total","X polarized","Y polarized");