#cbar.set_clim([-90,-50])
for radius_deg in [20, 40, 60, 80]:
r = np.sin(radius_deg * np.pi / 180.)
x = np.linspace(-r, r, 100)
beam_plot.plot(x, np.sqrt(r**2 - x**2), 'w-', linewidth=3)
beam_plot.plot(x, -np.sqrt(r**2 - x**2), 'w-', linewidth=3)
# Cuts plot initialization
cuts_plot = fig.add_subplot(gsr[1])
#receiver coordinates
ell = np.linspace(-np.pi / 2, np.pi / 2)
az = np.zeros_like(ell)
xticks = [-90, -60, -40, -20, 0, 20, 40, 60, 90]
beam_slice_E = get_interp_val(hpx_beam, ell, az)
beam_slice_E_err = get_interp_val(hpx_rms, ell, az)
beam_slice_H = get_interp_val(hpx_beam, ell, az + np.pi / 2)
beam_slice_H_err = get_interp_val(hpx_rms, ell, az + np.pi / 2)
cuts_E_line = cuts_plot.errorbar(ell * 180 / np.pi,
beam_slice_E,
beam_slice_E_err,
fmt='b.',
label='ECHO [E]')
cuts_H_line = cuts_plot.errorbar(ell * 180 / np.pi,
beam_slice_H,
beam_slice_H_err,
fmt='r.',
label='ECHO [H]')
cuts_plot.legend(loc='lower center')
#cbar.set_clim([-90,-50])
for radius_deg in [20,40,60,80]:
r = np.sin(radius_deg*np.pi/180.)
x = np.linspace(-r,r,100)
beam_plot.plot(x,np.sqrt(r**2-x**2),'w-',linewidth=3)
beam_plot.plot(x,-np.sqrt(r**2-x**2),'w-',linewidth=3)
# Cuts plot initialization
cuts_plot = fig.add_subplot(gsr[1])
#receiver coordinates
ell = np.linspace(-np.pi/2,np.pi/2)
az = np.zeros_like(ell)
xticks = [-90,-60,-40,-20,0,20,40,60,90]
beam_slice_E = get_interp_val(hpx_beam,ell,az)
beam_slice_E_err = get_interp_val(hpx_rms,ell,az)
beam_slice_H = get_interp_val(hpx_beam,ell,az+np.pi/2)
beam_slice_H_err = get_interp_val(hpx_rms,ell,az+np.pi/2)
cuts_E_line = cuts_plot.errorbar(ell*180/np.pi,
beam_slice_E,
beam_slice_E_err,
fmt='b.',
label='ECHO [E]')
cuts_H_line = cuts_plot.errorbar(ell*180/np.pi,
beam_slice_H,
beam_slice_H_err,
fmt='r.',
label='ECHO [H]')
cuts_plot.legend(loc='lower center')
ell = np.linspace(-np.pi/2,np.pi/2)
az = np.zeros_like(ell)
xticks = [-90,-60,-40,-20,0,20,40,60,90]
# Setup figure
fig = figure(figsize=(18,10))
gs = gridspec.GridSpec(2,4)
# Setup base_beam(rms)
base_beam = hp.read_map(opts.base_beam)
base_beam = hp.ma(base_beam,badval=badval)
base_rms = hp.read_map(opts.base_beam.replace('beam','rms'))
base_rms = hp.ma(base_rms,badval=badval)
base_slice_H = get_interp_val(base_beam,ell,az)
base_slice_H_err = get_interp_val(base_rms,ell,az)
base_slice_E = get_interp_val(base_beam,ell,az+np.pi/2)
base_slice_E_err = get_interp_val(base_rms,ell,az+np.pi/2)
base_coll = make_polycoll(base_beam,nsides=opts.nsides)
# Setup adj_beam(rms)
adj_beam = hp.read_map(opts.adjusted_beam)
adj_beam = hp.ma(adj_beam,badval=badval)
adj_rms = hp.read_map(opts.adjusted_beam.replace('beam','rms'))
adj_rms = hp.ma(adj_rms,badval=badval)
adj_slice_H = get_interp_val(adj_beam,ell,az)
adj_slice_H_err = get_interp_val(adj_rms,ell,az)
adj_slice_E = get_interp_val(adj_beam,ell,az+np.pi/2)
adj_slice_E_err = get_interp_val(adj_rms,ell,az+np.pi/2)
