GTO Debris Disk updates

JarronL · Jan 27, 2018 · 837f679 · 837f679
1 parent 2a393d1
commit 837f679
Show file tree

Hide file tree

Showing 25 changed files with 3,266 additions and 1,420 deletions.
diff --git a/notebooks/GTO/Debris/HD10647_contrast.pdf b/notebooks/GTO/Debris/HD10647_contrast.pdf
diff --git a/notebooks/GTO/Debris/HD10647_images.pdf b/notebooks/GTO/Debris/HD10647_images.pdf
diff --git a/notebooks/GTO/Debris/HD10647_spectra.pdf b/notebooks/GTO/Debris/HD10647_spectra.pdf
diff --git a/notebooks/GTO/Debris/HD107146_contrast.pdf b/notebooks/GTO/Debris/HD107146_contrast.pdf
diff --git a/notebooks/GTO/Debris/HD107146_images.pdf b/notebooks/GTO/Debris/HD107146_images.pdf
diff --git a/notebooks/GTO/Debris/HD107146_spectra.pdf b/notebooks/GTO/Debris/HD107146_spectra.pdf
diff --git a/notebooks/GTO/Debris/HD111398.vot b/notebooks/GTO/Debris/HD111398.vot
diff --git a/notebooks/GTO/Debris/HD181327_contrast.pdf b/notebooks/GTO/Debris/HD181327_contrast.pdf
diff --git a/notebooks/GTO/Debris/HD181327_images.pdf b/notebooks/GTO/Debris/HD181327_images.pdf
diff --git a/notebooks/GTO/Debris/HD181327_spectra.pdf b/notebooks/GTO/Debris/HD181327_spectra.pdf
diff --git a/notebooks/GTO/Debris/HD31411.vot b/notebooks/GTO/Debris/HD31411.vot
diff --git a/notebooks/GTO/Debris/HD32297_contrast.pdf b/notebooks/GTO/Debris/HD32297_contrast.pdf
diff --git a/notebooks/GTO/Debris/HD32297_images.pdf b/notebooks/GTO/Debris/HD32297_images.pdf
diff --git a/notebooks/GTO/Debris/HD32297_spectra.pdf b/notebooks/GTO/Debris/HD32297_spectra.pdf
diff --git a/notebooks/GTO/Debris/HD56161.vot b/notebooks/GTO/Debris/HD56161.vot
diff --git a/notebooks/GTO/Debris/HD61005_contrast.pdf b/notebooks/GTO/Debris/HD61005_contrast.pdf
diff --git a/notebooks/GTO/Debris/HD61005_images.pdf b/notebooks/GTO/Debris/HD61005_images.pdf
diff --git a/notebooks/GTO/Debris/HD61005_spectra.pdf b/notebooks/GTO/Debris/HD61005_spectra.pdf
diff --git a/notebooks/GTO/Debris/HR7297.vot b/notebooks/GTO/Debris/HR7297.vot
diff --git a/notebooks/GTO/Debris/iotHor.vot b/notebooks/GTO/Debris/iotHor.vot
diff --git a/notebooks/GTO/NRC_GTO_DebrisDisks-Copy1.ipynb b/notebooks/GTO/NRC_GTO_DebrisDisks-Copy1.ipynb
diff --git a/notebooks/GTO/NRC_GTO_DebrisDisks_LW.ipynb b/notebooks/GTO/NRC_GTO_DebrisDisks_LW.ipynb
diff --git a/notebooks/GTO/NRC_GTO_Debris_TargetAcq.ipynb b/notebooks/GTO/NRC_GTO_Debris_TargetAcq.ipynb
@@ -40,7 +40,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": 6,
    "metadata": {
     "collapsed": true
    },
