## The script will read the transmission at input 1
## for all three ONAs and plot them in the same figure

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


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



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

xlimlow = 199.949;
xlimhi = 200.051;
setplot("x min",xlimlow);
setplot("x max",xlimhi);
setplot("y min",0);
setplot("y max",0.6);

legend("ONA_1: scattering analysis", "ONA_2: impulse reponse, low sample rate", "ONA_3: impulse response, high sample rate");  
holdoff;