# get the pump time signal
tpump = getdata("pump","time");
spump = getdata("pump","time_signal");

# get the probe time signal
tprobe = getdata("probe","time");
sprobe = getdata("probe","time_signal");

# plot them on the same graph
plotxy(tpump*1e12,(spump),tprobe*1e12,(sprobe),"t (ps)","signal (scaled)");
legend("pump","probe");

# plot the level populations
t = getdata("medium","t");
N0 = getdata("medium","storage_Ez_4");
N1 = getdata("medium","storage_Ez_5");
N2 = getdata("medium","storage_Ez_6");
N3 = getdata("medium","storage_Ez_7");

# plot the level populations vs time
plot(t*1e12,N0,N1,N2,N3,"t (ps)","population");
legend("N0","N1","N2","N3");
 
# visualize the transmission of the T monitor
T = getresult("T","T");
visualize(T);