> > There is a perl module written by Andy (Starlink::PSF) that takes the PSF
> > output from the KAPPA PSF command and calculates the integral for it. It
That's the total intergral, not the integral out to a specified radius.
> > After the next release (or at least after Brad branches for the release) can
> > anyone see an argument for not adding this? The perl code is below.
No. It may be harder to compute total fluxes for other functions if we
choose to add them to KAPPA:PSF.
> > my $integral = $PI * $amp * $widthsq * Math::Gsl::Sf::sf_gamma($k)
> > / ( $gamma * $const );
>
> If anyone is interested I seem to have some RFWS code around from
> EDRS/DIPS that also does this job.
I have a sersicfit.f written for
Our galaxy-photometry papers that make extensive use of the Sersic
function; and sersicfit.f includes generating the total luminosity.
The above formula looks to be the same (Eqn 7 of Young, Dr. Moan et al.
Astron. Astrophys. Suppl. Ser. 130, 173) after switching from FWHM to
sigma.
Malcolm
|