function nodallthwp(tt,hwp)
t=size(tt);
n=t(3);
sec=(n-mod(n,9))/9;
ttttt=tt;
for
kkk=1:sec
ts=ttttt(:,:,(kkk-1)*9+1:kkk*9);
hwpsec=hwp((kkk-1)*54+1:kkk*54);
m4deg(ts);
load difff
load sigdif
load summ
load
sigsum2
load f
fall(:,:,kkk)=f;
sumall(kkk,:)=summ;
sigsum2all(kkk,:)=sigsum2;
stokehwp(difff,sigdif,summ,sigsum2,hwpsec)
load sq
load su
load
sigsq2
load sigsu2
load qq
load uu
load
sigqq2
load siguu2
tb1(kkk,:)=[sq(1:2) sq(4:9) su(1:2) su(4:9) ];
tb2(kkk,:)=[sigsq2(1:2) sigsq2(4:9)
sigsu2(1:2) sigsu2(4:9)];
btb1(kkk,:)=[qq(1:2) qq(4:9) uu(1:2)
uu(4:9) ];
btb2(kkk,:)=[sigqq2(1:2) sigqq2(4:9) siguu2(1:2) siguu2(4:9)];
end
%medd=fopen('m4.dat','w');
%fprintf(medd,'%10.3f',sq,sigsq2,);
bb=tb1(:,:)';
sb=tb2(:,:)';
bbb=btb1(:,:)';
bsb=btb2(:,:)';
for
j=1:16
for jj=1:sec
isb(jj)=1/sb(j,jj);
bisb(jj)=1/bsb(j,jj);
end
qu(j)=bb(j,1:sec)*isb(1:sec)'/sum(isb);
bqu(j)=bbb(j,1:sec)*bisb(1:sec)'/sum(bisb);
bb(j,sec+1)=qu(j);
bbb(j,sec+1)=bqu(j);
sb(j,sec+1)=1/sum(isb);
bsb(j,sec+1)=1/sum(bisb);
end
for
j=1:16
x=0;
bx=0;
for jj=1:sec
x=x+(bb(j,jj)-bb(j,sec+1))^2/sb(j,jj);
bx=bx+(bbb(j,jj)-bbb(j,sec+1))^2/bsb(j,jj);
end
X2(j)=x/(sec-1);
bX2(j)=bx/(sec-1);
bb(j,sec+2)=X2(j);
bbb(j,sec+2)=bX2(j);
sb(j,sec+2)= sb(j,sec+1)*X2(j);
bsb(j,sec+2)= bsb(j,sec+1)*bX2(j);
bb(j,sec+3)=sqrt(sb(j,sec+2));
bbb(j,sec+3)=sqrt(bsb(j,sec+2));
squ2(j)= sb(j,sec+2);
end
for
jj=1:sec
bbb(17,jj)=mean(fall(:,5,jj));
bbb(18,jj)= std(fall(:,5,jj));
bbb(19,jj)=sumall(jj,5);
bbb(20,jj)=sqrt(sigsum2all(jj,5));
bbb(21,jj)=(bbb(20,jj)/bbb(19,jj))/((sqrt(bsb(4,jj))/bbb(4,jj)+sqrt(bsb(12,jj))/bbb(12,jj))/2);
end
tab(1,:)=bbb(4,:);tab(2,1:sec+2)=sqrt(bsb(4,:));
tab(3,:)=bbb(12,:);tab(4,1:sec+2)=sqrt(bsb(12,:));
tab(5,:)=bbb(19,:);tab(6,:)=bbb(20,:);
tab(7,:)=bb(4,:);tab(8,1:sec+2)=sqrt(sb(4,:));
tab(9,:)=bb(12,:);tab(10,1:sec+2)=sqrt(sb(12,:));
tab(11,:)=bbb(21,:);
tab(12,:)=bbb(17,:);tab(13,:)=bbb(18,:);
tabb=tab';
save
fall
save bb
save sb
save sumall
save bbb
save
tabb
col1= textread('nameUT1','%18c');
col1(sec+1,:)='mean ';
col1(sec+2,:)='Xr2 ';
col1(sec+3,:)='sigmean ( *Xr ) ';
medd=fopen('tabb.dat','w');
fprintf(medd,'
Q5 sigQ5 U5
sigU5 I sigI
q5 sigq5 u5
sigu5 (sigI/I)/(sigU/U) f
stdf\n');
for
i=1:sec+3
fprintf(medd,'%1c',col1(i,:));
for j=1:13
fprintf(medd,'%7.2e ',tabb(i,j));
end
fprintf(medd,'\n');
end
fclose(medd);
for l=1:8
p(l)=sqrt(qu(l)^2+qu(l+8)^2);
perr(l)=sqrt(squ2(l)*qu(l)^2/(qu(l)^2+qu(l+8)^2)+squ2(l+8)*qu(l+8)^2/(qu(l)^2+qu(l+8)^2));
aatan=atan(qu(l)/qu(l+8));
if qu(l+8)<0
aatan=aatan+pi;
end
phi(l)=mod(aatan*90/pi-57.7+180,180);
phierr(l)=90*perr(l)/pi/p(l);
end
p
perr
phi
phierr
p=[p([1,2]),0,p([3:8])];
perr=[perr([1,2]),0,perr([3:8])];
phi=[phi([1,2]),0,phi([3:8])];
phierr=[phierr([1,2]),0,phierr([3:8])];
for m=1:9
v=mod((m-mod(m,3.1))/3.1+1+(m-1)*3,9);
if v==0
v=9;
end
p1(v)=p(m);
perr1(v)=perr(m);
phi1(v)=phi(m);
phierr1(v)=phierr(m);
end
zzzn=[p1'
perr1' phi1' phierr1']
save zzzn