Created
July 10, 2019 15:41
-
-
Save yuasatakayuki/830053ad2e0dc3764c0ded823d55d4b5 to your computer and use it in GitHub Desktop.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
{ | |
"cells": [ | |
{ | |
"cell_type": "code", | |
"execution_count": 29, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"%matplotlib inline\n", | |
"import astropy.io.fits as af\n", | |
"import matplotlib.pyplot as plt\n", | |
"import numpy as np" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 30, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"fits = af.open(\"20171124_184608.fits\")\n", | |
"event_list = fits[1].data\n", | |
"time = event_list[\"unixTime\"]-event_list[\"unixTime\"][0]\n", | |
"energy = np.array(event_list[\"energy\"])" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 32, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAD8CAYAAABQFVIjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAFuVJREFUeJzt3X+w3XWd3/Hnq2FhratLkDQTCWmijc6A00bJIM6uDtUVAjqinR0bZkeiskYrdNZ2Z7ah/oF1ywz+qltmLYhrKnRckBWVjMSykbFrO1OQsFJ+CXL5YUkmkgiudNYdV/TdP87n4pdwb+6P87333HPzfMycud/z/v76fPO997zy+X6+55xUFZIkDesfjLoBkqTlwUCRJPXCQJEk9cJAkST1wkCRJPXCQJEk9cJAkST1wkCRJPXCQJEk9eKYUTdgvk488cRav379qJshSWPlzjvv/FFVrVqIbY9toKxfv569e/eOuhmSNFaS/GChtu0lL0lSLwwUSVIvDBRJUi8MFElSLwwUSVIvDBRJUi8MFElSLwwUSVIvDBRJUi8MFEkaofU7bmb9jptH3YxeGCiSpF4YKJKkXswYKEl2JjmY5N5O7UtJ7mqPx5Lc1errk/xdZ95VnXVOS3JPkokkVyRJq5+QZE+Sh9rPlQtxoJK01I375a/Z9FC+AGzpFqrqX1bVpqraBNwIfKUz++HJeVX1gU79SuB9wMb2mNzmDuDWqtoI3NqeS5LGzIyBUlXfBp6aal7rZbwTuO5I20iyBnhxVd1WVQVcC7y9zT4PuKZNX9OpS5LGyLBjKK8Hnqiqhzq1DUm+m+Svkry+1U4C9nWW2ddqAKur6kCb/iGwesg2SZJGYNgv2Dqf5/ZODgDrqurJJKcBX0ty6mw3VlWVpKabn2Q7sB1g3bp182yyJGkhzLuHkuQY4F8AX5qsVdXPqurJNn0n8DDwCmA/sLaz+tpWA3iiXRKbvDR2cLp9VtXVVbW5qjavWrUg32ApSZqnYS55/Q7wQFU9eykryaokK9r0yxgMvj/SLmk9neSMNu5yAXBTW20XsK1Nb+vUJUljZDa3DV8H/G/glUn2JbmwzdrK8wfj3wDc3W4j/jLwgaqaHND/IPBnwASDnss3Wv1y4M1JHmIQUpcPcTySpBGZcQylqs6fpv7uKWo3MriNeKrl9wKvmqL+JPCmmdohSVrafKe8JKkXBookqRcGiiSpFwaKJKkXBookjcg4fxDkVAwUSVIvDBRJUi8MFElSLwwUSVIvDBRJUi8MFElSLwwUSVIvDBRJUi8MFElSLwwUSVIvDBRJUi8MFElSLwwUSVIvDBRJWkTL7ROGuwwUSVIvZgyUJDuTHExyb6f2kST7k9zVHud25l2SZCLJg0nO7tS3tNpEkh2d+oYkt7f6l5Ic2+cBSpIWx2x6KF8AtkxR/3RVbWqP3QBJTgG2Aqe2df5LkhVJVgCfAc4BTgHOb8sCfKxt658APwYuHOaAJEmjMWOgVNW3gadmub3zgOur6mdV9SgwAZzeHhNV9UhV/T1wPXBekgBvBL7c1r8GePscj0GStAQMM4ZycZK72yWxla12EvB4Z5l9rTZd/SXA31TVM4fVp5Rke5K9SfYeOnRoiKZLkvo230C5Eng5sAk4AHyqtxYdQVVdXVWbq2rzqlWrFmOXkqRZOmY+K1XVE5PTST4HfL093Q+c3Fl0basxTf1J4Pgkx7ReSnd5SdIYmVcPJcmaztN3AJN3gO0CtiY5LskGYCPwHeAOYGO7o+tYBgP3u6qqgG8Bv9vW3wbcNJ82SZJGa8YeSpLrgDOBE5PsAy4FzkyyCSjgMeD9AFV1X5IbgPuBZ4CLquoXbTsXA7cAK4CdVXVf28W/A65P8h+B7wKf7+3oJEmLZsZAqarzpyhP+6JfVZcBl01R3w3snqL+CIO7wCRp2VrO75Cf5DvlJUm9MFAkSb0wUCRJvTBQJEm9MFAkSb0wUCRJvTBQJEm9MFAkSb0wUCRJvTBQJEm9MFAkaYGs33HzUfGRK5MMFElSLwwUSVIvDBRJUi8MFElSLwwUSVIvDBRJUi8MFElSLwwUSVIvZgyUJDuTHExyb6f2iSQPJLk7yVeTHN/q65P8XZK72uOqzjqnJbknyUSSK5Kk1U9IsifJQ+3nyoU4UEnSwppND+ULwJbDanuAV1XVPwW+D1zSmfdwVW1qjw906lcC7wM2tsfkNncAt1bVRuDW9lySNGZmDJSq+jbw1GG1v6yqZ9rT24C1R9pGkjXAi6vqtqoq4Frg7W32ecA1bfqaTl2SNEb6GEN5L/CNzvMNSb6b5K+SvL7VTgL2dZbZ12oAq6vqQJv+IbC6hzZJ0kgdTZ/hNemYYVZO8mHgGeCLrXQAWFdVTyY5DfhaklNnu72qqiR1hP1tB7YDrFu3bv4NlyT1bt49lCTvBt4K/F67jEVV/ayqnmzTdwIPA68A9vPcy2JrWw3giXZJbPLS2MHp9llVV1fV5qravGrVqvk2XZK0AOYVKEm2AH8EvK2qftqpr0qyok2/jMHg+yPtktbTSc5od3ddANzUVtsFbGvT2zp1SdIYmfGSV5LrgDOBE5PsAy5lcFfXccCedvfvbe2OrjcAH03yc+CXwAeqanJA/4MM7hh7AYMxl8lxl8uBG5JcCPwAeGcvRyZJWlQzBkpVnT9F+fPTLHsjcOM08/YCr5qi/iTwppnaIUlL3fodN/PY5W8ZdTNGxnfKS5J6YaBIknphoEiSemGgSJJ6YaBIknphoEiSemGgSJJ6MdRneUmSjs4PgpyKPRRJUi8MFElSLwwUSVIvDBRJUi8MFElSLwwUSVIvDBRJUi8MFElSLwwUSVIvDBRJmgffHf98BookqRcGiiSpF7MKlCQ7kxxMcm+ndkKSPUkeaj9XtnqSXJFkIsndSV7TWWdbW/6hJNs69dOS3NPWuSJJ+jxISdLCm20P5QvAlsNqO4Bbq2ojcGt7DnAOsLE9tgNXwiCAgEuB1wKnA5dOhlBb5n2d9Q7flyRpiZtVoFTVt4GnDiufB1zTpq8B3t6pX1sDtwHHJ1kDnA3sqaqnqurHwB5gS5v34qq6raoKuLazLUnSmBhmDGV1VR1o0z8EVrfpk4DHO8vta7Uj1fdNUX+eJNuT7E2y99ChQ0M0XZLUt14G5VvPovrY1gz7ubqqNlfV5lWrVi307iRJczBMoDzRLlfRfh5s9f3AyZ3l1rbakeprp6hLksbIMIGyC5i8U2sbcFOnfkG72+sM4Cft0tgtwFlJVrbB+LOAW9q8p5Oc0e7uuqCzLUlaUtbvuNk3NU5jVt8pn+Q64EzgxCT7GNytdTlwQ5ILgR8A72yL7wbOBSaAnwLvAaiqp5L8MXBHW+6jVTU50P9BBneSvQD4RntIksbIrAKlqs6fZtabpli2gIum2c5OYOcU9b3Aq2bTFknS0uQ75SVJvTBQJEm9MFAkSb0wUCRJvTBQJGkG3io8OwaKJKkXBookqRcGiiSpFwaKJKkXBookqRcGiiQdgXd3zZ6BIknqhYEiSeqFgSJJ6oWBIknqhYEiSeqFgSJJ6oWBIknqhYEiSYfxvSfzM+9ASfLKJHd1Hk8n+VCSjyTZ36mf21nnkiQTSR5McnanvqXVJpLsGPagJEmL75j5rlhVDwKbAJKsAPYDXwXeA3y6qj7ZXT7JKcBW4FTgpcA3k7yizf4M8GZgH3BHkl1Vdf982yZJWnzzDpTDvAl4uKp+kGS6Zc4Drq+qnwGPJpkATm/zJqrqEYAk17dlDRRJGiN9BcpW4LrO84uTXADsBf6wqn4MnATc1llmX6sBPH5Y/bU9tUuSZs2xk+EMPSif5FjgbcBftNKVwMsZXA47AHxq2H109rU9yd4kew8dOtTXZiVJPejjLq9zgL+uqicAquqJqvpFVf0S+By/uqy1Hzi5s97aVpuu/jxVdXVVba6qzatWreqh6ZKkvvQRKOfTudyVZE1n3juAe9v0LmBrkuOSbAA2At8B7gA2JtnQejtb27KSpDEy1BhKkhcyuDvr/Z3yx5NsAgp4bHJeVd2X5AYGg+3PABdV1S/adi4GbgFWADur6r5h2iVJWnxDBUpV/S3wksNq7zrC8pcBl01R3w3sHqYtkjQfkwPxj13+lhG3ZPz5TnlJUi8MFElSLwwUSUct33fSLwNFktQLA0WS1AsDRdJRxctcC8dAkST1wkCRJPXCQJEk9cJAkST1oq/vQ5GkJc3B+IVnD0WS1AsDRdKytX7HzfZMFpGBIknqhYEiSeqFgSJJ6oWBIknqhYEiaVlyMH7xGSiSpF4YKJKkXgwdKEkeS3JPkruS7G21E5LsSfJQ+7my1ZPkiiQTSe5O8prOdra15R9Ksm3YdkmSFldfPZR/XlWbqmpze74DuLWqNgK3tucA5wAb22M7cCUMAgi4FHgtcDpw6WQISZLGw0Jd8joPuKZNXwO8vVO/tgZuA45PsgY4G9hTVU9V1Y+BPcCWBWqbpGXKgfjR6iNQCvjLJHcm2d5qq6vqQJv+IbC6TZ8EPN5Zd1+rTVeXJI2JPj5t+Leran+SfwTsSfJAd2ZVVZLqYT+0wNoOsG7duj42KUnqydCBUlX728+DSb7KYAzkiSRrqupAu6R1sC2+Hzi5s/raVtsPnHlY/X9Msa+rgasBNm/e3EtISRp/XupaGoa65JXkhUleNDkNnAXcC+wCJu/U2gbc1KZ3ARe0u73OAH7SLo3dApyVZGUbjD+r1SRJY2LYHspq4KtJJrf151X135PcAdyQ5ELgB8A72/K7gXOBCeCnwHsAquqpJH8M3NGW+2hVPTVk2yRJi2ioQKmqR4B/NkX9SeBNU9QLuGiabe0Edg7THknS6PhOeUljyS/PWnoMFElSLwwUSVIvDBRJUi/6eGOjJC0ax02WLnsokqReGCiSpF4YKJKkXhgokpY833MyHgwUSVIvDBRJS5o9k/FhoEiSeuH7UCQtGd3eyGOXv2WELdF82EORJPXCQJE0Ut7BtXwYKJJGxiBZXhxDkbQoHB9Z/uyhSHrWVJef5tKLsMdxdLOHIh3lJkOg22tYv+Pm5z2fapnJ2lRBYo/k6GMPRTrKdF/o59ujsCeiqcw7UJKcnORbSe5Pcl+SP2j1jyTZn+Su9ji3s84lSSaSPJjk7E59S6tNJNkx3CFJmop3U2mhDXPJ6xngD6vqr5O8CLgzyZ4279NV9cnuwklOAbYCpwIvBb6Z5BVt9meANwP7gDuS7Kqq+4domySmvlQlLZR5B0pVHQAOtOn/l+R7wElHWOU84Pqq+hnwaJIJ4PQ2b6KqHgFIcn1b1kCRhmBvRIutlzGUJOuBVwO3t9LFSe5OsjPJylY7CXi8s9q+VpuuPtV+tifZm2TvoUOH+mi6JKknQwdKkt8AbgQ+VFVPA1cCLwc2MejBfGrYfUyqqquranNVbV61alVfm5Uk9WCoQEnyawzC5ItV9RWAqnqiqn5RVb8EPsevLmvtB07urL621aarS5ojL3NplIa5yyvA54HvVdV/6tTXdBZ7B3Bvm94FbE1yXJINwEbgO8AdwMYkG5Icy2Dgftd82yVJGo1h7vL6LeBdwD1J7mq1fw+cn2QTUMBjwPsBquq+JDcwGGx/Brioqn4BkORi4BZgBbCzqu4bol3SUceeiZaCYe7y+l9Appi1+wjrXAZcNkV995HWk/R83hKspcZ3yktjyB6JliIDRRoThoiWOgNFktQLA0WS1AsDRVqCDv9EYC93aRz4fSjSiM30vSLSuLCHIo2QAaLlxB6KtICm63343hEtR/ZQpHmY7lsP+/g2RGlc2UPRUWOm3sLh36PeXWeq9aZbTjpaHZWBMtUfv5cglq7ZBEG3NtMys92XNCpT/edmHByVgTIdX0xGa5ggkDR6jqFMYar3APjCNnez+Xf031VaPuyhzNFyuVPn8Bfy6Y5lrpeTZvteCoNEWn4MlB6MOmQOv946mxf96bYz0zKSNB0DZQHM9gV92P/9TxUiBoGkUXEMRZLUCwNFktQLA0WS1IslEyhJtiR5MMlEkh2jbo8kaW6WRKAkWQF8BjgHOAU4P8kpo22VJGkulkSgAKcDE1X1SFX9PXA9cN6I2yRJmoOlEignAY93nu9rNUnSmEhVjboNJPldYEtV/X57/i7gtVV18WHLbQe2t6evBB6cZpMnAj9aoOaOksc1XpbrccHyPbaj4bj+cVWtWoidLJU3Nu4HTu48X9tqz1FVVwNXz7SxJHuranN/zVsaPK7xslyPC5bvsXlcw1kql7zuADYm2ZDkWGArsGvEbZIkzcGS6KFU1TNJLgZuAVYAO6vqvhE3S5I0B0siUACqajewu6fNzXhZbEx5XONluR4XLN9j87iGsCQG5SVJ42+pjKFIksbcsgqUcfj4liQnJ/lWkvuT3JfkD1r9hCR7kjzUfq5s9SS5oh3T3Ule09nWtrb8Q0m2deqnJbmnrXNFkizi8a1I8t0kX2/PNyS5vbXlS+2mC5Ic155PtPnrO9u4pNUfTHJ2pz6y85vk+CRfTvJAku8led1yOGdJ/k37Pbw3yXVJfn0cz1mSnUkOJrm3U1vw8zPdPhb4uD7Rfg/vTvLVJMd35s3pPMznXB9RVS2LB4PB/IeBlwHHAv8HOGXU7ZqinWuA17TpFwHfZ/BxMx8HdrT6DuBjbfpc4BtAgDOA21v9BOCR9nNlm17Z5n2nLZu27jmLeHz/Fvhz4Ovt+Q3A1jZ9FfCv2vQHgava9FbgS236lHbujgM2tHO6YtTnF7gG+P02fSxw/LifMwZvHn4UeEHnXL17HM8Z8AbgNcC9ndqCn5/p9rHAx3UWcEyb/ljnuOZ8HuZ6rmds72L9QS7CH8frgFs6zy8BLhl1u2bR7puANzN4k+aaVlsDPNimPwuc31n+wTb/fOCznfpnW20N8ECn/pzlFvhY1gK3Am8Evt7++H7U+eV/9hwxuKPvdW36mLZcDj9vk8uN8vwCv8nghTeH1cf6nPGrT6g4oZ2DrwNnj+s5A9bz3BfeBT8/0+1jIY/rsHnvAL441b/vTOdhPn+fM7V1OV3yGruPb2ndyFcDtwOrq+pAm/VDYHWbnu64jlTfN0V9MfwJ8EfAL9vzlwB/U1XPTNGWZ9vf5v+kLT/X410MG4BDwH/N4HLenyV5IWN+zqpqP/BJ4P8CBxicgztZHucMFuf8TLePxfJeBj0mmPtxzefv84iWU6CMlSS/AdwIfKiqnu7Oq8F/C8bq9rskbwUOVtWdo27LAjiGwWWHK6vq1cDfMri88awxPWcrGXwI6wbgpcALgS0jbdQCWYzzs9i/A0k+DDwDfHGx9jmT5RQos/r4lqUgya8xCJMvVtVXWvmJJGva/DXAwVaf7riOVF87RX2h/RbwtiSPMfi06DcC/xk4Psnk+526bXm2/W3+bwJPMvfjXQz7gH1VdXt7/mUGATPu5+x3gEer6lBV/Rz4CoPzuBzOGSzO+ZluHwsqybuBtwK/14IM5n5cTzL3c31kC3U9c7EfDP4X+QiD/21NDjydOup2TdHOANcCf3JY/RM8d3Dv4236LTx3APE7rX4Cg+v6K9vjUeCENu/wAcRzF/kYz+RXg/J/wXMH/T7Ypi/iuYN+N7TpU3nuwOIjDAYVR3p+gf8JvLJNf6Sdr7E+Z8BrgfuAf9j2ew3wr8f1nPH8MZQFPz/T7WOBj2sLcD+w6rDl5nwe5nquZ2zrYv1BLsaDwd0b32dwR8OHR92eadr42wy6xXcDd7XHuQyuT94KPAR8s/OLHAZfPvYwcA+wubOt9wIT7fGeTn0zcG9b50+ZxWBaz8d4Jr8KlJe1P8aJ9st7XKv/ens+0ea/rLP+h1vbH6Rzt9Mozy+wCdjbztvX2gvO2J8z4D8AD7R9/7f2YjR25wy4jsE40M8Z9CgvXIzzM90+Fvi4JhiMb0y+flw13/Mwn3N9pIfvlJck9WI5jaFIkkbIQJEk9cJAkST1wkCRJPXCQJEk9cJAkST1wkCRJPXCQJEk9eL/A6lIw/n6E74TAAAAAElFTkSuQmCC\n", | |
"text/plain": [ | |
"<Figure size 432x288 with 1 Axes>" | |
] | |
}, | |
"metadata": { | |
"needs_background": "light" | |
}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"ch = event_list[\"boardIndexAndChannel\"]\n", | |
"energy_filt = energy[condition]\n", | |
"_ = plt.hist(energy_filt, bins=200)" | |
] | |
} | |
], | |
"metadata": { | |
"kernelspec": { | |
"display_name": "Python 3", | |
"language": "python", | |
"name": "python3" | |
}, | |
"language_info": { | |
"codemirror_mode": { | |
"name": "ipython", | |
"version": 2 | |
}, | |
"file_extension": ".py", | |
"mimetype": "text/x-python", | |
"name": "python", | |
"nbconvert_exporter": "python", | |
"pygments_lexer": "ipython2", | |
"version": "2.7.10" | |
} | |
}, | |
"nbformat": 4, | |
"nbformat_minor": 2 | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment