Skip to content

DM-46582: Add the TEx PSF metric in summary stats. #803

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.

Sign up for GitHub

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

Open
wants to merge 1 commit into
base: main
Choose a base branch
from
Open

DM-46582: Add the TEx PSF metric in summary stats. #803

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
150 changes: 150 additions & 0 deletions python/lsst/afw/image/_exposureSummaryStats.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -162,6 +162,84 @@ class ExposureSummaryStats(Storable):
zeroPoint, and readNoise.
"""

psfTE1e1: float = float('nan')
"""Per-exposure TE1e1 ~ <de1 de1> of PSF residual ellipticity, averaged over
theta [0,1] arcmin via treecorr KK correlation. Dimensionless; used to form the
full-survey TE1 metric.
"""

psfTE1e2: float = float('nan')
"""Per-exposure TE1e2 ~ <de2 de2> of PSF residual ellipticity, averaged over
theta [0,1] arcmin via treecorr KK correlation. Dimensionless; used to form the
full-survey TE1 metric.
"""

psfTE1ex: float = float('nan')
"""Per-exposure TE1ex ~ <de1 de2> of PSF residual ellipticity, averaged over
theta [0,1] arcmin via treecorr KK correlation. Dimensionless; used to form the
full-survey TE1 metric.
"""

psfTE2e1: float = float('nan')
"""Per-exposure TE2e1 ~ <de1 de1> of PSF residual ellipticity, averaged over
theta [5,100] arcmin via treecorr KK correlation. Dimensionless; used to form the
full-survey TE2 metric.
"""

psfTE2e2: float = float('nan')
"""Per-exposure TE2e2 ~ <de2 de2> of PSF residual ellipticity, averaged over
theta [5,100] arcmin via treecorr KK correlation. Dimensionless; used to form the
full-survey TE2 metric.
"""

psfTE2ex: float = float('nan')
"""Per-exposure TE2ex ~ <de1 de2> of PSF residual ellipticity, averaged over
theta [5,100] arcmin via treecorr KK correlation. Dimensionless; used to form the
full-survey TE2 metric.
"""

psfTE3e1: float = float('nan')
"""Per-exposure median-over-CCDs of TE3e1 ~ <de1 de1> of PSF residual
ellipticity, where each CCD uses theta within [0,5] arcmin bins. Dimensionless;
downstream pipelines take the 85th percentile over
images to evaluate TE3.
"""

psfTE3e2: float = float('nan')
"""Per-exposure median-over-CCDs of TE3e2 ~ <de2 de2> of PSF residual
ellipticity, where each CCD uses theta within [0,5] arcmin bins. Dimensionless;
downstream pipelines take the 85th percentile over
images to evaluate TE3.
"""

psfTE3ex: float = float('nan')
"""Per-exposure median-over-CCDs of TE3ex ~ <de1 de2> of PSF residual
ellipticity, where each CCD uses theta within [0,5] arcmin bins. Dimensionless;
downstream pipelines take the 85th percentile over
images to evaluate TE3.
"""

psfTE4e1: float = float('nan')
"""Per-exposure median-over-CCDs of TE4e1 ~ <de1 de1> of PSF residual
ellipticity, where each CCD uses theta within [5,20] arcmin bins. Dimensionless;
downstream pipelines take the 85th percentile over
images to evaluate TE4.
"""

psfTE4e2: float = float('nan')
"""Per-exposure median-over-CCDs of TE4e2 ~ <de2 de2> of PSF residual
ellipticity, where each CCD uses theta within [5,20] arcmin bins. Dimensionless;
downstream pipelines take the 85th percentile over
images to evaluate TE4.
"""

psfTE4ex: float = float('nan')
"""Per-exposure median-over-CCDs of TE4ex ~ <de1 de2> of PSF residual
ellipticity, where each CCD uses theta within [5,20] arcmin bins. Dimensionless;
downstream pipelines take the 85th percentile over
images to evaluate TE4.
"""

def __post_init__(self):
Storable.__init__(self)

Expand Down Expand Up @@ -417,6 +495,78 @@ def update_schema(cls, schema: Schema) -> None:
"skyBg, zeroPoint, and readNoise.",
units="mag",
)
schema.addField(
"psfTE1e1",
type="F",
doc="Per-exposure E1e1 ~ <de1 de1> of PSF residual ellipticity "
"over theta within [0,1] arcmin. Dimensionless; contributes to TE1.",
)
schema.addField(
"psfTE1e2",
type="F",
doc="Per-exposure E1e2 ~ <de2 de2> of PSF residual ellipticity "
"over theta within [0,1] arcmin. Dimensionless; contributes to TE1.",
)
schema.addField(
"psfTE1ex",
type="F",
doc="Per-exposure E1ex ~ <de1 de2> of PSF residual ellipticity "
"over theta within [0,1] arcmin. Dimensionless; contributes to TE1.",
)
schema.addField(
"psfTE2e1",
type="F",
doc="Per-exposure E2e1 ~ <de1 de1> of PSF residual ellipticity "
"over theta within [5, 100] arcmin. Dimensionless; contributes to TE2.",
)
schema.addField(
"psfTE2e2",
type="F",
doc="Per-exposure E2e2 ~ <de2 de2> of PSF residual ellipticity "
"over theta within [5, 100] arcmin. Dimensionless; contributes to TE2.",
)
schema.addField(
"psfTE2ex",
type="F",
doc="Per-exposure E2ex ~ <de1 de2> of PSF residual ellipticity "
"over theta within [5, 100] arcmin. Dimensionless; contributes to TE2.",
)
schema.addField(
"psfTE3e1",
type="F",
doc="Per-exposure median-over-CCDs of TE3e1 ~ <de1 de1> with "
"per-CCD theta within [0,5] arcmin. Dimensionless; used for TE3.",
)
schema.addField(
"psfTE3e2",
type="F",
doc="Per-exposure median-over-CCDs of TE3e2 ~ <de2 de2> with "
"per-CCD theta within [0,5] arcmin. Dimensionless; used for TE3.",
)
schema.addField(
"psfTE3ex",
type="F",
doc="Per-exposure median-over-CCDs of TE3ex ~ <de1 de2> with "
"per-CCD theta within [0,5] arcmin. Dimensionless; used for TE3.",
)
schema.addField(
"psfTE4e1",
type="F",
doc="Per-exposure median-over-CCDs of TE4e1 ~ <de1 de1> with "
"per-CCD theta within [5, 20] arcmin. Dimensionless; used for TE4.",
)
schema.addField(
"psfTE4e2",
type="F",
doc="Per-exposure median-over-CCDs of TE4e2 ~ <de2 de2> with "
"per-CCD theta within [5, 20] arcmin. Dimensionless; used for TE4.",
)
schema.addField(
"psfTE4ex",
type="F",
doc="Per-exposure median-over-CCDs of TE4ex ~ <de1 de2> with "
"per-CCD theta within [5, 20] arcmin. Dimensionless; used for TE4.",
)

def update_record(self, record: BaseRecord) -> None:
"""Write summary-statistic columns into a record.
Expand Down