Loading [MathJax]/extensions/tex2jax.js
LSST Applications g04a91732dc+a3f7a6a005,g07dc498a13+5ab4d22ec3,g0fba68d861+870ee37b31,g1409bbee79+5ab4d22ec3,g1a7e361dbc+5ab4d22ec3,g1fd858c14a+11200c7927,g20f46db602+25d63fd678,g35bb328faa+fcb1d3bbc8,g4d2262a081+cc8af5cafb,g4d39ba7253+6b9d64fe03,g4e0f332c67+5d362be553,g53246c7159+fcb1d3bbc8,g60b5630c4e+6b9d64fe03,g78460c75b0+2f9a1b4bcd,g786e29fd12+cf7ec2a62a,g7b71ed6315+fcb1d3bbc8,g8048e755c2+a1301e4c20,g8852436030+a750987b4a,g89139ef638+5ab4d22ec3,g89e1512fd8+a86d53a4aa,g8d6b6b353c+6b9d64fe03,g9125e01d80+fcb1d3bbc8,g989de1cb63+5ab4d22ec3,g9f33ca652e+38ca901d1a,ga9baa6287d+6b9d64fe03,gaaedd4e678+5ab4d22ec3,gabe3b4be73+1e0a283bba,gb1101e3267+aa269f591c,gb58c049af0+f03b321e39,gb90eeb9370+af74afe682,gc741bbaa4f+7f5db660ea,gcf25f946ba+a750987b4a,gd315a588df+b78635c672,gd6cbbdb0b4+c8606af20c,gd9a9a58781+fcb1d3bbc8,gde0f65d7ad+5839af1903,ge278dab8ac+932305ba37,ge82c20c137+76d20ab76d,w.2025.11
LSST Data Management Base Package
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
sedterms.py
Go to the documentation of this file.
1# This file is part of fgcmcal.
2#
3# Developed for the LSST Data Management System.
4# This product includes software developed by the LSST Project
5# (https://www.lsst.org).
6# See the COPYRIGHT file at the top-level directory of this distribution
7# for details of code ownership.
8#
9# This program is free software: you can redistribute it and/or modify
10# it under the terms of the GNU General Public License as published by
11# the Free Software Foundation, either version 3 of the License, or
12# (at your option) any later version.
13#
14# This program is distributed in the hope that it will be useful,
15# but WITHOUT ANY WARRANTY; without even the implied warranty of
16# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
17# GNU General Public License for more details.
18#
19# You should have received a copy of the GNU General Public License
20# along with this program. If not, see <https://www.gnu.org/licenses/>.
21"""Configuration for SED terms in fgcmcal.
22Analogous to colorterms.
23"""
24
25from lsst.pex.config import Config, Field, ConfigDictField
26
27__all__ = ["Sedterm", "SedtermDict", "Sedboundaryterm", "SedboundarytermDict"]
28
29
31 """SED boundary term for a pair of bands.
32
33 The SED slope (in flux units) at the boundary between two bands is given by:
34
35 S = -0.921 * (primary - secondary) / (lambda_primary - lambda_secondary)
36
37 To construct a Sedboundaryterm, use keyword arguments:
38 Sedboundaryterm(primary=primaryBandName, secondary=secondaryBandName)
39
40 This is a subclass of Config. This follows the form of
41 `lsst.pipe.tasks.Colorterm`.
42 """
43 primary = Field(dtype=str, doc="name of primary band")
44 secondary = Field(dtype=str, doc="name of secondary band")
45
46
48 """A mapping of Sedboundaryterm name to Sedterm.
49
50 To construct a SedboundarytermDict use keyword arguments:
51 SedboundarytermDict(data=dataDict)
52 where dataDict is a Python dict of name: Sedterm
53 For example::
54
55 SedboundarytermDict(data={
56 'gr': Sedboundaryterm(primary="g", secondary="r"),
57 'ri': Sedboundaryterm(primary="r", secondary="i"),
58 })
59
60 This is a subclass of Config. This follows the form of
61 `lsst.pipe.tasks.ColortermDict`.
62 """
64 doc="Mapping of Sedboundary term name to Sedboundaryterm",
65 keytype=str,
66 itemtype=Sedboundaryterm,
67 default={},
68 )
69
70
72 """SED term for a single band.
73
74 The SED slope (in flux units) in the middle of a band is computed either
75 as an "interpolated" or "extrapolated" computation. See Burke et al. 2018
76 Appendix A (https://ui.adsabs.harvard.edu/abs/2018AJ....155...41B).
77
78 For interpolation, with a secondary term::
79
80 F'_nu ~ constant * (primaryTerm + secondaryTerm) / 2.0
81
82 For interpolation, without a secondary term::
83
84 F'_nu ~ constant * primaryTerm
85
86 For extrapolation::
87
88 F'_nu ~ primaryTerm + constant * (((lambda_primaryBand - lambda_secondaryBand) /
89 (lambda_primaryBand - lambda_tertiaryBand)) *
90 (primaryTerm - secondaryTerm))
91
92 where primaryTerm and secondaryTerm are names from a `SedboundarytermDict`, and
93 primaryBand, secondaryBand, and tertiaryBand are band names.
94
95 To construct a Sedterm, use keyword arguments::
96
97 Sedterm(primaryTerm=primaryTermName, secondaryTerm=secondaryTermName,
98 extrapolated=False, constant=1.0)
99
100 or::
101
102 Sedterm(primaryTerm=primaryTermName, secondaryTerm=secondaryTermName,
103 extrapolated=True, constant=1.0, primaryBand=primaryBandName,
104 secondaryBand=secondaryBandName, tertiaryBand=tertiaryBandName)
105
106 This is a subclass of Config. This follows the form of
107 `lsst.pipe.tasks.Colorterm`.
108 """
109 primaryTerm = Field(dtype=str, doc="Name of primary Sedboundaryterm")
110 secondaryTerm = Field(dtype=str, default=None, optional=True,
111 doc="Name of secondary Sedboundaryterm")
112 extrapolated = Field(dtype=bool, default=False, doc="Extrapolate to compute SED slope")
113 constant = Field(dtype=float, default=1.0, doc="Adjustment constant for SED slope")
114 primaryBand = Field(dtype=str, default=None, optional=True,
115 doc="Primary band name for extrapolation")
116 secondaryBand = Field(dtype=str, default=None, optional=True,
117 doc="Secondary band name for extrapolation")
118 tertiaryBand = Field(dtype=str, default=None, optional=True,
119 doc="Tertiary band name for extrapolation")
120
121 def validate(self):
122 Config.validate(self)
123 if self.extrapolated:
124 if self.primaryBand is None or \
125 self.secondaryBand is None or \
126 self.tertiaryBand is None:
127 raise RuntimeError("extrapolated requires primaryBand, secondaryBand, and "
128 "tertiaryBand are provided.")
129
130
132 """A mapping of bands to Sedterms.
133
134 To construct a SedtermDict use keyword arguments::
135
136 SedtermDict(data=dataDict)
137
138 where dataDict is a Python dict of band to Sedterm
139 For example::
140
141 SedtermDict(data={
142 'g': Sedterm(primaryTerm='gr', secondaryTerm='ri', extrapolated=True, constant=0.25,
143 primaryBand='g', secondaryBand='r', tertiaryBand='i'),
144 'r': Sedterm(primaryTerm='gr', secondaryTerm='ri', extrapolated=False)
145 })
146
147 This is a subclass of Config. This follows the form of
148 `lsst.pipe.tasks.ColortermDict`.
149 """
151 doc="Mapping of band name to Sedterm",
152 keytype=str,
153 itemtype=Sedterm,
154 default={},
155 )