PyBDSF gaussian fitting details

[Kincaidr]

Hello,

I have a simple question relating to the gaussian fitting procedure of PyBDSF. I am currently doing a radio simulation to derive source count corrections and for sanity check need to make sure my injected source properties (such as source size, peak flux, etc) match that with recovered. I ma using PyBDSF to extract the recovered sources properties.

The issue is that there is alot of scatter when comparing injected vs recovered particularly for major axis and total flux. When I locate these sources with large differences and do manual gaussian fitting with CASA viewer then the properties do match to my injected.

Si it seems the gaussian fitting PyBDSF is doing is not the same as with CASA. Do you know what it does different so that I can implement it into my simulation?

Thanks.