myr=getresult("ONA_1","input 2/mode 1/transmission");
x=myr.frequency;
y=myr.getattribute(myr.getattribute);
plot(x/1e12,abs(y)^2,"frequency (THz)","drop port transmission","");
holdon;

myr=getresult("ONA_2","input 2/mode 1/transmission");
x=myr.frequency;
y=myr.getattribute(myr.getattribute);
plot(x/1e12,abs(y)^2,"frequency (THz)","drop port transmission","");


myr=getresult("ONA_3","input 2/mode 1/transmission");
x=myr.frequency;
y=myr.getattribute(myr.getattribute);
plot(x/1e12,abs(y)^2,"frequency (THz)","drop port transmission)","");
setplot("x min",199.3);
setplot("x max",200.7);
setplot("y min",0);
setplot("y max",1.3);

legend("ONA_1: scattering analysis", "ONA_2: impulse reponse, inadequate time window", "ONA_3: impulse response, adequate time window");  
holdoff;