Here's some data mining/plotting tools for looking at Herschel Space Telescope line emissions!

Basic overview:

First, the herschel_correlations.py file computes correlations from the CDF_archive 1d lines files (either SPIRE, PACS, or the merged file I created myself). This produces a large CSV file of line emission pairs. It takes a CDF_archive 1d file and a CSV file with target properties (i.e. type, Lbol, etc ). The target properties file is also used for some other correlations. The list of acceptable types is hardcoded into the python file for now (that can be easily changed before re-running the script). An example of a properties file is included.

use: herschel_correlations.py (CDF lines filepath) (target property file) (output file path)

plot_lines.py also reads from the CDF_archive 1d lines file, and produces a scatter plot of one line vs another line with targets labeled.

use: plot_lines.py (CDF lines filepath) (target property file) (line to plot) (other line to plot)

The other functions, such as compare_prefixes.py, iterate through the large CSV file to produce subsets matching a specific pattern, such as all pairs where one is a CO line and the other is an OH line.

use: compare_prefixes.py (correlations file generated by herschel_correlations) (prefix 1, ex. OH or CO) (prefix 2) (correlation threshold, ex. everything with a correlation less than 0.9 will be omitted)

use: compare_like_lines.py (correlations file) (prefix for intra-prefix correlations) (threshold)

use line_correlations.py (correlations file) (specific line name, ex. CO12-11) (threshold)

CO/13CO vs J plotting is included in the herschel_mining.py, but there's no command-line interface yet. Feel free to import it into another script and use it!