sim_window_and_bbl.py
The second step of the simulation analysis pipeline is to define window functions and compute their associated mode coupling matrices.
Note that this bit of computation takes some time (for high lmax and high nside) so you can either use an interactive node:
salloc -N 1 -C haswell -q interactive -t 03:00:00
export OMP_NUM_THREADS=64
python sim_window_and_bbl.py global_sims_all.dict
or the submission script sufficient for lmax=5000 provided in
https://github.com/simonsobs/PSpipe/blob/master/project/AnalyseSims/NERSC_run/MCM_PREP.sb
Note that if you use the interactive node, don't forget to export: export OMP_NUM_THREADS=64 before running the script, otherwise the job will be single-threaded and will take forever.
the sim_window_and_bbl.py script is provided in :
https://github.com/simonsobs/PSpipe/blob/master/project/AnalyseSims/sim_window_and_bbl.py
and an example of dictionary file 'global_sims_all.dict' in:
You will have to edit the dictionary file and specify 'PATH_TO_COMBINE_MAPS' and 'PATH_TO_BEAM', and also provide a path to the binning file you want to use. The combined maps are the ones you just generated with combine_sims.py, the beam folder is available at:
https://github.com/simonsobs/PSpipe/tree/master/project/AnalyseSims/NERSC_run/beam
and an example of binning file (delta L=30) is provided in:
https://github.com/simonsobs/PSpipe/blob/master/data/binningFile_400_30.dat
If the code run successfully it will fill the mcm and window folder with the appropriate files. Here is an example of window function
