legend()
title('Residuals I')
subplot(3,1,2)
hist(aa[cut,1],bins=100,range=[-rr,rr],color='blue',alpha=0.2,label='All detectors:'+statstr(aa[cut,1]))
hist(bb[covmin > 40,1],bins=100,range=[-rr,rr],color='red',alpha=0.2,label='Per detectors:'+statstr(bb[cut,1]))
legend()
title('Residuals Q')
subplot(3,1,3)
hist(aa[cut,2],bins=100,range=[-rr,rr],color='blue',alpha=0.2,label='All detectors:'+statstr(aa[cut,2]))
hist(bb[cut,2],bins=100,range=[-rr,rr],color='red',alpha=0.2,label='Per detectors:'+statstr(bb[cut,2]))
legend()
title('Residuals U')
### Get Rho and epsilon
pars_all=mm.mixingmatrix_from_maps(output_maps_all,output_maps_all_spoiled,goodpix=covmin>40)
pars_det=mm.mixingmatrix_from_maps(output_maps_det,output_maps_det_spoiled,goodpix=covmin>40)
################### Make a MC
dtheta=15.
npointings=10000
calerror=1e-2
rhoT_all=[]
rhoQ_all=[]
rhoU_all=[]
epsTQ_all=[]
epsTU_all=[]
epsQU_all=[]
rhoT_det=[]
rhoQ_det=[]
rhoU_det=[]
