Add the ROSAT-to-eRASS1, FIR-radio and quiet-spotted searches and read HECATE's AGN flags correctly

LEDGER.md

@@ -270,7 +270,7 @@ is its docstring, and its calls are in the recipes of `README.md`.
 | Planets in step with their star's spin | NASA Exoplanet Archive composite table: 178 planets whose host carries a rotation period and whose orbital period is 1.5-20.5 times it, those at a whole multiple within 0.5 per cent against the half-integer windows and uniform fractions, each read for the source of the rotation period | 12 at whole multiples against 10.1 expected and 13 at half-integers (Poisson p 0.32) |
 | Nearby light bluer than any star | Gaia Catalogue of Nearby Stars: sources within 100 pc with BP-RP below -0.65, bluer than the Rayleigh-Jeans tail of any photosphere gives in Gaia's bands, each read for its BP/RP flux excess, image-peak and RUWE flags, white-dwarf probability and SIMBAD type | 1, a G 17 source with a BP/RP flux excess of 26 over the single-star relation, a blend or a corrupted BP and no bluer star |
 | Comets pushed toward the Sun | JPL Small-Body Database: 32 comets and 3 asteroids whose radial non-gravitational term A1 is negative, those negative at three sigma, each read for its arc and condition code, its transverse term A2 and the asymmetric outgassing parameter DT | 27 at three sigma: 1 on a short arc, 6 with a larger transverse term or a shifted outgassing law; 20 comets keep a radial term toward the Sun of 2e-10 to 6e-7 au per day squared on arcs of 3-75 years, which the symmetric outgassing law returns for a nucleus venting mostly after perihelion, and they await refits with an asymmetric law; the asteroid 4179 Toutatis's -3.1e-13 at 3.4 sigma is the lowest of the 22 asteroids SBDB fits a radial term for, whose median lies at +3.1 sigma, outward as sunlight pushes |
-| Galaxies warm without star formation | HECATE: 669 galaxies within 150 Mpc with GSWLC stellar masses and ultraviolet-optical star-formation rates and with WISE 22 um and IRAS far-infrared luminosities, the 38 quiescent ones (specific rate below 1e-11 per year) radiating more than a tenth of a solar luminosity per solar mass at 22 um within a factor of three of their far-infrared, against the quiescent population's own distribution; each read for HECATE's nuclear-activity flags and its inclination | 6, every one with an active nucleus; the quiescent galaxies' 22 um heat per mass runs 0.035 solar at the median and 0.55 at the 99th percentile |
+| Galaxies warm without star formation | HECATE: 669 galaxies within 150 Mpc with GSWLC stellar masses and ultraviolet-optical star-formation rates and with WISE 22 um and IRAS far-infrared luminosities, the 38 quiescent ones (specific rate below 1e-11 per year) radiating more than a tenth of a solar luminosity per solar mass at 22 um within a factor of three of their far-infrared, against the quiescent population's own distribution; each read for HECATE's two active-nucleus flags, its inclination and its homogenised star-formation rate, which counts the infrared | 6: 3 active nuclei, 1 edge-on disc, and CGCG 105-099 and UGC 12633, quiescent by their ultraviolet-optical rates and forming 10-15 solar masses a year once the infrared is counted; the quiescent galaxies' 22 um heat per mass runs 0.035 solar at the median and 0.55 at the 99th percentile (`hecate_warm_quiet`) |
 | Stars that dimmed since Hipparcos | Hipparcos new reduction and Gaia DR3: 17285 Hipparcos stars within 100 pc with a Gaia best neighbour, G against Hp through the sample's running median of G - Hp in BP-RP, those fainter by 0.5 mag and those brighter by as much, each read for Hipparcos's scatter and variability, Gaia's variability flag and RUWE, SIMBAD's type, a pair's name or component letter and a Gaia neighbour within 10 arcsec | 297 fainter against 4 brighter: 274 are pairs whose joint light Hipparcos measured and Gaia resolves, most by 0.75 mag, 23 are variables and 4 rest on Hipparcos records without a scatter; none unexplained |
 | Asteroids on polar orbits | JPL Small-Body Database: 51 bound orbits of asteroids and comets within 3 deg of 90 deg inclination, those within 0.5 deg of polar against a flat inclination density and the bands either side, each read for its arc and condition code | 9 against 8.5 expected (Poisson p 0.48): 8 long-period and Halley-type comets and one trans-Neptunian object |
 | Stars enriched evenly in every heavy element | GALAH DR4: 46868 stars with clean measurements of all eleven neutron-capture elements Sr, Y, Zr, Mo, Ru, Ba, La, Ce, Nd, Sm and Eu, those with every [X/Fe] above +0.4 dex and a spread below 0.1 dex across the eleven, against the stars depleted as evenly and the sample's own spread | none of the 41 stars enriched above +0.4 dex carries a flat pattern, and none is depleted evenly; the sample's spread across the elements runs 0.12 dex at the median against errors of 0.09 |
@@ -413,6 +413,9 @@ is its docstring, and its calls are in the recipes of `README.md`.
 | Deep temporary fades | Gaia Catalogue of Nearby Stars objects within 100 pc at G 13.5-19 north of -28 deg joined on IRSA to ZTF DR24 objects: 123,529 stars with 50 or more good points in both g and r | the faintest good point a magnitude below the median in both bands with the brightness skewed faint: 93; in their light curves 60 are faint on three nights or fewer, 3 are eclipsing binaries, 9 are fast movers (182-492 mas/yr) whose ZTF object steps down by 1.0-2.6 mag in both bands once the star has moved 1.1-1.8 arcsec from it, the newer object at the same star holding its brightness, and 21 hold faint nights in one band alone; no fade in both bands lasting more than a night (`ztf_deep_fades`, not recorded) |
 | Under our own radar fence at 217 MHz | the Gaia Catalogue of Nearby Stars within 10 pc at declinations 0 to +33 deg, which NAVSPASUR's 216.98 MHz fence has swept since 1961, carried to 2013.8; GLEAM EGC (VIII/100) within 2 arcmin and in the surrounding degree, the 216-223 MHz sub-band against the power law of the other nineteen | 3 GLEAM sources within 2 arcmin of the 86 stars, all 66-105 arcsec off and two of them one source beside the pair at 5.1 pc, against 2.0 from the fields' density (p 0.31), none above its power law at 216-223 MHz, where none of the 1596 field sources stands 3 sigma above either; limits 0.05-0.51 Jy at 3 sigma (median 0.12), 2e15 W isotropic over the 7.68 MHz sub-band at 5 pc for 0.1 Jy (`gleam_radar_fence`) |
 | Technetium in the nearest polluted white dwarf | van Maanen 2 (WD 0046+051, 4.3 pc): the ESO X-shooter UVB spectrum of 2016 (S/N 224), flat and free of metal lines across 4205-4312 A; NIST ASD lines and gf | Gaussians of the star's Fe I width (sigma 0.89 A) fixed at the Tc I resonance lines 4238.19, 4262.27 and 4297.06 A in the star's frame (-15.4 km/s from 15 Fe I lines): EW +2.0 +- 4.2, -1.3 +- 4.2 and -4.5 +- 4.5 mA against random places at robust sd 8.0 mA, 3-sigma 22.9 mA per line, the strength of a ground-term Fe I line of log gf -1.57 on the star's own curve of growth: N(Tc)/N(Fe) below about 0.03, partition functions and ionisation taken as equal; FEROS and UVES carry order-merging ripples there (`vma2_technetium`, not recorded) |
+| X-ray sources that went out | 2RXS (ROSAT All-Sky Survey, 1990-1991): 13,869 point sources at existence likelihood 10 or more and 0.05 ct/s or brighter, ten times eRASS1's limit, at Galactic longitude 181-359 deg, against the eRASS1 main catalogue (2019-2020) within 120 arcsec; SIMBAD within 90 arcsec, Gaia DR3 within 20 arcsec, AllWISE colours | 532 gone, all explained: 254 pieces of extended emission (galaxy clusters, SN 1006, three or more missing detections within half a degree), 161 variable nuclei, X-ray binaries and cataclysmic variables, 64 stars caught in a flare, 37 beside sources of 10 ct/s such as Sco X-1, 1 ROSAT variable, and 15 at likelihood 10-17, where 22 per cent of the sample lies and all 15 do (p 1.5e-10), the catalogue's spurious tail (`rosat_erass1_gone`) |
+| Heat without cosmic rays | HECATE: 820 galaxies within 200 Mpc north of -40 deg with IRAS 60 and 100 um fluxes bright enough for NVSS to reach q 3.3; NVSS components summed within each galaxy's D25 semi-major axis (45 arcsec to 3 arcmin); q = log of far-infrared over 1.4 GHz flux, 2.34 +- 0.26 in nature (Yun et al. 2001) | q median 2.34, scatter 0.17 over the 814 detected; 11 above 3.3, none unexplained: with NED's 1.4 GHz flux integrated over the galaxy's map M31, M33, M101 and IC 356 return to q 2.2-2.7 (IC 356's 30 mJy from Condon, Cotton and Broderick 2002, where NVSS's component list holds nothing, its cold dust at a 100/60 um ratio of 6.4 putting it where the coldest galaxies run, median q 2.65), 4 more discs are wider than the radio sum covers (NGC 300, NGC 247, NGC 6822, NGC 4753), 1 lies in the Galactic plane's cirrus, NGC 1317 has an active nucleus and NGC 1377 is the radio-deficient nascent starburst of Roussel et al. (2003) (`fir_radio_quiet`) |
+| Quiet but spotted | McQuillan et al. (2014) Kepler rotation periods and spot amplitudes: 104 dwarfs below 4000 K spinning faster than 10 d with amplitudes above their peers' median and a Gaia DR3 match, against the Kepler flare stars of Yang and Liu (2019) | 103 flare; the one quiet star is a binary at RUWE above 1.4, against 1.0 quiet star expected from the brightness bins' quiet shares; none unexplained (`kepler_quiet_spotted`) |
 
 ## Queued
 

