Skip to content

Commit

Permalink
commenting out bandpower windows
Browse files Browse the repository at this point in the history
  • Loading branch information
Javier Sanchez authored and Javier Sanchez committed Feb 15, 2024
1 parent d187a17 commit d8b6537
Showing 1 changed file with 7 additions and 7 deletions.
14 changes: 7 additions & 7 deletions augur/generate.py
Original file line number Diff line number Diff line change
Expand Up @@ -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])
Expand Down

0 comments on commit d8b6537

Please sign in to comment.