Package: pyspectools
PySpecToolsA Python library for analysis of rotational spectroscopy and beyondIntroductionPySpecToolsis a library written to help with analyzing rotational
spectroscopy data. The main functions of this library are:Wrapper for SPFIT and SPCAT programs of Herb Pickett, with YAML/JSON
interpretationGenerating specific figure types usingmatplotlib, such as polyads and
potential energy diagramsParsing and filtering of Fourier-transform centimeter-wave and
millimeter-wave absorption data. This includes:Fitting of lineshapes (e.g. Lorentizan second-derivative profiles)Fourier-filteringDouble resonance fittingAnalysis of broad band spectra with theAssignmentSessionandTransitionclasses.
These classes, combined with Jupyter notebooks, provide a way to assign spectra
reproducibly; astronomical and laboratory broadband spectra are supported.Autofit routines are available for a set of special cases, like linear/prolate
molecules. Eventually, SPFIT will be a backend option.Molecule identity inference (NEW!)—this uses a pre-trained probabilistic deep
learning model that allows users to perform inference on experimental constants
and expected composition to predict the most likely molecular formula and possible
functional groups present. Seeour paper on the development of the first generation of this model. An example of how to run this analysis
can be foundhere.The documentation for PySpecTools can be foundhere.If you use PySpecTools for research, please cite use the DOI badge below to cite the version
of the package; this is not to track adoption, but rather for the sake of reproducibility!Installationcondais the preferred way of maintaining software environments, andpoetryis used to manage Python package dependencies.As of PySpecTools 4.6.1, the installation is intended to be significantly more straightforward
with PyPI releases; in a given Python environment, just run:pip install PySpecToolsAlternatively, if you're having issues, we recommend creating a new Python environment
withinconda; with Python 3.7+:conda create -n pst python=3.7conda activate pstpip install poetrypoetry installInstallation on Windows is less straightforward. The following instructions avoid
issues originating from virtual environments created by poetry and include a workaround
for a known issue with poetry in Windows.conda create -n pst python=3.7conda activate pstpip install poetrypoetry config virtualenvs.in-project falsepoetry config virtualenvs.create falseNavigate to the folderC:\Users\user\AppData\Local\pypoetry\Cacheand delete all contents of this folder.Navigate to the folder containing PySpecToolspoetry installPyPickettPySpecToolsincludes a set of routines for wrapping SPFIT/SPCAT. The design
philosophy behind these functions is that the formatting and running of
SPFIT/SPCAT can be a little tricky, in terms of reproducibility, parameter
coding, and visualization. These problems are solved by wrapping and managing
input files in an object-oriented fashion:Able to serialize SPFIT/SPCAT files from more human-friendly formats like
YAML and JSON.Automatic file/folder management, allowing the user to go back to an earlier
fit/parameters. Ability to "finalize" the fit so the final parameter set is
known.Display the predicted spectrum usingmatplotlibin a Jupyter notebook,
which could be useful for analysis and publication.A parameter scan mode, allowing the RMS to be explored as a function of
whatever parameter.There is still much to do for this module, including a way of managing experimental lines.Notes on releasePySpecToolsis currently being released on a irregular schedule, using a sequence-based version numbering system.
The numbering works as X.Y.Z, where X denotes huge changes that are backwards incompatible, Y are significant changes
(typically new features) and Z are minor bug fixes. A freeze and release will typically occur when
a new version with potentially backwards breaking changes are about to be made. The large changes typically occur once a year (based on the trend so far).Currently,PySpecToolsis under the MIT license, which allows anyone to freely use and modify as you wish!Planned featuresIntegration of deep learning tools for molecule identifiction and spectral assignmentProbability-based assignment routines - rather than single assignments.Revamp of codebase - needs a substantial re-organization that will likely result in backwards compatibility breaking.Additional Cython routines - many functions withinPySpecToolsare fast enough, but we can always go faster 😀Better abstraction in thespectra.assignmentmodules - need to move a lot of the complicated routines into subclasses (especially for transitions and line lists), although there is a case to be made for a simpler user interface (only have to deal withLineList, instead of three subclasses ofLineList)ContributingIf you have features you think would benefit other spectroscopists, you can raise an issue in the repo. Alternatively (and even better) would be to fork the repo, and submit a pull request!The only comments on coding style are:Documentation is written in NumPy styleObject-oriented PythonFormatted withblackThere are a set of unit tests that can be run to ensure the most complicated routines in the
library are working as intended. Right now coverage is poor, and so the tests I've written
focus on theassignmentmodule. There is a script contained in thetestsfolder that will
generate a synthetic spectrum to test functionality on. To run these tests:cdtestspythongenerate_test_spectrum.pypytestYou will need to havepytestinstalled. These tests are designed to raise errors when there
are huge errors; some tolerance is included for imperfect peak detection, for example. These
are defined as constants within thetest_assignment.pytesting script.Questions? Comments?If you have features you would like to have added, please raise an issue on the repo, or
feel free to send me an email at kinlee_at_cfa.harvard.edu.Also, please feel free to fork and contribute! The code is being formatted withblack,
and uses NumPy-style docstrings. If you have any questions about contributing, drop me an
email!