Skip to content

sim_window_and_bbl.py

Jump to bottom
thibautlouis edited this page Jul 3, 2019 · 15 revisions

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:

https://github.com/simonsobs/PSpipe/blob/master/project/AnalyseSims/NERSC_run/global_sims_all.dict

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

Clone this wiki locally