CalMresfdaH.m

function [normI,phi,firstX,MresH,MresVH] = CalMresfdaH(n_and_dateX,Xamp,qemres)
%uses calculated sine and cosine instead of fda Cor96XZstuf signals with a shift by 22 samples for the Q Cor signal
% to get the Iand Q data
%set_path
clockstart=clock
%load lastMresV
[Xsnif, Zsnif] =getFDAreaderXZ(n_and_dateX);
[XCor, ZCor] =getFDAreaderXZCor(n_and_dateX);
siz=size(Xsnif)
MresH=zeros(224,96);
MresV=zeros(224,96);
MresVH=zeros(224,96);
MresHI=zeros(224,96);
MresVI=zeros(224,96);

MresHQ=zeros(224,96);
MresVQ=zeros(224,96);
MresHVQ=zeros(224,96);



%[Cor96Xstuf,Cor96Zstuf ] = SortCorrXZstuf( XZCor );

for i=1:96
    Cor96Xstuf(i,1:siz(2))=XCor(i,1:siz(2))-mean(XCor(i,:));
    
end

[ CXI,CXQ,firstX,lastX,lengthX]  = Corfast_fda2(Cor96Xstuf,Xamp,Xamp/20);
%[ CZI,CZQ,firstZ,lastZ,lengthZ]  = Corfast_fda2(Cor96Zstuf,Zamp,Zamp/20);


Xres=zeros(1,siz(2));
Zres=zeros(1,siz(2));
CX=zeros(1,siz(2));
CZ=zeros(1,siz(2));


for j=1:96
    j
    
    
    for k=1:224
        
        
        Xres(1:siz(2))=Xsnif(k,1:siz(2));
        Zres(1:siz(2))=Zsnif(k,1:siz(2));
        CX(1:siz(2))=CXI(j,1:siz(2));
        %CZ(1:siz(2))=CZI(j,1:siz(2));
        
        C=CX/sum(abs(CX));
        A=max(abs(CX));
        
        
        
        MresHI(k,j)=MresH(k,j)+Xres*C'/A;
        MresVI(k,j)=MresV(k,j)+Zres*C'/A;
        
        
        [ CX] = CXQ(j,1:siz(2));
        
        
        C=CX/sum(abs(CX));
        A=max(abs(CX));
        
        
        MresHQ(k,j)=MresH(k,j)+Xres*C'/A;
        MresVQ(k,j)=MresV(k,j)+Zres*C'/A;
        
        
    end
    
end

stdHI=std(MresHI);
stdHQ=std(MresHQ);
stdVI=std(MresVI);
stdVQ=std(MresVQ);
%calibration factor OK for.1A=32768
G=1e9*.2/pi*2*.5/32768;

normI=sqrt(stdHI.^2+stdHQ.^2)./stdHI;
normQ=sqrt(stdHI.^2+stdHQ.^2)./stdHQ;
phi=stdHQ./sqrt(stdHI.^2+stdHQ.^2);



for i=1:96
    if phi(i)< .8
        MresH(:,i)=MresHI(:,i)*normI(i)/G;
        MresVH(:,i)=MresVI(:,i)*normI(i)/G;
    else
        MresH(:,i)=MresHQ(:,i)*normQ(i)/G;
        MresVH(:,i)=MresVQ(:,i)*normQ(i)/G;
    end
end

time=time_label;
file=['TimeLabelMresX',time];

save(fullfile(qemres.datadir,file),'n_and_dateX' , 'MresH','MresVH')

clock-clockstart