for subdirname, subdir in zip(subdirnames, subdirs): figure() ccount = 0 patches = [] labels = [] for simname in simnames: simdat = data[simname] if simname == 'BCC': zlow = simdat[subdir]['redshifts'] < 0.36 massbin, ratiobin, ratioerr = bbb.summary2DMass(simdat[subdir], selection=zlow) fill_between(massbin, ratiobin + ratioerr, ratiobin - ratioerr, alpha=0.3, color=c[ccount], label='{0} z < 0.36'.format(simname)) patches.append(Rectangle((0, 0), 1, 1, fc=c[ccount], alpha=0.3)) labels.append('{0} z < 0.36'.format(simname)) ccount += 1 zlow = np.logical_and(simdat[subdir]['redshifts'] > 0.36, simdat[subdir]['redshifts'] < 0.5) massbin, ratiobin, ratioerr = bbb.summary2DMass(simdat[subdir], selection=zlow) fill_between(massbin,
##### ## Bias # for mcr in mcrs: figure() linecount = 0 ccount = 0 for strategy in strategies: axhline(1.0, c='DimGray', linestyle='--') simdat = mxxldata[mcr][strategy] massbin, ratiobin, ratioerr = bbb.summary2DMass(simdat) plot(massbin, ratiobin, color=c[ccount], linestyle=linestyle[linecount], label=strategy) ccount += 1 xlabel(r'Mass $M_{200}$', fontsize=16) ylabel(r'Ratio $M_{\mathrm{lensing}} / M_{\mathrm{true}}$', fontsize=16) axis([14.1, 15.2, 0.7, 1.3]) legend(loc='upper left') title('MC Relation: {0}'.format(mcr)) tight_layout() savefig('piscoproposal_{0}.png'.format(mcr))
if key not in data: print 'Missing %s' % ' '.join(key) continue simdat = data[key] print print key if simname == 'BCC': zlow = simdat['redshifts'] < 0.36 massbin, massedge, ratiobin, ratioerr = bbb.summary2DMass(simdat, selection = zlow) print len(simdat['redshifts'][zlow]), len(massbin) # errorbar(massbin, ratiobin, ratioerr, massedge, marker='None', color=c[ccount], linestyle='None') xpoints, y_low, y_high = fillerrs(massbin, massedge, ratiobin, ratioerr) # fill_between(massbin, ratiobin + ratioerr, ratiobin - ratioerr, alpha=0.3, color = c[ccount], label='{0} z < 0.36'.format(simname)) fill_between(xpoints, y_low, y_high, alpha=0.8, color = c[ccount], label='{0} z < 0.36'.format(simname)) patches.append(Rectangle((0, 0), 1, 1, fc=c[ccount], alpha=0.8)) labels.append('{0} z < 0.36'.format(simname)) ccount += 1