README.md

@@ -428,9 +428,11 @@ on any vet that is not a code.
 | Planets in step with their star's spin | planets of the composite table whose host carries a rotation period; an orbital period at a whole multiple of 2-20 of the rotation within 0.5 per cent, against the half-integer windows and uniform fractions, each with the source of the rotation period | `seti.py transits spinsteps [record] [--tol-rel 0.005] [--n-max 20]` |
 | Nearby light bluer than any star | Gaia Catalogue of Nearby Stars objects with BP-RP below -0.65; the BP/RP flux excess, the image-peak and RUWE flags, the white-dwarf probability and SIMBAD's type as vets | `seti.py ultraviolet bluest [record] [--colour-max -0.65] [--excess-max 0.2]` |
 | Comets pushed toward the Sun | comets and asteroids of the Small-Body Database with a negative radial non-gravitational term A1 at three sigma; the arc, condition code, transverse term and asymmetric outgassing parameter as vets | `seti.py solar_system sunward [record] [--sigma 3]` |
-| Galaxies warm without star formation | quiescent HECATE galaxies within 150 Mpc with GSWLC masses and rates; a 22 um luminosity above a tenth of a solar luminosity per solar mass and within a factor of three of the far-infrared, against the quiescent population, with the nuclear-activity flags and the inclination as vets | `seti.py infrared warmgalaxies [record] [--d-max-mpc 150] [--ssfr-max -11] [--heat-min -1.0] [--warm-min -0.5] [--edge-on-deg 70]` |
+| Galaxies warm without star formation | quiescent HECATE galaxies within 150 Mpc with GSWLC masses and rates; a 22 um luminosity above a tenth of a solar luminosity per solar mass and within a factor of three of the far-infrared, against the quiescent population, with the active-nucleus flags, the inclination and the infrared-inclusive star-formation rate as vets | `seti.py infrared warmgalaxies [record] [--d-max-mpc 150] [--ssfr-max -11] [--heat-min -1.0] [--warm-min -0.5] [--edge-on-deg 70]` |
+| Heat without cosmic rays | HECATE galaxies with IRAS fluxes against NVSS summed over the disc; the far-infrared to radio ratio q above 3.3, with NED's flux integrated over the map, the Galactic plane, discs wider than the sum, active nuclei, SIMBAD types and published explanations as vets | `seti.py infrared radioquiet [record] [--d-max-mpc 200] [--q-flag 3.3] [--nvss-limit-mjy 2.5] [--b-min 10] [--extent-arcmin 6]` |
 | Stars that dimmed since Hipparcos | Hipparcos stars within 100 pc with a Gaia DR3 best neighbour; G - Hp against the running median in BP-RP, a step of 0.5 mag either way, with Hipparcos's scatter and variability, Gaia's variability flag and RUWE, SIMBAD's type, pair names and Gaia neighbours as vets | `seti.py dimming hipdim [record] [--dmag 0.5] [--plx-min 10] [--bright-g 3] [--pair-arcsec 10]` |
 | Single-transit flashes in Hipparcos | Hipparcos epoch photometry of quiet stars within 100 pc; one transit 0.3 mag and five robust deviations brighter than the median against the transits as much fainter, with catalogue-wide events, recurrence, variability and RUWE as vets and the pass partner as the test of a real flash | `seti.py flashes hipflash [record] [--mag 0.3] [--sig 5] [--quiet 0.02] [--pair-days 0.02] [--plx-min 10]` |
+| Quiet but spotted | Kepler M dwarfs spinning faster than 10 days with large spot amplitudes against the Kepler flare list, in brightness bins where most peers flare; Gaia neighbours, RUWE and colour as vets | `seti.py flashes quietspots [record] [--teff-max 4000] [--prot-max 10] [--frac-min 0.7] [--g-step 1] [--beam-arcsec 4]` |
 | Asteroids on polar orbits | asteroids and comets of the Small-Body Database on bound orbits within 3 deg of 90 deg inclination; those within 0.5 deg of polar against a flat density and the bands either side, each read for its arc and condition code | `seti.py solar_system polar [record] [--half-width-deg 0.5] [--band-deg 3] [--condition-max 6]` |
 | Stars enriched evenly in every heavy element | GALAH DR4 stars with all eleven neutron-capture elements clean; every [X/Fe] above +0.4 dex with a spread below 0.1 dex, against the stars depleted as evenly and the sample's own spread | `seti.py composition evenheavy [record] [--mean-min 0.4] [--spread-max 0.1]` |
 | Transits longer than the star allows | transiting planets of the composite table; a duration beyond 1.25 times the central circular maximum at the catalogued geometry, with the duration error, the sources of duration and radius and the eccentricity as vets | `seti.py transits longtransits [record] [--excess 1.25] [--err-max 0.2]` |
@@ -504,6 +506,7 @@ on any vet that is not a code.
 | Retrograde co-orbitals of the giant planets | SBDB bodies sharing a giant planet's semi-major axis on retrograde orbits of modest eccentricity, read against the published cases and the arcs | `seti.py solar_system retroco [record] [--a-tol 0.05] [--i-min 90] [--e-max 0.5] [--arc-min-d 365]` |
 | Radio sources that switched off | Compact NVSS sources of 1995 absent from the VLASS quick-look components of 2019, read against SIMBAD, bright neighbours, RACS-low and the catalogue's holes, then the later surveys, the lobes of brighter sources and the VLASS images themselves | `seti.py radio_sky switchoff [record] [--flux-min-mjy 100] [--major-max-arcsec 25] [--radius-arcsec 20] [--dec-min -40] [--hole-deg 0.25] [--b-min 10] [--bright-jy 5] [--bright-deg 0.5] [--racs-dec-max 41] [--later-arcsec 30] [--lobe-arcsec 90] [--image-sigma 10]` |
 | Under our own radar fence at 217 MHz | Gaia Catalogue of Nearby Stars within 10 pc in the strip NAVSPASUR's 216.98 MHz fence sweeps; GLEAM EGC within 2 arcmin against the density of the surrounding degree, the 216-223 MHz sub-band against the power law of the other nineteen, and three times the sub-band's local noise as the limit where nothing is found | `seti.py radio_sky radarfence [record] [--dist-pc 10] [--dec-min 0] [--dec-max 33] [--near-arcmin 2] [--field-deg 1] [--sigma 3] [--epoch 2013.8]` |
+| X-ray sources that went out | 2RXS point sources ten times eRASS1's limit in eRASS1's half of the sky against the eRASS1 main catalogue; SIMBAD types, Gaia stars, WISE AGN colours, bright neighbours, clumps of missing detections and the likelihood of what is left as vets | `seti.py highenergy rosatgone [record] [--exi-min 10] [--crate-min 0.05] [--match-arcsec 120] [--edge-deg 1] [--simbad-arcsec 90] [--gaia-arcsec 20] [--bright-ct 10] [--bright-deg 0.5] [--clump-deg 0.5] [--spur-exi 20]` |
 | Bodies pushed out of their orbital planes | SBDB bodies whose normal non-gravitational term rivals the radial, read against jets, dark comets and arcs | `seti.py solar_system outofplane [record] [--sigma 3] [--ratio-min 0.5] [--arc-min-d 1000]` |
 | Gamma rays harder than any shock | 4FGL sources whose photon index is harder than any accelerator gives, read against their associations and extension | `seti.py highenergy hardgamma [record] [--index-max 1.3] [--sig-min 10]` |
 | Stars whose distance changed since Hipparcos | Hipparcos against Gaia parallaxes with scaled errors, outliers read for misidentification, binarity, solution quality and the Gaussian tail | `seti.py neighbourhood receded [record] [--sigma 5] [--min-mas 2] [--g-sat 6] [--f2-max 3] [--match-mag 2]` |

RESULTS.md

