From b7da4dbcebbc211496e56d34b0dbe6b3fd293ce3 Mon Sep 17 00:00:00 2001
From: Kristopher Cooper <kris@infra04-wg012.eduroam.wireless.umn.edu>
Date: Tue, 23 Sep 2025 17:04:09 -0500
Subject: [PATCH 1/3] Removed build folder

---
 .../build/lib/response_tools_py/__init__.py   |  34 -
 .../build/lib/response_tools_py/att.py        | 242 -------
 .../build/lib/response_tools_py/det_resp.py   | 310 ---------
 .../build/lib/response_tools_py/eff_area.py   | 320 ---------
 .../lib/response_tools_py/io/__init__.py      |   0
 .../lib/response_tools_py/io/load_yaml.py     |  19 -
 .../response_tools_py/io/tests/__init__.py    |   0
 .../build/lib/response_tools_py/ir-plots.py   | 641 ------------------
 .../build/lib/response_tools_py/phot_spec.py  | 100 ---
 .../lib/response_tools_py/quantum_eff.py      |  50 --
 .../build/lib/response_tools_py/version.py    |   3 -
 11 files changed, 1719 deletions(-)
 delete mode 100644 response-tools-py/build/lib/response_tools_py/__init__.py
 delete mode 100644 response-tools-py/build/lib/response_tools_py/att.py
 delete mode 100644 response-tools-py/build/lib/response_tools_py/det_resp.py
 delete mode 100644 response-tools-py/build/lib/response_tools_py/eff_area.py
 delete mode 100644 response-tools-py/build/lib/response_tools_py/io/__init__.py
 delete mode 100644 response-tools-py/build/lib/response_tools_py/io/load_yaml.py
 delete mode 100644 response-tools-py/build/lib/response_tools_py/io/tests/__init__.py
 delete mode 100644 response-tools-py/build/lib/response_tools_py/ir-plots.py
 delete mode 100644 response-tools-py/build/lib/response_tools_py/phot_spec.py
 delete mode 100644 response-tools-py/build/lib/response_tools_py/quantum_eff.py
 delete mode 100644 response-tools-py/build/lib/response_tools_py/version.py

