lfast-etc 0.1.0__tar.gz

lfast_etc-0.1.0/LICENSE

@@ -0,0 +1,28 @@
+BSD 3-Clause License
+
+Copyright (c) 2026, Schmidt Observatories
+
+Redistribution and use in source and binary forms, with or without
+modification, are permitted provided that the following conditions are met:
+
+1. Redistributions of source code must retain the above copyright notice, this
+   list of conditions and the following disclaimer.
+
+2. Redistributions in binary form must reproduce the above copyright notice,
+   this list of conditions and the following disclaimer in the documentation
+   and/or other materials provided with the distribution.
+
+3. Neither the name of the copyright holder nor the names of its
+   contributors may be used to endorse or promote products derived from
+   this software without specific prior written permission.
+
+THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
+DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE
+FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
+SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
+CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
+OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.

lfast_etc-0.1.0/PKG-INFO

@@ -0,0 +1,135 @@
+Metadata-Version: 2.4
+Name: lfast-etc
+Version: 0.1.0
+Summary: Exposure Time Calculator for the LFAST telescope array (WIYN Bench spectrograph and CSHELL).
+Author-email: Gudmundur Stefansson <gstefansson@schmidtsciences.org>
+License-Expression: BSD-3-Clause
+Project-URL: Homepage, https://github.com/schmidt-observatories/lfast-etc
+Requires-Python: >=3.11
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: numpy>=2.0
+Requires-Dist: scipy>=1.10
+Requires-Dist: astropy>=5.3
+Requires-Dist: matplotlib>=3.7
+Provides-Extra: web
+Requires-Dist: flask>=3; extra == "web"
+Requires-Dist: bokeh>=3.4; extra == "web"
+Provides-Extra: docs
+Requires-Dist: sphinx>=7; extra == "docs"
+Requires-Dist: pydata-sphinx-theme>=0.15; extra == "docs"
+Requires-Dist: nbsphinx; extra == "docs"
+Requires-Dist: sphinx-design; extra == "docs"
+Requires-Dist: sphinx-copybutton; extra == "docs"
+Requires-Dist: myst-parser; extra == "docs"
+Requires-Dist: ipykernel; extra == "docs"
+Provides-Extra: dev
+Requires-Dist: pytest; extra == "dev"
+Requires-Dist: nbformat; extra == "dev"
+Requires-Dist: nbclient; extra == "dev"
+Requires-Dist: ipykernel; extra == "dev"
+Requires-Dist: jupyter; extra == "dev"
+Dynamic: license-file
+
+![Logo](docs/logo_schmidt_observatory_system.png)
+
+# lfast-etc
+
+Exposure Time Calculator for the **LFAST** telescope array at Kitt Peak (nominally 20 x 0.726 m
+unit telescopes feeding one spectrograph through an 18 um fiber combiner), with two instruments:
+the **WIYN Bench** spectrograph (visible; the 15 grating/order modes of the
+[WIYN configuration table](https://www.wiyn.org/Instruments/wiynbench.html) that reach 400 nm, scaled from the
+reference 600 @ 10.1 deg mode by `scripts/build_bench_presets.py`) and **CSHELL** (1.08-1.75 um, R up to 43,000).
+
+The number of telescopes, the fiber, and every instrument parameter are plain inputs, so the
+calculator follows the array design as it evolves. The physics is the per-resolution-element model
+of the reference LFAST ETC notebook (2026-10-12) with the Choi et al. (2026, SPIE 14147)
+throughput model; the test-suite reproduces that notebook to 1e-6.
+
+## Install
+
+Requires Python >= 3.11 and numpy >= 2.0.
+
+```sh
+git clone git@github.com:schmidt-observatories/lfast-etc.git
+cd lfast-etc
+pip install -e .            # core
+pip install -e ".[web]"     # + browser ETC
+pip install -e ".[docs]"    # + Sphinx docs
+```
+
+## Quick start
+
+```python
+from lfast_etc import LFASTETC, Source, Telescope
+
+etc = LFASTETC("bench_r1800", n_telescopes=20)          # presets: list_instruments()
+star = Source.blackbody(5800, mag=15, band="V")
+res = etc.exptime(star, snr=50)                         # per resolution element
+res.plot(); res.summary(); res.at(550.0); res.to_table()
+
+res = etc.snr(star, exptime=600, n_exp=3)               # SNR for 600 s in 3 reads
+res.plot_noise_budget(); res.plot_spectrum()            # what limits you; what the spectrum looks like
+etc.plot_throughput(); etc.plot_sky(); star.plot()
+etc.limiting_magnitude(star, exptime=900, snr=10, wavelength=(500, 650))
+etc.snr_grid(star, mags, exptimes, wavelength=(500, 650))   # sensitivity map
+
+# scale the array
+LFASTETC("cshell_r43000", telescope=Telescope(n_telescopes=40, fiber="18um"),
+         wavelength_min_nm=1081, wavelength_max_nm=1086)
+```
+
+Everything is per resolution element; SNR per pixel is lower by about `sqrt(pixels_per_resel)`.
+The bundled sky is grey time (moon up), so faint-end numbers are 0.3-0.5 mag pessimistic for
+dark time. The docs (`docs/`, build with `make html`) have a quick start, the concepts, every
+preset, the physics and its known limitations, and a plot gallery.
+
+## Browser ETC
+
+```sh
+pip install -e ".[web]"
+lfast-etc-web --port 5000     # http://127.0.0.1:5000  (JSON at /api/calc?...)
+```
+
+## Notebooks
+
+`notebooks/` holds executed examples, written as a walkthrough for someone new to the package.
+Start with 01 (the whole API, with checks and exercises at the end). 01-05 are rebuilt with
+`python scripts/build_notebooks.py`: quick start, sensitivity maps and array scaling, stellar
+metallicity, transient classification, exoplanet atmospheres (He 10830) with CSHELL. 06-10
+(`python scripts/build_science_notebooks.py`) compare science cases versus the array size
+(20 telescopes, 10 m equivalent = 200, 30 m equivalent = 1720): array-scaling overview,
+LSST-era time domain, Galactic archaeology, stellar activity and flares, CSHELL NIR science.
+The notebooks are generated from those scripts; edit the scripts, not the `.ipynb` files.
+
+## Docs
+
+```sh
+pip install -e ".[docs]"
+cd docs && make html         # needs pandoc on PATH for the notebooks
+```
+
+## Tests
+
+```sh
+pytest
+```
+
+CI (GitHub Actions) runs the tests on Python 3.11-3.13 and builds the docs on every push and
+pull request. Bugs and planned features are tracked in the
+[issues](https://github.com/schmidt-observatories/lfast-etc/issues), labelled by area and priority.
+
+## Layout
+
+```
+src/lfast_etc/        package (etc.py, source.py, telescope.py, instrument.py, sky.py, plotting.py)
+src/lfast_etc/data/   throughput curves, SkyCalc table, instrument presets (TOML)
+src/lfast_etc/web/    Flask + Bokeh browser ETC
+notebooks/            executed example notebooks
+docs/                 Sphinx site
+tests/                pytest, including the regression test against the reference notebook
+```
+
+## License
+
+BSD 3-Clause. See [LICENSE](LICENSE).

lfast_etc-0.1.0/README.md

@@ -0,0 +1,102 @@
+![Logo](docs/logo_schmidt_observatory_system.png)
+
+# lfast-etc
+
+Exposure Time Calculator for the **LFAST** telescope array at Kitt Peak (nominally 20 x 0.726 m
+unit telescopes feeding one spectrograph through an 18 um fiber combiner), with two instruments:
+the **WIYN Bench** spectrograph (visible; the 15 grating/order modes of the
+[WIYN configuration table](https://www.wiyn.org/Instruments/wiynbench.html) that reach 400 nm, scaled from the
+reference 600 @ 10.1 deg mode by `scripts/build_bench_presets.py`) and **CSHELL** (1.08-1.75 um, R up to 43,000).
+
+The number of telescopes, the fiber, and every instrument parameter are plain inputs, so the
+calculator follows the array design as it evolves. The physics is the per-resolution-element model
+of the reference LFAST ETC notebook (2026-10-12) with the Choi et al. (2026, SPIE 14147)
+throughput model; the test-suite reproduces that notebook to 1e-6.
+
+## Install
+
+Requires Python >= 3.11 and numpy >= 2.0.
+
+```sh
+git clone git@github.com:schmidt-observatories/lfast-etc.git
+cd lfast-etc
+pip install -e .            # core
+pip install -e ".[web]"     # + browser ETC
+pip install -e ".[docs]"    # + Sphinx docs
+```
+
+## Quick start
+
+```python
+from lfast_etc import LFASTETC, Source, Telescope
+
+etc = LFASTETC("bench_r1800", n_telescopes=20)          # presets: list_instruments()
+star = Source.blackbody(5800, mag=15, band="V")
+res = etc.exptime(star, snr=50)                         # per resolution element
+res.plot(); res.summary(); res.at(550.0); res.to_table()
+
+res = etc.snr(star, exptime=600, n_exp=3)               # SNR for 600 s in 3 reads
+res.plot_noise_budget(); res.plot_spectrum()            # what limits you; what the spectrum looks like
+etc.plot_throughput(); etc.plot_sky(); star.plot()
+etc.limiting_magnitude(star, exptime=900, snr=10, wavelength=(500, 650))
+etc.snr_grid(star, mags, exptimes, wavelength=(500, 650))   # sensitivity map
+
+# scale the array
+LFASTETC("cshell_r43000", telescope=Telescope(n_telescopes=40, fiber="18um"),
+         wavelength_min_nm=1081, wavelength_max_nm=1086)
+```
+
+Everything is per resolution element; SNR per pixel is lower by about `sqrt(pixels_per_resel)`.
+The bundled sky is grey time (moon up), so faint-end numbers are 0.3-0.5 mag pessimistic for
+dark time. The docs (`docs/`, build with `make html`) have a quick start, the concepts, every
+preset, the physics and its known limitations, and a plot gallery.
+
+## Browser ETC
+
+```sh
+pip install -e ".[web]"
+lfast-etc-web --port 5000     # http://127.0.0.1:5000  (JSON at /api/calc?...)
+```
+
+## Notebooks
+
+`notebooks/` holds executed examples, written as a walkthrough for someone new to the package.
+Start with 01 (the whole API, with checks and exercises at the end). 01-05 are rebuilt with
+`python scripts/build_notebooks.py`: quick start, sensitivity maps and array scaling, stellar
+metallicity, transient classification, exoplanet atmospheres (He 10830) with CSHELL. 06-10
+(`python scripts/build_science_notebooks.py`) compare science cases versus the array size
+(20 telescopes, 10 m equivalent = 200, 30 m equivalent = 1720): array-scaling overview,
+LSST-era time domain, Galactic archaeology, stellar activity and flares, CSHELL NIR science.
+The notebooks are generated from those scripts; edit the scripts, not the `.ipynb` files.
+
+## Docs
+
+```sh
+pip install -e ".[docs]"
+cd docs && make html         # needs pandoc on PATH for the notebooks
+```
+
+## Tests
+
+```sh
+pytest
+```
+
+CI (GitHub Actions) runs the tests on Python 3.11-3.13 and builds the docs on every push and
+pull request. Bugs and planned features are tracked in the
+[issues](https://github.com/schmidt-observatories/lfast-etc/issues), labelled by area and priority.
+
+## Layout
+
+```
+src/lfast_etc/        package (etc.py, source.py, telescope.py, instrument.py, sky.py, plotting.py)
+src/lfast_etc/data/   throughput curves, SkyCalc table, instrument presets (TOML)
+src/lfast_etc/web/    Flask + Bokeh browser ETC
+notebooks/            executed example notebooks
+docs/                 Sphinx site
+tests/                pytest, including the regression test against the reference notebook
+```
+
+## License
+
+BSD 3-Clause. See [LICENSE](LICENSE).

lfast_etc-0.1.0/pyproject.toml

@@ -0,0 +1,51 @@
+[build-system]
+requires = ["setuptools>=77", "wheel"]
+build-backend = "setuptools.build_meta"
+
+[project]
+name = "lfast-etc"
+version = "0.1.0"
+authors = [{ name = "Gudmundur Stefansson", email = "gstefansson@schmidtsciences.org" }]
+description = "Exposure Time Calculator for the LFAST telescope array (WIYN Bench spectrograph and CSHELL)."
+readme = "README.md"
+license = "BSD-3-Clause"
+requires-python = ">=3.11"
+dependencies = [
+  "numpy>=2.0",
+  "scipy>=1.10",
+  "astropy>=5.3",
+  "matplotlib>=3.7",
+]
+
+[project.optional-dependencies]
+web = ["flask>=3", "bokeh>=3.4"]
+docs = [
+  "sphinx>=7", "pydata-sphinx-theme>=0.15", "nbsphinx", "sphinx-design",
+  "sphinx-copybutton", "myst-parser", "ipykernel",
+]   # plus the pandoc binary on PATH
+dev = ["pytest", "nbformat", "nbclient", "ipykernel", "jupyter"]
+
+[project.scripts]
+lfast-etc-web = "lfast_etc.web.app:main"
+
+[project.urls]
+Homepage = "https://github.com/schmidt-observatories/lfast-etc"
+
+[tool.setuptools]
+package-dir = { "" = "src" }
+include-package-data = true
+
+[tool.setuptools.packages.find]
+where = ["src"]
+
+[tool.setuptools.package-data]
+"lfast_etc" = [
+  "data/throughput/*.npy",
+  "data/sky/*.fits",
+  "data/config/*.toml",
+  "web/templates/*.html",
+]
+
+[tool.pytest.ini_options]
+testpaths = ["tests"]
+pythonpath = ["tests"]

lfast_etc-0.1.0/setup.cfg

@@ -0,0 +1,4 @@
+[egg_info]
+tag_build = 
+tag_date = 0
+

lfast_etc-0.1.0/src/lfast_etc/__init__.py

@@ -0,0 +1,27 @@
+"""lfast_etc: Exposure Time Calculator for the LFAST telescope array.
+
+Quick start::
+
+    from lfast_etc import LFASTETC, Source
+    etc = LFASTETC("bench_r1800", n_telescopes=20)
+    src = Source.blackbody(5800, mag=15, band="V")
+    res = etc.exptime(src, snr=50)
+    res.plot()
+"""
+import importlib.metadata
+
+try:
+    __version__ = importlib.metadata.version("lfast-etc")
+except importlib.metadata.PackageNotFoundError:  # pragma: no cover
+    __version__ = "0.0.0"
+
+from .instrument import Instrument, list_instruments
+from .telescope import Telescope
+from .source import Source, BANDS
+from .sky import Sky
+from .etc import LFASTETC, ETCResult
+
+__all__ = [
+    "__version__", "LFASTETC", "ETCResult", "Instrument", "list_instruments",
+    "Telescope", "Source", "Sky", "BANDS",
+]

lfast_etc-0.1.0/src/lfast_etc/constants.py

@@ -0,0 +1,23 @@
+"""Physical constants (SI) and package data paths."""
+from pathlib import Path
+import numpy as np
+
+H_PLANCK = 6.62607015e-34   # J s
+C_LIGHT = 2.99792458e8      # m/s
+K_BOLTZMANN = 1.380649e-23  # J/K
+
+DATA_DIR = Path(__file__).parent / "data"
+
+
+def positive(name, value, allow_zero=False):
+    """Validate a scalar parameter at the public boundary; returns it as float."""
+    try:
+        v = float(value)
+    except (TypeError, ValueError):
+        raise ValueError(f"{name} must be a number, got {value!r}") from None
+    if not np.isfinite(v) or v < 0 or (v == 0 and not allow_zero):
+        raise ValueError(f"{name} must be {'>= 0' if allow_zero else '> 0'}, got {value!r}")
+    return v
+THROUGHPUT_DIR = DATA_DIR / "throughput"
+SKY_DIR = DATA_DIR / "sky"
+CONFIG_DIR = DATA_DIR / "config"

lfast_etc-0.1.0/src/lfast_etc/data/config/bench_1200_28p7_o1.toml

@@ -0,0 +1,18 @@
+# WIYN Bench spectrograph, 1200 @ 28.7 grating, order 1, blaze peak 7500 A.
+# Row of https://www.wiyn.org/Instruments/wiynbench.html; generated by
+# scripts/build_bench_presets.py from the bench_r1800 anchor (see that script).
+# Table coverage 1231 A, red-cable resolution 1.19 A at 2x2 binning.
+[instrument]
+name = "bench_1200_28p7_o1"
+description = "WIYN Bench, 1200 @ 28.7 deg, order 1, 688-812 nm, R~6910"
+wavelength_min_nm = 688.5
+wavelength_max_nm = 811.5
+resolving_power = 6910
+pixels_per_resel = 2.2    # dispersion direction only, as bench_r1800
+spectrograph_efficiency = 0.229  # grating peak efficiency 55% vs 72% for the 600 @ 10.1 anchor
+slit_throughput = 1.0
+qe = 0.75              # placeholder flat QE, as bench_r1800
+dark_current = 3.0e-4  # e-/s/pix
+read_noise = 3.7       # e-/pix
+full_well = 100000.0   # placeholder
+min_exptime = 1.0

lfast_etc-0.1.0/src/lfast_etc/data/config/bench_1200_28p7_o2.toml

@@ -0,0 +1,18 @@
+# WIYN Bench spectrograph, 1200 @ 28.7 grating, order 2, blaze peak 4000 A.
+# Row of https://www.wiyn.org/Instruments/wiynbench.html; generated by
+# scripts/build_bench_presets.py from the bench_r1800 anchor (see that script).
+# Table coverage 615 A, red-cable resolution 0.58 A at 2x2 binning.
+[instrument]
+name = "bench_1200_28p7_o2"
+description = "WIYN Bench, 1200 @ 28.7 deg, order 2, 400-431 nm, R~7560 (clipped at 400 nm)"
+wavelength_min_nm = 400.0
+wavelength_max_nm = 430.8
+resolving_power = 7561
+pixels_per_resel = 2.2    # dispersion direction only, as bench_r1800
+spectrograph_efficiency = 0.229  # grating peak efficiency 55% vs 72% for the 600 @ 10.1 anchor
+slit_throughput = 1.0
+qe = 0.75              # placeholder flat QE, as bench_r1800
+dark_current = 3.0e-4  # e-/s/pix
+read_noise = 3.7       # e-/pix
+full_well = 100000.0   # placeholder
+min_exptime = 1.0

lfast_etc-0.1.0/src/lfast_etc/data/config/bench_316_63p4_o9.toml

@@ -0,0 +1,18 @@
+# WIYN Bench spectrograph, 316 @ 63.4 grating, order 9, blaze peak 6563 A.
+# Row of https://www.wiyn.org/Instruments/wiynbench.html; generated by
+# scripts/build_bench_presets.py from the bench_r1800 anchor (see that script).
+# Table coverage 250 A, red-cable resolution 0.23 A at 2x2 binning.
+[instrument]
+name = "bench_316_63p4_o9"
+description = "WIYN Bench, 316 @ 63.4 deg, order 9, 644-669 nm, R~31280"
+wavelength_min_nm = 643.8
+wavelength_max_nm = 668.8
+resolving_power = 31284
+pixels_per_resel = 2.2    # dispersion direction only, as bench_r1800
+spectrograph_efficiency = 0.208  # grating peak efficiency 50% vs 72% for the 600 @ 10.1 anchor
+slit_throughput = 1.0
+qe = 0.75              # placeholder flat QE, as bench_r1800
+dark_current = 3.0e-4  # e-/s/pix
+read_noise = 3.7       # e-/pix
+full_well = 100000.0   # placeholder
+min_exptime = 1.0

lfast_etc-0.1.0/src/lfast_etc/data/config/bench_316_7p0_o1.toml

@@ -0,0 +1,18 @@
+# WIYN Bench spectrograph, 316 @ 7.0 grating, order 1, blaze peak 7500 A.
+# Row of https://www.wiyn.org/Instruments/wiynbench.html; generated by
+# scripts/build_bench_presets.py from the bench_r1800 anchor (see that script).
+# Table coverage 5290 A, red-cable resolution 5.72 A at 2x2 binning.
+[instrument]
+name = "bench_316_7p0_o1"
+description = "WIYN Bench, 316 @ 7.0 deg, order 1, 486-1014 nm, R~1440"
+wavelength_min_nm = 485.5
+wavelength_max_nm = 1014.5
+resolving_power = 1438
+pixels_per_resel = 2.2    # dispersion direction only, as bench_r1800
+spectrograph_efficiency = 0.292  # grating peak efficiency 70% vs 72% for the 600 @ 10.1 anchor
+slit_throughput = 1.0
+qe = 0.75              # placeholder flat QE, as bench_r1800
+dark_current = 3.0e-4  # e-/s/pix
+read_noise = 3.7       # e-/pix
+full_well = 100000.0   # placeholder
+min_exptime = 1.0

lfast_etc-0.1.0/src/lfast_etc/data/config/bench_316_7p0_o2.toml

@@ -0,0 +1,18 @@
+# WIYN Bench spectrograph, 316 @ 7.0 grating, order 2, blaze peak 3750 A.
+# Row of https://www.wiyn.org/Instruments/wiynbench.html; generated by
+# scripts/build_bench_presets.py from the bench_r1800 anchor (see that script).
+# Table coverage 2645 A, red-cable resolution 2.86 A at 2x2 binning.
+[instrument]
+name = "bench_316_7p0_o2"
+description = "WIYN Bench, 316 @ 7.0 deg, order 2, 400-507 nm, R~1440 (clipped at 400 nm)"
+wavelength_min_nm = 400.0
+wavelength_max_nm = 507.2
+resolving_power = 1438
+pixels_per_resel = 2.2    # dispersion direction only, as bench_r1800
+spectrograph_efficiency = 0.292  # grating peak efficiency 70% vs 72% for the 600 @ 10.1 anchor
+slit_throughput = 1.0
+qe = 0.75              # placeholder flat QE, as bench_r1800
+dark_current = 3.0e-4  # e-/s/pix
+read_noise = 3.7       # e-/pix
+full_well = 100000.0   # placeholder
+min_exptime = 1.0

lfast_etc-0.1.0/src/lfast_etc/data/config/bench_3300vph_o1.toml

@@ -0,0 +1,18 @@
+# WIYN Bench spectrograph, 3300 VPH grating, order 1, blaze peak 5120 A.
+# Row of https://www.wiyn.org/Instruments/wiynbench.html; generated by
+# scripts/build_bench_presets.py from the bench_r1800 anchor (see that script).
+# Table coverage 273 A, red-cable resolution 0.31 A at 2x2 binning.
+[instrument]
+name = "bench_3300vph_o1"
+description = "WIYN Bench, 3300 VPH, order 1, 498-526 nm, R~18110"
+wavelength_min_nm = 498.4
+wavelength_max_nm = 525.6
+resolving_power = 18108
+pixels_per_resel = 2.2    # dispersion direction only, as bench_r1800
+spectrograph_efficiency = 0.300  # no peak efficiency listed for the VPH gratings: anchor value kept
+slit_throughput = 1.0
+qe = 0.75              # placeholder flat QE, as bench_r1800
+dark_current = 3.0e-4  # e-/s/pix
+read_noise = 3.7       # e-/pix
+full_well = 100000.0   # placeholder
+min_exptime = 1.0

lfast_etc-0.1.0/src/lfast_etc/data/config/bench_400_4p2_o1.toml

@@ -0,0 +1,18 @@
+# WIYN Bench spectrograph, 400 @ 4.2 grating, order 1, blaze peak 5000 A.
+# Row of https://www.wiyn.org/Instruments/wiynbench.html; generated by
+# scripts/build_bench_presets.py from the bench_r1800 anchor (see that script).
+# Table coverage 4199 A, red-cable resolution 4.56 A at 2x2 binning.
+[instrument]
+name = "bench_400_4p2_o1"
+description = "WIYN Bench, 400 @ 4.2 deg, order 1, 400-710 nm, R~1200 (clipped at 400 nm) (not recommended by WIYN: hard to focus over the range)"
+wavelength_min_nm = 400.0
+wavelength_max_nm = 710.0
+resolving_power = 1202
+pixels_per_resel = 2.2    # dispersion direction only, as bench_r1800
+spectrograph_efficiency = 0.292  # grating peak efficiency 70% vs 72% for the 600 @ 10.1 anchor
+slit_throughput = 1.0
+qe = 0.75              # placeholder flat QE, as bench_r1800
+dark_current = 3.0e-4  # e-/s/pix
+read_noise = 3.7       # e-/pix
+full_well = 100000.0   # placeholder
+min_exptime = 1.0

lfast_etc-0.1.0/src/lfast_etc/data/config/bench_600_13p9_o1.toml

@@ -0,0 +1,18 @@
+# WIYN Bench spectrograph, 600 @ 13.9 grating, order 1, blaze peak 7500 A.
+# Row of https://www.wiyn.org/Instruments/wiynbench.html; generated by
+# scripts/build_bench_presets.py from the bench_r1800 anchor (see that script).
+# Table coverage 2724 A, red-cable resolution 2.86 A at 2x2 binning.
+[instrument]
+name = "bench_600_13p9_o1"
+description = "WIYN Bench, 600 @ 13.9 deg, order 1, 614-886 nm, R~2880"
+wavelength_min_nm = 613.8
+wavelength_max_nm = 886.2
+resolving_power = 2875
+pixels_per_resel = 2.2    # dispersion direction only, as bench_r1800
+spectrograph_efficiency = 0.292  # grating peak efficiency 70% vs 72% for the 600 @ 10.1 anchor
+slit_throughput = 1.0
+qe = 0.75              # placeholder flat QE, as bench_r1800
+dark_current = 3.0e-4  # e-/s/pix
+read_noise = 3.7       # e-/pix
+full_well = 100000.0   # placeholder
+min_exptime = 1.0

lfast_etc-0.1.0/src/lfast_etc/data/config/bench_600_13p9_o2.toml

@@ -0,0 +1,18 @@
+# WIYN Bench spectrograph, 600 @ 13.9 grating, order 2, blaze peak 3800 A.
+# Row of https://www.wiyn.org/Instruments/wiynbench.html; generated by
+# scripts/build_bench_presets.py from the bench_r1800 anchor (see that script).
+# Table coverage 1362 A, red-cable resolution 1.43 A at 2x2 binning.
+[instrument]
+name = "bench_600_13p9_o2"
+description = "WIYN Bench, 600 @ 13.9 deg, order 2, 400-448 nm, R~2910 (clipped at 400 nm)"
+wavelength_min_nm = 400.0
+wavelength_max_nm = 448.1
+resolving_power = 2913
+pixels_per_resel = 2.2    # dispersion direction only, as bench_r1800
+spectrograph_efficiency = 0.292  # grating peak efficiency 70% vs 72% for the 600 @ 10.1 anchor
+slit_throughput = 1.0
+qe = 0.75              # placeholder flat QE, as bench_r1800
+dark_current = 3.0e-4  # e-/s/pix
+read_noise = 3.7       # e-/pix
+full_well = 100000.0   # placeholder
+min_exptime = 1.0

lfast_etc-0.1.0/src/lfast_etc/data/config/bench_740vph_o1.toml

@@ -0,0 +1,18 @@
+# WIYN Bench spectrograph, 740 VPH grating, order 1, blaze peak 8798 A.
+# Row of https://www.wiyn.org/Instruments/wiynbench.html; generated by
+# scripts/build_bench_presets.py from the bench_r1800 anchor (see that script).
+# Table coverage 2152 A, red-cable resolution 2.68 A at 2x2 binning.
+[instrument]
+name = "bench_740vph_o1"
+description = "WIYN Bench, 740 VPH, order 1, 772-987 nm, R~3600"
+wavelength_min_nm = 772.2
+wavelength_max_nm = 987.4
+resolving_power = 3599
+pixels_per_resel = 2.2    # dispersion direction only, as bench_r1800
+spectrograph_efficiency = 0.300  # no peak efficiency listed for the VPH gratings: anchor value kept
+slit_throughput = 1.0
+qe = 0.75              # placeholder flat QE, as bench_r1800
+dark_current = 3.0e-4  # e-/s/pix
+read_noise = 3.7       # e-/pix
+full_well = 100000.0   # placeholder
+min_exptime = 1.0

lfast_etc-0.1.0/src/lfast_etc/data/config/bench_740vph_o2.toml

@@ -0,0 +1,18 @@
+# WIYN Bench spectrograph, 740 VPH grating, order 2, blaze peak 5000 A.
+# Row of https://www.wiyn.org/Instruments/wiynbench.html; generated by
+# scripts/build_bench_presets.py from the bench_r1800 anchor (see that script).
+# Table coverage 1076 A, red-cable resolution 1.34 A at 2x2 binning.
+[instrument]
+name = "bench_740vph_o2"
+description = "WIYN Bench, 740 VPH, order 2, 446-554 nm, R~4090"
+wavelength_min_nm = 446.2
+wavelength_max_nm = 553.8
+resolving_power = 4091
+pixels_per_resel = 2.2    # dispersion direction only, as bench_r1800
+spectrograph_efficiency = 0.300  # no peak efficiency listed for the VPH gratings: anchor value kept
+slit_throughput = 1.0
+qe = 0.75              # placeholder flat QE, as bench_r1800
+dark_current = 3.0e-4  # e-/s/pix
+read_noise = 3.7       # e-/pix
+full_well = 100000.0   # placeholder
+min_exptime = 1.0

lfast_etc-0.1.0/src/lfast_etc/data/config/bench_860_30p9_o1.toml

@@ -0,0 +1,18 @@
+# WIYN Bench spectrograph, 860 @ 30.9 grating, order 1, blaze peak 10800 A.
+# Row of https://www.wiyn.org/Instruments/wiynbench.html; generated by
+# scripts/build_bench_presets.py from the bench_r1800 anchor (see that script).
+# Table coverage 1680 A, red-cable resolution 1.64 A at 2x2 binning.
+[instrument]
+name = "bench_860_30p9_o1"
+description = "WIYN Bench, 860 @ 30.9 deg, order 1, 996-1164 nm, R~7220"
+wavelength_min_nm = 996.0
+wavelength_max_nm = 1164.0
+resolving_power = 7220
+pixels_per_resel = 2.2    # dispersion direction only, as bench_r1800
+spectrograph_efficiency = 0.271  # grating peak efficiency 65% vs 72% for the 600 @ 10.1 anchor
+slit_throughput = 1.0
+qe = 0.75              # placeholder flat QE, as bench_r1800
+dark_current = 3.0e-4  # e-/s/pix
+read_noise = 3.7       # e-/pix
+full_well = 100000.0   # placeholder
+min_exptime = 1.0

lfast_etc-0.1.0/src/lfast_etc/data/config/bench_860_30p9_o2.toml

@@ -0,0 +1,18 @@
+# WIYN Bench spectrograph, 860 @ 30.9 grating, order 2, blaze peak 6000 A.
+# Row of https://www.wiyn.org/Instruments/wiynbench.html; generated by
+# scripts/build_bench_presets.py from the bench_r1800 anchor (see that script).
+# Table coverage 840 A, red-cable resolution 0.77 A at 2x2 binning.
+[instrument]
+name = "bench_860_30p9_o2"
+description = "WIYN Bench, 860 @ 30.9 deg, order 2, 558-642 nm, R~8540"
+wavelength_min_nm = 558.0
+wavelength_max_nm = 642.0
+resolving_power = 8543
+pixels_per_resel = 2.2    # dispersion direction only, as bench_r1800
+spectrograph_efficiency = 0.271  # grating peak efficiency 65% vs 72% for the 600 @ 10.1 anchor
+slit_throughput = 1.0
+qe = 0.75              # placeholder flat QE, as bench_r1800
+dark_current = 3.0e-4  # e-/s/pix
+read_noise = 3.7       # e-/pix
+full_well = 100000.0   # placeholder
+min_exptime = 1.0

lfast_etc-0.1.0/src/lfast_etc/data/config/bench_860_30p9_o3.toml

@@ -0,0 +1,18 @@
+# WIYN Bench spectrograph, 860 @ 30.9 grating, order 3, blaze peak 4000 A.
+# Row of https://www.wiyn.org/Instruments/wiynbench.html; generated by
+# scripts/build_bench_presets.py from the bench_r1800 anchor (see that script).
+# Table coverage 560 A, red-cable resolution 0.52 A at 2x2 binning.
+[instrument]
+name = "bench_860_30p9_o3"
+description = "WIYN Bench, 860 @ 30.9 deg, order 3, 400-428 nm, R~8430 (clipped at 400 nm)"
+wavelength_min_nm = 400.0
+wavelength_max_nm = 428.0
+resolving_power = 8434
+pixels_per_resel = 2.2    # dispersion direction only, as bench_r1800
+spectrograph_efficiency = 0.271  # grating peak efficiency 65% vs 72% for the 600 @ 10.1 anchor
+slit_throughput = 1.0
+qe = 0.75              # placeholder flat QE, as bench_r1800
+dark_current = 3.0e-4  # e-/s/pix
+read_noise = 3.7       # e-/pix
+full_well = 100000.0   # placeholder
+min_exptime = 1.0