@@ -4,7 +4,7 @@ Generated by `python -I seti.py store table` from the Searched table of `LEDGER.
 
 A run is one row of the store's `searches` table. Its sample is the number of things it searched (`n_targets`) in the unit named (`n_unit`), and samples in different units are never added. Its rows are its rows in the `coverage` table, one for each target for which it recorded a statistic, a limit or a vet; many runs record only their flagged targets, so a run's rows are not its sample. A vetted row carries the codes of `README.md` for the explanations that removed it.
 
-518 runs of 406 approaches; 514356 per-target rows over 432573 distinct target identifiers.
+521 runs of 409 approaches; 514900 per-target rows over 433116 distinct target identifiers.
 
 ## Curved drift tracks, radio
 
@@ -2112,11 +2112,11 @@ A run is one row of the store's `searches` table. Its sample is the number of th
 
 ## Galaxies warm without star formation
 
-**Data.** HECATE: 669 galaxies within 150 Mpc with GSWLC stellar masses and ultraviolet-optical star-formation rates and with WISE 22 um and IRAS far-infrared luminosities, the 38 quiescent ones (specific rate below 1e-11 per year) radiating more than a tenth of a solar luminosity per solar mass at 22 um within a factor of three of their far-infrared, against the quiescent population's own distribution; each read for HECATE's nuclear-activity flags and its inclination
+**Data.** HECATE: 669 galaxies within 150 Mpc with GSWLC stellar masses and ultraviolet-optical star-formation rates and with WISE 22 um and IRAS far-infrared luminosities, the 38 quiescent ones (specific rate below 1e-11 per year) radiating more than a tenth of a solar luminosity per solar mass at 22 um within a factor of three of their far-infrared, against the quiescent population's own distribution; each read for HECATE's two active-nucleus flags, its inclination and its homogenised star-formation rate, which counts the infrared
 
 **Runs.** `hecate_warm_quiet`: 38 galaxies, 6 per-target rows, 6 vetted
 
-**Result.** 6, every one with an active nucleus; the quiescent galaxies' 22 um heat per mass runs 0.035 solar at the median and 0.55 at the 99th percentile
+**Result.** 6: 3 active nuclei, 1 edge-on disc, and CGCG 105-099 and UGC 12633, quiescent by their ultraviolet-optical rates and forming 10-15 solar masses a year once the infrared is counted; the quiescent galaxies' 22 um heat per mass runs 0.035 solar at the median and 0.55 at the 99th percentile (`hecate_warm_quiet`)
 
 ## Stars that dimmed since Hipparcos
 
@@ -3253,3 +3253,27 @@ A run is one row of the store's `searches` table. Its sample is the number of th
 **Runs.** none in the store
 
 **Result.** Gaussians of the star's Fe I width (sigma 0.89 A) fixed at the Tc I resonance lines 4238.19, 4262.27 and 4297.06 A in the star's frame (-15.4 km/s from 15 Fe I lines): EW +2.0 +- 4.2, -1.3 +- 4.2 and -4.5 +- 4.5 mA against random places at robust sd 8.0 mA, 3-sigma 22.9 mA per line, the strength of a ground-term Fe I line of log gf -1.57 on the star's own curve of growth: N(Tc)/N(Fe) below about 0.03, partition functions and ionisation taken as equal; FEROS and UVES carry order-merging ripples there (`vma2_technetium`, not recorded)
+
+## X-ray sources that went out
+
+**Data.** 2RXS (ROSAT All-Sky Survey, 1990-1991): 13,869 point sources at existence likelihood 10 or more and 0.05 ct/s or brighter, ten times eRASS1's limit, at Galactic longitude 181-359 deg, against the eRASS1 main catalogue (2019-2020) within 120 arcsec; SIMBAD within 90 arcsec, Gaia DR3 within 20 arcsec, AllWISE colours
+
+**Runs.** `rosat_erass1_gone`: 13,869 sources, 532 per-target rows, 532 vetted
+
+**Result.** 532 gone, all explained: 254 pieces of extended emission (galaxy clusters, SN 1006, three or more missing detections within half a degree), 161 variable nuclei, X-ray binaries and cataclysmic variables, 64 stars caught in a flare, 37 beside sources of 10 ct/s such as Sco X-1, 1 ROSAT variable, and 15 at likelihood 10-17, where 22 per cent of the sample lies and all 15 do (p 1.5e-10), the catalogue's spurious tail (`rosat_erass1_gone`)
+
+## Heat without cosmic rays
+
+**Data.** HECATE: 820 galaxies within 200 Mpc north of -40 deg with IRAS 60 and 100 um fluxes bright enough for NVSS to reach q 3.3; NVSS components summed within each galaxy's D25 semi-major axis (45 arcsec to 3 arcmin); q = log of far-infrared over 1.4 GHz flux, 2.34 +- 0.26 in nature (Yun et al. 2001)
+
+**Runs.** `fir_radio_quiet`: 820 galaxies, 11 per-target rows, 11 vetted
+
+**Result.** q median 2.34, scatter 0.17 over the 814 detected; 11 above 3.3, none unexplained: with NED's 1.4 GHz flux integrated over the galaxy's map M31, M33, M101 and IC 356 return to q 2.2-2.7 (IC 356's 30 mJy from Condon, Cotton and Broderick 2002, where NVSS's component list holds nothing, its cold dust at a 100/60 um ratio of 6.4 putting it where the coldest galaxies run, median q 2.65), 4 more discs are wider than the radio sum covers (NGC 300, NGC 247, NGC 6822, NGC 4753), 1 lies in the Galactic plane's cirrus, NGC 1317 has an active nucleus and NGC 1377 is the radio-deficient nascent starburst of Roussel et al. (2003) (`fir_radio_quiet`)
+
+## Quiet but spotted
+
+**Data.** McQuillan et al. (2014) Kepler rotation periods and spot amplitudes: 104 dwarfs below 4000 K spinning faster than 10 d with amplitudes above their peers' median and a Gaia DR3 match, against the Kepler flare stars of Yang and Liu (2019)
+
+**Runs.** `kepler_quiet_spotted`: 104 stars, 1 per-target row, 1 vetted
+
+**Result.** 103 flare; the one quiet star is a binary at RUWE above 1.4, against 1.0 quiet star expected from the brightness bins' quiet shares; none unexplained (`kepler_quiet_spotted`)

coverage.db

@@ -1,3 +1,3 @@
 version https://git-lfs.github.com/spec/v1
-oid sha256:3e4ac89c6fdf68b1deeeb0cd0aac966058a55c3c924db87bf2546e83357f3e73
-size 110522368
+oid sha256:36338eb3e6cdf3a9c0b100540ba5357a9e4ba5d1a21a5ec90bc38e853ec716bc
+size 110686208

setitools/composition.py

@@ -3387,6 +3387,12 @@ def giants_residuals(spectra, vsys, whiten=True, band="optical"):
     master = np.nanmedian(F[~ins], axis=0)
     with np.errstate(all="ignore"):
         R, ER = F / master, E / master
+    if band == "nir":
+        # in the near infrared the median spectrum falls toward zero in the water bands, where the ratio would swamp each row's cubic, so the
+        # telluric bands and the samples where the median lies at 0.2 or below are set aside before the fit
+        R[:, ~(master > 0.2)] = np.nan
+        for lo, hi in giants_BANDS[band]["telluric"]:
+            R[:, (g > lo) & (g < hi)] = np.nan
     x = (np.arange(len(g)) - len(g) / 2.0) / len(g)
     for i in range(len(R)):
         good = np.isfinite(R[i])
@@ -3427,7 +3433,8 @@ def giants_ccf(R, T, sd=None, vmax=300):
     dev = "cuda" if torch.cuda.is_available() else "cpu"
     n = R.shape[1]
     w = 1.0 / (sd if sd is not None else np.ones(n))
-    Rt = torch.as_tensor(np.nan_to_num(R * w[None, :], nan=0.0), dtype=torch.float32, device=dev)
+    # a non-finite sample, infinite as well as missing, counts as zero, since an infinite one would meet its zero weight as a NaN
+    Rt = torch.as_tensor(np.where(np.isfinite(R), R * w[None, :], 0.0), dtype=torch.float32, device=dev)
     Wt = torch.as_tensor(np.where(np.isfinite(R).all(0), w, 0.0), dtype=torch.float32, device=dev)
     Tt = torch.as_tensor(T, dtype=torch.float32, device=dev)
     out = torch.zeros((R.shape[0], 2 * vmax + 1), device=dev)
@@ -3446,7 +3453,7 @@ def giants_ccf_fft(R, T, sd=None, vmax=300):
     n = R.shape[1]
     w = 1.0 / (sd if sd is not None else np.ones(n))
     Wt = np.where(np.isfinite(R).all(0), w, 0.0)
-    X = np.nan_to_num(R * w[None, :], nan=0.0) * Wt[None, :]
+    X = np.where(np.isfinite(R), R * w[None, :], 0.0) * Wt[None, :]
     Tt = np.asarray(T, float)
     m = 1 << int(np.ceil(np.log2(2 * n)))
     ft = np.conj(np.fft.rfft(Tt, m))
@@ -3813,6 +3820,12 @@ def companions_residuals(spectra, vsys, band, align=True, info=None):
             info["vstar"] = np.array([(meta["rv"] if meta.get("rv") is not None else vsys) + (dv if mv else 0.0) for (_, _, _, meta), dv, mv in zip(spectra, shifts, moved)])
     with np.errstate(all="ignore"):
         R = F / master
+    if band == "nir":
+        # in the near infrared the median spectrum falls toward zero in the water bands, where the ratio would swamp each row's cubic, so the
+        # telluric bands and the samples where the median lies at 0.2 or below are set aside before the fit
+        R[:, ~(master > 0.2)] = np.nan
+        for lo, hi in giants_BANDS[band]["telluric"]:
+            R[:, (g > lo) & (g < hi)] = np.nan
     x = (np.arange(len(g)) - len(g) / 2.0) / len(g)
     for i in range(len(R)):
         good = np.isfinite(R[i])
@@ -4253,7 +4266,7 @@ binaries_DOC = r"""Companions of single-lined binaries along their known orbits:
     python -I seti.py composition binaries template NAME [--inst HARPS] [--out FILE.npy]       a primary's line depths from its own spectra at its catalogue orbit's velocities
     python -I seti.py composition binaries template --spectra DIR --object NAME [NAME ...] --out FILE.npy   a star's line depths from its spectra in DIR in its own frame
     python -I seti.py composition binaries build NAME [NAME ...] [--inst HARPS] [--min-snr 10] [--per-star 300] [--keep] [--template FILE.npy] [--variant V] [--band LO HI]
-                                                 [--inject RATIO --inject-k2 KMS [--inject-shift 0.25] [--inject-template FILE.npy]]   each binary's rows: ROOT/rows or ROOT/rows_V
+                                                 [--inject RATIO --inject-k2 KMS [--inject-shift 0.25] [--inject-template FILE.npy]] [--orbit-frames]   each binary's rows: ROOT/rows or ROOT/rows_V
     python -I seti.py composition binaries search [NAME ...] [--inst HARPS] [--dv 1.0]          the summit over K2 at a fixed systemic offset: ROOT/results
     python -I seti.py composition binaries calibrate [NAME ...] [--inst HARPS] [--draws 10000] [--window 5]   the summits' false-alarm probabilities: ROOT/calibrate
     python -I seti.py composition binaries calibrate --aggregate [--inst HARPS]                the survey's binaries taken together: ROOT/calibrate/survey_<inst>.json
@@ -4479,8 +4492,8 @@ def binaries_planck(w, t2, t1):
 
 
 def binaries_build(names, inst="HARPS", min_snr=10.0, per_star=300, streams=4, rv_template=None, template=None, keep=False, variant=None, band=None,
-                   inject=0.0, inject_k2=None, inject_shift=0.25, inject_dv=1.0, inject_template=None, teff1=6000.0, teff2=3400.0):
-    """Each named binary's rows: its spectra's velocities measured against RV_TEMPLATE and the orbit fitted to them, each spectrum put in the primary's frame at its own velocity, the residuals cross-correlated with the M-dwarf TEMPLATE (zero outside BAND where given) and the star's pattern taken out, kept as ROOT/rows/<tag>.npz, or ROOT/rows_<VARIANT>, with a note; INJECT lays the lines of INJECT_TEMPLATE (the M-dwarf template unless given) into every spectrum at that V-band flux ratio, scaled as a TEFF2 black body against a TEFF1 one, along the orbit moved by INJECT_SHIFT of a period at INJECT_K2 and offset INJECT_DV."""
+                   inject=0.0, inject_k2=None, inject_shift=0.25, inject_dv=1.0, inject_template=None, teff1=6000.0, teff2=3400.0, orbit_frames=False):
+    """Each named binary's rows: its spectra's velocities measured against RV_TEMPLATE and the orbit fitted to them, each spectrum put in the primary's frame at its own velocity, the residuals cross-correlated with the M-dwarf TEMPLATE (zero outside BAND where given) and the star's pattern taken out, kept as ROOT/rows/<tag>.npz, or ROOT/rows_<VARIANT>, with a note; with ORBIT_FRAMES each spectrum is put at the plan orbit's velocity instead and left unaligned, as in the near infrared, where a solar-type primary's own lines are too few to place it; INJECT lays the lines of INJECT_TEMPLATE (the M-dwarf template unless given) into every spectrum at that V-band flux ratio, scaled as a TEFF2 black body against a TEFF1 one, along the orbit moved by INJECT_SHIFT of a period at INJECT_K2 and offset INJECT_DV."""
     rv_template = rv_template or os.path.join(binaries_ROOT, "templates", "rv.npy")
     template = template or os.path.join(binaries_ROOT, "templates", "mdwarf.npy")
     Trv, T = np.load(rv_template), np.load(template)
@@ -4497,8 +4510,12 @@ def binaries_build(names, inst="HARPS", min_snr=10.0, per_star=300, streams=4, r
             print("%s: %d spectra loaded of %d listed; too few" % (target["name"], len(spectra), listed), flush=True)
             continue
         t = np.array([s[3]["mjd"] for s in spectra])
-        v = binaries_velocities(spectra, Trv, band=giants_band(inst))
-        P, tp, e, w, k1, gamma, rms = binaries_refit(t, v, *binaries_orbit(target)[:5])
+        if orbit_frames:
+            P, tp, e, w, k1, gamma = binaries_orbit(target)
+            v, rms = gamma + k1 * binaries_radial(t, P, tp, e, w), 0.0
+        else:
+            v = binaries_velocities(spectra, Trv, band=giants_band(inst))
+            P, tp, e, w, k1, gamma, rms = binaries_refit(t, v, *binaries_orbit(target)[:5])
         for s, x in zip(spectra, v):
             s[3]["rv"] = float(x)
         if inject:
@@ -4508,7 +4525,9 @@ def binaries_build(names, inst="HARPS", min_snr=10.0, per_star=300, streams=4, r
             spectra = [(wl, fl * (1.0 - inject * scale(wl) * np.interp(wl, g * (1.0 + x / 299792.458), Ti, left=0.0, right=0.0)), er, meta)
                        for (wl, fl, er, meta), x in zip(spectra, v2)]
         info = {}
-        R, sd, _ = companions_residuals(spectra, gamma, giants_band(inst), info=info)
+        R, sd, _ = companions_residuals(spectra, gamma, giants_band(inst), align=not orbit_frames, info=info)
+        if "vstar" not in info:
+            info["vstar"] = np.asarray(v, float)
         ccf = companions_unpattern(giants_ccf(R, Tc, sd))
         tag = binaries_tag(target)
         np.savez_compressed(os.path.join(folder, tag + ".npz"), ccf=ccf.astype(np.float32), mjd=t, berv=np.array([s[3]["berv"] for s in spectra], float),
@@ -4920,6 +4939,7 @@ def binaries_cli():
     p.add_argument("--dv-fixed", type=float, default=1.0)
     p.add_argument("--components", type=int, default=None)
     p.add_argument("--shift", type=float, default=0.0, help="measure: the track moved by this fraction of a period")
+    p.add_argument("--orbit-frames", action="store_true", help="build: each spectrum at the plan orbit's velocity, unaligned, as for the near infrared")
     a = p.parse_args()
     if a.mode == "plan":
         binaries_plan(a.min)
@@ -4927,7 +4947,7 @@ def binaries_cli():
         binaries_template(a.names[0] if a.names else None, a.inst, a.out, a.spectra, a.object)
     elif a.mode == "build":
         binaries_build(a.names, a.inst, a.min_snr, a.per_star, a.streams, a.rv_template, a.template, a.keep, a.variant, a.band, a.inject, a.inject_k2, a.inject_shift,
-                       a.inject_dv, a.inject_template)
+                       a.inject_dv, a.inject_template, orbit_frames=a.orbit_frames)
     elif a.mode == "search":
         binaries_search(a.names, a.inst, a.dv, a.draws or 300, a.detect)
     elif a.mode == "calibrate":

setitools/flashes.py

@@ -3807,7 +3807,8 @@ def sndeficit_main(record=False, d_max_mpc=20.0, sn_per_msun=0.01, year_min=1980
         tir, uv, ltir = core.num(g["logSFRTIR"]), core.num(g["logSFRGSW"]), core.num(g["logLTIR"])
         # a luminous infrared galaxy hides its supernovae in dust; an active nucleus lends its infrared to the rate
         dusty = (tir is not None and uv is not None and tir - uv > math.log10(dust_ratio)) or (ltir is not None and ltir > ltir_max)
-        agn = bool((g["AGNs17"] or "").strip() or (g["AGNh"] or "").strip() or (g["Nuclear"] or "").strip() not in ("", "0"))
+        # HECATE writes its two active-nucleus flags as Y, N or ?, and Nuclear as -1 where no emission-line class exists
+        agn = "Y" in (str(g["AGNs17"]).strip(), str(g["AGNh"]).strip())
         out.append(dict(pgc=g["PGC"], name=(g["NED"] or "").strip(), d=float(g["D"]), sfr=round(sfr, 2), expected=round(exp, 1), seen=n, p=p, dusty=dusty, agn=agn,
                         vet="dust" if dusty else "background-galaxy" if agn else None))
     short = sorted([o for o in out if o["p"] < p_max], key=lambda o: o["p"])
@@ -4907,6 +4908,74 @@ def chirps_main(record=False, rate_max=60.0, rate_min=1.0, p_max=1e-12, tmag_max
         print("recorded")
 
 
+quietspots_DOC = r"""Quiet but spotted: a fast-spinning, heavily spotted M dwarf runs a saturated dynamo, and its peers flare, so one that never flares while its spots show the dynamo at work is a star whose flaring something suppresses; McQuillan et al. (2014) Kepler rotation periods and spot amplitudes of the dwarfs below 4000 K spinning faster than 10 days with an amplitude above the median of their peers, against the Kepler stars Yang and Liu (2019) found flaring, matched to Gaia DR3 for brightness, colour, distance and RUWE; a quiet star is one in a brightness bin where most of its peers flare, the expected count of quiet stars is the sum of each bin's quiet share, and each quiet star is read for a Gaia neighbour within 4 arcsec and 3 mag (spots or light from another star), RUWE above 1.4 (a binary whose companion carries the spots), and a Gaia colour or absolute magnitude that is not an M dwarf's.
+
+    python -I seti.py flashes quietspots [record] [--teff-max 4000] [--prot-max 10] [--frac-min 0.7] [--g-step 1] [--beam-arcsec 4]"""
+
+
+def quietspots_main(record=False, teff_max=4000.0, prot_max=10.0, frac_min=0.7, g_step=1.0, beam_arcsec=4.0):
+    rot = core.tap("vizier", 'SELECT KIC, Teff, "log(g)", Prot, Rper, _RA, _DE FROM "J/ApJS/211/24/table1" WHERE Teff < %g AND Prot < %g' % (teff_max, prot_max), timeout=600)
+    flaring = {str(r["KIC"]).strip(): num(r["Nf"]) for r in core.tap("vizier", 'SELECT KIC, Nf FROM "J/ApJS/241/29/table1"', timeout=600)}
+    st = [dict(source_id=str(r["KIC"]).strip(), ra=float(r["_RA"]), dec=float(r["_DE"]), teff=num(r["Teff"]), logg=num(r["log(g)"]), prot=num(r["Prot"]), amp=num(r["Rper"]))
+          for r in rot if num(r["log(g)"]) and num(r["log(g)"]) > 4.0 and num(r["Rper"])]
+    med = float(np.median([s["amp"] for s in st]))
+    st = [s for s in st if s["amp"] >= med]
+    gaia = {}
+    for m in core.xmatch(st, "I/355/gaiadr3", 2.0):
+        if m["source_id"] not in gaia or float(m["angDist"]) < float(gaia[m["source_id"]]["angDist"]):
+            gaia[m["source_id"]] = m
+    st = [s for s in st if s["source_id"] in gaia and num(gaia[s["source_id"]].get("Gmag")) is not None]
+    for s in st:
+        g = gaia[s["source_id"]]
+        s["g"], s["bprp"], s["plx"], s["ruwe"] = num(g["Gmag"]), num(g.get("BP-RP")), num(g.get("Plx")), num(g.get("RUWE"))
+        s["nf"] = flaring.get(s["source_id"])
+        s["bin"] = int(math.floor(s["g"] / g_step))
+    bins = {}
+    for s in st:
+        bins.setdefault(s["bin"], []).append(s)
+    share = {b: sum(x["nf"] is not None for x in v) / len(v) for b, v in bins.items()}
+    quiet = [s for s in st if s["nf"] is None and share[s["bin"]] >= frac_min]
+    expected = sum(1 - share[b] for b, v in bins.items() if share[b] >= frac_min for _ in v)
+    near = {}
+    for k in range(0, len(quiet), 50):
+        chunk = quiet[k:k + 50]
+        where = " OR ".join("1 = CONTAINS(POINT('ICRS', ra, dec), CIRCLE('ICRS', %.6f, %.6f, %.6f))" % (s["ra"], s["dec"], beam_arcsec / 3600) for s in chunk)
+        for r in core.tap("gaia", "SELECT source_id, ra, dec, phot_g_mean_mag FROM gaiadr3.gaia_source WHERE " + where):
+            for s in chunk:
+                if core.sep_arcsec(float(r["ra"]), float(r["dec"]), s["ra"], s["dec"]) < beam_arcsec and num(r["phot_g_mean_mag"]) is not None \
+                        and abs(num(r["phot_g_mean_mag"]) - s["g"]) > 0.01 and num(r["phot_g_mean_mag"]) < s["g"] + 3:
+                    near.setdefault(s["source_id"], []).append(round(num(r["phot_g_mean_mag"]), 1))
+    for s in quiet:
+        mg = s["g"] + 5 * math.log10(s["plx"] / 100) if s["plx"] and s["plx"] > 0 else None
+        s["mg"] = mg
+        s["near"] = near.get(s["source_id"], [])
+        s["vet"] = ("blend" if s["near"] else "binary" if s["ruwe"] and s["ruwe"] > 1.4 else
+                    "misidentified" if (s["bprp"] is not None and s["bprp"] < 1.5) or mg is None or mg < 5.0 else None)
+    left = [s for s in quiet if s["vet"] is None]
+    print("%d Kepler dwarfs below %g K spinning faster than %g d with spot amplitudes above their median (%.0f), %d with Gaia; flaring in Yang and Liu (2019): %d; "
+          "quiet in brightness bins where %.0f per cent or more flare: %d against %.1f expected from those bins' quiet shares; %d blends, %d binaries, %d not M dwarfs; %d left" % (
+              len(rot), teff_max, prot_max, med, len(st), sum(s["nf"] is not None for s in st), 100 * frac_min, len(quiet), expected,
+              sum(s["vet"] == "blend" for s in quiet), sum(s["vet"] == "binary" for s in quiet), sum(s["vet"] == "misidentified" for s in quiet), len(left)))
+    for b in sorted(bins):
+        print("  G %.0f-%.0f: %d stars, %.0f per cent flaring" % (b * g_step, (b + 1) * g_step, len(bins[b]), 100 * share[b]))
+    for s in left:
+        print("  KIC %s: G %.2f, BP-RP %s, M_G %s, Teff %.0f, P %.2f d, amplitude %.0f, RUWE %s" % (
+            s["source_id"], s["g"], "%.2f" % s["bprp"] if s["bprp"] is not None else "-", "%.2f" % s["mg"] if s["mg"] is not None else "-", s["teff"], s["prot"], s["amp"], s["ruwe"]))
+    if record:
+        p_left = 1.0 - sum(math.exp(-expected) * expected ** i / math.factorial(i) for i in range(len(quiet)))
+        core.record_run(
+            "kepler_quiet_spotted", "Quiet but spotted", "McQuillan et al. (2014) Kepler rotation (J/ApJS/211/24); Yang and Liu (2019) Kepler flares (J/ApJS/241/29); Gaia DR3 by CDS XMatch",
+            "%d Kepler dwarfs below %g K faster than %g d with spot amplitudes above their median" % (len(st), teff_max, prot_max), os.path.abspath(__file__), "quietspots",
+            len(st), "stars", dict(teff_max=teff_max, prot_max=prot_max, frac_min=frac_min, g_step=g_step, beam_arcsec=beam_arcsec),
+            dict(flaring=sum(s["nf"] is not None for s in st), quiet=len(quiet), expected=expected, p=p_left),
+            "a spotted fast rotator absent from the flare list in a brightness bin where %.0f per cent of its peers flare" % (100 * frac_min),
+            "%d quiet against %.1f expected (p %.2f): %d blends, %d binaries, %d not M dwarfs; %d left" % (len(quiet), expected, p_left, sum(s["vet"] == "blend" for s in quiet),
+                                                                                                   sum(s["vet"] == "binary" for s in quiet), sum(s["vet"] == "misidentified" for s in quiet), len(left)),
+            [dict(target="kic:%s" % s["source_id"], stat=s["prot"], unit="rotation period, d, no flare", vet=s["vet"] or ("chance-rate" if p_left > 0.01 else None),
+                  note="G %.2f, BP-RP %s, Teff %.0f, amplitude %.0f, RUWE %s, neighbours %s" % (s["g"], s["bprp"], s["teff"], s["amp"], s["ruwe"], s["near"] or "none")) for s in quiet])
+        print("recorded")
+
+
 hipflash_DOC = r"""Single-transit brightenings in Hipparcos epoch photometry: stars within 100 pc with a Gaia DR3 best neighbour and a robust scatter below 0.02 mag, a clean transit brighter than the star's median by 0.3 mag and five robust deviations, against the transits as much fainter, the mirror; each read for a catalogue-wide event (two or more other stars within half an hour, the satellite's own), a recurrence on the same star, Hipparcos's variability type, SIMBAD's type and RUWE, and for the pass partner 20.6 minutes away, which sees the same star through the other field of view and shares a real flash.
 
     python -I seti.py flashes hipflash [record] [--mag 0.3] [--sig 5] [--quiet 0.02] [--pair-days 0.02] [--plx-min 10]"""
@@ -4985,6 +5054,7 @@ def hipflash_main(record=False, mag=0.3, sig=5.0, quiet=0.02, pair_days=0.02, pl
 
 APPROACHES = {
     "hipflash": (lambda: core.run_cli(hipflash_main), hipflash_DOC),
+    "quietspots": (lambda: core.run_cli(quietspots_main), quietspots_DOC),
     "chirps": (lambda: core.run_cli(chirps_main), chirps_DOC),
     "tones": (lambda: core.run_cli(tones_main), tones_DOC),
     "flareecho": (lambda: core.run_cli(flareecho_main), flareecho_DOC),

setitools/highenergy.py

@@ -3363,7 +3363,109 @@ def gasunder_main(record=False, l_deg=312.07, b_deg=0.0, radius_deg=2.0):
         print("recorded")
 
 
+rosatgone_DOC = r"""X-ray sources that went out: the 2RXS catalogue of the ROSAT All-Sky Survey (1990-1991) in the half of the sky eROSITA-DE's first survey covers (Galactic longitude 180-360 deg), point sources at existence likelihood 10 or more and 0.05 ct/s or brighter, ten times eRASS1's limit, a degree clear of the half's edges, against the eRASS1 main catalogue (2019-2020) within 120 arcsec, the reach of a faint 2RXS position; eRASS1 covers its half whole, so a source with no counterpart faded by an order of magnitude in thirty years; each read for SIMBAD's type within 90 arcsec (flare stars caught in a flare, galaxy clusters whose gas eRASS1 catalogues as one extended source, variable AGN, X-ray binaries and cataclysmic variables in quiescence, transients), a Gaia DR3 star brighter than G 13 or within 500 pc within 20 arcsec, WISE AGN colours (W1-W2 above 0.8), a source of 10 ct/s within half a degree, three or more missing detections within half a degree (one extended source cut up), extent, and the ROSAT light curve's own variability; what is left is the catalogue's spurious tail when every such source lies below likelihood 20, against the share of the sample there.
+
+    python -I seti.py highenergy rosatgone [record] [--exi-min 10] [--crate-min 0.05] [--match-arcsec 120] [--edge-deg 1] [--simbad-arcsec 90] [--gaia-arcsec 20] [--bright-ct 10] [--bright-deg 0.5] [--clump-deg 0.5] [--spur-exi 20]"""
+
+# SIMBAD's types as CDS XMatch writes them, long names and short codes alike
+ROSATGONE_CLUSTER = ("ClG", "GtowardsCl", "BrightestCG", "GroupG", "GrG", "CGG", "SCG", "C?G", "SuperClG", "SNRemnant", "SNR")
+ROSATGONE_VARIABLE = ("AGN", "QSO", "Seyfert", "Sy1", "Sy2", "SyG", "BLLac", "Blazar", "BLL", "Bla", "LINER", "LIN", "XrayB", "XB", "LXB", "HXB", "CataclyV", "CV", "Nova", "No*",
+                      "SN", "GRB", "TDE", "Pulsar", "Psr")
+
+
+def rosatgone_kind(t):
+    """A SIMBAD type read as a galaxy cluster or remnant, a star, a variable extragalactic or accreting source, or none of these."""
+    if any(t.startswith(k) for k in ROSATGONE_CLUSTER):
+        return "cluster"
+    if any(k in t for k in ROSATGONE_VARIABLE):
+        return "variable"
+    if t.endswith("*") or t.startswith("Star") or t in ("**", "*", "Star?"):
+        return "star"
+    return "other" if t else None
+
+
+def rosatgone_main(record=False, exi_min=10.0, crate_min=0.05, match_arcsec=120.0, edge_deg=1.0, simbad_arcsec=90.0, gaia_arcsec=20.0, bright_ct=10.0, bright_deg=0.5, clump_deg=0.5, spur_exi=20.0):
+    path = core.data("rosatgone", "cat2rxs.tsv")
+    if not os.path.exists(path):
+        os.makedirs(os.path.dirname(path), exist_ok=True)
+        cols = "2RXS,ExiML,CRate,e_CRate,RAJ2000,DEJ2000,GLON,GLAT,Ext,FAll,Fluxp,sigExc"
+        open(path, "wb").write(core.get("https://vizier.cds.unistra.fr/viz-bin/asu-tsv?-source=J/A%2BA/588/A103/cat2rxs&-out.max=unlimited&-out=" + cols, timeout=900))
+    lines = [l for l in open(path, encoding="utf-8", errors="replace").read().splitlines() if l and not l.startswith("#")]
+    head = lines[0].split("\t")
+    rows = [dict(zip(head, l.split("\t"))) for l in lines[3:]]
+    src = []
+    for r in rows:
+        ra, dec, l, exi, cr = (num(r.get(k)) for k in ("RAJ2000", "DEJ2000", "GLON", "ExiML", "CRate"))
+        if None in (ra, dec, l, exi, cr) or not 180.0 + edge_deg <= l < 360.0 - edge_deg or exi < exi_min or cr < crate_min:
+            continue
+        src.append(dict(source_id=r["2RXS"].strip(), ra=ra, dec=dec, cr=cr, ecr=num(r.get("e_CRate")), exi=exi, ext=num(r.get("Ext")) or 0.0,
+                        flux=num(r.get("Fluxp")), sigexc=num(r.get("sigExc"))))
+    near = {}
+    for k in range(0, len(src), 20000):
+        for m in core.xmatch(src[k:k + 20000], "J/A+A/682/A34/erass1-m", match_arcsec):
+            d = float(m["angDist"])
+            near[m["source_id"]] = min(near.get(m["source_id"], 1e9), d)
+    gone = [s for s in src if near.get(s["source_id"], 1e9) > match_arcsec]
+    covered = gone
+    print("%d 2RXS point sources in the eRASS1 half at likelihood %g and %g ct/s or brighter; %d without an eRASS1 source within %g arcsec" % (
+        len(src), exi_min, crate_min, len(gone), match_arcsec), flush=True)
+    sim, gaia = {}, {}
+    for m in core.xmatch(covered, "simbad", simbad_arcsec):
+        if m["source_id"] not in sim or float(m["angDist"]) < sim[m["source_id"]][0]:
+            sim[m["source_id"]] = (float(m["angDist"]), (m.get("main_id") or "").strip(), (m.get("main_type") or m.get("otype") or "").strip())
+    for m in core.xmatch(covered, "I/355/gaiadr3", gaia_arcsec):
+        g, plx = num(m.get("Gmag")), num(m.get("Plx"))
+        # a bright star, or any star within 500 pc, of which M dwarfs flare by an order of magnitude
+        if g is not None and (g < 13 or (plx or 0) > 2.0) and (m["source_id"] not in gaia or g < gaia[m["source_id"]]):
+            gaia[m["source_id"]] = g
+    agn = {m["source_id"] for m in core.xmatch(covered, "II/328/allwise", gaia_arcsec)
+           if num(m.get("W1mag")) is not None and num(m.get("W2mag")) is not None and num(m["W1mag"]) - num(m["W2mag"]) > 0.8}
+    # a source ROSAT saw at bright_ct or more spoils its surroundings; three or more missing detections within clump_deg are one extended source cut up
+    bright = [x for x in src if x["cr"] >= bright_ct]
+    for s in covered:
+        s["bright"] = next((b["source_id"] for b in bright if b is not s and sep_arcsec(b["ra"], b["dec"], s["ra"], s["dec"]) < bright_deg * 3600), None)
+        s["clump"] = sum(1 for x in covered if sep_arcsec(x["ra"], x["dec"], s["ra"], s["dec"]) < clump_deg * 3600)
+    for s in covered:
+        s["simbad"] = sim.get(s["source_id"])
+        t = s["simbad"][2] if s["simbad"] else ""
+        s["gaia_g"] = gaia.get(s["source_id"])
+        s["vet"] = ("bright-neighbour" if s["bright"] else
+                    "diffuse-emission" if s["ext"] > 0 or s["clump"] >= 3 or rosatgone_kind(t) == "cluster" else
+                    "flare" if rosatgone_kind(t) == "star" or s["gaia_g"] is not None else
+                    "natural-type" if rosatgone_kind(t) == "variable" or s["source_id"] in agn else
+                    "variable-star" if s["sigexc"] and s["sigexc"] > 3 else
+                    "natural-type" if t else None)
+    # the rest are spurious detections when they all sit at the catalogue's lowest likelihoods, which few of the sample do
+    low = sum(x["exi"] < spur_exi for x in src) / max(len(src), 1)
+    rest = [s for s in covered if s["vet"] is None]
+    p_low = low ** sum(s["exi"] < spur_exi for s in rest) if rest else 1.0
+    if rest and all(s["exi"] < spur_exi for s in rest) and p_low < 1e-3:
+        for s in rest:
+            s["vet"] = "artefact"
+    print("%d left by the other vets, %d of them below likelihood %g, where %.0f per cent of the sample lies (p %.1e)" % (
+        len(rest), sum(s["exi"] < spur_exi for s in rest), spur_exi, 100 * low, p_low))
+    left = [s for s in covered if s["vet"] is None]
+    counts = {}
+    for s in covered:
+        counts[s["vet"]] = counts.get(s["vet"], 0) + 1
+    print("vets %s; %d unexplained" % (counts, len(left)))
+    for s in sorted(left, key=lambda s: -s["cr"])[:30]:
+        print("  2RXS %s: %.3f +- %.3f ct/s at RA %.4f Dec %+.4f, likelihood %.0f" % (s["source_id"], s["cr"], s["ecr"] or 0, s["ra"], s["dec"], s["exi"]))
+    if record:
+        core.record_run(
+            "rosat_erass1_gone", "X-ray sources that went out", "2RXS (VizieR J/A+A/588/A103) and eRASS1 main (J/A+A/682/A34/erass1-m) by CDS XMatch; SIMBAD; Gaia DR3",
+            "%d 2RXS point sources at likelihood %g and %g ct/s or brighter at Galactic longitude 180-360 deg" % (len(src), exi_min, crate_min), os.path.abspath(__file__), "rosatgone",
+            len(src), "sources", dict(exi_min=exi_min, crate_min=crate_min, match_arcsec=match_arcsec, edge_deg=edge_deg, simbad_arcsec=simbad_arcsec, gaia_arcsec=gaia_arcsec, bright_ct=bright_ct, bright_deg=bright_deg, clump_deg=clump_deg, spur_exi=spur_exi),
+            dict(gone=len(gone), vets=counts), "a 2RXS source ten times eRASS1's limit with no eRASS1 source within %g arcsec" % match_arcsec,
+            "%d of %d without an eRASS1 counterpart: %s; %d unexplained" % (len(covered), len(src), ", ".join("%d %s" % (v, k) for k, v in sorted(counts.items(), key=lambda x: -x[1]) if k), len(left)),
+            [dict(target="2rxs:%s" % s["source_id"], stat=round(s["cr"], 3), unit="ROSAT ct/s, absent from eRASS1", vet=s["vet"],
+                  note="likelihood %.0f; SIMBAD %s; Gaia G %s" % (s["exi"], s["simbad"], "%.1f" % s["gaia_g"] if s["gaia_g"] is not None else "-"))
+             for s in covered])
+        print("recorded")
+
+
 APPROACHES = {
+    "rosatgone": (lambda: core.run_cli(rosatgone_main), rosatgone_DOC),
     "gasunder": (lambda: core.run_cli(gasunder_main), gasunder_DOC),
     "hgpsoverlap": (lambda: core.run_cli(hgpsoverlap_main), hgpsoverlap_DOC),
     "unasclass": (lambda: core.run_cli(unasclass_main), unasclass_DOC),

setitools/infrared.py

@@ -3502,13 +3502,13 @@ def ksexcess_main(record=False, widths=5.0, errors=5.0, w1_widths=2.0, colour_mi
         print("recorded")
 
 
-warmgalaxies_DOC = r"""Galaxies warm without star formation: HECATE galaxies within 150 Mpc with a GSWLC stellar mass and ultraviolet-optical star-formation rate and with WISE 22 um and IRAS far-infrared luminosities, the quiescent ones (specific star-formation rate below 1e-11 per year) whose 22 um luminosity exceeds a tenth of a solar luminosity per solar mass and lies within a factor of three of their far-infrared luminosity; star-forming dust at 20-40 K peaks at 60-100 um and old stars warm little dust, so mid-infrared heat without young stars or cold dust is what structures at room temperature around a galaxy's stars would radiate. Each read for an active nucleus in HECATE's three flags and for an edge-on disc, whose dust lanes hide the young stars.
+warmgalaxies_DOC = r"""Galaxies warm without star formation: HECATE galaxies within 150 Mpc with a GSWLC stellar mass and ultraviolet-optical star-formation rate and with WISE 22 um and IRAS far-infrared luminosities, the quiescent ones (specific star-formation rate below 1e-11 per year) whose 22 um luminosity exceeds a tenth of a solar luminosity per solar mass and lies within a factor of three of their far-infrared luminosity; star-forming dust at 20-40 K peaks at 60-100 um and old stars warm little dust, so mid-infrared heat without young stars or cold dust is what structures at room temperature around a galaxy's stars would radiate. Each read for an active nucleus in HECATE's two flags, for an edge-on disc, whose dust lanes hide the young stars, and for HECATE's homogenised star-formation rate, which counts the infrared and puts a galaxy whose young stars dust hides from the ultraviolet and optical back among the star-forming.
 
     python -I seti.py infrared warmgalaxies [record] [--d-max-mpc 150] [--ssfr-max -11] [--heat-min -1.0] [--warm-min -0.5] [--edge-on-deg 70]"""
 
 
 def warmgalaxies_main(record=False, d_max_mpc=150.0, ssfr_max=-11.0, heat_min=-1.0, warm_min=-0.5, edge_on_deg=70.0):
-    rows = core.tap("vizier", 'SELECT PGC, NED, RAJ2000, DEJ2000, D, i, logLTIR, logLFIR, logL22um, logL12um, logMGSW, logSFRGSW, logSFRTIR, Nuclear, AGNs17, AGNh, W3mag, W4mag '
+    rows = core.tap("vizier", 'SELECT PGC, NED, RAJ2000, DEJ2000, D, i, logLTIR, logLFIR, logL22um, logL12um, logMGSW, logSFRGSW, logSFRTIR, logSFRh, Nuclear, AGNs17, AGNh, W3mag, W4mag '
                               'FROM "J/MNRAS/506/1896/hecate" WHERE D < %g AND logSFRGSW IS NOT NULL AND logMGSW IS NOT NULL AND logL22um IS NOT NULL AND logLFIR IS NOT NULL'
                     % d_max_mpc, timeout=900)
     gal = []
@@ -3518,8 +3518,8 @@ def warmgalaxies_main(record=False, d_max_mpc=150.0, ssfr_max=-11.0, heat_min=-1
             continue
         # HECATE's luminosities are in erg/s; a solar luminosity is 3.828e33 erg/s
         l22, lfir = l22 - 33.583, lfir - 33.583
-        gal.append(dict(pgc=r["PGC"], name=(r["NED"] or "").strip(), d=num(r["D"]), inc=num(r["i"]), m=m, ssfr=sfr - m, heat=l22 - m, warm=l22 - lfir, l22=l22, lfir=lfir,
-                        ltir=num(r["logLTIR"]), agn=bool((r["AGNs17"] or "").strip() or (r["AGNh"] or "").strip() or (r["Nuclear"] or "").strip() not in ("", "0")),
+        gal.append(dict(pgc=r["PGC"], name=(r["NED"] or "").strip(), d=num(r["D"]), inc=num(r["i"]), m=m, ssfr=sfr - m, ssfr_h=(num(r["logSFRh"]) - m) if num(r["logSFRh"]) is not None else None, heat=l22 - m, warm=l22 - lfir, l22=l22, lfir=lfir,
+                        ltir=num(r["logLTIR"]), agn="Y" in (str(r["AGNs17"]).strip(), str(r["AGNh"]).strip()),
                         w34=(num(r["W3mag"]) - num(r["W4mag"])) if num(r["W3mag"]) is not None and num(r["W4mag"]) is not None else None))
     quiet = [g for g in gal if g["ssfr"] < ssfr_max]
     forming = [g for g in gal if g["ssfr"] >= ssfr_max]
@@ -3527,14 +3527,15 @@ def warmgalaxies_main(record=False, d_max_mpc=150.0, ssfr_max=-11.0, heat_min=-1
     warm_q = np.array([g["warm"] for g in quiet])
     hits = [g for g in quiet if g["heat"] > heat_min and g["warm"] > warm_min]
     for g in hits:
-        g["vet"] = "background-galaxy" if g["agn"] else "dust-lane" if g["inc"] is not None and g["inc"] > edge_on_deg else None
+        g["vet"] = ("background-galaxy" if g["agn"] else "dust-lane" if g["inc"] is not None and g["inc"] > edge_on_deg else
+                    "dust" if g["ssfr_h"] is not None and g["ssfr_h"] > ssfr_max + 1 else None)
     left = [g for g in hits if g["vet"] is None]
     print("%d galaxies within %g Mpc with masses, ultraviolet-optical star-formation rates, 22 um and far-infrared luminosities; quiescent (specific rate below 1e%g): %d, "
           "whose 22 um heat per stellar mass runs median %.2f (99th percentile %.2f) and 22 um over far-infrared median %.2f (99th %.2f); above %g per solar mass and within "
-          "%g dex of the far-infrared: %d, %d with active nuclei, %d edge-on; unexplained %d" % (
+          "%g dex of the far-infrared: %d, %d with active nuclei, %d edge-on, %d star-forming by the infrared rate; unexplained %d" % (
               len(gal), d_max_mpc, ssfr_max, len(quiet), float(np.median(heat_q)) if len(heat_q) else 0, float(np.percentile(heat_q, 99)) if len(heat_q) else 0,
               float(np.median(warm_q)) if len(warm_q) else 0, float(np.percentile(warm_q, 99)) if len(warm_q) else 0, heat_min, warm_min, len(hits),
-              sum(g["vet"] == "background-galaxy" for g in hits), sum(g["vet"] == "dust-lane" for g in hits), len(left)))
+              sum(g["vet"] == "background-galaxy" for g in hits), sum(g["vet"] == "dust-lane" for g in hits), sum(g["vet"] == "dust" for g in hits), len(left)))
     for g in sorted(hits, key=lambda g: -g["heat"]):
         print("  %s (PGC %s, %s Mpc): log M %.2f, specific rate %.2f, 22 um per mass %.2f, 22 um over far-infrared %+.2f, W3-W4 %s, inclination %s; %s" % (
             g["name"], g["pgc"], g["d"], g["m"], g["ssfr"], g["heat"], g["warm"], "%.2f" % g["w34"] if g["w34"] is not None else "-", g["inc"], g["vet"] or "unexplained"))
@@ -3548,8 +3549,8 @@ def warmgalaxies_main(record=False, d_max_mpc=150.0, ssfr_max=-11.0, heat_min=-1
             "galaxies apart" % (float(np.median(heat_q)) if len(heat_q) else 0, float(np.percentile(heat_q, 99)) if len(heat_q) else 0,
                                 float(np.median(warm_q)) if len(warm_q) else 0, float(np.percentile(warm_q, 99)) if len(warm_q) else 0, len(forming)),
             "a quiescent galaxy radiating a tenth of a solar luminosity per solar mass at 22 um with little cold dust",
-            "%d: %d with active nuclei, %d edge-on; %d unexplained" % (len(hits), sum(g["vet"] == "background-galaxy" for g in hits), sum(g["vet"] == "dust-lane" for g in hits),
-                                                                     len(left)),
+            "%d: %d with active nuclei, %d edge-on, %d star-forming once the infrared is counted; %d unexplained" % (
+                len(hits), sum(g["vet"] == "background-galaxy" for g in hits), sum(g["vet"] == "dust-lane" for g in hits), sum(g["vet"] == "dust" for g in hits), len(left)),
             [dict(target="pgc:%s" % g["pgc"], stat=g["heat"], unit="log 22 um luminosity per stellar mass, solar units", vet=g["vet"],
                   note="%s, %s Mpc, log M %.2f, specific rate %.2f, 22 um over far-infrared %+.2f, inclination %s" % (g["name"], g["d"], g["m"], g["ssfr"], g["warm"], g["inc"]))
              for g in hits])
@@ -4644,7 +4645,105 @@ def pvedge_main(record=False, widths=5.0, jk_widths=2.0, emax=0.03, beam_arcsec=
         print("recorded")
 
 
+radioquiet_DOC = r"""Heat without cosmic rays: star formation heats dust and, through its supernovae, lights synchrotron, so the far-infrared to 1.4 GHz ratio q holds at 2.34 with 0.26 dex of scatter over four decades of luminosity (Yun, Reddy and Condon 2001), and waste heat from engines carries no synchrotron with it; HECATE galaxies within 200 Mpc north of -40 deg with IRAS 60 and 100 um fluxes, the far-infrared flux 1.26e-14 (2.58 S60 + S100) W/m^2, against every NVSS component within the galaxy's D25 semi-major axis, at least 45 arcsec and at most 3 arcmin, q above 3.3 (four times the scatter), or NVSS's 2.5 mJy giving a lower limit above it; each read for an IRAS source in the Galactic plane's cirrus (|b| below 10 deg), a disc wider than 6 arcmin, whose synchrotron the 3 arcmin sum misses, HECATE's active-nucleus flags or WISE AGN colours (W1-W2 above 0.8), NED's largest 1.4 GHz flux integrated over the galaxy's map, which holds a faint disc NVSS's component list breaks up, SIMBAD's active-nucleus types, NGC 1377's published nascent starburst, and a 60 um flux IRAS flags as an upper limit.
+
+    python -I seti.py infrared radioquiet [record] [--d-max-mpc 200] [--q-flag 3.3] [--nvss-limit-mjy 2.5] [--b-min 10] [--extent-arcmin 6]"""
+
+
+# galaxies the literature has already explained: PGC to reading
+RADIOQUIET_KNOWN = {"13324": "a nascent starburst too young for its supernovae to have lit synchrotron (Roussel et al. 2003)"}
+
+
+def radioquiet_ned_1p4(name):
+    """The largest integrated 1.4 GHz flux density (mJy) NED holds for a galaxy, or None."""
+    t = core.get("https://ned.ipac.caltech.edu/cgi-bin/datasearch?" + urllib.parse.urlencode(dict(search_type="Photometry", objname=name, of="ascii_bar", meas_type="bot")),
+                 timeout=120, soft=True)
+    best = None
+    for line in (t or b"").decode(errors="replace").splitlines():
+        c = [x.strip() for x in line.split("|")]
+        # columns: number, band, value, uncertainty, units, frequency (Hz), flux in Jy
+        # the continuum alone: the 21 cm line sits at 1.42 GHz in Jy Hz
+        if len(c) > 6 and c[4] in ("mJy", "Jy") and "HI" not in c[1] and "21 cm" not in c[1] and num(c[5]) and abs(num(c[5]) - 1.4e9) < 0.03e9 and num(c[6])                 and "<" not in c[3]:
+            best = max(best or 0.0, 1000.0 * num(c[6]))
+    return best
+
+
+def radioquiet_main(record=False, d_max_mpc=200.0, q_flag=3.3, nvss_limit_mjy=2.5, b_min=10.0, extent_arcmin=6.0):
+    from astropy.coordinates import SkyCoord
+    import astropy.units as u
+    rows = core.tap("vizier", 'SELECT PGC, NED, RAJ2000, DEJ2000, D, MajAxis, F60um, F100um, q_F60um, q_F100um, W1mag, W2mag, logSFRh, Nuclear, AGNs17, AGNh '
+                              'FROM "J/MNRAS/506/1896/hecate" WHERE D < %g AND DEJ2000 > -40 AND F60um > 0 AND F100um > 0' % d_max_mpc, timeout=900)
+    gal = []
+    for r in rows:
+        s60, s100 = num(r["F60um"]), num(r["F100um"])
+        fir = 1.26e-14 * (2.58 * s60 + s100)
+        # testable only where NVSS's limit would put q above the flag
+        if math.log10(fir / 3.75e12 / (nvss_limit_mjy * 1e-29)) <= q_flag:
+            continue
+        maj = num(r["MajAxis"]) or 0.0
+        gal.append(dict(source_id=str(r["PGC"]), name=(r["NED"] or "").strip(), ra=float(r["RAJ2000"]), dec=float(r["DEJ2000"]), d=num(r["D"]), maj=maj, fir=fir,
+                        upper=str(r["q_F60um"]).strip() == "1", agn="Y" in (str(r["AGNs17"]).strip(), str(r["AGNh"]).strip()),
+                        w12=(num(r["W1mag"]) - num(r["W2mag"])) if num(r["W1mag"]) is not None and num(r["W2mag"]) is not None else None, s1=0.0))
+    # HECATE's MajAxis is the D25 semi-major axis in arcmin
+    radius = {g["source_id"]: max(45.0, min(60.0 * g["maj"], 180.0)) for g in gal}
+    byid = {g["source_id"]: g for g in gal}
+    for k in range(0, len(gal), 20000):
+        for m in core.xmatch(gal[k:k + 20000], "VIII/65/nvss", 180.0):
+            if float(m["angDist"]) <= radius[m["source_id"]]:
+                byid[m["source_id"]]["s1"] += num(m.get("S1.4")) or 0.0
+    b = SkyCoord([g["ra"] for g in gal] * u.deg, [g["dec"] for g in gal] * u.deg).galactic.b.deg
+    for g, bb in zip(gal, b):
+        g["b"] = float(bb)
+        g["limit"] = g["s1"] <= 0
+        g["q"] = math.log10(g["fir"] / 3.75e12 / (max(g["s1"], nvss_limit_mjy) * 1e-29))
+    # NVSS's component list breaks a faint extended disc into nothing; NED holds the flux integrated over the map (Condon, Cotton and Broderick 2002)
+    for g in [g for g in gal if g["q"] > q_flag]:
+        g["ned"] = radioquiet_ned_1p4(g["name"]) if g["name"] else None
+        if g["ned"] and g["ned"] > g["s1"]:
+            g["q_ned"] = math.log10(g["fir"] / 3.75e12 / (g["ned"] * 1e-29))
+    qs = np.array([g["q"] for g in gal if not g["limit"]])
+    high = [g for g in gal if g["q"] > q_flag]
+    simbad = {m["source_id"]: (m.get("main_type") or "").strip() for m in core.xmatch(high, "simbad", 30.0)} if high else {}
+    for g in high:
+        # NVSS's components are summed within 3 arcmin of the centre, so a disc wider than twice that loses its extended synchrotron
+        g["vet"] = ("catalogue-dropped" if g.get("q_ned") is not None and g["q_ned"] <= q_flag else
+                    "diffuse-emission" if abs(g["b"]) < b_min else
+                    "out-of-reach" if 2 * g["maj"] > extent_arcmin else
+                    "published-classification" if g["agn"] or (g["w12"] is not None and g["w12"] > 0.8) else
+                    "natural-type" if any(k in simbad.get(g["source_id"], "") for k in ("AGN", "Seyfert", "LINER", "QSO", "Sy")) else
+                    "young-star" if g["source_id"] in RADIOQUIET_KNOWN else
+                    "single-record" if g["upper"] else None)
+    left = [g for g in high if g["vet"] is None]
+    print("%d HECATE galaxies within %g Mpc north of -40 deg bright enough in the far infrared that NVSS reaches q %g; %d detected, q median %.2f, robust scatter %.2f; "
+          "above %g: %d (%d upper limits), %d active nuclei, %d in the plane, %d larger than %g arcmin, %d with a 60 um limit; unexplained %d" % (
+              len(gal), d_max_mpc, q_flag, len(qs), float(np.median(qs)), 1.4826 * float(np.median(np.abs(qs - np.median(qs)))), q_flag, len(high), sum(g["limit"] for g in high),
+              sum(g["vet"] == "published-classification" for g in high), sum(g["vet"] == "diffuse-emission" for g in high), sum(g["vet"] == "out-of-reach" for g in high), extent_arcmin,
+              sum(g["vet"] == "single-record" for g in high), len(left)))
+    for g in sorted(high, key=lambda g: -g["q"])[:40]:
+        print("  %s (PGC %s, %.0f Mpc): q %s%.2f, NVSS %.1f mJy within %.0f arcsec, NED %s mJy%s, major axis %.1f arcmin, b %+.1f; %s" % (
+            g["name"], g["source_id"], g["d"] or 0, ">" if g["limit"] else "", g["q"], g["s1"], radius[g["source_id"]], "%.1f" % g["ned"] if g.get("ned") else "-",
+            " (q %.2f)" % g["q_ned"] if g.get("q_ned") is not None else "", g["maj"], g["b"], g["vet"] or "unexplained"))
+    if record:
+        core.record_run(
+            "fir_radio_quiet", "Heat without cosmic rays", "HECATE (J/MNRAS/506/1896) with its IRAS fluxes; NVSS (VIII/65) by CDS XMatch",
+            "%d HECATE galaxies within %g Mpc north of -40 deg with IRAS 60 and 100 um fluxes bright enough for NVSS to reach q %g" % (len(gal), d_max_mpc, q_flag),
+            os.path.abspath(__file__), "radioquiet", len(gal), "galaxies", dict(d_max_mpc=d_max_mpc, q_flag=q_flag, nvss_limit_mjy=nvss_limit_mjy, b_min=b_min, extent_arcmin=extent_arcmin),
+            dict(detected=len(qs), q_median=float(np.median(qs)), q_scatter=1.4826 * float(np.median(np.abs(qs - np.median(qs))))),
+            "q above %g, four times the natural relation's 0.26 dex above its 2.34" % q_flag,
+            "%d above q %g of %d (%d of them NVSS limits): %d back under it with NED's flux integrated over the map%s, %d larger than NVSS recovers, %d in the plane's cirrus, %d active nuclei, %d nascent starburst (NGC 1377), %d on IRAS limits; "
+            "%d unexplained%s" % (
+                len(high), q_flag, len(gal), sum(g["limit"] for g in high), sum(g["vet"] == "catalogue-dropped" for g in high),
+                " (%s)" % ", ".join("%s, %.0f mJy, q %.2f" % (g["name"], g["ned"], g["q_ned"]) for g in high if g["vet"] == "catalogue-dropped") if any(g["vet"] == "catalogue-dropped" for g in high) else "",
+                sum(g["vet"] == "out-of-reach" for g in high), sum(g["vet"] == "diffuse-emission" for g in high),
+                sum(g["vet"] in ("published-classification", "natural-type") for g in high), sum(g["vet"] == "young-star" for g in high), sum(g["vet"] == "single-record" for g in high),
+                len(left), " (%s)" % ", ".join("%s, q above %.2f" % (g["name"], g["q"]) for g in left) if left else ""),
+            [dict(target="pgc:%s" % g["source_id"], stat=round(g["q"], 2), unit="q, far-infrared over 1.4 GHz" + (" (lower limit)" if g["limit"] else ""), vet=g["vet"],
+                  note="%s, %s Mpc; NVSS %.1f mJy; major axis %.1f arcmin; b %+.1f" % (g["name"], g["d"], g["s1"], g["maj"], g["b"])) for g in high])
+        print("recorded")
+
+
 APPROACHES = {
+    "radioquiet": (lambda: core.run_cli(radioquiet_main), radioquiet_DOC),
     "pvedge": (lambda: core.run_cli(pvedge_main), pvedge_DOC),
     "coldcomp": (lambda: core.run_cli(coldcomp_main), coldcomp_DOC),
     "diskvoids": (lambda: core.run_cli(diskvoids_main), diskvoids_DOC),