@@ -62,13 +62,13 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 4,
+   "execution_count": 7,
    "metadata": {
     "collapsed": true
    },
    "outputs": [],
    "source": [
-    "bp_k = S.ObsBandpass('k')"
+    "bp_k = pynrc.bp_2mass('k')"
    ]
   },
   {
@@ -80,10 +80,8 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 5,
-   "metadata": {
-    "collapsed": false
-   },
+   "execution_count": 19,
+   "metadata": {},
    "outputs": [
     {
      "name": "stdout",
@@ -92,46 +90,46 @@
       "HD10647\n",
       " Pattern   NGRP NINT   t_int     t_exp     t_acq     SNR      Well     eff   \n",
       "---------- ---- ---- --------- --------- --------- -------- -------- --------\n",
-      "RAPID         9    1      1.64      1.64      1.83     85.7    0.028   63.390\n",
+      "RAPID        65    1      3.26      3.26      3.32    152.0    0.055   83.472\n",
       " Pattern   NGRP NINT   t_int     t_exp     t_acq     SNR      Well     eff   \n",
       "---------- ---- ---- --------- --------- --------- -------- -------- --------\n",
-      "RAPID         9    1      1.64      1.64      1.83     95.3    0.033   70.472\n",
+      "RAPID        65    1      3.26      3.26      3.32    167.1    0.067   91.749\n",
       "\n",
       "HD107146\n",
       " Pattern   NGRP NINT   t_int     t_exp     t_acq     SNR      Well     eff   \n",
       "---------- ---- ---- --------- --------- --------- -------- -------- --------\n",
-      "RAPID        17    1      3.09      3.09      3.28     70.1    0.018   38.714\n",
+      "RAPID        65    1      3.26      3.26      3.32     84.5    0.018   46.393\n",
       " Pattern   NGRP NINT   t_int     t_exp     t_acq     SNR      Well     eff   \n",
       "---------- ---- ---- --------- --------- --------- -------- -------- --------\n",
-      "RAPID        17    1      3.09      3.09      3.28     70.1    0.018   38.714\n",
+      "RAPID        65    1      3.26      3.26      3.32     84.5    0.018   46.393\n",
       "\n",
       "HD181327\n",
       " Pattern   NGRP NINT   t_int     t_exp     t_acq     SNR      Well     eff   \n",
       "---------- ---- ---- --------- --------- --------- -------- -------- --------\n",
-      "RAPID        17    1      3.09      3.09      3.28     56.0    0.012   30.910\n",
+      "RAPID        65    1      3.26      3.26      3.32     68.4    0.013   37.552\n",
       " Pattern   NGRP NINT   t_int     t_exp     t_acq     SNR      Well     eff   \n",
       "---------- ---- ---- --------- --------- --------- -------- -------- --------\n",
-      "RAPID         9    1      1.64      1.64      1.83     54.4    0.013   40.233\n",
+      "RAPID        65    1      3.26      3.26      3.32    103.0    0.027   56.543\n",
       "\n",
       "HD61005\n",
       " Pattern   NGRP NINT   t_int     t_exp     t_acq     SNR      Well     eff   \n",
       "---------- ---- ---- --------- --------- --------- -------- -------- --------\n",
-      "RAPID        33    1      6.01      6.01      6.20     63.8    0.014   25.650\n",
+      "RAPID        65    1      3.26      3.26      3.32     49.5    0.007   27.190\n",
       " Pattern   NGRP NINT   t_int     t_exp     t_acq     SNR      Well     eff   \n",
       "---------- ---- ---- --------- --------- --------- -------- -------- --------\n",
-      "RAPID         9    1      1.64      1.64      1.83     63.7    0.017   47.137\n",
+      "RAPID        65    1      3.26      3.26      3.32    116.3    0.033   63.886\n",
       "\n",
       "HD32297\n",
       " Pattern   NGRP NINT   t_int     t_exp     t_acq     SNR      Well     eff   \n",
       "---------- ---- ---- --------- --------- --------- -------- -------- --------\n",
-      "RAPID        33    1      6.01      6.01      6.20     66.4    0.015   26.678\n",
+      "RAPID        65    1      3.26      3.26      3.32     51.7    0.008   28.415\n",
       "\n"
      ]
     }
    ],
    "source": [
     "nrc = pynrc.NIRCam(filter='F335M', pupil='CIRCLYOT', ND_acq=True,\n",
-    "                   wind_mode='WINDOW', xpix=128, ypix=128)\n",
+    "                   wind_mode='WINDOW', xpix=64, ypix=64)\n",
     "\n",
     "for i in range(len(args_sources)):\n",
     "\n",
@@ -147,7 +145,7 @@
     "    \n",
     "    print(name_sci)\n",
     "\n",
-    "    ng_levels = [9,17,33]\n",
+    "    ng_levels = [65]\n",
     "    for sp in [sp_sci, sp_ref]:\n",
     "        for ng in ng_levels: \n",
     "            pynrc.setup_logging('ERROR', verbose=False)\n",
@@ -163,6 +161,56 @@
     "pynrc.setup_logging('WARNING', verbose=False)"
    ]
   },
