%example using Wayne script
[angr,pscat,beamc,snr, ang1,ang2,rings] = readpscat('DRE01.log', 'factory_zsc_1183.asc', 'ringarea_1183.asc');


c_lisst=median(beamc);
VSF=median(pscat);
loglog(angr*180/pi,VSF)
xlabel('angle [degrees]');
ylabel('VSF [m^{-1} sr^{-1}]');