{ "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-16\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.5.0+9.g83a85f8'" ] }, "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', 'massinit', 'age', 'feh', 'helium', 'carbon', 'nitrogen', 'oxygen', 'neon', 'magnesium', 'silicon', 'sulphur', 'calcium', 'alpha', 'parentid', 'partitionid', 'dform', 'id', 'log10_NH', 'mass'])" ] }, "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 `massinit`: particle initial 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 phase space kernel sizes, but one can include pre-computed kernel sizes with the following entry:\n", "- key `kernels`: phase-space kernel radii in $kpc$ and $km/s$ (Nx2)\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, `k_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 kernels 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_kernels_docs()` for `k_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/main/lib/python3.10/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/main/lib/python3.10/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_gens` – 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", "- `caching` - Flag to activate checkpoint caching of survey production. Default to False\n", "- `append_hash` - Flag to append a unique hash to survey file names. Default to True\n", "- `no_post_processing` - Flag to skip post-processing entirely and only produce `Galaxia`'s outputs. Default to False\n", "- `partitioning_rule` - Function that returns an array of partition id from a given dataframe of `Galaxia`'s catalog properties. Default to None\n", "- `max_pp_workers` - In post-processing, parallelization is done per survey partition, for which this integer parameter control the max number of independent workers. Default to None\n", "- `pp_auto_flush` - Flag to control auto-flushing of survey data during post-processing steps. Default to True (it is strongly recommended to keep it True)\n", "\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 GADGET format file: to_enbid \n", "Allocated 3.8147 MByte for particle storage.\n", "....done with reading.\n", "Type = 1 Particles = 100000\n", "Total particles = 100000\n", "WARNING! particles with multiple mass results might not be accurate\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.162309 s \n", "\n", "Density Calculation. Smoothing .....\n", "Evaluated = 1 % Time Left = 1.273820 s of 1.286700 s Par no = 94545 Density = 2.073172e+02 \n", "Evaluated = 2 % Time Left = 1.254485 s of 1.280100 s Par no = 7661 Density = 9.490810e+02 \n", "Evaluated = 3 % Time Left = 1.210483 s of 1.247933 s Par no = 45728 Density = 4.426955e+02 \n", "Evaluated = 4 % Time Left = 1.213859 s of 1.264450 s Par no = 50673 Density = 7.158524e+02 \n", "Evaluated = 5 % Time Left = 1.191154 s of 1.253860 s Par no = 84687 Density = 1.401414e+02 \n", "Evaluated = 6 % Time Left = 1.177604 s of 1.252783 s Par no = 55217 Density = 6.157888e+02 \n", "Evaluated = 7 % Time Left = 1.179493 s of 1.268286 s Par no = 76523 Density = 4.584541e+02 \n", "Evaluated = 8 % Time Left = 1.175414 s of 1.277637 s Par no = 17728 Density = 3.759916e+02 \n", "Evaluated = 9 % Time Left = 1.172017 s of 1.287944 s Par no = 24571 Density = 3.816299e+02 \n", "Evaluated = 10 % Time Left = 1.164506 s of 1.293910 s Par no = 56796 Density = 2.710653e+02 \n", "Evaluated = 11 % Time Left = 1.134972 s of 1.275264 s Par no = 62092 Density = 2.468291e+02 \n", "Evaluated = 12 % Time Left = 1.106697 s of 1.257625 s Par no = 92210 Density = 3.054839e+02 \n", "Evaluated = 13 % Time Left = 1.094662 s of 1.258246 s Par no = 52517 Density = 4.655019e+02 \n", "Evaluated = 14 % Time Left = 1.082070 s of 1.258236 s Par no = 53473 Density = 6.649532e+02 \n", "Evaluated = 15 % Time Left = 1.068001 s of 1.256487 s Par no = 36690 Density = 3.385491e+02 \n", "Evaluated = 16 % Time Left = 1.052145 s of 1.252569 s Par no = 52644 Density = 2.297310e+02 \n", "Evaluated = 17 % Time Left = 1.038625 s of 1.251371 s Par no = 39271 Density = 3.912303e+02 \n", "Evaluated = 18 % Time Left = 1.034923 s of 1.262117 s Par no = 72159 Density = 5.531728e+02 \n", "Evaluated = 19 % Time Left = 1.025307 s of 1.265826 s Par no = 70783 Density = 3.475029e+02 \n", "Evaluated = 20 % Time Left = 1.013295 s of 1.266635 s Par no = 4946 Density = 3.470684e+02 \n", "Evaluated = 21 % Time Left = 1.001576 s of 1.267833 s Par no = 55704 Density = 4.690976e+02 \n", "Evaluated = 22 % Time Left = 0.987762 s of 1.266377 s Par no = 65148 Density = 5.908168e+02 \n", "Evaluated = 23 % Time Left = 0.971373 s of 1.261539 s Par no = 66239 Density = 7.375975e+01 \n", "Evaluated = 24 % Time Left = 0.959820 s of 1.262937 s Par no = 5808 Density = 2.837905e+02 \n", "Evaluated = 25 % Time Left = 0.948050 s of 1.264084 s Par no = 82119 Density = 3.643292e+02 \n", "Evaluated = 26 % Time Left = 0.936733 s of 1.265873 s Par no = 18234 Density = 5.550048e+02 \n", "Evaluated = 27 % Time Left = 0.922580 s of 1.263826 s Par no = 17562 Density = 3.339526e+02 \n", "Evaluated = 28 % Time Left = 0.915550 s of 1.271614 s Par no = 82032 Density = 5.470850e+02 \n", "Evaluated = 29 % Time Left = 0.901753 s of 1.270093 s Par no = 49435 Density = 4.979107e+02 \n", "Evaluated = 30 % Time Left = 0.887849 s of 1.268373 s Par no = 55540 Density = 2.183091e+02 \n", "Evaluated = 31 % Time Left = 0.877106 s of 1.271187 s Par no = 89348 Density = 6.990241e+02 \n", "Evaluated = 32 % Time Left = 0.863549 s of 1.269944 s Par no = 71900 Density = 5.763120e+02 \n", "Evaluated = 33 % Time Left = 0.851243 s of 1.270530 s Par no = 19060 Density = 3.338016e+02 \n", "Evaluated = 34 % Time Left = 0.837537 s of 1.269015 s Par no = 48097 Density = 5.252333e+02 \n", "Evaluated = 35 % Time Left = 0.822455 s of 1.265334 s Par no = 92512 Density = 3.829829e+02 \n", "Evaluated = 36 % Time Left = 0.808691 s of 1.263600 s Par no = 11034 Density = 3.941921e+01 \n", "Evaluated = 37 % Time Left = 0.794852 s of 1.261689 s Par no = 92858 Density = 3.082108e+02 \n", "Evaluated = 38 % Time Left = 0.783785 s of 1.264189 s Par no = 64470 Density = 9.409036e+02 \n", "Evaluated = 39 % Time Left = 0.770327 s of 1.262851 s Par no = 93035 Density = 5.775200e+01 \n", "Evaluated = 40 % Time Left = 0.757684 s of 1.262827 s Par no = 66380 Density = 4.528047e+02 \n", "Evaluated = 41 % Time Left = 0.746865 s of 1.265895 s Par no = 47276 Density = 6.639922e+02 \n", "Evaluated = 42 % Time Left = 0.733574 s of 1.264805 s Par no = 98619 Density = 3.078024e+02 \n", "Evaluated = 43 % Time Left = 0.722016 s of 1.266716 s Par no = 83819 Density = 5.742248e+02 \n", "Evaluated = 44 % Time Left = 0.711135 s of 1.269907 s Par no = 40506 Density = 2.339407e+02 \n", "Evaluated = 45 % Time Left = 0.695202 s of 1.264027 s Par no = 76662 Density = 1.204026e+02 \n", "Evaluated = 46 % Time Left = 0.684348 s of 1.267335 s Par no = 64120 Density = 5.082795e+02 \n", "Evaluated = 47 % Time Left = 0.671974 s of 1.267900 s Par no = 3645 Density = 6.558582e+02 \n", "Evaluated = 48 % Time Left = 0.658162 s of 1.265721 s Par no = 7894 Density = 8.308378e+02 \n", "Evaluated = 49 % Time Left = 0.646334 s of 1.267347 s Par no = 4532 Density = 9.386861e+01 \n", "Evaluated = 50 % Time Left = 0.632788 s of 1.265602 s Par no = 22331 Density = 2.882854e+02 \n", "Evaluated = 51 % Time Left = 0.620248 s of 1.265837 s Par no = 35966 Density = 1.752636e+02 \n", "Evaluated = 52 % Time Left = 0.607555 s of 1.265765 s Par no = 57105 Density = 3.234318e+02 \n", "Evaluated = 53 % Time Left = 0.593450 s of 1.262687 s Par no = 89855 Density = 1.073572e+02 \n", "Evaluated = 54 % Time Left = 0.579516 s of 1.259844 s Par no = 8210 Density = 4.065989e+01 \n", "Evaluated = 55 % Time Left = 0.566915 s of 1.259840 s Par no = 61142 Density = 3.700869e+02 \n", "Evaluated = 56 % Time Left = 0.555144 s of 1.261720 s Par no = 77421 Density = 7.643566e+02 \n", "Evaluated = 57 % Time Left = 0.541872 s of 1.260196 s Par no = 94858 Density = 9.271213e+01 \n", "Evaluated = 58 % Time Left = 0.527967 s of 1.257095 s Par no = 73283 Density = 2.528433e+02 \n", "Evaluated = 59 % Time Left = 0.514749 s of 1.255517 s Par no = 41030 Density = 1.077854e+02 \n", "Evaluated = 60 % Time Left = 0.502717 s of 1.256823 s Par no = 1288 Density = 1.555944e+02 \n", "Evaluated = 61 % Time Left = 0.489859 s of 1.256080 s Par no = 60528 Density = 1.972560e+02 \n", "Evaluated = 62 % Time Left = 0.476699 s of 1.254503 s Par no = 35369 Density = 2.816853e+02 \n", "Evaluated = 63 % Time Left = 0.464544 s of 1.255557 s Par no = 81571 Density = 1.742148e+02 \n", "Evaluated = 64 % Time Left = 0.451223 s of 1.253431 s Par no = 74571 Density = 2.872113e+02 \n", "Evaluated = 65 % Time Left = 0.438006 s of 1.251482 s Par no = 66808 Density = 2.310972e+02 \n", "Evaluated = 66 % Time Left = 0.424732 s of 1.249248 s Par no = 66795 Density = 9.541369e+01 \n", "Evaluated = 67 % Time Left = 0.413082 s of 1.251803 s Par no = 27749 Density = 4.961322e+02 \n", "Evaluated = 68 % Time Left = 0.400098 s of 1.250346 s Par no = 76853 Density = 3.065290e+01 \n", "Evaluated = 69 % Time Left = 0.387386 s of 1.249671 s Par no = 27545 Density = 4.113157e+02 \n", "Evaluated = 70 % Time Left = 0.376180 s of 1.253976 s Par no = 94330 Density = 9.211449e+02 \n", "Evaluated = 71 % Time Left = 0.362905 s of 1.251441 s Par no = 60575 Density = 6.903663e+01 \n", "Evaluated = 72 % Time Left = 0.350146 s of 1.250567 s Par no = 69730 Density = 1.650653e+01 \n", "Evaluated = 73 % Time Left = 0.337283 s of 1.249241 s Par no = 25092 Density = 1.591412e+00 \n", "Evaluated = 74 % Time Left = 0.325312 s of 1.251247 s Par no = 92614 Density = 1.146144e+02 \n", "Evaluated = 75 % Time Left = 0.313020 s of 1.252131 s Par no = 78117 Density = 3.205864e+02 \n", "Evaluated = 76 % Time Left = 0.300023 s of 1.250150 s Par no = 34177 Density = 1.334262e+02 \n", "Evaluated = 77 % Time Left = 0.287393 s of 1.249590 s Par no = 37747 Density = 2.525646e+01 \n", "Evaluated = 78 % Time Left = 0.274951 s of 1.249833 s Par no = 71398 Density = 2.129484e+02 \n", "Evaluated = 79 % Time Left = 0.262358 s of 1.249385 s Par no = 85573 Density = 3.863819e+01 \n", "Evaluated = 80 % Time Left = 0.249632 s of 1.248224 s Par no = 2558 Density = 3.960752e+01 \n", "Evaluated = 81 % Time Left = 0.236996 s of 1.247412 s Par no = 49255 Density = 2.495420e+01 \n", "Evaluated = 82 % Time Left = 0.224137 s of 1.245277 s Par no = 64542 Density = 3.184729e+02 \n", "Evaluated = 83 % Time Left = 0.211627 s of 1.244939 s Par no = 32410 Density = 6.581928e+01 \n", "Evaluated = 84 % Time Left = 0.199318 s of 1.245815 s Par no = 72779 Density = 5.631017e+02 \n", "Evaluated = 85 % Time Left = 0.186454 s of 1.243112 s Par no = 31876 Density = 1.751137e+02 \n", "Evaluated = 86 % Time Left = 0.173909 s of 1.242298 s Par no = 40987 Density = 2.338241e+02 \n", "Evaluated = 87 % Time Left = 0.161622 s of 1.243345 s Par no = 51332 Density = 4.748823e+02 \n", "Evaluated = 88 % Time Left = 0.149264 s of 1.243970 s Par no = 28773 Density = 2.243994e+02 \n", "Evaluated = 89 % Time Left = 0.136789 s of 1.243649 s Par no = 78440 Density = 4.169291e+01 \n", "Evaluated = 90 % Time Left = 0.124370 s of 1.243820 s Par no = 83393 Density = 6.845449e+02 \n", "Evaluated = 91 % Time Left = 0.111948 s of 1.244003 s Par no = 45461 Density = 2.227390e+01 \n", "Evaluated = 92 % Time Left = 0.099556 s of 1.244602 s Par no = 59958 Density = 1.922444e+01 \n", "Evaluated = 93 % Time Left = 0.087227 s of 1.246272 s Par no = 71816 Density = 8.541215e+00 \n", "Evaluated = 94 % Time Left = 0.074864 s of 1.247948 s Par no = 28276 Density = 8.799874e+02 \n", "Evaluated = 95 % Time Left = 0.062546 s of 1.251176 s Par no = 80170 Density = 1.667844e+02 \n", "Evaluated = 96 % Time Left = 0.050143 s of 1.253883 s Par no = 45447 Density = 6.512313e+02 \n", "Evaluated = 97 % Time Left = 0.037625 s of 1.254595 s Par no = 39654 Density = 2.771305e+02 \n", "Evaluated = 98 % Time Left = 0.025100 s of 1.255604 s Par no = 73295 Density = 4.006583e+00 \n", "Evaluated = 99 % Time Left = 0.012554 s of 1.256668 s Par no = 43505 Density = 9.523295e+01 \n", "Evaluated = 100 % Time Left = -0.000013 s of 1.258743 s Par no = 96564 Density = 3.785476e+00 \n", "\n", "Total Smoothing Time = 1.2588 s\n", "1 3\n", "\n", "writing snapshot file ....\n", "to_enbid_d3n64.est\n", "done with snapshot.\n", "\n", "Total Time = 1.44277 s \n", "\n", "Dimensions = 3\n", "Normalization constant of Kernel type 3: 0.596831\n", "Reading GADGET format file: to_enbid \n", "Allocated 3.8147 MByte for particle storage.\n", "....done with reading.\n", "Type = 1 Particles = 100000\n", "Total particles = 100000\n", "WARNING! particles with multiple mass results might not be accurate\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.163697 s \n", "\n", "Density Calculation. Smoothing .....\n", "Evaluated = 1 % Time Left = 1.217886 s of 1.230200 s Par no = 27998 Density = 2.171451e+02 \n", "Evaluated = 2 % Time Left = 1.184465 s of 1.208650 s Par no = 55969 Density = 8.945573e+01 \n", "Evaluated = 3 % Time Left = 1.211259 s of 1.248733 s Par no = 23539 Density = 1.659165e+02 \n", "Evaluated = 4 % Time Left = 1.243283 s of 1.295100 s Par no = 42127 Density = 6.482698e+01 \n", "Evaluated = 5 % Time Left = 1.235595 s of 1.300640 s Par no = 90717 Density = 3.819016e+01 \n", "Evaluated = 6 % Time Left = 1.191704 s of 1.267783 s Par no = 32585 Density = 4.963705e+01 \n", "Evaluated = 7 % Time Left = 1.190321 s of 1.279929 s Par no = 63192 Density = 8.197108e+01 \n", "Evaluated = 8 % Time Left = 1.177093 s of 1.279463 s Par no = 54953 Density = 1.460569e+01 \n", "Evaluated = 9 % Time Left = 1.155242 s of 1.269511 s Par no = 29569 Density = 5.293639e+01 \n", "Evaluated = 10 % Time Left = 1.154921 s of 1.283260 s Par no = 36769 Density = 2.903416e+01 \n", "Evaluated = 11 % Time Left = 1.141080 s of 1.282127 s Par no = 48903 Density = 3.880676e+01 \n", "Evaluated = 12 % Time Left = 1.125456 s of 1.278942 s Par no = 87768 Density = 2.989370e+01 \n", "Evaluated = 13 % Time Left = 1.107002 s of 1.272431 s Par no = 68373 Density = 1.613992e+02 \n", "Evaluated = 14 % Time Left = 1.097249 s of 1.275886 s Par no = 35219 Density = 5.155452e+01 \n", "Evaluated = 15 % Time Left = 1.077204 s of 1.267313 s Par no = 59873 Density = 1.214169e+02 \n", "Evaluated = 16 % Time Left = 1.065554 s of 1.268531 s Par no = 68594 Density = 3.922213e+01 \n", "Evaluated = 17 % Time Left = 1.048844 s of 1.263682 s Par no = 70126 Density = 1.702955e+01 \n", "Evaluated = 18 % Time Left = 1.043255 s of 1.272278 s Par no = 87334 Density = 8.252219e+01 \n", "Evaluated = 19 % Time Left = 1.029493 s of 1.270995 s Par no = 40178 Density = 6.439714e+01 \n", "Evaluated = 20 % Time Left = 1.016243 s of 1.270320 s Par no = 5913 Density = 1.354404e+01 \n", "Evaluated = 21 % Time Left = 0.999755 s of 1.265529 s Par no = 33195 Density = 2.604624e+00 \n", "Evaluated = 22 % Time Left = 0.987272 s of 1.265750 s Par no = 69263 Density = 1.359985e+02 \n", "Evaluated = 23 % Time Left = 0.972280 s of 1.262717 s Par no = 7749 Density = 6.080642e+01 \n", "Evaluated = 24 % Time Left = 0.953217 s of 1.254250 s Par no = 71558 Density = 1.210650e+02 \n", "Evaluated = 25 % Time Left = 0.937835 s of 1.250464 s Par no = 45848 Density = 1.070795e+02 \n", "Evaluated = 26 % Time Left = 0.926883 s of 1.252562 s Par no = 36853 Density = 3.120538e+01 \n", "Evaluated = 27 % Time Left = 0.916849 s of 1.255974 s Par no = 80804 Density = 1.585974e+02 \n", "Evaluated = 28 % Time Left = 0.910782 s of 1.264993 s Par no = 89314 Density = 4.844725e+01 \n", "Evaluated = 29 % Time Left = 0.899829 s of 1.267383 s Par no = 24521 Density = 1.661036e+02 \n", "Evaluated = 30 % Time Left = 0.887534 s of 1.267923 s Par no = 94700 Density = 1.208696e+02 \n", "Evaluated = 31 % Time Left = 0.875936 s of 1.269490 s Par no = 54669 Density = 1.868439e+02 \n", "Evaluated = 32 % Time Left = 0.861792 s of 1.267359 s Par no = 79174 Density = 4.603831e+01 \n", "Evaluated = 33 % Time Left = 0.847405 s of 1.264803 s Par no = 68045 Density = 1.456620e+02 \n", "Evaluated = 34 % Time Left = 0.835415 s of 1.265800 s Par no = 99153 Density = 3.661007e+01 \n", "Evaluated = 35 % Time Left = 0.823417 s of 1.266814 s Par no = 89524 Density = 8.644379e+01 \n", "Evaluated = 36 % Time Left = 0.811068 s of 1.267314 s Par no = 12288 Density = 5.322282e+01 \n", "Evaluated = 37 % Time Left = 0.798744 s of 1.267868 s Par no = 59478 Density = 1.532470e+02 \n", "Evaluated = 38 % Time Left = 0.787255 s of 1.269787 s Par no = 16941 Density = 8.638396e+01 \n", "Evaluated = 39 % Time Left = 0.777012 s of 1.273810 s Par no = 42784 Density = 1.020941e+02 \n", "Evaluated = 40 % Time Left = 0.766589 s of 1.277670 s Par no = 84531 Density = 1.737449e+02 \n", "Evaluated = 41 % Time Left = 0.756595 s of 1.282385 s Par no = 66470 Density = 9.489573e+01 \n", "Evaluated = 42 % Time Left = 0.744873 s of 1.284286 s Par no = 32233 Density = 7.794218e+01 \n", "Evaluated = 43 % Time Left = 0.731106 s of 1.282665 s Par no = 14257 Density = 1.429659e+01 \n", "Evaluated = 44 % Time Left = 0.719546 s of 1.284927 s Par no = 83299 Density = 1.562222e+01 \n", "Evaluated = 45 % Time Left = 0.705614 s of 1.282958 s Par no = 31159 Density = 2.521319e+01 \n", "Evaluated = 46 % Time Left = 0.694300 s of 1.285765 s Par no = 49702 Density = 1.460362e+02 \n", "Evaluated = 47 % Time Left = 0.681405 s of 1.285694 s Par no = 12255 Density = 5.861372e+01 \n", "Evaluated = 48 % Time Left = 0.669702 s of 1.287913 s Par no = 99186 Density = 1.738853e+02 \n", "Evaluated = 49 % Time Left = 0.657144 s of 1.288543 s Par no = 77170 Density = 1.840628e+02 \n", "Evaluated = 50 % Time Left = 0.643249 s of 1.286524 s Par no = 44414 Density = 4.307461e+00 \n", "Evaluated = 51 % Time Left = 0.631599 s of 1.289004 s Par no = 44617 Density = 1.121407e+02 \n", "Evaluated = 52 % Time Left = 0.621462 s of 1.294738 s Par no = 133 Density = 4.638072e+01 \n", "Evaluated = 53 % Time Left = 0.607019 s of 1.291557 s Par no = 71723 Density = 1.046514e+02 \n", "Evaluated = 54 % Time Left = 0.592164 s of 1.287341 s Par no = 39551 Density = 3.995141e+01 \n", "Evaluated = 55 % Time Left = 0.578321 s of 1.285187 s Par no = 93018 Density = 1.550717e+02 \n", "Evaluated = 56 % Time Left = 0.565430 s of 1.285096 s Par no = 86261 Density = 2.966780e+01 \n", "Evaluated = 57 % Time Left = 0.551268 s of 1.282047 s Par no = 35473 Density = 7.820364e+00 \n", "Evaluated = 58 % Time Left = 0.537366 s of 1.279474 s Par no = 5858 Density = 1.709624e+02 \n", "Evaluated = 59 % Time Left = 0.523300 s of 1.276373 s Par no = 72956 Density = 2.428934e+01 \n", "Evaluated = 60 % Time Left = 0.510206 s of 1.275547 s Par no = 35202 Density = 1.415912e+02 \n", "Evaluated = 61 % Time Left = 0.496994 s of 1.274377 s Par no = 36040 Density = 6.637925e+01 \n", "Evaluated = 62 % Time Left = 0.483388 s of 1.272108 s Par no = 16293 Density = 5.401261e+00 \n", "Evaluated = 63 % Time Left = 0.470853 s of 1.272611 s Par no = 80591 Density = 2.856942e+01 \n", "Evaluated = 64 % Time Left = 0.457367 s of 1.270498 s Par no = 31354 Density = 5.230644e+01 \n", "Evaluated = 65 % Time Left = 0.445038 s of 1.271574 s Par no = 60568 Density = 1.536678e+01 \n", "Evaluated = 66 % Time Left = 0.431423 s of 1.268929 s Par no = 73768 Density = 1.115788e+01 \n", "Evaluated = 67 % Time Left = 0.420478 s of 1.274213 s Par no = 59537 Density = 1.131722e+02 \n", "Evaluated = 68 % Time Left = 0.408279 s of 1.275912 s Par no = 88477 Density = 1.674836e+01 \n", "Evaluated = 69 % Time Left = 0.396154 s of 1.277957 s Par no = 54406 Density = 2.433154e+01 \n", "Evaluated = 70 % Time Left = 0.383890 s of 1.279677 s Par no = 14370 Density = 3.662899e+00 \n", "Evaluated = 71 % Time Left = 0.371719 s of 1.281832 s Par no = 31474 Density = 9.792092e-01 \n", "Evaluated = 72 % Time Left = 0.359379 s of 1.283542 s Par no = 46820 Density = 1.158320e+02 \n", "Evaluated = 73 % Time Left = 0.346205 s of 1.282289 s Par no = 62606 Density = 1.650805e+01 \n", "Evaluated = 74 % Time Left = 0.333304 s of 1.281988 s Par no = 54937 Density = 8.059295e+01 \n", "Evaluated = 75 % Time Left = 0.319762 s of 1.279099 s Par no = 86775 Density = 4.669335e+00 \n", "Evaluated = 76 % Time Left = 0.306463 s of 1.276982 s Par no = 82481 Density = 1.206733e+01 \n", "Evaluated = 77 % Time Left = 0.293332 s of 1.275412 s Par no = 52308 Density = 8.028543e+01 \n", "Evaluated = 78 % Time Left = 0.280839 s of 1.276599 s Par no = 4377 Density = 6.577021e+01 \n", "Evaluated = 79 % Time Left = 0.268479 s of 1.278532 s Par no = 40260 Density = 1.177500e+02 \n", "Evaluated = 80 % Time Left = 0.255872 s of 1.279424 s Par no = 53065 Density = 1.968351e+02 \n", "Evaluated = 81 % Time Left = 0.243246 s of 1.280310 s Par no = 26435 Density = 2.804098e+00 \n", "Evaluated = 82 % Time Left = 0.230164 s of 1.278762 s Par no = 86123 Density = 1.680275e+01 \n", "Evaluated = 83 % Time Left = 0.217131 s of 1.277314 s Par no = 41271 Density = 1.177749e+01 \n", "Evaluated = 84 % Time Left = 0.204722 s of 1.279590 s Par no = 12204 Density = 5.935811e+01 \n", "Evaluated = 85 % Time Left = 0.191933 s of 1.279638 s Par no = 67785 Density = 1.943354e+01 \n", "Evaluated = 86 % Time Left = 0.179082 s of 1.279250 s Par no = 23611 Density = 9.835917e+01 \n", "Evaluated = 87 % Time Left = 0.166544 s of 1.281206 s Par no = 34587 Density = 1.865564e+01 \n", "Evaluated = 88 % Time Left = 0.153683 s of 1.280795 s Par no = 34487 Density = 1.615047e+02 \n", "Evaluated = 89 % Time Left = 0.141111 s of 1.282942 s Par no = 98454 Density = 3.704479e+00 \n", "Evaluated = 90 % Time Left = 0.128477 s of 1.284894 s Par no = 13623 Density = 1.174157e+02 \n", "Evaluated = 91 % Time Left = 0.115927 s of 1.288221 s Par no = 16259 Density = 9.483730e+01 \n", "Evaluated = 92 % Time Left = 0.103060 s of 1.288409 s Par no = 71674 Density = 1.076167e+00 \n", "Evaluated = 93 % Time Left = 0.090085 s of 1.287110 s Par no = 4103 Density = 1.704676e+02 \n", "Evaluated = 94 % Time Left = 0.077316 s of 1.288821 s Par no = 62978 Density = 1.411068e+02 \n", "Evaluated = 95 % Time Left = 0.064585 s of 1.291953 s Par no = 79016 Density = 1.412207e+02 \n", "Evaluated = 96 % Time Left = 0.051694 s of 1.292664 s Par no = 69795 Density = 3.356945e+01 \n", "Evaluated = 97 % Time Left = 0.038739 s of 1.291739 s Par no = 63171 Density = 6.152881e+01 \n", "Evaluated = 98 % Time Left = 0.025842 s of 1.292722 s Par no = 84522 Density = 1.332426e+02 \n", "Evaluated = 99 % Time Left = 0.012936 s of 1.294907 s Par no = 84660 Density = 1.301089e+02 \n", "Evaluated = 100 % Time Left = -0.000013 s of 1.297500 s Par no = 37622 Density = 1.020922e-01 \n", "\n", "Total Smoothing Time = 1.29756 s\n", "1 3\n", "\n", "writing snapshot file ....\n", "to_enbid_d3n64.est\n", "done with snapshot.\n", "\n", "Total Time = 1.48319 s \n", "\n", "Executing JOB 1/1 = /home/athob/miniconda3/envs/main/.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/main/.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] -72.25056128347201 \n", "JOB 1/1 | rSun[1] 180.23501013398354 \n", "JOB 1/1 | rSun[2] 47.908219714638804 \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/main/.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 -72.2506 180.235 47.9082 12.9 245.6 7.78\n", "JOB 1/1 | Reading Isochrones from dir- /home/athob/miniconda3/envs/main/.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.69504 \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=7768601 accepted=1236527 rejected=6532074\n", "JOB 1/1 | -----------Done---------------\n", "JOB 1/1 | Total stars written 1236527 \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= 12.6233 \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', 'minit', 'mtip', 'neon', 'nitrogen', 'oxygen', 'parentid', 'partid', 'partitionid', 'px', 'py', 'pz', 'rad', 'satid', 'silicon', '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', 'glon', 'mub', 'vr']\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", "['mura', 'mudec', 'ra', 'dec']\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/16/26 22:48:23] 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/16/26 22:48:23]\u001b[0m\u001b[2;36m \u001b[0m\u001b[33mWARNING \u001b[0m could not close memmap for column \u001b]8;id=8662479;file:///home/athob/miniconda3/envs/main/lib/python3.10/site-packages/vaex/dataset_mmap.py\u001b\\\u001b[2mdataset_mmap.py\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=8662480;file:///home/athob/miniconda3/envs/main/lib/python3.10/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_grp_Intrinsic', 'gaia__dr2_g_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_grp', 'gaia__dr2_g', 'gaia__dr2_gbp']\n" ] }, { "data": { "text/html": [ "
[04/16/26 22:48:24] 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/16/26 22:48:24]\u001b[0m\u001b[2;36m \u001b[0m\u001b[33mWARNING \u001b[0m could not close memmap for column \u001b]8;id=8662485;file:///home/athob/miniconda3/envs/main/lib/python3.10/site-packages/vaex/dataset_mmap.py\u001b\\\u001b[2mdataset_mmap.py\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=8662486;file:///home/athob/miniconda3/envs/main/lib/python3.10/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', 'log10_NH', 'A_gaia__dr2_g', 'E(B-V)', 'A_gaia__dr2_gbp', 'A_0']\n", "Overwritten the following quantities to /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.h5\n", "['gaia__dr2_grp', 'gaia__dr2_g', 'gaia__dr2_gbp']\n" ] }, { "data": { "text/html": [ "
[04/16/26 22:48:46] 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/16/26 22:48:46]\u001b[0m\u001b[2;36m \u001b[0m\u001b[33mWARNING \u001b[0m could not close memmap for column \u001b]8;id=8662491;file:///home/athob/miniconda3/envs/main/lib/python3.10/site-packages/vaex/dataset_mmap.py\u001b\\\u001b[2mdataset_mmap.py\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=8662492;file:///home/athob/miniconda3/envs/main/lib/python3.10/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 error_modeling post-processing pipeline\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/athob/miniconda3/envs/main/lib/python3.10/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/main/lib/python3.10/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/main/lib/python3.10/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", "/home/athob/miniconda3/envs/main/lib/python3.10/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", "['pi_Clean', 'pi_Sig', 'dec_Clean', 'mudec_Clean', 'gaia__dr2_gbp_Sig', 'gaia__dr2_grp_Err', 'mudec_Err', 'gaia__dr2_gbp_Err', 'mudec_Sig', 'dec_Err', 'vr_Sig', 'mura_Sig', 'gaia__dr2_g_Err', 'pi_Err', 'ra_Clean', 'vr_Err', 'vr_Clean', 'gaia__dr2_g_Sig', 'mura_Err', 'dec_Sig', 'ra_Sig', 'ra_Err', 'mura_Clean', 'gaia__dr2_grp_Sig']\n", "Overwritten the following quantities to /home/athob/Software/repos/py-ananke/jupyter/survey.sim.0.h5\n", "['pi', 'gaia__dr2_g', 'mura', 'ra', 'gaia__dr2_grp', 'mudec', 'vr', 'gaia__dr2_gbp', 'dec']\n" ] }, { "data": { "text/html": [ "
[04/16/26 22:48:48] 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/16/26 22:48:48]\u001b[0m\u001b[2;36m \u001b[0m\u001b[33mWARNING \u001b[0m could not close memmap for column \u001b]8;id=8662497;file:///home/athob/miniconda3/envs/main/lib/python3.10/site-packages/vaex/dataset_mmap.py\u001b\\\u001b[2mdataset_mmap.py\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=8662498;file:///home/athob/miniconda3/envs/main/lib/python3.10/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", "['glon', 'mul', 'mub', 'glat']\n" ] }, { "data": { "text/html": [ "
[04/16/26 22:48:50] 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/16/26 22:48:50]\u001b[0m\u001b[2;36m \u001b[0m\u001b[33mWARNING \u001b[0m could not close memmap for column \u001b]8;id=8662503;file:///home/athob/miniconda3/envs/main/lib/python3.10/site-packages/vaex/dataset_mmap.py\u001b\\\u001b[2mdataset_mmap.py\u001b[0m\u001b]8;;\u001b\\\u001b[2m:\u001b[0m\u001b]8;id=8662504;file:///home/athob/miniconda3/envs/main/lib/python3.10/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`, `minit` 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 minit 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 sulphur teff vr vr_Clean vr_Err vr_Sig vx vy vz
0 0.7093977911593942 0.5527585393775436 0.7284652298464588 0.41253265460200594 0.22883799714819167 10.0163259506225590.3380742371082306-1.1830222606658936-0.39430123567581177nan -67.64321024253071 nan nan 0.0 21.61633894249464 -0.8867433071136475nan nan 6.7480845 nan nan nan 7.214833 nan nan nan 6.124478 nan nan nan 4.695984 -0.7315406799316406 5.720949928704791e+21 0.16360928 0.59283054-0.5486690402030945 0.5929527191533049 0.886202 nan nan 0.059893632115547144 nan nan nan nan -0.13614990460575904 nan nan -0.7326117157936096-0.652482807636261-0.66972786188125610 0 0 nan 0.004750422218261355 nan nan 125.17213110769076-168.2302942144889 -18.54608096646986 nan 202.63994832632554nan nan 210.507604177970930 -0.7563512325286865-1.09243309497833254863.5093275.22938484682527 275.22938484682527 nan nan 67.38685760498046 -300.3354850769043 55.13861343383789
1 0.03675363349150529 0.0308575008248499270.0391577021009943540.0218023930335774530.01185601080371138410.0163259506225590.3380742371082306-1.1830222606658936-0.39430123567581177nan -8.417383029591718 nan nan 0.0 22.89004704298827 -0.8867433071136475nan nan 5.6813674 nan nan nan 6.011444 nan nan nan 5.179238 nan nan nan 4.604853 -0.7315406799316406 2.9640027009278457e+200.41011643 0.69291294-0.5486690402030945 0.693232385460927 0.886202 nan nan -0.17179515341704701 nan nan nan nan -0.0962338436840866 nan nan -0.7326117157936096-0.652482807636261-0.66972786188125610 1 0 nan 0.0026423515125431757nan nan 150.0012221090916 -163.99231761527685306.3187042969622 nan 198.4922231483621 nan nan 378.4507834226541 0 -0.7563512325286865-1.09243309497833255579.93 73.8975744999208 73.8975744999208 nan nan 21.952759213430355 -344.88078691848267-104.08814365114452
2 3.5918606960695265 2.1166452934477245 3.3890797523932945 1.8598608456075434 1.1586647406675892 10.0163259506225590.3380742371082306-1.1830222606658936-0.39430123567581177nan -83.40311673017392 nan nan 0.0 17.57349665301026 -0.8867433071136475nan nan 11.699921 nan nan nan 12.905635 nan nan nan 10.669007 nan nan nan 5.164279 -0.7315406799316406 2.896661851668973e+22 0.00312047360.14802878-0.5486690402030945 0.148171417781676580.886202 nan nan 1.6732193473257568 nan nan nan nan 0.5649156080988607 nan nan -0.7326117157936096-0.652482807636261-0.66972786188125610 1 0 nan 0.030570368189346952 nan nan 12.434140598656441-25.78037465608216 -15.836701054884445nan 85.11047732952994 nan nan 32.71141498218776 0 -0.7563512325286865-1.09243309497833253348.2634407.1541354487033 407.1541354487033 nan nan -86.89760936346514 -448.4464746099155 -179.20189567325986
3 66.93287715606313 284.3907929593648 163.41445726792637 21.934311636022855 21.591250695504236 10.0163259506225590.3380742371082306-1.1830222606658936-0.39430123567581177nan 25.116892055647813 nan nan 0.0 21.819560732239374-0.8867433071136475nan nan 5.3372197 nan nan nan 5.6315928 nan nan nan 4.869831 nan nan nan 4.5574493-0.7315406799316406 5.397812673876058e+23 0.5614562 0.726579 -0.5486690402030945 0.7268972973128587 0.886202 nan nan -0.029328045437368597nan nan nan nan 0.3642257246132188 nan nan -0.7326117157936096-0.652482807636261-0.66972786188125610 1 0 nan 0.004326013332288946 nan nan 60.94197332301033 36.55344834119549 219.96530222118827 nan 212.6952906695271 nan nan 231.159712924621970 -0.7563512325286865-1.09243309497833255804.0195334.43550222465706 334.43550222465706 nan nan 390.8553665397576 284.69060588601087 195.85859030817775
4 0.17676035492413658 0.14116341642877317 0.18347508509139218 0.10354375093408806 0.05701946933036664 10.0163259506225590.3380742371082306-1.1830222606658936-0.39430123567581177nan -33.806606767270964nan nan 0.0 22.61096317412976 -0.8867433071136475nan nan 6.672208 nan nan nan 7.129779 nan nan nan 6.0562134 nan nan nan 4.68971 -0.7315406799316406 1.425486733259166e+21 0.17565803 0.60063034-0.5486690402030945 0.6007358197592728 0.886202 nan nan -0.13198050254089727 nan nan nan nan 0.19442476536089973 nan nan -0.7326117157936096-0.652482807636261-0.66972786188125610 1 0 nan 0.0030047432288688994nan nan 181.68874977198072-228.103924192566 160.3695429719637 nan 198.77232806002579nan nan 332.8071398554872 0 -0.7563512325286865-1.09243309497833254916.587 317.9421695907456 317.9421695907456 nan nan 466.9477690545565 -131.57910665429952-56.366331595593635
... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...
1,236,5220.0303838129888442980.02416082274495299 0.0314535457891423760.0177809992001699350.0098012299964013879.99021053314209 0.7036365270614624-0.89319157600402830.053123194724321365nan -70.40179275209826 nan nan 0.0 21.834277500479267-0.8060703873634338nan nan 6.846554 nan nan nan 7.335413 nan nan nan 6.2058263 nan nan nan 4.69534 -0.051267735660076142.4503074991003464e+200.15197778 0.5896 -0.102433867752552030.5905590121481706 0.8955961nan nan -0.39623652997028513 nan nan nan nan -0.1811989925496919 nan nan -0.56209683418273930.0669776126742363-1.214622616767883399999 1 0 nan 0.00429679367331618 nan nan 164.23541287200152-114.9885546206074 -118.18796569550966nan 298.05796999049454nan nan 232.7316776251488 0 -1.247308373451233 -0.99559468030929574777.316 142.47080473264475 142.47080473264475 nan nan -43.50419760572824 -482.63820143200155128.57036538752644
1,236,5232.8420045694329636 2.072874078803592 2.842470614227543 1.6171073022156517 0.9167756675590205 9.99021053314209 0.7036365270614624-0.89319157600402830.053123194724321365nan -8.743879258564638 nan nan 0.0 21.544367049554253-0.8060703873634338nan nan 6.3442564 nan nan nan 6.7639413 nan nan nan 5.762414 nan nan nan 4.661165 -0.051267735660076142.291939168897551e+22 0.2312045 0.6365434 -0.102433867752552030.6367230372102219 0.8955961nan nan -0.3826783814589959 nan nan nan nan 0.2714363153702788 nan nan -0.56209683418273930.0669776126742363-1.214622616767883399999 1 0 nan 0.00491050989022538 nan nan -152.9011299896323-134.10169435546433-10.449918743993287nan 104.03220288582199nan nan 203.644839813997920 -1.247308373451233 -0.99559468030929575102.662 405.6882227116895 405.6882227116895 nan nan -11.650050224520553-606.2818977234989 44.825464050083866
1,236,5240.14304244168711758 0.1208889298137058 0.15303129080388636 0.0849543001198953 0.04614272312487663 9.99021053314209 0.7036365270614624-0.89319157600402830.053123194724321365nan 27.20632557200984 nan nan 0.0 22.376794337207283-0.8060703873634338nan nan 5.32091 nan nan nan 5.621333 nan nan nan 4.849912 nan nan nan 4.5497336-0.051267735660076141.1535680781219158e+210.5675964 0.73500913-0.102433867752552030.7355422150028529 0.8955961nan nan -0.009462294142350725nan nan nan nan -0.057672877536177315nan nan -0.56209683418273930.0669776126742363-1.214622616767883399999 1 0 nan 0.003346887638984693 nan nan -231.1958588399615-73.37442862346133 174.46253428367396 nan 131.17425918434174nan nan 298.785052821001950 -1.247308373451233 -0.99559468030929575780.352 -233.94048014316843-233.94048014316843nan nan 132.63361502828016 50.96038956307919 -203.44990893403917
1,236,5250.17676035492413658 0.15172366208631505 0.190760602368683 0.10532176282817238 0.05701946933036664 9.99021053314209 0.7036365270614624-0.89319157600402830.053123194724321365nan -37.04737836059199 nan nan 0.0 22.10515257719927 -0.8060703873634338nan nan 3.429879 nan nan nan 3.6858425 nan nan nan 3.007895 nan nan nan 3.955146 -0.051267735660076141.425486733259166e+21 3.2290587 0.8733245 -0.102433867752552030.875238606457467 0.8955961nan nan -0.10573725549747043 nan nan nan nan 0.11806638977856074 nan nan -0.56209683418273930.0669776126742363-1.214622616767883399999 1 0 nan 0.0037928833357135567nan nan 145.02732041818132-188.33496882134352114.0579475068687 nan 198.1319965929415 nan nan 263.6516632568319 0 -1.247308373451233 -0.99559468030929576069.9087-114.39019236072886-114.39019236072886nan nan 82.28321437494888 114.36010574478351 -180.21077944859172
1,236,5260.0 0.0 0.0 0.0 0.0 9.99021053314209 0.7036365270614624-0.89319157600402830.053123194724321365nan 66.34501065058804 nan nan 0.0 22.365319084848508-0.8060703873634338nan nan 7.145454 nan nan nan 7.683845 nan nan nan 6.4657936 nan nan nan 4.7130475-0.051267735660076140.0 0.12192643 0.56545436-0.102433867752552030.5665924329936798 0.8955961nan nan 0.02282984335249972 nan nan nan nan -0.2731995547858111 nan nan -0.56209683418273930.0669776126742363-1.214622616767883399999 1 0 nan 0.0033646212465128435nan nan -44.87257181094431209.88325093720366 205.5952815364384 nan 232.72612631240358nan nan 297.2102732325128 0 -1.247308373451233 -0.99559468030929574615.3115-152.21224150477337-152.21224150477337nan nan -191.52291961963877-354.50122000134684100.05434225104445
" ], "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 minit 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 sulphur teff vr vr_Clean vr_Err vr_Sig vx vy vz\n", "0 0.7093977911593942 0.5527585393775436 0.7284652298464588 0.41253265460200594 0.22883799714819167 10.016325950622559 0.3380742371082306 -1.1830222606658936 -0.39430123567581177 nan -67.64321024253071 nan nan 0.0 21.61633894249464 -0.8867433071136475 nan nan 6.7480845 nan nan nan 7.214833 nan nan nan 6.124478 nan nan nan 4.695984 -0.7315406799316406 5.720949928704791e+21 0.16360928 0.59283054 -0.5486690402030945 0.5929527191533049 0.886202 nan nan 0.059893632115547144 nan nan nan nan -0.13614990460575904 nan nan -0.7326117157936096 -0.652482807636261 -0.6697278618812561 0 0 0 nan 0.004750422218261355 nan nan 125.17213110769076 -168.2302942144889 -18.54608096646986 nan 202.63994832632554 nan nan 210.50760417797093 0 -0.7563512325286865 -1.0924330949783325 4863.5093 275.22938484682527 275.22938484682527 nan nan 67.38685760498046 -300.3354850769043 55.13861343383789\n", "1 0.03675363349150529 0.030857500824849927 0.039157702100994354 0.021802393033577453 0.011856010803711384 10.016325950622559 0.3380742371082306 -1.1830222606658936 -0.39430123567581177 nan -8.417383029591718 nan nan 0.0 22.89004704298827 -0.8867433071136475 nan nan 5.6813674 nan nan nan 6.011444 nan nan nan 5.179238 nan nan nan 4.604853 -0.7315406799316406 2.9640027009278457e+20 0.41011643 0.69291294 -0.5486690402030945 0.693232385460927 0.886202 nan nan -0.17179515341704701 nan nan nan nan -0.0962338436840866 nan nan -0.7326117157936096 -0.652482807636261 -0.6697278618812561 0 1 0 nan 0.0026423515125431757 nan nan 150.0012221090916 -163.99231761527685 306.3187042969622 nan 198.4922231483621 nan nan 378.4507834226541 0 -0.7563512325286865 -1.0924330949783325 5579.93 73.8975744999208 73.8975744999208 nan nan 21.952759213430355 -344.88078691848267 -104.08814365114452\n", "2 3.5918606960695265 2.1166452934477245 3.3890797523932945 1.8598608456075434 1.1586647406675892 10.016325950622559 0.3380742371082306 -1.1830222606658936 -0.39430123567581177 nan -83.40311673017392 nan nan 0.0 17.57349665301026 -0.8867433071136475 nan nan 11.699921 nan nan nan 12.905635 nan nan nan 10.669007 nan nan nan 5.164279 -0.7315406799316406 2.896661851668973e+22 0.0031204736 0.14802878 -0.5486690402030945 0.14817141778167658 0.886202 nan nan 1.6732193473257568 nan nan nan nan 0.5649156080988607 nan nan -0.7326117157936096 -0.652482807636261 -0.6697278618812561 0 1 0 nan 0.030570368189346952 nan nan 12.434140598656441 -25.78037465608216 -15.836701054884445 nan 85.11047732952994 nan nan 32.71141498218776 0 -0.7563512325286865 -1.0924330949783325 3348.2634 407.1541354487033 407.1541354487033 nan nan -86.89760936346514 -448.4464746099155 -179.20189567325986\n", "3 66.93287715606313 284.3907929593648 163.41445726792637 21.934311636022855 21.591250695504236 10.016325950622559 0.3380742371082306 -1.1830222606658936 -0.39430123567581177 nan 25.116892055647813 nan nan 0.0 21.819560732239374 -0.8867433071136475 nan nan 5.3372197 nan nan nan 5.6315928 nan nan nan 4.869831 nan nan nan 4.5574493 -0.7315406799316406 5.397812673876058e+23 0.5614562 0.726579 -0.5486690402030945 0.7268972973128587 0.886202 nan nan -0.029328045437368597 nan nan nan nan 0.3642257246132188 nan nan -0.7326117157936096 -0.652482807636261 -0.6697278618812561 0 1 0 nan 0.004326013332288946 nan nan 60.94197332301033 36.55344834119549 219.96530222118827 nan 212.6952906695271 nan nan 231.15971292462197 0 -0.7563512325286865 -1.0924330949783325 5804.0195 334.43550222465706 334.43550222465706 nan nan 390.8553665397576 284.69060588601087 195.85859030817775\n", "4 0.17676035492413658 0.14116341642877317 0.18347508509139218 0.10354375093408806 0.05701946933036664 10.016325950622559 0.3380742371082306 -1.1830222606658936 -0.39430123567581177 nan -33.806606767270964 nan nan 0.0 22.61096317412976 -0.8867433071136475 nan nan 6.672208 nan nan nan 7.129779 nan nan nan 6.0562134 nan nan nan 4.68971 -0.7315406799316406 1.425486733259166e+21 0.17565803 0.60063034 -0.5486690402030945 0.6007358197592728 0.886202 nan nan -0.13198050254089727 nan nan nan nan 0.19442476536089973 nan nan -0.7326117157936096 -0.652482807636261 -0.6697278618812561 0 1 0 nan 0.0030047432288688994 nan nan 181.68874977198072 -228.103924192566 160.3695429719637 nan 198.77232806002579 nan nan 332.8071398554872 0 -0.7563512325286865 -1.0924330949783325 4916.587 317.9421695907456 317.9421695907456 nan nan 466.9477690545565 -131.57910665429952 -56.366331595593635\n", "... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...\n", "1,236,522 0.030383812988844298 0.02416082274495299 0.031453545789142376 0.017780999200169935 0.009801229996401387 9.99021053314209 0.7036365270614624 -0.8931915760040283 0.053123194724321365 nan -70.40179275209826 nan nan 0.0 21.834277500479267 -0.8060703873634338 nan nan 6.846554 nan nan nan 7.335413 nan nan nan 6.2058263 nan nan nan 4.69534 -0.05126773566007614 2.4503074991003464e+20 0.15197778 0.5896 -0.10243386775255203 0.5905590121481706 0.8955961 nan nan -0.39623652997028513 nan nan nan nan -0.1811989925496919 nan nan -0.5620968341827393 0.0669776126742363 -1.2146226167678833 99999 1 0 nan 0.00429679367331618 nan nan 164.23541287200152 -114.9885546206074 -118.18796569550966 nan 298.05796999049454 nan nan 232.7316776251488 0 -1.247308373451233 -0.9955946803092957 4777.316 142.47080473264475 142.47080473264475 nan nan -43.50419760572824 -482.63820143200155 128.57036538752644\n", "1,236,523 2.8420045694329636 2.072874078803592 2.842470614227543 1.6171073022156517 0.9167756675590205 9.99021053314209 0.7036365270614624 -0.8931915760040283 0.053123194724321365 nan -8.743879258564638 nan nan 0.0 21.544367049554253 -0.8060703873634338 nan nan 6.3442564 nan nan nan 6.7639413 nan nan nan 5.762414 nan nan nan 4.661165 -0.05126773566007614 2.291939168897551e+22 0.2312045 0.6365434 -0.10243386775255203 0.6367230372102219 0.8955961 nan nan -0.3826783814589959 nan nan nan nan 0.2714363153702788 nan nan -0.5620968341827393 0.0669776126742363 -1.2146226167678833 99999 1 0 nan 0.00491050989022538 nan nan -152.9011299896323 -134.10169435546433 -10.449918743993287 nan 104.03220288582199 nan nan 203.64483981399792 0 -1.247308373451233 -0.9955946803092957 5102.662 405.6882227116895 405.6882227116895 nan nan -11.650050224520553 -606.2818977234989 44.825464050083866\n", "1,236,524 0.14304244168711758 0.1208889298137058 0.15303129080388636 0.0849543001198953 0.04614272312487663 9.99021053314209 0.7036365270614624 -0.8931915760040283 0.053123194724321365 nan 27.20632557200984 nan nan 0.0 22.376794337207283 -0.8060703873634338 nan nan 5.32091 nan nan nan 5.621333 nan nan nan 4.849912 nan nan nan 4.5497336 -0.05126773566007614 1.1535680781219158e+21 0.5675964 0.73500913 -0.10243386775255203 0.7355422150028529 0.8955961 nan nan -0.009462294142350725 nan nan nan nan -0.057672877536177315 nan nan -0.5620968341827393 0.0669776126742363 -1.2146226167678833 99999 1 0 nan 0.003346887638984693 nan nan -231.1958588399615 -73.37442862346133 174.46253428367396 nan 131.17425918434174 nan nan 298.78505282100195 0 -1.247308373451233 -0.9955946803092957 5780.352 -233.94048014316843 -233.94048014316843 nan nan 132.63361502828016 50.96038956307919 -203.44990893403917\n", "1,236,525 0.17676035492413658 0.15172366208631505 0.190760602368683 0.10532176282817238 0.05701946933036664 9.99021053314209 0.7036365270614624 -0.8931915760040283 0.053123194724321365 nan -37.04737836059199 nan nan 0.0 22.10515257719927 -0.8060703873634338 nan nan 3.429879 nan nan nan 3.6858425 nan nan nan 3.007895 nan nan nan 3.955146 -0.05126773566007614 1.425486733259166e+21 3.2290587 0.8733245 -0.10243386775255203 0.875238606457467 0.8955961 nan nan -0.10573725549747043 nan nan nan nan 0.11806638977856074 nan nan -0.5620968341827393 0.0669776126742363 -1.2146226167678833 99999 1 0 nan 0.0037928833357135567 nan nan 145.02732041818132 -188.33496882134352 114.0579475068687 nan 198.1319965929415 nan nan 263.6516632568319 0 -1.247308373451233 -0.9955946803092957 6069.9087 -114.39019236072886 -114.39019236072886 nan nan 82.28321437494888 114.36010574478351 -180.21077944859172\n", "1,236,526 0.0 0.0 0.0 0.0 0.0 9.99021053314209 0.7036365270614624 -0.8931915760040283 0.053123194724321365 nan 66.34501065058804 nan nan 0.0 22.365319084848508 -0.8060703873634338 nan nan 7.145454 nan nan nan 7.683845 nan nan nan 6.4657936 nan nan nan 4.7130475 -0.05126773566007614 0.0 0.12192643 0.56545436 -0.10243386775255203 0.5665924329936798 0.8955961 nan nan 0.02282984335249972 nan nan nan nan -0.2731995547858111 nan nan -0.5620968341827393 0.0669776126742363 -1.2146226167678833 99999 1 0 nan 0.0033646212465128435 nan nan -44.87257181094431 209.88325093720366 205.5952815364384 nan 232.72612631240358 nan nan 297.2102732325128 0 -1.247308373451233 -0.9955946803092957 4615.3115 -152.21224150477337 -152.21224150477337 nan nan -191.52291961963877 -354.50122000134684 100.05434225104445" ] }, "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 minit 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 sulphur teff vr vr_Clean vr_Err vr_Sig vx vy vz
0 0.0 0.0 0.0 0.0 0.0 9.6704635620117190.9966661930084229 -1.2210134267807007 -0.866862952709198 9.379641502642725 9.37964184359441 -3.409516847471325e-07 4.0032654e-070.0 20.98134159368987 -1.4514390230178833 20.699148 0.008917205539981635 -0.29110986 0.014730277 20.60625 -0.5002811976175674 0.12518968 0.34471595 20.423008 0.3151946119749979 -0.8735288 0.34471595 21.218313926482217 32.789977594364814 2.210598 -0.8765081167221069 0.0 106.547681.074 -0.454772859811782841.07800815552214971.07894952.8218627555487608 4.390176195197285 -0.038668291846166214.428844487043452 2.2539077 3.5166128676565345 -1.02753513187815360.3161678891983761 -1.34370302107652972.2539077 -0.28822585940361023-0.952394425868988 -1.123855471611023 41 1 0 0.19151221663947868 0.006364022131264957 0.1851481945082137 1.4411755 123.14089310556496 79.32840839786846 56.87009735146566 263.7400068221982 263.7400069780021 -1.5580387532406174e-074.0032654e-07157.1333316217166 0 -0.6097681522369385 -0.326174199581146245063.5337-73.96488052073 -73.96488052073 nan nan -30.52740044242802570.35033465341655 -236.39764039403698
1 0.0726896631811602 0.0577266114969466060.07520941214860483 0.04252197804621342 0.023448278445535545 9.225929260253906-0.01319162081927061-1.3375555276870728 -0.29447129368782043 14.2927087364468 14.29270868966921 4.6777590423784194e-08 3.968168e-07 0.0 21.436982367710744-0.6304842233657837 20.683773 0.001921076973055727 -0.812856 0.014661002 21.105263 -0.08356368900664732 -0.32336566 0.34244165 20.407858 0.37898734021979 -1.4506352 0.34244165 75.53030028095948 -85.98998161429085 2.069631 -0.568984866142273 5.8620696113838857e+20177.871981.5886616-0.6436758637428284 1.59280976488731871.5488037-0.5145322630470786-0.0230524908400771 -0.214918292304910660.19186580146483356 2.2366896 -0.9859950123745546 -1.1119344390789918-0.15111201784831976-0.960822421230672 2.2366896 -0.626400351524353 -0.38219770789146423-0.7955443859100342 92 1 0 -3.1099475548599145 0.005159451406320136 -3.1151070062662347 1.4285405 3.3866957395069335 -48.31061257443122 187.67109898551485 185.69603212533264185.696032072546985.2785654222662604e-08 3.968168e-07 193.819055796326950 -1.6153771877288818 -1.2225399017333984 4808.5986-178.8177574749873 -178.8177574749873 nan nan -23.566758992497327-188.44689556565487-232.76072165508538
2 0.0 0.0 0.0 0.0 0.0 9.3016777038574220.415499746799469 -0.8815332651138306 -1.2968095541000366 41.934818492910786 41.93481806239753 4.305132521827234e-07 4.4025873e-070.0 17.452539050806816-0.8148350119590759 20.78876 0.008206054293774 3.3280153 0.015498512 20.906279 -0.0100527573335505463.463793 0.37021703 20.830341 0.3051409998318023 3.0726616 0.37021703 31.048341890855873 178.31954593582773 4.334337 -0.159046620130538940.0 3.607051 1.1337712-0.399335265159606931.13487866439809661.46496628.373033956499727 2.367282624058058 -1.6702253484205312 4.037507972478589 2.448507 -0.9190682013280174 8.08383306193433 4.818271968610433 3.265561093323897 2.448507 -0.7606052756309509 -0.6690134406089783 -0.3620690107345581 139 1 0 -0.9454505736129999 0.03232155077215023 -0.9777721243851502 1.5849314 -26.4951391573380870.7773107074618792 15.957189242705711 121.44684273727863121.446842581137631.561409923332064e-07 4.4025873e-0730.93910954488134 0 -1.03615140914917 -0.5998236536979675 7274.878 -367.46402582128326-367.46402582128326nan nan 640.65426316444 -383.57268749821895369.9504450871109
3 0.14868789260284618 0.12135829983653816 0.15609333325867777 0.08758908040221669 0.047963836323498765 8.982826232910156-0.3807957172393799 -0.7434932589530945 -0.042149968445301056-41.150652478427745-41.1506522476057 -2.308220473060798e-07 2.4726e-07 0.0 19.86556003372781 -0.3753343820571899420.22587 -0.0019696531326764540.24092163 0.011362129 20.75327 0.09069456594160508 0.64092165 0.23931459 19.820011 0.18694107591979076 -0.32007837 0.23931459 -60.16596236787171 -3.42740173719806762.7828066 -1.0834593772888184 1.199095908087469e+21 62.1827352.0512393-0.7561300992965698 2.051319213870563 1.96755831.7245934347929401 -2.7258924834137592 -0.400478742683286 -2.325413740730473 1.4821514 -2.3083182247726475 -0.9338445861152858-0.4090677587592231 -0.52477682735606271.4821514 0.11303850263357162 -0.17889174818992615-0.32079774141311646180 1 0 1.0922577799679152 0.01063868605740242 1.0816190939105128 0.8901359446.67865443015671 -2.7956262522421014-81.53921074635058 340.71437023828565340.7143702901641 -5.187844878001224e-08 2.4726e-07 93.9965701219464 0 -0.8397403359413147 -0.138868972659111025230.8975-388.19351059930756-388.19351059930756nan nan -1.9812087468119364-114.05739650504492450.27711082009785
4 0.0 0.0 0.0 0.0 0.0 9.473405838012695-0.348300039768219 -0.39555588364601135-0.8425897359848022 55.798487008887825 55.798487318897415-3.1000958964835713e-074.1201386e-070.0 20.73028154197521 -0.7897023558616638 20.727884 0.010290117189420935 -0.012687408 0.014958844 21.11538 0.006784309276688641 0.3783126 0.35224897 20.338593 0.1779984986904901 -0.5696874 0.35224897 58.734925197266485 104.14834870936562 2.476772 -0.217520266771316530.0 79.58692 1.298 -1.1380023956298828 1.32226961595118731.306232 1.4714652089146298 -0.9201217086903218 -0.5021471434900779 -0.4179745652002439 2.3111072 0.6781918589887956 -1.3336154242665539-0.41506939162943274-0.91854603263712112.3111072 -0.28932100534439087-0.4842122793197632 -0.3538527190685272 216 1 0 -3.7493629724747346 0.007144036941517213 -3.756507009416252 1.4832499 -17.75753198260054 70.44405532971666 119.6487860170618 210.15504747383238210.1550472993074 1.745249841360272e-07 4.1201386e-07139.976879765073730 -0.7041721343994141 -0.3010624349117279 5229.513 -416.9178793700232 -416.9178793700232 nan nan 77.08563863202212 -580.1798241091653 -134.72592647993883
... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...
1,6060.0 0.0 0.0 0.0 0.0 9.38451099395752 0.4432217478752136 -0.5602891445159912 -0.15777307748794556 -47.519042777241765-47.51904164012377-1.1371179972705115e-064.3653733e-070.0 18.27586522733421 -0.6387149095535278 20.749012 -0.02348805575141002 2.496635 0.015428401 20.775738 -0.1409814731332613 2.640854 0.36786798 20.43804 -0.068203521488147932.2303782 0.36786798 5.8527532895933625 -29.3844196386559134.0334554 -0.292838811874389650.0 7.60750771.2815708-0.195493146777153021.28214773869931051.3825252-3.2740209610150734-3.1106336961705394 -1.8346307809055318 -1.2760029152650076 2.4304676 -0.8525601334330479 1.330424760192722 -0.084664935980061861.4150896961727837 2.4304676 -0.8938306570053101 -0.7070809602737427 -0.6943075060844421 99656 1 0 2.650887127847258 0.022122130734428303 2.62876499711283 1.5715344 39.18271529980879 -22.0642892847665734.609514293273968 236.04343116372632236.043430770889423.928368963040253e-07 4.3653733e-0745.2036023114046940 -0.19964775443077087-0.595932126045227 7147.856 -46.578151367606665-46.578151367606665nan nan -138.8826103990196 -214.49165905196995-302.91878343717974
1,6071.319065373073896 0.9657225294100228 1.3160882286128504 0.752903994550527 0.4255049590560955 9.6946840286254880.3144649267196655 0.08044461905956268 -0.4594436585903168 -59.04860574167606 -59.04860565493767-8.673839898486664e-08 4.5912518e-070.0 21.60467767336341 -0.9855790138244629 20.81061 -0.009727128963152181-1.7500637 0.015849648 21.972767 0.21179698850423126 -1.1597953 0.3820421 19.682194 -0.20888238370122916-2.4665062 0.3820421 -32.30080556630673 -21.7667173392555361.332253 -0.9605326056480408 1.0637623976402386e+22469.386141.031 -0.6711140871047974 1.080351738633246 1.07197121.6902510513774607 3.7215264225566607 0.07520679762366758 3.646319624932993 2.5396707 3.7304781895975694 -1.7098705290045788-0.40183603058781336-1.30803449841676562.5396707 -0.4108820855617523 0.4820970594882965 -0.1632679104804992799678 1 0 0.638995171979518 0.004776001618846268 0.6342191703606718 1.6528506 164.36109754557685 -65.62902415411033 -111.88527018494862301.1050140034034 301.105013737890032.6551336875137766e-07 4.5912518e-07209.3801635355326 0 -0.8625441193580627 -1.1022558212280273 4467.9136539.8676637102133 539.8676637102133 nan nan 649.9855591531561 -176.8994063585329348.303906133631614
1,6080.16358290083015492 0.12371660081515007 0.1654368400221287 0.09432478442134579 0.05276867768714675 9.427410125732422-0.2979351878166199 -0.6779230237007141 -0.7961729168891907 -72.7101843268879 -72.71018421071417-1.1617372723931937e-073.7475579e-070.0 21.361563734776887-0.2609029710292816 20.668995 0.04161743326636356 -0.85790354 0.014218576 21.416006 0.11707942196517243 -0.22807339 0.32801843 19.521221 -0.3335890842186817 -1.6010762 0.32801843 -26.40912035420234 -76.43554906345355 1.7518698 -0.4149706959724426 1.3192169421786686e+21217.475561.415 -0.5588381290435791 1.4169202383032 1.408751 -1.5996339505161272-0.110694374543010840.2381306993703463 -0.348825073913357152.1280324 -0.07796421241912085-1.597702709989097 0.13583834735295117 -1.73354105734204822.1280324 -0.9929465651512146 -0.5758118033409119 -1.362631916999817 99690 0 0 -0.378467493181161630.005341795440224328 -0.383809288621385971.3491207 39.32437186086459 -162.9898795950431 -83.26371898221693 100.82329117453118100.823290855514643.190165371759531e-07 3.7475579e-07187.2029753273377 0 -0.20092321932315826-1.201646089553833 4337.5713187.42927928665588 187.42927928665588 nan nan -135.7710708618164 -271.5865379333496346.11163589477539
1,6090.15232777350090693 0.1089682081976194 0.15096974896059798 0.08613405815641825 0.049137991451905456 9.8438262939453120.12904362380504608 -0.6916096210479736 -0.6358122229576111 -64.27468659646041 -64.27468649583398-1.0062644001080644e-071.4157862e-070.0 21.565664005942352-0.4644570052623749 19.689575 0.011780971512784738 -1.9968392 0.008357561 20.560171 0.021424787517824746 -1.1778891 0.15522619 19.097677 0.3026562245371609 -2.8567772 0.15522619 -13.405302205946317-33.77399564406335 0.8676923 -0.260052889585495 1.2284497862976364e+21784.4015 0.972 -0.3354133665561676 1.03970750490846281.0328658-1.5401791476241433-1.0591716943137626 -0.1815668600789209 -0.8776048342348416 0.91319853 -0.194531372798605021.134966776493329 -0.212719467707689241.3476862442010185 0.91319853-0.7727305889129639 -0.3794154226779938 -0.6729522943496704 99839 1 0 0.12550882027637808 0.00486258491898170550.12064623535739638 0.509683 166.28796316364597 -111.21082703477833-47.6779175481562 255.26327329261466255.263273247606834.5007809269193156e-08 1.4157862e-07205.651935474972480 -1.2591712474822998 -0.5824271440505981 3945.6924144.6540682929069 144.6540682929069 nan nan -17.856087991760113-306.7362385631388429.253566910702535
1,6100.0063332968966788820.0042692271488931560.00622090421555935350.00348760810566237130.00204299899892867149.2706127166748050.5517753958702087 -0.5592101812362671 -0.7457828521728516 -46.350927080086294-46.350927116704223.661792446971074e-08 1.2829528e-070.0 22.621986534395614-0.46490278840065 19.584358 0.0062754730196515635-3.0481722 0.00790691 20.579638 -0.05893251145537814 -1.9896367 0.1435379 18.640377 0.039363469610319696-4.0244627 0.1435379 4.928791669064881 -87.72635226024357 0.42134246-0.6085554361343384 5.107497497321678e+19 2744.25591.60901370.0868726298213005 1.63866386410213651.58655880.6013688432415865 0.45123193933600814 0.0292358914625752970.4219960478734328 0.8383238 0.0028622142128001490.3975454877306359 0.0113519261329458680.38619356159769 0.8383238 -0.42113539576530457-1.5036002397537231 -0.668512225151062 99950 1 0 0.48053825214363877 0.00298952846624161950.47754872367739715 0.4618629813.221324967569547 -333.0016480281500428.73951046519076 145.46588050834504145.465880494448211.3896815611021882e-08 1.2829528e-07334.5009125325981 0 -0.9105316400527954 -0.228339418768882753680.93 293.80275593643427 293.80275593643427 nan nan -5.450666446110489 -289.1422501173651 71.83522477022365
" ], "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 minit 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 sulphur teff vr vr_Clean vr_Err vr_Sig vx vy vz\n", "0 0.0 0.0 0.0 0.0 0.0 9.670463562011719 0.9966661930084229 -1.2210134267807007 -0.866862952709198 9.379641502642725 9.37964184359441 -3.409516847471325e-07 4.0032654e-07 0.0 20.98134159368987 -1.4514390230178833 20.699148 0.008917205539981635 -0.29110986 0.014730277 20.60625 -0.5002811976175674 0.12518968 0.34471595 20.423008 0.3151946119749979 -0.8735288 0.34471595 21.218313926482217 32.789977594364814 2.210598 -0.8765081167221069 0.0 106.54768 1.074 -0.45477285981178284 1.0780081555221497 1.0789495 2.8218627555487608 4.390176195197285 -0.03866829184616621 4.428844487043452 2.2539077 3.5166128676565345 -1.0275351318781536 0.3161678891983761 -1.3437030210765297 2.2539077 -0.28822585940361023 -0.952394425868988 -1.123855471611023 41 1 0 0.19151221663947868 0.006364022131264957 0.1851481945082137 1.4411755 123.14089310556496 79.32840839786846 56.87009735146566 263.7400068221982 263.7400069780021 -1.5580387532406174e-07 4.0032654e-07 157.1333316217166 0 -0.6097681522369385 -0.32617419958114624 5063.5337 -73.96488052073 -73.96488052073 nan nan -30.527400442428025 70.35033465341655 -236.39764039403698\n", "1 0.0726896631811602 0.057726611496946606 0.07520941214860483 0.04252197804621342 0.023448278445535545 9.225929260253906 -0.01319162081927061 -1.3375555276870728 -0.29447129368782043 14.2927087364468 14.29270868966921 4.6777590423784194e-08 3.968168e-07 0.0 21.436982367710744 -0.6304842233657837 20.683773 0.001921076973055727 -0.812856 0.014661002 21.105263 -0.08356368900664732 -0.32336566 0.34244165 20.407858 0.37898734021979 -1.4506352 0.34244165 75.53030028095948 -85.98998161429085 2.069631 -0.568984866142273 5.8620696113838857e+20 177.87198 1.5886616 -0.6436758637428284 1.5928097648873187 1.5488037 -0.5145322630470786 -0.0230524908400771 -0.21491829230491066 0.19186580146483356 2.2366896 -0.9859950123745546 -1.1119344390789918 -0.15111201784831976 -0.960822421230672 2.2366896 -0.626400351524353 -0.38219770789146423 -0.7955443859100342 92 1 0 -3.1099475548599145 0.005159451406320136 -3.1151070062662347 1.4285405 3.3866957395069335 -48.31061257443122 187.67109898551485 185.69603212533264 185.69603207254698 5.2785654222662604e-08 3.968168e-07 193.81905579632695 0 -1.6153771877288818 -1.2225399017333984 4808.5986 -178.8177574749873 -178.8177574749873 nan nan -23.566758992497327 -188.44689556565487 -232.76072165508538\n", "2 0.0 0.0 0.0 0.0 0.0 9.301677703857422 0.415499746799469 -0.8815332651138306 -1.2968095541000366 41.934818492910786 41.93481806239753 4.305132521827234e-07 4.4025873e-07 0.0 17.452539050806816 -0.8148350119590759 20.78876 0.008206054293774 3.3280153 0.015498512 20.906279 -0.010052757333550546 3.463793 0.37021703 20.830341 0.3051409998318023 3.0726616 0.37021703 31.048341890855873 178.31954593582773 4.334337 -0.15904662013053894 0.0 3.607051 1.1337712 -0.39933526515960693 1.1348786643980966 1.4649662 8.373033956499727 2.367282624058058 -1.6702253484205312 4.037507972478589 2.448507 -0.9190682013280174 8.08383306193433 4.818271968610433 3.265561093323897 2.448507 -0.7606052756309509 -0.6690134406089783 -0.3620690107345581 139 1 0 -0.9454505736129999 0.03232155077215023 -0.9777721243851502 1.5849314 -26.495139157338087 0.7773107074618792 15.957189242705711 121.44684273727863 121.44684258113763 1.561409923332064e-07 4.4025873e-07 30.93910954488134 0 -1.03615140914917 -0.5998236536979675 7274.878 -367.46402582128326 -367.46402582128326 nan nan 640.65426316444 -383.57268749821895 369.9504450871109\n", "3 0.14868789260284618 0.12135829983653816 0.15609333325867777 0.08758908040221669 0.047963836323498765 8.982826232910156 -0.3807957172393799 -0.7434932589530945 -0.042149968445301056 -41.150652478427745 -41.1506522476057 -2.308220473060798e-07 2.4726e-07 0.0 19.86556003372781 -0.37533438205718994 20.22587 -0.001969653132676454 0.24092163 0.011362129 20.75327 0.09069456594160508 0.64092165 0.23931459 19.820011 0.18694107591979076 -0.32007837 0.23931459 -60.16596236787171 -3.4274017371980676 2.7828066 -1.0834593772888184 1.199095908087469e+21 62.182735 2.0512393 -0.7561300992965698 2.051319213870563 1.9675583 1.7245934347929401 -2.7258924834137592 -0.400478742683286 -2.325413740730473 1.4821514 -2.3083182247726475 -0.9338445861152858 -0.4090677587592231 -0.5247768273560627 1.4821514 0.11303850263357162 -0.17889174818992615 -0.32079774141311646 180 1 0 1.0922577799679152 0.01063868605740242 1.0816190939105128 0.89013594 46.67865443015671 -2.7956262522421014 -81.53921074635058 340.71437023828565 340.7143702901641 -5.187844878001224e-08 2.4726e-07 93.9965701219464 0 -0.8397403359413147 -0.13886897265911102 5230.8975 -388.19351059930756 -388.19351059930756 nan nan -1.9812087468119364 -114.05739650504492 450.27711082009785\n", "4 0.0 0.0 0.0 0.0 0.0 9.473405838012695 -0.348300039768219 -0.39555588364601135 -0.8425897359848022 55.798487008887825 55.798487318897415 -3.1000958964835713e-07 4.1201386e-07 0.0 20.73028154197521 -0.7897023558616638 20.727884 0.010290117189420935 -0.012687408 0.014958844 21.11538 0.006784309276688641 0.3783126 0.35224897 20.338593 0.1779984986904901 -0.5696874 0.35224897 58.734925197266485 104.14834870936562 2.476772 -0.21752026677131653 0.0 79.58692 1.298 -1.1380023956298828 1.3222696159511873 1.306232 1.4714652089146298 -0.9201217086903218 -0.5021471434900779 -0.4179745652002439 2.3111072 0.6781918589887956 -1.3336154242665539 -0.41506939162943274 -0.9185460326371211 2.3111072 -0.28932100534439087 -0.4842122793197632 -0.3538527190685272 216 1 0 -3.7493629724747346 0.007144036941517213 -3.756507009416252 1.4832499 -17.75753198260054 70.44405532971666 119.6487860170618 210.15504747383238 210.1550472993074 1.745249841360272e-07 4.1201386e-07 139.97687976507373 0 -0.7041721343994141 -0.3010624349117279 5229.513 -416.9178793700232 -416.9178793700232 nan nan 77.08563863202212 -580.1798241091653 -134.72592647993883\n", "... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...\n", "1,606 0.0 0.0 0.0 0.0 0.0 9.38451099395752 0.4432217478752136 -0.5602891445159912 -0.15777307748794556 -47.519042777241765 -47.51904164012377 -1.1371179972705115e-06 4.3653733e-07 0.0 18.27586522733421 -0.6387149095535278 20.749012 -0.02348805575141002 2.496635 0.015428401 20.775738 -0.1409814731332613 2.640854 0.36786798 20.43804 -0.06820352148814793 2.2303782 0.36786798 5.8527532895933625 -29.384419638655913 4.0334554 -0.29283881187438965 0.0 7.6075077 1.2815708 -0.19549314677715302 1.2821477386993105 1.3825252 -3.2740209610150734 -3.1106336961705394 -1.8346307809055318 -1.2760029152650076 2.4304676 -0.8525601334330479 1.330424760192722 -0.08466493598006186 1.4150896961727837 2.4304676 -0.8938306570053101 -0.7070809602737427 -0.6943075060844421 99656 1 0 2.650887127847258 0.022122130734428303 2.62876499711283 1.5715344 39.18271529980879 -22.064289284766573 4.609514293273968 236.04343116372632 236.04343077088942 3.928368963040253e-07 4.3653733e-07 45.203602311404694 0 -0.19964775443077087 -0.595932126045227 7147.856 -46.578151367606665 -46.578151367606665 nan nan -138.8826103990196 -214.49165905196995 -302.91878343717974\n", "1,607 1.319065373073896 0.9657225294100228 1.3160882286128504 0.752903994550527 0.4255049590560955 9.694684028625488 0.3144649267196655 0.08044461905956268 -0.4594436585903168 -59.04860574167606 -59.04860565493767 -8.673839898486664e-08 4.5912518e-07 0.0 21.60467767336341 -0.9855790138244629 20.81061 -0.009727128963152181 -1.7500637 0.015849648 21.972767 0.21179698850423126 -1.1597953 0.3820421 19.682194 -0.20888238370122916 -2.4665062 0.3820421 -32.30080556630673 -21.766717339255536 1.332253 -0.9605326056480408 1.0637623976402386e+22 469.38614 1.031 -0.6711140871047974 1.080351738633246 1.0719712 1.6902510513774607 3.7215264225566607 0.07520679762366758 3.646319624932993 2.5396707 3.7304781895975694 -1.7098705290045788 -0.40183603058781336 -1.3080344984167656 2.5396707 -0.4108820855617523 0.4820970594882965 -0.16326791048049927 99678 1 0 0.638995171979518 0.004776001618846268 0.6342191703606718 1.6528506 164.36109754557685 -65.62902415411033 -111.88527018494862 301.1050140034034 301.10501373789003 2.6551336875137766e-07 4.5912518e-07 209.3801635355326 0 -0.8625441193580627 -1.1022558212280273 4467.9136 539.8676637102133 539.8676637102133 nan nan 649.9855591531561 -176.89940635853293 48.303906133631614\n", "1,608 0.16358290083015492 0.12371660081515007 0.1654368400221287 0.09432478442134579 0.05276867768714675 9.427410125732422 -0.2979351878166199 -0.6779230237007141 -0.7961729168891907 -72.7101843268879 -72.71018421071417 -1.1617372723931937e-07 3.7475579e-07 0.0 21.361563734776887 -0.2609029710292816 20.668995 0.04161743326636356 -0.85790354 0.014218576 21.416006 0.11707942196517243 -0.22807339 0.32801843 19.521221 -0.3335890842186817 -1.6010762 0.32801843 -26.40912035420234 -76.43554906345355 1.7518698 -0.4149706959724426 1.3192169421786686e+21 217.47556 1.415 -0.5588381290435791 1.4169202383032 1.408751 -1.5996339505161272 -0.11069437454301084 0.2381306993703463 -0.34882507391335715 2.1280324 -0.07796421241912085 -1.597702709989097 0.13583834735295117 -1.7335410573420482 2.1280324 -0.9929465651512146 -0.5758118033409119 -1.362631916999817 99690 0 0 -0.37846749318116163 0.005341795440224328 -0.38380928862138597 1.3491207 39.32437186086459 -162.9898795950431 -83.26371898221693 100.82329117453118 100.82329085551464 3.190165371759531e-07 3.7475579e-07 187.2029753273377 0 -0.20092321932315826 -1.201646089553833 4337.5713 187.42927928665588 187.42927928665588 nan nan -135.7710708618164 -271.58653793334963 46.11163589477539\n", "1,609 0.15232777350090693 0.1089682081976194 0.15096974896059798 0.08613405815641825 0.049137991451905456 9.843826293945312 0.12904362380504608 -0.6916096210479736 -0.6358122229576111 -64.27468659646041 -64.27468649583398 -1.0062644001080644e-07 1.4157862e-07 0.0 21.565664005942352 -0.4644570052623749 19.689575 0.011780971512784738 -1.9968392 0.008357561 20.560171 0.021424787517824746 -1.1778891 0.15522619 19.097677 0.3026562245371609 -2.8567772 0.15522619 -13.405302205946317 -33.77399564406335 0.8676923 -0.260052889585495 1.2284497862976364e+21 784.4015 0.972 -0.3354133665561676 1.0397075049084628 1.0328658 -1.5401791476241433 -1.0591716943137626 -0.1815668600789209 -0.8776048342348416 0.91319853 -0.19453137279860502 1.134966776493329 -0.21271946770768924 1.3476862442010185 0.91319853 -0.7727305889129639 -0.3794154226779938 -0.6729522943496704 99839 1 0 0.12550882027637808 0.0048625849189817055 0.12064623535739638 0.509683 166.28796316364597 -111.21082703477833 -47.6779175481562 255.26327329261466 255.26327324760683 4.5007809269193156e-08 1.4157862e-07 205.65193547497248 0 -1.2591712474822998 -0.5824271440505981 3945.6924 144.6540682929069 144.6540682929069 nan nan -17.856087991760113 -306.73623856313884 29.253566910702535\n", "1,610 0.006333296896678882 0.004269227148893156 0.0062209042155593535 0.0034876081056623713 0.0020429989989286714 9.270612716674805 0.5517753958702087 -0.5592101812362671 -0.7457828521728516 -46.350927080086294 -46.35092711670422 3.661792446971074e-08 1.2829528e-07 0.0 22.621986534395614 -0.46490278840065 19.584358 0.0062754730196515635 -3.0481722 0.00790691 20.579638 -0.05893251145537814 -1.9896367 0.1435379 18.640377 0.039363469610319696 -4.0244627 0.1435379 4.928791669064881 -87.72635226024357 0.42134246 -0.6085554361343384 5.107497497321678e+19 2744.2559 1.6090137 0.0868726298213005 1.6386638641021365 1.5865588 0.6013688432415865 0.45123193933600814 0.029235891462575297 0.4219960478734328 0.8383238 0.002862214212800149 0.3975454877306359 0.011351926132945868 0.38619356159769 0.8383238 -0.42113539576530457 -1.5036002397537231 -0.668512225151062 99950 1 0 0.48053825214363877 0.0029895284662416195 0.47754872367739715 0.46186298 13.221324967569547 -333.00164802815004 28.73951046519076 145.46588050834504 145.46588049444821 1.3896815611021882e-08 1.2829528e-07 334.5009125325981 0 -0.9105316400527954 -0.22833941876888275 3680.93 293.80275593643427 293.80275593643427 nan nan -5.450666446110489 -289.1422501173651 71.83522477022365" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "survey[~survey.gaia__dr2_g.isna()]" ] } ], "metadata": { "kernelspec": { "display_name": "main", "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.10.13" } }, "nbformat": 4, "nbformat_minor": 5 }