jupyter/testing_ananke.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "f130bf6c",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'0.0.2b2'"
      ]
     },
     "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. Py-ananke has a method to produce such data."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "f4dafbb1",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "dict_keys(['pos3', 'vel3', 'mass', 'age', 'feh', 'helium', 'carbon', 'nitrogen', 'oxygen', 'neon', 'magnesium', 'silicon', 'sulphur', 'calcium', 'alpha', 'parentid', 'dform', 'log10_NH'])"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.random.seed(0)\n",
    "p = an.Ananke.make_dummy_particles_input()\n",
    "p.keys()"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "55c50b11",
   "metadata": {},
   "source": [
    "We can define here 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 choise `photo_sys`\n",
    " - the CMD `cmd_magnames` and its box limits `cmd_box`"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "1ecc74b8",
   "metadata": {},
   "outputs": [],
   "source": [
    "D = 5000 # *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'\n",
    "\n",
    "cmd_magnames = {'magnitude': 'Gmag',\n",
    "                'color_minuend': 'G_BPmag',\n",
    "                'color_subtrahend': 'G_RPmag'}\n",
    "\n",
    "cmd_box = {'abs_mag_lim_lo': -1000,\n",
    "           'abs_mag_lim_hi': 1000,\n",
    "           'app_mag_lim_hi': 30}"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "1f132dcb",
   "metadata": {},
   "source": [
    "Preparing now the pipeline with GAIA photometry"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "f56bb9c8",
   "metadata": {},
   "outputs": [],
   "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": [
    "This runs the pipeline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "f04f19ca",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Dimensions = 3\n",
      "Normalization constant of Kernel type 3: 0.596831\n",
      "Reading ASCII format file: to_enbid \n",
      "Read 300000 records\n",
      "Allocated 3.8147 MByte for particle storage.\n",
      "Type = 1 Particles = 100000\n",
      "Total particles = 100000\n",
      "Allocated 22.1252 MByte for Binary-Tree and Entropy bins.\n",
      "\n",
      "Scaling Co-ordinates as x[i]=x[i]/h[i] with h[i]->\n",
      "1 1 1 \n",
      "\n",
      "Starting to Build Tree .......\n",
      "Particle Type = 1 First node number = 1\n",
      "                  Last node number  = 199999\n",
      "Total number of nodes = 199999\n",
      "Treebuild time = 0.183057 s \n",
      "\n",
      "Density Calculation. Smoothing .....\n",
      "Evaluated =   1 % Time Left = 0.616863 s of 0.623100 s Par no = 96805 Density = 4.540412e-02 \n",
      "Evaluated =   2 % Time Left = 0.581330 s of 0.593200 s Par no = 55454 Density = 4.964532e-03 \n",
      "Evaluated =   3 % Time Left = 0.597449 s of 0.615933 s Par no = 27717 Density = 1.827217e-01 \n",
      "Evaluated =   4 % Time Left = 0.599490 s of 0.624475 s Par no = 42234 Density = 2.712376e-02 \n",
      "Evaluated =   5 % Time Left = 0.601838 s of 0.633520 s Par no = 99785 Density = 2.988042e-01 \n",
      "Evaluated =   6 % Time Left = 0.612027 s of 0.651100 s Par no = 88021 Density = 1.791133e-01 \n",
      "Evaluated =   7 % Time Left = 0.623957 s of 0.670929 s Par no = 75708 Density = 1.165057e-01 \n",
      "Evaluated =   8 % Time Left = 0.613150 s of 0.666475 s Par no = 45579 Density = 7.255409e-02 \n",
      "Evaluated =   9 % Time Left = 0.599815 s of 0.659144 s Par no = 30351 Density = 8.828936e-02 \n",
      "Evaluated =  10 % Time Left = 0.597764 s of 0.664190 s Par no = 89176 Density = 1.775999e-01 \n",
      "Evaluated =  11 % Time Left = 0.596010 s of 0.669682 s Par no = 60361 Density = 8.299880e-02 \n",
      "Evaluated =  12 % Time Left = 0.594521 s of 0.675600 s Par no = 90082 Density = 2.203585e-01 \n",
      "Evaluated =  13 % Time Left = 0.588562 s of 0.676515 s Par no = 26814 Density = 1.910994e-02 \n",
      "Evaluated =  14 % Time Left = 0.579781 s of 0.674171 s Par no = 74297 Density = 2.157651e-02 \n",
      "Evaluated =  15 % Time Left = 0.574707 s of 0.676133 s Par no = 94905 Density = 9.119636e-02 \n",
      "Evaluated =  16 % Time Left = 0.571382 s of 0.680225 s Par no = 51653 Density = 9.666828e-02 \n",
      "Evaluated =  17 % Time Left = 0.563256 s of 0.678629 s Par no = 39308 Density = 1.287775e-01 \n",
      "Evaluated =  18 % Time Left = 0.559425 s of 0.682233 s Par no = 7800 Density = 4.130107e-02 \n",
      "Evaluated =  19 % Time Left = 0.550226 s of 0.679300 s Par no = 42838 Density = 5.914273e-02 \n",
      "Evaluated =  20 % Time Left = 0.542649 s of 0.678320 s Par no = 49016 Density = 1.554842e-01 \n",
      "Evaluated =  21 % Time Left = 0.534601 s of 0.676719 s Par no = 29173 Density = 5.981782e-02 \n",
      "Evaluated =  22 % Time Left = 0.523926 s of 0.671709 s Par no = 14018 Density = 4.713365e-02 \n",
      "Evaluated =  23 % Time Left = 0.515465 s of 0.669443 s Par no = 90617 Density = 4.088211e-02 \n",
      "Evaluated =  24 % Time Left = 0.511362 s of 0.672854 s Par no = 20588 Density = 1.273881e-01 \n",
      "Evaluated =  25 % Time Left = 0.506177 s of 0.674912 s Par no = 15465 Density = 1.274245e-01 \n",
      "Evaluated =  26 % Time Left = 0.501995 s of 0.678381 s Par no = 99360 Density = 2.044210e-01 \n",
      "Evaluated =  27 % Time Left = 0.495990 s of 0.679448 s Par no = 71132 Density = 2.079416e-02 \n",
      "Evaluated =  28 % Time Left = 0.489498 s of 0.679868 s Par no = 6665 Density = 4.366050e-03 \n",
      "Evaluated =  29 % Time Left = 0.485476 s of 0.683779 s Par no = 44439 Density = 1.463732e-01 \n",
      "Evaluated =  30 % Time Left = 0.478289 s of 0.683280 s Par no = 6386 Density = 1.825201e-01 \n",
      "Evaluated =  31 % Time Left = 0.472015 s of 0.684090 s Par no = 2716 Density = 1.495581e-01 \n",
      "Evaluated =  32 % Time Left = 0.465253 s of 0.684206 s Par no = 69470 Density = 6.785187e-02 \n",
      "Evaluated =  33 % Time Left = 0.459497 s of 0.685827 s Par no = 95884 Density = 1.793627e-01 \n",
      "Evaluated =  34 % Time Left = 0.456265 s of 0.691321 s Par no = 73768 Density = 1.295312e-02 \n",
      "Evaluated =  35 % Time Left = 0.451578 s of 0.694746 s Par no = 47898 Density = 1.088661e-01 \n",
      "Evaluated =  36 % Time Left = 0.448923 s of 0.701453 s Par no = 66396 Density = 3.959185e-02 \n",
      "Evaluated =  37 % Time Left = 0.445665 s of 0.707416 s Par no = 16463 Density = 1.572218e-01 \n",
      "Evaluated =  38 % Time Left = 0.441071 s of 0.711416 s Par no = 55669 Density = 1.712036e-01 \n",
      "Evaluated =  39 % Time Left = 0.435561 s of 0.714046 s Par no = 97705 Density = 1.234432e-01 \n",
      "Evaluated =  40 % Time Left = 0.428115 s of 0.713538 s Par no = 43442 Density = 6.719715e-02 \n",
      "Evaluated =  41 % Time Left = 0.420981 s of 0.713539 s Par no = 10242 Density = 2.102901e-01 \n",
      "Evaluated =  42 % Time Left = 0.414709 s of 0.715029 s Par no = 40929 Density = 1.391802e-01 \n",
      "Evaluated =  43 % Time Left = 0.407640 s of 0.715170 s Par no = 5720 Density = 6.551556e-02 \n",
      "Evaluated =  44 % Time Left = 0.400033 s of 0.714357 s Par no = 66730 Density = 9.093652e-02 \n",
      "Evaluated =  45 % Time Left = 0.393590 s of 0.715631 s Par no = 48819 Density = 2.248661e-01 \n",
      "Evaluated =  46 % Time Left = 0.386977 s of 0.716637 s Par no = 32014 Density = 9.211309e-02 \n",
      "Evaluated =  47 % Time Left = 0.379441 s of 0.715940 s Par no = 31533 Density = 2.932794e-02 \n",
      "Evaluated =  48 % Time Left = 0.371361 s of 0.714169 s Par no = 86677 Density = 1.263195e-01 \n",
      "Evaluated =  49 % Time Left = 0.364454 s of 0.714631 s Par no = 81951 Density = 8.802209e-02 \n",
      "Evaluated =  50 % Time Left = 0.357381 s of 0.714776 s Par no = 80149 Density = 1.559105e-01 \n",
      "Evaluated =  51 % Time Left = 0.350363 s of 0.715041 s Par no = 25229 Density = 5.093941e-02 \n",
      "Evaluated =  52 % Time Left = 0.343227 s of 0.715071 s Par no = 89731 Density = 4.987880e-02 \n",
      "Evaluated =  53 % Time Left = 0.334934 s of 0.712642 s Par no = 76366 Density = 1.511493e-01 \n",
      "Evaluated =  54 % Time Left = 0.327808 s of 0.712641 s Par no = 73159 Density = 1.548675e-02 \n",
      "Evaluated =  55 % Time Left = 0.320680 s of 0.712638 s Par no = 54007 Density = 8.763472e-02 \n",
      "Evaluated =  56 % Time Left = 0.313058 s of 0.711511 s Par no = 36235 Density = 5.465963e-02 \n",
      "Evaluated =  57 % Time Left = 0.305663 s of 0.710861 s Par no = 62068 Density = 4.050222e-02 \n",
      "Evaluated =  58 % Time Left = 0.297614 s of 0.708621 s Par no = 37661 Density = 2.656842e-02 \n",
      "Evaluated =  59 % Time Left = 0.290793 s of 0.709268 s Par no = 67914 Density = 1.940649e-01 \n",
      "Evaluated =  60 % Time Left = 0.283896 s of 0.709757 s Par no = 64732 Density = 5.859732e-02 \n",
      "Evaluated =  61 % Time Left = 0.276151 s of 0.708098 s Par no = 88788 Density = 7.161248e-03 \n",
      "Evaluated =  62 % Time Left = 0.270871 s of 0.712837 s Par no = 44083 Density = 2.877788e-02 \n",
      "Evaluated =  63 % Time Left = 0.264371 s of 0.714537 s Par no = 14045 Density = 3.265437e-02 \n",
      "Evaluated =  64 % Time Left = 0.257069 s of 0.714100 s Par no = 93621 Density = 5.216300e-02 \n",
      "Evaluated =  65 % Time Left = 0.251367 s of 0.718211 s Par no = 9132 Density = 8.411559e-02 \n",
      "Evaluated =  66 % Time Left = 0.246055 s of 0.723712 s Par no = 702 Density = 1.110411e-02 \n",
      "Evaluated =  67 % Time Left = 0.239801 s of 0.726693 s Par no = 45670 Density = 2.465221e-02 \n",
      "Evaluated =  68 % Time Left = 0.232487 s of 0.726546 s Par no = 53860 Density = 1.580593e-02 \n",
      "Evaluated =  69 % Time Left = 0.225279 s of 0.726729 s Par no = 79719 Density = 4.454535e-02 \n",
      "Evaluated =  70 % Time Left = 0.217550 s of 0.725191 s Par no = 96401 Density = 1.035200e-02 \n",
      "Evaluated =  71 % Time Left = 0.210009 s of 0.724194 s Par no = 48105 Density = 3.704986e-02 \n",
      "Evaluated =  72 % Time Left = 0.202603 s of 0.723608 s Par no = 75056 Density = 2.642116e-02 \n",
      "Evaluated =  73 % Time Left = 0.195117 s of 0.722682 s Par no = 97912 Density = 1.775965e-01 \n",
      "Evaluated =  74 % Time Left = 0.187544 s of 0.721353 s Par no = 87879 Density = 7.451770e-02 \n",
      "Evaluated =  75 % Time Left = 0.180503 s of 0.722043 s Par no = 85159 Density = 1.462781e-02 \n",
      "Evaluated =  76 % Time Left = 0.173254 s of 0.721921 s Par no = 78132 Density = 5.302112e-02 \n",
      "Evaluated =  77 % Time Left = 0.165699 s of 0.720461 s Par no = 21990 Density = 3.497640e-03 \n",
      "Evaluated =  78 % Time Left = 0.158175 s of 0.719012 s Par no = 16462 Density = 2.482335e-02 \n",
      "Evaluated =  79 % Time Left = 0.151059 s of 0.719365 s Par no = 22311 Density = 3.571432e-02 \n",
      "Evaluated =  80 % Time Left = 0.143648 s of 0.718275 s Par no = 43269 Density = 1.776266e-02 \n",
      "Evaluated =  81 % Time Left = 0.136401 s of 0.717940 s Par no = 43142 Density = 1.532302e-02 \n",
      "Evaluated =  82 % Time Left = 0.129060 s of 0.717038 s Par no = 52197 Density = 3.049043e-02 \n",
      "Evaluated =  83 % Time Left = 0.121641 s of 0.715576 s Par no = 74476 Density = 1.449572e-01 \n",
      "Evaluated =  84 % Time Left = 0.114357 s of 0.714774 s Par no = 8891 Density = 1.050831e-02 \n",
      "Evaluated =  85 % Time Left = 0.107208 s of 0.714768 s Par no = 1559 Density = 1.452041e-03 \n",
      "Evaluated =  86 % Time Left = 0.100040 s of 0.714626 s Par no = 20602 Density = 2.028736e-01 \n",
      "Evaluated =  87 % Time Left = 0.093209 s of 0.717046 s Par no = 54383 Density = 2.137116e-01 \n",
      "Evaluated =  88 % Time Left = 0.086471 s of 0.720651 s Par no = 70246 Density = 2.757094e-03 \n",
      "Evaluated =  89 % Time Left = 0.079635 s of 0.724019 s Par no = 51472 Density = 5.058417e-04 \n",
      "Evaluated =  90 % Time Left = 0.072548 s of 0.725549 s Par no = 66269 Density = 3.630537e-02 \n",
      "Evaluated =  91 % Time Left = 0.065477 s of 0.727601 s Par no = 62858 Density = 1.776032e-01 \n",
      "Evaluated =  92 % Time Left = 0.058494 s of 0.731264 s Par no = 84061 Density = 1.439861e-01 \n",
      "Evaluated =  93 % Time Left = 0.051510 s of 0.735969 s Par no = 78610 Density = 4.029672e-03 \n",
      "Evaluated =  94 % Time Left = 0.044449 s of 0.740934 s Par no = 3871 Density = 1.597447e-01 \n",
      "Evaluated =  95 % Time Left = 0.037172 s of 0.743594 s Par no = 68736 Density = 1.582606e-01 \n",
      "Evaluated =  96 % Time Left = 0.029821 s of 0.745711 s Par no = 30472 Density = 6.743144e-02 \n",
      "Evaluated =  97 % Time Left = 0.022374 s of 0.746060 s Par no = 86838 Density = 1.353215e-03 \n",
      "Evaluated =  98 % Time Left = 0.014940 s of 0.747395 s Par no = 92913 Density = 5.211163e-02 \n",
      "Evaluated =  99 % Time Left = 0.007479 s of 0.748625 s Par no = 30015 Density = 9.179294e-04 \n",
      "Evaluated = 100 % Time Left = -0.000007 s of 0.749219 s Par no = 86525 Density = 1.579017e-04 \n",
      "\n",
      "Total Smoothing Time = 0.749258 s\n",
      "1 3\n",
      "\n",
      "writing snapshot file ....\n",
      "to_enbid_d3n64.est\n",
      "done with snapshot.\n",
      "\n",
      "Total Time = 1.1589 s \n",
      "\n",
      "Dimensions = 3\n",
      "Normalization constant of Kernel type 3: 0.596831\n",
      "Reading ASCII format file: to_enbid \n",
      "Read 300000 records\n",
      "Allocated 3.8147 MByte for particle storage.\n",
      "Type = 1 Particles = 100000\n",
      "Total particles = 100000\n",
      "Allocated 22.1252 MByte for Binary-Tree and Entropy bins.\n",
      "\n",
      "Scaling Co-ordinates as x[i]=x[i]/h[i] with h[i]->\n",
      "1 1 1 \n",
      "\n",
      "Starting to Build Tree .......\n",
      "Particle Type = 1 First node number = 1\n",
      "                  Last node number  = 199999\n",
      "Total number of nodes = 199999\n",
      "Treebuild time = 0.181858 s \n",
      "\n",
      "Density Calculation. Smoothing .....\n",
      "Evaluated =   1 % Time Left = 0.903564 s of 0.912700 s Par no = 3379 Density = 1.586355e-02 \n",
      "Evaluated =   2 % Time Left = 0.926287 s of 0.945200 s Par no = 75725 Density = 4.756319e-02 \n",
      "Evaluated =   3 % Time Left = 0.924724 s of 0.953333 s Par no = 95160 Density = 6.928091e-03 \n",
      "Evaluated =   4 % Time Left = 0.869343 s of 0.905575 s Par no = 40983 Density = 2.235913e-02 \n",
      "Evaluated =   5 % Time Left = 0.839905 s of 0.884120 s Par no = 61196 Density = 3.542700e-02 \n",
      "Evaluated =   6 % Time Left = 0.835635 s of 0.888983 s Par no = 68808 Density = 4.706392e-02 \n",
      "Evaluated =   7 % Time Left = 0.863761 s of 0.928786 s Par no = 54967 Density = 2.581413e-03 \n",
      "Evaluated =   8 % Time Left = 0.828531 s of 0.900588 s Par no = 92704 Density = 1.288175e-02 \n",
      "Evaluated =   9 % Time Left = 0.831498 s of 0.913744 s Par no = 40727 Density = 3.485199e-02 \n",
      "Evaluated =  10 % Time Left = 0.829917 s of 0.922140 s Par no = 56891 Density = 2.531235e-02 \n",
      "Evaluated =  11 % Time Left = 0.824721 s of 0.926664 s Par no = 5777 Density = 3.864953e-02 \n",
      "Evaluated =  12 % Time Left = 0.812003 s of 0.922742 s Par no = 4039 Density = 3.381603e-02 \n",
      "Evaluated =  13 % Time Left = 0.837211 s of 0.962323 s Par no = 99176 Density = 2.977309e-02 \n",
      "Evaluated =  14 % Time Left = 0.838982 s of 0.975571 s Par no = 23426 Density = 2.395511e-02 \n",
      "Evaluated =  15 % Time Left = 0.820864 s of 0.965733 s Par no = 18363 Density = 1.992155e-02 \n",
      "Evaluated =  16 % Time Left = 0.794263 s of 0.945563 s Par no = 87665 Density = 1.640359e-03 \n",
      "Evaluated =  17 % Time Left = 0.776529 s of 0.935588 s Par no = 14016 Density = 7.559446e-03 \n",
      "Evaluated =  18 % Time Left = 0.752546 s of 0.917750 s Par no = 54692 Density = 2.281061e-03 \n",
      "Evaluated =  19 % Time Left = 0.736746 s of 0.909574 s Par no = 17861 Density = 3.635737e-02 \n",
      "Evaluated =  20 % Time Left = 0.728015 s of 0.910030 s Par no = 67475 Density = 4.535319e-02 \n",
      "Evaluated =  21 % Time Left = 0.717815 s of 0.908638 s Par no = 8904 Density = 1.495583e-02 \n",
      "Evaluated =  22 % Time Left = 0.713063 s of 0.914195 s Par no = 83588 Density = 2.024327e-02 \n",
      "Evaluated =  23 % Time Left = 0.701384 s of 0.910900 s Par no = 31865 Density = 3.190889e-02 \n",
      "Evaluated =  24 % Time Left = 0.686423 s of 0.903200 s Par no = 48390 Density = 2.340200e-02 \n",
      "Evaluated =  25 % Time Left = 0.676164 s of 0.901564 s Par no = 62455 Density = 6.559847e-03 \n",
      "Evaluated =  26 % Time Left = 0.673160 s of 0.909688 s Par no = 40761 Density = 4.604229e-02 \n",
      "Evaluated =  27 % Time Left = 0.664905 s of 0.910841 s Par no = 76737 Density = 5.852960e-03 \n",
      "Evaluated =  28 % Time Left = 0.650706 s of 0.903771 s Par no = 54978 Density = 3.478805e-03 \n",
      "Evaluated =  29 % Time Left = 0.642810 s of 0.905379 s Par no = 96011 Density = 1.743980e-02 \n",
      "Evaluated =  30 % Time Left = 0.638242 s of 0.911787 s Par no = 36915 Density = 2.532993e-02 \n",
      "Evaluated =  31 % Time Left = 0.640128 s of 0.927735 s Par no = 588 Density = 2.564653e-02 \n",
      "Evaluated =  32 % Time Left = 0.631726 s of 0.929022 s Par no = 31595 Density = 5.423155e-03 \n",
      "Evaluated =  33 % Time Left = 0.624636 s of 0.932306 s Par no = 65882 Density = 9.776751e-03 \n",
      "Evaluated =  34 % Time Left = 0.612962 s of 0.928744 s Par no = 36652 Density = 4.575477e-03 \n",
      "Evaluated =  35 % Time Left = 0.601842 s of 0.925926 s Par no = 78212 Density = 1.707879e-02 \n",
      "Evaluated =  36 % Time Left = 0.591497 s of 0.924228 s Par no = 98062 Density = 2.259785e-02 \n",
      "Evaluated =  37 % Time Left = 0.578486 s of 0.918246 s Par no = 99847 Density = 1.371562e-02 \n",
      "Evaluated =  38 % Time Left = 0.566103 s of 0.913084 s Par no = 40933 Density = 3.523312e-02 \n",
      "Evaluated =  39 % Time Left = 0.555006 s of 0.909862 s Par no = 2518 Density = 3.364817e-02 \n",
      "Evaluated =  40 % Time Left = 0.544866 s of 0.908125 s Par no = 78254 Density = 6.255044e-02 \n",
      "Evaluated =  41 % Time Left = 0.537533 s of 0.911088 s Par no = 88151 Density = 4.831521e-02 \n",
      "Evaluated =  42 % Time Left = 0.526229 s of 0.907307 s Par no = 1224 Density = 3.085156e-02 \n",
      "Evaluated =  43 % Time Left = 0.519343 s of 0.911144 s Par no = 85166 Density = 1.231266e-02 \n",
      "Evaluated =  44 % Time Left = 0.508700 s of 0.908409 s Par no = 90763 Density = 2.646071e-02 \n",
      "Evaluated =  45 % Time Left = 0.497701 s of 0.904927 s Par no = 76867 Density = 5.733746e-03 \n",
      "Evaluated =  46 % Time Left = 0.487112 s of 0.902076 s Par no = 65993 Density = 2.850425e-02 \n",
      "Evaluated =  47 % Time Left = 0.476793 s of 0.899626 s Par no = 37327 Density = 7.533906e-03 \n",
      "Evaluated =  48 % Time Left = 0.466131 s of 0.896423 s Par no = 23978 Density = 7.339587e-03 \n",
      "Evaluated =  49 % Time Left = 0.455841 s of 0.893822 s Par no = 98943 Density = 1.799383e-03 \n",
      "Evaluated =  50 % Time Left = 0.445164 s of 0.890346 s Par no = 54150 Density = 2.290839e-02 \n",
      "Evaluated =  51 % Time Left = 0.437084 s of 0.892025 s Par no = 7580 Density = 1.406420e-02 \n",
      "Evaluated =  52 % Time Left = 0.430000 s of 0.895852 s Par no = 65512 Density = 4.249471e-03 \n",
      "Evaluated =  53 % Time Left = 0.420528 s of 0.894760 s Par no = 63408 Density = 3.415667e-03 \n",
      "Evaluated =  54 % Time Left = 0.411903 s of 0.895461 s Par no = 56448 Density = 1.298454e-02 \n",
      "Evaluated =  55 % Time Left = 0.401227 s of 0.891636 s Par no = 63701 Density = 1.132553e-03 \n",
      "Evaluated =  56 % Time Left = 0.391437 s of 0.889650 s Par no = 95209 Density = 1.302888e-02 \n",
      "Evaluated =  57 % Time Left = 0.382989 s of 0.890693 s Par no = 57672 Density = 4.024602e-02 \n",
      "Evaluated =  58 % Time Left = 0.372475 s of 0.886867 s Par no = 21270 Density = 1.373904e-02 \n",
      "Evaluated =  59 % Time Left = 0.361748 s of 0.882334 s Par no = 7187 Density = 2.129810e-02 \n",
      "Evaluated =  60 % Time Left = 0.351558 s of 0.878917 s Par no = 57532 Density = 2.673977e-02 \n",
      "Evaluated =  61 % Time Left = 0.341650 s of 0.876048 s Par no = 20643 Density = 3.174773e-02 \n",
      "Evaluated =  62 % Time Left = 0.332333 s of 0.874582 s Par no = 8171 Density = 3.570614e-02 \n",
      "Evaluated =  63 % Time Left = 0.323123 s of 0.873329 s Par no = 70708 Density = 2.276234e-02 \n",
      "Evaluated =  64 % Time Left = 0.314735 s of 0.874287 s Par no = 86426 Density = 3.718382e-04 \n",
      "Evaluated =  65 % Time Left = 0.305075 s of 0.871668 s Par no = 20171 Density = 9.558010e-03 \n",
      "Evaluated =  66 % Time Left = 0.296208 s of 0.871224 s Par no = 79285 Density = 7.091322e-03 \n",
      "Evaluated =  67 % Time Left = 0.287110 s of 0.870057 s Par no = 76449 Density = 4.521368e-03 \n",
      "Evaluated =  68 % Time Left = 0.277699 s of 0.867837 s Par no = 90892 Density = 4.311953e-03 \n",
      "Evaluated =  69 % Time Left = 0.269504 s of 0.869397 s Par no = 77427 Density = 1.956905e-02 \n",
      "Evaluated =  70 % Time Left = 0.259969 s of 0.866591 s Par no = 63619 Density = 1.020839e-03 \n",
      "Evaluated =  71 % Time Left = 0.250305 s of 0.863152 s Par no = 23141 Density = 1.147707e-02 \n",
      "Evaluated =  72 % Time Left = 0.240902 s of 0.860394 s Par no = 97062 Density = 2.525904e-03 \n",
      "Evaluated =  73 % Time Left = 0.231420 s of 0.857142 s Par no = 20151 Density = 3.074151e-03 \n",
      "Evaluated =  74 % Time Left = 0.222578 s of 0.856101 s Par no = 16518 Density = 1.823752e-03 \n",
      "Evaluated =  75 % Time Left = 0.213493 s of 0.854007 s Par no = 30015 Density = 1.477548e-03 \n",
      "Evaluated =  76 % Time Left = 0.204995 s of 0.854182 s Par no = 69971 Density = 3.932012e-02 \n",
      "Evaluated =  77 % Time Left = 0.196378 s of 0.853853 s Par no = 92629 Density = 2.091704e-02 \n",
      "Evaluated =  78 % Time Left = 0.187942 s of 0.854321 s Par no = 83958 Density = 3.231096e-02 \n",
      "Evaluated =  79 % Time Left = 0.179384 s of 0.854251 s Par no = 1207 Density = 6.786271e-03 \n",
      "Evaluated =  80 % Time Left = 0.170778 s of 0.853935 s Par no = 75504 Density = 8.697630e-03 \n",
      "Evaluated =  81 % Time Left = 0.162165 s of 0.853543 s Par no = 9731 Density = 4.173785e-03 \n",
      "Evaluated =  82 % Time Left = 0.154118 s of 0.856257 s Par no = 51510 Density = 2.677194e-03 \n",
      "Evaluated =  83 % Time Left = 0.145800 s of 0.857696 s Par no = 36411 Density = 2.679306e-02 \n",
      "Evaluated =  84 % Time Left = 0.137023 s of 0.856448 s Par no = 32995 Density = 2.042019e-02 \n",
      "Evaluated =  85 % Time Left = 0.128157 s of 0.854438 s Par no = 1296 Density = 1.403352e-02 \n",
      "Evaluated =  86 % Time Left = 0.119589 s of 0.854270 s Par no = 15953 Density = 4.376418e-02 \n",
      "Evaluated =  87 % Time Left = 0.111154 s of 0.855100 s Par no = 67362 Density = 3.897928e-02 \n",
      "Evaluated =  88 % Time Left = 0.102504 s of 0.854274 s Par no = 1175 Density = 5.999018e-03 \n",
      "Evaluated =  89 % Time Left = 0.093843 s of 0.853196 s Par no = 35183 Density = 1.744862e-02 \n",
      "Evaluated =  90 % Time Left = 0.085318 s of 0.853268 s Par no = 57201 Density = 2.958058e-02 \n",
      "Evaluated =  91 % Time Left = 0.076902 s of 0.854562 s Par no = 35155 Density = 1.742160e-02 \n",
      "Evaluated =  92 % Time Left = 0.068342 s of 0.854379 s Par no = 98988 Density = 3.247329e-03 \n",
      "Evaluated =  93 % Time Left = 0.059719 s of 0.853247 s Par no = 7049 Density = 1.366043e-02 \n",
      "Evaluated =  94 % Time Left = 0.051279 s of 0.854797 s Par no = 63956 Density = 1.613105e-02 \n",
      "Evaluated =  95 % Time Left = 0.042791 s of 0.855998 s Par no = 22220 Density = 3.321449e-02 \n",
      "Evaluated =  96 % Time Left = 0.034418 s of 0.860671 s Par no = 53538 Density = 8.815587e-04 \n",
      "Evaluated =  97 % Time Left = 0.025905 s of 0.863784 s Par no = 54890 Density = 2.194311e-02 \n",
      "Evaluated =  98 % Time Left = 0.017521 s of 0.876478 s Par no = 10360 Density = 3.505276e-03 \n",
      "Evaluated =  99 % Time Left = 0.008956 s of 0.896463 s Par no = 62553 Density = 2.918437e-02 \n",
      "Evaluated = 100 % Time Left = -0.000009 s of 0.910781 s Par no = 98678 Density = 2.092983e-05 \n",
      "\n",
      "Total Smoothing Time = 0.910877 s\n",
      "1 3\n",
      "\n",
      "writing snapshot file ....\n",
      "to_enbid_d3n64.est\n",
      "done with snapshot.\n",
      "\n",
      "Total Time = 1.34495 s \n",
      "\n",
      "/home/athob/miniconda3/envs/main/.cache/Galaxia_ananke/bin/galaxia -r --hdim=6 --nfile=sim /home/athob/Software/repos/py-ananke/jupyter/survey_params\n",
      "CODEDATAPATH=/home/athob/miniconda3/envs/main/.cache/Galaxia_ananke/GalaxiaData/\n",
      "Reading Parameter file-             /home/athob/Software/repos/py-ananke/jupyter/survey_params\n",
      "--------------------------------------------------------\n",
      "outputFile               survey                  \n",
      "outputDir                /home/athob/Software/repos/py-ananke/jupyter\n",
      "photoCateg               padova                  \n",
      "photoSys                 GAIA                    \n",
      "magcolorNames            Gmag,G_BPmag-G_RPmag    \n",
      "appMagLimits[0]          -1000                   \n",
      "appMagLimits[1]          30                      \n",
      "absMagLimits[0]          -1000                   \n",
      "absMagLimits[1]          1000                    \n",
      "colorLimits[0]           -1000                   \n",
      "colorLimits[1]           1000                    \n",
      "geometryOption           0                       \n",
      "surveyArea               207.455                 \n",
      "fSample                  0.01                    \n",
      "popID                    10                      \n",
      "warpFlareOn              0                       \n",
      "longitude                76.273                  \n",
      "latitude                 13.4725                 \n",
      "starType                 0                       \n",
      "photoError               0                       \n",
      "seed                     17052                   \n",
      "r_max                    10000                   \n",
      "r_min                    0                       \n",
      "nres                     64                      \n",
      "nstart                   0                       \n",
      "rSun[0]                  4873.920645095638       \n",
      "rSun[1]                  -855.0925212732859      \n",
      "rSun[2]                  716.7386729994488       \n",
      "vSun[0]                  0.0                     \n",
      "vSun[1]                  0.0                     \n",
      "vSun[2]                  0.0                     \n",
      "--------------------------------------------------------\n",
      "Reading Halo Sat File=/home/athob/miniconda3/envs/main/.cache/Galaxia_ananke/GalaxiaData/nbody1/filenames/sim.txt\n",
      "Halo 1 Sat 0 fname nbody1/sim/sim.ebf 0\n",
      "nbody1/sim/sim.ebf 0\n",
      "No of Satellites   =1\n",
      "Generating catalog. Max allowed stars is 18446744073709551615\n",
      "Source numbering will start at 0\n",
      "Writing to /home/athob/Software/repos/py-ananke/jupyter/survey.sim.ebf\n",
      "Using geometry:                     All Sky\n",
      "setting center to 4873.92 -855.093 716.739 0 0 0\n",
      "Reading Isochrones from dir-        /home/athob/miniconda3/envs/main/.cache/Galaxia_ananke/GalaxiaData/Isochrones/padova/GAIA\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "nbody1/sim/\n",
      "\n",
      "\t1\t1\n",
      "\n",
      "path nbody1/sim/ sats 1 x[0] 1\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Isochrone Grid Size:                (Age bins=71,Feh bins=34,Alpha bins=1)\n",
      "Time Isochrone Reading              3.0945      \n",
      "------------------------------\n",
      "nbody1/sim/sim.ebf  Sat No=0\n",
      "Particles=100000\n",
      "Satellite Info\n",
      "Particles=100000 Mass=5.49989e+08 0.379706\n",
      "Total Stars=48429 accepted=14062 rejected=34367\n",
      "-----------Done---------------\n",
      "Total stars written                 14062                   \n",
      "30\n",
      "File written-                       /home/athob/Software/repos/py-ananke/jupyter/survey.sim.ebf\n",
      "Calculating magnitudes in GAIA system................\n",
      "initializing isochrone data\n",
      "interpolating on isochrone tables\n",
      "Appending spherical coordinates................\n",
      "Total Time=                         3.21774     \n",
      "Exported the following quantities to /home/athob/Software/repos/py-ananke/jupyter/survey.sim.h5\n",
      "['age', 'alpha', 'dec', 'dmod', 'feh', 'gaia_g_bpmag', 'gaia_g_rpmag', 'gaia_gmag', 'glat', 'glon', 'grav', 'lum', 'mact', 'mtip', 'parentid', 'partid', 'px', 'py', 'pz', 'ra', 'rad', 'smass', 'teff', 'vx', 'vy', 'vz']\n"
     ]
    }
   ],
   "source": [
    "survey = ananke.run()"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "1664586c",
   "metadata": {},
   "source": [
    "The output is saved as a `vaex` dataframe"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "42e393f6",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>.vaex-description pre {\n",
       "          max-width : 450px;\n",
       "          white-space : nowrap;\n",
       "          overflow : hidden;\n",
       "          text-overflow: ellipsis;\n",
       "        }\n",
       "\n",
       "        .vex-description pre:hover {\n",
       "          max-width : initial;\n",
       "          white-space: pre;\n",
       "        }</style>"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<thead>\n",
       "<tr><th>#                                 </th><th>age              </th><th>alpha                </th><th>dec      </th><th>dmod              </th><th>feh                   </th><th>gaia_g_bpmag  </th><th>gaia_g_rpmag  </th><th>gaia_gmag  </th><th>glat      </th><th>glon     </th><th>grav     </th><th>lum      </th><th>mact     </th><th>mtip     </th><th>parentid  </th><th>partid  </th><th>px                </th><th>py               </th><th>pz                </th><th>ra       </th><th>rad               </th><th>smass             </th><th>teff     </th><th>vx                 </th><th>vy                 </th><th>vz                 </th><th>log10_NH          </th><th>E(B-V)               </th><th>A_0                 </th><th>A_gaia_gmag          </th><th>A_gaia_g_bpmag      </th><th>A_gaia_g_rpmag      </th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "<tr><td><i style='opacity: 0.6'>0</i>     </td><td>9.308772087097168</td><td>-1.2332874536514282  </td><td>32.418007</td><td>28.484609519084557</td><td>-0.0037263932172209024</td><td>1.0253423     </td><td>-0.20729075   </td><td>0.4835399  </td><td>-8.125952 </td><td>169.96152</td><td>2.554411 </td><td>1.7408993</td><td>1.686    </td><td>1.6848724</td><td>4         </td><td>0       </td><td>-4851.089514114193</td><td>858.7427717081614</td><td>-703.4227760115338</td><td>71.96125 </td><td>4976.475883968641 </td><td>1.6863790627254682</td><td>3.6696734</td><td>-9.245811462402344 </td><td>20.93867301940918  </td><td>54.28696823120117  </td><td>19.64856970249834 </td><td>0.0017808596494772331</td><td>0.005520664913379423</td><td>0.0043117042588483685</td><td>0.005664448771016373</td><td>0.003214528588895797</td></tr>\n",
       "<tr><td><i style='opacity: 0.6'>1</i>     </td><td>8.795965194702148</td><td>-0.14946995675563812 </td><td>31.928228</td><td>28.462504480824826</td><td>-0.9277363419532776   </td><td>1.4407223     </td><td>1.5200416     </td><td>1.4697223  </td><td>-8.3500395</td><td>170.41609</td><td>4.3145313</td><td>1.5043409</td><td>1.8545854</td><td>2.1818192</td><td>8         </td><td>0       </td><td>-4805.82800898724 </td><td>811.4615512720652</td><td>-715.3659174719525</td><td>72.09454 </td><td>4926.073598418888 </td><td>1.8555609252927683</td><td>4.040267 </td><td>-23.906299591064453</td><td>-23.05091667175293 </td><td>-119.33078002929688</td><td>21.72862339032904 </td><td>0.21413289160411542  </td><td>0.6638119639727578  </td><td>0.6376884618807591   </td><td>0.7588439124834473  </td><td>0.40263227038951077 </td></tr>\n",
       "<tr><td><i style='opacity: 0.6'>2</i>     </td><td>8.795965194702148</td><td>-0.14946995675563812 </td><td>31.747633</td><td>28.461467119434154</td><td>-0.9277363419532776   </td><td>1.3381635     </td><td>1.4195217     </td><td>1.3675218  </td><td>-8.421482 </td><td>170.59296</td><td>4.285985 </td><td>1.5501645</td><td>1.8897315</td><td>2.1818192</td><td>8         </td><td>1       </td><td>-4805.130089160716</td><td>796.0965347754252</td><td>-721.0979676465176</td><td>72.1601  </td><td>4923.720863877546 </td><td>1.8902745034999806</td><td>4.042415 </td><td>-20.556483367955813</td><td>-11.19114403311702 </td><td>-109.55960273134058</td><td>21.8529690185789  </td><td>0.28512087155530075  </td><td>0.8838747018214324  </td><td>0.8411038085359314   </td><td>1.0061914597172579  </td><td>0.5346279194339567  </td></tr>\n",
       "<tr><td><i style='opacity: 0.6'>3</i>     </td><td>9.381908416748047</td><td>0.1484881341457367   </td><td>32.212658</td><td>28.51779980536087 </td><td>-0.3162284791469574   </td><td>0.67582935    </td><td>-0.41144708   </td><td>0.21355292 </td><td>-8.292467 </td><td>170.09169</td><td>2.541    </td><td>1.8177235</td><td>1.5157236</td><td>1.50785  </td><td>21        </td><td>0       </td><td>-4925.709123946956</td><td>860.415305697114 </td><td>-728.7921006956768</td><td>71.90629 </td><td>5053.1240631958135</td><td>1.5248914165021865</td><td>3.6970446</td><td>-11.731164932250977</td><td>-8.347715377807617 </td><td>-1.6173608303070068</td><td>21.6300966176637  </td><td>0.17066977226862987  </td><td>0.5290762940327526  </td><td>0.415480433490774    </td><td>0.5451291338087634  </td><td>0.30839324719142563 </td></tr>\n",
       "<tr><td><i style='opacity: 0.6'>4</i>     </td><td>9.095987319946289</td><td>-0.017641598358750343</td><td>33.066193</td><td>28.478124799440803</td><td>-0.6019725203514099   </td><td>1.559597      </td><td>1.4334579     </td><td>1.5148839  </td><td>-7.8933043</td><td>169.30615</td><td>4.0220127</td><td>1.2962763</td><td>1.6351371</td><td>1.7589755</td><td>28        </td><td>0       </td><td>-4829.273081161312</td><td>911.9691959253855</td><td>-681.3752849379254</td><td>71.68203 </td><td>4961.636684166859 </td><td>1.6357240475969435</td><td>3.9286296</td><td>-61.722835540771484</td><td>17.274423599243164 </td><td>-79.17394256591797 </td><td>21.760445674655806</td><td>0.2304123036393029   </td><td>0.714278141281839   </td><td>0.6607178790936152   </td><td>0.8031938276712772  </td><td>0.4305090931187769  </td></tr>\n",
       "<tr><td>...                               </td><td>...              </td><td>...                  </td><td>...      </td><td>...               </td><td>...                   </td><td>...           </td><td>...           </td><td>...        </td><td>...       </td><td>...      </td><td>...      </td><td>...      </td><td>...      </td><td>...      </td><td>...       </td><td>...     </td><td>...               </td><td>...              </td><td>...               </td><td>...      </td><td>...               </td><td>...               </td><td>...      </td><td>...                </td><td>...                </td><td>...                </td><td>...               </td><td>...                  </td><td>...                 </td><td>...                  </td><td>...                 </td><td>...                 </td></tr>\n",
       "<tr><td><i style='opacity: 0.6'>14,057</i></td><td>8.878317832946777</td><td>-0.02291102707386017 </td><td>32.20407 </td><td>28.49246106311801 </td><td>-1.0532594919204712   </td><td>1.4030721     </td><td>1.4818171     </td><td>1.4319022  </td><td>-8.007031 </td><td>170.3397 </td><td>4.2904644</td><td>1.5196012</td><td>1.8255359</td><td>2.06986  </td><td>99961     </td><td>0       </td><td>-4875.678212501339</td><td>829.9435230127879</td><td>-695.7071121016216</td><td>72.351135</td><td>4994.502244376301 </td><td>1.826051196048781 </td><td>4.0398893</td><td>-94.0284652709961  </td><td>26.59738540649414  </td><td>-36.70679473876953 </td><td>21.999959794531037</td><td>0.3999629711086328   </td><td>1.2398852104367617  </td><td>1.1608994747421792   </td><td>1.4014781998642274  </td><td>0.7464984366283208  </td></tr>\n",
       "<tr><td><i style='opacity: 0.6'>14,058</i></td><td>8.878317832946777</td><td>-0.02291102707386017 </td><td>32.17042 </td><td>28.493987005649632</td><td>-1.0532594919204712   </td><td>1.4676036     </td><td>1.5443108     </td><td>1.495075   </td><td>-8.012118 </td><td>170.37938</td><td>4.310163 </td><td>1.4924301</td><td>1.8058368</td><td>2.06986  </td><td>99961     </td><td>1       </td><td>-4879.618764690416</td><td>827.137164014496 </td><td>-696.6355146700624</td><td>72.375854</td><td>4998.013226984538 </td><td>1.8067118855478792</td><td>4.0390067</td><td>-100.84668284186372</td><td>31.843649395250658 </td><td>-39.45628194423663 </td><td>22.08752176398854 </td><td>0.48930736787882817  </td><td>1.5168528404243673  </td><td>1.4025157137536441   </td><td>1.705200374928969   </td><td>0.9099650029954219  </td></tr>\n",
       "<tr><td><i style='opacity: 0.6'>14,059</i></td><td>9.533735275268555</td><td>0.4945138394832611   </td><td>31.501265</td><td>28.469588757430913</td><td>-0.8962724804878235   </td><td>-0.5359952    </td><td>-1.7004024    </td><td>-1.0372192 </td><td>-8.640716 </td><td>170.73337</td><td>1.855697 </td><td>2.3409772</td><td>1.233    </td><td>1.2443578</td><td>99972     </td><td>0       </td><td>-4822.310854293636</td><td>786.8050914637156</td><td>-742.5012519636822</td><td>72.06434 </td><td>4942.1708121592765</td><td>1.2627260865250525</td><td>3.6789184</td><td>-86.41449737548828 </td><td>-1.34444260597229  </td><td>72.12390899658203  </td><td>21.802687379019993</td><td>0.2539495048316513   </td><td>0.787243464978119   </td><td>0.6034074981093309   </td><td>0.8019255460851893  </td><td>0.4556852491475808  </td></tr>\n",
       "<tr><td><i style='opacity: 0.6'>14,060</i></td><td>9.103982925415039</td><td>0.24086332321166992  </td><td>32.159424</td><td>28.49281689796909 </td><td>-1.0517905950546265   </td><td>1.209412      </td><td>1.1521823     </td><td>1.1902008  </td><td>-8.776904 </td><td>169.75446</td><td>3.954495 </td><td>1.4586608</td><td>1.6370276</td><td>1.678095 </td><td>99978     </td><td>0       </td><td>-4858.104566909603</td><td>878.1032882563181</td><td>-762.2230861464215</td><td>71.21421 </td><td>4995.320750555371 </td><td>1.6375016392882022</td><td>3.9523072</td><td>-22.41616439819336 </td><td>-24.103288650512695</td><td>-21.061344146728516</td><td>22.418101990482803</td><td>1.0475191756565996   </td><td>3.2473094445354587  </td><td>2.726742013004422    </td><td>3.5040009480813974  </td><td>1.8987928076857015  </td></tr>\n",
       "<tr><td><i style='opacity: 0.6'>14,061</i></td><td>9.475979804992676</td><td>0.5997751951217651   </td><td>31.852282</td><td>28.503839968745147</td><td>-0.9664602875709534   </td><td>0.13237825    </td><td>-0.83745486   </td><td>-0.26653275</td><td>-8.732972 </td><td>170.19609</td><td>2.3693337</td><td>2.0007885</td><td>1.283    </td><td>1.2862073</td><td>99994     </td><td>0       </td><td>-4890.06365237408 </td><td>845.0093150965281</td><td>-762.2973659315777</td><td>71.583595</td><td>5020.74302678162  </td><td>1.3062057771381834</td><td>3.7179322</td><td>27.791748046875    </td><td>113.67877197265625 </td><td>-6.7907280921936035</td><td>22.30888045775735 </td><td>0.8145925785740111   </td><td>2.5252369935794348  </td><td>1.8733800066502095   </td><td>2.545348166726454   </td><td>1.4444557989660336  </td></tr>\n",
       "</tbody>\n",
       "</table>"
      ],
      "text/plain": [
       "#       age                alpha                  dec        dmod                feh                     gaia_g_bpmag    gaia_g_rpmag    gaia_gmag    glat        glon       grav       lum        mact       mtip       parentid    partid    px                  py                 pz                  ra         rad                 smass               teff       vx                   vy                   vz                   log10_NH            E(B-V)                 A_0                   A_gaia_gmag            A_gaia_g_bpmag        A_gaia_g_rpmag\n",
       "0       9.308772087097168  -1.2332874536514282    32.418007  28.484609519084557  -0.0037263932172209024  1.0253423       -0.20729075     0.4835399    -8.125952   169.96152  2.554411   1.7408993  1.686      1.6848724  4           0         -4851.089514114193  858.7427717081614  -703.4227760115338  71.96125   4976.475883968641   1.6863790627254682  3.6696734  -9.245811462402344   20.93867301940918    54.28696823120117    19.64856970249834   0.0017808596494772331  0.005520664913379423  0.0043117042588483685  0.005664448771016373  0.003214528588895797\n",
       "1       8.795965194702148  -0.14946995675563812   31.928228  28.462504480824826  -0.9277363419532776     1.4407223       1.5200416       1.4697223    -8.3500395  170.41609  4.3145313  1.5043409  1.8545854  2.1818192  8           0         -4805.82800898724   811.4615512720652  -715.3659174719525  72.09454   4926.073598418888   1.8555609252927683  4.040267   -23.906299591064453  -23.05091667175293   -119.33078002929688  21.72862339032904   0.21413289160411542    0.6638119639727578    0.6376884618807591     0.7588439124834473    0.40263227038951077\n",
       "2       8.795965194702148  -0.14946995675563812   31.747633  28.461467119434154  -0.9277363419532776     1.3381635       1.4195217       1.3675218    -8.421482   170.59296  4.285985   1.5501645  1.8897315  2.1818192  8           1         -4805.130089160716  796.0965347754252  -721.0979676465176  72.1601    4923.720863877546   1.8902745034999806  4.042415   -20.556483367955813  -11.19114403311702   -109.55960273134058  21.8529690185789    0.28512087155530075    0.8838747018214324    0.8411038085359314     1.0061914597172579    0.5346279194339567\n",
       "3       9.381908416748047  0.1484881341457367     32.212658  28.51779980536087   -0.3162284791469574     0.67582935      -0.41144708     0.21355292   -8.292467   170.09169  2.541      1.8177235  1.5157236  1.50785    21          0         -4925.709123946956  860.415305697114   -728.7921006956768  71.90629   5053.1240631958135  1.5248914165021865  3.6970446  -11.731164932250977  -8.347715377807617   -1.6173608303070068  21.6300966176637    0.17066977226862987    0.5290762940327526    0.415480433490774      0.5451291338087634    0.30839324719142563\n",
       "4       9.095987319946289  -0.017641598358750343  33.066193  28.478124799440803  -0.6019725203514099     1.559597        1.4334579       1.5148839    -7.8933043  169.30615  4.0220127  1.2962763  1.6351371  1.7589755  28          0         -4829.273081161312  911.9691959253855  -681.3752849379254  71.68203   4961.636684166859   1.6357240475969435  3.9286296  -61.722835540771484  17.274423599243164   -79.17394256591797   21.760445674655806  0.2304123036393029     0.714278141281839     0.6607178790936152     0.8031938276712772    0.4305090931187769\n",
       "...     ...                ...                    ...        ...                 ...                     ...             ...             ...          ...         ...        ...        ...        ...        ...        ...         ...       ...                 ...                ...                 ...        ...                 ...                 ...        ...                  ...                  ...                  ...                 ...                    ...                   ...                    ...                   ...\n",
       "14,057  8.878317832946777  -0.02291102707386017   32.20407   28.49246106311801   -1.0532594919204712     1.4030721       1.4818171       1.4319022    -8.007031   170.3397   4.2904644  1.5196012  1.8255359  2.06986    99961       0         -4875.678212501339  829.9435230127879  -695.7071121016216  72.351135  4994.502244376301   1.826051196048781   4.0398893  -94.0284652709961    26.59738540649414    -36.70679473876953   21.999959794531037  0.3999629711086328     1.2398852104367617    1.1608994747421792     1.4014781998642274    0.7464984366283208\n",
       "14,058  8.878317832946777  -0.02291102707386017   32.17042   28.493987005649632  -1.0532594919204712     1.4676036       1.5443108       1.495075     -8.012118   170.37938  4.310163   1.4924301  1.8058368  2.06986    99961       1         -4879.618764690416  827.137164014496   -696.6355146700624  72.375854  4998.013226984538   1.8067118855478792  4.0390067  -100.84668284186372  31.843649395250658   -39.45628194423663   22.08752176398854   0.48930736787882817    1.5168528404243673    1.4025157137536441     1.705200374928969     0.9099650029954219\n",
       "14,059  9.533735275268555  0.4945138394832611     31.501265  28.469588757430913  -0.8962724804878235     -0.5359952      -1.7004024      -1.0372192   -8.640716   170.73337  1.855697   2.3409772  1.233      1.2443578  99972       0         -4822.310854293636  786.8050914637156  -742.5012519636822  72.06434   4942.1708121592765  1.2627260865250525  3.6789184  -86.41449737548828   -1.34444260597229    72.12390899658203    21.802687379019993  0.2539495048316513     0.787243464978119     0.6034074981093309     0.8019255460851893    0.4556852491475808\n",
       "14,060  9.103982925415039  0.24086332321166992    32.159424  28.49281689796909   -1.0517905950546265     1.209412        1.1521823       1.1902008    -8.776904   169.75446  3.954495   1.4586608  1.6370276  1.678095   99978       0         -4858.104566909603  878.1032882563181  -762.2230861464215  71.21421   4995.320750555371   1.6375016392882022  3.9523072  -22.41616439819336   -24.103288650512695  -21.061344146728516  22.418101990482803  1.0475191756565996     3.2473094445354587    2.726742013004422      3.5040009480813974    1.8987928076857015\n",
       "14,061  9.475979804992676  0.5997751951217651     31.852282  28.503839968745147  -0.9664602875709534     0.13237825      -0.83745486     -0.26653275  -8.732972   170.19609  2.3693337  2.0007885  1.283      1.2862073  99994       0         -4890.06365237408   845.0093150965281  -762.2973659315777  71.583595  5020.74302678162    1.3062057771381834  3.7179322  27.791748046875      113.67877197265625   -6.7907280921936035  22.30888045775735   0.8145925785740111     2.5252369935794348    1.8733800066502095     2.545348166726454     1.4444557989660336"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "survey"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "base",
   "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.9"
  },
  "vscode": {
   "interpreter": {
    "hash": "93aa337d3be5595eae1049f286a1f242db5cd4d40fd613efbe29fb4c796799bf"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}