-
Notifications
You must be signed in to change notification settings - Fork 0
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
DM-35404: Add tests for focus analyses #18
Merged
Merged
Changes from all commits
Commits
Show all changes
3 commits
Select commit
Hold shift + click to select a range
File filter
Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,132 @@ | ||
# This file is part of summit_extras. | ||
# | ||
# Developed for the LSST Data Management System. | ||
# This product includes software developed by the LSST Project | ||
# (https://www.lsst.org). | ||
# See the COPYRIGHT file at the top-level directory of this distribution | ||
# for details of code ownership. | ||
# | ||
# This program is free software: you can redistribute it and/or modify | ||
# it under the terms of the GNU General Public License as published by | ||
# the Free Software Foundation, either version 3 of the License, or | ||
# (at your option) any later version. | ||
# | ||
# This program is distributed in the hope that it will be useful, | ||
# but WITHOUT ANY WARRANTY; without even the implied warranty of | ||
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the | ||
# GNU General Public License for more details. | ||
# | ||
# You should have received a copy of the GNU General Public License | ||
# along with this program. If not, see <https://www.gnu.org/licenses/>. | ||
|
||
import unittest | ||
from typing import Iterable | ||
|
||
import lsst.utils.tests | ||
|
||
import matplotlib as mpl | ||
mpl.use('Agg') | ||
|
||
from lsst.summit.extras import SpectralFocusAnalyzer, NonSpectralFocusAnalyzer # noqa: E402 | ||
import lsst.summit.utils.butlerUtils as butlerUtils # noqa: E402 | ||
|
||
|
||
class FocusAnalysisTestCase(lsst.utils.tests.TestCase): | ||
|
||
@classmethod | ||
def setUpClass(cls): | ||
try: | ||
cls.butler = butlerUtils.makeDefaultLatissButler() | ||
except FileNotFoundError: | ||
raise unittest.SkipTest("Skipping tests that require the LATISS butler repo.") | ||
|
||
# dataIds have been chosen match those of a spectral focus sweep which | ||
# is in the LATISS-test-data collection at the summit, TTS and USDF | ||
cls.dayObs = 20220628 | ||
cls.seqNums = range(280, 288+1) | ||
cls.focusAnalyzer = SpectralFocusAnalyzer() | ||
|
||
# default is 3 boxes, so setting four tests the generality of the code. | ||
# The values have been picked to land on the pretty small spectrum we | ||
# have in the test data. | ||
cls.focusAnalyzer.setSpectrumBoxOffsets([1600, 1700, 1800, 1900]) | ||
|
||
def test_run(self): | ||
# we don't check the plots, but set doDisplay to True to check the | ||
# plots are generated without error | ||
self.focusAnalyzer.getFocusData(self.dayObs, self.seqNums, doDisplay=True) | ||
result = self.focusAnalyzer.fitDataAndPlot() | ||
self.assertIsInstance(result, Iterable) | ||
self.assertEqual(len(result), len(self.focusAnalyzer.getSpectrumBoxOffsets())) | ||
|
||
for number in result: | ||
# check they're all numbers, non-nan, and vaguely sensible. These | ||
# are values in mm of the hexapod offset from nominal focus and are | ||
# usually very small (around 0.01mm) so this is just testing that | ||
# something hasn't gone wildly wrong with the fit | ||
self.assertGreater(number, -2.0) | ||
self.assertLess(number, 2.0) | ||
|
||
|
||
class NonSpectralFocusAnalysisTestCase(lsst.utils.tests.TestCase): | ||
|
||
@classmethod | ||
def setUpClass(cls): | ||
try: | ||
cls.butler = butlerUtils.makeDefaultLatissButler() | ||
except FileNotFoundError: | ||
raise unittest.SkipTest("Skipping tests that require the LATISS butler repo.") | ||
|
||
# dataIds have been chosen match those of a non-spectral focus sweep | ||
# which is in the LATISS-test-data collection at the summit, TTS and | ||
# USDF | ||
cls.dayObs = 20220405 | ||
cls.seqNums = range(523, 531+1) | ||
cls.focusAnalyzer = NonSpectralFocusAnalyzer() | ||
|
||
def test_run(self): | ||
# we don't check the plots, but set doDisplay to True to check the | ||
# plots are generated without error | ||
self.focusAnalyzer.getFocusData(self.dayObs, self.seqNums, doDisplay=True) | ||
result = self.focusAnalyzer.fitDataAndPlot() | ||
self.assertIsInstance(result, dict) | ||
|
||
# result is a dict which looks like this: | ||
# {'fwhmFitMin': 0.029221417454391708, | ||
# 'ee90FitMin': 0.0754762884665358, | ||
# 'ee80FitMin': 0.07188778363981346, | ||
# 'ee50FitMin': 0.13998855716378267} | ||
self.assertEqual(len(result), 4) | ||
|
||
for k, v in result.items(): | ||
# check they're all numbers, non-nan, and vaguely sensible. The | ||
# values are the position of the fwhm which should be close to zero | ||
# on the x-axis assuming we were roughly symmetric around focus | ||
# (which is the case for the current test dataset) and encircled | ||
# energy radii, which also should be small as the data in question | ||
# was roughly in focus. This is really just to check that the code | ||
# isn't horribly broken, so not much attention should be paid | ||
# to these numerical values. | ||
self.assertIsInstance(k, str) | ||
self.assertIsInstance(v, float) | ||
|
||
if k == 'fwhmFitMin': | ||
# everything else is strictly positive but the fwhm position | ||
# can be negative. | ||
v = abs(v) | ||
|
||
self.assertGreater(v, 0.0) | ||
self.assertLess(v, 1.0) | ||
|
||
|
||
class TestMemory(lsst.utils.tests.MemoryTestCase): | ||
pass | ||
|
||
|
||
def setup_module(module): | ||
lsst.utils.tests.init() | ||
|
||
|
||
if __name__ == "__main__": | ||
lsst.utils.tests.init() | ||
unittest.main() |
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Minor point, but shouldn't this import happen after
lsst.summit.utils
?There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Is this referring to fact I'm importing the tested code after the
mpl.use('Agg')
line? This is because you have to set the matplotlib backend before any other code which imports matplotlib does (hence the#noqa: E402
line decorators) or people have plots pop up.There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
No, this is a very minor point. Shouldn't the order of imports be:
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Ohhhh, there's no
isort
over here, and I've made minimal effort anywhere to keep things alphabetical, so that didn't even occur to me! Look elsewhere in the package and you'll see the imports are a free-for-all! I'm considering going pre-commit + black + isort etc at some point though, but that's definitely another ticket 🙂