RIN = getsweepdata("RIN_sweep", "source_RIN");
input_SNR = getsweepdata("RIN_sweep", "input_SNR");
reference_SNR = getsweepdata("RIN_sweep", "reference_SNR");
measured_NF = reference_SNR - input_SNR;
amplifier_NF = getnamed("AMP_1", "noise figure");

plot(RIN, measured_NF, "Source RIN (dB/Hz)", "Measured NF (dB)", "");
holdon;
plot(RIN, amplifier_NF*ones(length(RIN)), "Source RIN (dB/Hz)", "Measured NF (dB)", "", "pen=--");

# calculate shot noise RIN
f = getnamed("CWL_1", "frequency");
Pref = getnamed("CWL_1", "reference power");
RIN_sn = 10*log10(2*h*f/Pref);
plot(RIN_sn*ones(50), linspace(min(measured_NF), max(measured_NF), 50),  "Source RIN (dB/Hz)", "Measured NF (dB)", "", "pen=--");
holdoff;

legend("Measured NF", "Amplifier NF", "Shot noise RIN");