= psphot Extended Source Fit 20130712 tag =

This page discusses some experiments where we ran psphot on simulated images using the current ipp tag ipp-20130712

This can be compared to the results with the ipp trunk discussed [wiki:psphot_extended_source_fit here]

The 40 simulated images were provided last year by Peter Draper for use in testing the psphot kron radius and magnitude
code. They were constructed using E Bertin's SkyMaker program. 


=== Sersic Index ===
If we look at the sersic index for all galaxies we get the following histogram. The hisogram on the right shows the counts for pure disks and bulges. We see that the fits for bulges are clustered around 4 as expected while disks mostly have values near 1.

[[Image(sersic_index.jpg)]] [[Image(sersic_index_pure.jpg)]]

=== Sersic Fit Magnitudes ===

The difference between the simulated galaxys' magnitudes and the sersic fit magnitudes are shown in the following plot. 91% of the sources are contained in the band near zero. The reason for the other bands needs to be investigated.

[[Image(dmag_sersic.jpg)]]

If we consider only the pure bulge sources the fits are better. psphot under estimates the flux somewhat. If we look at pure disks we see the banding.

[[Image(dmag_sersic_bulges.jpg)]] [[Image(dmag_sersic_disks.jpg)]]

The bands away from dmag == 0 are primarily objects with measured sersic index > 4.5 and small effective disk length. 

[[Image(dmag_disks_v_dl.jpg)]]

=== Radii ===
If we compare the measured widths with the source sizes we see effects that are not understood.

==== disks ====
The following plots compare the measured semi-major axis with the disk scale length. The plot on the left shows all pure disks while the plot on the
right zooms in on the linear regime. We find that the measured width is a factor of 2 larger.

[[Image(disk_widths.jpg)]] [[Image(disk_widths_zoom.jpg)]]

==== bulges ====

If we compare the measured width to the bulge effective radius we find that the measured widths are a fraction smaller
than expected.

[[Image(bulge_widths.jpg)]] 

== Exponential Model ==

=== Exponential magnitudes ===
In general the magnitude fit with exponential model is slightly less than the expected value. However when the galaxy
sample is restricted to pure disks the magnitude differences are clustered around zero.

[[Image(dmag_exp_20130712.jpg)]] [[Image(dmag_exp_disks_20130712.jpg)]]

=== Radii ===

As with the sersice models, the measured width is about 2 x the modeled disk scale length.

[[Image(disks_widths_exp_20130712.jpg)]]

== de Vacouleur Model ==
The de Vacouleur model provides a good fit to elliptical galaxies. In our sample these are represented by the galaxies
with ht == -6 pure disks

{{{

    I(r) = I(0) exp[ -k * R ** (1 / 4) ]

psphot extended source parameters measured for PS_MODE_DEV
   EXT_INST_MAG
   EXT_WIDTH_MAJ
   EXT_WIDTH_MIN
   EXT_THETA
}}}

It looks like there is something wrong with the DEV models.The magnitude differences are all over the place and there was a problem
matching up the models with the pure bulge objects.

[[Image(dmag_dev.jpg)]] [[Image(dmag_dev_bulges.jpg)]]

[[Image(bulge_widths_dev.jpg)]]