{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "Spectra\n", "---\n", "\n", "Download all the Jupyter notebooks from: https://github.com/HeloiseS/hoki/tree/master/tutorials" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Initial Imports" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from hoki import load\n", "from hoki.spec import dopcor\n", "import pandas as pd\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "from sklearn.preprocessing import MinMaxScaler\n", "\n", "%matplotlib inline\n", "plt.style.use('tuto.mplstyle')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### In this tutorial you will:\n", "\n", "- Load BPASS synthetic spectra through `hoki`\n", "- Manipulated an plot these synthetic spectra\n", "- Compare them to observational data and find best matches. \n", "- Get more familiar with `pandas` if you are not already." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Loading Spectra in `hoki`\n", "\n", "The synthetic spectra returned by BPASS are given in ASCII files, grouped by metallicity (`zXXX`) and whether the population is made up of just single stars (`-sin-`) or also contains binary stars (`-bin`).\n", "Each file contains the wavelenght bin in its first column and the following 51 columns are the synthetic spectra of the population at each BPASS time bins. The wavelength ranges from 1 to 100,000 Angstrom.\n", "\n", "We can load each file into a `pandas.DataFrame` using the `hoki.load.model_output()` function.\n", "**This may take a few seconds.**\n", "For the purposes of this tutorial we are only going to load two sets of synthetic spectra: solar metallicity and a tenth solar - both including binary stars. " ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "ss_z002 = load.model_output('./data/spectra/spectra-bin-imf135_300.z002.dat')\n", "ss_z020 = load.model_output('./data/spectra/spectra-bin-imf135_300.z020.dat')" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", " | WL | \n", "6.0 | \n", "6.1 | \n", "6.2 | \n", "6.3 | \n", "6.4 | \n", "6.5 | \n", "6.6 | \n", "6.7 | \n", "6.8 | \n", "... | \n", "10.1 | \n", "10.2 | \n", "10.3 | \n", "10.4 | \n", "10.5 | \n", "10.6 | \n", "10.7 | \n", "10.8 | \n", "10.9 | \n", "11.0 | \n", "
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
0 | \n", "1.0 | \n", "2.745775e-40 | \n", "2.702434e-40 | \n", "3.047202e-40 | \n", "4.195247e-38 | \n", "7.912627e-38 | \n", "2.884525e-38 | \n", "9.237552e-39 | \n", "1.976199e-38 | \n", "5.478071e-38 | \n", "... | \n", "3.143051e-40 | \n", "4.732241e-40 | \n", "2.407628e-40 | \n", "1.854254e-40 | \n", "1.770071e-40 | \n", "3.505892e-40 | \n", "1.656219e-40 | \n", "1.040945e-40 | \n", "2.193011e-40 | \n", "3.402488e-41 | \n", "
1 | \n", "2.0 | \n", "2.745775e-40 | \n", "2.702434e-40 | \n", "3.047202e-40 | \n", "4.195247e-38 | \n", "7.912627e-38 | \n", "2.884525e-38 | \n", "9.237552e-39 | \n", "1.976199e-38 | \n", "5.478071e-38 | \n", "... | \n", "3.143051e-40 | \n", "4.732241e-40 | \n", "2.407628e-40 | \n", "1.854254e-40 | \n", "1.770071e-40 | \n", "3.505892e-40 | \n", "1.656219e-40 | \n", "1.040945e-40 | \n", "2.193011e-40 | \n", "3.402488e-41 | \n", "
2 | \n", "3.0 | \n", "2.745775e-40 | \n", "2.702434e-40 | \n", "3.047202e-40 | \n", "4.195247e-38 | \n", "7.912627e-38 | \n", "2.884525e-38 | \n", "9.237552e-39 | \n", "1.976199e-38 | \n", "5.478071e-38 | \n", "... | \n", "3.143051e-40 | \n", "4.732241e-40 | \n", "2.407628e-40 | \n", "1.854254e-40 | \n", "1.770071e-40 | \n", "3.505892e-40 | \n", "1.656219e-40 | \n", "1.040945e-40 | \n", "2.193011e-40 | \n", "3.402488e-41 | \n", "
3 | \n", "4.0 | \n", "2.745775e-40 | \n", "2.702434e-40 | \n", "3.047202e-40 | \n", "4.195247e-38 | \n", "7.912627e-38 | \n", "2.884525e-38 | \n", "9.237552e-39 | \n", "1.976199e-38 | \n", "5.478071e-38 | \n", "... | \n", "3.143051e-40 | \n", "4.732241e-40 | \n", "2.407628e-40 | \n", "1.854254e-40 | \n", "1.770071e-40 | \n", "3.505892e-40 | \n", "1.656219e-40 | \n", "1.040945e-40 | \n", "2.193011e-40 | \n", "3.402488e-41 | \n", "
4 | \n", "5.0 | \n", "2.745775e-40 | \n", "2.702434e-40 | \n", "3.047202e-40 | \n", "4.195247e-38 | \n", "7.912627e-38 | \n", "2.884525e-38 | \n", "9.237552e-39 | \n", "1.976199e-38 | \n", "5.478071e-38 | \n", "... | \n", "3.143051e-40 | \n", "4.732241e-40 | \n", "2.407628e-40 | \n", "1.854254e-40 | \n", "1.770071e-40 | \n", "3.505892e-40 | \n", "1.656219e-40 | \n", "1.040945e-40 | \n", "2.193011e-40 | \n", "3.402488e-41 | \n", "
5 rows × 52 columns
\n", "\n", " | WL | \n", "Counts | \n", "
---|---|---|
0 | \n", "4036.82 | \n", "5437.012 | \n", "
1 | \n", "4038.36 | \n", "5506.620 | \n", "
2 | \n", "4039.90 | \n", "5747.366 | \n", "
3 | \n", "4041.44 | \n", "5748.782 | \n", "
4 | \n", "4042.99 | \n", "5761.701 | \n", "