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=[]