{ "cells": [ { "cell_type": "markdown", "id": "f28a6e79", "metadata": {}, "source": [ "# Example with dummy data\n", "\n", "---\n", "\n", "**Author:** Adrien CR Thob\n", "\n", "*© 2022 Adrien CR Thob*\n", "\n", "*This [notebook](https://github.com/athob/py-ananke/blob/main/jupyter/example_notebook.ipynb) is part of the [py-ananke](https://github.com/athob/py-ananke) project, which is licensed under the GNU Affero General Public License v3.0 (AGPL-3.0). The full copyright notice, including terms governing use, modification, and redistribution, is contained in the files [LICENSE](https://github.com/athob/py-ananke/blob/main/LICENSE) and [COPYRIGHT](https://github.com/athob/py-ananke/blob/main/COPYRIGHT), which can be found at the root of the source code distribution tree*\n", "\n", "**Updated On:** 2026-04-11\n", "\n", "---\n", "\n", "This notebook shows a simplified use case utilizing a dummy set of randomly generated particle data. It is by default set to simulate a catalog of stars with Gaia DR2 astrometry & photometry, but you may change the parameters defined below to your liking. In short, the pipeline is accessible via the class `Ananke`, which method `run` can eventually be called after creation of an `Ananke` object to run the full ananke pipeline. The method `run` is designed to returns the catalog in the form of an `Output` object of the submodule `galaxia_ananke` which can be used as a `vaex` dataframe.\n", "\n", "Below are the packages that this notebook must import:" ] }, { "cell_type": "code", "execution_count": 1, "id": "f130bf6c", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "'0.4.0+14.ge1134b7.dirty'" ] }, "execution_count": 1, "metadata": {}, "output_type": "execute_result" } ], "source": [ "import numpy as np\n", "import ananke as an\n", "an.__version__" ] }, { "attachments": {}, "cell_type": "markdown", "id": "d0809afe", "metadata": {}, "source": [ "We define here some dummy input data. `Ananke` has a method `make_dummy_particles_input` to produce such data in a dictionary formatted output. Below we show how to call it and the keys that compose that dictionary." ] }, { "cell_type": "code", "execution_count": 2, "id": "f4dafbb1", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "dict_keys(['pos3', 'vel3', 'mass', 'age', 'feh', 'helium', 'carbon', 'nitrogen', 'oxygen', 'neon', 'magnesium', 'silicon', 'sulphur', 'calcium', 'alpha', 'parentid', 'partitionid', 'dform', 'id', 'log10_NH'])" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "np.random.seed(0)\n", "p = an.Ananke.make_dummy_particles_input()\n", "p.keys()" ] }, { "cell_type": "markdown", "id": "51f0f1ae", "metadata": {}, "source": [ "The input data must be formatted as a dictionary of equal-length arrays. The dictionary must have the following entries:\n", "- key `pos3`: particle position coordinates in $kpc$ (shape Nx3)\n", "- key `vel3`: particle velocity coordinates in $km.s^{-1}$ (shape Nx3)\n", "- key `mass`: particle stellar mass in solar masses\n", "- key `age`: particle log10 stellar age in years\n", "- key `feh`: particle stellar metallicity \\[Fe/H\\] in dex relative to solar\n", "\n", "Additionally, the following entries can optionally be added:\n", "- key `parentid`: index to give to the parent particle\n", "- key `id`: additional index to classify the parent particle\n", "- key `log10_NH`: log10 hydrogen column densities between Observer position and particle in $cm^{-2}$ - must be provided to estimate extinctions\n", "- key `dform`: particle formation distance\n", "- keys `helium`, `carbon`, `nitrogen`, `oxygen`, `neon`, `magnesium`, `silicon`, `sulphur`, `calcium`: particle various chemical abundances \\[X/H\\]\n", "- key `alpha`: particle alpha chemical abundances \\[Mg/Fe\\]\n", "\n", "Ananke will compute the phase space densities that are used to determine particle smoothing lengths, but one can include pre-computed densities with the following entries:\n", "- key `rho_pos`: particle density in position space in $kpc^{-3}$\n", "- key `rho_vel`: particle density in velocity space in $km^{-3}.s^{3}$\n", "\n", "At any time, you may access this format via the helper of the `Ananke.make_dummy_particles_input` method using\n", "\n", "`help(an.Ananke.make_dummy_particles_input)`" ] }, { "attachments": {}, "cell_type": "markdown", "id": "55c50b11", "metadata": {}, "source": [ "We define below some parameters for Ananke such as\n", " - the observer position `observer`\n", " - the shell of particles to mask `rshell`\n", " - the sampling factor `fsample`\n", " - the photometric system of choice `photo_sys` (in our case Gaia DR2)\n", " - the CMD `cmd_magnames` and its box limits `cmd_box`" ] }, { "cell_type": "code", "execution_count": 3, "id": "1ecc74b8", "metadata": {}, "outputs": [], "source": [ "D = 200 # *units.kpc\n", "\n", "observer = np.nan*np.ones(3)\n", "while not np.linalg.norm(observer)<1:\n", " observer = 2*np.random.rand(3)-1\n", "\n", "observer *= D/np.linalg.norm(observer)\n", "\n", "rshell = [0, 2*D]\n", "\n", "fsample = 0.01\n", "\n", "photo_sys = 'padova/GAIA__DR2'\n", "\n", "cmd_magnames = {'magnitude': 'G',\n", " 'color_minuend': 'Gbp',\n", " 'color_subtrahend': 'Grp'}\n", "\n", "cmd_box = {\n", " 'abs_mag_lim_lo': -1000,\n", " 'abs_mag_lim_hi': 1000,\n", " # 'app_mag_lim_lo' : -1000,\n", " 'app_mag_lim_hi': 30,\n", " # 'color_lim_lo' : -1000,\n", " # 'color_lim_hi' : 1000\n", " }" ] }, { "cell_type": "markdown", "id": "8e3b700c", "metadata": {}, "source": [ "For more details regarding these parameters and more, you may consult the docstring associated to the class `Ananke` via the lines:\n", "\n", "`help(an.Ananke.__init__)`\n", "\n", "or\n", "\n", "`help(an.Ananke)`\n", "\n", "Note that among those are 3 special parameters, `d_params`, `e_params` and `err_params`, which receive dictionaries of subparameters to configure the sub-steps of the pipeline corresponding respectively to the phase space density estimation, the extinction mapping and the instrument error modeling. Each have further documentation that can be accessed calling the following methods:\n", "\n", "`an.display_density_docs()` for `d_params`\n", "\n", "`an.display_extinction_docs()` for `e_params`\n", "\n", "`an.display_errormodel_docs()` for `err_params`" ] }, { "attachments": {}, "cell_type": "markdown", "id": "1f132dcb", "metadata": {}, "source": [ "We show below how to create the `Ananke` pipeline object before calling `run`, using the parameters we defined above." ] }, { "cell_type": "code", "execution_count": 4, "id": "f56bb9c8", "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/athob/miniconda3/envs/galaxia_kernels/lib/python3.12/site-packages/ananke/utils.py:45: UserWarning: Pandas doesn't allow columns to be created via a new attribute name - see https://pandas.pydata.org/pandas-docs/stable/indexing.html#attribute-access\n", " self._record_of_all_used_keys = set()\n", "/home/athob/miniconda3/envs/galaxia_kernels/lib/python3.12/site-packages/ananke/utils.py:45: UserWarning: Pandas doesn't allow columns to be created via a new attribute name - see https://pandas.pydata.org/pandas-docs/stable/indexing.html#attribute-access\n", " self._record_of_all_used_keys = set()\n" ] } ], "source": [ "name = 'sim'\n", "ananke = an.Ananke(p, name, fsample=fsample,\n", " observer=observer, rshell=rshell,\n", " photo_sys=photo_sys, cmd_magnames=cmd_magnames,\n", " **cmd_box)" ] }, { "attachments": {}, "cell_type": "markdown", "id": "fb2ba696", "metadata": {}, "source": [ "As mentioned above, the `Ananke` method `run` runs the pipeline once it has been created. Some additional parameters can be provided to this method, which are listed in its docstring (accessible by calling `help(ananke.run)`). Those include:\n", "- `input_dir`, `output_dir`, `i_o_dir` - Optional arguments to specify paths for the directories where `ananke` should generate input and output data. If the `i_o_dir` keyword argument is provided, it overrides any path given to the `input_dir` and `output_dir` keyword arguments.\n", "- `k_factor` – Scaling factor applied to the kernels lengths to adjust all the kernels sizes uniformly. Lower values reduces the kernels extents, while higher values increases them. Default to 1 (no adjustment).\n", "- `surveyname` – Optional name `Galaxia` should use for the output files. Default to `'survey'`.\n", "- `n_jobs` – Number of independent catalog generations ran in parallel. Default to 1.\n", "- `verbose` – Verbose boolean flag to allow pipeline to print what it's doing to stdout. Default to `True`.\n", "- `rand_seed` – Seed to be used by `Galaxia`'s pseudorandom number generator. Default to 17052\n", "- `nstart` – Index at which to start indexing synthetic stars. Default to 0\n", "For our example however, we will use the method as it is.\n", "\n", "\n", "\n", "
\n", "
Warning
\n", "

Some of the `run` method parameters are currently not properly implemented, using them may lead to unexpected results

\n", "
" ] }, { "cell_type": "code", "execution_count": 5, "id": "f04f19ca", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Dimensions = 3\n", "Normalization constant of Kernel type 3: 0.596831\n", "Reading ASCII format file: to_enbid \n", "Read 300000 records\n", "Allocated 3.8147 MByte for particle storage.\n", "Type = 1 Particles = 100000\n", "Total particles = 100000\n", "Allocated 22.1252 MByte for Binary-Tree and Entropy bins.\n", "\n", "Scaling Co-ordinates as x[i]=x[i]/h[i] with h[i]->\n", "1 1 1 \n", "\n", "Starting to Build Tree .......\n", "Particle Type = 1 First node number = 1\n", " Last node number = 199999\n", "Total number of nodes = 199999\n", "Treebuild time = 0.176679 s \n", "\n", "Density Calculation. Smoothing .....\n", "Evaluated = 1 % Time Left = 1.146210 s of 1.157800 s Par no = 96805 Density = 4.540411e-02 \n", "Evaluated = 2 % Time Left = 1.141394 s of 1.164700 s Par no = 55454 Density = 4.964532e-03 \n", "Evaluated = 3 % Time Left = 1.201074 s of 1.238233 s Par no = 27717 Density = 1.827217e-01 \n", "Evaluated = 4 % Time Left = 1.267739 s of 1.320575 s Par no = 42234 Density = 2.712376e-02 \n", "Evaluated = 5 % Time Left = 1.285868 s of 1.353560 s Par no = 99785 Density = 2.988043e-01 \n", "Evaluated = 6 % Time Left = 1.299552 s of 1.382517 s Par no = 88021 Density = 1.791133e-01 \n", "Evaluated = 7 % Time Left = 1.327773 s of 1.427729 s Par no = 75708 Density = 1.165058e-01 \n", "Evaluated = 8 % Time Left = 1.306328 s of 1.419938 s Par no = 45579 Density = 7.255407e-02 \n", "Evaluated = 9 % Time Left = 1.288758 s of 1.416233 s Par no = 30351 Density = 8.828938e-02 \n", "Evaluated = 10 % Time Left = 1.275331 s of 1.417050 s Par no = 89176 Density = 1.775999e-01 \n", "Evaluated = 11 % Time Left = 1.253535 s of 1.408482 s Par no = 60361 Density = 8.299880e-02 \n", "Evaluated = 12 % Time Left = 1.248420 s of 1.418675 s Par no = 90082 Density = 2.203585e-01 \n", "Evaluated = 13 % Time Left = 1.230641 s of 1.414546 s Par no = 26814 Density = 1.910994e-02 \n", "Evaluated = 14 % Time Left = 1.209232 s of 1.406100 s Par no = 74297 Density = 2.157650e-02 \n", "Evaluated = 15 % Time Left = 1.199432 s of 1.411113 s Par no = 94905 Density = 9.119637e-02 \n", "Evaluated = 16 % Time Left = 1.183997 s of 1.409538 s Par no = 51653 Density = 9.666827e-02 \n", "Evaluated = 17 % Time Left = 1.161888 s of 1.399882 s Par no = 39308 Density = 1.287774e-01 \n", "Evaluated = 18 % Time Left = 1.142816 s of 1.393694 s Par no = 7800 Density = 4.130107e-02 \n", "Evaluated = 19 % Time Left = 1.120583 s of 1.383453 s Par no = 42838 Density = 5.914271e-02 \n", "Evaluated = 20 % Time Left = 1.101118 s of 1.376415 s Par no = 49016 Density = 1.554842e-01 \n", "Evaluated = 21 % Time Left = 1.084502 s of 1.372805 s Par no = 29173 Density = 5.981785e-02 \n", "Evaluated = 22 % Time Left = 1.063286 s of 1.363205 s Par no = 14018 Density = 4.713366e-02 \n", "Evaluated = 23 % Time Left = 1.046336 s of 1.358896 s Par no = 90617 Density = 4.088210e-02 \n", "Evaluated = 24 % Time Left = 1.033124 s of 1.359392 s Par no = 20588 Density = 1.273882e-01 \n", "Evaluated = 25 % Time Left = 1.017022 s of 1.356048 s Par no = 15465 Density = 1.274245e-01 \n", "Evaluated = 26 % Time Left = 1.005453 s of 1.358738 s Par no = 99360 Density = 2.044210e-01 \n", "Evaluated = 27 % Time Left = 0.989583 s of 1.355611 s Par no = 71132 Density = 2.079416e-02 \n", "Evaluated = 28 % Time Left = 0.971526 s of 1.349361 s Par no = 6665 Density = 4.366049e-03 \n", "Evaluated = 29 % Time Left = 0.959243 s of 1.351066 s Par no = 44439 Density = 1.463732e-01 \n", "Evaluated = 30 % Time Left = 0.942306 s of 1.346170 s Par no = 6386 Density = 1.825200e-01 \n", "Evaluated = 31 % Time Left = 0.928796 s of 1.346100 s Par no = 2716 Density = 1.495580e-01 \n", "Evaluated = 32 % Time Left = 0.913660 s of 1.343638 s Par no = 69470 Density = 6.785184e-02 \n", "Evaluated = 33 % Time Left = 0.901914 s of 1.346161 s Par no = 95884 Density = 1.793627e-01 \n", "Evaluated = 34 % Time Left = 0.888447 s of 1.346153 s Par no = 73768 Density = 1.295312e-02 \n", "Evaluated = 35 % Time Left = 0.875124 s of 1.346366 s Par no = 47898 Density = 1.088661e-01 \n", "Evaluated = 36 % Time Left = 0.857914 s of 1.340511 s Par no = 66396 Density = 3.959183e-02 \n", "Evaluated = 37 % Time Left = 0.842889 s of 1.337941 s Par no = 16463 Density = 1.572218e-01 \n", "Evaluated = 38 % Time Left = 0.831748 s of 1.341550 s Par no = 55669 Density = 1.712036e-01 \n", "Evaluated = 39 % Time Left = 0.820568 s of 1.345215 s Par no = 97705 Density = 1.234432e-01 \n", "Evaluated = 40 % Time Left = 0.806258 s of 1.343785 s Par no = 43442 Density = 6.719712e-02 \n", "Evaluated = 41 % Time Left = 0.792800 s of 1.343751 s Par no = 10242 Density = 2.102901e-01 \n", "Evaluated = 42 % Time Left = 0.781737 s of 1.347845 s Par no = 40929 Density = 1.391802e-01 \n", "Evaluated = 43 % Time Left = 0.768943 s of 1.349047 s Par no = 5720 Density = 6.551559e-02 \n", "Evaluated = 44 % Time Left = 0.754300 s of 1.346989 s Par no = 66730 Density = 9.093650e-02 \n", "Evaluated = 45 % Time Left = 0.740453 s of 1.346302 s Par no = 48819 Density = 2.248662e-01 \n", "Evaluated = 46 % Time Left = 0.728266 s of 1.348665 s Par no = 32014 Density = 9.211309e-02 \n", "Evaluated = 47 % Time Left = 0.714427 s of 1.348000 s Par no = 31533 Density = 2.932795e-02 \n", "Evaluated = 48 % Time Left = 0.699242 s of 1.344723 s Par no = 86677 Density = 1.263195e-01 \n", "Evaluated = 49 % Time Left = 0.687264 s of 1.347602 s Par no = 81951 Density = 8.802204e-02 \n", "Evaluated = 50 % Time Left = 0.675195 s of 1.350418 s Par no = 80149 Density = 1.559104e-01 \n", "Evaluated = 51 % Time Left = 0.660847 s of 1.348694 s Par no = 25229 Density = 5.093939e-02 \n", "Evaluated = 52 % Time Left = 0.647404 s of 1.348787 s Par no = 89731 Density = 4.987880e-02 \n", "Evaluated = 53 % Time Left = 0.632298 s of 1.345343 s Par no = 76366 Density = 1.511494e-01 \n", "Evaluated = 54 % Time Left = 0.618909 s of 1.345483 s Par no = 73159 Density = 1.548675e-02 \n", "Evaluated = 55 % Time Left = 0.605247 s of 1.345024 s Par no = 54007 Density = 8.763470e-02 \n", "Evaluated = 56 % Time Left = 0.591500 s of 1.344348 s Par no = 36235 Density = 5.465964e-02 \n", "Evaluated = 57 % Time Left = 0.577849 s of 1.343867 s Par no = 62068 Density = 4.050222e-02 \n", "Evaluated = 58 % Time Left = 0.562545 s of 1.339424 s Par no = 37661 Density = 2.656841e-02 \n", "Evaluated = 59 % Time Left = 0.549200 s of 1.339544 s Par no = 67914 Density = 1.940649e-01 \n", "Evaluated = 60 % Time Left = 0.536683 s of 1.341742 s Par no = 64732 Density = 5.859731e-02 \n", "Evaluated = 61 % Time Left = 0.521806 s of 1.337998 s Par no = 88788 Density = 7.161247e-03 \n", "Evaluated = 62 % Time Left = 0.507941 s of 1.336721 s Par no = 44083 Density = 2.877787e-02 \n", "Evaluated = 63 % Time Left = 0.493786 s of 1.334592 s Par no = 14045 Density = 3.265437e-02 \n", "Evaluated = 64 % Time Left = 0.479127 s of 1.330945 s Par no = 93621 Density = 5.216302e-02 \n", "Evaluated = 65 % Time Left = 0.465022 s of 1.328672 s Par no = 9132 Density = 8.411557e-02 \n", "Evaluated = 66 % Time Left = 0.452096 s of 1.329732 s Par no = 702 Density = 1.110411e-02 \n", "Evaluated = 67 % Time Left = 0.439246 s of 1.331090 s Par no = 45670 Density = 2.465222e-02 \n", "Evaluated = 68 % Time Left = 0.425309 s of 1.329131 s Par no = 53860 Density = 1.580593e-02 \n", "Evaluated = 69 % Time Left = 0.411221 s of 1.326562 s Par no = 79719 Density = 4.454536e-02 \n", "Evaluated = 70 % Time Left = 0.397054 s of 1.323559 s Par no = 96401 Density = 1.035200e-02 \n", "Evaluated = 71 % Time Left = 0.382895 s of 1.320372 s Par no = 48105 Density = 3.704987e-02 \n", "Evaluated = 72 % Time Left = 0.368985 s of 1.317851 s Par no = 75056 Density = 2.642116e-02 \n", "Evaluated = 73 % Time Left = 0.355416 s of 1.316405 s Par no = 97912 Density = 1.775965e-01 \n", "Evaluated = 74 % Time Left = 0.341894 s of 1.315027 s Par no = 87879 Density = 7.451773e-02 \n", "Evaluated = 75 % Time Left = 0.328814 s of 1.315307 s Par no = 85159 Density = 1.462781e-02 \n", "Evaluated = 76 % Time Left = 0.315834 s of 1.316030 s Par no = 78132 Density = 5.302110e-02 \n", "Evaluated = 77 % Time Left = 0.302310 s of 1.314448 s Par no = 21990 Density = 3.497640e-03 \n", "Evaluated = 78 % Time Left = 0.288448 s of 1.311186 s Par no = 16462 Density = 2.482335e-02 \n", "Evaluated = 79 % Time Left = 0.275033 s of 1.309743 s Par no = 22311 Density = 3.571432e-02 \n", "Evaluated = 80 % Time Left = 0.261456 s of 1.307345 s Par no = 43269 Density = 1.776266e-02 \n", "Evaluated = 81 % Time Left = 0.248436 s of 1.307628 s Par no = 43142 Density = 1.532303e-02 \n", "Evaluated = 82 % Time Left = 0.235222 s of 1.306861 s Par no = 52197 Density = 3.049044e-02 \n", "Evaluated = 83 % Time Left = 0.221922 s of 1.305500 s Par no = 74476 Density = 1.449572e-01 \n", "Evaluated = 84 % Time Left = 0.208860 s of 1.305456 s Par no = 8891 Density = 1.050831e-02 \n", "Evaluated = 85 % Time Left = 0.195778 s of 1.305273 s Par no = 1559 Density = 1.452041e-03 \n", "Evaluated = 86 % Time Left = 0.182713 s of 1.305185 s Par no = 20602 Density = 2.028735e-01 \n", "Evaluated = 87 % Time Left = 0.169825 s of 1.306448 s Par no = 54383 Density = 2.137116e-01 \n", "Evaluated = 88 % Time Left = 0.157068 s of 1.309008 s Par no = 70246 Density = 2.757094e-03 \n", "Evaluated = 89 % Time Left = 0.144046 s of 1.309625 s Par no = 51472 Density = 5.058417e-04 \n", "Evaluated = 90 % Time Left = 0.130866 s of 1.308787 s Par no = 66269 Density = 3.630537e-02 \n", "Evaluated = 91 % Time Left = 0.117969 s of 1.310915 s Par no = 62858 Density = 1.776033e-01 \n", "Evaluated = 92 % Time Left = 0.105036 s of 1.313115 s Par no = 84061 Density = 1.439860e-01 \n", "Evaluated = 93 % Time Left = 0.092021 s of 1.314775 s Par no = 78610 Density = 4.029672e-03 \n", "Evaluated = 94 % Time Left = 0.078931 s of 1.315732 s Par no = 3871 Density = 1.597447e-01 \n", "Evaluated = 95 % Time Left = 0.065875 s of 1.317764 s Par no = 68736 Density = 1.582606e-01 \n", "Evaluated = 96 % Time Left = 0.052799 s of 1.320299 s Par no = 30472 Density = 6.743145e-02 \n", "Evaluated = 97 % Time Left = 0.039605 s of 1.320600 s Par no = 86838 Density = 1.353215e-03 \n", "Evaluated = 98 % Time Left = 0.026424 s of 1.321883 s Par no = 92913 Density = 5.211162e-02 \n", "Evaluated = 99 % Time Left = 0.013232 s of 1.324565 s Par no = 30015 Density = 9.179291e-04 \n", "Evaluated = 100 % Time Left = -0.000013 s of 1.327918 s Par no = 86525 Density = 1.579017e-04 \n", "\n", "Total Smoothing Time = 1.32805 s\n", "1 3\n", "\n", "writing snapshot file ....\n", "to_enbid_d3n64.est\n", "done with snapshot.\n", "\n", "Total Time = 1.8399 s \n", "\n", "Dimensions = 3\n", "Normalization constant of Kernel type 3: 0.596831\n", "Reading ASCII format file: to_enbid \n", "Read 300000 records\n", "Allocated 3.8147 MByte for particle storage.\n", "Type = 1 Particles = 100000\n", "Total particles = 100000\n", "Allocated 22.1252 MByte for Binary-Tree and Entropy bins.\n", "\n", "Scaling Co-ordinates as x[i]=x[i]/h[i] with h[i]->\n", "1 1 1 \n", "\n", "Starting to Build Tree .......\n", "Particle Type = 1 First node number = 1\n", " Last node number = 199999\n", "Total number of nodes = 199999\n", "Treebuild time = 0.186382 s \n", "\n", "Density Calculation. Smoothing .....\n", "Evaluated = 1 % Time Left = 1.292036 s of 1.305100 s Par no = 3379 Density = 1.586355e-02 \n", "Evaluated = 2 % Time Left = 1.232680 s of 1.257850 s Par no = 75725 Density = 4.756319e-02 \n", "Evaluated = 3 % Time Left = 1.252451 s of 1.291200 s Par no = 95160 Density = 6.928090e-03 \n", "Evaluated = 4 % Time Left = 1.232075 s of 1.283425 s Par no = 40983 Density = 2.235913e-02 \n", "Evaluated = 5 % Time Left = 1.232232 s of 1.297100 s Par no = 79587 Density = 2.321368e-02 \n", "Evaluated = 6 % Time Left = 1.252834 s of 1.332817 s Par no = 62143 Density = 2.106821e-02 \n", "Evaluated = 7 % Time Left = 1.230669 s of 1.323314 s Par no = 48098 Density = 2.322558e-02 \n", "Evaluated = 8 % Time Left = 1.221666 s of 1.327913 s Par no = 73424 Density = 1.566558e-02 \n", "Evaluated = 9 % Time Left = 1.210135 s of 1.329833 s Par no = 41738 Density = 1.527212e-02 \n", "Evaluated = 10 % Time Left = 1.199624 s of 1.332930 s Par no = 34821 Density = 3.462067e-02 \n", "Evaluated = 11 % Time Left = 1.200071 s of 1.348409 s Par no = 99566 Density = 5.806165e-02 \n", "Evaluated = 12 % Time Left = 1.184511 s of 1.346050 s Par no = 59443 Density = 1.071449e-02 \n", "Evaluated = 13 % Time Left = 1.159094 s of 1.332308 s Par no = 7663 Density = 2.551240e-03 \n", "Evaluated = 14 % Time Left = 1.159131 s of 1.347843 s Par no = 69540 Density = 2.717029e-02 \n", "Evaluated = 15 % Time Left = 1.140262 s of 1.341500 s Par no = 28859 Density = 1.270314e-02 \n", "Evaluated = 16 % Time Left = 1.115339 s of 1.327800 s Par no = 77774 Density = 2.211896e-02 \n", "Evaluated = 17 % Time Left = 1.109648 s of 1.336941 s Par no = 65109 Density = 3.644556e-03 \n", "Evaluated = 18 % Time Left = 1.090527 s of 1.329928 s Par no = 17146 Density = 6.157028e-03 \n", "Evaluated = 19 % Time Left = 1.084244 s of 1.338589 s Par no = 79145 Density = 2.474514e-02 \n", "Evaluated = 20 % Time Left = 1.069751 s of 1.337205 s Par no = 84563 Density = 3.054690e-02 \n", "Evaluated = 21 % Time Left = 1.055863 s of 1.336552 s Par no = 35101 Density = 3.370794e-02 \n", "Evaluated = 22 % Time Left = 1.041737 s of 1.335577 s Par no = 31911 Density = 1.401559e-03 \n", "Evaluated = 23 % Time Left = 1.030752 s of 1.338657 s Par no = 96514 Density = 1.869492e-02 \n", "Evaluated = 24 % Time Left = 1.015157 s of 1.335750 s Par no = 41208 Density = 4.063115e-02 \n", "Evaluated = 25 % Time Left = 0.999020 s of 1.332044 s Par no = 50899 Density = 1.130457e-03 \n", "Evaluated = 26 % Time Left = 0.982192 s of 1.327304 s Par no = 25354 Density = 2.573581e-02 \n", "Evaluated = 27 % Time Left = 0.968507 s of 1.326741 s Par no = 38397 Density = 1.824382e-02 \n", "Evaluated = 28 % Time Left = 0.960030 s of 1.333393 s Par no = 99925 Density = 5.098955e-02 \n", "Evaluated = 29 % Time Left = 0.949878 s of 1.337876 s Par no = 67233 Density = 3.835664e-02 \n", "Evaluated = 30 % Time Left = 0.937886 s of 1.339857 s Par no = 21415 Density = 1.821073e-02 \n", "Evaluated = 31 % Time Left = 0.922078 s of 1.336365 s Par no = 39830 Density = 5.986182e-03 \n", "Evaluated = 32 % Time Left = 0.904372 s of 1.329978 s Par no = 62656 Density = 1.851001e-03 \n", "Evaluated = 33 % Time Left = 0.893253 s of 1.333233 s Par no = 48049 Density = 1.832367e-02 \n", "Evaluated = 34 % Time Left = 0.881638 s of 1.335835 s Par no = 71467 Density = 3.566542e-02 \n", "Evaluated = 35 % Time Left = 0.869493 s of 1.337703 s Par no = 94467 Density = 8.080338e-03 \n", "Evaluated = 36 % Time Left = 0.852868 s of 1.332628 s Par no = 24663 Density = 3.556317e-02 \n", "Evaluated = 37 % Time Left = 0.837934 s of 1.330076 s Par no = 55355 Density = 2.436876e-02 \n", "Evaluated = 38 % Time Left = 0.826399 s of 1.332924 s Par no = 77628 Density = 3.192216e-02 \n", "Evaluated = 39 % Time Left = 0.813989 s of 1.334431 s Par no = 6428 Density = 3.955819e-02 \n", "Evaluated = 40 % Time Left = 0.800352 s of 1.333943 s Par no = 29219 Density = 3.686058e-04 \n", "Evaluated = 41 % Time Left = 0.786888 s of 1.333732 s Par no = 11475 Density = 1.075123e-02 \n", "Evaluated = 42 % Time Left = 0.776366 s of 1.338586 s Par no = 25064 Density = 3.927625e-02 \n", "Evaluated = 43 % Time Left = 0.762786 s of 1.338244 s Par no = 14672 Density = 1.066699e-02 \n", "Evaluated = 44 % Time Left = 0.747779 s of 1.335343 s Par no = 46917 Density = 2.114631e-02 \n", "Evaluated = 45 % Time Left = 0.733245 s of 1.333198 s Par no = 35522 Density = 9.714150e-03 \n", "Evaluated = 46 % Time Left = 0.720657 s of 1.334574 s Par no = 10564 Density = 2.478037e-03 \n", "Evaluated = 47 % Time Left = 0.706405 s of 1.332864 s Par no = 90836 Density = 3.490644e-02 \n", "Evaluated = 48 % Time Left = 0.692669 s of 1.332081 s Par no = 91231 Density = 8.647473e-03 \n", "Evaluated = 49 % Time Left = 0.679982 s of 1.333324 s Par no = 16819 Density = 9.508127e-03 \n", "Evaluated = 50 % Time Left = 0.664736 s of 1.329498 s Par no = 30766 Density = 4.927925e-03 \n", "Evaluated = 51 % Time Left = 0.651281 s of 1.329173 s Par no = 9459 Density = 7.940475e-03 \n", "Evaluated = 52 % Time Left = 0.637009 s of 1.327131 s Par no = 76505 Density = 8.072356e-03 \n", "Evaluated = 53 % Time Left = 0.622958 s of 1.325472 s Par no = 52019 Density = 3.462386e-03 \n", "Evaluated = 54 % Time Left = 0.609822 s of 1.325728 s Par no = 56938 Density = 1.337130e-02 \n", "Evaluated = 55 % Time Left = 0.595117 s of 1.322511 s Par no = 33321 Density = 2.006831e-02 \n", "Evaluated = 56 % Time Left = 0.581560 s of 1.321757 s Par no = 4046 Density = 2.271462e-02 \n", "Evaluated = 57 % Time Left = 0.569096 s of 1.323511 s Par no = 68529 Density = 3.147267e-02 \n", "Evaluated = 58 % Time Left = 0.557718 s of 1.327933 s Par no = 23845 Density = 3.357711e-02 \n", "Evaluated = 59 % Time Left = 0.545196 s of 1.329778 s Par no = 57420 Density = 3.986822e-02 \n", "Evaluated = 60 % Time Left = 0.531483 s of 1.328742 s Par no = 89896 Density = 3.369183e-03 \n", "Evaluated = 61 % Time Left = 0.518158 s of 1.328644 s Par no = 39610 Density = 1.535302e-03 \n", "Evaluated = 62 % Time Left = 0.503581 s of 1.325247 s Par no = 26107 Density = 4.015532e-03 \n", "Evaluated = 63 % Time Left = 0.490500 s of 1.325713 s Par no = 98943 Density = 1.799383e-03 \n", "Evaluated = 64 % Time Left = 0.476865 s of 1.324662 s Par no = 99085 Density = 1.412553e-02 \n", "Evaluated = 65 % Time Left = 0.463286 s of 1.323712 s Par no = 17 Density = 1.179375e-02 \n", "Evaluated = 66 % Time Left = 0.448900 s of 1.320332 s Par no = 3087 Density = 1.429490e-02 \n", "Evaluated = 67 % Time Left = 0.436383 s of 1.322412 s Par no = 86468 Density = 2.256122e-02 \n", "Evaluated = 68 % Time Left = 0.422731 s of 1.321076 s Par no = 52080 Density = 4.635496e-03 \n", "Evaluated = 69 % Time Left = 0.408457 s of 1.317646 s Par no = 26417 Density = 1.433997e-02 \n", "Evaluated = 70 % Time Left = 0.395156 s of 1.317230 s Par no = 79634 Density = 7.429452e-03 \n", "Evaluated = 71 % Time Left = 0.382125 s of 1.317717 s Par no = 97643 Density = 2.184668e-04 \n", "Evaluated = 72 % Time Left = 0.369273 s of 1.318878 s Par no = 78794 Density = 5.369595e-03 \n", "Evaluated = 73 % Time Left = 0.355545 s of 1.316884 s Par no = 45708 Density = 1.294336e-02 \n", "Evaluated = 74 % Time Left = 0.341985 s of 1.315378 s Par no = 66525 Density = 2.125199e-03 \n", "Evaluated = 75 % Time Left = 0.328374 s of 1.313548 s Par no = 21984 Density = 8.619656e-04 \n", "Evaluated = 76 % Time Left = 0.315002 s of 1.312563 s Par no = 43799 Density = 2.740289e-03 \n", "Evaluated = 77 % Time Left = 0.301307 s of 1.310087 s Par no = 63854 Density = 3.854665e-03 \n", "Evaluated = 78 % Time Left = 0.287887 s of 1.308636 s Par no = 91160 Density = 2.256671e-03 \n", "Evaluated = 79 % Time Left = 0.274637 s of 1.307859 s Par no = 1891 Density = 1.555006e-04 \n", "Evaluated = 80 % Time Left = 0.261657 s of 1.308353 s Par no = 17306 Density = 9.342227e-04 \n", "Evaluated = 81 % Time Left = 0.248832 s of 1.309714 s Par no = 83119 Density = 2.469888e-02 \n", "Evaluated = 82 % Time Left = 0.235883 s of 1.310537 s Par no = 73346 Density = 1.956386e-03 \n", "Evaluated = 83 % Time Left = 0.222614 s of 1.309570 s Par no = 98412 Density = 4.288590e-03 \n", "Evaluated = 84 % Time Left = 0.209807 s of 1.311374 s Par no = 52509 Density = 8.264857e-03 \n", "Evaluated = 85 % Time Left = 0.197188 s of 1.314676 s Par no = 76716 Density = 3.547589e-02 \n", "Evaluated = 86 % Time Left = 0.184462 s of 1.317680 s Par no = 95957 Density = 2.923764e-02 \n", "Evaluated = 87 % Time Left = 0.171235 s of 1.317291 s Par no = 3295 Density = 1.361911e-02 \n", "Evaluated = 88 % Time Left = 0.158085 s of 1.317487 s Par no = 4026 Density = 3.781490e-02 \n", "Evaluated = 89 % Time Left = 0.144992 s of 1.318226 s Par no = 78656 Density = 6.492402e-02 \n", "Evaluated = 90 % Time Left = 0.132167 s of 1.321801 s Par no = 9255 Density = 4.661644e-02 \n", "Evaluated = 91 % Time Left = 0.119129 s of 1.323800 s Par no = 89389 Density = 2.612382e-03 \n", "Evaluated = 92 % Time Left = 0.105705 s of 1.321476 s Par no = 91357 Density = 1.136847e-02 \n", "Evaluated = 93 % Time Left = 0.092644 s of 1.323670 s Par no = 75212 Density = 1.186542e-02 \n", "Evaluated = 94 % Time Left = 0.079625 s of 1.327313 s Par no = 79246 Density = 4.605766e-02 \n", "Evaluated = 95 % Time Left = 0.066454 s of 1.329338 s Par no = 93038 Density = 2.164083e-02 \n", "Evaluated = 96 % Time Left = 0.053176 s of 1.329741 s Par no = 65384 Density = 2.308048e-02 \n", "Evaluated = 97 % Time Left = 0.039831 s of 1.328131 s Par no = 32352 Density = 1.225704e-03 \n", "Evaluated = 98 % Time Left = 0.026585 s of 1.329904 s Par no = 16971 Density = 3.436726e-02 \n", "Evaluated = 99 % Time Left = 0.013303 s of 1.331632 s Par no = 22517 Density = 2.430708e-02 \n", "Evaluated = 100 % Time Left = -0.000013 s of 1.334523 s Par no = 42877 Density = 2.940734e-06 \n", "\n", "Total Smoothing Time = 1.33457 s\n", "1 3\n", "\n", "writing snapshot file ....\n", "to_enbid_d3n64.est\n", "done with snapshot.\n", "\n", "Total Time = 1.86428 s \n", "\n", "Executing JOB 1/1 = /home/athob/miniconda3/envs/galaxia_kernels/.cache/galaxia_ananke/bin/galaxia -r --hdim=6 --nfile=sim --ngen=0 /home/athob/Software/repos/py-ananke/jupyter/survey_params\n", "JOB 1/1 | CODEDATAPATH=/home/athob/miniconda3/envs/galaxia_kernels/.cache/galaxia_ananke/GalaxiaData/\n", "JOB 1/1 | Reading Parameter file- /home/athob/Software/repos/py-ananke/jupyter/survey_params\n", "JOB 1/1 | --------------------------------------------------------\n", "JOB 1/1 | outputFile survey \n", "JOB 1/1 | outputDir /home/athob/Software/repos/py-ananke/jupyter\n", "JOB 1/1 | photoCateg padova \n", "JOB 1/1 | photoSys GAIA__DR2 \n", "JOB 1/1 | magcolorNames G,Gbp-Grp \n", "JOB 1/1 | appMagLimits[0] -1000 \n", "JOB 1/1 | appMagLimits[1] 30 \n", "JOB 1/1 | absMagLimits[0] -1000 \n", "JOB 1/1 | absMagLimits[1] 1000 \n", "JOB 1/1 | colorLimits[0] -1000 \n", "JOB 1/1 | colorLimits[1] 1000 \n", "JOB 1/1 | geometryOption 0 \n", "JOB 1/1 | surveyArea 207.455 \n", "JOB 1/1 | fSample 0.01 \n", "JOB 1/1 | popID 10 \n", "JOB 1/1 | warpFlareOn 0 \n", "JOB 1/1 | longitude 76.273 \n", "JOB 1/1 | latitude 13.4725 \n", "JOB 1/1 | starType 0 \n", "JOB 1/1 | photoError 0 \n", "JOB 1/1 | seed 17052 \n", "JOB 1/1 | r_max 400 \n", "JOB 1/1 | r_min 0 \n", "JOB 1/1 | nres 64 \n", "JOB 1/1 | nstart 0 \n", "JOB 1/1 | rSun[0] 194.95682580382555 \n", "JOB 1/1 | rSun[1] -34.20370085093144 \n", "JOB 1/1 | rSun[2] 28.66954691997796 \n", "JOB 1/1 | vSun[0] 12.9 \n", "JOB 1/1 | vSun[1] 245.6 \n", "JOB 1/1 | vSun[2] 7.78 \n", "JOB 1/1 | --------------------------------------------------------\n", "JOB 1/1 | Reading Halo Sat File=/home/athob/miniconda3/envs/galaxia_kernels/.cache/galaxia_ananke/GalaxiaData/nbody1/filenames/sim.txt\n", "JOB 1/1 | nbody1/sim/\n", "JOB 1/1 | \n", "JOB 1/1 | \t1\t1\n", "JOB 1/1 | \n", "JOB 1/1 | path nbody1/sim/ sats 1 x[0] 1\n", "JOB 1/1 | Halo 1 Sat 0 fname nbody1/sim/sim.ebf 0\n", "JOB 1/1 | nbody1/sim/sim.ebf 0\n", "JOB 1/1 | No of Satellites =1\n", "JOB 1/1 | Generating catalog. Max allowed stars is 18446744073709551615\n", "JOB 1/1 | Source numbering will start at 0\n", "JOB 1/1 | Writing to /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.ebf\n", "JOB 1/1 | Using geometry: All Sky\n", "JOB 1/1 | setting center to 194.957 -34.2037 28.6695 12.9 245.6 7.78\n", "JOB 1/1 | Reading Isochrones from dir- /home/athob/miniconda3/envs/galaxia_kernels/.cache/galaxia_ananke/GalaxiaData/Isochrones/padova/GAIA__DR2\n", "JOB 1/1 | Isochrone Grid Size: (Age bins=71,Feh bins=34,Alpha bins=1)\n", "JOB 1/1 | Time Isochrone Reading 3.71064 \n", "JOB 1/1 | ------------------------------\n", "JOB 1/1 | nbody1/sim/sim.ebf Sat No=0\n", "JOB 1/1 | Particles=100000\n", "JOB 1/1 | Satellite Info\n", "JOB 1/1 | Particles=100000 Mass=5.49989e+08 0.379706\n", "JOB 1/1 | Total Stars=2213240 accepted=1611712 rejected=601528\n", "JOB 1/1 | -----------Done---------------\n", "JOB 1/1 | Total stars written 1611712 \n", "JOB 1/1 | 31\n", "JOB 1/1 | File written- /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.ebf\n", "JOB 1/1 | Calculating magnitudes in GAIA__DR2 system................\n", "JOB 1/1 | initializing isochrone data\n", "JOB 1/1 | interpolating on isochrone tables\n", "JOB 1/1 | Total Time= 9.21257 \n", "Exported the following quantities from /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.ebf to /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.h5 for partition 0\n", "['age', 'alpha', 'calcium', 'carbon', 'dform', 'dmod', 'feh', 'gaia__dr2_g', 'gaia__dr2_gbp', 'gaia__dr2_grp', 'grav', 'helium', 'lum', 'mact', 'magnesium', 'mtip', 'neon', 'nitrogen', 'oxygen', 'parentid', 'partid', 'partitionid', 'px', 'py', 'pz', 'rad', 'satid', 'silicon', 'smass', 'sulphur', 'teff', 'vx', 'vy', 'vz']\n", "Running convert_cartesian_to_galactic post-processing pipeline\n", "Exported the following quantities to /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.h5\n", "['mul', 'glat', 'vr', 'mub', 'glon']\n", "Overwritten the following quantities to /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.h5\n", "['rad']\n", "Running convert_galactic_to_icrs post-processing pipeline\n", "Exported the following quantities to /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.h5\n", "['dec', 'mudec', 'ra', 'mura']\n", "Running last_conversions post-processing pipeline\n", "Exported the following quantities to /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.h5\n", "['pi']\n", "Overwritten the following quantities to /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.h5\n", "['teff', 'lum']\n" ] }, { "data": { "text/html": [ "
[04/11/26 21:22:18] WARNING  could not close memmap for column                                   dataset_mmap.py:94\n",
       "                             /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.h5                          \n",
       "
\n" ], "text/plain": [ "\u001b[2;36m[04/11/26 21:22:18]\u001b[0m\u001b[2;36m \u001b[0m\u001b[33mWARNING \u001b[0m could not close memmap for column \u001b]8;id=4642755;file:///home/athob/miniconda3/envs/galaxia_kernels/lib/python3.12/site-packages/vaex/dataset_mmap.py\u001b\\\u001b[2mdataset_mmap.py\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=4642756;file:///home/athob/miniconda3/envs/galaxia_kernels/lib/python3.12/site-packages/vaex/dataset_mmap.py#94\u001b\\\u001b[2m94\u001b[0m\u001b]8;;\u001b\\\n", "\u001b[2;36m \u001b[0m \u001b[35m/home/athob/Software/repos/py-ananke/jupyter/\u001b[0m\u001b[95msurvey.sim.0.h5\u001b[0m \u001b[2m \u001b[0m\n" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "Running observed_magnitudes post-processing pipeline\n", "Exported the following quantities to /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.h5\n", "['gaia__dr2_g_Intrinsic', 'gaia__dr2_grp_Intrinsic', 'gaia__dr2_gbp_Intrinsic']\n", "Overwritten the following quantities to /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.h5\n", "['gaia__dr2_gbp', 'gaia__dr2_grp', 'gaia__dr2_g']\n" ] }, { "data": { "text/html": [ "
                    WARNING  could not close memmap for column                                   dataset_mmap.py:94\n",
       "                             /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.h5                          \n",
       "
\n" ], "text/plain": [ "\u001b[2;36m \u001b[0m\u001b[2;36m \u001b[0m\u001b[33mWARNING \u001b[0m could not close memmap for column \u001b]8;id=4642761;file:///home/athob/miniconda3/envs/galaxia_kernels/lib/python3.12/site-packages/vaex/dataset_mmap.py\u001b\\\u001b[2mdataset_mmap.py\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=4642762;file:///home/athob/miniconda3/envs/galaxia_kernels/lib/python3.12/site-packages/vaex/dataset_mmap.py#94\u001b\\\u001b[2m94\u001b[0m\u001b]8;;\u001b\\\n", "\u001b[2;36m \u001b[0m \u001b[35m/home/athob/Software/repos/py-ananke/jupyter/\u001b[0m\u001b[95msurvey.sim.0.h5\u001b[0m \u001b[2m \u001b[0m\n" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "Running extinctions post-processing pipeline\n", "Preparing interpolator\n", "Now parallelizing extinctions pipeline\n", "Exported the following quantities to /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.h5\n", "['A_gaia__dr2_grp', 'E(B-V)', 'A_0', 'A_gaia__dr2_gbp', 'A_gaia__dr2_g', 'log10_NH']\n", "Overwritten the following quantities to /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.h5\n", "['gaia__dr2_gbp', 'gaia__dr2_grp', 'gaia__dr2_g']\n" ] }, { "data": { "text/html": [ "
[04/11/26 21:22:29] WARNING  could not close memmap for column                                   dataset_mmap.py:94\n",
       "                             /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.h5                          \n",
       "
\n" ], "text/plain": [ "\u001b[2;36m[04/11/26 21:22:29]\u001b[0m\u001b[2;36m \u001b[0m\u001b[33mWARNING \u001b[0m could not close memmap for column \u001b]8;id=4642767;file:///home/athob/miniconda3/envs/galaxia_kernels/lib/python3.12/site-packages/vaex/dataset_mmap.py\u001b\\\u001b[2mdataset_mmap.py\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=4642768;file:///home/athob/miniconda3/envs/galaxia_kernels/lib/python3.12/site-packages/vaex/dataset_mmap.py#94\u001b\\\u001b[2m94\u001b[0m\u001b]8;;\u001b\\\n", "\u001b[2;36m \u001b[0m \u001b[35m/home/athob/Software/repos/py-ananke/jupyter/\u001b[0m\u001b[95msurvey.sim.0.h5\u001b[0m \u001b[2m \u001b[0m\n" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stderr", "output_type": "stream", "text": [ "/home/athob/miniconda3/envs/galaxia_kernels/lib/python3.12/site-packages/ananke/_default_error_model.py:51: RuntimeWarning: overflow encountered in power\n", " grvs = rpmag + 132.32 - 377.28*ggrp + 402.32*ggrp**2 - 190.97*ggrp**3 + 34.026*ggrp**4\n", "/home/athob/miniconda3/envs/galaxia_kernels/lib/python3.12/site-packages/ananke/_default_error_model.py:51: RuntimeWarning: overflow encountered in multiply\n", " grvs = rpmag + 132.32 - 377.28*ggrp + 402.32*ggrp**2 - 190.97*ggrp**3 + 34.026*ggrp**4\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "Running error_modeling post-processing pipeline\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/athob/miniconda3/envs/galaxia_kernels/lib/python3.12/site-packages/ananke/utils.py:45: UserWarning: Pandas doesn't allow columns to be created via a new attribute name - see https://pandas.pydata.org/pandas-docs/stable/indexing.html#attribute-access\n", " self._record_of_all_used_keys = set()\n", "/home/athob/miniconda3/envs/galaxia_kernels/lib/python3.12/site-packages/ananke/ErrorModelDriver.py:137: FutureWarning: The behavior of Series.argmax/argmin with skipna=False and NAs, or with all-NAs is deprecated. In a future version this will raise ValueError.\n", " i_max_err = np.abs(df[prop_err_name] if isinstance(df, pd.DataFrame) else df[prop_err_name].to_pandas_series()).argmax()\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "Exported the following quantities to /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.h5\n", "['mudec_Clean', 'gaia__dr2_gbp_Sig', 'ra_Sig', 'mudec_Sig', 'vr_Err', 'pi_Sig', 'pi_Clean', 'gaia__dr2_g_Err', 'gaia__dr2_grp_Sig', 'dec_Sig', 'mura_Err', 'dec_Clean', 'ra_Err', 'gaia__dr2_gbp_Err', 'mura_Sig', 'vr_Clean', 'dec_Err', 'gaia__dr2_grp_Err', 'pi_Err', 'ra_Clean', 'gaia__dr2_g_Sig', 'vr_Sig', 'mudec_Err', 'mura_Clean']\n", "Overwritten the following quantities to /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.h5\n", "['gaia__dr2_gbp', 'mura', 'gaia__dr2_grp', 'pi', 'dec', 'gaia__dr2_g', 'ra', 'vr', 'mudec']\n" ] }, { "data": { "text/html": [ "
[04/11/26 21:22:30] WARNING  could not close memmap for column                                   dataset_mmap.py:94\n",
       "                             /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.h5                          \n",
       "
\n" ], "text/plain": [ "\u001b[2;36m[04/11/26 21:22:30]\u001b[0m\u001b[2;36m \u001b[0m\u001b[33mWARNING \u001b[0m could not close memmap for column \u001b]8;id=4642773;file:///home/athob/miniconda3/envs/galaxia_kernels/lib/python3.12/site-packages/vaex/dataset_mmap.py\u001b\\\u001b[2mdataset_mmap.py\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=4642774;file:///home/athob/miniconda3/envs/galaxia_kernels/lib/python3.12/site-packages/vaex/dataset_mmap.py#94\u001b\\\u001b[2m94\u001b[0m\u001b]8;;\u001b\\\n", "\u001b[2;36m \u001b[0m \u001b[35m/home/athob/Software/repos/py-ananke/jupyter/\u001b[0m\u001b[95msurvey.sim.0.h5\u001b[0m \u001b[2m \u001b[0m\n" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "Running convert_icrs_to_galactic post-processing pipeline\n", "Overwritten the following quantities to /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.h5\n", "['glat', 'mub', 'mul', 'glon']\n" ] }, { "data": { "text/html": [ "
[04/11/26 21:22:35] WARNING  could not close memmap for column                                   dataset_mmap.py:94\n",
       "                             /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.h5                          \n",
       "
\n" ], "text/plain": [ "\u001b[2;36m[04/11/26 21:22:35]\u001b[0m\u001b[2;36m \u001b[0m\u001b[33mWARNING \u001b[0m could not close memmap for column \u001b]8;id=4642779;file:///home/athob/miniconda3/envs/galaxia_kernels/lib/python3.12/site-packages/vaex/dataset_mmap.py\u001b\\\u001b[2mdataset_mmap.py\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=4642780;file:///home/athob/miniconda3/envs/galaxia_kernels/lib/python3.12/site-packages/vaex/dataset_mmap.py#94\u001b\\\u001b[2m94\u001b[0m\u001b]8;;\u001b\\\n", "\u001b[2;36m \u001b[0m \u001b[35m/home/athob/Software/repos/py-ananke/jupyter/\u001b[0m\u001b[95msurvey.sim.0.h5\u001b[0m \u001b[2m \u001b[0m\n" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "survey = ananke.run()" ] }, { "attachments": {}, "cell_type": "markdown", "id": "1664586c", "metadata": {}, "source": [ "The output is saved as a `vaex` dataframe, with its columns organized in alphabetical order. The photometry is stored in columns that are named according to their associated photometric system and filter name. The keys are in the lowercase format `photosys_filtername` where `photo_sys` is the photometric system and `filtername` is the filter name. In this example, the apparent photometry in filters `gbp`, `grp` & `g` of the Gaia DR2 system (identified as `GAIA__DR2` in Galaxia) are respectively under keys `gaia__dr2_gbp`, `gaia__dr2_grp` & `gaia__dr2_g`. Beside these, the columns include:\n", "\n", "- key `A_0` for the reference extinction which extinction coefficients are based on (at $\\lambda = 550 \\, nm$ in the case of Gaia DR2)\n", "- key `A_{filter_name}` for the extinction value in each filter designated by `filter_name` (here `A_gaia__dr2_g`, `A_gaia__dr2_gbp` & `A_gaia__dr2_grp`)\n", "- key `E(B-V)` for the reddening index\n", "- key `age` for the $log_{10}$ stellar age in years\n", "- key `alpha`, `calcium`, `carbon`, `helium`, `magnesium`, `neon`, `nitrogen`, `oxygen`, `silicon`, `sulphur` for the various chemical abundances as given as input\n", "- key `dec`, `ra` for the astrometric declination and right ascension celestial coordinates in degrees\n", "- key `dform` for the formation distance as given as input\n", "- key `dmod` for the distance modulus\n", "- key `feh` for the stellar metallicity \\[Fe/H\\] in dex relative to solar\n", "- key `{filter_name}_Intrinsic` for the intrinsic stellar photometry in each filter designated by `filter_name` (here `gaia__dr2_g_Intrinsic`, `gaia__dr2_gbp_Intrinsic` & `gaia__dr2_grp_Intrinsic`)\n", "- key `glat`, `glon` for the astrometric galactic latitude and longitude celestrial coordinates in degrees\n", "- key `grav` for the $log_{10}$ surface gravity in CGS units\n", "- key `log10_NH` for the $log_{10}$ hydrogen column density between Observer position and star in $cm^{-2}$\n", "- key `lum` for the stellar luminosity in solar luminosities\n", "- key `mact`, `mtip`, `smass` for respectively the current stellar mass, the mass of that same star at tip of giant branch for its given age \\& metallicity and its stellar mass on zero-age main sequence, all in solar masses\n", "- key `mub`, `mudec`, `mul`, `mura` for the astrometric proper motions, respectively in the direction of the galactic latitude, declination, galactic longitude and right ascension, all in milliarcseconds per year\n", "- key `parentid` for the parent particle index as given as input\n", "- key `partid` for the flag that identifies stars that are *not* central relatively to their parent particle\n", "- key `pi` for the star parallax in milliarcseconds\n", "- key `px`, `py`, `pz` for the star position cartesian coordinates in $kpc$ relative to the Observer's position\n", "- key `rad` for the star distance to the Observer in $kpc$\n", "- key `satid` for the index ``id`` the population the parent particle is part of\n", "- key `teff` for the star effective temperature in Kelvin\n", "- key `vr` for the star astrometric radial velocity in $km.s^{-1}$\n", "- key `vx`, `vy`, `vz` for the star velocity cartesian coordinates in $km.s^{-1}$ relative to the Observer's velocity\n", "\n", "Additionally, astrometric and photometric quantities `X` all have associated columns identified as:\n", "- key `{X}_Sig` for the standard error on the quantity `X`\n", "- key `{X}_Err` for the actual drawn gaussian error on the quantity `X`\n", "\n", "Taking the photometry as an example of the latter, the standard errors are stored under `gaia__dr2_g_Sig`, `gaia__dr2_gbp_Sig` & `gaia__dr2_grp_Sig`, while the drawn errors are under `gaia__dr2_g_Err`, `gaia__dr2_gbp_Err` & `gaia__dr2_grp_Err`." ] }, { "cell_type": "code", "execution_count": 6, "id": "42e393f6", "metadata": {}, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
# A_0 A_gaia__dr2_g A_gaia__dr2_gbp A_gaia__dr2_grp E(B-V) age alpha calcium carbon dec dec_Clean dec_Err dec_Sig dform dmod feh gaia__dr2_g gaia__dr2_g_Err gaia__dr2_g_Intrinsic gaia__dr2_g_Sig gaia__dr2_gbp gaia__dr2_gbp_Err gaia__dr2_gbp_Intrinsic gaia__dr2_gbp_Sig gaia__dr2_grp gaia__dr2_grp_Err gaia__dr2_grp_Intrinsic gaia__dr2_grp_Sig glat glon grav helium log10_NH lum mact magnesium mtip mub mudec mudec_Clean mudec_Err mudec_Sig mul mura mura_Clean mura_Err mura_Sig neon nitrogen oxygen parentid partid partitionid pi pi_Clean pi_Err pi_Sig px py pz ra ra_Clean ra_Err ra_Sig rad satid silicon smass sulphur teff vr vr_Clean vr_Err vr_Sig vx vy vz
0 0.709397791159409 0.52206985445396150.70919313390714510.4055614050521911 0.2288379971481964410.0163259506225590.3380742371082306-1.1830222606658936-0.39430123567581177nan 43.89190561692206 nan nan 0.0 20.87123488088669 -0.8867433071136475nan nan 7.8221183 nan nan nan 8.470699 nan nan nan 7.0639534 nan nan nan 4.769691 -0.7315406799316406 5.72094992870491e+21 0.0710139650.50516033-0.5486690402030945 0.886202 nan nan -0.2413280326361944 nan nan nan nan 0.29630942610455613nan nan -0.7326117157936096-0.652482807636261-0.66972786188125610 0 0 nan 0.006695037652823928 nan nan -142.0352559796068 46.20841677042607 0.6925918281909844 nan 73.8027984762779 nan nan 149.3643578805275 0 -0.75635123252868650.5054699889829319 -1.09243309497833254285.4614-156.7385192430598 -156.7385192430598 nan nan 67.38685760498046 -300.3354850769043 55.13861343383789
1 0.75668220632304230.55145941111718840.75344249911562890.431195670717934460.2440910342977555810.0163259506225590.3380742371082306-1.1830222606658936-0.39430123567581177nan 43.735186299167324nan nan 0.0 20.87527071994942 -0.8867433071136475nan nan 7.993824 nan nan nan 8.668457 nan nan nan 7.2174444 nan nan nan 4.7857795-0.7315406799316406 6.102275857443889e+21 0.06093182 0.48871416-0.5486690402030945 0.886202 nan nan -0.24008076703265266nan nan nan nan 0.29130663692538156nan nan -0.7326117157936096-0.652482807636261-0.66972786188125610 1 0 nan 0.006682605995616205 nan nan -142.2608617629229346.417971376138425-0.11534282804241658nan 73.42607530130886nan nan 149.6422205133747 0 -0.75635123252868650.48905900859965035-1.09243309497833254198.5107-156.97767104891037-156.97767104891037nan nan 67.7484200428451 -298.2995731452656452.69162146715175
2 136.6464782405509 2904.23239167338441370.444979191365778.02675978260835 44.079509109855124 10.0163259506225590.3380742371082306-1.1830222606658936-0.39430123567581177nan 44.364568278584706nan nan 0.0 20.888157433594138-0.8867433071136475nan nan 8.71385 nan nan nan 9.517461 nan nan nan 7.8597684 nan nan nan 4.8554068-0.7315406799316406 1.101987727746378e+24 0.03486588 0.42435014-0.5486690402030945 0.886202 nan nan -0.23863717453488725nan nan nan nan 0.2958006234265731 nan nan -0.7326117157936096-0.652482807636261-0.66972786188125610 1 0 nan 0.006643065174432654 nan nan -142.9270461828407447.20842184009032 1.8396875936103667 nan 74.0498933920498 nan nan 150.532920232173440 -0.75635123252868650.4250665825385847 -1.09243309497833253935.5205-158.4535061010785 -158.4535061010785 nan nan 68.03093236871236 -301.4758806839234 56.09200627357448
3 0.71509160840002480.59382870248778010.75988277212997570.422318035751309160.2306747123871047710.0163259506225590.3380742371082306-1.1830222606658936-0.39430123567581177nan 43.89621033939353 nan nan 0.0 20.86868071849355 -0.8867433071136475nan nan 5.1929526 nan nan nan 5.4745564 nan nan nan 4.7391353 nan nan nan 4.5336547-0.7315406799316406 5.766867809677618e+21 0.6409676 0.74067265-0.5486690402030945 0.886202 nan nan -0.24429902665121833nan nan nan nan 0.2940592689423017 nan nan -0.7326117157936096-0.652482807636261-0.66972786188125610 1 0 nan 0.006702917225320304 nan nan -141.8805554620098 46.1152954631301 0.7600375956495602 nan 73.8447277872187 nan nan 149.188773542136970 -0.75635123252868650.7407077517512485 -1.09243309497833255891.9644-156.7353025299949 -156.7353025299949 nan nan 67.38594209547328 -300.6020691055117 52.577100775679746
4 1.18695442936753710.81358271517667041.15634729390904580.6614930768091819 0.3828885256024313 10.0163259506225590.3380742371082306-1.1830222606658936-0.39430123567581177nan 43.93297434426074 nan nan 0.0 20.879530213898438-0.8867433071136475nan nan 8.90345 nan nan nan 9.740426 nan nan nan 8.03157 nan nan nan 4.877773 -0.7315406799316406 9.572213140060782e+21 0.02983189 0.4045688 -0.5486690402030945 0.886202 nan nan -0.24092819622063041nan nan nan nan 0.296377068576234 nan nan -0.7326117157936096-0.652482807636261-0.66972786188125610 1 0 nan 0.006669510447933342 nan nan -142.5205780120727246.5650255517057160.6324477943182331 nan 73.713875054752 nan nan 149.9360422037972 0 -0.75635123252868650.4050935914883408 -1.09243309497833253882.164 -156.2199566593424 -156.2199566593424 nan nan 66.21987999393174 -301.0910757394999655.35190950080388
... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...
1,611,7079.849335135566943 5.905034450988529 8.998823684092358 5.115782383513528 3.177204882440949 9.99021053314209 0.7036365270614624-0.89319157600402830.053123194724321365nan 52.47591875749927 nan nan 0.0 21.61297084448185 -0.8060703873634338nan nan 6.9566293 nan nan nan 7.4624496 nan nan nan 6.302246 nan nan nan 4.7021213-0.051267735660076147.943012206102372e+22 0.1396583 0.5802463 -0.102433867752552030.8955961nan nan -0.1151919563604325 nan nan nan nan 0.08138204404640034nan nan -0.56209683418273930.0669776126742363-1.214622616767883399999 1 0 nan 0.004757796156067555 nan nan -186.3597503860403 96.74107735172117 9.349095600452483 nan 67.0515975296681 nan nan 210.181345984046230 -1.247308373451233 0.5813063179170757 -0.99559468030929574714.454 -90.2409497740284 -90.2409497740284 nan nan 15.159623414379855-164.46277627244382-24.763303315808546
1,611,70813.3052868969287657.972572810389155 12.1990440944498646.734497001850922 4.292028031267344 9.99021053314209 0.7036365270614624-0.89319157600402830.053123194724321365nan 52.10502299677406 nan nan 0.0 21.608894668241703-0.8060703873634338nan nan 5.3323083 nan nan nan 5.6339135 nan nan nan 4.8603163 nan nan nan 4.551513 -0.051267735660076141.0730070078168358e+230.5618677 0.73399246-0.102433867752552030.8955961nan nan -0.11464220238473403nan nan nan nan 0.08162295523810104nan nan -0.56209683418273930.0669776126742363-1.214622616767883399999 1 0 nan 0.0047667356338638904nan nan -186.0710412860123 96.61308355256259 7.357869720071434 nan 66.41106254567475nan nan 209.787174454523980 -1.247308373451233 0.7344738879028191 -0.99559468030929575773.816 -90.25301122452501 -90.25301122452501 nan nan 15.649286825866614-163.9671077947139 -24.555644142144015
1,611,70918.18659863039707511.42519172980502916.50377702486021 8.639310462536832 5.866644719482927 9.99021053314209 0.7036365270614624-0.89319157600402830.053123194724321365nan 52.281639500974386nan nan 0.0 21.614814425414835-0.8060703873634338nan nan 6.256601 nan nan nan 6.6638975 nan nan nan 5.6850376 nan nan nan 4.654846 -0.051267735660076141.4666611798707317e+230.2484262 0.6444424 -0.102433867752552030.8955961nan nan -0.11541291785371721nan nan nan nan 0.08405587899877405nan nan -0.56209683418273930.0669776126742363-1.214622616767883399999 1 0 nan 0.0047537584994812 nan nan -186.401187850548 97.17315223354157 7.952910951827339 nan 66.48740383323538nan nan 210.359865800741630 -1.247308373451233 0.6446346203563409 -0.99559468030929575160.5283-89.2744102255874 -89.2744102255874 nan nan 13.274847794294757-165.87604067238016-23.46408831673764
1,611,71017.80165910608414 10.99406360503620316.2662064130328558.558298042503163 5.742470679381981 9.99021053314209 0.7036365270614624-0.89319157600402830.053123194724321365nan 52.60308949195131 nan nan 0.0 21.6255893004592 -0.8060703873634338nan nan 5.077172 nan nan nan 5.3556504 nan nan nan 4.6276927 nan nan nan 4.507982 -0.051267735660076141.435617669845495e+23 0.71100533 0.7585091 -0.102433867752552030.8955961nan nan -0.11182195701593674nan nan nan nan 0.08124741721522437nan nan -0.56209683418273930.0669776126742363-1.214622616767883399999 1 0 nan 0.004730228712270713 nan nan -186.3828534065042799.48057937118065 7.593190154300342 nan 65.72063877969379nan nan 211.406268243625160 -1.247308373451233 0.7591386913124265 -0.99559468030929575922.2593-91.18185681851935 -91.18185681851935 nan nan 15.256838439138757-163.3252054589539 -24.379336029180863
1,611,7118.816373078046517 5.566795304062381 8.328230058759424 4.706684597638085 2.8439913154988763 9.99021053314209 0.7036365270614624-0.89319157600402830.053123194724321365nan 52.47916141482985 nan nan 0.0 21.613667682089748-0.8060703873634338nan nan 5.4040017 nan nan nan 5.7124815 nan nan nan 4.925138 nan nan nan 4.5628715-0.051267735660076147.1099782887471904e+220.52572614 0.7271273 -0.102433867752552030.8955961nan nan -0.1148307999308442 nan nan nan nan 0.08279535335029264nan nan -0.56209683418273930.0669776126742363-1.214622616767883399999 1 0 nan 0.004756269597697463 nan nan -186.4613449414160396.6796112442321 9.475214471500284 nan 67.12279484473645nan nan 210.2488051737239 0 -1.247308373451233 0.7273869656881858 -0.99559468030929575729.2954-88.7516720457417 -88.7516720457417 nan nan 13.603362480412988-164.48282279687962-23.355289452640978
" ], "text/plain": [ "# A_0 A_gaia__dr2_g A_gaia__dr2_gbp A_gaia__dr2_grp E(B-V) age alpha calcium carbon dec dec_Clean dec_Err dec_Sig dform dmod feh gaia__dr2_g gaia__dr2_g_Err gaia__dr2_g_Intrinsic gaia__dr2_g_Sig gaia__dr2_gbp gaia__dr2_gbp_Err gaia__dr2_gbp_Intrinsic gaia__dr2_gbp_Sig gaia__dr2_grp gaia__dr2_grp_Err gaia__dr2_grp_Intrinsic gaia__dr2_grp_Sig glat glon grav helium log10_NH lum mact magnesium mtip mub mudec mudec_Clean mudec_Err mudec_Sig mul mura mura_Clean mura_Err mura_Sig neon nitrogen oxygen parentid partid partitionid pi pi_Clean pi_Err pi_Sig px py pz ra ra_Clean ra_Err ra_Sig rad satid silicon smass sulphur teff vr vr_Clean vr_Err vr_Sig vx vy vz\n", "0 0.709397791159409 0.5220698544539615 0.7091931339071451 0.4055614050521911 0.22883799714819644 10.016325950622559 0.3380742371082306 -1.1830222606658936 -0.39430123567581177 nan 43.89190561692206 nan nan 0.0 20.87123488088669 -0.8867433071136475 nan nan 7.8221183 nan nan nan 8.470699 nan nan nan 7.0639534 nan nan nan 4.769691 -0.7315406799316406 5.72094992870491e+21 0.071013965 0.50516033 -0.5486690402030945 0.886202 nan nan -0.2413280326361944 nan nan nan nan 0.29630942610455613 nan nan -0.7326117157936096 -0.652482807636261 -0.6697278618812561 0 0 0 nan 0.006695037652823928 nan nan -142.0352559796068 46.20841677042607 0.6925918281909844 nan 73.8027984762779 nan nan 149.3643578805275 0 -0.7563512325286865 0.5054699889829319 -1.0924330949783325 4285.4614 -156.7385192430598 -156.7385192430598 nan nan 67.38685760498046 -300.3354850769043 55.13861343383789\n", "1 0.7566822063230423 0.5514594111171884 0.7534424991156289 0.43119567071793446 0.24409103429775558 10.016325950622559 0.3380742371082306 -1.1830222606658936 -0.39430123567581177 nan 43.735186299167324 nan nan 0.0 20.87527071994942 -0.8867433071136475 nan nan 7.993824 nan nan nan 8.668457 nan nan nan 7.2174444 nan nan nan 4.7857795 -0.7315406799316406 6.102275857443889e+21 0.06093182 0.48871416 -0.5486690402030945 0.886202 nan nan -0.24008076703265266 nan nan nan nan 0.29130663692538156 nan nan -0.7326117157936096 -0.652482807636261 -0.6697278618812561 0 1 0 nan 0.006682605995616205 nan nan -142.26086176292293 46.417971376138425 -0.11534282804241658 nan 73.42607530130886 nan nan 149.6422205133747 0 -0.7563512325286865 0.48905900859965035 -1.0924330949783325 4198.5107 -156.97767104891037 -156.97767104891037 nan nan 67.7484200428451 -298.29957314526564 52.69162146715175\n", "2 136.6464782405509 2904.2323916733844 1370.4449791913657 78.02675978260835 44.079509109855124 10.016325950622559 0.3380742371082306 -1.1830222606658936 -0.39430123567581177 nan 44.364568278584706 nan nan 0.0 20.888157433594138 -0.8867433071136475 nan nan 8.71385 nan nan nan 9.517461 nan nan nan 7.8597684 nan nan nan 4.8554068 -0.7315406799316406 1.101987727746378e+24 0.03486588 0.42435014 -0.5486690402030945 0.886202 nan nan -0.23863717453488725 nan nan nan nan 0.2958006234265731 nan nan -0.7326117157936096 -0.652482807636261 -0.6697278618812561 0 1 0 nan 0.006643065174432654 nan nan -142.92704618284074 47.20842184009032 1.8396875936103667 nan 74.0498933920498 nan nan 150.53292023217344 0 -0.7563512325286865 0.4250665825385847 -1.0924330949783325 3935.5205 -158.4535061010785 -158.4535061010785 nan nan 68.03093236871236 -301.4758806839234 56.09200627357448\n", "3 0.7150916084000248 0.5938287024877801 0.7598827721299757 0.42231803575130916 0.23067471238710477 10.016325950622559 0.3380742371082306 -1.1830222606658936 -0.39430123567581177 nan 43.89621033939353 nan nan 0.0 20.86868071849355 -0.8867433071136475 nan nan 5.1929526 nan nan nan 5.4745564 nan nan nan 4.7391353 nan nan nan 4.5336547 -0.7315406799316406 5.766867809677618e+21 0.6409676 0.74067265 -0.5486690402030945 0.886202 nan nan -0.24429902665121833 nan nan nan nan 0.2940592689423017 nan nan -0.7326117157936096 -0.652482807636261 -0.6697278618812561 0 1 0 nan 0.006702917225320304 nan nan -141.8805554620098 46.1152954631301 0.7600375956495602 nan 73.8447277872187 nan nan 149.18877354213697 0 -0.7563512325286865 0.7407077517512485 -1.0924330949783325 5891.9644 -156.7353025299949 -156.7353025299949 nan nan 67.38594209547328 -300.6020691055117 52.577100775679746\n", "4 1.1869544293675371 0.8135827151766704 1.1563472939090458 0.6614930768091819 0.3828885256024313 10.016325950622559 0.3380742371082306 -1.1830222606658936 -0.39430123567581177 nan 43.93297434426074 nan nan 0.0 20.879530213898438 -0.8867433071136475 nan nan 8.90345 nan nan nan 9.740426 nan nan nan 8.03157 nan nan nan 4.877773 -0.7315406799316406 9.572213140060782e+21 0.02983189 0.4045688 -0.5486690402030945 0.886202 nan nan -0.24092819622063041 nan nan nan nan 0.296377068576234 nan nan -0.7326117157936096 -0.652482807636261 -0.6697278618812561 0 1 0 nan 0.006669510447933342 nan nan -142.52057801207272 46.565025551705716 0.6324477943182331 nan 73.713875054752 nan nan 149.9360422037972 0 -0.7563512325286865 0.4050935914883408 -1.0924330949783325 3882.164 -156.2199566593424 -156.2199566593424 nan nan 66.21987999393174 -301.09107573949996 55.35190950080388\n", "... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...\n", "1,611,707 9.849335135566943 5.905034450988529 8.998823684092358 5.115782383513528 3.177204882440949 9.99021053314209 0.7036365270614624 -0.8931915760040283 0.053123194724321365 nan 52.47591875749927 nan nan 0.0 21.61297084448185 -0.8060703873634338 nan nan 6.9566293 nan nan nan 7.4624496 nan nan nan 6.302246 nan nan nan 4.7021213 -0.05126773566007614 7.943012206102372e+22 0.1396583 0.5802463 -0.10243386775255203 0.8955961 nan nan -0.1151919563604325 nan nan nan nan 0.08138204404640034 nan nan -0.5620968341827393 0.0669776126742363 -1.2146226167678833 99999 1 0 nan 0.004757796156067555 nan nan -186.3597503860403 96.74107735172117 9.349095600452483 nan 67.0515975296681 nan nan 210.18134598404623 0 -1.247308373451233 0.5813063179170757 -0.9955946803092957 4714.454 -90.2409497740284 -90.2409497740284 nan nan 15.159623414379855 -164.46277627244382 -24.763303315808546\n", "1,611,708 13.305286896928765 7.972572810389155 12.199044094449864 6.734497001850922 4.292028031267344 9.99021053314209 0.7036365270614624 -0.8931915760040283 0.053123194724321365 nan 52.10502299677406 nan nan 0.0 21.608894668241703 -0.8060703873634338 nan nan 5.3323083 nan nan nan 5.6339135 nan nan nan 4.8603163 nan nan nan 4.551513 -0.05126773566007614 1.0730070078168358e+23 0.5618677 0.73399246 -0.10243386775255203 0.8955961 nan nan -0.11464220238473403 nan nan nan nan 0.08162295523810104 nan nan -0.5620968341827393 0.0669776126742363 -1.2146226167678833 99999 1 0 nan 0.0047667356338638904 nan nan -186.0710412860123 96.61308355256259 7.357869720071434 nan 66.41106254567475 nan nan 209.78717445452398 0 -1.247308373451233 0.7344738879028191 -0.9955946803092957 5773.816 -90.25301122452501 -90.25301122452501 nan nan 15.649286825866614 -163.9671077947139 -24.555644142144015\n", "1,611,709 18.186598630397075 11.425191729805029 16.50377702486021 8.639310462536832 5.866644719482927 9.99021053314209 0.7036365270614624 -0.8931915760040283 0.053123194724321365 nan 52.281639500974386 nan nan 0.0 21.614814425414835 -0.8060703873634338 nan nan 6.256601 nan nan nan 6.6638975 nan nan nan 5.6850376 nan nan nan 4.654846 -0.05126773566007614 1.4666611798707317e+23 0.2484262 0.6444424 -0.10243386775255203 0.8955961 nan nan -0.11541291785371721 nan nan nan nan 0.08405587899877405 nan nan -0.5620968341827393 0.0669776126742363 -1.2146226167678833 99999 1 0 nan 0.0047537584994812 nan nan -186.401187850548 97.17315223354157 7.952910951827339 nan 66.48740383323538 nan nan 210.35986580074163 0 -1.247308373451233 0.6446346203563409 -0.9955946803092957 5160.5283 -89.2744102255874 -89.2744102255874 nan nan 13.274847794294757 -165.87604067238016 -23.46408831673764\n", "1,611,710 17.80165910608414 10.994063605036203 16.266206413032855 8.558298042503163 5.742470679381981 9.99021053314209 0.7036365270614624 -0.8931915760040283 0.053123194724321365 nan 52.60308949195131 nan nan 0.0 21.6255893004592 -0.8060703873634338 nan nan 5.077172 nan nan nan 5.3556504 nan nan nan 4.6276927 nan nan nan 4.507982 -0.05126773566007614 1.435617669845495e+23 0.71100533 0.7585091 -0.10243386775255203 0.8955961 nan nan -0.11182195701593674 nan nan nan nan 0.08124741721522437 nan nan -0.5620968341827393 0.0669776126742363 -1.2146226167678833 99999 1 0 nan 0.004730228712270713 nan nan -186.38285340650427 99.48057937118065 7.593190154300342 nan 65.72063877969379 nan nan 211.40626824362516 0 -1.247308373451233 0.7591386913124265 -0.9955946803092957 5922.2593 -91.18185681851935 -91.18185681851935 nan nan 15.256838439138757 -163.3252054589539 -24.379336029180863\n", "1,611,711 8.816373078046517 5.566795304062381 8.328230058759424 4.706684597638085 2.8439913154988763 9.99021053314209 0.7036365270614624 -0.8931915760040283 0.053123194724321365 nan 52.47916141482985 nan nan 0.0 21.613667682089748 -0.8060703873634338 nan nan 5.4040017 nan nan nan 5.7124815 nan nan nan 4.925138 nan nan nan 4.5628715 -0.05126773566007614 7.1099782887471904e+22 0.52572614 0.7271273 -0.10243386775255203 0.8955961 nan nan -0.1148307999308442 nan nan nan nan 0.08279535335029264 nan nan -0.5620968341827393 0.0669776126742363 -1.2146226167678833 99999 1 0 nan 0.004756269597697463 nan nan -186.46134494141603 96.6796112442321 9.475214471500284 nan 67.12279484473645 nan nan 210.2488051737239 0 -1.247308373451233 0.7273869656881858 -0.9955946803092957 5729.2954 -88.7516720457417 -88.7516720457417 nan nan 13.603362480412988 -164.48282279687962 -23.355289452640978" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "survey" ] }, { "cell_type": "markdown", "id": "a9a28a5d", "metadata": {}, "source": [ "Please refer to [`vaex`'s documentation](https://vaex.io/docs/tutorial.html) for further help on how to use `vaex` dataframes: the following line for example isolate only the rows with non-NaN photometry." ] }, { "cell_type": "code", "execution_count": 7, "id": "09f840d5", "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
# A_0 A_gaia__dr2_g A_gaia__dr2_gbp A_gaia__dr2_grp E(B-V) age alpha calcium carbon dec dec_Clean dec_Err dec_Sig dform dmod feh gaia__dr2_g gaia__dr2_g_Err gaia__dr2_g_Intrinsic gaia__dr2_g_Sig gaia__dr2_gbp gaia__dr2_gbp_Err gaia__dr2_gbp_Intrinsic gaia__dr2_gbp_Sig gaia__dr2_grp gaia__dr2_grp_Err gaia__dr2_grp_Intrinsic gaia__dr2_grp_Sig glat glon grav helium log10_NH lum mact magnesium mtip mub mudec mudec_Clean mudec_Err mudec_Sig mul mura mura_Clean mura_Err mura_Sig neon nitrogen oxygen parentid partid partitionid pi pi_Clean pi_Err pi_Sig px py pz ra ra_Clean ra_Err ra_Sig rad satid silicon smass sulphur teff vr vr_Clean vr_Err vr_Sig vx vy vz
0 0.358846324020237750.262481850769453460.3580179075913724 0.204661728531405350.115756878716205739.6928701400756840.3099161982536316 -1.1548902988433838-0.4665657579898834 18.86801579630566318.8680157785651871.7740477296505672e-08 1.691274e-07 0.0 21.118883008502078-0.4829272329807281519.868202 0.012658170630218355 -1.5258217 0.009227399 20.504751 -0.16132866071418359 -0.8108196 0.17850117 18.7721 -0.22652792657603063 -2.3249161 0.17850117 -14.991976836287705-177.956902949767941.2548169 -0.278090775012969972.893921967905143e+21 448.282441.1287807 -0.173011064529418951.1379396 -0.15521054578348675-1.3589967238651148 -0.3270482743773745 -1.0319484494877402 1.0656426 1.570216570810097 0.8017470282407725 0.1632684267596954 0.638478601481077 1.0656426 -0.7034189105033875 -0.09953558444976807-0.4723181426525116 49 1 0 -0.6145466549364363 0.005973424767124436 -0.6205200797035608 0.60885864-161.60712291393943-5.765154922998849-43.30577433967765 74.02626514013622 74.026265027250491.1288572264651892e-07 1.691274e-07 167.408151769758180 -0.60188096761703491.1408413190663447-0.7073405385017395 4129.926 69.84271348583249 69.84271348583249 nan nan -44.649975087957735-287.75369665394226-65.06155285732873
1 0.5999316419840481 0.3951491437418697 0.5843014320686075 0.327567946448567 0.1935263361238865 9.9881792068481450.117806412279605870.09089594334363937-0.2674128711223602328.32392700064622228.323927151103447-1.5045722692676327e-077.903525e-08 0.0 21.313550824736375-0.6772944927215576 19.078148 0.007681412416285324 -2.6382327 0.005971381 20.259933 -0.07032803519170927 -1.5675896 0.09647448 17.950785 -0.07180883013983862 -3.6185248 0.09647448 -18.536442761608946165.70363972881054 0.3074857 -0.213550269603729254.838158403097162e+21 1852.92180.85137147-0.55948805809021 0.912569340.3131789409812662 -0.1038156537874877 -0.149861869114884920.04604621532739723 0.55013824 0.48286076265325345 0.5660899886174509 0.172859030109131340.393230958508319460.55013824-0.3997773826122284 -0.7002913355827332 -0.08155464380979538144 1 0 -0.152748911588277470.005461222050076329 -0.1582101336383538 0.2845269 -168.2333465832280442.8707615557611 -58.21183715231835 59.86462922169792559.864629159298656.239927571995072e-08 7.903525e-08 183.1091998879671 0 -0.31031048297882080.914106032084624 -0.380739361047744753662.81 -25.762507382985724-25.762507382985724nan nan -27.923177997805908-197.6293933960381316.189820511232778
2 0.176081952679942480.137703394273812170.1808537019520938 0.102556633208450040.056800629896755639.579007148742676-0.6709350943565369-0.9936923980712891-0.9404328465461731 33.04906634892173433.04906634871991 2.0182171922849153e-10 3.6569543e-070.0 21.38615087842112 -0.9346992373466492 20.56803 -0.03600744374280186 -0.9198183 0.014033196 21.098215 -0.07126815067421044 -0.3975206 0.32202673 19.32279 -0.5803682522135298 -1.585548 0.32202673 -8.484585766867529 168.82659578503788 1.813334 -0.152460306882858281.4200157474188907e+21204.612841.146 -1.605634331703186 1.1506394 -0.3324971332482145 -1.1621686215766727 -0.12728573316172523-1.0348828884149475 2.0831733 1.242114141035326 0.550241745054646 0.1341992658590463 0.4160424791955997 2.0831733 -0.9467113018035889 -0.7777135372161865 -0.36034896969795227355 1 0 -2.8237769876652146 0.005281652504745023 -2.8290586401699596 1.3165035 -183.7130045148962636.287543243958545-27.93507061765807 70.77774203758443 70.777741676888433.6069600030252607e-07 3.6569543e-07189.334682488407280 -0.71507573127746581.1526367829812323-0.128577053546905524663.5015-75.65222985360518 -75.65222985360518 nan nan 38.83221747966123 -175.9006125758071 28.874331454692708
3 0.178611749883683950.128920333499283950.177478230315429920.101347475205828720.057616693510865799.481425285339355-1.1020725965499878-1.0413920879364014-0.5413622856140137 28.23936555420345628.23936538065812 1.735453331636329e-07 2.606068e-07 0.0 21.781771243448272-0.2500181496143341 20.272312 -0.006433955642751986 -1.6319455 0.0116908895 21.282545 0.17120013764702185 -0.8479041 0.24911553 19.201895 -0.20884711491488975 -2.4723763 0.24911553 -4.592810387657218 177.86633782029145 1.2065276 -0.8023523092269897 1.4404173377716447e+21537.386 1.3710679 -1.3520907163619995 1.369177 0.8879077655799169 0.7300725836482986 -0.217375011185357930.9474475948336566 1.551452 -0.27902170477556787 0.5772584642238547 0.150781561917897540.426476902305957251.551452 -0.6431713700294495 -0.4919672906398773 -0.6521893739700317 359 1 0 0.4967736437882954 0.004401956548212814 0.4923716872400826 0.93818444-226.285253498574928.430630112877303 -18.19051182758227780.6929868604775 80.69298721775455-3.572770414100613e-07 2.606068e-07 227.171710817090680 -0.75073283910751341.3773316962276039-0.6594099998474121 4003.228 27.33035105089606 27.33035105089606 nan nan -38.03903369359221 -283.6631352491681 0.4139964206713862
4 0.945820473694447 0.6678443096074201 0.930577483029963 0.5330434047235866 0.3051033786111119 9.6038694381713870.19408565759658813-0.6169882416725159-0.7964903116226196 32.35306966051205532.35306979227988 -1.3176782507459781e-071.0227677e-070.0 21.17193983212506 -0.869421124458313 19.34399 -0.0006781485567268116-2.495115 0.0069470913 20.299686 -0.07008269042370177 -1.7327482 0.11953725 18.383387 -3.715706189872538e-05-3.3215604 0.11953725 -6.795830400723491 171.12720912421034 0.82894164-1.036791443824768 7.627584465277797e+21 1142.33591.1365215 -0.6753354668617249 1.1560693 -0.1418502376652223 -0.09612761524104929-0.395926026957800550.29979841171675126 0.6884891 0.010582078436602305 -0.104847278924096340.23107472902875997-0.33592200795285630.6884891 -0.7131249904632568 -1.3325629234313965 -0.8150839805603027 490 1 0 0.33585443984875474 0.005829241304978611 0.33002519854377615 0.36819637-168.3051656639284526.27400588210036 -20.29967534295857274.04671686055353 74.046716797872016.268152637233299e-08 1.0227677e-07171.548911373067480 -0.327743262052536 1.159396049264403 -0.090744525194168094076.2021-19.64436147888908 -19.64436147888908 nan nan -31.620475053741032-368.72625625434716-49.06766686013123
... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...
6740.151018003458936150.118293240257224670.155223255710963570.088000205664986190.048715484986753599.5345153808593751.2677661180496216 -0.4057428538799286-0.59430992603302 21.34353415994132421.34353477135828 -6.11416954624569e-07 3.1462778e-070.0 21.592642702411716-1.1074150800704956 20.454937 -0.005222778643958658 -1.2507769 0.012943381 21.053337 0.03625308061370761 -0.73078275 0.2874478 19.719486 -0.045624436526205674 -1.9155316 0.2874478 -17.503740519827847176.7477689805897 1.6993186 -0.453119516372680661.2178871246688396e+21276.454681.174 0.16035111248493195 1.1788915 -0.20376656041788468-0.11886161920922561-0.2466959595835509 0.12783434037432528 1.8276072 -0.014066707641402212-0.166104180824460880.1330427667679313 -0.29914694759239221.8276072 -0.48284637928009033-0.3429383337497711 -0.361402690410614 99192 1 0 1.1484148018479465 0.004802545183696615 1.14361225666425 1.13266 -198.2618205019605 11.265869808201476-62.62680511726111468.59347486776855 68.593474359692085.080764713807229e-07 3.1462778e-07208.222923835206870 -1.429317593574524 1.1812293030115915-0.138396084308624274680.429 35.33961960826008 35.33961960826008 nan nan -32.91202385958913 -269.96195551855874-61.869312861074555
6750.9635104765596588 0.7031919691501173 0.9598374375111064 0.5493862070334198 0.310809831148277049.3014926910400390.18354283273220062-0.7873296141624451-0.5832746028900146 36.29390846857701536.29390899617845 -5.276014372685639e-07 5.21441e-07 0.0 21.986108149698197-0.6020106673240662 20.938156 -0.00258938495380322 -1.7485539 0.01696236 21.846714 0.0033400161006579407-1.1025709 0.42019537 20.104744 0.07144425513196052 -2.5021954 0.42019537 -18.73755156006982 151.75488224208473 1.3968682 -0.5497167110443115 7.770245778706926e+21 504.3772 1.508 -0.418467819690704351.501947 3.299747393217878 -0.28967745228217223-0.11117319854321678-0.178504253738955452.8375878 5.791327156440196 6.6591208177663495 0.076938645695560446.582182172070789 2.8375878 -0.7871666550636292 -0.6889442801475525 -0.6478202939033508 99272 1 0 -0.6687868125300893 0.004006621987939852 -0.6727934345180292 1.8771876 -208.2157182786129 111.85547750288215-80.17569728370337 47.36247265711304 47.36247304202983-3.849167924330201e-07 5.21441e-07 249.586809788907940 -0.32534617185592651.510235532950821 -0.967689037322998 4293.1353-70.86896352290194 -70.86896352290194 nan nan 7.312308989325453 -170.89071803084343-36.78972701008118
6760.375003121188620840.283381137761839760.3791604979470969 0.216191012780430990.1209687487705228510.16224193572998-0.5301219820976257-0.6149452328681946-0.9359602928161621 20.23864675894622420.23864607329147 6.856547555591803e-07 4.3608154e-070.0 21.28507619772096 -0.466217964887619 20.76201 -0.009498002225838745 -0.7969467 0.015419783 20.97355 -0.4993928131066979 -0.19129108 0.3675803 19.7897 -0.1881952233286074 -1.5233725 0.3675803 -25.702886990754198170.24573776609319 1.5625119 -0.6190789937973022 3.024218719263071e+21 199.889330.792 -0.9963399171829224 0.877864060.42368958053650185 0.9520933879339335 -0.206467772701862951.1585611606357964 2.4282553 -0.9299044207050995 -0.3711512804161766 0.11890487188601179-0.49005615230218842.4282553 -0.6988312005996704 -0.2878437340259552 -0.5345746278762817 99878 0 0 -1.7825581982250396 0.00553330692308911 -1.7880915051481288 1.5698935 -160.487938474724 27.589162658182417-78.38070567241937 57.76421066762154 57.764210030561916.370596314964964e-07 4.3608154e-07180.7237541491236 0 -1.31891858577728270.8813156093858641-0.7903018593788147 4401.795 -4.640823433510391 -4.640823433510391 nan nan -5.309618854522705 -198.32230911254882-48.23520782470703
6773.257310648137865 2.2137026339358314 3.1487449676409014 1.8124928839213121 1.0507453703670533 9.7130212783813481.1312025785446167 -0.9191203713417053-0.0349808745086193132.19865945514827 32.19865976920392 -3.1405565121668e-07 5.04319e-07 0.0 21.54974334945573 -1.7587264776229858 20.905123 -0.004074007588229828 -2.8542504 0.016663374 22.090984 -0.36655122087409825 -2.2409532 0.4098431 19.830896 0.05489917004070952 -3.5862386 0.4098431 -8.142774803767256 170.2343085995188 0.909421 -0.8930411338806152 2.626863425917633e+22 1282.17521.0054736 -0.6275238394737244 1.034234 3.76935254840918 3.7578791566953673 -0.2190059146437332 3.9768850713391006 2.7561986 -1.7677753073125324 1.7920358297247334 0.177648139364831341.6143876903599021 2.7561986 -0.22640863060951233-0.433468759059906 -0.0626202598214149599897 1 0 0.512121779005605 0.00489836710700053850.5072234118986044 1.8155484 -199.1630888478725334.27863511741595 -28.91581930639344572.14905734016035 72.14905771278723-3.7262687352547246e-075.04319e-07 204.1496641954096 0 -0.43520781397819521.035315803776429 0.16094648838043213 4533.121 -51.247183737712234-51.247183737712234nan nan 3.9261128581479183 -277.590976406493275.696494128336444
6780.2641105417770973 0.245390353058514750.297365905861833260.159682139985030560.085196948960353978.8293781280517580.945320188999176 -1.3487151861190796-0.3256357312202453635.85495989711826 35.8549602904271 -3.933088386264988e-07 2.846336e-07 0.0 21.455110431946277-1.5585298538208008 20.364681 0.0008218263761016395 -1.3366417 0.012262382 20.276108 -0.21644250444128993 -1.2599249 0.26641363 20.328588 0.18773063435298895 -1.4739338 0.26641363 -13.655983988490425159.63295240619402 2.855254 -0.069797731935977942.1299237240088491e+21256.322662.0916967 -0.6132096648216248 2.002991 0.5167124356330226 0.27332072910096816 -0.226788478189218740.5001092072901869 1.6750921 0.17761734906933685 0.4731125055121081 0.226260593005451820.246851912506656281.6750921 -0.40418189764022827-0.7111912369728088 -0.6041487455368042 99991 1 0 -0.284595303589919860.005116558142176688 -0.289711861732096571.0246809 -178.0455509686296 66.09801073444817 -46.14267931512884558.62774848228099 58.62774870405757-2.2177657977211457e-072.846336e-07 195.443884778094140 -0.51304239034652712.091994004496485 -0.105139151215553287735.15 -73.44707433898377 -73.44707433898377 nan nan -31.489219320865587-304.1184333153507 -3.04114398973808
" ], "text/plain": [ "# A_0 A_gaia__dr2_g A_gaia__dr2_gbp A_gaia__dr2_grp E(B-V) age alpha calcium carbon dec dec_Clean dec_Err dec_Sig dform dmod feh gaia__dr2_g gaia__dr2_g_Err gaia__dr2_g_Intrinsic gaia__dr2_g_Sig gaia__dr2_gbp gaia__dr2_gbp_Err gaia__dr2_gbp_Intrinsic gaia__dr2_gbp_Sig gaia__dr2_grp gaia__dr2_grp_Err gaia__dr2_grp_Intrinsic gaia__dr2_grp_Sig glat glon grav helium log10_NH lum mact magnesium mtip mub mudec mudec_Clean mudec_Err mudec_Sig mul mura mura_Clean mura_Err mura_Sig neon nitrogen oxygen parentid partid partitionid pi pi_Clean pi_Err pi_Sig px py pz ra ra_Clean ra_Err ra_Sig rad satid silicon smass sulphur teff vr vr_Clean vr_Err vr_Sig vx vy vz\n", "0 0.35884632402023775 0.26248185076945346 0.3580179075913724 0.20466172853140535 0.11575687871620573 9.692870140075684 0.3099161982536316 -1.1548902988433838 -0.4665657579898834 18.868015796305663 18.868015778565187 1.7740477296505672e-08 1.691274e-07 0.0 21.118883008502078 -0.48292723298072815 19.868202 0.012658170630218355 -1.5258217 0.009227399 20.504751 -0.16132866071418359 -0.8108196 0.17850117 18.7721 -0.22652792657603063 -2.3249161 0.17850117 -14.991976836287705 -177.95690294976794 1.2548169 -0.27809077501296997 2.893921967905143e+21 448.28244 1.1287807 -0.17301106452941895 1.1379396 -0.15521054578348675 -1.3589967238651148 -0.3270482743773745 -1.0319484494877402 1.0656426 1.570216570810097 0.8017470282407725 0.1632684267596954 0.638478601481077 1.0656426 -0.7034189105033875 -0.09953558444976807 -0.4723181426525116 49 1 0 -0.6145466549364363 0.005973424767124436 -0.6205200797035608 0.60885864 -161.60712291393943 -5.765154922998849 -43.30577433967765 74.02626514013622 74.02626502725049 1.1288572264651892e-07 1.691274e-07 167.40815176975818 0 -0.6018809676170349 1.1408413190663447 -0.7073405385017395 4129.926 69.84271348583249 69.84271348583249 nan nan -44.649975087957735 -287.75369665394226 -65.06155285732873\n", "1 0.5999316419840481 0.3951491437418697 0.5843014320686075 0.327567946448567 0.1935263361238865 9.988179206848145 0.11780641227960587 0.09089594334363937 -0.26741287112236023 28.323927000646222 28.323927151103447 -1.5045722692676327e-07 7.903525e-08 0.0 21.313550824736375 -0.6772944927215576 19.078148 0.007681412416285324 -2.6382327 0.005971381 20.259933 -0.07032803519170927 -1.5675896 0.09647448 17.950785 -0.07180883013983862 -3.6185248 0.09647448 -18.536442761608946 165.70363972881054 0.3074857 -0.21355026960372925 4.838158403097162e+21 1852.9218 0.85137147 -0.55948805809021 0.91256934 0.3131789409812662 -0.1038156537874877 -0.14986186911488492 0.04604621532739723 0.55013824 0.48286076265325345 0.5660899886174509 0.17285903010913134 0.39323095850831946 0.55013824 -0.3997773826122284 -0.7002913355827332 -0.08155464380979538 144 1 0 -0.15274891158827747 0.005461222050076329 -0.1582101336383538 0.2845269 -168.23334658322804 42.8707615557611 -58.21183715231835 59.864629221697925 59.86462915929865 6.239927571995072e-08 7.903525e-08 183.1091998879671 0 -0.3103104829788208 0.914106032084624 -0.38073936104774475 3662.81 -25.762507382985724 -25.762507382985724 nan nan -27.923177997805908 -197.62939339603813 16.189820511232778\n", "2 0.17608195267994248 0.13770339427381217 0.1808537019520938 0.10255663320845004 0.05680062989675563 9.579007148742676 -0.6709350943565369 -0.9936923980712891 -0.9404328465461731 33.049066348921734 33.04906634871991 2.0182171922849153e-10 3.6569543e-07 0.0 21.38615087842112 -0.9346992373466492 20.56803 -0.03600744374280186 -0.9198183 0.014033196 21.098215 -0.07126815067421044 -0.3975206 0.32202673 19.32279 -0.5803682522135298 -1.585548 0.32202673 -8.484585766867529 168.82659578503788 1.813334 -0.15246030688285828 1.4200157474188907e+21 204.61284 1.146 -1.605634331703186 1.1506394 -0.3324971332482145 -1.1621686215766727 -0.12728573316172523 -1.0348828884149475 2.0831733 1.242114141035326 0.550241745054646 0.1341992658590463 0.4160424791955997 2.0831733 -0.9467113018035889 -0.7777135372161865 -0.36034896969795227 355 1 0 -2.8237769876652146 0.005281652504745023 -2.8290586401699596 1.3165035 -183.71300451489626 36.287543243958545 -27.93507061765807 70.77774203758443 70.77774167688843 3.6069600030252607e-07 3.6569543e-07 189.33468248840728 0 -0.7150757312774658 1.1526367829812323 -0.12857705354690552 4663.5015 -75.65222985360518 -75.65222985360518 nan nan 38.83221747966123 -175.9006125758071 28.874331454692708\n", "3 0.17861174988368395 0.12892033349928395 0.17747823031542992 0.10134747520582872 0.05761669351086579 9.481425285339355 -1.1020725965499878 -1.0413920879364014 -0.5413622856140137 28.239365554203456 28.23936538065812 1.735453331636329e-07 2.606068e-07 0.0 21.781771243448272 -0.2500181496143341 20.272312 -0.006433955642751986 -1.6319455 0.0116908895 21.282545 0.17120013764702185 -0.8479041 0.24911553 19.201895 -0.20884711491488975 -2.4723763 0.24911553 -4.592810387657218 177.86633782029145 1.2065276 -0.8023523092269897 1.4404173377716447e+21 537.386 1.3710679 -1.3520907163619995 1.369177 0.8879077655799169 0.7300725836482986 -0.21737501118535793 0.9474475948336566 1.551452 -0.27902170477556787 0.5772584642238547 0.15078156191789754 0.42647690230595725 1.551452 -0.6431713700294495 -0.4919672906398773 -0.6521893739700317 359 1 0 0.4967736437882954 0.004401956548212814 0.4923716872400826 0.93818444 -226.28525349857492 8.430630112877303 -18.190511827582277 80.6929868604775 80.69298721775455 -3.572770414100613e-07 2.606068e-07 227.17171081709068 0 -0.7507328391075134 1.3773316962276039 -0.6594099998474121 4003.228 27.33035105089606 27.33035105089606 nan nan -38.03903369359221 -283.6631352491681 0.4139964206713862\n", "4 0.945820473694447 0.6678443096074201 0.930577483029963 0.5330434047235866 0.3051033786111119 9.603869438171387 0.19408565759658813 -0.6169882416725159 -0.7964903116226196 32.353069660512055 32.35306979227988 -1.3176782507459781e-07 1.0227677e-07 0.0 21.17193983212506 -0.869421124458313 19.34399 -0.0006781485567268116 -2.495115 0.0069470913 20.299686 -0.07008269042370177 -1.7327482 0.11953725 18.383387 -3.715706189872538e-05 -3.3215604 0.11953725 -6.795830400723491 171.12720912421034 0.82894164 -1.036791443824768 7.627584465277797e+21 1142.3359 1.1365215 -0.6753354668617249 1.1560693 -0.1418502376652223 -0.09612761524104929 -0.39592602695780055 0.29979841171675126 0.6884891 0.010582078436602305 -0.10484727892409634 0.23107472902875997 -0.3359220079528563 0.6884891 -0.7131249904632568 -1.3325629234313965 -0.8150839805603027 490 1 0 0.33585443984875474 0.005829241304978611 0.33002519854377615 0.36819637 -168.30516566392845 26.27400588210036 -20.299675342958572 74.04671686055353 74.04671679787201 6.268152637233299e-08 1.0227677e-07 171.54891137306748 0 -0.327743262052536 1.159396049264403 -0.09074452519416809 4076.2021 -19.64436147888908 -19.64436147888908 nan nan -31.620475053741032 -368.72625625434716 -49.06766686013123\n", "... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...\n", "674 0.15101800345893615 0.11829324025722467 0.15522325571096357 0.08800020566498619 0.04871548498675359 9.534515380859375 1.2677661180496216 -0.4057428538799286 -0.59430992603302 21.343534159941324 21.34353477135828 -6.11416954624569e-07 3.1462778e-07 0.0 21.592642702411716 -1.1074150800704956 20.454937 -0.005222778643958658 -1.2507769 0.012943381 21.053337 0.03625308061370761 -0.73078275 0.2874478 19.719486 -0.045624436526205674 -1.9155316 0.2874478 -17.503740519827847 176.7477689805897 1.6993186 -0.45311951637268066 1.2178871246688396e+21 276.45468 1.174 0.16035111248493195 1.1788915 -0.20376656041788468 -0.11886161920922561 -0.2466959595835509 0.12783434037432528 1.8276072 -0.014066707641402212 -0.16610418082446088 0.1330427667679313 -0.2991469475923922 1.8276072 -0.48284637928009033 -0.3429383337497711 -0.361402690410614 99192 1 0 1.1484148018479465 0.004802545183696615 1.14361225666425 1.13266 -198.2618205019605 11.265869808201476 -62.626805117261114 68.59347486776855 68.59347435969208 5.080764713807229e-07 3.1462778e-07 208.22292383520687 0 -1.429317593574524 1.1812293030115915 -0.13839608430862427 4680.429 35.33961960826008 35.33961960826008 nan nan -32.91202385958913 -269.96195551855874 -61.869312861074555\n", "675 0.9635104765596588 0.7031919691501173 0.9598374375111064 0.5493862070334198 0.31080983114827704 9.301492691040039 0.18354283273220062 -0.7873296141624451 -0.5832746028900146 36.293908468577015 36.29390899617845 -5.276014372685639e-07 5.21441e-07 0.0 21.986108149698197 -0.6020106673240662 20.938156 -0.00258938495380322 -1.7485539 0.01696236 21.846714 0.0033400161006579407 -1.1025709 0.42019537 20.104744 0.07144425513196052 -2.5021954 0.42019537 -18.73755156006982 151.75488224208473 1.3968682 -0.5497167110443115 7.770245778706926e+21 504.3772 1.508 -0.41846781969070435 1.501947 3.299747393217878 -0.28967745228217223 -0.11117319854321678 -0.17850425373895545 2.8375878 5.791327156440196 6.6591208177663495 0.07693864569556044 6.582182172070789 2.8375878 -0.7871666550636292 -0.6889442801475525 -0.6478202939033508 99272 1 0 -0.6687868125300893 0.004006621987939852 -0.6727934345180292 1.8771876 -208.2157182786129 111.85547750288215 -80.17569728370337 47.36247265711304 47.36247304202983 -3.849167924330201e-07 5.21441e-07 249.58680978890794 0 -0.3253461718559265 1.510235532950821 -0.967689037322998 4293.1353 -70.86896352290194 -70.86896352290194 nan nan 7.312308989325453 -170.89071803084343 -36.78972701008118\n", "676 0.37500312118862084 0.28338113776183976 0.3791604979470969 0.21619101278043099 0.12096874877052285 10.16224193572998 -0.5301219820976257 -0.6149452328681946 -0.9359602928161621 20.238646758946224 20.23864607329147 6.856547555591803e-07 4.3608154e-07 0.0 21.28507619772096 -0.466217964887619 20.76201 -0.009498002225838745 -0.7969467 0.015419783 20.97355 -0.4993928131066979 -0.19129108 0.3675803 19.7897 -0.1881952233286074 -1.5233725 0.3675803 -25.702886990754198 170.24573776609319 1.5625119 -0.6190789937973022 3.024218719263071e+21 199.88933 0.792 -0.9963399171829224 0.87786406 0.42368958053650185 0.9520933879339335 -0.20646777270186295 1.1585611606357964 2.4282553 -0.9299044207050995 -0.3711512804161766 0.11890487188601179 -0.4900561523021884 2.4282553 -0.6988312005996704 -0.2878437340259552 -0.5345746278762817 99878 0 0 -1.7825581982250396 0.00553330692308911 -1.7880915051481288 1.5698935 -160.487938474724 27.589162658182417 -78.38070567241937 57.76421066762154 57.76421003056191 6.370596314964964e-07 4.3608154e-07 180.7237541491236 0 -1.3189185857772827 0.8813156093858641 -0.7903018593788147 4401.795 -4.640823433510391 -4.640823433510391 nan nan -5.309618854522705 -198.32230911254882 -48.23520782470703\n", "677 3.257310648137865 2.2137026339358314 3.1487449676409014 1.8124928839213121 1.0507453703670533 9.713021278381348 1.1312025785446167 -0.9191203713417053 -0.03498087450861931 32.19865945514827 32.19865976920392 -3.1405565121668e-07 5.04319e-07 0.0 21.54974334945573 -1.7587264776229858 20.905123 -0.004074007588229828 -2.8542504 0.016663374 22.090984 -0.36655122087409825 -2.2409532 0.4098431 19.830896 0.05489917004070952 -3.5862386 0.4098431 -8.142774803767256 170.2343085995188 0.909421 -0.8930411338806152 2.626863425917633e+22 1282.1752 1.0054736 -0.6275238394737244 1.034234 3.76935254840918 3.7578791566953673 -0.2190059146437332 3.9768850713391006 2.7561986 -1.7677753073125324 1.7920358297247334 0.17764813936483134 1.6143876903599021 2.7561986 -0.22640863060951233 -0.433468759059906 -0.06262025982141495 99897 1 0 0.512121779005605 0.0048983671070005385 0.5072234118986044 1.8155484 -199.16308884787253 34.27863511741595 -28.915819306393445 72.14905734016035 72.14905771278723 -3.7262687352547246e-07 5.04319e-07 204.1496641954096 0 -0.4352078139781952 1.035315803776429 0.16094648838043213 4533.121 -51.247183737712234 -51.247183737712234 nan nan 3.9261128581479183 -277.59097640649327 5.696494128336444\n", "678 0.2641105417770973 0.24539035305851475 0.29736590586183326 0.15968213998503056 0.08519694896035397 8.829378128051758 0.945320188999176 -1.3487151861190796 -0.32563573122024536 35.85495989711826 35.8549602904271 -3.933088386264988e-07 2.846336e-07 0.0 21.455110431946277 -1.5585298538208008 20.364681 0.0008218263761016395 -1.3366417 0.012262382 20.276108 -0.21644250444128993 -1.2599249 0.26641363 20.328588 0.18773063435298895 -1.4739338 0.26641363 -13.655983988490425 159.63295240619402 2.855254 -0.06979773193597794 2.1299237240088491e+21 256.32266 2.0916967 -0.6132096648216248 2.002991 0.5167124356330226 0.27332072910096816 -0.22678847818921874 0.5001092072901869 1.6750921 0.17761734906933685 0.4731125055121081 0.22626059300545182 0.24685191250665628 1.6750921 -0.40418189764022827 -0.7111912369728088 -0.6041487455368042 99991 1 0 -0.28459530358991986 0.005116558142176688 -0.28971186173209657 1.0246809 -178.0455509686296 66.09801073444817 -46.142679315128845 58.62774848228099 58.62774870405757 -2.2177657977211457e-07 2.846336e-07 195.44388477809414 0 -0.5130423903465271 2.091994004496485 -0.10513915121555328 7735.15 -73.44707433898377 -73.44707433898377 nan nan -31.489219320865587 -304.1184333153507 -3.04114398973808" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "survey[~survey.gaia__dr2_g.isna()]" ] } ], "metadata": { "kernelspec": { "display_name": "galaxia_kernels", "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.12.13" } }, "nbformat": 4, "nbformat_minor": 5 }