Description

The function allows one to construct a library of spectrum-weighted ULTRASAT PSFs based on
a) Pickles' stellar spectra determined by a grid of T_eff and log(g) OR
b) blackbody spectra on a grid of T OR
c) mixed library, containing any types of spectra in a linear sequence, determined by ultrasat.weightedPSFindex.

The output is two matlab objects containing a) spectrum-weighted PSFs and b) their flux containment radii for a range of two spectral parameters -- T_{eff} and log(g) [or just T for the black body case] -- and a range of radial distances of the source from the detector tile's inner corner. All the output PSF stamps are normalized to unity. The flux containment radius is calculated for a given level of flux containment.

Examples

% A cube of PSF stamps weighted with the Pickles' stellar spectra is produced and the PSF containment radius 
% is calculated for a 90% flux containment:
ultrasat.weightedPSF('ContainmentLevel',0.9);
% A cube of PSF stamps weighted with the black body spectra is produced and the PSF containment radius 
% is calculated for an 80% flux containment:
ultrasat.weightedPSF('Type','BB','ContainmentLevel',0.8);
% A cube of PSF stamps weighted with arbitrary spectra determined by ultrasat.weightedPSFindex is produced 
% and the PSF containment radius is calculated for a 90% flux containment:
ultrasat.weightedPSF('Type','all');
% A single 2D PSF stamp is produced for the given stellar class and source distance [in deg] from the tile's inner corner:
P = ultrasat.weightedPSF('Type','stellarclass','Class',{'o9','v'}, 'RDist', 5.5);