Difference between revisions of "SMICA propagation code"
Revision as of 08:35, 30 November 2018
The script "smica_coadd" allows users to combine a set of Planck frequency maps into a single component-separated map using SMICA weights. It has been used by the Planck team to propagate CMB and noise Monte Carlo realizations to obtain the set of SMICA simulations.
The script is written in Python2 and depends on several Python modules:
- healpy >= 1.7.4 (https://healpy.readthedocs.io/en/latest/) ;
- pyfits >= 3.4 (https://pythonhosted.org/pyfits/) ;
- numpy .
Ancillary data provided along with the script corresponds to:
- "transition_mask.fits.gz": definition of the Xhigh area (see Planck Collaboration IV. 2018 for a description of the SMICA pipeline);
- "weights_T_smica_R3.00_<X>.txt" with X in `["Xhigh", "Xfull"]`: harmonic weights to be used in temperature;
- "weights_EB_smica_R3.00.txt": harmonic weights to be used in polarization.
Input data consist of a set of nine Planck frequency maps, having the same properties as the Planck simulation maps.
- Units, maps should be in "K_CMB" units, except for the 545-GHz and 857-GHz maps, which should be in "MJy/sr" (conversion factors are explicitly given in the code).
- Beam convolution. the weights include a beam correction necessary to put all the maps at the same resolution. The beam window functions that we used are provided within the PLA, and can be used to generate new simulations.
- Calibration: Since LFI beam window functions have been de-normalized to account for calibration correction of real data, our weights also include a "re-normalization" of the LFI beam window function. Therefore, this can be ignored by users when generating new simulations.
The script provides a basic routine (named "my_freq_map_finder" L45) to retrieve input files from a given directory, which should be customized by users to match their needs.
smica_coadd input_directory output_directory 0 1000 smica_coadd input_directory output_directory 0 1000 --polar
will process sequencially the realizations numbered 0 to 1000 from the input directory. One can call the script in a parallel mode, where realization indices are split over processes or threads.
The outputs of a temperature run are FITS files of one column corresponding to intensity pixel values. The outputs of a polarization run are FITS files of two columns for Q and U stokes parameters. Units of the output maps are given in the UNITS header keyword.
Cosmic Microwave background
Planck Legacy Archive
(Planck) Low Frequency Instrument
Flexible Image Transfer Specification