#=====================================================
# SCRIPT reproduces figure 3
#=====================================================
# S. Lin and S. Wang, 
# "High-throughput GaAs PIN electrooptic modulator with a 3-dB bandwidth of 9.6 GHz at 1.3 ?m,"
# Appl. Opt.  26, 1696-1700 (1987).
#======================================================
##
runsweep;
##
freq=linspace(3,24,40);
r3=getsweepresult("frequency_sweep","Out");
x=r3.getparameter("time");
y=r3.getattribute(r3.getattribute);
S21=matrix(40);
for (i=1:40)
{
    Vpp_optical=max(y(121200:end,i))-min(y(121200:end,i));
    S21(i)=Vpp_optical/0.001;     
}
##
S21_normalized=S21/S21(1);
plot(freq,S21_normalized,"Frequency (GHz)","Normalized |S21| ","Frequency response","precision=6,enhanced");
setplot("x min", 3);
setplot("x max", 23);
setplot("y max", 1);
setplot("y min", 0);
##
BW=1.9*c/(pi*8e-3*(6.6-3.41));
?"TWOMM 3dB BW = " + num2str(BW*1e-9) + " GHz";