This analysis tool runs a bow-tie analysis (Van Allen et al. 1974) for the energy channels of a generic particle instrument. The input is a csv table of channel responses indexed by the incident energy. The results of analysis are the geometric factor (with errors) and the effective energy of the channel.
Tested in Ubuntu 20.04.6 LTS with Python version 3.12.8.
- This tool requires a recent Python (>=3.10) installation. Following SunPy's approach, we recommend installing Python via miniforge (click for instructions).
- Download this file and extract to a folder of your choice (or clone the repository https://github.com/spearhead-he/bowtie if you know how to use
git). - Open a terminal or the miniforge prompt and move to the directory where the code is.
- (Optional, but recommended) Create a new virtual environment (e.g.,
conda create --name bowtie, orpython -m venv venv_bowtie_toolif you don't use miniforge/conda) and activate it (e.g.,conda activate bowtie, orsource venv_bowtie_tool/bin/activateif you don't use miniforge/conda). - Install the Python dependencies from the requirements.txt file with
pip install -r requirements.txt - Open the Jupyter Notebook by running
jupyter-lab bowtie_example.ipynb
The Notebook is a simple example that also acts as a tutorial to teach the user on how to run the bow-tie analysis with this tool.
The tool operates with two main classes, which are called Bowtie and Spectra. Bowtie stores response functions and contains the methods to run bow-tie analysis, while Spectra contains the information of the spectral indices and the amount of different spectra that are used in the bow-tie calculation.
The Bowtie class contains the data that the bow-tie analysis is applied on, and the energy range to be considered in the calculations. Its methods make running analysis easy and straightforward.
Methods:
set_energy_range(energy_min, energy_max):
energy_min : {float} The minimum energy in MeV to consider
energy_max : {float} See energy_min.
bowtie_analysis(channel, spectra, plot):
channel : {str} The channel name as it appears in the csv table.
spectra : {Spectra} The Spectra class object, introduced in this package. Contains the
spectral indices and the power law spectra used in the bow-tie analysis.
plot : {bool} A boolean switch to produce a plot visualizing the analysis.
bowtie_analysis_full_stack(spectra, plot):
A wrapper for bowtie_analysis(). Runs the analysis on all channels that appear in the input file.
spectra : {Spectra} See bowtie_analysis().
plot : {bool} See bowtie_analysis()
The Spectra class contains the range of spectra that are applied on the response function to run bow-tie analysis.
Methods:
set_spectral_indices(gamma_min, gamma_max):
gamma_min : {float} The minimum spectral index to consider in the calculation.
gamma_max : {float} See gamma_min.
produce_power_law_spectra(response_df):
response_df : {pandas.DataFrame} The input csv table read in to a pandas DataFrame. Contains
the channel responses as a function of incident energy.
Contributions to this tool are very much welcome and encouraged! Contributions can take the form of issues to report bugs and request new features or pull requests to submit new code.
If you don't have a GitHub account, you can sign-up for free here, or you can also reach out to us with feedback by sending an email to [email protected].
This tool is developed within the SPEARHEAD (SPEcification, Analysis & Re-calibration of High Energy pArticle Data) project. SPEARHEAD has received funding from the European Union’s Horizon Europe programme under grant agreement No 101135044.
The tool reflects only the authors’ view and the European Commission is not responsible for any use that may be made of the information it contains.