sim_spectra.py
Now all the preparatory work is done and we can compute some spectra !
This is done with sim_spectra.py
https://github.com/simonsobs/PSpipe/blob/master/project/AnalyseSims/sim_spectra.py
and should be ran in the folder you used for the window function and mode coupling computation.
This bit is pretty fast and can be done in interactive mode
salloc -N 1 -C haswell -q interactive -t 01:00:00
export OMP_NUM_THREADS=64
python sim_spectra.py global_sims_all.dict
Using the same dictionnary file as used for the mode coupling computation, you will get a bunch of signal only spectra. This will produce a bunch of spectra:
You can plot the spectra using:
python sim_spectra_plot.py global_sims_all.dict
We will have to edit 'PATH_TO_THEORY_FILE' to the location of
https://github.com/simonsobs/PSpipe/blob/master/data/lensed_scalar_cls.npy
It will produce plot looking like
Producing the signal only spectra was important because first we want to make sure they look sensible, second we want to use them for the construction of the covariance matrix. We will go back to this in the sim_covariance wiki page. For now we want to continue and now compute the cross spectra between the two noisy splits that we have generated in combine_sims.py
For this we follow the same procedure as for the signal only spectra just replacing 'global_sims_all.dict' by 'global_sims_all_noise.dict' provided in:
salloc -N 1 -C haswell -q interactive -t 01:00:00
export OMP_NUM_THREADS=64
python sim_spectra.py global_sims_all_noise.dict
and produce plot like
