Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

DM-48591: Add multiprofit two-Gaussian PSF and Sersic columns #309

Merged
merged 1 commit into from
Feb 21, 2025

Conversation

taranu
Copy link
Contributor

@taranu taranu commented Feb 12, 2025

Checklist

When making changes to YAML files in the schemas directory:

  • If applicable, incremented the schema version number, following the guidelines in the contribution guide
  • Referred to the documentation on specific schemas for additional versioning information, change constraints, or tasks that may need to be performed, based on which schema is being updated

- name: g_psfModel_TwoGaussian_cen_x
"@id": "#Object.g_psfModel_TwoGaussian_cen_x"
datatype: double
description: Centroid (x-axis) of the two-Gaussian PSF model (g-band).
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This is relative to the declared center of the original PSF model? If so, it might be clearer to call it a delta_x or something.

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

No, it's in (multiprofit, not converted to afw/numpy convention) PSF image coordinates, so it's a lot of 17.5 +/- 0.005 values. They're not really useful for anything except as a diagnostic in case it's far off x.5, and I guess you can multiply them by 2 and round to get the PSF stamp size.

- name: z_psfModel_TwoGaussian_is_parent_flag
"@id": "#Object.z_psfModel_TwoGaussian_is_parent_flag"
datatype: boolean
description: Flag set for objects not fit because they are parents for the two-Gaussian PSF model (z-band).
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

If this is derivable from some other object table flag, I don't think we need to replicate it here, as long as we document that the PSF model Gaussian approximations and galaxy model fits don't run on parent objects.

description: Flag set for objects not fit because they are parents for the two-Gaussian PSF model (z-band).
fits:tunit:
# MultiProFit Sersic model
- name: sersic_cen_x
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I would be inclined to drop the pixel-coordinate centers from the standardized tables, since are redundant with the (ra, dec).

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I'm a little hesitant to do this because I derive the ra, dec uncertainties in a fairly stupid way. I just convert x + x_err, y + err to ra, dec and subtract the original ra_dec. This is probably fine in most cases but I don't think I was careful enough with edge cases like RA=360 or dec=-90. I suppose I should use whatever Clare implemented for the ra, dec errors instead.

datatype: double
description: Centroid (tract, y-axis) from the multiband Sersic model fit.
fits:tunit: pixel
- name: sersic_cen_ra
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think cen_ here is redundant; I'd drop it.

description: Centroid (declination) from the multiband Sersic model fit.
fits:tunit: deg
- name: sersic_reff_x
"@id": "#Object.sersic_reff_x"
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

As I think you're aware, I'm hoping to switch all ellipse-parameter fields to angular coordinates (DM-48425; @ctslater may have another related ticket) and standardize on one ellipse parameterization (probably {a, b, theta} or {r_det, q, theta}?). Since we don't know what we'll switching to yet, I don't think that ought to be part of this ticket, but it's coming.

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Yeah, I haven't polled galaxies about this yet. I'm a bit concerned that transforming the errors is not trivial, even ignoring covariances.

datatype: double
description: Sersic profile index parameter from the multiband Sersic model fit.
fits:tunit:
- name: sersic_cen_xErr
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I recommend dropping these pixel-coordinate uncertainties along with the pixel-coordinate values.

"@id": "#Object.sersic_rhoErr"
datatype: double
description: Error on ellipse rho (correlation coefficient) from the multiband Sersic model fit.
fits:tunit:
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Do you have any covariance uncertainties for the ellipse parameters + sersic index + flux?

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Not yet, unfortunately.

datatype: int
description: Number of iterations in the non-linear fit for the two-Gaussian PSF model (g-band).
fits:tunit:
- name: g_psfModel_TwoGaussian_chisq_red
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I'd just spell out _reduced; _red is too overloaded.

datatype: boolean
description: Flag set for failures with an unexpected or unknown cause for the multiband Sersic model.
fits:tunit:
- name: sersic_is_parent_flag
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

If this and sersic_not_primary_flag are redundant with the corresponding deblend_ and detect_ flags, I think we can drop them.

@taranu taranu force-pushed the tickets/DM-48591 branch 4 times, most recently from a55ce7e to 6e1ffda Compare February 20, 2025 02:20
@taranu taranu merged commit 8fa1191 into main Feb 21, 2025
10 checks passed
@taranu taranu deleted the tickets/DM-48591 branch February 21, 2025 00:06
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants