snr_check.py

Computes the SNRs on the secondary eclipses of all TRAPPIST-1 planets in the 15 micron JWST MIRI filter. This is used for benchmarking the code. The benchmarked values are

    planetplanet    analytic     ratio
b:     4.250          4.488      0.947
c:     2.752          2.950      0.933
d:     0.931          0.978      0.953
e:     0.863          0.900      0.958
f:     0.652          0.691      0.942
g:     0.502          0.539      0.931
h:     0.117          0.122      0.958

Note that the planetplanet estimates are lower than the analytic estimates by a few percent because the latter neglects thermal/background noise.