From b9592bd4052dcce80d3eb0d5f85e2c4315fb64c5 Mon Sep 17 00:00:00 2001 From: TristanCantatGaudin Date: Mon, 6 May 2024 11:11:42 +0200 Subject: [PATCH] Add files via upload --- docs/notebooks/SubsampleSF_Tutorial.ipynb | 40 ++++++------- docs/notebooks/dr3-unresolved-binaries.ipynb | 62 ++++++++++---------- docs/notebooks/rave_ratios.ipynb | 28 ++++----- 3 files changed, 65 insertions(+), 65 deletions(-) diff --git a/docs/notebooks/SubsampleSF_Tutorial.ipynb b/docs/notebooks/SubsampleSF_Tutorial.ipynb index d9a0180..3a38ecf 100644 --- a/docs/notebooks/SubsampleSF_Tutorial.ipynb +++ b/docs/notebooks/SubsampleSF_Tutorial.ipynb @@ -2,10 +2,10 @@ "cells": [ { "cell_type": "markdown", - "id": "58d9e9e6", + "id": "08d746b7", "metadata": {}, "source": [ - "# Selection function for a subset of the Gaia catalogue\n", + "# 📊 Selection function for a subset of the Gaia catalogue\n", "\n", "The GaiaUnlimited module **subsample** provides a simple way of counting sources satisfying user-defined criteria, out of the whole Gaia catalogue. The result is given per bin of colour and magnitude, and per healpix region on the sky. The coulour and magnitude binning can be chosen by the user, as well as the spatial resolution (order of the healpix tessellation).\n", "\n", @@ -17,7 +17,7 @@ { "cell_type": "code", "execution_count": 1, - "id": "5d307655", + "id": "6fc778ba", "metadata": {}, "outputs": [], "source": [ @@ -35,7 +35,7 @@ { "cell_type": "code", "execution_count": 2, - "id": "9ea2d42b", + "id": "a05f7ef2", "metadata": {}, "outputs": [], "source": [ @@ -48,7 +48,7 @@ }, { "cell_type": "markdown", - "id": "583b0c9e", + "id": "54bb32f6", "metadata": {}, "source": [ "Define dependencies of the selection function:\n", @@ -63,7 +63,7 @@ { "cell_type": "code", "execution_count": 3, - "id": "c9576b2d", + "id": "f439a7ab", "metadata": {}, "outputs": [], "source": [ @@ -72,7 +72,7 @@ }, { "cell_type": "markdown", - "id": "c831785e", + "id": "8759897c", "metadata": {}, "source": [ "Launch the query to the Gaia DR3 catalogue to determine which fraction of sources have a valid radial velocity (**radial_velocity is not null**). Depending on the exact query, this operation can take ~40 minutes.\n", @@ -85,7 +85,7 @@ { "cell_type": "code", "execution_count": 4, - "id": "012878a9", + "id": "80b09818", "metadata": {}, "outputs": [ { @@ -105,7 +105,7 @@ }, { "cell_type": "markdown", - "id": "43a3d105", + "id": "afd1f19a", "metadata": {}, "source": [ "Now we want to visualise the results for the entire sky, so we generate a list of coordinates of the centers of all healpix regions of order 4.\n", @@ -118,7 +118,7 @@ { "cell_type": "code", "execution_count": 5, - "id": "50d49621", + "id": "558ff246", "metadata": {}, "outputs": [ { @@ -168,7 +168,7 @@ }, { "cell_type": "markdown", - "id": "7a50ac0e", + "id": "03ea0986", "metadata": {}, "source": [ "## Selection function for a user-updated table (requires authenticated query)\n", @@ -181,7 +181,7 @@ { "cell_type": "code", "execution_count": 9, - "id": "d80a54f2", + "id": "655de081", "metadata": {}, "outputs": [ { @@ -203,7 +203,7 @@ }, { "cell_type": "markdown", - "id": "33f1b328", + "id": "488b3867", "metadata": {}, "source": [ "In this example we construct the selection function for the Gaia-Enceladus Sausage sample. In Castro-Ginard et al. (2023, Sect. 4) we reproduced the cuts from Helmi et al. (2018) and uploaded the table to the Gaia archive as:\n", @@ -218,7 +218,7 @@ { "cell_type": "code", "execution_count": 10, - "id": "9350b3d3", + "id": "18495cc9", "metadata": {}, "outputs": [], "source": [ @@ -230,7 +230,7 @@ }, { "cell_type": "markdown", - "id": "8cf4cd4b", + "id": "0ca17b9d", "metadata": {}, "source": [ "which takes about 45 minutes to complete. The result can then be queried and displayed:" @@ -239,7 +239,7 @@ { "cell_type": "code", "execution_count": 11, - "id": "964f85ac", + "id": "e0f0c574", "metadata": {}, "outputs": [ { @@ -290,7 +290,7 @@ { "cell_type": "code", "execution_count": null, - "id": "a81d187e", + "id": "78310f53", "metadata": {}, "outputs": [], "source": [] @@ -298,7 +298,7 @@ { "cell_type": "code", "execution_count": null, - "id": "8088a229", + "id": "3f246b2f", "metadata": {}, "outputs": [], "source": [] @@ -306,7 +306,7 @@ ], "metadata": { "kernelspec": { - "display_name": "Python 3 (ipykernel)", + "display_name": "Python 3", "language": "python", "name": "python3" }, @@ -320,7 +320,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.9.6" + "version": "3.7.6" } }, "nbformat": 4, diff --git a/docs/notebooks/dr3-unresolved-binaries.ipynb b/docs/notebooks/dr3-unresolved-binaries.ipynb index cd77f08..5f596a2 100644 --- a/docs/notebooks/dr3-unresolved-binaries.ipynb +++ b/docs/notebooks/dr3-unresolved-binaries.ipynb @@ -2,17 +2,17 @@ "cells": [ { "cell_type": "markdown", - "id": "3a213997", + "id": "880e039b", "metadata": {}, "source": [ - "# 💫 Gaia DR3 detectability of unresolved binary systems\n", + "# Gaia DR3 detectability of unresolved binary systems\n", "\n", "updated 23 Apr. 2024" ] }, { "cell_type": "markdown", - "id": "7e92c986", + "id": "742e1aff", "metadata": {}, "source": [ "In this notebook, we discuss how to estimate the detectability of unresolved binary systems based on RUWE. \n", @@ -24,7 +24,7 @@ }, { "cell_type": "markdown", - "id": "4332eee3", + "id": "a465afe3", "metadata": {}, "source": [ "### Sky-varying RUWE threshold" @@ -32,7 +32,7 @@ }, { "cell_type": "markdown", - "id": "6be0eb6c", + "id": "d440f844", "metadata": {}, "source": [ "One of the main results in [Castro-Ginard et al. (2024)](https://arxiv.org/abs/2404.14127) is Fig. 3, which shows a sky-varying RUWE thresdhold above which soucers show potential signs of binarity. \n", @@ -43,7 +43,7 @@ { "cell_type": "code", "execution_count": 1, - "id": "50dd2942", + "id": "77c623e1", "metadata": {}, "outputs": [ { @@ -62,7 +62,7 @@ }, { "cell_type": "markdown", - "id": "0ad1fde3", + "id": "3837dca0", "metadata": {}, "source": [ "This object provides a pre-computed table at HEALPix level 5, which includes for each pixel the RUWE thresholds together with the Gaia scanning law details needed to compute them and the probability to select sources with a RUWE higher than the threshold (estimated as in [Castro-Ginard et al., 2023](https://ui.adsabs.harvard.edu/abs/2023A%26A...677A..37C/abstract)).\n", @@ -73,7 +73,7 @@ { "cell_type": "code", "execution_count": 5, - "id": "61bfc12c", + "id": "e59366bd", "metadata": {}, "outputs": [ { @@ -98,7 +98,7 @@ }, { "cell_type": "markdown", - "id": "e400db45", + "id": "d9587bdd", "metadata": {}, "source": [ "Then, we query the values for the RUWE threshold with:" @@ -107,7 +107,7 @@ { "cell_type": "code", "execution_count": 6, - "id": "efcf413c", + "id": "f7f89852", "metadata": {}, "outputs": [ { @@ -128,7 +128,7 @@ }, { "cell_type": "markdown", - "id": "0c13ad54", + "id": "c0ad55ae", "metadata": {}, "source": [ "The argument \"crowding = True\" gives us the contribution to a increased RUWE from sources other than binarity, gauged directly from Gaia data. Setting \"crowding = False\" will provide a RUWE threshold where values above are only increased by binarity. The default and recommended setting for \"crowding\" is True.\n", @@ -139,7 +139,7 @@ { "cell_type": "code", "execution_count": 14, - "id": "e475dfc6", + "id": "968b3918", "metadata": {}, "outputs": [ { @@ -175,7 +175,7 @@ }, { "cell_type": "markdown", - "id": "ca4475c5", + "id": "adb0d7d0", "metadata": {}, "source": [ "The Gaia scanning law details to estimate the above values are also accessible through the selection function object we created. This will provide the observation times, scanning angles or the along-scan parallax factors for each CCD observation of a certain pixel (usually 9 CCD observations per source transit). The number of observations can be derived from the previous. These data are collected through the Gaia Observation Forecast Tool ([GOST](https://gaia.esac.esa.int/gost/))." @@ -184,7 +184,7 @@ { "cell_type": "code", "execution_count": 17, - "id": "386f4928", + "id": "246b7802", "metadata": {}, "outputs": [ { @@ -233,7 +233,7 @@ }, { "cell_type": "markdown", - "id": "9ad2401d", + "id": "43fbabe0", "metadata": {}, "source": [ "### Selection function of potential binary systems based on RUWE" @@ -241,7 +241,7 @@ }, { "cell_type": "markdown", - "id": "2560b2a2", + "id": "7bbbd66d", "metadata": {}, "source": [ "Based on this sky-varyin RUWE threshold we can estimate the probability to select sources with RUWE higher than the threshold, which can be seen as potential binary systems (note that there are other contributions that will lead to an increased RUWE which are not taken into account).\n", @@ -252,7 +252,7 @@ { "cell_type": "code", "execution_count": 31, - "id": "99e89a6d", + "id": "b43b2d0c", "metadata": {}, "outputs": [ { @@ -308,7 +308,7 @@ }, { "cell_type": "markdown", - "id": "530f90f9", + "id": "39e993d5", "metadata": {}, "source": [ "### Forward model to estimate the RUWE for a simulated source" @@ -316,7 +316,7 @@ }, { "cell_type": "markdown", - "id": "8e7c4efc", + "id": "53b88a7f", "metadata": {}, "source": [ "This functionality shows how to estimate a RUWE value for a simulated source observed by Gaia, being the source either a single source or a binary." @@ -325,7 +325,7 @@ { "cell_type": "code", "execution_count": 1, - "id": "1cf862c8", + "id": "936b8084", "metadata": {}, "outputs": [], "source": [ @@ -334,7 +334,7 @@ }, { "cell_type": "markdown", - "id": "cf9e4513", + "id": "c04eed1a", "metadata": {}, "source": [ "We simulate both a single source and a binary at a given location in the sky, providing $\\alpha^*$ and $\\delta$. \n", @@ -347,7 +347,7 @@ { "cell_type": "code", "execution_count": 7, - "id": "02b0f0d1", + "id": "a0980d88", "metadata": {}, "outputs": [ { @@ -376,7 +376,7 @@ }, { "cell_type": "markdown", - "id": "2971d54a", + "id": "1d50d19a", "metadata": {}, "source": [ "We need the G magnitude of the source to compure its along-scan astrometric error and the parallax of the source.\n", @@ -387,7 +387,7 @@ { "cell_type": "code", "execution_count": 9, - "id": "8f980d07", + "id": "e387779e", "metadata": {}, "outputs": [], "source": [ @@ -414,7 +414,7 @@ }, { "cell_type": "markdown", - "id": "389cbb52", + "id": "77318641", "metadata": {}, "source": [ "We now estimate the ruwes using the following function: " @@ -423,7 +423,7 @@ { "cell_type": "code", "execution_count": 11, - "id": "91c7f3b3", + "id": "ba41834d", "metadata": {}, "outputs": [ { @@ -445,7 +445,7 @@ }, { "cell_type": "markdown", - "id": "a6c90347", + "id": "7fce77c8", "metadata": {}, "source": [ "The signs of binarity on RUWE may be negligible depending on the properties of the binary system. So, the RUWE values for both the single and binary source may be very similar. \n", @@ -458,7 +458,7 @@ { "cell_type": "code", "execution_count": 24, - "id": "416dd665", + "id": "50e33335", "metadata": {}, "outputs": [], "source": [ @@ -478,7 +478,7 @@ { "cell_type": "code", "execution_count": 25, - "id": "b533e768", + "id": "2de61e9a", "metadata": {}, "outputs": [ { @@ -501,7 +501,7 @@ { "cell_type": "code", "execution_count": 38, - "id": "2d6ce5e8", + "id": "3eb5009e", "metadata": {}, "outputs": [ { @@ -531,7 +531,7 @@ { "cell_type": "code", "execution_count": null, - "id": "feb7d92b", + "id": "8b36b2a1", "metadata": {}, "outputs": [], "source": [] diff --git a/docs/notebooks/rave_ratios.ipynb b/docs/notebooks/rave_ratios.ipynb index c907cfb..03306aa 100644 --- a/docs/notebooks/rave_ratios.ipynb +++ b/docs/notebooks/rave_ratios.ipynb @@ -2,10 +2,10 @@ "cells": [ { "cell_type": "markdown", - "id": "7a74c602", + "id": "2b6f8b3e", "metadata": {}, "source": [ - "# RAVE DR6 selection function\n", + "# 💫 RAVE DR6 selection function\n", "\n", "The RAVE (RAdial Velocity Experiment) target selection and associated selection function is described in [Wojno et al. (2017)](https://ui.adsabs.harvard.edu/abs/2017MNRAS.468.3368W/abstract), who construct it as count ratios relative to the 2MASS catalogue. The RAVE target selection was initially performed in the $I$ band (9 to 12). The survey in fact relies on predicted $I$ estimated from the 2MASS $J$ and $K_s$ bands. In some fields ($5^{\\circ} < |b| < 25^{\\circ}$) an additional requirement of $(J-K_s) > 0.5$ was imposed.\n", "\n", @@ -67,7 +67,7 @@ { "cell_type": "code", "execution_count": 1, - "id": "a2805647", + "id": "4060fb3c", "metadata": {}, "outputs": [ { @@ -115,7 +115,7 @@ { "cell_type": "code", "execution_count": 2, - "id": "276c5d24", + "id": "0d0333c9", "metadata": {}, "outputs": [], "source": [ @@ -129,7 +129,7 @@ { "cell_type": "code", "execution_count": 3, - "id": "0b4954e2", + "id": "b71c9b64", "metadata": {}, "outputs": [ { @@ -167,7 +167,7 @@ { "cell_type": "code", "execution_count": 4, - "id": "fab63adf", + "id": "f79eb22b", "metadata": {}, "outputs": [ { @@ -210,7 +210,7 @@ { "cell_type": "code", "execution_count": 5, - "id": "07a5519f", + "id": "467d0963", "metadata": {}, "outputs": [ { @@ -253,7 +253,7 @@ { "cell_type": "code", "execution_count": 6, - "id": "c129edb4", + "id": "02fe4d83", "metadata": {}, "outputs": [ { @@ -296,7 +296,7 @@ { "cell_type": "code", "execution_count": 7, - "id": "f9ab25a2", + "id": "4b202bf1", "metadata": {}, "outputs": [ { @@ -338,7 +338,7 @@ }, { "cell_type": "markdown", - "id": "0c7246e0", + "id": "e603ac92", "metadata": {}, "source": [ "## Subset: RAVE DR6 stars with small radial velocity errors \n", @@ -370,7 +370,7 @@ { "cell_type": "code", "execution_count": 8, - "id": "f7e78c16", + "id": "0fac397d", "metadata": {}, "outputs": [ { @@ -418,7 +418,7 @@ { "cell_type": "code", "execution_count": 9, - "id": "06bc3a25", + "id": "0deb91e7", "metadata": {}, "outputs": [ { @@ -461,7 +461,7 @@ ], "metadata": { "kernelspec": { - "display_name": "Python 3 (ipykernel)", + "display_name": "Python 3", "language": "python", "name": "python3" }, @@ -475,7 +475,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.10.9" + "version": "3.7.6" } }, "nbformat": 4,