diff --git a/augur/generate.py b/augur/generate.py index 8c0e924..fc94158 100644 --- a/augur/generate.py +++ b/augur/generate.py @@ -232,14 +232,14 @@ def generate_sacc_and_stats(config): else: ells_here = ells # Trying to add bandpower windows - ells_aux = np.arange(0, np.max(ells_here)+1) - wgt = np.zeros((len(ells_aux), len(ells_here))) - for i in range(len(ells_here)): - in_win = (ells_aux > ell_edges[i]) & (ells_aux < ell_edges[i+1]) - wgt[in_win, i] = 1.0 - win = sacc.BandpowerWindow(ells_aux, wgt) + #ells_aux = np.arange(0, np.max(ells_here)+1) + #wgt = np.zeros((len(ells_aux), len(ells_here))) + #for i in range(len(ells_here)): + # in_win = (ells_aux > ell_edges[i]) & (ells_aux < ell_edges[i+1]) + # wgt[in_win, i] = 1.0 + #win = sacc.BandpowerWindow(ells_aux, wgt) S.add_ell_cl(key, tr1, tr2, - ells_here, np.zeros(len(ells_here)), window=win) + ells_here, np.zeros(len(ells_here))) #, window=win) # # For some reason add_ell_cl does not update the sacc file properly on the fly... # for i, ell in enumerate(ells_here): # S.add_data_point(key, (tr1, tr2), 0.0, ell=ell, error=1e30, window=win[i])