Skip to content

Instantly share code, notes, and snippets.

@pllim
Created October 6, 2017 14:56
Show Gist options
  • Save pllim/ba822cb4d2b37a34185c889a227e1f6a to your computer and use it in GitHub Desktop.
Save pllim/ba822cb4d2b37a34185c889a227e1f6a to your computer and use it in GitHub Desktop.
PySynphot example: Calculation in ABMAG
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"import os\n",
"import numpy as np\n",
"import pysynphot as S"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"# Where you installed the data files\n",
"cdbs = os.environ['PYSYN_CDBS']"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
"# Change this to your source spectrum\n",
"sp = S.FileSpectrum(cdbs + 'calspec/alpha_lyr_stis_008.fits')"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"# Bandpass\n",
"bp = S.FileBandpass(cdbs + 'comp/nonhst/sdss_g_005_syn.fits')"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"# Define an array for optimal binning of your bandpass.\n",
"# This is just an example.\n",
"binset = np.arange(3000, 7000)\n",
"\n",
"# Create an observation to convolve them\n",
"obs = S.Observation(sp, bp, binset=binset)"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"-0.09533066034714466"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Effective stimulus in ABMAG\n",
"obs.effstim('abmag')"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.5.4"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment