{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# KAGRA PRX with Mirror maps" ] }, { "cell_type": "code", "execution_count": 74, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Once deleted, variables cannot be recovered. Proceed (y/[n])? y\n" ] } ], "source": [ "%reset" ] }, { "cell_type": "code", "execution_count": 17, "metadata": {}, "outputs": [], "source": [ "import pykat\n", "from pykat import finesse\n", "\n", "import matplotlib.pyplot as plt\n", "import matplotlib.cm as cm\n", "import numpy as np\n", "%matplotlib inline\n", "pykat.init_pykat_plotting(dpi=90)" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [], "source": [ "base = finesse.kat()\n", "base.verbose=False\n", "base.parse(\"\"\"\n", "# clight 299792458.0 \n", "\n", "l i1 83.74 0 nL #I0=780W\n", "const fmod1 16.880962M #f1\n", "const fmod2 45.0159M #f2\n", "const fmod3 56.2699M #f3\n", "\n", "mod eo1 $fmod1 0.1 1 pm 0 nL neo1\n", "s seo1 0 neo1 neo2\n", "mod eo2 $fmod2 0.1 1 pm 0 neo2 neo3\n", "s seo2 0 neo3 neo4\n", "mod eo3 $fmod3 0.05 1 am 0 neo4 neo5\n", "\n", "s eo1refl 0 neo5 neo22\n", "bs refl 0.0001 0.9999 0 45 neo22 dump npr nREFL #REFL\n", "s bsm 0 npr nprb\n", "m prm 0.90 0.10 -13 nprb nf0\n", "s Lp1 14.7615 nf0 nf1 #14.7615 \n", "bs PR2 1 0 0 0.6860 nf1 nf2 dump dump # tilted\n", "s Lp2 11.0661 nf2 nf3\n", "bs PR3 1 0 0 0.6860 nf3 nf4 dump dump # tilted\n", "s Lp3 15.7638 nf4 n1\n", "\n", "bs bs1 0.5 0.5 0 45 n1 n2 n3hr n4hr #BS\n", "s bs1bsAR1 0 n3hr n3hr2\n", "m bsAR1 0 1 0 n3hr2 n3ar\n", "s subBS1 0.0 1.45 n3ar n3ar2 # BS thickness is ignored \n", "m bsAR3 0 1 0 n3ar2 n3\n", "s bs1bsAR2 0 n4hr n4hr2\n", "m bsAR2 0 1 0 n4hr2 n4ar\n", "s subBS2 0.0 1.45 n4ar n4ar2 # BS thickness is ignored \n", "m bsAR4 0 1 0 n4ar2 n4\n", "\n", "s Lmx 26.6649 n3 n6\n", "s Lmy 23.3351 n2 n5\n", "\n", "bs ARx 200e-6 0.9998 0 0.025 n6 nPOX n6ar dump # wedged ITM\n", "bs ARy 200e-6 0.9998 0 0.025 n5 dump n5ar dump # wedged ITM\n", "\n", "s subx 0.15 1.754 n6ar n6hr # thickness 150mm\n", "s suby 0.15 1.754 n5ar n5hr # thickness 150mm\n", "\n", "m ITMx 0.996 4000e-6 90 n6hr dump\n", "#m ITMy 0.996 4000e-6 0 n5hr dump\n", "m ITMy 0.996 0e-6 0 n5hr dump\n", "\n", "\n", "attr ITMx Rc -1900 % no error\n", "attr ITMy Rc -1900 % no error\n", "attr ARx Rc 0\n", "attr ARy Rc 0\n", "attr prm Rc -458.1285\n", "attr PR2 Rc -3.0764\n", "attr PR3 Rc 24.9165\n", "\"\"\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Field Checks" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Fields in PRC - Sideband Lock" ] }, { "cell_type": "code", "execution_count": 19, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 19, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/svg+xml": [ "\n", "\n", "\n", "\n" ], "text/plain": [ "

" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "kat0 = base.deepcopy()\n", "kat0.parse(\"\"\"\n", "ad CRPRC 0 nf0\n", "ad f1PRC 16.880962M nf0\n", "ad f2PRC 45.0159M nf0\n", "ad f3PRC 56.2699M nf0\n", "xaxis* prm phi lin -180 180 1000\n", "\n", "cav prcN prm nf0 ITMx n6hr\n", "cav prcS prm nf0 ITMy n5hr\n", "\n", "maxtem 4\n", "phase 2\n", "\"\"\")\n", "\n", "out = kat0.run()\n", "plt.semilogy(out.x, out[\"CRPRC\"],label=\"Carrier\")\n", "plt.semilogy(out.x, out[\"f1PRC\"],label=\"f1\")\n", "plt.semilogy(out.x, out[\"f2PRC\"],label=\"f2\")\n", "plt.semilogy(out.x, out[\"f3PRC\"],label=\"f3\")\n", "plt.xlabel(\"PRM tuning [deg]\")\n", "plt.ylabel(\"Amplitude [a.u.]\")\n", "plt.legend(loc=0)\n", "plt.show" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### With Mirror Map" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "## Cavity scan with all ITM mirror maps \n", "\n", "kat0 = base.deepcopy()\n", "kat0.prm.remove()\n", "kat0.parse(\"\"\"\n", "## -45 deg is a tuning for SB lock\n", "m prm 0.90 0.10 -45 nprb nf0\n", "\n", "ad CRPRC 0 nf0\n", "ad f1PRC 16.880962M nf0\n", "ad f2PRC 45.0159M nf0\n", "ad f3PRC 56.2699M nf0\n", "xaxis* prm phi lin -180 180 100\n", "\n", "map ITMy ITMYhr.txt\n", "map ITMx ITMXhr.txt\n", "map ARx ITMXsub.txt\n", "map ARx ITMXsub.txt\n", "map ITMx Apperture.txt\n", "map ITMy Apperture.txt\n", "\n", "maxtem 8\n", "phase 2\n", "\"\"\")\n", "\n", "out = kat0.run()\n", "plt.semilogy(out.x, out[\"CRPRC\"],label=\"Carrier\")\n", "plt.semilogy(out.x, out[\"f1PRC\"],label=\"f1\")\n", "plt.semilogy(out.x, out[\"f2PRC\"],label=\"f2\")\n", "plt.semilogy(out.x, out[\"f3PRC\"],label=\"f3\")\n", "plt.xlabel(\"PRM tuning [deg]\")\n", "plt.ylabel(\"Amplitude [a.u.]\")\n", "plt.legend(loc=0)\n", "plt.show" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/svg+xml": [ "\n", "\n", "\n", "\n" ], "text/plain": [ "

" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "## Cavity scan with only ITM Reflection phase mirror maps \n", "\n", "katp = base.deepcopy()\n", "katp.prm.remove()\n", "katp.parse(\"\"\"\n", "## -45 deg is a tuning for SB lock\n", "m prm 0.90 0.10 -45 nprb nf0\n", "\n", "ad CRPRC 0 nf0\n", "ad f1PRC 16.880962M nf0\n", "ad f2PRC 45.0159M nf0\n", "ad f3PRC 56.2699M nf0\n", "xaxis* prm phi lin -180 180 100\n", "\n", "map ITMy ITMYhr.txt\n", "map ITMx ITMXhr.txt\n", "map ITMx Apperture.txt\n", "map ITMy Apperture.txt\n", "\n", "maxtem 6\n", "phase 2\n", "\"\"\")\n", "\n", "out = katp.run()\n", "plt.semilogy(out.x, out[\"CRPRC\"],label=\"Carrier\")\n", "plt.semilogy(out.x, out[\"f1PRC\"],label=\"f1\")\n", "plt.semilogy(out.x, out[\"f2PRC\"],label=\"f2\")\n", "plt.semilogy(out.x, out[\"f3PRC\"],label=\"f3\")\n", "plt.xlabel(\"PRM tuning [deg]\")\n", "plt.ylabel(\"Amplitude [a.u.]\")\n", "plt.legend(loc=0)\n", "plt.show" ] }, { "cell_type": "code", "execution_count": 79, "metadata": {}, "outputs": [ { "data": { "image/svg+xml": [ "\n", "\n", "\n", "\n" ], "text/plain": [ "

" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "text/plain": [ "

" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "## AS spot without mirror maps\n", "\n", "kat1 = base.deepcopy()\n", "kat1.parse(\"\"\"\n", "beam beam n4 # AS port\n", "gauss input i1 nL 4.4088792e-3 7.3067531\n", "xaxis beam x lin -2 2 300\n", "x2axis beam y lin -2 2 300\n", "\n", "maxtem 8\n", "phase 2\n", "\"\"\")\n", "\n", "beamprof = kat1.run()\n", "\n", "fig = plt.figure()\n", "fig = plt.contour(beamprof.x, beamprof.y, beamprof[\"beam\"], 25, linewidths=0.5)\n", "plt.axis(\"equal\")\n", "plt.xlim(-5, 5)\n", "plt.ylim(-5, 5)\n", "plt.title('AS spot in PRMI, with NO Map')\n", "plt.show()\n" ] }, { "cell_type": "code", "execution_count": 80, "metadata": {}, "outputs": [], "source": [ "## AS spot with both AR and HR mirror maps of ITMs\n", "\n", "kat2 = base.deepcopy()\n", "kat2.parse(\"\"\"\n", "\n", "map ITMx ITMXhr.txt\n", "map ARx ITMXsub.txt\n", "map ITMx Apperture.txt\n", "\n", "map ITMy ITMYhr.txt\n", "map ARy ITMYsub.txt\n", "map ITMy Apperture.txt\n", "\n", "beam beam n4 # AS port\n", "gauss input i1 nL 4.4088792e-3 7.3067531\n", "xaxis beam x lin -2 2 300\n", "x2axis beam y lin -2 2 300\n", "maxtem 8\n", "phase 2\n", "\"\"\")\n", "beam_mapall = kat2.run()" ] }, { "cell_type": "code", "execution_count": 81, "metadata": {}, "outputs": [ { "data": { "image/svg+xml": [ "\n", "\n", "\n", "\n" ], "text/plain": [ "

" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "text/plain": [ "

" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig = plt.figure()\n", "fig = plt.contour(beam_mapall.x, beam_mapall.y, beam_mapall[\"beam\"], 25, linewidths=0.5)\n", "plt.axis(\"equal\")\n", "plt.xlim(-5, 5)\n", "plt.ylim(-5, 5)\n", "plt.title('AS spot in PRMI, with all maps')\n", "\n", "plt.show()\n", "plt.savefig('AS_all_map.png')" ] }, { "cell_type": "code", "execution_count": 82, "metadata": {}, "outputs": [ { "data": { "image/svg+xml": [ "\n", "\n", "\n", "\n" ], "text/plain": [ "

" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "text/plain": [ "

" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "kat3 = base.deepcopy()\n", "kat3.parse(\"\"\"\n", "\n", "map ARx ITMXsub.txt\n", "map ARx Apperture.txt\n", "\n", "map ARy ITMYsub.txt\n", "map ARy Apperture.txt\n", "\n", "beam beam n4 # AS port\n", "gauss input i1 nL 4.4088792e-3 7.3067531\n", "xaxis beam x lin -2 2 300\n", "x2axis beam y lin -2 2 300\n", "maxtem 8\n", "phase 2\n", "\"\"\")\n", "\n", "beam_submap = kat3.run()\n", "\n", "fig = plt.figure()\n", "fig = plt.contour(beam_submap.x, beam_submap.y, beam_submap[\"beam\"], 25, linewidths=0.5)\n", "plt.axis(\"equal\")\n", "plt.xlim(-5, 5)\n", "plt.ylim(-5, 5)\n", "plt.title('AS spot in PRMI, with only transmitting maps')\n", "\n", "plt.show()\n", "plt.savefig('AS_TRANS_map.png')" ] }, { "cell_type": "code", "execution_count": 83, "metadata": {}, "outputs": [], "source": [ "kat4 = base.deepcopy()\n", "kat4.parse(\"\"\"\n", "\n", "map ITMx ITMXhr.txt\n", "map ITMx Apperture.txt\n", "\n", "map ITMy ITMYhr.txt\n", "map ITMy Apperture.txt\n", "\n", "beam beam n4 # AS port\n", "gauss input i1 nL 4.4088792e-3 7.3067531\n", "xaxis beam x lin -2 2 300\n", "x2axis beam y lin -2 2 300\n", "maxtem 8\n", "phase 2\n", "\"\"\")\n", "beamp_hrmap = kat4.run()" ] }, { "cell_type": "code", "execution_count": 84, "metadata": {}, "outputs": [ { "data": { "image/svg+xml": [ "\n", "\n", "\n", "\n" ], "text/plain": [ "

" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "text/plain": [ "

" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig = plt.figure()\n", "fig = plt.contour(beamp_hrmap.x, beamp_hrmap.y, beamp_hrmap[\"beam\"], 25, linewidths=0.5)\n", "plt.axis(\"equal\")\n", "plt.xlim(-5, 5)\n", "plt.ylim(-5, 5)\n", "plt.title('AS spot in PRMI, only HR maps')\n", "\n", "plt.show()\n", "plt.savefig('AS_HR_map.png')" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.16" } }, "nbformat": 4, "nbformat_minor": 2 }