+  {
+   "cell_type": "code",
+   "execution_count": 23,
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "HD32297\n",
+      " Pattern   NGRP NINT   t_int     t_exp     t_acq     SNR      Well     eff   \n",
+      "---------- ---- ---- --------- --------- --------- -------- -------- --------\n",
+      "RAPID         5    1      0.25      0.25      0.31    185.6    0.112  335.849\n",
+      " Pattern   NGRP NINT   t_int     t_exp     t_acq     SNR      Well     eff   \n",
+      "---------- ---- ---- --------- --------- --------- -------- -------- --------\n",
+      "RAPID         5    1      0.25      0.25      0.31    332.1    0.338  600.879\n"
+     ]
+    }
+   ],
+   "source": [
+    "# HD 32297 and reference star are too faint\n",
+    "# Turn off ND_acq for faint TA position\n",
+    "nrc = pynrc.NIRCam(filter='F335M', pupil='CIRCLYOT', ND_acq=False,\n",
+    "                   wind_mode='WINDOW', xpix=64, ypix=64)\n",
+    "\n",
+    "i=4\n",
+    "name_sci, dist_sci, fov, age, spt_sci, kmag_sci = args_sources[i]\n",
+    "name_ref, spt_ref, kmag_ref = ref_sources[i]\n",
+    "\n",
+    "# Science target and reference star spectra\n",
+    "sp_sci = pynrc.stellar_spectrum(spt_sci, kmag_sci, 'vegamag', bp_k)\n",
+    "sp_ref = pynrc.stellar_spectrum(spt_ref, kmag_ref, 'vegamag', bp_k)\n",
+    "sp_sci.name = name_sci\n",
+    "sp_ref.name = name_ref\n",
+    "\n",
+    "\n",
+    "print(name_sci)\n",
+    "\n",
+    "ng_levels = [5,9,17,33,65]\n",
+    "for sp in [sp_sci, sp_ref]:\n",
+    "    for ng in ng_levels: \n",
+    "        pynrc.setup_logging('WARNING', verbose=False)\n",
+    "        tbl = nrc.ramp_optimize(sp, snr_goal=100, snr_frac=0.01, well_frac_max=0.5, verbose=False,\n",
+    "                                ng_min=ng, ng_max=ng, nint_max=1, patterns=['RAPID'])\n",
+    "        nrow = len(tbl)\n",
+    "        if nrow>0:\n",
+    "            print(tbl)\n",
+    "            break\n"
+   ]
+  },
   {
    "cell_type": "markdown",
    "metadata": {},
@@ -174,7 +222,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 6,
+   "execution_count": 24,
    "metadata": {
     "collapsed": true
    },
@@ -185,34 +233,34 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 7,
-   "metadata": {
-    "collapsed": false
-   },
+   "execution_count": 33,
+   "metadata": {},
    "outputs": [
     {
      "name": "stdout",
      "output_type": "stream",
      "text": [
       "BRIGHT1\n",
+      "BRIGHT2\n",
+      "DEEP2\n",
+      "DEEP8\n",
+      "MEDIUM2\n",
+      "MEDIUM8\n",
       "RAPID\n",
       "SHALLOW2\n",
+      "SHALLOW4\n",
       " Pattern   NGRP NINT   t_int     t_exp     t_acq     SNR      Well     eff   \n",
       "---------- ---- ---- --------- --------- --------- -------- -------- --------\n",
-      "BRIGHT1       3    3     53.68    161.05    193.27     31.5    1.040    2.263\n",
-      "BRIGHT1       3    4     53.68    214.74    257.69     36.3    1.040    2.263\n",
-      "RAPID         4    4     42.95    171.79    214.74     29.3    0.832    2.001\n",
-      "RAPID         4    5     42.95    214.74    268.42     32.8    0.832    2.001\n",
-      "RAPID         4    6     42.95    257.68    322.11     35.9    0.832    2.001\n"
+      "BRIGHT1       3    3     53.68    161.05    193.27     64.8    1.043    4.657\n",
+      "RAPID         3    3     32.21     96.63    128.85     36.2    0.626    3.185\n"
      ]
     }
    ],
    "source": [
-    "sp_dim = pynrc.stellar_spectrum('flat', 17, 'vegamag', nrc.bandpass)\n",
+    "sp_dim = pynrc.stellar_spectrum('flat', 16, 'vegamag', nrc.bandpass)\n",
     "sp_bright = pynrc.stellar_spectrum('flat', 11, 'vegamag', nrc.bandpass)\n",
-    "nrc = pynrc.NIRCam(filter='F335M', pupil='CIRCLYOT', wind_mode='FULL')\n",
-    "tbl = nrc.ramp_optimize(sp_dim, sp_bright, snr_goal=30, snr_frac=0.1, ng_min=2, ng_max=10, nint_min=3,\n",
-    "                        well_frac_max=1.2, patterns=['RAPID', 'BRIGHT1', 'SHALLOW2'], verbose=True)\n"
+    "tbl = nrc.ramp_optimize(sp_dim, sp_bright, snr_goal=30, snr_frac=0.1, ng_min=3, ng_max=10, nint_min=3,\n",
+    "                        well_frac_max=1.2, verbose=True)\n"
    ]
   },
   {
@@ -228,23 +276,23 @@
  "metadata": {
   "anaconda-cloud": {},
   "kernelspec": {
-   "display_name": "Python [default]",
+   "display_name": "Python 3",
    "language": "python",
-   "name": "python2"
+   "name": "python3"
   },
   "language_info": {
    "codemirror_mode": {
     "name": "ipython",
-    "version": 2
+    "version": 3
    },
    "file_extension": ".py",
    "mimetype": "text/x-python",
    "name": "python",
    "nbconvert_exporter": "python",
-   "pygments_lexer": "ipython2",
-   "version": "2.7.12"
+   "pygments_lexer": "ipython3",
+   "version": "3.6.1"
   }
  },
  "nbformat": 4,
- "nbformat_minor": 0
+ "nbformat_minor": 1
 }
diff --git a/notebooks/GTO/NRC_GTO_YSO_Imaging.ipynb b/notebooks/GTO/NRC_GTO_YSO_Imaging.ipynb
diff --git a/notebooks/GTO/YSOs/TWHya_spectra.pdf b/notebooks/GTO/YSOs/TWHya_spectra.pdf