diff --git a/response-tools-py/build/lib/response_tools_py/__init__.py b/response-tools-py/build/lib/response_tools_py/__init__.py
deleted file mode 100644
index 3f5cdf0..0000000
--- a/response-tools-py/build/lib/response_tools_py/__init__.py
+++ /dev/null
@@ -1,34 +0,0 @@
-"""
-Top level of the `response-tools-py` package.
-
-Namespace
----------
-
-The base `response_tools_py` namespace includes:
-
-`~response_tools_py.contextResponseInfo`
-    The information stored in a YAML file that includes information from
-    the flight.
-
-`~response_tools_py.__version__`
-    The current version of the code as stated in the setup script.
-
-Examples
---------
-# importing the module is as easy as:
->>> import response_tools_py
-
-# then accessing, e.g., the context information for the flight/flare
->>> print(response_tools_py.contextResponseInfo)
-...
-
-# that can be accessed like a native Python dictionary
->>> print(response_tools_py.contextResponseInfo[""])
-...
-"""
-
-from .version import __version__
-from response_tools_py.io.load_yaml import load_response_context
-
-# for global context info
-contextResponseInfo = load_response_context()
diff --git a/response-tools-py/build/lib/response_tools_py/att.py b/response-tools-py/build/lib/response_tools_py/att.py
deleted file mode 100644
index e7ed8e8..0000000
--- a/response-tools-py/build/lib/response_tools_py/att.py
+++ /dev/null
@@ -1,242 +0,0 @@
-"""Code to load different attenuators. """
-
-import logging
-import os
-import pathlib
-
-from astropy.io import fits
-import astropy.units as u
-import numpy as np
-import pandas
-import scipy
-
-# thermal blanket attenuation
-@u.quantity_input(mid_energies=u.keV)
-def att_thermal_blanket(mid_energies, file=None):
-    """Return thermal blanket transmittance interpolated to the given energies."""
-    _f = os.path.join(pathlib.Path(__file__).parent, 
-                            "att-data", 
-                            "F4_Blanket_transmission.dat") if file is None else file
-    att = scipy.io.readsav(_f)
-    att_es, att_values = att["energy_kev"] << u.keV, att["f4_transmission"] << u.dimensionless_unscaled
-    return np.interp(mid_energies.value, att_es.value, att_values.value, left=0, right=1) << u.dimensionless_unscaled
-
-# thermal blanket attenuation
-@u.quantity_input(mid_energies=u.keV)
-def att_uniform_al_cdte(mid_energies, file=None, position=None):
-    """Return thermal blanket transmittance interpolated to the given energies."""
-    if position is None:
-        logging.info("`position` must be 2 or 4.")
-    _f = os.path.join(pathlib.Path(__file__).parent, 
-                            "att-data", 
-                            f"unif_att_p{position}_theoretical.csv") if file is None else file
-    att = pandas.read_csv(_f)
-    att_es, att_values = att["energy[keV]"] << u.keV, att["transmission"] << u.dimensionless_unscaled
-    return np.interp(mid_energies.value, att_es.value, att_values.value, left=0, right=1) << u.dimensionless_unscaled
-
-# pixelated attenuator attenuation
-@u.quantity_input(mid_energies=u.keV)
-def att_pixelated(mid_energies, file=None, values=None):
-    """Return pixelated attenuator transmittance interpolated to the given energies."""
-    _f = os.path.join(pathlib.Path(__file__).parent, 
-                            "att-data", 
-                            "20240607_fosxi4_transmission.csv") if file is None else file
-    att = pandas.read_csv(_f)
-    att_es, att_values_measured, att_values_modelled = att["energy"] << u.keV, att["measured_transmission"] << u.dimensionless_unscaled, att["modeled_transmission"] << u.dimensionless_unscaled
-
-    if values=="measured":
-        return np.interp(mid_energies.value, att_es.value, att_values_measured.value, left=0, right=1) << u.dimensionless_unscaled
-    elif values in ["modelled", "modeled"]:
-        return np.interp(mid_energies.value, att_es.value, att_values_modelled.value, left=0, right=1) << u.dimensionless_unscaled
-    else:
-        logging.info("`values` must be \"measured\" or \"modelled\" (or \"modeled\" for you Americans).")
-
-# aluminized mylar attenuation
-@u.quantity_input(mid_energies=u.keV)
-def att_al_mylar(mid_energies, file=None):
-    """Return aluminized mylar transmittance interpolated to the given energies."""
-    _f = os.path.join(pathlib.Path(__file__).parent, 
-                            "att-data", 
-                            "thin_mylar_p3_p5_theoretical.csv") if file is None else file
-    att = pandas.read_csv(_f)
-    att_es, att_values = att["energy[keV]"] << u.keV, att["transmission"] << u.dimensionless_unscaled
-    return np.interp(mid_energies.value, att_es.value, att_values.value, left=0, right=1) << u.dimensionless_unscaled
-
-# OBF attenuation
-@u.quantity_input(mid_energies=u.keV)
-def att_obf0(mid_energies, file=None):
-    """Return OBF0 transmittance interpolated to the given energies."""
-    logging.warning("This might not be the function you are looking for, please see `att_cmos_obfilter`.")
-    _f = os.path.join(pathlib.Path(__file__).parent, 
-                            "att-data", 
-                            "CMOST_Prefilter_transmission.dat") if file is None else file
-    att = scipy.io.readsav(_f)
-    att_es, att_values = att["cmos_prefilter_transmission"]["energy_kev"][0], att["cmos_prefilter_transmission"]["position0"][0]
-
-    # this attenuation acts weird because it stops noticably above 0 at 1 keV (the lower energy), so...
-    # assume same drop-off as OBF1
-    _obf1_numbers = att_obf1(mid_energies, file=None).value
-    _ex_obf0_inds = np.nonzero(_obf1_numbers<att_values[0])
-    _energy_ext = mid_energies[_ex_obf0_inds].value
-    att_es = np.concatenate((_energy_ext-np.max(_energy_ext)+att_es[0], att_es)) << u.keV
-    att_values = np.concatenate((_obf1_numbers[_ex_obf0_inds], att_values)) << u.dimensionless_unscaled
-
-    return np.interp(mid_energies.value, att_es.value, att_values.value, left=0, right=1) << u.dimensionless_unscaled
-    
-# OBF attenuation
-@u.quantity_input(mid_energies=u.keV)
-def att_obf1(mid_energies, file=None):
-    """Return OBF1 transmittance interpolated to the given energies."""
-    logging.warning("This might not be the function you are looking for, please see `att_cmos_obfilter`.")
-    _f = os.path.join(pathlib.Path(__file__).parent, 
-                            "att-data", 
-                            "CMOST_Prefilter_transmission.dat") if file is None else file
-    att = scipy.io.readsav(_f)
-    att_es, att_values = att["cmos_prefilter_transmission"]["energy_kev"][0] << u.keV, att["cmos_prefilter_transmission"]["position1"][0] << u.dimensionless_unscaled
-    return np.interp(mid_energies.value, att_es.value, att_values.value, left=0, right=1) << u.dimensionless_unscaled
-
-@u.quantity_input(mid_energies=u.keV)
-def att_sigmoid(mid_energies, l ,x0, k, b):
-    """Sigmoid model for a general/fake attenuator."""
-    return (l / (1 + np.exp(-k*(mid_energies.value-x0))) + b) << u.dimensionless_unscaled
-
-@u.quantity_input(mid_energies=u.keV)
-def att_cmos_filter(mid_energies, file=None, telescope=None):
-    if telescope is None:
-        logging.warning("`telescope` input in `att_cmos_filter()` must be 0 or 1.")
-        return
-    _f = os.path.join(pathlib.Path(__file__).parent, 
-                            "att-data", 
-                            f"foxsi4_telescope-{telescope}_BASIC_attenuation_filter_transmittance_V25APR13.fits") if file is None else file
-    with fits.open(_f) as hdul:
-        es, t = hdul[2].data << u.keV, hdul[1].data << u.dimensionless_unscaled
-    return np.interp(mid_energies.value, es.value, t.value, left=0, right=1) << u.dimensionless_unscaled
-
-@u.quantity_input(mid_energies=u.keV)
-def att_cmos_obfilter(mid_energies, file=None, telescope=None):
-    if telescope is None:
-        logging.warning("`telescope` input in `att_cmos_obfilter()` must be 0 or 1.")
-        return
-    _f = os.path.join(pathlib.Path(__file__).parent, 
-                            "att-data", 
-                            f"foxsi4_telescope-{telescope}_BASIC_optical_blocking_filter_transmittance_V25APR13.fits") if file is None else file
-    with fits.open(_f) as hdul:
-        es, t = hdul[2].data << u.keV, hdul[1].data << u.dimensionless_unscaled
-    return np.interp(mid_energies.value, es.value, t.value, left=0, right=1) << u.dimensionless_unscaled
-
-@u.quantity_input(off_axis_angle=u.arcsec)
-def att_cmos_collimator_ratio(off_axis_angle, file=None, telescope=None):
-    if telescope is None:
-        logging.warning("`telescope` input in `att_cmos_collimator_ratio()` must be 0 or 1.")
-        return
-    _f = os.path.join(pathlib.Path(__file__).parent, 
-                            "att-data", 
-                            f"foxsi4_telescope-{telescope}_BASIC_collimator_aperture_ratio_V25APR13.fits") if file is None else file
-    with fits.open(_f) as hdul:
-        oa_angles, aperture_ratio = hdul[2].data << u.arcsec, hdul[1].data << u.dimensionless_unscaled
-    return np.interp(off_axis_angle.value, oa_angles.value, aperture_ratio.value, left=0, right=0) << u.dimensionless_unscaled
-
-if __name__=="__main__":
-    import matplotlib.pyplot as plt
-
-    from phot_spec import create_energy_midpoints, zeroes2nans
-
-    SAVE_ASSETS = False
-    assets_dir = os.path.join(pathlib.Path(__file__).parent, "..", "..", "assets", "response-tools-py-figs", "att-figs")
-    pathlib.Path(assets_dir).mkdir(parents=True, exist_ok=True)
-
-    mid_energies = create_energy_midpoints()
-
-    tb_col, obf0_col, obf1_col, cdte_fixed2_col, cdte_fixed4_col = plt.cm.viridis([0, 0.2, 0.4, 0.6, 0.8])
-    pix_att_meas_col, pix_att_mod_col, al_mylar_mod_col, cmost0_col, cmost1_col = plt.cm.plasma([0, 0.2, 0.4, 0.6, 0.8])
-    cmost2_col, cmost3_col = plt.cm.cividis([0.1, 0.9])
-
-    # all attenuators
-    plt.figure(figsize=(10,8))
-    att_therm_bl = zeroes2nans(att_thermal_blanket(mid_energies))
-    plt.ylabel(f"Transmission [{att_therm_bl.unit:latex}]")
-    plt.xlabel(f"Energy [{mid_energies.unit:latex}]")
-    p1 = plt.plot(mid_energies, att_therm_bl, color=tb_col, ls="-", label="Thermal blanket")
-
-    obf0 = zeroes2nans(att_obf0(mid_energies))
-    p2 = plt.plot(mid_energies, obf0, color=obf0_col, ls="--", label="OBF0", lw=3)
-
-    obf1 = zeroes2nans(att_obf1(mid_energies))
-    p3 = plt.plot(mid_energies, obf1, color=obf1_col, ls="-.", label="OBF1", lw=3)
-
-    cdte_fixed4 = zeroes2nans(att_uniform_al_cdte(mid_energies, position=4))
-    p4 = plt.plot(mid_energies, cdte_fixed4, color=cdte_fixed4_col, ls="-.", label="CdTe fixed p4")
-
-    cdte_fixed2 = zeroes2nans(att_uniform_al_cdte(mid_energies, position=2))
-    p5 = plt.plot(mid_energies, cdte_fixed2, color=cdte_fixed2_col, ls="-.", label="CdTe fixed p2")
-
-    pix_att_meas = zeroes2nans(att_pixelated(mid_energies, values="measured"))
-    p6 = plt.plot(mid_energies, pix_att_meas, color=pix_att_meas_col, ls="--", label="Pix. Att. Measured")
-
-    pix_att_mod = zeroes2nans(att_pixelated(mid_energies, values="modelled"))
-    p7 = plt.plot(mid_energies, pix_att_mod, color=pix_att_mod_col, ls="--", label="Pix. Att. Modelled")
-
-    al_mylar_mod = zeroes2nans(att_al_mylar(mid_energies))
-    p8 = plt.plot(mid_energies, al_mylar_mod, color=al_mylar_mod_col, ls="-", label="Al-Mylar")
-
-    cmos_t0 = att_cmos_filter(mid_energies, telescope=0)
-    p9 = plt.plot(mid_energies, cmos_t0, color=cmost0_col, ls="-", label="CMOS0-filter", lw=2)
-
-    cmos_t1 = att_cmos_filter(mid_energies, telescope=1)
-    p10 = plt.plot(mid_energies, cmos_t1, color=cmost1_col, ls="-.", label="CMOS1-filter", lw=2)
-
-    cmos_t2 = att_cmos_obfilter(mid_energies, telescope=0)
-    p11 = plt.plot(mid_energies, cmos_t2, color=cmost2_col, ls="--", label="CMOS0-OBF", lw=1)
-
-    cmos_t3 = att_cmos_obfilter(mid_energies, telescope=1)
-    p12 = plt.plot(mid_energies, cmos_t3, color=cmost3_col, ls=":", label="CMOS1-OBF", lw=1)
-
-    plt.ylim([0,1.05])
-    plt.xlim([0.5, 100])
-    plt.xscale("log")
-
-    plt.title("Attenuators")
-    
-    plt.legend(handles=p1+p2+p3+p4+p5+p6+p7+p8+p9+p10+p11+p12)
-    plt.tight_layout()
-    if SAVE_ASSETS:
-        plt.savefig(os.path.join(assets_dir,"transmissions.png"), dpi=200, bbox_inches="tight")
-    plt.show()
-
-    ## other
-    # collimator (so far, only have value for on-axis)
-    print(att_cmos_collimator_ratio(0<<u.arcsec, file=None, telescope=0))
-    print(att_cmos_collimator_ratio(0<<u.arcsec, file=None, telescope=1))
-
-    # any fake/functionaly model attenuators
-    # sigmoid
-    l ,x0, k, b = 1, 0, 1, 0
-    example = att_sigmoid(mid_energies, l ,x0, k, b)
-    p1 = plt.plot(mid_energies, example, ls="-", label=f"l:{l} ,x0:{x0}, k:{k}, b:{b}")
-    l ,x0, k, b = 1.04, 0, 1, 0
-    p2 = plt.plot(mid_energies, att_sigmoid(mid_energies, l, x0, k, b), ls="-", label=f"l:{l} ,x0:{x0}, k:{k}, b:{b}")
-    l ,x0, k, b = 1, 3, 1, 0
-    p3 = plt.plot(mid_energies, att_sigmoid(mid_energies, l, x0, k, b), ls="-", label=f"l:{l} ,x0:{x0}, k:{k}, b:{b}")
-    l ,x0, k, b = 1, 3, 1, 0.04
-    p4 = plt.plot(mid_energies, att_sigmoid(mid_energies, l, x0, k, b), ls="-", label=f"l:{l} ,x0:{x0}, k:{k}, b:{b}")
-    l ,x0, k, b = 1, 10, 0.5, 0
-    p5 = plt.plot(mid_energies, att_sigmoid(mid_energies, l, x0, k, b), ls="-", label=f"l:{l} ,x0:{x0}, k:{k}, b:{b}")
-    l ,x0, k, b = 1, 10, 1, 0
-    p6 = plt.plot(mid_energies, att_sigmoid(mid_energies, l, x0, k, b), ls="-", label=f"l:{l} ,x0:{x0}, k:{k}, b:{b}")
-    l ,x0, k, b = 1, 10, 2, 0
-    p7 = plt.plot(mid_energies, att_sigmoid(mid_energies, l, x0, k, b), ls="-", label=f"l:{l} ,x0:{x0}, k:{k}, b:{b}")
-    l ,x0, k, b = 1, 25, 1, 0
-    p8 = plt.plot(mid_energies, att_sigmoid(mid_energies, l, x0, k, b), ls="-", label=f"spectrometer-like:\nl:{l} ,x0:{x0}, k:{k}, b:{b}")
-
-    plt.ylabel(f"Transmission [{example.unit:latex}]")
-    plt.xlabel(f"Energy [{mid_energies.unit:latex}]")
-    plt.title("Fake Attenuator, Sigmoid$=l/(1+exp^{-k*(x-x0)})+b")
-    plt.ylim([0,1.05])
-    plt.xlim([0.5, 100])
-    plt.xscale("log")
-
-    plt.legend(handles=p1+p2+p3+p4+p5+p6+p7+p8)
-    plt.tight_layout()
-    if SAVE_ASSETS:
-        plt.savefig(os.path.join(assets_dir,"model-transmissions.png"), dpi=200, bbox_inches="tight")
-    plt.show()
\ No newline at end of file
diff --git a/response-tools-py/build/lib/response_tools_py/det_resp.py b/response-tools-py/build/lib/response_tools_py/det_resp.py
deleted file mode 100644
index a12b49d..0000000
--- a/response-tools-py/build/lib/response_tools_py/det_resp.py
+++ /dev/null
@@ -1,310 +0,0 @@
-"""Code to load different detectro responses. """
-
-import logging
-import os
-import pathlib
-
-from astropy.io import fits
-import astropy.units as u
-from matplotlib.colors import LogNorm, Normalize
-import matplotlib.gridspec as gridspec
-import matplotlib.pyplot as plt
-import numpy as np
-
-def cdte_det_resp_rmf(file):
-    """Return the redistribution matrix from a given file."""
-    e_lo, e_hi, ngrp, fchan, nchan, matrix = _read_rmf(file)
-
-    fchan_array = col2arr_py(fchan)
-    nchan_array = col2arr_py(nchan)
-
-    return e_lo<<u.keV, e_hi<<u.keV, vrmf2arr_py(data=matrix,
-                                                 n_grp_list=ngrp,
-                                                 f_chan_array=fchan_array,
-                                                 n_chan_array=nchan_array)<<(u.ct/u.ph)
-
-def cmos_det_resp(file=None, telescope=None):
-
-    if telescope is None:
-        logging.warning("`telescope` input in cmos_det_resp()` must be 0 or 1.")
-        return
-        
-    _f = os.path.join(pathlib.Path(__file__).parent, 
-                            "det-resp-data", 
-                            "cmos", 
-                            f"foxsi4_telescope-{telescope}_BASIC_RESPONSE_MATRIX_V25APR13.fits") if file is None else file
-    
-    with fits.open(_f) as hdul:
-        matrix, counts, energy = hdul[1].data<<(u.ct/u.ph), hdul[2].data<<u.dimensionless_unscaled, hdul[3].data<<u.keV # units?
-    return counts, energy, matrix 
-
-def _read_rmf(file):
-    """
-    Read a .rmf file and extract useful information from it.
-
-    Parameters
-    ----------
-    file :  `str`, `file-like` or `pathlib.Path`
-        A .rmf file (see `~astropy.fits.io.open` for details).
-
-    Returns
-    -------
-    `tuple`
-        The low and high boundary of energy bins (data['energ_lo'], 
-        data['energ_hi']), number of sub-set channels in the energybin 
-        (data['n_grp']), starting index of each sub-set of channels 
-        (data['f_chan']),number of channels in each sub-set 
-        (data['n_chan']), redistribution matrix [counts per photon] 
-        (data['matrix']).
-    """
-
-    with fits.open(file) as hdul:
-        data = hdul[2].data
-
-    return data["energ_lo"], data["energ_hi"], data["n_grp"], data["f_chan"], data["n_chan"], data["matrix"]
-
-def col2arr_py(data):
-    """Takes a list of parameters for each energy channel from a .rmf 
-    file and returns it in the correct format.
-
-    From: https://lost-contact.mit.edu/afs/physics.wisc.edu/home/craigm/lib/idl/util/vcol2arr.pro
-
-    **This original function is taken from Sunkit-spex**
-    - https://github.com/sunpy/sunkit-spex
-    - Tweaked to change `r` below to `[r]` since the values in data 
-    aren't lists (e.g., NuSTAR)
-
-    Parameters
-    ----------
-    data : array/list-like object
-            One parameter's array/list from the .rmf file.
-
-    Returns
-    -------
-    A 2D numpy array of the correctly ordered input data.
-
-    Example
-    -------
-    data = FITS_rec([(  1.6 ,   1.64,   1, [0]   , [18]  , [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]),
-                     (  1.64,   1.68,   1, [0]   , [20]  , [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]),
-                     (  1.68,   1.72,   2, [0,22], [20,1], [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]),
-                     dtype=(numpy.record, [('ENERG_LO', '>f4'), ('ENERG_HI', '>f4'), ('N_GRP', '>i2'),
-                                           ('F_CHAN', '>i4', (2,)), ('N_CHAN', '>i4', (2,)), ('MATRIX', '>i4', (2,))]))
-
-     col2arr_py(data['F_CHAN'])
-    array([[  0.,   0.],
-           [  0.,   0.],
-           [  0.,  22.]])
-    ## max row length of 2 so 2 columns, each row is an energy channel.
-    """
-
-    # this is the quicker way I have chosen to do in Python (this may be revised later but is ~30x faster than way below in Python)
-    max_len = np.max([len([r]) for r in data])  # find max row length
-    chan_array_py = np.array(
-        [[*[r], *(max_len - len([r])) * [0]] for r in data]
-    )  # make each row that length (padding with 0)
-
-    return chan_array_py
-
-def vrmf2arr_py(data=None, n_grp_list=None, f_chan_array=None, n_chan_array=None):
-    """Takes redistribution parameters for each energy channel from a 
-    .rmf file and returns it in the correct format.
-
-    From: https://lost-contact.mit.edu/afs/physics.wisc.edu/home/craigm/lib/idl/spectral/vrmf2arr.pro
-
-    **This original function is taken from Sunkit-spex**
-    - https://github.com/sunpy/sunkit-spex
-
-    Parameters
-    ----------
-    data : array/list-like object
-            Redistribution matrix parameter array/list from the .rmf file. Units are counts per photon.
-            Default : None
-
-    no_of_channels : int
-            Number of entries is the total number of photon channels, the entries themselves show the total number
-            of count channels to which that photon channel contributes.
-            Default : None
-
-    f_chan_array : numpy.array
-            The index of each sub-set channel from each energy bin from the .rmf file run through col2arr_py().
-            Default : None
-
-    n_chan_array : numpy.array
-            The number of sub-set channels in each index for each energy bin from the .rmf file run through col2arr_py().
-            Default : None
-
-    Returns
-    -------
-    A 2D numpy array of the correctly ordered input data with dimensions of energy in the rows and channels in
-    the columns.
-
-    Example
-    -------
-     d_rmf = 'directory/'
-     f_rmf = 'file.rmf'
-     e_lo, e_hi, ngrp, fchan, nchan, matrix = det_io._read_rmf(d_rmf+f_rmf)
-
-     fchan_array = nu_spec.col2arr_py(fchan)
-     nchan_array = nu_spec.col2arr_py(nchan)
-
-     rmf = nu_spec.vrmf2arr_py(data=matrix,
-                                  n_grp_list=ngrp,
-                                  f_chan_array=fchan_array,
-                                  n_chan_array=nchan_array)
-     rmf
-
-    array([[0.00033627, 0.0007369 , 0.00113175, ..., 0.        , 0.        , 0.        ],
-           [0.00039195, 0.00079259, 0.00138341, ..., 0.        , 0.        , 0.        ],
-           [0.00042811, 0.00083381, 0.00157794, ..., 0.        , 0.        , 0.        ],
-                                                ...,
-           [0.        , 0.        , 0.        , ..., 0.00408081, 0.00409889, 0.00403308],
-           [0.        , 0.        , 0.        , ..., 0.00405333, 0.00413722, 0.00413216],
-           [0.        , 0.        , 0.        , ..., 0.        , 0.        , 0.        ]])
-    ## rows = photon/energy channels, columns = counts channels
-
-    What's Going On?
-    ----------------
-    The RMF file has the photon-to-counts conversion information in it.
-    The martix has the photon-to-count conversion value for each count channel (columns) that is involved with theach photon channel (rows).
-            E.g., matrix = [ [a, b, c, d, e, f, ...] ,
-                             [        ...          ] ,
-                             [        ...          ] ,
-                                      ...             ]
-    F_chan is the starting index of contiguous counts channels that are involved with the photon channel.
-            E.g., f_chan = [ [0, 5, 0, 0, 0, ...] ,
-                             [       ...        ] ,
-                             [       ...        ] ,
-                                     ...           ]
-                            For the first photon channel, there are rows of counts channels starting at index 0 and 5
-    N_chan is the corresponding number of counts channels from each index in the f_chan array.
-            E.g., n_chan = [ [2, 3, 0, 0, 0, ...] ,
-                             [        ...        ] ,
-                             [        ...        ] ,
-                                      ...           ]
-                            Starting at index 0 for the first photon channel we have the first 2 matrix values, then at index 5 we have the next 3.
-                            The total of each row is the same as the n_grp_list and the number of entries in each row of the matrix entry.
-    Putting all this together, the rmf matrix is:
-            rmf_matrix = [ [a, b, 0, 0, 0, c , d , e, 0 , 0 , ...] ,   #<-- index 0 (f_chan) with 2 entries (n_chan) with photon-to-counts conversion (matrix)
-                         [                 ...                   ] ,
-                         [                 ...                   ] ,
-                                           ...                      ]
-    """
-    # this was is about >6x quicker in than the IDL code written in Python
-
-    n_grp_list = n_grp_list.astype("<i2")  # change from ‘big-endian’ (">i2") to ‘little-endian’ ("<i2")
-
-    # find the non-zero entries in Nchan, this is the number to counts channels
-    #  in a row that contribute so will have a value if it is useful
-    b = np.nonzero(n_chan_array)
-
-    # now only want the useful entries from the pre-formatted Nchan and Fchan arrays
-    c = f_chan_array[b]
-    d = n_chan_array[b]
-
-    # to help with indexing, this provides a running sum of the number of counts
-    #  channels that a single photon channel contributes to
-    e = np.cumsum(n_chan_array, axis=1)
-
-    # these entries will give the final indices in the row on counts channels
-    final_inds = e[b]
-
-    # need to find the starting index so -1, but that means any entry that is
-    #  -1 will be where a zero is needed
-    starting_inds = b[1] - 1
-
-    # get the  starting indices but the ones that should be 0 are replaced with
-    #  the final one in the list at the minute (-1 in starting_inds)
-    start_inds = np.cumsum(n_chan_array, axis=1)[(b[0], starting_inds)]
-
-    # where starting_inds==-1 that value should be 0, i.e. starting from the first
-    #  value in the rmf matrix
-    new_e = np.where(starting_inds != -1, start_inds, 0)
-
-    # initialise the rmf matrix
-    mat_array_py = np.zeros((len(data), len(n_grp_list)))
-
-    # now go through row by row (this is the slowest part and needs to be made faster).
-    #  Here we go through each photon channel's number of discrete rows of counts channels.
-    for r in range(len(c)):
-        mat_array_py[b[0][r], c[r] : c[r] + d[r]] = data[b[0][r]][new_e[r] : final_inds[r]]
-
-    return mat_array_py
-
-
-if __name__=="__main__":
-
-    SAVE_ASSETS = False
-    assets_dir = os.path.join(pathlib.Path(__file__).parent, "..", "..", "assets", "response-tools-py-figs", "det-resp-figs")
-    pathlib.Path(assets_dir).mkdir(parents=True, exist_ok=True)
-
-    # CMOS
-    fig = plt.figure(figsize=(12,7))
-
-    gs = gridspec.GridSpec(1, 2)
-    gs_ax0 = fig.add_subplot(gs[0, 0])
-    telescope = 0
-    c, e, m = cmos_det_resp(file=None, telescope=telescope)
-    extent = [np.min(c.value), np.max(c.value), np.min(e.value), np.max(e.value)]
-    i = gs_ax0.imshow(m.value, origin="lower", aspect=(extent[1]-extent[0])/(extent[3]-extent[2]), extent=extent, norm=LogNorm())
-    gs_ax0.set_ylim([0,10])
-    gs_ax0.set_xlabel("Counts?")
-    gs_ax0.set_ylabel(f"Energy? [{e.unit:latex}]")
-    gs_ax0.set_title(f"CMOS{telescope}")
-    cax = gs_ax0.inset_axes([0.1, 0.08, 0.8, 0.05],)
-    cbar = fig.colorbar(i, cax=cax, orientation='horizontal')
-    cbar.set_label(f"{m.unit:latex}?", size=8, labelpad=0)
-    cax.tick_params(axis='both', which='major', labelsize=6)
-
-
-    gs_ax1 = fig.add_subplot(gs[0, 1])
-    telescope = 1
-    c, e, m = cmos_det_resp(file=None, telescope=telescope)
-    extent = [np.min(c.value), np.max(c.value), np.min(e.value), np.max(e.value)]
-    i = gs_ax1.imshow(m.value, origin="lower", aspect=(extent[1]-extent[0])/(extent[3]-extent[2]), extent=extent, norm=LogNorm())
-    gs_ax1.set_ylim([0,10])
-    gs_ax1.set_xlabel("Counts?")
-    gs_ax1.set_ylabel(f"Energy? [{e.unit:latex}]")
-    gs_ax1.set_title(f"CMOS{telescope}")
-    cax = gs_ax1.inset_axes([0.1, 0.08, 0.8, 0.05],)
-    cbar = fig.colorbar(i, cax=cax, orientation='horizontal')
-    cbar.set_label(f"{m.unit:latex}?", size=8, labelpad=0)
-    cax.tick_params(axis='both', which='major', labelsize=6)
-
-    plt.tight_layout()
-    if SAVE_ASSETS:
-        plt.savefig(os.path.join(assets_dir,"cmos-response-matrices.png"), dpi=200, bbox_inches="tight")
-    plt.show()
-
-    # CdTe
-    d_rmf = "./det-resp-data/cdte/merged/"
-    d_rmf = "./det-resp-data/cdte/pt/"
-    f_rmf = "Resp_3keVto30keV_CdTe1_reg0_1hit.rmf"
-
-
-    e_lo, e_hi, rmf = cdte_det_resp_rmf(os.path.join(pathlib.Path(__file__).parent, d_rmf, f_rmf))
-
-    fig = plt.figure(figsize=(12, 5))
-    gs = gridspec.GridSpec(1, 2)
-
-    gs_ax0 = fig.add_subplot(gs[0, 0])
-    r = gs_ax0.imshow(rmf.value, origin="lower", norm=Normalize(vmin=0.001, vmax=0.12), extent=[np.min(e_lo.value), np.max(e_lo.value), np.min(e_lo.value), np.max(e_lo.value)])
-    cbar = plt.colorbar(r)
-    cbar.ax.set_ylabel('Counts photon$^{-1}$')
-    fig.suptitle(d_rmf+f_rmf)
-    gs_ax0.set_xlabel("Count Energy [keV]")
-    gs_ax0.set_ylabel("Photon Energy [keV]")
-    gs_ax0.set_title("Linear Scale")
-
-    gs_ax0 = fig.add_subplot(gs[0, 1])
-    r = gs_ax0.imshow(rmf.value, origin="lower", norm=LogNorm(vmin=0.001, vmax=0.12), extent=[np.min(e_lo.value), np.max(e_lo.value), np.min(e_lo.value), np.max(e_lo.value)])
-    cbar = plt.colorbar(r)
-    cbar.ax.set_ylabel('Counts photon$^{-1}$')
-    fig.suptitle(d_rmf+f_rmf)
-    gs_ax0.set_xlabel("Count Energy [keV]")
-    gs_ax0.set_ylabel("Photon Energy [keV]")
-    gs_ax0.set_title("Log Scale")
-
-    if SAVE_ASSETS:
-        plt.savefig(os.path.join(assets_dir,"cdte-response-matrix.png"), dpi=200, bbox_inches="tight")
-    plt.show()
\ No newline at end of file
diff --git a/response-tools-py/build/lib/response_tools_py/eff_area.py b/response-tools-py/build/lib/response_tools_py/eff_area.py
deleted file mode 100644
index 951020e..0000000
--- a/response-tools-py/build/lib/response_tools_py/eff_area.py
+++ /dev/null
@@ -1,320 +0,0 @@
-"""Code to load different effective areas. """
-
-import logging
-import os
-import pathlib
-
-from astropy.io import fits
-import astropy.units as u
-import numpy as np
-from scipy.interpolate import CloughTocher2DInterpolator, LinearNDInterpolator 
-import pandas
-
-@u.quantity_input(mid_energies=u.keV, off_axis=u.arcmin)
-def eff_area_msfc_10shell(mid_energies, off_axis=0, file_tilt=None, file_pan=None, optic_id=None):
-    """ 
-    Tilt is the off-axis angle when the optics up coincides with Solar North (“vertical”). Pan is the other angle (“horizontal”).
-    """
-    _id = optic_id if optic_id in ["X-7", "X-8"] else None
-    if _id is None:
-        logging.warning("Please provide a MSFC heritage optic ID from [\'X-7\', \'X-8\'].")
-        
-    vals_tilt = _get_ea_file_info(_id, "tilt", file=file_tilt)
-    vals_pan = _get_ea_file_info(_id, "pan", file=file_pan)
-
-    # from Milo on 17/3/2025 @ 13:47 (Slack)
-    ea_energies = [4.5,  5.5,  6.5,  7.5,  8.5,  9.5, 11. , 13. , 15. , 17. , 19. , 22.5, 27.5] << u.keV
-
-    off_axis_angles_tilt, eff_areas_tilt = _get_oa_and_ea_msfc_10shell(vals_tilt)
-    _, eff_areas_pan = _get_oa_and_ea_msfc_10shell(vals_pan)
-
-    mean_eff_areas_areas = (eff_areas_tilt+eff_areas_pan)/2
-
-    x, y = np.meshgrid(ea_energies, off_axis_angles_tilt)
-    i = CloughTocher2DInterpolator(list(zip(x.flatten().value, y.flatten().value)), mean_eff_areas_areas.T.flatten().value)
-
-    return i(mid_energies, off_axis) << u.cm**2
-
-def _get_ea_file_info(optics_id, axis, file=None):
-    _f = os.path.join(pathlib.Path(__file__).parent, 
-                            "eff-area-data", 
-                            f"FOXSI3_Module_{optics_id}_EA_{axis}.txt") if file is None else file
-    return np.loadtxt(_f, delimiter=",")
-
-def _get_oa_and_ea_msfc_10shell(grid):
-    off_axis_angles = grid[0]
-    eff_areas = grid[1:]
-    return off_axis_angles << u.arcmin, eff_areas << u.cm**2
-
-@u.quantity_input(mid_energies=u.keV)
-def eff_area_msfc_hi_res(mid_energies, file=None, position=None, use_model=False):
-    """Return MSCF hi-res effective areas interpolated to the given energies.
-    
-    Latest from Wayne.
-    - older function is `eff_area_msfc`
-    
-    Position 0: X10/FM2
-    Position 3: X09/FM1
-    Position 6: X11/FM3
-    """
-    # msfc_hi_res effective areas
-    if use_model:
-        _f = os.path.join(pathlib.Path(__file__).parent, 
-                                "eff-area-data", 
-                                "FOXSI4_Module_MSFC_HiRes_EA_v1.txt") if file is None else file
-        e, fm1, fm2, fm3 = np.loadtxt(_f).T
-    else:
-        _f = os.path.join(pathlib.Path(__file__).parent, 
-                                "eff-area-data", 
-                                "FOXSI4_Module_MSFC_HiRes_EA_with models_v1.txt") if file is None else file
-        e, _, _, _, fm1, fm2, fm3 = np.loadtxt(_f).T
-
-    e <<= u.keV
-    fm1 <<= u.cm**2
-    fm2 <<= u.cm**2
-    fm3 <<= u.cm**2
-    if position==0:
-        return np.interp(mid_energies.value, e.value, fm2.value, left=fm2.value[0], right=0) << u.cm**2
-    elif position==3:
-        return np.interp(mid_energies.value, e.value, fm1.value, left=fm1.value[0], right=0) << u.cm**2
-    elif position==6:
-        return np.interp(mid_energies.value, e.value, fm3.value, left=fm3.value[0], right=0) << u.cm**2
-    else:
-        logging.warning("`position` must be 0 (X10/FM2), 3 (X09/FM1), or 6 (X11/FM3).")
-
-@u.quantity_input(mid_energies=u.keV)
-def eff_area_msfc(mid_energies, file=None):
-    """Return MSCF hi-res effective areas interpolated to the given energies."""
-    # msfc_hi_res effective areas
-    logging.warning("This might not be the function you are looking for, please see `eff_area_msfc_hi_res`.")
-    _f = os.path.join(pathlib.Path(__file__).parent, 
-                            "eff-area-data", 
-                            "3Inner_EA_EPDL97_14AA.csv") if file is None else file
-    msfc_hi_res = pandas.read_csv(_f).to_numpy()[:,1:] # remove the first column that only indexes
-    # in cm2 ; we use the innermost and the 3rd innermost shells (S10 and S08) [from Yixian]
-    msfc_hi_res_es, msfc_hi_res_effas08, msfc_hi_res_effas10 = msfc_hi_res[:,0] << u.keV, msfc_hi_res[:,1] << u.cm**2, msfc_hi_res[:,3] << u.cm**2 
-    msfc_hi_res_effa = msfc_hi_res_effas08 + msfc_hi_res_effas10
-    return np.interp(mid_energies.value, msfc_hi_res_es.value, msfc_hi_res_effa.value, left=msfc_hi_res_effa.value[0], right=0) << u.cm**2
-
-@u.quantity_input(mid_energies=u.keV)
-def eff_area_nagoya(mid_energies, file=None):
-    """Return Nagoya SXR hi-res effective areas interpolated to the given energies."""
-    # nagoya sxr effective areas
-    _f = os.path.join(pathlib.Path(__file__).parent, 
-                            "eff-area-data", 
-                            "effective-area_raytracing_soft-xray-optic_on-axis.txt") if file is None else file
-    nagoya_sxr = np.loadtxt(_f)
-    nagoya_sxr_es, nagoya_sxr_effa = nagoya_sxr[:,0] << u.keV, nagoya_sxr[:,1]/100 << u.cm**2
-    return np.interp(mid_energies.value, nagoya_sxr_es.value, nagoya_sxr_effa.value, left=0, right=0) << u.cm**2
-
-@u.quantity_input(mid_energies=u.keV)
-def eff_area_cmos(mid_energies, file=None, telescope=None):
-    """Return optics paired with CMOS effective areas interpolated to 
-    the given energies.
-    
-    Telescope 0: position 0, X10/FM2(?)
-    Telescope 1: position 1, Nagoya(?)
-    """
-    if telescope is None:
-        logging.warning("`telescope` input in `eff_area_cmos()` must be 0 or 1.")
-        return
-        
-    _f = os.path.join(pathlib.Path(__file__).parent, 
-                            "eff-area-data", 
-                            f"foxsi4_telescope-{telescope}_BASIC_mirror_effective_area_V25APR13.fits") if file is None else file
-    with fits.open(_f) as hdul:
-        es, effas = hdul[2].data << u.keV, hdul[1].data << u.cm**2
-    return np.interp(mid_energies.value, es.value, effas.value, left=0, right=0) << u.cm**2
-
-@u.quantity_input(mid_energies=u.keV)
-def eff_area_cmos_telescope(mid_energies, file=None, telescope=None):
-    """Return full telescope(?) with CMOS effective areas interpolated to 
-    the given energies.
-    
-    Telescope 0: position 0, X10/FM2(?)
-    Telescope 1: position 1, Nagoya(?)
-
-    **Note**
-    Will come back to this when there are off-axis angles and do 2D 
-    interp. Currently only on-axis angle so just interp across energies.
-    """
-    if telescope is None:
-        logging.warning("`telescope` input in `eff_area_cmos_telescope()` must be 0 or 1.")
-        return
-        
-    _f = os.path.join(pathlib.Path(__file__).parent, 
-                            "eff-area-data", 
-                            f"foxsi4_telescope-{telescope}_BASIC_TELESCOPE_RESPONSE_V25APR13.fits") if file is None else file
-    with fits.open(_f) as hdul:
-        # _ is the off-axis angle but it's just [0] at the minute
-        ea_energies, _, effas = hdul[2].data << u.keV, hdul[3].data << u.arcsec, hdul[1].data << u.cm**2
-
-    return np.interp(mid_energies.value, ea_energies.value, effas.value, left=0, right=0) << u.cm**2
-
-if __name__=="__main__":
-    import matplotlib.pyplot as plt
-    import matplotlib.gridspec as gridspec
-
-    SAVE_ASSETS = False
-    assets_dir = os.path.join(pathlib.Path(__file__).parent, "..", "..", "assets", "response-tools-py-figs", "eff-area-figs")
-    pathlib.Path(assets_dir).mkdir(parents=True, exist_ok=True)
-
-    mid_energies = np.linspace(0, 20, 1000)<<u.keV
-    
-    #cmos optics/telescopes
-    fig = plt.figure(figsize=(14,7))
-    gs = gridspec.GridSpec(1, 2)
-
-
-    gs_ax0 = fig.add_subplot(gs[0, 0])
-    a0 = eff_area_cmos(mid_energies, file=None, telescope=0)
-    a1 = eff_area_cmos(mid_energies, file=None, telescope=1)
-    msfcp0 = eff_area_msfc_hi_res(mid_energies, file=None, position=0)
-    msfcp0m = eff_area_msfc_hi_res(mid_energies, file=None, position=0, use_model=True)
-    gs_ax0.plot(mid_energies, a0, label="CMOS telescope 0, position 0")
-    gs_ax0.plot(mid_energies, a1, label="CMOS telescope 1, position 1")
-    gs_ax0.plot(mid_energies, msfcp0, label="MSFC (meas.) position 0")
-    gs_ax0.plot(mid_energies, msfcp0m, label="MSFC (mod.) position 0")
-    gs_ax0.set_title("CMOS Optics")
-    gs_ax0.set_ylabel(f"Effective Area [{a0.unit:latex}]")
-    gs_ax0.set_xlabel(f"Energy [{mid_energies.unit:latex}]")
-    plt.legend()
-
-    # CMOS full telescope ones
-    gs_ax1 = fig.add_subplot(gs[0, 1])
-    a0ft = eff_area_cmos_telescope(mid_energies, file=None, telescope=0)
-    a1ft = eff_area_cmos_telescope(mid_energies, file=None, telescope=1)
-    gs_ax1.plot(mid_energies, a0ft, label="CMOS telescope 0, position 0")
-    gs_ax1.plot(mid_energies, a1ft, label="CMOS telescope 1, position 1")
-    gs_ax1.set_title("CMOS Telescope?")
-    gs_ax1.set_ylabel(f"Effective Area [{a0ft.unit:latex}]")
-    gs_ax1.set_xlabel(f"Energy [{mid_energies.unit:latex}]")
-    plt.legend()
-
-    plt.tight_layout()
-    if SAVE_ASSETS:
-        plt.savefig(os.path.join(assets_dir,"cmos-optics-resp.png"), dpi=200, bbox_inches="tight")
-    plt.show()
-
-    # most optics
-    fig = plt.figure(figsize=(9.5, 9.5))
-    gs = gridspec.GridSpec(3, 3)
-
-    ## OPTIC: X-7
-    optic = "X-7" 
-
-    vals_tilt = _get_ea_file_info(optic, "tilt")
-    vals_pan = _get_ea_file_info(optic, "pan")
-    off_axis_angles_tilt, eff_areas_tilt = _get_oa_and_ea_msfc_10shell(vals_tilt)
-    off_axis_angles_pan, eff_areas_pan = _get_oa_and_ea_msfc_10shell(vals_pan)
-    ea_energies = [4.5,  5.5,  6.5,  7.5,  8.5,  9.5, 11. , 13. , 15. , 17. , 19. , 22.5, 27.5] << u.keV
-    ea_energies_binned = ea_energies
-
-    off_axis_angles_tilt = off_axis_angles_tilt << u.arcmin
-    off_axis_angles_pan = off_axis_angles_pan << u.arcmin
-
-    e_i = 0
-    gs_ax0 = fig.add_subplot(gs[0, 0])
-    p00 = gs_ax0.plot(off_axis_angles_tilt, eff_areas_tilt[e_i], label="tilt", lw=4)
-    p01 = gs_ax0.plot(off_axis_angles_pan, eff_areas_pan[e_i], label="pan")
-    gs_ax0.set_title(f"@ Energy: {ea_energies[e_i]}")
-    gs_ax0.set_ylabel(f"Heritage MSFC {optic} [{eff_areas_pan.unit:latex}]")
-    gs_ax0.set_xlabel(f"off-axis angle [{off_axis_angles_tilt.unit:latex}]")
-    plt.legend(handles=p00+p01)
-    gs_ax0.set_xlim([off_axis_angles_tilt[0].value, off_axis_angles_tilt[-1].value])
-
-    oa_i = 5
-    gs_ax1 = fig.add_subplot(gs[0, 1])
-    p10 = gs_ax1.plot(ea_energies, eff_areas_tilt[:,oa_i], label="tilt", lw=4)
-    p11 = gs_ax1.plot(ea_energies, eff_areas_pan[:,oa_i], label="pan", lw=3)
-    p12 = gs_ax1.plot(ea_energies, (eff_areas_pan[:,oa_i]+eff_areas_tilt[:,oa_i])/2, label="ave")
-    gs_ax1.set_title(f"@ Tilt:{off_axis_angles_tilt[oa_i]:latex}, Pan:{off_axis_angles_pan[oa_i]:latex}")
-    gs_ax1.set_ylabel(f"{optic} [{eff_areas_pan.unit:latex}]")
-    gs_ax1.set_xlabel(f"Energy [{ea_energies.unit:latex}]")
-    plt.legend(handles=p10+p11+p12)
-    gs_ax1.set_xlim([ea_energies[0].value, ea_energies[-1].value])
-    # gs_ax1.set_ylim([0, np.nanmax(eff_areas_tilt[:,oa_i]).value*1.01])
-
-    gs_ax2 = fig.add_subplot(gs[0, 2])
-    _ps = []
-    for oaa in off_axis_angles_tilt:
-        _lw = 1 if oaa >=0 else 2
-        _p = gs_ax2.plot(ea_energies, eff_area_msfc_10shell(ea_energies, off_axis=oaa<<u.arcmin, optic_id=optic), label=f"{oaa:latex}", lw=_lw)
-        _ps += _p
-    gs_ax2.set_ylabel(f"{optic} [{eff_areas_pan.unit:latex}]")
-    gs_ax2.set_xlabel(f"Energy [{ea_energies.unit:latex}]")
-    plt.legend(handles=_ps, prop={"size": 9})
-    gs_ax2.set_xlim([ea_energies[0].value, ea_energies[-1].value])
-
-    for gsax in [gs_ax0, gs_ax1, gs_ax2]:
-        gsax.set_ylim([0, np.nanmax(eff_area_msfc_10shell(ea_energies, off_axis=0<<u.arcmin, optic_id=optic)).value*1.01])
-
-
-    ## OPTIC: X-8
-    optic = "X-8"
-
-    vals_tilt = _get_ea_file_info(optic, "tilt")
-    vals_pan = _get_ea_file_info(optic, "pan")
-    off_axis_angles_tilt, eff_areas_tilt = _get_oa_and_ea_msfc_10shell(vals_tilt)
-    off_axis_angles_pan, eff_areas_pan = _get_oa_and_ea_msfc_10shell(vals_pan)
-    ea_energies = [4.5,  5.5,  6.5,  7.5,  8.5,  9.5, 11. , 13. , 15. , 17. , 19. , 22.5, 27.5] << u.keV
-    ea_energies_binned = ea_energies
-    e_i = 0
-    gs_ax3 = fig.add_subplot(gs[1, 0])
-    p00 = gs_ax3.plot(off_axis_angles_tilt, eff_areas_tilt[e_i], label="tilt", lw=4)
-    p01 = gs_ax3.plot(off_axis_angles_pan, eff_areas_pan[e_i], label="pan")
-    gs_ax3.set_title(f"@ Energy: {ea_energies[e_i]}")
-    gs_ax3.set_ylabel(f"Heritage MSFC {optic} [{eff_areas_pan.unit:latex}]")
-    gs_ax3.set_xlabel(f"off-axis angle [{off_axis_angles_tilt.unit:latex}]")
-    plt.legend(handles=p00+p01)
-    gs_ax3.set_xlim([off_axis_angles_tilt[0].value, off_axis_angles_tilt[-1].value])
-
-    oa_i = 5
-    gs_ax4 = fig.add_subplot(gs[1, 1])
-    p10 = gs_ax4.plot(ea_energies, eff_areas_tilt[:,oa_i], label="tilt", lw=4)
-    p11 = gs_ax4.plot(ea_energies, eff_areas_pan[:,oa_i], label="pan", lw=3)
-    p12 = gs_ax4.plot(ea_energies, (eff_areas_pan[:,oa_i]+eff_areas_tilt[:,oa_i])/2, label="ave")
-    gs_ax4.set_title(f"@ Tilt:{off_axis_angles_tilt[oa_i]:latex}, Pan:{off_axis_angles_pan[oa_i]:latex}")
-    gs_ax4.set_ylabel(f"{optic} [{eff_areas_pan.unit:latex}]")
-    gs_ax4.set_xlabel(f"Energy [{ea_energies.unit:latex}]")
-    plt.legend(handles=p10+p11+p12)
-    gs_ax4.set_xlim([ea_energies[0].value, ea_energies[-1].value])
-    # gs_ax4.set_ylim([0, np.nanmax(eff_areas_tilt[:,oa_i]).value*1.01])
-
-    gs_ax5 = fig.add_subplot(gs[1, 2])
-    _ps = []
-    for oaa in off_axis_angles_tilt:
-        _lw = 1 if oaa >=0 else 2
-        _p = gs_ax5.plot(ea_energies, eff_area_msfc_10shell(ea_energies, off_axis=oaa<<u.arcmin, optic_id=optic), label=f"{oaa:latex}", lw=_lw)
-        _ps += _p
-    gs_ax5.set_ylabel(f"{optic} [{eff_areas_pan.unit:latex}]")
-    gs_ax5.set_xlabel(f"Energy [{ea_energies.unit:latex}]")
-    plt.legend(handles=_ps, prop={"size": 9})
-    gs_ax5.set_xlim([ea_energies[0].value, ea_energies[-1].value])
-
-    for gsax in [gs_ax3, gs_ax4, gs_ax5]:
-        gsax.set_ylim([0, np.nanmax(eff_area_msfc_10shell(ea_energies, off_axis=0<<u.arcmin, optic_id=optic)).value*1.01])
-
-    gs_ax5 = fig.add_subplot(gs[2, 1])
-    _f = os.path.join(pathlib.Path(__file__).parent, 
-                            "eff-area-data", 
-                            "FOXSI4_Module_MSFC_HiRes_EA_v1.txt")
-    e, *_ = np.loadtxt(_f).T
-    e <<= u.keV
-    msfc_hi_res_p0 = eff_area_msfc_hi_res(e, position=0)
-    msfc_hi_res_p3 = eff_area_msfc_hi_res(e, position=3)
-    msfc_hi_res_p6 = eff_area_msfc_hi_res(e, position=6)
-    p1 = gs_ax5.plot(e, msfc_hi_res_p0, label="Pos. 0 (X10/FM2)")
-    p2 = gs_ax5.plot(e, msfc_hi_res_p3, label="Pos. 3 (X09/FM1)")
-    p3 = gs_ax5.plot(e, msfc_hi_res_p6, label="Pos. 6 (X11/FM3)")
-    gs_ax5.set_title("MSFC High-res.")
-    gs_ax5.set_xlabel(f"Energy {e.unit:latex}")
-    gs_ax5.set_ylabel(f"Eff. Area {msfc_hi_res_p0.unit:latex}")
-    plt.legend(handles=p1+p2+p3)
-
-    gs_ax5.set_xlim([e[0].value, e[-1].value])
-
-    plt.tight_layout()
-    if SAVE_ASSETS:
-        plt.savefig(os.path.join(assets_dir,"heritage-and-msfc-optics.png"), dpi=200, bbox_inches="tight")
-    plt.show()
\ No newline at end of file
diff --git a/response-tools-py/build/lib/response_tools_py/io/__init__.py b/response-tools-py/build/lib/response_tools_py/io/__init__.py
deleted file mode 100644
index e69de29..0000000
diff --git a/response-tools-py/build/lib/response_tools_py/io/load_yaml.py b/response-tools-py/build/lib/response_tools_py/io/load_yaml.py
deleted file mode 100644
index 12c9175..0000000
--- a/response-tools-py/build/lib/response_tools_py/io/load_yaml.py
+++ /dev/null
@@ -1,19 +0,0 @@
-"""
-Module containing methods to load in specific or general YAML files.
-"""
-
-import os
-import pathlib
-
-import yaml
-
-def load_yaml(filename):
-    """Function to load in a general YAML file. """
-    with open(filename, "r") as file:
-        file_info = yaml.safe_load(file)
-
-    return file_info
-
-def load_cdte_context():
-    """Function to load in the context information from file. """
-    return load_yaml(os.path.join(pathlib.Path(__file__).parent, "..", "..", "..", "response-information", "info.yaml"))
\ No newline at end of file
diff --git a/response-tools-py/build/lib/response_tools_py/io/tests/__init__.py b/response-tools-py/build/lib/response_tools_py/io/tests/__init__.py
deleted file mode 100644
index e69de29..0000000
diff --git a/response-tools-py/build/lib/response_tools_py/ir-plots.py b/response-tools-py/build/lib/response_tools_py/ir-plots.py
deleted file mode 100644
index 5d74fce..0000000
--- a/response-tools-py/build/lib/response_tools_py/ir-plots.py
+++ /dev/null
@@ -1,641 +0,0 @@
-
-import os
-import pathlib
-
-import astropy.units as u
-from matplotlib.colors import LogNorm
-import matplotlib.pyplot as plt
-import matplotlib.gridspec as gridspec
-import numpy as np
-
-from att import (att_thermal_blanket, 
-                    att_uniform_al_cdte, 
-                    att_pixelated, 
-                    att_al_mylar,
-                    att_cmos_obfilter,
-                    att_cmos_collimator_ratio,
-                    att_cmos_filter,)
-from det_resp import (cdte_det_resp_rmf, 
-                        cmos_det_resp,)
-from eff_area import (eff_area_msfc_10shell, 
-                        eff_area_msfc_hi_res,
-                        eff_area_cmos,)
-from phot_spec import zeroes2nans
-from quantum_eff import qe_cmos
-
-def foxsi4_response_paths(save_assets=False):
-    
-    assets_dir = os.path.join(pathlib.Path(__file__).parent, "..", "..", "assets", "response-tools-py-figs", "response-paths")
-    pathlib.Path(assets_dir).mkdir(parents=True, exist_ok=True)
-
-    plt.rcParams.update({'font.size': 7})
-
-    tb_col = plt.cm.viridis([.5])
-
-    FIGSIZE = (18,4)
-    GS = (1, 12)
-
-    ###########################################################################
-    ########################## CdTe 4 (pos. 2) ################################
-    ###########################################################################
-
-    fig = plt.figure(figsize=FIGSIZE)
-    gs = gridspec.GridSpec(*GS)
-
-    # define det response first so we know what energies we actually care about
-    f_rmf = os.path.join(pathlib.Path(__file__).parent, "det-resp-data", "cdte", "merged", "Resp_3keVto30keV_CdTe1_reg0_1hit.rmf")
-    e_lo, e_hi, rmf = cdte_det_resp_rmf(f_rmf)
-    mid_energies = (e_lo+e_hi)/2
-    photon_energy_lims = [np.nanmin(mid_energies.value), np.nanmax(mid_energies.value)]
-
-    ## thermal blanket
-    gs_ax0 = fig.add_subplot(gs[0, 1:3])
-    att_therm_bl = zeroes2nans(att_thermal_blanket(mid_energies))
-    gs_ax0.set_xlabel(f"Transmission [{att_therm_bl.unit:latex}]")
-    gs_ax0.set_ylabel(f"Photon Energy [{e_lo.unit:latex}]")
-    gs_ax0.set_title("Thermal blanket")
-    gs_ax0.plot(att_therm_bl, mid_energies, color=tb_col, ls="-")
-
-    y_lims = [0,1.01]
-    gs_ax0.set_xlim(y_lims)
-    gs_ax0.set_ylim(photon_energy_lims)
-    gs_ax0.set_aspect(np.diff(y_lims)/np.diff(photon_energy_lims))
-
-    ## Marshall 10-shell
-    optic = "X-7"
-    gs_ax1 = fig.add_subplot(gs[0, 3:5])
-    gs_ax1.set_title(f"Marshall 10-shell ({optic})")
-    gs_ax1.set_aspect('equal', 'box')
-    gs_ax1.set_yticks([])
-    _ps = []
-    for oaa in [-8.64, -6.72, -4.8, -2.88, -0.96, 0., 0.96, 2.88, 4.8, 6.72, 8.64]<<u.arcmin:
-        _lw = 1 if oaa >=0 else 2
-        _eff_areas = eff_area_msfc_10shell(mid_energies, off_axis=oaa<<u.arcmin, optic_id=optic)
-        _p = gs_ax1.plot(_eff_areas, mid_energies, label=f"{oaa:latex}", lw=_lw)
-        _ps += _p
-    gs_ax1.set_xlabel(f"{optic} [{_eff_areas.unit:latex}]")
-    plt.legend(handles=_ps, prop={"size": 6})
-    oa_eff_area = eff_area_msfc_10shell(mid_energies, off_axis=0<<u.arcmin, optic_id=optic)
-    y_lims = [0, np.nanmax(oa_eff_area).value*1.01]
-    gs_ax1.set_xlim(y_lims)
-    gs_ax1.set_ylim(photon_energy_lims)
-    gs_ax1.set_aspect(np.diff(y_lims)/np.diff(photon_energy_lims))
-    gs_ax1.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## Al filter
-    gs_ax2 = fig.add_subplot(gs[0, 5:7])
-    gs_ax2.set_title("Al (0.015\")")
-    cdte_fixed2 = zeroes2nans(att_uniform_al_cdte(mid_energies, file=None, position=2))
-    plt.plot(cdte_fixed2, mid_energies, label="CdTe fixed p2")
-    gs_ax2.set_xlabel(f"Transmission [{cdte_fixed2.unit:latex}]")
-    y_lims = [0,1.01]
-    gs_ax2.set_xlim(y_lims)
-    gs_ax2.set_ylim(photon_energy_lims)
-    gs_ax2.set_aspect(np.diff(y_lims)/np.diff(photon_energy_lims))
-    gs_ax2.set_yticks([])
-    gs_ax2.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## no2021_07
-    gs_ax3 = fig.add_subplot(gs[0, 7:9])
-    gs_ax3.imshow(rmf.value, origin="lower", norm=LogNorm(vmin=0.001, vmax=0.12), extent=[np.min(e_lo.value), np.max(e_lo.value), np.min(e_lo.value), np.max(e_lo.value)])
-    gs_ax3.set_xlabel(f"Count Energy [{e_lo.unit:latex}]")
-    gs_ax3.set_title(f"{f_rmf.split('/')[-1]}", size=8)
-    gs_ax3.set_aspect('equal', 'box')
-    gs_ax3.set_yticks([])
-    gs_ax3.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## total response
-    gs_ax4 = fig.add_subplot(gs[0, 9:11])
-    title = "Pos. 2, CdTe 4 Response"
-    gs_ax4.set_title(title)
-    gs_ax4.set_aspect('equal', 'box')
-    gs_ax4.set_yticks([])
-    pos2_resp = rmf * att_therm_bl[:,None] * oa_eff_area[:,None] * cdte_fixed2[:,None]
-    gs_ax4.imshow(pos2_resp.value, origin="lower", norm=LogNorm(vmin=0.001, vmax=0.12), extent=[np.min(e_lo.value), np.max(e_lo.value), np.min(e_lo.value), np.max(e_lo.value)])
-    gs_ax4.set_xlabel(f"Count Energy [{e_lo.unit:latex}]")
-    gs_ax4.annotate("=", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    plt.tight_layout()
-    plt.subplots_adjust(left=None, bottom=None, right=None, top=None, wspace=0.01, hspace=None)
-
-    if save_assets:
-        plt.savefig(os.path.join(assets_dir,f"{title}.png"), dpi=200, bbox_inches="tight")
-    plt.show()
-
-    ###########################################################################
-    ###########################################################################
-
-    ###########################################################################
-    ########################## CdTe 2 (pos. 3) ################################
-    ###########################################################################
-
-    fig = plt.figure(figsize=FIGSIZE)
-    gs = gridspec.GridSpec(*GS)
-
-    # define det response first so we know what energies we actually care about
-    f_rmf = os.path.join(pathlib.Path(__file__).parent, "det-resp-data", "cdte", "merged", "Resp_3keVto30keV_CdTe1_reg0_1hit.rmf")
-    e_lo, e_hi, rmf = cdte_det_resp_rmf(f_rmf)
-    mid_energies = (e_lo+e_hi)/2
-    photon_energy_lims = [np.nanmin(mid_energies.value), np.nanmax(mid_energies.value)]
-
-    ## thermal blanket
-    gs_ax0 = fig.add_subplot(gs[0, :2])
-    att_therm_bl = zeroes2nans(att_thermal_blanket(mid_energies))
-    gs_ax0.set_xlabel(f"Transmission [{att_therm_bl.unit:latex}]")
-    gs_ax0.set_ylabel(f"Photon Energy [{e_lo.unit:latex}]")
-    gs_ax0.set_title("Thermal blanket")
-    gs_ax0.plot(att_therm_bl, mid_energies, color=tb_col, ls="-")
-
-    y_lims = [0,1.01]
-    gs_ax0.set_xlim(y_lims)
-    gs_ax0.set_ylim(photon_energy_lims)
-    gs_ax0.set_aspect(np.diff(y_lims)/np.diff(photon_energy_lims))
-
-    ## Marshall 2-shell, high-res
-    optic = "X09/FM2"
-    gs_ax1 = fig.add_subplot(gs[0, 2:4])
-    gs_ax1.set_title(f"Marshall 2-shell ({optic})")
-    gs_ax1.set_aspect('equal', 'box')
-    gs_ax1.set_yticks([])
-    _eff_areas = eff_area_msfc_hi_res(mid_energies, position=3)
-    gs_ax1.plot(_eff_areas, mid_energies, label=f"{oaa:latex}", lw=_lw)
-    gs_ax1.set_xlabel(f"{optic} [{_eff_areas.unit:latex}]")
-    y_lims = [0, np.nanmax(_eff_areas).value*1.01]
-    gs_ax1.set_xlim(y_lims)
-    gs_ax1.set_ylim(photon_energy_lims)
-    gs_ax1.set_aspect(np.diff(y_lims)/np.diff(photon_energy_lims))
-    gs_ax1.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## Al Mylar
-    gs_ax2 = fig.add_subplot(gs[0, 4:6])
-    gs_ax2.set_title("Al Mylar")
-    al_mylar_att = zeroes2nans(att_al_mylar(mid_energies))
-    plt.plot(al_mylar_att, mid_energies, label="CdTe fixed p2")
-    gs_ax2.set_xlabel(f"Transmission [{al_mylar_att.unit:latex}]")
-    y_lims = [0,1.01]
-    gs_ax2.set_xlim(y_lims)
-    gs_ax2.set_ylim(photon_energy_lims)
-    gs_ax2.set_aspect(np.diff(y_lims)/np.diff(photon_energy_lims))
-    gs_ax2.set_yticks([])
-    gs_ax2.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## Pix. Att.
-    gs_ax3 = fig.add_subplot(gs[0, 6:8])
-    gs_ax3.set_title("Pixelated Attenuator")
-    pix_att = zeroes2nans(att_pixelated(mid_energies, values="modelled"))
-    plt.plot(pix_att, mid_energies, label="Pix. Att., pos. 5")
-    gs_ax3.set_xlabel(f"Transmission [{pix_att.unit:latex}]")
-    y_lims = [0,1.01]
-    gs_ax3.set_xlim(y_lims)
-    gs_ax3.set_ylim(photon_energy_lims)
-    gs_ax3.set_aspect(np.diff(y_lims)/np.diff(photon_energy_lims))
-    gs_ax3.set_yticks([])
-    gs_ax3.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## no2021_07
-    gs_ax4 = fig.add_subplot(gs[0, 8:10])
-    gs_ax4.imshow(rmf.value, origin="lower", norm=LogNorm(vmin=0.001, vmax=0.12), extent=[np.min(e_lo.value), np.max(e_lo.value), np.min(e_lo.value), np.max(e_lo.value)])
-    gs_ax4.set_xlabel(f"Count Energy [{e_lo.unit:latex}]")
-    gs_ax4.set_title(f"{f_rmf.split('/')[-1]}", size=8)
-    gs_ax4.set_aspect('equal', 'box')
-    gs_ax4.set_yticks([])
-    gs_ax4.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## total response
-    gs_ax5 = fig.add_subplot(gs[0, 10:])
-    title = "Pos. 3, CdTe 2 Response"
-    gs_ax5.set_title(title)
-    gs_ax5.set_aspect('equal', 'box')
-    gs_ax5.set_yticks([])
-    pos2_resp = rmf * att_therm_bl[:,None] * _eff_areas[:,None] * al_mylar_att[:,None] * pix_att[:,None]
-    gs_ax5.imshow(pos2_resp.value, origin="lower", norm=LogNorm(vmin=0.001, vmax=0.12), extent=[np.min(e_lo.value), np.max(e_lo.value), np.min(e_lo.value), np.max(e_lo.value)])
-    gs_ax5.set_xlabel(f"Count Energy [{e_lo.unit:latex}]")
-    gs_ax5.annotate("=", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    plt.tight_layout()
-    plt.subplots_adjust(left=None, bottom=None, right=None, top=None, wspace=0.01, hspace=None)
-
-    if save_assets:
-        plt.savefig(os.path.join(assets_dir,f"{title}.png"), dpi=200, bbox_inches="tight")
-    plt.show()
-
-    ###########################################################################
-    ###########################################################################
-
-    ###########################################################################
-    ########################## CdTe 3 (pos. 4) ################################
-    ###########################################################################
-
-    fig = plt.figure(figsize=FIGSIZE)
-    gs = gridspec.GridSpec(*GS)
-
-    # define det response first so we know what energies we actually care about
-    f_rmf = os.path.join(pathlib.Path(__file__).parent, "det-resp-data", "cdte", "merged", "Resp_3keVto30keV_CdTe1_reg0_1hit.rmf")
-    e_lo, e_hi, rmf = cdte_det_resp_rmf(f_rmf)
-    mid_energies = (e_lo+e_hi)/2
-    photon_energy_lims = [np.nanmin(mid_energies.value), np.nanmax(mid_energies.value)]
-
-    ## thermal blanket
-    gs_ax0 = fig.add_subplot(gs[0, 1:3])
-    att_therm_bl = zeroes2nans(att_thermal_blanket(mid_energies))
-    gs_ax0.set_xlabel(f"Transmission [{att_therm_bl.unit:latex}]")
-    gs_ax0.set_ylabel(f"Photon Energy [{e_lo.unit:latex}]")
-    gs_ax0.set_title("Thermal blanket")
-    gs_ax0.plot(att_therm_bl, mid_energies, color=tb_col, ls="-")
-
-    y_lims = [0,1.01]
-    gs_ax0.set_xlim(y_lims)
-    gs_ax0.set_ylim(photon_energy_lims)
-    gs_ax0.set_aspect(np.diff(y_lims)/np.diff(photon_energy_lims))
-
-    ## Nagoya 1-shell
-    gs_ax1 = fig.add_subplot(gs[0, 3:5])
-    gs_ax1.set_title("Nagoya 1-Shell HXT")
-    gs_ax1.set_aspect('equal', 'box')
-    gs_ax1.set_yticks([])
-    gs_ax1.annotate("?", (0.5, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-    y_lims = [0, np.nanmax(oa_eff_area).value*1.01]
-    gs_ax1.set_xlim(y_lims)
-    gs_ax1.set_ylim(photon_energy_lims)
-    gs_ax1.set_aspect(np.diff(y_lims)/np.diff(photon_energy_lims))
-    gs_ax1.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## Al filter
-    gs_ax2 = fig.add_subplot(gs[0, 5:7])
-    gs_ax2.set_title("Al (0.005\")")
-    cdte_fixed2 = zeroes2nans(att_uniform_al_cdte(mid_energies, file=None, position=4))
-    plt.plot(cdte_fixed2, mid_energies, label="CdTe fixed p4")
-    gs_ax2.set_xlabel(f"Transmission [{cdte_fixed2.unit:latex}]")
-    y_lims = [0,1.01]
-    gs_ax2.set_xlim(y_lims)
-    gs_ax2.set_ylim(photon_energy_lims)
-    gs_ax2.set_aspect(np.diff(y_lims)/np.diff(photon_energy_lims))
-    gs_ax2.set_yticks([])
-    gs_ax2.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## no2021_07
-    gs_ax3 = fig.add_subplot(gs[0, 7:9])
-    gs_ax3.imshow(rmf.value, origin="lower", norm=LogNorm(vmin=0.001, vmax=0.12), extent=[np.min(e_lo.value), np.max(e_lo.value), np.min(e_lo.value), np.max(e_lo.value)])
-    gs_ax3.set_xlabel(f"Count Energy [{e_lo.unit:latex}]")
-    gs_ax3.set_title(f"{f_rmf.split('/')[-1]}", size=8)
-    gs_ax3.set_aspect('equal', 'box')
-    gs_ax3.set_yticks([])
-    gs_ax3.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## total response
-    gs_ax4 = fig.add_subplot(gs[0, 9:11])
-    title = "Pos. 4, CdTe 3 Response"
-    gs_ax4.set_title(title)
-    gs_ax4.set_aspect('equal', 'box')
-    gs_ax4.set_yticks([])
-    gs_ax4.annotate("?", (0.5, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-    gs_ax4.set_xlabel(f"Count Energy [{e_lo.unit:latex}]")
-    gs_ax4.annotate("=", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    plt.tight_layout()
-    plt.subplots_adjust(left=None, bottom=None, right=None, top=None, wspace=0.01, hspace=None)
-    
-    if save_assets:
-        plt.savefig(os.path.join(assets_dir,f"{title}.png"), dpi=200, bbox_inches="tight")
-    plt.show()
-
-
-    ###########################################################################
-    ###########################################################################
-
-    ###########################################################################
-    ########################## CdTe 1 (pos. 5) ################################
-    ###########################################################################
-
-    fig = plt.figure(figsize=FIGSIZE)
-    gs = gridspec.GridSpec(*GS)
-
-    # define det response first so we know what energies we actually care about
-    f_rmf = os.path.join(pathlib.Path(__file__).parent, "det-resp-data", "cdte", "merged", "Resp_3keVto30keV_CdTe1_reg0_1hit.rmf")
-    e_lo, e_hi, rmf = cdte_det_resp_rmf(f_rmf)
-    mid_energies = (e_lo+e_hi)/2
-    photon_energy_lims = [np.nanmin(mid_energies.value), np.nanmax(mid_energies.value)]
-
-    ## thermal blanket
-    gs_ax0 = fig.add_subplot(gs[0, :2])
-    att_therm_bl = zeroes2nans(att_thermal_blanket(mid_energies))
-    gs_ax0.set_xlabel(f"Transmission [{att_therm_bl.unit:latex}]")
-    gs_ax0.set_ylabel(f"Photon Energy [{e_lo.unit:latex}]")
-    gs_ax0.set_title("Thermal blanket")
-    gs_ax0.plot(att_therm_bl, mid_energies, color=tb_col, ls="-")
-
-    y_lims = [0,1.01]
-    gs_ax0.set_xlim(y_lims)
-    gs_ax0.set_ylim(photon_energy_lims)
-    gs_ax0.set_aspect(np.diff(y_lims)/np.diff(photon_energy_lims))
-
-    ## Marshall 10-shell
-    optic = "X-8"
-    gs_ax1 = fig.add_subplot(gs[0, 2:4])
-    gs_ax1.set_title(f"Marshall 10-shell ({optic})")
-    gs_ax1.set_aspect('equal', 'box')
-    gs_ax1.set_yticks([])
-    _ps = []
-    for oaa in [-8.64, -6.72, -4.8, -2.88, -0.96, 0., 0.96, 2.88, 4.8, 6.72, 8.64]<<u.arcmin:
-        _lw = 1 if oaa >=0 else 2
-        _eff_areas = eff_area_msfc_10shell(mid_energies, off_axis=oaa<<u.arcmin, optic_id=optic)
-        _p = gs_ax1.plot(_eff_areas, mid_energies, label=f"{oaa:latex}", lw=_lw)
-        _ps += _p
-    gs_ax1.set_xlabel(f"{optic} [{_eff_areas.unit:latex}]")
-    plt.legend(handles=_ps, prop={"size": 6})
-    oa_eff_area = eff_area_msfc_10shell(mid_energies, off_axis=0<<u.arcmin, optic_id=optic)
-    y_lims = [0, np.nanmax(oa_eff_area).value*1.01]
-    gs_ax1.set_xlim(y_lims)
-    gs_ax1.set_ylim(photon_energy_lims)
-    gs_ax1.set_aspect(np.diff(y_lims)/np.diff(photon_energy_lims))
-    gs_ax1.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## Al Mylar
-    gs_ax2 = fig.add_subplot(gs[0, 4:6])
-    gs_ax2.set_title("Al Mylar")
-    al_mylar_att = zeroes2nans(att_al_mylar(mid_energies))
-    plt.plot(al_mylar_att, mid_energies, label="CdTe fixed p2")
-    gs_ax2.set_xlabel(f"Transmission [{al_mylar_att.unit:latex}]")
-    y_lims = [0,1.01]
-    gs_ax2.set_xlim(y_lims)
-    gs_ax2.set_ylim(photon_energy_lims)
-    gs_ax2.set_aspect(np.diff(y_lims)/np.diff(photon_energy_lims))
-    gs_ax2.set_yticks([])
-    gs_ax2.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## Pix. Att.
-    gs_ax3 = fig.add_subplot(gs[0, 6:8])
-    gs_ax3.set_title("Pixelated Attenuator")
-    pix_att = zeroes2nans(att_pixelated(mid_energies, values="modelled"))
-    plt.plot(pix_att, mid_energies, label="Pix. Att., pos. 5")
-    gs_ax3.set_xlabel(f"Transmission [{pix_att.unit:latex}]")
-    y_lims = [0,1.01]
-    gs_ax3.set_xlim(y_lims)
-    gs_ax3.set_ylim(photon_energy_lims)
-    gs_ax3.set_aspect(np.diff(y_lims)/np.diff(photon_energy_lims))
-    gs_ax3.set_yticks([])
-    gs_ax3.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## no2021_07
-    gs_ax4 = fig.add_subplot(gs[0, 8:10])
-    gs_ax4.imshow(rmf.value, origin="lower", norm=LogNorm(vmin=0.001, vmax=0.12), extent=[np.min(e_lo.value), np.max(e_lo.value), np.min(e_lo.value), np.max(e_lo.value)])
-    gs_ax4.set_xlabel(f"Count Energy [{e_lo.unit:latex}]")
-    gs_ax4.set_title(f"{f_rmf.split('/')[-1]}", size=8)
-    gs_ax4.set_aspect('equal', 'box')
-    gs_ax4.set_yticks([])
-    gs_ax4.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## total response
-    gs_ax5 = fig.add_subplot(gs[0, 10:])
-    title = "Pos. 5, CdTe 1 Response"
-    gs_ax5.set_title(title)
-    gs_ax5.set_aspect('equal', 'box')
-    gs_ax5.set_yticks([])
-    pos2_resp = rmf * att_therm_bl[:,None] * oa_eff_area[:,None] * al_mylar_att[:,None] * pix_att[:, None]
-    gs_ax5.imshow(pos2_resp.value, origin="lower", norm=LogNorm(vmin=0.001, vmax=0.12), extent=[np.min(e_lo.value), np.max(e_lo.value), np.min(e_lo.value), np.max(e_lo.value)])
-    gs_ax5.set_xlabel(f"Count Energy [{e_lo.unit:latex}]")
-    gs_ax5.annotate("=", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    plt.tight_layout()
-    plt.subplots_adjust(left=None, bottom=None, right=None, top=None, wspace=0.01, hspace=None)
-    
-    if save_assets:
-        plt.savefig(os.path.join(assets_dir,f"{title}.png"), dpi=200, bbox_inches="tight")
-    plt.show()
-
-
-    ###########################################################################
-    ###########################################################################
-
-
-    ###########################################################################
-    ###########################################################################
-
-    ###########################################################################
-    ########################## CMOS 1 (pos. 0) ################################
-    ###########################################################################
-
-    fig = plt.figure(figsize=FIGSIZE)
-    GS = (1, 7)
-    gs = gridspec.GridSpec(*GS)
-
-    # PATH: optical_blocking_filter_transmittance * collimator_aperture_ratio * mirror_effective_area * attenuation_filter_transmittance * sensor_quantum_efficiency @ response_matrix
-
-    # define det response first so we know what energies we actually care about
-    telescope = 0
-    counts, mid_energies, rmf = cmos_det_resp(telescope=telescope)
-    photon_energy_lims = [np.nanmin(mid_energies.value), np.nanmax(mid_energies.value)]
-
-    ## optical_blocking_filter_transmittance
-    gs_ax0 = fig.add_subplot(gs[0, 0])
-    obf1 = zeroes2nans(att_cmos_obfilter(mid_energies, telescope=telescope))
-    gs_ax0.set_xlabel(f"Transmission [{obf1.unit:latex}]")
-    gs_ax0.set_ylabel(f"Photon Energy [{mid_energies.unit:latex}]")
-    gs_ax0.set_title(f"OBF{telescope}")
-    gs_ax0.plot(obf1, mid_energies, color=tb_col, ls="-")
-
-    y_lims = [0,1.01]
-    gs_ax0.set_xlim(y_lims)
-    gs_ax0.set_ylim(photon_energy_lims)
-    gs_ax0.set_aspect(np.diff(y_lims)/np.diff(photon_energy_lims))
-
-    ## collimator_aperture_ratio
-    gs_ax1 = fig.add_subplot(gs[0, 1])
-    col1 = att_cmos_collimator_ratio(0<<u.arcsec, telescope=telescope)
-    gs_ax1.set_yticks([])
-    gs_ax1.set_xticks([])
-    gs_ax1.annotate(f"Collimator\naperture\nratio\n{col1}", (0.5, 0.5), xycoords="axes fraction", color="k", size=25, weight="bold", ha="center", va="center")
-    gs_ax1.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-    gs_ax1.set_aspect("equal", "box")
-
-    ## mirror_effective_area
-    gs_ax2 = fig.add_subplot(gs[0, 2])
-    cmos_ea1 = zeroes2nans(eff_area_cmos(mid_energies, telescope=telescope))
-    gs_ax2.set_xlabel(f"Effective Area [{cmos_ea1.unit:latex}]")
-    gs_ax2.set_yticks([])
-    gs_ax2.set_title(f"CMOS{telescope+1} Mirror")
-    gs_ax2.plot(cmos_ea1, mid_energies, color=tb_col, ls="-")
-
-    x_lims = [0, np.nanmax(cmos_ea1).value*1.01]
-    gs_ax2.set_xlim(x_lims)
-    gs_ax2.set_ylim(photon_energy_lims)
-    gs_ax2.set_aspect(np.diff(x_lims)/np.diff(photon_energy_lims))
-    gs_ax2.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## attenuation_filter_transmittance
-    gs_ax3 = fig.add_subplot(gs[0, 3])
-    filter1 = zeroes2nans(att_cmos_filter(mid_energies, telescope=telescope))
-    gs_ax3.set_xlabel(f"Transmission [{filter1.unit:latex}]")
-    gs_ax3.set_yticks([])
-    gs_ax3.set_title(f"Filter{telescope}")
-    gs_ax3.plot(filter1, mid_energies, color=tb_col, ls="-")
-
-    y_lims = [0,1.01]
-    gs_ax3.set_xlim(y_lims)
-    gs_ax3.set_ylim(photon_energy_lims)
-    gs_ax3.set_aspect(np.diff(y_lims)/np.diff(photon_energy_lims))
-    gs_ax3.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## sensor_quantum_efficiency
-    gs_ax4 = fig.add_subplot(gs[0, 4])
-    qe1 = zeroes2nans(qe_cmos(mid_energies, telescope=telescope))
-    gs_ax4.set_xlabel(f"Quantum Efficiency [{qe1.unit:latex}]")
-    gs_ax4.set_yticks([])
-    gs_ax4.set_title(f"CMOS{telescope+1} QE")
-    gs_ax4.plot(qe1, mid_energies, color=tb_col, ls="-")
-
-    y_lims = [0,1.01]
-    gs_ax4.set_xlim(y_lims)
-    gs_ax4.set_ylim(photon_energy_lims)
-    gs_ax4.set_aspect(np.diff(y_lims)/np.diff(photon_energy_lims))
-    gs_ax4.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## response_matrix
-    gs_ax5 = fig.add_subplot(gs[0, 5])
-    extent = [np.min(counts.value), np.max(counts.value), np.min(mid_energies.value), np.max(mid_energies.value)]
-    gs_ax5.set_ylim([0,10])
-    gs_ax5.set_yticks([])
-    gs_ax5.set_title(f"CMOS{telescope+1} RMF")
-    gs_ax5.set_aspect('equal', 'box')
-    gs_ax5.imshow(rmf.value, origin="lower", aspect=(extent[1]-extent[0])/(extent[3]-0), extent=extent, norm=LogNorm())
-    gs_ax5.set_xlabel(f"Counts? [{counts.unit:latex}]")
-    gs_ax5.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## total response
-    gs_ax6 = fig.add_subplot(gs[0, 6])
-    title = f"Pos. {telescope}, CMOS {telescope+1} Response"
-    gs_ax6.set_title(title)
-    gs_ax6.set_aspect('equal', 'box')
-    gs_ax6.set_yticks([])
-    gs_ax6.set_ylim([0,10])
-    pos0_resp = rmf * obf1[:,None] * col1 * cmos_ea1[:,None] * filter1[:,None] * qe1[:, None]
-    gs_ax6.imshow(pos0_resp.value, origin="lower", aspect=(extent[1]-extent[0])/(extent[3]-0), extent=extent, norm=LogNorm())
-    gs_ax6.set_xlabel(f"Counts? [{counts.unit:latex}]")
-    gs_ax6.annotate("=", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    plt.tight_layout()
-    plt.subplots_adjust(left=None, bottom=None, right=None, top=None, wspace=0.01, hspace=None)
-    
-    if save_assets:
-        plt.savefig(os.path.join(assets_dir,f"{title}.png"), dpi=200, bbox_inches="tight")
-    plt.show()
-
-
-    ###########################################################################
-    ###########################################################################
-
-
-    ###########################################################################
-    ###########################################################################
-
-    ###########################################################################
-    ########################## CMOS 2 (pos. 1) ################################
-    ###########################################################################
-
-    fig = plt.figure(figsize=FIGSIZE)
-    GS = (1, 7)
-    gs = gridspec.GridSpec(*GS)
-
-    # PATH: optical_blocking_filter_transmittance * collimator_aperture_ratio * mirror_effective_area * attenuation_filter_transmittance * sensor_quantum_efficiency @ response_matrix
-
-    # define det response first so we know what energies we actually care about
-    telescope = 1
-    counts, mid_energies, rmf = cmos_det_resp(telescope=telescope)
-    photon_energy_lims = [np.nanmin(mid_energies.value), np.nanmax(mid_energies.value)]
-
-    ## optical_blocking_filter_transmittance
-    gs_ax0 = fig.add_subplot(gs[0, 0])
-    obf2 = zeroes2nans(att_cmos_obfilter(mid_energies, telescope=telescope))
-    gs_ax0.set_xlabel(f"Transmission [{obf2.unit:latex}]")
-    gs_ax0.set_ylabel(f"Photon Energy [{mid_energies.unit:latex}]")
-    gs_ax0.set_title(f"OBF{telescope}")
-    gs_ax0.plot(obf2, mid_energies, color=tb_col, ls="-")
-
-    y_lims = [0,1.01]
-    gs_ax0.set_xlim(y_lims)
-    gs_ax0.set_ylim(photon_energy_lims)
-    gs_ax0.set_aspect(np.diff(y_lims)/np.diff(photon_energy_lims))
-
-    ## collimator_aperture_ratio
-    gs_ax1 = fig.add_subplot(gs[0, 1])
-    col2 = att_cmos_collimator_ratio(0<<u.arcsec, telescope=telescope)
-    gs_ax1.set_yticks([])
-    gs_ax1.set_xticks([])
-    gs_ax1.annotate(f"Collimator\naperture\nratio\n{col2}", (0.5, 0.5), xycoords="axes fraction", color="k", size=25, weight="bold", ha="center", va="center")
-    gs_ax1.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-    gs_ax1.set_aspect("equal", "box")
-
-    ## mirror_effective_area
-    gs_ax2 = fig.add_subplot(gs[0, 2])
-    cmos_ea2 = zeroes2nans(eff_area_cmos(mid_energies, telescope=telescope))
-    gs_ax2.set_xlabel(f"Effective Area [{cmos_ea2.unit:latex}]")
-    gs_ax2.set_yticks([])
-    gs_ax2.set_title(f"CMOS{telescope+1} Mirror")
-    gs_ax2.plot(cmos_ea2, mid_energies, color=tb_col, ls="-")
-
-    x_lims = [0, np.nanmax(cmos_ea2).value*1.01]
-    gs_ax2.set_xlim(x_lims)
-    gs_ax2.set_ylim(photon_energy_lims)
-    gs_ax2.set_aspect(np.diff(x_lims)/np.diff(photon_energy_lims))
-    gs_ax2.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## attenuation_filter_transmittance
-    gs_ax3 = fig.add_subplot(gs[0, 3])
-    filter2 = zeroes2nans(att_cmos_filter(mid_energies, telescope=telescope))
-    gs_ax3.set_xlabel(f"Transmission [{filter2.unit:latex}]")
-    gs_ax3.set_yticks([])
-    gs_ax3.set_title(f"Filter{telescope}")
-    gs_ax3.plot(filter2, mid_energies, color=tb_col, ls="-")
-
-    y_lims = [0,1.01]
-    gs_ax3.set_xlim(y_lims)
-    gs_ax3.set_ylim(photon_energy_lims)
-    gs_ax3.set_aspect(np.diff(y_lims)/np.diff(photon_energy_lims))
-    gs_ax3.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## sensor_quantum_efficiency
-    gs_ax4 = fig.add_subplot(gs[0, 4])
-    qe2 = zeroes2nans(qe_cmos(mid_energies, telescope=telescope))
-    gs_ax4.set_xlabel(f"Quantum Efficiency [{qe2.unit:latex}]")
-    gs_ax4.set_yticks([])
-    gs_ax4.set_title(f"CMOS{telescope+1} QE")
-    gs_ax4.plot(qe2, mid_energies, color=tb_col, ls="-")
-
-    y_lims = [0,1.01]
-    gs_ax4.set_xlim(y_lims)
-    gs_ax4.set_ylim(photon_energy_lims)
-    gs_ax4.set_aspect(np.diff(y_lims)/np.diff(photon_energy_lims))
-    gs_ax4.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## response_matrix
-    gs_ax5 = fig.add_subplot(gs[0, 5])
-    extent = [np.min(counts.value), np.max(counts.value), np.min(mid_energies.value), np.max(mid_energies.value)]
-    gs_ax5.set_ylim([0,10])
-    gs_ax5.set_yticks([])
-    gs_ax5.set_title(f"CMOS{telescope+1} RMF")
-    gs_ax5.set_aspect('equal', 'box')
-    gs_ax5.imshow(rmf.value, origin="lower", aspect=(extent[1]-extent[0])/(extent[3]-0), extent=extent, norm=LogNorm())
-    gs_ax5.set_xlabel(f"Counts? [{counts.unit:latex}]")
-    gs_ax5.annotate("x", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    ## total response
-    gs_ax6 = fig.add_subplot(gs[0, 6])
-    title = f"Pos. {telescope}, CMOS {telescope+1} Response"
-    gs_ax6.set_title(title)
-    gs_ax6.set_aspect('equal', 'box')
-    gs_ax6.set_yticks([])
-    gs_ax6.set_ylim([0,10])
-    pos1_resp = rmf * obf2[:,None] * col2 * cmos_ea2[:,None] * filter2[:,None] * qe2[:, None]
-    gs_ax6.imshow(pos1_resp.value, origin="lower", aspect=(extent[1]-extent[0])/(extent[3]-0), extent=extent, norm=LogNorm())
-    gs_ax6.set_xlabel(f"Counts? [{counts.unit:latex}]")
-    gs_ax6.annotate("=", (0, 0.5), xycoords="axes fraction", color="red", size=25, weight="bold", ha="center", va="center")
-
-    plt.tight_layout()
-    plt.subplots_adjust(left=None, bottom=None, right=None, top=None, wspace=0.01, hspace=None)
-    
-    if save_assets:
-        plt.savefig(os.path.join(assets_dir,f"{title}.png"), dpi=200, bbox_inches="tight")
-    plt.show()
-
-if __name__=="__main__":
-    foxsi4_response_paths(save_assets=False)
\ No newline at end of file
diff --git a/response-tools-py/build/lib/response_tools_py/phot_spec.py b/response-tools-py/build/lib/response_tools_py/phot_spec.py
deleted file mode 100644
index 3597df0..0000000
--- a/response-tools-py/build/lib/response_tools_py/phot_spec.py
+++ /dev/null
@@ -1,100 +0,0 @@
-"""Code to create photon models. """
-
-import astropy.units as u
-from copy import deepcopy
-import numpy as np
-from scipy.ndimage import gaussian_filter1d
-from scipy.io import readsav
-
-def get_energy_delta():
-    """Return energy bin width."""
-    return 0.0445 << u.keV
-
-def create_energy_edges():
-    """One function to handle giving out the same energies."""
-    # no longer needed with above files `lower_e = 1.0002920302956426` # from lower bounds thermal() error, in KeV
-    lower_e = 0.5 # from lower bounds thermal() error, in KeV
-    return np.arange(lower_e, 100, get_energy_delta().value) << u.keV
-
-def create_energy_midpoints(edges=None):
-    """Tied to `create_energy_edges`."""
-    _edges = create_energy_edges() if edges is None else edges
-    return (_edges[1:]+_edges[:-1])/2
-
-def zeroes2nans(array):
-    """Replace 0s in the array with np.nans."""
-    _array_copy = deepcopy(array)
-    _array_copy[np.nonzero(array==0)] = np.nan
-    return _array_copy
-
-def nans2zeroes(array):
-    """Replace Nans in the array with Os."""
-    _array_copy = deepcopy(array)
-    _array_copy[np.isnan(array)] = 0
-    return _array_copy
-
-def gaussian_blur(array, fwhm):
-    """Apply a Gaussian blur (with FWHM) to an array. 
-    
-    Parameters
-    ----------
-    array : `numpy.ndarray`
-        Array to which apply blurring.
-        
-    fwhm : `int`, `float`
-        Full-width half-max of the Gaussian. must be in units of "number
-        of bins" for the function,  not anything else like energy, etc.
-    """
-    sigma = fwhm/2.355 # std for kernal in units of "number of bins"
-    return gaussian_filter1d(array, sigma)
-
-def sxr_res():
-    """Return the SXR energy resolution (from Albert).
-    
-    FOXSI-4: 0.4 keV @ 10 keV
-    Albert: 0.2 keV
-    """
-    return 0.2 << u.keV
-
-def sxr_blur(array):
-    """Apply the FOXSI-4 SXR energy resolution to the incoming array. """
-    return gaussian_blur(array, sxr_res()/get_energy_delta()) << array.unit
-
-def hxr_res():
-    """Return the HXR energy resolution (from FOXSI-4).
-    
-    0.8 keV @ 14 keV
-    """
-    return 0.8 << u.keV
-
-def hxr_blur(array):
-    """Apply the FOXSI-4 HXR energy resolution to the incoming array. """
-    return gaussian_blur(array, hxr_res()/get_energy_delta()) << array.unit
-
-def _load_idl_sim(filename):
-    """Function to return simulated spectrum."""
-    results = readsav(filename)
-
-    energy_edges = results["e_hist_bins"]
-    full_loop = results["spatially_resolved_I_photon"][:,0]
-    loop_top = results["spatially_resolved_I_photon"][:,1]
-    mean_footpoint = results["spatially_resolved_I_photon"][:,2]
-
-    return {"model":{"energy_edges":energy_edges << u.keV, "full_loop":full_loop << (u.ph / (u.keV * u.s * u.cm**2))}, 
-            "submodels":{"loop_top":loop_top << (u.ph / (u.keV * u.s * u.cm**2)), "footpoints":mean_footpoint << (u.ph / (u.keV * u.s * u.cm**2))}}
-
-def sim_energy_edges(*_, filename=None, **__):
-    """Access energy bin edges from Morgan's simulations."""
-    return _load_idl_sim(filename)["model"]["energy_edges"]
-
-def sim_full_loop(*_, filename=None, **__):
-    """Access full loop spectrum from Morgan's simulations."""
-    return _load_idl_sim(filename)["model"]["full_loop"]
-
-def sim_loop_top(*_, filename=None, **__):
-    """Access loop-top spectrum from Morgan's simulations."""
-    return _load_idl_sim(filename)["submodels"]["loop_top"]
-
-def sim_footpoints(*_, filename=None, **__):
-    """Access ONE footpoint spectrum from Morgan's simulations."""
-    return _load_idl_sim(filename)["submodels"]["footpoints"]
\ No newline at end of file
diff --git a/response-tools-py/build/lib/response_tools_py/quantum_eff.py b/response-tools-py/build/lib/response_tools_py/quantum_eff.py
deleted file mode 100644
index fe2cd0c..0000000
--- a/response-tools-py/build/lib/response_tools_py/quantum_eff.py
+++ /dev/null
@@ -1,50 +0,0 @@
-"""Code to load different quantum efficiencies. """
-
-import logging
-import os
-import pathlib
-
-from astropy.io import fits
-import astropy.units as u
-import numpy as np
-
-@u.quantity_input(mid_energies=u.keV)
-def qe_cmos(mid_energies, file=None, telescope=None):
-    """Return quantum efficency of CMOS detectors.
-    
-    Telescope 0: position 0, X10/FM2(?)
-    Telescope 1: position 1, Nagoya(?)
-    """
-    if telescope is None:
-        logging.warning("`telescope` input in `qe_cmos()` must be 0 or 1.")
-        return
-        
-    _f = os.path.join(pathlib.Path(__file__).parent, 
-                            "quantum-eff-data", 
-                            f"foxsi4_telescope-{telescope}_BASIC_sensor_quantum_efficiency_V25APR13.fits") if file is None else file
-    with fits.open(_f) as hdul:
-        es, qe = hdul[2].data << u.keV, hdul[1].data << u.dimensionless_unscaled
-    return np.interp(mid_energies.value, es.value, qe.value, left=0, right=0) << u.dimensionless_unscaled
-
-if __name__=="__main__":
-    import matplotlib.pyplot as plt
-
-    SAVE_ASSETS = False
-    assets_dir = os.path.join(pathlib.Path(__file__).parent, "..", "..", "assets", "response-tools-py-figs", "quantum-eff-figs")
-    pathlib.Path(assets_dir).mkdir(parents=True, exist_ok=True)
-
-    mid_energies = np.linspace(0, 20, 1000)<<u.keV
-    plt.figure()
-    qe0 = qe_cmos(mid_energies, telescope=0)
-    qe1 = qe_cmos(mid_energies, telescope=1)
-    plt.plot(mid_energies, qe0, lw=2, label="CMOS telescope 0, position 0")
-    plt.plot(mid_energies, qe1, lw=1, ls=":", label="CMOS telescope 1, position 1")
-    plt.title("CMOS Detectors")
-    plt.ylabel(f"Quantum efficiency [{qe0.unit:latex}]")
-    plt.xlabel(f"Energy [{mid_energies.unit:latex}]")
-    plt.xlim([np.min(mid_energies.value), np.max(mid_energies.value)])
-    plt.ylim([0, 1.01])
-    plt.legend()
-    if SAVE_ASSETS:
-        plt.savefig(os.path.join(assets_dir,"cmos-qunatum-eff.png"), dpi=200, bbox_inches="tight")
-    plt.show()
diff --git a/response-tools-py/build/lib/response_tools_py/version.py b/response-tools-py/build/lib/response_tools_py/version.py
deleted file mode 100644
index ef37232..0000000
--- a/response-tools-py/build/lib/response_tools_py/version.py
+++ /dev/null
@@ -1,3 +0,0 @@
-import importlib.metadata
-
-__version__ = importlib.metadata.version("response_tools_py")

From a3c1d8ada94e9220af42591d04aed8fa1c006897 Mon Sep 17 00:00:00 2001
From: Kristopher Cooper <kris@infra04-wg012.eduroam.wireless.umn.edu>
Date: Wed, 24 Sep 2025 16:05:20 -0500
Subject: [PATCH 2/3] fixed some typos

---
 .../examples/functions_and_outputs.py         | 46 +++++++++++++------
 .../response_tools_py/responses.py            | 36 +++++++++------
 2 files changed, 54 insertions(+), 28 deletions(-)

diff --git a/response-tools-py/examples/functions_and_outputs.py b/response-tools-py/examples/functions_and_outputs.py
index 7e40302..ba90ea6 100644
--- a/response-tools-py/examples/functions_and_outputs.py
+++ b/response-tools-py/examples/functions_and_outputs.py
@@ -8,7 +8,7 @@
 - Thermal blanket -> Marshall 10-shell X-7 -> Al (0.015") -> CdTe4
 
 Mainly a quick example on unit aware objects and how to access the 
-returned dataclass objects.
+returned data-class objects.
 
 This example shows the use of nice high level functions that are tied to
 FOXSI-4 telescopes. 
@@ -128,13 +128,13 @@
 when necessary.
 
 The output of `pos2_optics` (an `pos2_optics_new`) isn't just an array 
-now, it's a dataclass. The dataclass contains the effective areas of the 
-optics but also the energy, file, off-axis angle information used to 
+now, it's a data-class. The data-class contains the effective areas of 
+the optics but also the energy, file, off-axis angle information used to 
 produce it. This is crucial to track when there are a lot of files 
 flying around.
 
 As suggested in the "Output" section of the helpful documentation 
-earlier, we can see the contents of the dataclass and how to access the 
+earlier, we can see the contents of the data-class and how to access the 
 information within it:
 """
 
@@ -164,11 +164,28 @@
  'model': False}
 ```
 
+The above might be a bit cumbersome and you may only wish to look at the
+data fields contained in the output. For this we can simply run:
+"""
+print(pos2_optics.fields)
+"""
+yielding something like:
+```
+['filename',
+ 'function_path',
+ 'mid_energies',
+ 'off_axis_angle',
+ 'effective_areas',
+ 'optic_id',
+ 'model']
+```
+which might be a bit easier to read.
+
 Note: there is a method called `print_contents()` you can use on the 
 function output that might format the contents a little nicer than the
 above you may wish to use.
 
-Each field can be accessed with the displayed name. For excample, to get
+Each field can be accessed with the displayed name. For example, to get
 the effective areas of the optics, simply:
 """
 
@@ -184,10 +201,11 @@
 Redistribution Matrix Function (RMF), and/or Detector Response Matrix 
 (DRM).
 
-First, we can get the RMF for a telescope, say, Telescope 2:
+First, we can get the RMF for a telescope, say, Telescope 2 to be 
+consistent with using position 2's components previously:
 """
 
-pos2_rmf = responses.foxsi4_telescope2_rmf(region=0)
+tel2_rmf = responses.foxsi4_telescope2_rmf(region=0)
 
 """
 The `region` input refers to the different pitch regions across the CdTe
@@ -200,9 +218,9 @@
 the ARF values for:
 """
 
-mid_energies = (pos2_rmf.input_energy_edges[:-1]\
-                +pos2_rmf.input_energy_edges[1:])/2
-pos2_arf = responses.foxsi4_telescope2_arf(mid_energies=mid_energies, 
+mid_energies = (tel2_rmf.input_energy_edges[:-1]\
+                +tel2_rmf.input_energy_edges[1:])/2
+tel2_arf = responses.foxsi4_telescope2_arf(mid_energies=mid_energies, 
                                            off_axis_angle=0<<u.arcmin)
 
 """
@@ -211,20 +229,20 @@
 RMF to the general `responses.foxsi4_telescope_response` function:
 """
 
-pos2_drm = responses.foxsi4_telescope_response(pos2_arf, pos2_rmf)
+tel2_drm = responses.foxsi4_telescope_response(tel2_arf, tel2_rmf)
 
 """
 Note that with the FOXSI-4 Telescope fucntions (`foxsi4_telescope*`), 
 there exists a field in the class called `elements`. This field contains
-all the dataclasses used to produce the objects `response` field.
+all the data-classes used to produce the objects `response` field.
 
 Checking the `elements` field for, say, the ARF object:
 """
 
-print(pos2_arf.elements)
+print(tel2_arf.elements)
 
 """
-We find that the `elements` field contains all the dataclasses that 
+We find that the `elements` field contains all the data-classes that 
 produced the Telescope 2 ARF such as the thermal blanket transmission 
 (dimensionless), the Marshall 10-shell X-7 optics effective areas 
 (cm**2), and the Al (0.015") attenuator transmissions (dimensionless).
diff --git a/response-tools-py/response_tools_py/responses.py b/response-tools-py/response_tools_py/responses.py
index e785f4c..435f380 100644
--- a/response-tools-py/response_tools_py/responses.py
+++ b/response-tools-py/response_tools_py/responses.py
@@ -51,6 +51,11 @@ def foxsi4_telescope_response(arf_response, rmf_response):
     combining both the ARF and RMF into the Detector Response Matrix 
     (DRM).
 
+    This function will check both the ARF and RMF have come from the 
+    same telescope using the "telescope" field. A warning will be 
+    produce if they are different to inform the user but the function 
+    will still output the product if possible.
+
     Parameters
     ----------
     arf_response : `responses.Response1DOutput`
@@ -70,7 +75,10 @@ def foxsi4_telescope_response(arf_response, rmf_response):
     # check compatibility
     rmf_mids = (rmf_response.input_energy_edges[:-1]+rmf_response.input_energy_edges[1:])/2
     if not np.all(arf_response.mid_energies==rmf_mids):
-        raise ValueError("The `arf_response.mid_energies` do not match the bin centers of `rmf_response.input_energy_edges`.\nDRM product cannot be calculated")
+        raise ValueError(f"In {sys._getframe().f_code.co_name}, the `arf_response.mid_energies` do not match the bin centers of `rmf_response.input_energy_edges`.\nDRM product cannot be calculated")
+    
+    if arf_response.telescope!=rmf_response.telescope:
+        logging.warning(f"In {sys._getframe().f_code.co_name}, the \"telescope\" fields from `arf_response` ({arf_response.telescope}) and `rmf_response` ({rmf_response.telescope}) do not match.\nAn output shall be produced if possible (check the returned objects\"telescope\" field) but the user should be cautious.")
 
     total_response = arf_response.response[:,None] * rmf_response.response
 
@@ -84,7 +92,7 @@ def foxsi4_telescope_response(arf_response, rmf_response):
                             output_energy_edges=rmf_response.output_energy_edges,
                             response=total_response,
                             response_type="DRM",
-                            telescope="foxsi4-2",
+                            telescope=f"ARF:{arf_response.telescope},RMF{arf_response.telescope}",
                             elements=(arf_response,
                                       rmf_response,
                                       ),
@@ -132,7 +140,7 @@ def foxsi4_telescope2_arf(mid_energies, off_axis_angle):
                             mid_energies=mid_energies,
                             response=arf,
                             response_type="ARF",
-                            telescope="foxsi4-2",
+                            telescope="foxsi4-telescope2",
                             elements=(tb, 
                                       opt, 
                                       uni_al,
@@ -186,7 +194,7 @@ def foxsi4_telescope2_flight_arf(mid_energies, off_axis_angle, time_range):
                             mid_energies=mid_energies,
                             response=flight_arf,
                             response_type="ARF-flight",
-                            telescope=f"foxsi4-2",
+                            telescope=f"foxsi4-telescope2",
                             elements=(atm, 
                                       arf,
                                       ),
@@ -248,7 +256,7 @@ def foxsi4_telescope2_rmf(region:int=None, pitch=None, _side:str="merged", _even
                             output_energy_edges=rmf.output_energy_edges,
                             response=rmf.detector_response,
                             response_type="RMF",
-                            telescope="foxsi4-2",
+                            telescope="foxsi4-telescope2",
                             elements=(rmf,
                                       ),
                             )
@@ -299,7 +307,7 @@ def foxsi4_telescope3_arf(mid_energies, off_axis_angle=None):
                             mid_energies=mid_energies,
                             response=arf,
                             response_type="ARF",
-                            telescope="foxsi4-2",
+                            telescope="foxsi4-telescope3",
                             elements=(tb, 
                                       opt,
                                       mylar, 
@@ -357,7 +365,7 @@ def foxsi4_telescope3_flight_arf(mid_energies, time_range, off_axis_angle=None):
                             mid_energies=mid_energies,
                             response=flight_arf,
                             response_type="ARF-flight",
-                            telescope=f"foxsi4-3",
+                            telescope=f"foxsi4-telescope3",
                             elements=(atm, 
                                       arf,
                                       ),
@@ -419,7 +427,7 @@ def foxsi4_telescope3_rmf(region:int=None, pitch=None, _side:str="merged", _even
                             output_energy_edges=rmf.output_energy_edges,
                             response=rmf.detector_response,
                             response_type="RMF",
-                            telescope="foxsi4-3",
+                            telescope="foxsi4-telescope3",
                             elements=(rmf,
                                       ),
                             )
@@ -468,7 +476,7 @@ def foxsi4_telescope4_arf(mid_energies, off_axis_angle=None):
                             mid_energies=mid_energies,
                             response=arf,
                             response_type="ARF",
-                            telescope="foxsi4-4",
+                            telescope="foxsi4-telescope4",
                             elements=(tb, 
                                       opt, 
                                       uni_al,
@@ -525,7 +533,7 @@ def foxsi4_telescope4_flight_arf(mid_energies, time_range, off_axis_angle=None):
                             mid_energies=mid_energies,
                             response=flight_arf,
                             response_type="ARF-flight",
-                            telescope=f"foxsi4-4",
+                            telescope=f"foxsi4-telescope4",
                             elements=(atm, 
                                       arf,
                                       ),
@@ -587,7 +595,7 @@ def foxsi4_telescope4_rmf(region:int=None, pitch=None, _side:str="merged", _even
                             output_energy_edges=rmf.output_energy_edges,
                             response=rmf.detector_response,
                             response_type="RMF",
-                            telescope="foxsi4-4",
+                            telescope="foxsi4-telescope4",
                             elements=(rmf,
                                       ),
                             )
@@ -636,7 +644,7 @@ def foxsi4_telescope5_arf(mid_energies, off_axis_angle):
                             mid_energies=mid_energies,
                             response=arf,
                             response_type="ARF",
-                            telescope="foxsi4-5",
+                            telescope="foxsi4-telescope5",
                             elements=(tb, 
                                       opt, 
                                       mylar,
@@ -691,7 +699,7 @@ def foxsi4_telescope5_flight_arf(mid_energies, off_axis_angle, time_range):
                             mid_energies=mid_energies,
                             response=flight_arf,
                             response_type="ARF-flight",
-                            telescope=f"foxsi4-5",
+                            telescope=f"foxsi4-telescope5",
                             elements=(atm, 
                                       arf,
                                       ),
@@ -753,7 +761,7 @@ def foxsi4_telescope5_rmf(region:int=None, pitch=None, _side:str="merged", _even
                             output_energy_edges=rmf.output_energy_edges,
                             response=rmf.detector_response,
                             response_type="RMF",
-                            telescope="foxsi4-5",
+                            telescope="foxsi4-telescope5",
                             elements=(rmf,
                                       ),
                             )

From 33232461a60b9ce9900ce638febece239e5b58a1 Mon Sep 17 00:00:00 2001
From: Kristopher Cooper <kris@infra04-wg012.eduroam.wireless.umn.edu>
Date: Wed, 24 Sep 2025 16:28:37 -0500
Subject: [PATCH 3/3] restructure to Python-only package

---
 .gitignore                                    |   2 +-
 README.md                                     |  98 ++++++++++++------
 .../att-figs/atmospheric-transmissions.png    | Bin
 .../att-figs/model-transmissions.png          | Bin
 .../att-figs/transmissions.png                | Bin
 .../det-resp-figs/cdte-response-matrix.png    | Bin
 .../det-resp-figs/cmos-response-matrices.png  | Bin
 .../eff-area-figs/cmos-optics-resp.png        | Bin
 .../heritage-and-msfc-optics.png              | Bin
 .../quantum-eff-figs/cmos-qunatum-eff.png     | Bin
 .../response-figs/response-chain.png          | Bin
 .../response-hit-combinations.png             | Bin
 .../response-paths/Pos_0_CMOS_1_Response.png  | Bin
 .../response-paths/Pos_1_CMOS_2_Response.png  | Bin
 .../response-paths/Pos_2_CdTe_4_Response.png  | Bin
 .../response-paths/Pos_3_CdTe_2_Response.png  | Bin
 .../response-paths/Pos_4_CdTe_3_Response.png  | Bin
 .../response-paths/Pos_5_CdTe_1_Response.png  | Bin
 examples/README.md                            |  14 ++-
 .../functions_and_outputs.py                  |   6 +-
 .../examples => examples}/plot_arf_rmf_drm.py |   2 +-
 response-information/README.md                |   2 +-
 response-tools-py/README.md                   |  58 -----------
 response-tools-py/examples/README.md          |  13 ---
 .../response_tools_py/version.py              |   3 -
 .../__init__.py                               |  16 +--
 .../attenuation.py                            |   8 +-
 .../aux/README.md                             |   2 +-
 .../aux/format-msfc-eff-areas.py              |   0
 .../detector_response.py                      |   6 +-
 .../effective_area.py                         |   6 +-
 .../io/__init__.py                            |   0
 .../io/load_yaml.py                           |   2 +-
 .../io/tests/README.md                        |   2 +-
 .../io/tests/__init__.py                      |   0
 .../ir-plots.py                               |   2 +-
 .../phot_spec.py                              |   0
 .../quantum_efficiency.py                     |   6 +-
 .../responses.py                              |   8 +-
 .../telescope_parts.py                        |  22 ++--
 .../util.py                                   |   0
 response_tools/version.py                     |   3 +
 response-tools-py/setup.py => setup.py        |   2 +-
 tests/README.md                               |   2 +-
 44 files changed, 129 insertions(+), 156 deletions(-)
 rename assets/{response-tools-py-figs => response-tools-figs}/att-figs/atmospheric-transmissions.png (100%)
 rename assets/{response-tools-py-figs => response-tools-figs}/att-figs/model-transmissions.png (100%)
 rename assets/{response-tools-py-figs => response-tools-figs}/att-figs/transmissions.png (100%)
 rename assets/{response-tools-py-figs => response-tools-figs}/det-resp-figs/cdte-response-matrix.png (100%)
 rename assets/{response-tools-py-figs => response-tools-figs}/det-resp-figs/cmos-response-matrices.png (100%)
 rename assets/{response-tools-py-figs => response-tools-figs}/eff-area-figs/cmos-optics-resp.png (100%)
 rename assets/{response-tools-py-figs => response-tools-figs}/eff-area-figs/heritage-and-msfc-optics.png (100%)
 rename assets/{response-tools-py-figs => response-tools-figs}/quantum-eff-figs/cmos-qunatum-eff.png (100%)
 rename assets/{response-tools-py-figs => response-tools-figs}/response-figs/response-chain.png (100%)
 rename assets/{response-tools-py-figs => response-tools-figs}/response-figs/response-hit-combinations.png (100%)
 rename assets/{response-tools-py-figs => response-tools-figs}/response-paths/Pos_0_CMOS_1_Response.png (100%)
 rename assets/{response-tools-py-figs => response-tools-figs}/response-paths/Pos_1_CMOS_2_Response.png (100%)
 rename assets/{response-tools-py-figs => response-tools-figs}/response-paths/Pos_2_CdTe_4_Response.png (100%)
 rename assets/{response-tools-py-figs => response-tools-figs}/response-paths/Pos_3_CdTe_2_Response.png (100%)
 rename assets/{response-tools-py-figs => response-tools-figs}/response-paths/Pos_4_CdTe_3_Response.png (100%)
 rename assets/{response-tools-py-figs => response-tools-figs}/response-paths/Pos_5_CdTe_1_Response.png (100%)
 rename {response-tools-py/examples => examples}/functions_and_outputs.py (98%)
 rename {response-tools-py/examples => examples}/plot_arf_rmf_drm.py (98%)
 delete mode 100644 response-tools-py/README.md
 delete mode 100644 response-tools-py/examples/README.md
 delete mode 100644 response-tools-py/response_tools_py/version.py
 rename {response-tools-py/response_tools_py => response_tools}/__init__.py (57%)
 rename {response-tools-py/response_tools_py => response_tools}/attenuation.py (99%)
 rename {response-tools-py/response_tools_py => response_tools}/aux/README.md (97%)
 rename {response-tools-py/response_tools_py => response_tools}/aux/format-msfc-eff-areas.py (100%)
 rename {response-tools-py/response_tools_py => response_tools}/detector_response.py (99%)
 rename {response-tools-py/response_tools_py => response_tools}/effective_area.py (99%)
 rename {response-tools-py/response_tools_py => response_tools}/io/__init__.py (100%)
 rename {response-tools-py/response_tools_py => response_tools}/io/load_yaml.py (89%)
 rename {response-tools-py/response_tools_py => response_tools}/io/tests/README.md (56%)
 rename {response-tools-py/response_tools_py => response_tools}/io/tests/__init__.py (100%)
 rename {response-tools-py/response_tools_py => response_tools}/ir-plots.py (99%)
 rename {response-tools-py/response_tools_py => response_tools}/phot_spec.py (100%)
 rename {response-tools-py/response_tools_py => response_tools}/quantum_efficiency.py (94%)
 rename {response-tools-py/response_tools_py => response_tools}/responses.py (99%)
 rename {response-tools-py/response_tools_py => response_tools}/telescope_parts.py (97%)
 rename {response-tools-py/response_tools_py => response_tools}/util.py (100%)
 create mode 100644 response_tools/version.py
 rename response-tools-py/setup.py => setup.py (93%)

diff --git a/.gitignore b/.gitignore
index bd9ed0b..e4f75e7 100644
--- a/.gitignore
+++ b/.gitignore
@@ -6,7 +6,7 @@ __pycache__/
 .ipynb_checkpoints
 
 # Python
-response_tools_py.egg-info/
+response_tools.egg-info/
 build/
 
 # data directories
diff --git a/README.md b/README.md
index d28110c..4cfedff 100644
--- a/README.md
+++ b/README.md
@@ -6,6 +6,65 @@ Repository to manage response elements, originally for the FOXSI-4 sounding rock
 
 **Note:** Although this repositiory can be extanded to any other instrument (e.g., FOXSI-1, -2, and -3), this repository is only currently aimed for FOXSI-4 use.
 
+## `response-tools`
+
+A repository to software tools for the response system used for the FOXSI mission in Python (woo, Python). We will include helpful loader functions for the all the response elements and download clients to obtain the files.
+
+The tools being developed in Python should be placed in the `response_tools` folder (note the underscores and not dashes).
+
+There is an "examples" and a "tests" folder. The "examples" folder is a great place to include scripts that show how some of the code in the repository works and the "tests" folder is a fantastic place to put code that tests the tools that have been created.
+
+More information will be placed here with regards as to how this package is recommended to be used.
+
+## Install tips
+
+In order to work with some preliminary data, it is instructive to set up a virtual environment (more information below) and install necessary packages. The suggested way to do this is to create an environment with:
+
+- `conda create -n response-tools-env python==3.12`
+
+(recommend using Python 3.12 just because this has been proven to be stable with the software).
+
+The Python code can then be installed with:
+
+- `pip install -e .`
+
+while in this directory (the `response-tools` directory) containing the `setup.py` file.
+
+Any time the code is updated and, say, you pull it from Github, make sure to perform the `pip install -e .` line from above. This ensures any new changes, updated versions, etc. take effect.
+
+## Namespace
+
+The base `response_tools` namespace includes:
+
+- `~response_tools.contextResponseInfo`
+  - The information stored in a [YAML file](../response-information/info.yaml) that includes file version information (e.g., the current versions to be used) and can be accessed like a Python dictionary.
+- `~response_tools.__version__`
+  - The version of the code determined by `response-tools/setup.py`.
+    - `0.0.1`: First rendition of the code.
+
+### Examples using the namespace
+
+```python
+# importing the module is as easy as:
+>>> import response_tools
+
+# can access the YAML file contents in the code as 
+>>> print(response_tools.contextResponseInfo[""])
+...
+
+# then accessing the version as:
+>>> print(response_tools.__version__)
+'0.0.1'
+```
+
+## Useful Python tips
+
+Some useful things to keep in mind while using the Python code in this repository.
+
+### Virtual environments
+
+It might be a good idea to look into ([conda](https://conda.io/projects/conda/en/latest/user-guide/tasks/manage-environments.html)) virtual environments if you are not familiar, this includes looking into them yourself or getting in touch with someone to help explain. This requires downloading [miniconda](https://docs.anaconda.com/miniconda/install/#quick-command-line-install) (or just conda if you prefer).
+
 ## Response files
 
 There are many response elements for FOXSI-4. Once downloaded, they can be places or found [here](./response-information/).
@@ -26,7 +85,7 @@ For a historical description of the response files below (old versions, etc.), s
 - **Detector response:** ""
   - **Quantum efficiency:** ""
 
-![Position 0 photon path](./assets/response-tools-py-figs/response-paths/Pos_0_CMOS_1_Response.png "Position 0 photon path")
+![Position 0 photon path](./assets/response-tools-figs/response-paths/Pos_0_CMOS_1_Response.png "Position 0 photon path")
 
 ### Position 1 (SXR, CMOS2)
 
@@ -37,7 +96,7 @@ For a historical description of the response files below (old versions, etc.), s
 - **Detector response:** ""
   - **Quantum efficiency:** ""
 
-![Position 1 photon path](./assets/response-tools-py-figs/response-paths/Pos_1_CMOS_2_Response.png "Position 1 photon path")
+![Position 1 photon path](./assets/response-tools-figs/response-paths/Pos_1_CMOS_2_Response.png "Position 1 photon path")
 
 ### Position 2 (HXR, CdTe4)
 
@@ -46,7 +105,7 @@ For a historical description of the response files below (old versions, etc.), s
 - **Filter:** ""
 - **Detector response:** ""
 
-![Position 2 photon path](./assets/response-tools-py-figs/response-paths/Pos_2_CdTe_4_Response.png "Position 2 photon path")
+![Position 2 photon path](./assets/response-tools-figs/response-paths/Pos_2_CdTe_4_Response.png "Position 2 photon path")
 
 ### Position 3 (HXR, CdTe2)
 
@@ -56,7 +115,7 @@ For a historical description of the response files below (old versions, etc.), s
 - **Pixelated attenuator:** ""
 - **Detector response:** ""
 
-![Position 3 photon path](./assets/response-tools-py-figs/response-paths/Pos_3_CdTe_2_Response.png "Position 3 photon path")
+![Position 3 photon path](./assets/response-tools-figs/response-paths/Pos_3_CdTe_2_Response.png "Position 3 photon path")
 
 ### Position 4 (HXR, CdTe3)
 
@@ -65,7 +124,7 @@ For a historical description of the response files below (old versions, etc.), s
 - **Filter:** ""
 - **Detector response:** ""
 
-![Position 4 photon path](./assets/response-tools-py-figs/response-paths/Pos_4_CdTe_3_Response.png "Position 4 photon path")
+![Position 4 photon path](./assets/response-tools-figs/response-paths/Pos_4_CdTe_3_Response.png "Position 4 photon path")
 
 ### Position 5 (HXR, CdTe1)
 
@@ -75,34 +134,17 @@ For a historical description of the response files below (old versions, etc.), s
 - **Pixelated attenuator:** ""
 - **Detector response:** ""
 
-![Position 5 photon path](./assets/response-tools-py-figs/response-paths/Pos_5_CdTe_1_Response.png "Position 5 photon path")
+![Position 5 photon path](./assets/response-tools-figs/response-paths/Pos_5_CdTe_1_Response.png "Position 5 photon path")
 
 ### Position 6 (Timepix)
 
 - N/A
 
-## Software and software help
-
-The repository may make use of multiple languages and each folder like `response-tools-???` hosts code in their respective languages relating to response tools. These include:
-
-- [`response-tools-py`](response-tools-py/) which contains Python code to download and read in the response files.
-  - The README in this directory also contains Python specific information & help
-
 ## Example code
 
 There are a few existing example scripts showing how to use a lot of tools in the repository.
 
-**Python examples:** Python example scripts can be found in response-tools-py [examples](response-tools-py/examples/) folder which has an associated [README file](response-tools-py/examples/README.md).
-
-## Repository Aim
-
-This repository will contain all code that proves useful in loading FOXSI response elements (at least from FOXSI-4). Several languages might be used and so the top level will be to contain language specific packages.
-
-For example, the `response-tools-py` folder will be a Python package containing all the necessary `.py` files. If other languages are to be included, like `IDL`, then a folder called `response-tools-idl` should be created and used to contain all the `IDL`-ness of the repository. This standard can be applied to the inclusion of other languages (e.g., C++ as `response-tools-cpp`, Shell code as `response-tools-shell`, etc.).
-
-Every `response-tools-<?>` folder should have an "examples" and a "tests" folder. The "examples" folder is a great place to include scripts that show how some of the code in the repository works and the "tests" folder is a fantastic place to put code that tests the tools that have been created.
-
-Additionally, there is also an "examples" and "tests" folder in the top level of the repository so there is a place for anything that fits these folders that spans across code from multiple languages.
+**Python examples:** Python example scripts can be found in response-tools [examples](response-tools/examples/) folder which has an associated [README file](response-tools/examples/README.md).
 
 ## Contributing to the repository
 
@@ -112,14 +154,6 @@ In order to contribute, we ask that you first create your own fork of the reposi
 
 **Note:** We aim to _never_ `push` from a local machine to this repository directly. If this happens then it can be very difficult for other contributers to understand what changes are being made and how it affects their own PRs. _If the repository is pushed to directly, in order to help track changes and make them visible to other contributers, the repository will be reverted back to it's state before the push and the undone changes will be asked to be proposed via a PR to then be merged._
 
-## The `external` directory
-
-The `external` directory is a place for software external to the repository. There is no guarentee that it will follow any specific coding language and so it would perhaps not be ideal to place it in a specific coding language `tools` directory. Some may be [`git` submodules](https://git-scm.com/book/en/v2/Git-Tools-Submodules).
-
-To clone submodules when cloning the main repository, try the following:
-
-- `git clone --recurse-submodules https://github.com/foxsi/response-tools.git`
-
 ## The `response-information` directory
 
 This directory contains context information about the response files that can be shared between a lot of different scripts. E.g., the latest versions of the files the code should be using.
diff --git a/assets/response-tools-py-figs/att-figs/atmospheric-transmissions.png b/assets/response-tools-figs/att-figs/atmospheric-transmissions.png
similarity index 100%
rename from assets/response-tools-py-figs/att-figs/atmospheric-transmissions.png
rename to assets/response-tools-figs/att-figs/atmospheric-transmissions.png
diff --git a/assets/response-tools-py-figs/att-figs/model-transmissions.png b/assets/response-tools-figs/att-figs/model-transmissions.png
similarity index 100%
rename from assets/response-tools-py-figs/att-figs/model-transmissions.png
rename to assets/response-tools-figs/att-figs/model-transmissions.png
diff --git a/assets/response-tools-py-figs/att-figs/transmissions.png b/assets/response-tools-figs/att-figs/transmissions.png
similarity index 100%
rename from assets/response-tools-py-figs/att-figs/transmissions.png
rename to assets/response-tools-figs/att-figs/transmissions.png
diff --git a/assets/response-tools-py-figs/det-resp-figs/cdte-response-matrix.png b/assets/response-tools-figs/det-resp-figs/cdte-response-matrix.png
similarity index 100%
rename from assets/response-tools-py-figs/det-resp-figs/cdte-response-matrix.png
rename to assets/response-tools-figs/det-resp-figs/cdte-response-matrix.png
diff --git a/assets/response-tools-py-figs/det-resp-figs/cmos-response-matrices.png b/assets/response-tools-figs/det-resp-figs/cmos-response-matrices.png
similarity index 100%
rename from assets/response-tools-py-figs/det-resp-figs/cmos-response-matrices.png
rename to assets/response-tools-figs/det-resp-figs/cmos-response-matrices.png
diff --git a/assets/response-tools-py-figs/eff-area-figs/cmos-optics-resp.png b/assets/response-tools-figs/eff-area-figs/cmos-optics-resp.png
similarity index 100%
rename from assets/response-tools-py-figs/eff-area-figs/cmos-optics-resp.png
rename to assets/response-tools-figs/eff-area-figs/cmos-optics-resp.png
diff --git a/assets/response-tools-py-figs/eff-area-figs/heritage-and-msfc-optics.png b/assets/response-tools-figs/eff-area-figs/heritage-and-msfc-optics.png
similarity index 100%
rename from assets/response-tools-py-figs/eff-area-figs/heritage-and-msfc-optics.png
rename to assets/response-tools-figs/eff-area-figs/heritage-and-msfc-optics.png
diff --git a/assets/response-tools-py-figs/quantum-eff-figs/cmos-qunatum-eff.png b/assets/response-tools-figs/quantum-eff-figs/cmos-qunatum-eff.png
similarity index 100%
rename from assets/response-tools-py-figs/quantum-eff-figs/cmos-qunatum-eff.png
rename to assets/response-tools-figs/quantum-eff-figs/cmos-qunatum-eff.png
diff --git a/assets/response-tools-py-figs/response-figs/response-chain.png b/assets/response-tools-figs/response-figs/response-chain.png
similarity index 100%
rename from assets/response-tools-py-figs/response-figs/response-chain.png
rename to assets/response-tools-figs/response-figs/response-chain.png
diff --git a/assets/response-tools-py-figs/response-figs/response-hit-combinations.png b/assets/response-tools-figs/response-figs/response-hit-combinations.png
similarity index 100%
rename from assets/response-tools-py-figs/response-figs/response-hit-combinations.png
rename to assets/response-tools-figs/response-figs/response-hit-combinations.png
diff --git a/assets/response-tools-py-figs/response-paths/Pos_0_CMOS_1_Response.png b/assets/response-tools-figs/response-paths/Pos_0_CMOS_1_Response.png
similarity index 100%
rename from assets/response-tools-py-figs/response-paths/Pos_0_CMOS_1_Response.png
rename to assets/response-tools-figs/response-paths/Pos_0_CMOS_1_Response.png
diff --git a/assets/response-tools-py-figs/response-paths/Pos_1_CMOS_2_Response.png b/assets/response-tools-figs/response-paths/Pos_1_CMOS_2_Response.png
similarity index 100%
rename from assets/response-tools-py-figs/response-paths/Pos_1_CMOS_2_Response.png
rename to assets/response-tools-figs/response-paths/Pos_1_CMOS_2_Response.png
diff --git a/assets/response-tools-py-figs/response-paths/Pos_2_CdTe_4_Response.png b/assets/response-tools-figs/response-paths/Pos_2_CdTe_4_Response.png
similarity index 100%
rename from assets/response-tools-py-figs/response-paths/Pos_2_CdTe_4_Response.png
rename to assets/response-tools-figs/response-paths/Pos_2_CdTe_4_Response.png
diff --git a/assets/response-tools-py-figs/response-paths/Pos_3_CdTe_2_Response.png b/assets/response-tools-figs/response-paths/Pos_3_CdTe_2_Response.png
similarity index 100%
rename from assets/response-tools-py-figs/response-paths/Pos_3_CdTe_2_Response.png
rename to assets/response-tools-figs/response-paths/Pos_3_CdTe_2_Response.png
diff --git a/assets/response-tools-py-figs/response-paths/Pos_4_CdTe_3_Response.png b/assets/response-tools-figs/response-paths/Pos_4_CdTe_3_Response.png
similarity index 100%
rename from assets/response-tools-py-figs/response-paths/Pos_4_CdTe_3_Response.png
rename to assets/response-tools-figs/response-paths/Pos_4_CdTe_3_Response.png
diff --git a/assets/response-tools-py-figs/response-paths/Pos_5_CdTe_1_Response.png b/assets/response-tools-figs/response-paths/Pos_5_CdTe_1_Response.png
similarity index 100%
rename from assets/response-tools-py-figs/response-paths/Pos_5_CdTe_1_Response.png
rename to assets/response-tools-figs/response-paths/Pos_5_CdTe_1_Response.png
diff --git a/examples/README.md b/examples/README.md
index 62925cd..5637cc9 100644
--- a/examples/README.md
+++ b/examples/README.md
@@ -1,3 +1,13 @@
-# `cdte-tools`
+# `response-tools` Examples
 
-Please put examples spanning multiple languages in here.
+This module contains functions to handle the visualization of the various response data products.
+
+## 1. Data handling examples
+
+- [using_the_functions.py](./functions_and_outputs.py)
+  - Shows how to use the functions and their outputs in the package.
+
+## 2. Plotting examples
+
+- [plot_arf_rmf_drm.py](./plot_arf_rmf_drm.py)
+  - Shows how to obtain and plot the ARF, RMF, and DRM for Telescope 2.
diff --git a/response-tools-py/examples/functions_and_outputs.py b/examples/functions_and_outputs.py
similarity index 98%
rename from response-tools-py/examples/functions_and_outputs.py
rename to examples/functions_and_outputs.py
index ba90ea6..1a42225 100644
--- a/response-tools-py/examples/functions_and_outputs.py
+++ b/examples/functions_and_outputs.py
@@ -46,8 +46,8 @@
 (DRM).
 """
 
-import response_tools_py.responses as responses
-import response_tools_py.telescope_parts as telescope_parts
+import response_tools.responses as responses
+import response_tools.telescope_parts as telescope_parts
 
 """
 Let's look at the `foxsi4_position2_optics` function since we mentioned
@@ -57,7 +57,7 @@
 ```
 >>> help(telescope_parts.foxsi4_position2_optics)
 Help on function foxsi4_position2_optics in module 
-response_tools_py.telescope_parts:
+response_tools.telescope_parts:
 
 foxsi4_position2_optics(mid_energies, off_axis_angle)
     Position 2 MSFC heritage X-7 optic effective areas.
diff --git a/response-tools-py/examples/plot_arf_rmf_drm.py b/examples/plot_arf_rmf_drm.py
similarity index 98%
rename from response-tools-py/examples/plot_arf_rmf_drm.py
rename to examples/plot_arf_rmf_drm.py
index 3bee15f..6bbd1e1 100644
--- a/response-tools-py/examples/plot_arf_rmf_drm.py
+++ b/examples/plot_arf_rmf_drm.py
@@ -23,7 +23,7 @@
 import matplotlib.pyplot as plt
 import numpy as np
 
-import response_tools_py.responses as responses
+import response_tools.responses as responses
 
 fig = plt.figure(figsize=(18, 5))
 gs = gridspec.GridSpec(1, 3)
diff --git a/response-information/README.md b/response-information/README.md
index a3af295..afee3cc 100644
--- a/response-information/README.md
+++ b/response-information/README.md
@@ -85,7 +85,7 @@ When a newer version of a file comes along, the text describing any old file can
   - _Origin:_ Uses [see Morgan]
   - _Version description:_
 - **Optics:** `"FOXSI4_Module_MSFC_HiRes_EA_with_models_v1.txt"` (Marshall 2-shell X09/FM2)
-  - _Original file:_ `"FOXSI4_Effective_Areas_2025April11.txt"`[*](../response-tools-py/response_tools_py/aux)
+  - _Original file:_ `"FOXSI4_Effective_Areas_2025April11.txt"`[*](../response_tools/aux)
   - _Origin:_ `FOXSI`>`FOXSI-4 Flight 36.370`>`Flight Science`>`Response Work`>`Optics`>`msfc-hi-res`>`model-data-comparison`
   - _Version description:_
 - **Filter:** `"thin_mylar_p3_p5_theoretical_v1.csv"`
diff --git a/response-tools-py/README.md b/response-tools-py/README.md
deleted file mode 100644
index 551d482..0000000
--- a/response-tools-py/README.md
+++ /dev/null
@@ -1,58 +0,0 @@
-# `response-tools-py`
-
-A repository to software tools for the response system used for the FOXSI mission in Python (woo, Python). We will include helpful loader functions for the all the response elements and download clients to obtain the files.
-
-The tools being developed in Python should be placed in the `response_tools_py` folder (note the underscores and not dashes).
-
-There is an "examples" and a "tests" folder. The "examples" folder is a great place to include scripts that show how some of the code in the repository works and the "tests" folder is a fantastic place to put code that tests the tools that have been created.
-
-More information will be placed here with regards as to how this package is recommended to be used.
-
-## Install tips
-
-In order to work with some preliminary data, it is instructive to set up a virtual environment (more information below) and install necessary packages. The suggested way to do this is to create an environment with:
-
-- `conda create -n response-tools-env python==3.12`
-
-(recommend using Python 3.12 just because this has been proven to be stable with the software).
-
-The Python code can then be installed with:
-
-- `pip install -e .`
-
-while in this directory (the `response-tools-py` directory) containing the `setup.py` file.
-
-Any time the code is updated and, say, you pull it from Github, make sure to perform the `pip install -e .` line from above. This ensures any new changes, updated versions, etc. take effect.
-
-## Namespace
-
-The base `response_tools_py` namespace includes:
-
-- `~response_tools_py.contextResponseInfo`
-  - The information stored in a [YAML file](../response-information/info.yaml) that includes file version information (e.g., the current versions to be used) and can be accessed like a Python dictionary.
-- `~response_tools_py.__version__`
-  - The version of the code determined by `response-tools-py/setup.py`.
-    - `0.0.1`: First rendition of the code.
-
-### Examples using the namespace
-
-```python
-# importing the module is as easy as:
->>> import response_tools_py
-
-# can access the YAML file contents in the code as 
->>> print(response_tools_py.contextResponseInfo[""])
-...
-
-# then accessing the version as:
->>> print(response_tools_py.__version__)
-'0.0.1'
-```
-
-## Useful Python tips
-
-Some useful things to keep in mind while using the Python code in this repository.
-
-### Virtual environments
-
-It might be a good idea to look into ([conda](https://conda.io/projects/conda/en/latest/user-guide/tasks/manage-environments.html)) virtual environments if you are not familiar, this includes looking into them yourself or getting in touch with someone to help explain. This requires downloading [miniconda](https://docs.anaconda.com/miniconda/install/#quick-command-line-install) (or just conda if you prefer).
diff --git a/response-tools-py/examples/README.md b/response-tools-py/examples/README.md
deleted file mode 100644
index 82e351f..0000000
--- a/response-tools-py/examples/README.md
+++ /dev/null
@@ -1,13 +0,0 @@
-# `response-tools-py` Examples
-
-This module contains functions to handle the visualization of the various response data products.
-
-## 1. Data handling examples
-
-- [using_the_functions.py](./functions_and_outputs.py)
-  - Shows how to use the functions and their outputs in the package.
-
-## 2. Plotting examples
-
-- [plot_arf_rmf_drm.py](./plot_arf_rmf_drm.py)
-  - Shows how to obtain and plot the ARF, RMF, and DRM for Telescope 2.
diff --git a/response-tools-py/response_tools_py/version.py b/response-tools-py/response_tools_py/version.py
deleted file mode 100644
index ef37232..0000000
--- a/response-tools-py/response_tools_py/version.py
+++ /dev/null
@@ -1,3 +0,0 @@
-import importlib.metadata
-
-__version__ = importlib.metadata.version("response_tools_py")
diff --git a/response-tools-py/response_tools_py/__init__.py b/response_tools/__init__.py
similarity index 57%
rename from response-tools-py/response_tools_py/__init__.py
rename to response_tools/__init__.py
index 3f5cdf0..fea117f 100644
--- a/response-tools-py/response_tools_py/__init__.py
+++ b/response_tools/__init__.py
@@ -1,34 +1,34 @@
 """
-Top level of the `response-tools-py` package.
+Top level of the `response-tools` package.
 
 Namespace
 ---------
 
-The base `response_tools_py` namespace includes:
+The base `response_tools` namespace includes:
 
-`~response_tools_py.contextResponseInfo`
+`~response_tools.contextResponseInfo`
     The information stored in a YAML file that includes information from
     the flight.
 
-`~response_tools_py.__version__`
+`~response_tools.__version__`
     The current version of the code as stated in the setup script.
 
 Examples
 --------
 # importing the module is as easy as:
->>> import response_tools_py
+>>> import response_tools
 
 # then accessing, e.g., the context information for the flight/flare
->>> print(response_tools_py.contextResponseInfo)
+>>> print(response_tools.contextResponseInfo)
 ...
 
 # that can be accessed like a native Python dictionary
->>> print(response_tools_py.contextResponseInfo[""])
+>>> print(response_tools.contextResponseInfo["..."])
 ...
 """
 
 from .version import __version__
-from response_tools_py.io.load_yaml import load_response_context
+from response_tools.io.load_yaml import load_response_context
 
 # for global context info
 contextResponseInfo = load_response_context()
diff --git a/response-tools-py/response_tools_py/attenuation.py b/response_tools/attenuation.py
similarity index 99%
rename from response-tools-py/response_tools_py/attenuation.py
rename to response_tools/attenuation.py
index 85b5c7c..076f259 100644
--- a/response-tools-py/response_tools_py/attenuation.py
+++ b/response_tools/attenuation.py
@@ -13,11 +13,11 @@
 import pandas
 import scipy
 
-from response_tools_py.util import BaseOutput, native_resolution
+from response_tools.util import BaseOutput, native_resolution
 
-ATT_PATH = os.path.join(pathlib.Path(__file__).parent, "..", "..", "response-information", "attenuation-data")
-ATM_PATH = os.path.join(pathlib.Path(__file__).parent, "..", "..", "response-information", "atmospheric-data")
-ASSETS_PATH = os.path.join(pathlib.Path(__file__).parent, "..", "..", "assets", "response-tools-py-figs", "att-figs")
+ATT_PATH = os.path.join(pathlib.Path(__file__).parent, "..", "response-information", "attenuation-data")
+ATM_PATH = os.path.join(pathlib.Path(__file__).parent, "..", "response-information", "atmospheric-data")
+ASSETS_PATH = os.path.join(pathlib.Path(__file__).parent, "..", "assets", "response-tools-figs", "att-figs")
 
 @dataclass
 class AttOutput(BaseOutput):
diff --git a/response-tools-py/response_tools_py/aux/README.md b/response_tools/aux/README.md
similarity index 97%
rename from response-tools-py/response_tools_py/aux/README.md
rename to response_tools/aux/README.md
index b4dc63b..15adf54 100644
--- a/response-tools-py/response_tools_py/aux/README.md
+++ b/response_tools/aux/README.md
@@ -1,4 +1,4 @@
-# `response_tools_py.aux`
+# `response_tools.aux`
 
 A folder to contain any extra code/processes used for, say, formatting the files we have. Somewhere to document this stuff.
 
diff --git a/response-tools-py/response_tools_py/aux/format-msfc-eff-areas.py b/response_tools/aux/format-msfc-eff-areas.py
similarity index 100%
rename from response-tools-py/response_tools_py/aux/format-msfc-eff-areas.py
rename to response_tools/aux/format-msfc-eff-areas.py
diff --git a/response-tools-py/response_tools_py/detector_response.py b/response_tools/detector_response.py
similarity index 99%
rename from response-tools-py/response_tools_py/detector_response.py
rename to response_tools/detector_response.py
index a76398d..60d7486 100644
--- a/response-tools-py/response_tools_py/detector_response.py
+++ b/response_tools/detector_response.py
@@ -11,10 +11,10 @@
 import astropy.units as u
 import numpy as np
 
-from response_tools_py.util import BaseOutput
+from response_tools.util import BaseOutput
 
-DET_RESP_PATH = os.path.join(pathlib.Path(__file__).parent, "..", "..", "response-information", "detector-response-data")
-ASSETS_PATH = os.path.join(pathlib.Path(__file__).parent, "..", "..", "assets", "response-tools-py-figs", "det-resp-figs")
+DET_RESP_PATH = os.path.join(pathlib.Path(__file__).parent, "..", "response-information", "detector-response-data")
+ASSETS_PATH = os.path.join(pathlib.Path(__file__).parent, "..", "assets", "response-tools-figs", "det-resp-figs")
 
 @dataclass
 class DetectorResponseOutput(BaseOutput):
diff --git a/response-tools-py/response_tools_py/effective_area.py b/response_tools/effective_area.py
similarity index 99%
rename from response-tools-py/response_tools_py/effective_area.py
rename to response_tools/effective_area.py
index c33c48f..d9cd4b2 100644
--- a/response-tools-py/response_tools_py/effective_area.py
+++ b/response_tools/effective_area.py
@@ -14,10 +14,10 @@
 from scipy.interpolate import CloughTocher2DInterpolator
 import pandas
 
-from response_tools_py.util import BaseOutput, native_resolution
+from response_tools.util import BaseOutput, native_resolution
 
-EFF_PATH = os.path.join(pathlib.Path(__file__).parent, "..", "..", "response-information", "effective-area-data")
-ASSETS_PATH = os.path.join(pathlib.Path(__file__).parent, "..", "..", "assets", "response-tools-py-figs", "eff-area-figs")
+EFF_PATH = os.path.join(pathlib.Path(__file__).parent, "..", "response-information", "effective-area-data")
+ASSETS_PATH = os.path.join(pathlib.Path(__file__).parent, "..", "assets", "response-tools-figs", "eff-area-figs")
 
 @dataclass
 class EffAreaOutput(BaseOutput):
diff --git a/response-tools-py/response_tools_py/io/__init__.py b/response_tools/io/__init__.py
similarity index 100%
rename from response-tools-py/response_tools_py/io/__init__.py
rename to response_tools/io/__init__.py
diff --git a/response-tools-py/response_tools_py/io/load_yaml.py b/response_tools/io/load_yaml.py
similarity index 89%
rename from response-tools-py/response_tools_py/io/load_yaml.py
rename to response_tools/io/load_yaml.py
index eda0c4c..2bda7ad 100644
--- a/response-tools-py/response_tools_py/io/load_yaml.py
+++ b/response_tools/io/load_yaml.py
@@ -16,4 +16,4 @@ def load_yaml(filename):
 
 def load_response_context():
     """Function to load in the context information from file. """
-    return load_yaml(os.path.join(pathlib.Path(__file__).parent, "..", "..", "..", "response-information", "info.yaml"))
\ No newline at end of file
+    return load_yaml(os.path.join(pathlib.Path(__file__).parent, "..", "..", "response-information", "info.yaml"))
\ No newline at end of file
diff --git a/response-tools-py/response_tools_py/io/tests/README.md b/response_tools/io/tests/README.md
similarity index 56%
rename from response-tools-py/response_tools_py/io/tests/README.md
rename to response_tools/io/tests/README.md
index f5f1f48..5a36b57 100644
--- a/response-tools-py/response_tools_py/io/tests/README.md
+++ b/response_tools/io/tests/README.md
@@ -1,3 +1,3 @@
-# `response-tools-py.io` Tests
+# `response-tools.io` Tests
 
 Please put Python test scripts in here.
diff --git a/response-tools-py/response_tools_py/io/tests/__init__.py b/response_tools/io/tests/__init__.py
similarity index 100%
rename from response-tools-py/response_tools_py/io/tests/__init__.py
rename to response_tools/io/tests/__init__.py
diff --git a/response-tools-py/response_tools_py/ir-plots.py b/response_tools/ir-plots.py
similarity index 99%
rename from response-tools-py/response_tools_py/ir-plots.py
rename to response_tools/ir-plots.py
index 54d7ee1..40ede16 100644
--- a/response-tools-py/response_tools_py/ir-plots.py
+++ b/response_tools/ir-plots.py
@@ -26,7 +26,7 @@
 
 def foxsi4_response_paths(save_assets=False):
     
-    assets_dir = os.path.join(pathlib.Path(__file__).parent, "..", "..", "assets", "response-tools-py-figs", "response-paths")
+    assets_dir = os.path.join(pathlib.Path(__file__).parent, "..", "..", "assets", "response-tools-figs", "response-paths")
     pathlib.Path(assets_dir).mkdir(parents=True, exist_ok=True)
 
     plt.rcParams.update({'font.size': 7})
diff --git a/response-tools-py/response_tools_py/phot_spec.py b/response_tools/phot_spec.py
similarity index 100%
rename from response-tools-py/response_tools_py/phot_spec.py
rename to response_tools/phot_spec.py
diff --git a/response-tools-py/response_tools_py/quantum_efficiency.py b/response_tools/quantum_efficiency.py
similarity index 94%
rename from response-tools-py/response_tools_py/quantum_efficiency.py
rename to response_tools/quantum_efficiency.py
index 819cbc6..7ec72b7 100644
--- a/response-tools-py/response_tools_py/quantum_efficiency.py
+++ b/response_tools/quantum_efficiency.py
@@ -10,10 +10,10 @@
 import astropy.units as u
 import numpy as np
 
-from response_tools_py.util import BaseOutput, native_resolution
+from response_tools.util import BaseOutput, native_resolution
 
-Q_PATH = os.path.join(pathlib.Path(__file__).parent, "..", "..", "response-information", "quantum-efficiency-data")
-ASSETS_PATH = os.path.join(pathlib.Path(__file__).parent, "..", "..", "assets", "response-tools-py-figs", "quantum-eff-figs")
+Q_PATH = os.path.join(pathlib.Path(__file__).parent, "..", "response-information", "quantum-efficiency-data")
+ASSETS_PATH = os.path.join(pathlib.Path(__file__).parent, "..", "assets", "response-tools-figs", "quantum-eff-figs")
 
 @dataclass
 class QuantumEffOutput(BaseOutput):
diff --git a/response-tools-py/response_tools_py/responses.py b/response_tools/responses.py
similarity index 99%
rename from response-tools-py/response_tools_py/responses.py
rename to response_tools/responses.py
index 435f380..28188f9 100644
--- a/response-tools-py/response_tools_py/responses.py
+++ b/response_tools/responses.py
@@ -9,11 +9,11 @@
 import astropy.units as u
 import numpy as np
 
-from response_tools_py.attenuation import att_foxsi4_atmosphere
-import response_tools_py.telescope_parts as tp
-from response_tools_py.util import BaseOutput
+from response_tools.attenuation import att_foxsi4_atmosphere
+import response_tools.telescope_parts as tp
+from response_tools.util import BaseOutput
 
-ASSETS_PATH = os.path.join(pathlib.Path(__file__).parent, "..", "..", "assets", "response-tools-py-figs", "response-figs")
+ASSETS_PATH = os.path.join(pathlib.Path(__file__).parent, "..", "assets", "response-tools-figs", "response-figs")
 
 @dataclass
 class Response1DOutput(BaseOutput):
diff --git a/response-tools-py/response_tools_py/telescope_parts.py b/response_tools/telescope_parts.py
similarity index 97%
rename from response-tools-py/response_tools_py/telescope_parts.py
rename to response_tools/telescope_parts.py
index fb17526..f7513e0 100644
--- a/response-tools-py/response_tools_py/telescope_parts.py
+++ b/response_tools/telescope_parts.py
@@ -5,18 +5,18 @@
 
 import astropy.units as u
 
-from response_tools_py.attenuation import (att_thermal_blanket,
-                                           att_pixelated, 
-                                           att_al_mylar,
-                                           att_uniform_al_cdte
-                                           )
-from response_tools_py.detector_response import (cdte_det_resp, 
-                                                 cmos_det_resp,
-                                                 )
-from response_tools_py.effective_area import (eff_area_msfc_10shell,
-                                              eff_area_msfc_hi_res,
-                                              eff_area_nagoya_hxt,
+from response_tools.attenuation import (att_thermal_blanket,
+                                        att_pixelated, 
+                                        att_al_mylar,
+                                        att_uniform_al_cdte
+                                        )
+from response_tools.detector_response import (cdte_det_resp, 
+                                              cmos_det_resp,
                                               )
+from response_tools.effective_area import (eff_area_msfc_10shell,
+                                           eff_area_msfc_hi_res,
+                                           eff_area_nagoya_hxt,
+                                           )
 
 # position 2
 @u.quantity_input(mid_energies=u.keV)
diff --git a/response-tools-py/response_tools_py/util.py b/response_tools/util.py
similarity index 100%
rename from response-tools-py/response_tools_py/util.py
rename to response_tools/util.py
diff --git a/response_tools/version.py b/response_tools/version.py
new file mode 100644
index 0000000..ee5165f
--- /dev/null
+++ b/response_tools/version.py
@@ -0,0 +1,3 @@
+import importlib.metadata
+
+__version__ = importlib.metadata.version("response_tools")
diff --git a/response-tools-py/setup.py b/setup.py
similarity index 93%
rename from response-tools-py/setup.py
rename to setup.py
index eaabd7f..7c8254c 100644
--- a/response-tools-py/setup.py
+++ b/setup.py
@@ -1,7 +1,7 @@
 import setuptools
 
 setuptools.setup(
-    name="response_tools_py",
+    name="response_tools",
     version="0.0.1",
     description="Software used for the FOXSI response files.",
     url="https://github.com/foxsi/response-tools",
diff --git a/tests/README.md b/tests/README.md
index 87d6817..28d2a44 100644
--- a/tests/README.md
+++ b/tests/README.md
@@ -1,3 +1,3 @@
 # `response-tools`
 
-Please put tests spanning multiple languages in here.
+Please put Python tests in here.