0001 if nargin<3, sLinLog='log'; end
0002
0003 lambda=stEsurf.lambda;
0004 theta=stEsurf.theta;
0005
0006 indLambda=zeros(1,length(lambda0));
0007 sLegend=cell(1,length(lambda0));
0008 for ll=1:length(lambda0)
0009 ind=find(lambda==lambda0(ll),1);
0010 if isempty(ind)
0011 disp(['PwePlotThetaDep: lambda=' num2str(lambda0(ll)) ...
0012 ' cannot be found. Using lambda=' num2str(lambda(1)) ' instead']);
0013 ind=1;
0014 end
0015 indLambda(ll)=ind;
0016
0017 sLegend{2*ll-1}=['F_{E4}^0(\theta,\phi=0) - \lambda=' num2str(lambda(indLambda(ll)))];
0018 sLegend{2*ll} =['F_{E4}^0(\theta,\phi=90^\circ) - \lambda=' num2str(lambda(indLambda(ll)))];
0019 end
0020
0021
0022 F0E4phi0=(abs(stEsurf.Ecr(indLambda,:)).^2+abs(stEsurf.Ect(indLambda,:)).^2).^2;
0023 F0E4phi90=abs(stEsurf.Esf(indLambda,:)).^4;
0024
0025 figure('Name','Angle-dependent SERS EFs');
0026
0027 plot(theta/pi*180,F0E4phi0,theta/pi*180,F0E4phi90,':');
0028 if ~strcmpi(sLinLog,'linear')
0029 set(gca,'Yscale','log')
0030 end
0031 title(['Angle-dependent SERS EFs for r=' num2str(stEsurf.r0) 'nm']);
0032 xlabel('Theta [degrees]');
0033 ylabel('SERS EF');
0034 legend(sLegend);
0035