source: ETALON/THZ_filters/2017_THz_FILTER/cross_filter_2d_g.m @ 733

 clear all; close all; clc
%% Double Cross filter 
n1=1;
n2=1;
n=1;%ref. indx of media between filters
Thetp=0;

a1=1e-4;
A1=0.53;


Thet=pi/180*0;
darr=(0:0.05:8)*1e-3;
for cnt2=1:length(darr)
g=2.9e-3;
a=1e-4;
b=1e-4;
d=darr(cnt2);
loss=0;
c=3e+8;
FreqArr=1e+9:1e+8:1e+11;
for cnt=1:length(FreqArr)

freq=FreqArr(cnt);

 Stot= S(freq, Thet, g, a, b, n1, n2)*T(freq,loss,d,n)* S(freq, Thet, g, a, b, n1, n2);
%  Stot=  S(freq, Thet, g, a, b, n1, n2);

TTE(cnt2,cnt)=abs(1./Stot(2,2))^2;
TTM(cnt2,cnt)=abs(1./Stot(4,4))^2;
end
end
% Ttot=sin(Thetp).^2.*TOP+cos(Thetp).^2.*TIP;
% figure
% 
% semilogy(FreqArr,TTE,'k')
% hold on
% 
% semilogy(FreqArr,TTM,'r')
figure
[hc,hc]=contourf(FreqArr*1e-9,darr/g,TTE);
set(hc,'LineStyle','none');
colorbar

set(gca,'fontsize',16)
ylabel('d/g, [1]')
xlabel('Frequency, [GHz]')
title('Double cross filter transmissivity')

% 3e+8./(6e-3*(1-cos((48:7:125)/180*pi)))*1e-9
% 
% 
% for cnt2=1:length(darr)
% 
% Px=(TTE(cnt2,:));
% fwm(cnt2)=length(Px(Px>0.5*max(Px)));
% Px2=(TTE(cnt2,1:600));
% 
% fwhm(cnt2)=length(Px2(Px2>0.5));
% 
% end
% figure
% hold on
% plot(darr,fwm)
% plot(darr,fwhm,'k')