myr=getresult("ONA_1","input 1/mode 1/gain");
x=myr.frequency;
y=myr.getattribute(myr.getattribute);
plot(x/1e12,y,"frequency (THz)","Transmission (dB)","");
holdon;

myr=getresult("ONA_2","input 1/mode 1/gain");
x=myr.frequency;
y=myr.getattribute(myr.getattribute);
plot(x/1e12,y,"frequency (THz)","Transmission (dB)","");


myr=getresult("ONA_3","input 1/mode 1/gain");
x=myr.frequency;
y=myr.getattribute(myr.getattribute);
plot(x/1e12,y,"frequency (THz)","Transmission (dB))","");
legend("ONA_1 scattering analysis","ONA_2 impulse reponse, no delay comp,", "ONA_3 impulse reponse, with delay comp,");  
holdoff;