notebooks/paper-figures.ipynb

{
 "cells": [
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "f7bdaf49",
   "metadata": {},
   "source": [
    "<img width='120px' style='float: left; padding: 5px 20px 20px 0px;' src='data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAJYAAACWCAYAAAA8AXHiAAAAAXNSR0IArs4c6QAAAARnQU1BAACxjwv8YQUAAAAJcEhZcwAACxMAAAsTAQCanBgAAAAhdEVYdENyZWF0aW9uIFRpbWUAMjAyMjowODozMSAxMTozOToxMIfy3oUAAHwgSURBVHhe3X0HfF1Hlf73en9PvcuWe+8lduw4cSoptCQECAsLAZZACKGEHlgDu8vuUpZespRQE0JIIN1JHMeOe7fcLVuyJKu313v5f2eeZLWnatnJ/o9/15Luu+/emTPnfOc7M3NnNC+cfzFVbqtAMpVEj0STYTxR+11kmQpwbeH7kUgluj/JICmg0v04GkOVmOa4FuXWldDwXzjpwzH3U/DFGzHZthqznG/hWS20Gh3qAyfweusz6Ajn8AYJFJrKcNrdgbVFi/Hpee+DUWdI3/sSSzAeRkOgFae9tdjdeoTlalHnc0wuTHeWY37OdBRZctXfLqMdepb9jZauaBe+c/TbLHsALQE9wnEjKlxNyLP6sCT7A2gNHcP50D62ZwLFlvmY4bgBVl0ukoh332GwSHtJm29q/h0sOgduKvkgdBojP2HjDiFajRZnfKfwTP3fcGPxrVhTeDViyVj3p4Dui1/96oZUnxvIQ+oCJ1HtP4iZjhXIN09ioXqNrr9ocDbwGq/fhWLrIky1X9N9HjBoTCgwz4E/3oKm8BEE4q3IM82ARqOHTe9EW6QavlgIsYQJ+RYjpjtmY0vTfuh0OszOqmDFJr4RY8k4OiMeHOk6g7+f24wnzr2Cv9RsxN7W4yxLEFMcpXQmJ5pD7TjhrsFrTfvwcsMu7GytxK7Wo6jxN8BhsNLIHEqxb4Q0B5vwQsNz6ApbEU1Ql8Yoiu1NKLEsRKF5HnKM05BnnolI0oOW8HE0BA/QyOIscwl1Kg6b2Vj0WiPbqpNteRwVtgWwso2GMyz5zKyzorLrEM0ghflZi3mm105093/5gQ0DDeeEZxc6oy1Ymntjd2EGixhgW+QETnqfh91QgDnOt1LZ+u5P0yJ/55tmIpoKoDl8FO5YPXLN0+gVLj4xgbZwNcIJM8+7cVv5OlYsycbci2pfgzIAE5HLrDOqRtRqNN13Hb0ICkeSUdQGmvFq4x78o+41/L7qWWxq2I1zgUYauAWrCxbibRXX4O5pb8E1xcvV3zeVXolVhQsxJ2sKSm2Fyvsbg6041HEKLzfuRmu4E5PsRbDrrXSUsZdrvCIAsLd9Lw39DB3SDIM2Aac5gBxzkBHhFv5tUY1r1Frp1PPZLkXUaSujyWGW+QRchhLqPlvdaaBQw2yRJM4HT9B5coh2U/oZSibRa/W8byPbsRVzsxbwuYJyadF97Isf39D9uzKWcCKAvR0vshB5hNHlGRUn13mJRAe7/szPtZjvupMGIBY+WOTzHGO6kC00xK5IDT1qGrKNpWgKHWWoTCCS0FEpSbx7yjt5ro1hsZbKO4zNjXtpEHvVOW8swLul4GRIGsnIuqJeHO44jVea9uCR08/g77WvYm/bcSo5iKmOMlxfuhLvm3Yrbi5bi4W5s4i2eUpJ0nDpf1BGU2orwGzXFKzIn4cbSldhZf4Cde/tzYewreUAw32UKFdCB+hV6KWUBB3lz2efQ2MgCIteA502ToTtIv2oIGIt5RV9DSbFEJhDFJurDK0tWsWo0cG6LlDtN1Ck1ha9gw53hO0RJK1Z3P3J0KIncHhjXpzyHsM05yy2c273J2z3ys7jF0ojRtAaqsPT539MaFuHRdnXsTIDY7MG8WSIxvcbNlQbLfUdyCX8jiRS8PrAbtQEttAbFmBx9ntx3L0N+zr2wB1xEXqBz859gIhWwLBEhAu2Y1/7cRWSOohevmhA8Zwsk4NGqGdZwedmURmmtD75d545W32vyltH/uZW9SkmR5rhmqyMo8xaCLvRqr6PVNqAxiJaPiRG9DrSeRp/qX6J3t2CSbYivGvKDbiycLFC10slKZZ3Ex3l58cfJwclL7XpiRRhVDjqsDL3TuqmvPvKwSJR55jnKRUar8p/gLrO49nBtRcOvLfzOZwnB35r2X28Zz6fOzRqyfWNofN4svZRrMhbRW52G+Ld9qL7+Bfvu4BYEnIOdb0KT6yNRrVeEbnB6tfguPdpVqoKU+xrSW4XSpuOKOIlWcbJRMQuVvAEjWI6+VuFQi1/TINAIg6H0YS5rnkw6Y3kXdlYQjRZV7yUhZ6HRbkzqSAtAuRlXREvQvGICpfNoQ42MOGeBnXGW4+2UCfLlIe3lF+JO6dcj9smXY0rChagkAZmphGKcYxXRBOCliVEsisLF9GQs3C06yw2McxWeeqUYWeT6Es5J1JCiQj2tB3FwyefZFiOYla2mQ5vIcXwsA2cqGA7ZDKUvqLXmqjrSqKXWXHdTG4l7R9LRsglK9k2ZbyuNON1vSI8y0yAqCQaBrA4d3n3ebZ3D2JJw/vjbvzj/E8YBnOxvvAD/Fp/axUEqPZvxSnPC2yoBZjtuJXnRt9Q8gyfCqF/JAeYxSzmvWyYzUS/w/BEnDQmM7604IswsfJ9RRpTvEMqMrCiSXqyN+ZXHi2fCBo5DDZen1Z1+vxwyhmfSIn0Wh3R1I0X6raTe+1SnG5p3hwVYudmT1EcbrwiZfZE/Qq1X2vcx4TjLEMPkG0JYH3xVXjh/E7W040bS69ThH0kPiSyj1FG2vCKvI/yXkT6ASJ1iiSD2Nj4a+SaivGWkg9niFiDZVPTiwyHx/H+aR9mVGCyRz1cQCzJwk54d6HWfwwLGQKzjAU8279BQozRle4nyD8KyKveyQYfe+ZmJnGX+wiZzCaC5ZoqCKdHafF6ptARFFhymJVUUE29ilIGcuFfz9+9pRMkkkPCohhW77U9V1waEQXKcxflzCS6zlbIsq/tGDYzm6xkuJSyCKpJ94k02mhEunZaibrP123DL0/9jYnGLuomwpBnZ8bcjCyzUAGXCsMVDgtmuVbyWyPfWxKpGI2mKVxJJJrOcNjLh/qKkQlAa6SO0aCRCcocmMnPhhNBObGd4+4j/KnHTHIt0bsyrB5L3dr6N2YS2ViacyO/0r9R5OIT3meJDk3MQG6mN0qcHo8IfJILsYJS0QrbSqJYKz2/jaTRxBTaz0ZaPC6jfaMkyTplGR1Eq7lYwxAp7XzSU6uQZi8RJ5FM0IHIIw2WQSYgeo3z82gyhhpmw0/XbcFvTv0dO1oreU87bi5fg/dMux6N4W3M1IOkEPl0wJTqz7qiYB6jS7G6y2jEqLPRoY+SB0VIF+Z3n+0vEoEkC67xH2ZUEccvUWUcTmx6O+qD58jNarEge5GKOMqwpBFrmA2c9u7BvKyraNHlA26mUdlclW8TUWYqptjWjviw4cSotasMpTVyHDmmKSz8ZNQzJQ7GzVSYm2FkBs/nXtQzLrf0lNXGbFL44BXMIG0Gs+J9rzcfwG5ypJZQOw4TyaRP7HVmlc/WbcXLDbsZ1rbjpYad6u+jndVE7Wy8b/otiiMuyJnBMLgVh9x76YhGlNkmqf40qy6JZblXKsQYnaSodxvcsVqG2AYUmeepbHGQ0LAEeaQ/S1xmsm2OnFQfDSWCWELyT/mOk/BnY7J9StqwRCVbWx5n+mogWt00AC3EgqM44nmSiOInuX6bKuDFiCCkRZ9FEn9MIdRU+xq4o/XKqCIJA0OKBytyVygk+L8oUmob0Wlu1lSszl+IMnsBzvrOY3/HCZWxVvP3en8L2pl8SJORBqqMcnbWFNw59TrcXnG96qw1sT26om14vvER0ockkcZAFMni99vZDjNRbi9VzxutiBFK2ZrDRxSdcRrL0h/0kxSfa2PbVKtwOMO5XHWejiQ2g4N1PM1svJ3UYCk0R7tOpYRXvdT8CHnTOizJuUFlBj1C+2U834Mj7qdQbluB6fYb+OiRieJIIsZ1yvcCs7lDJJMfYTjQ4uWmJ9AVzlH9Ne+bejuuLrqKITqirv2/KlJ2nVaLILPZTmazKVZFjMhBZJPfdSop6b1OGl9Co4gkHk83/C+Ode3muRx+boJGOFhYh5tLV6PEWsy2SPGbco+e5Gbo4Td5htCPvZ2/VeCwOv9efmsw5RDEOuHdicquV7Gu4C7VpyXhcTiRe29vfQ172nfivVP/GVph/dLTbtJYVYdof6MS7uVBlX8TESYbk2xrJsSoeqTMshwGnZUZxUbkMr0tthTAZfKwYhr85dzT2Nm2hwq4POOGl0qkqcVQjDSmEls+Sq35auxROmQNbEAxJDnS3CZ9rYicO8dwdMKzF3ZjIdshye+EiexxhkEjOVKhurc4fjDeQBT7LQHgaRrj0FmcXG/UOdRoiAy1eRglMnFZMaISy3RldPXBUyMalYgg4YLsJWxPA/a07YC2LVyPhmAVM4C5hDNX+qoLomHldiAc92KSdRUVMf70eaBIJW2E40LTXMJ9Lb35DKbZl8GkDyLbHFOV+d2Zx/BK0+uqh/f/BxHDkUOQaDhJO3SYydTTqnG1cKWNQiNDVBoU0TjTBpFCW2grzvh+xaTqFNrD29EW3qHuMJSkqNciywL1u1CRTB2g8izpw8wxFRMda9VojJRpeEkRhZ0EpzmoDVRDt/pfZm/oijXjqoJ38AN6EhtR0mSjTk++04lDnX8jipRiftat6rx42kQdBp2OxswMMXQMgUQHw/AtNLIWKqmRutGTU0gaW6UKvTBnjho7lNT9//dD6nnMvYuIvRHTGYZaQn7qK06j9BFl7FhdsBhOQwLn/H+lIW3n9TmosL+d3LSViHaMyZeARD7bUjdY51odo4+djiwcqg4znVfy+5YL7d5zWHgumPCiLniCBH42uV2pIul9rxl4yNCWjd873LUfujkfzN1Qbp2JRdnraJMpFYbkkEIccv+dWUQDFmW9DQ5DnoJr6aycyMPGEBukUTXSuArMUzDHdSXh9wRC5ALilTqNGUfdJwn3QczLnkV+YlIV/P/5CDBCPH7upzQyMxObpeS3R0i0Q/DFDGwDB4l9EPWBJxCI1aPQcgWmOu+gHktpTCXoDB+GL17N8yvYhoZB+pZDBovjTMik/zDHWEpDnMz79rZ9+tDyOhNOe/er+8xyLpGEUZ0f+tAgy5SNU54TwrESbMwr+gEdH8/MoQq1/v0oNs9mIWcreLwUIlM6ptivZOHNJIybVIfcNYV30iOFjEZQYLUwds/Fs/Uv45cnf0+P9aa/+P+pCIfd3voCkagRK3OvJWKfYgPrmKkFaXAWIoOHIe85GoiBGdv7mNq/lWjDRICk3WYoR6n9ehpcA875/j5k+JI2L2E4NFDX9YFDbIPeeVQ9IuVwGnKQTX7XHKphRPENeb8eEQsRMFidvxa6d3zqxg0rc2/kbXpjbYLWvLf9UZVBLMt9N/T0nImTwYWTaTThhIdc77CaN1Rum0e4tjJWn4I3GsJVhauIZsXY3LKdfOIcIddEZcj0EPHINPkdjSiOw3+SdcpgaTgRRhdT/lamyOeDTSTAtWgNtfPzuIJ2UaSg9EgKnUhpDtXhb7UPo8w2jWFwOV5r2Qy7IaxIeZcMe5k9mOqYRtR6F0m99Df25Ugp2PWTGBLbWKd91GE+/5YuhcGgIGOH/pjMlTvBpGyJ6mIYKJIdhpMBnPUdVv1ZWTIoneFeA8VudED32a8+uEFmivY8XNDqXGAvIXArK3AlG3mJInwXK3JfkWQqwsYanInY9Dls3MPMVjpQZl1EgjmZRNVH1DzPcOzBB2d8gBU14HDnUWxpEf6xl7H8OOr8DenKsu0lfEsoibERgnEJHX7VNyYdk+f8dajsPM50+CBebd6G15q3YxMTg5catqSPxs3Y3XaQDblDXXOg8wiqvDXkLB5liOLlCu6JHpfK0CQj//v53yrjurnsbvLbY2gI1RI1/KyLlUhkxrriZeRT6+hM0rc0mHhLyRyGKSz3SRqidEDPo5PY0x/2EUEWrVarEMuiczIqTR9kNDJHS+Jfte8I72FicrVAOeZIIhMANZWdh1J6NpjcVBQWJuS90vwDdYOrCu4lKow/E0wbE9NoxOgd9egIV/LnOWSb5tMjr+MTewspFa3segbV/h1Ym3+P8qI4Ifrp+r+wgWtxS6nMn3q7MpBzvvPqZ5W3Gu2RTtXwMuA73TGFCsqHJ+aDO+JGJ42iiz/FKBTuCPrwEKKaZXASHeWQaccONSTjMNgRioeVMdcFztPrO9R35chixjzZXs7wMwWrCpYROXLp5RM3TUb0f9JzEI+c+W/Mci3GtUXvwnePfYeljtCgoyTaTtbRjrsq1iuCP5yQIdGhjuOM98+s13TMdn1E9YENRC7JCF9r+QmzfTOuK35AOe7AawS9n6r/ueJl7614kGdGdiql55PuEylRnDrBAlW6n2Fcf4W85i1ErLVEmLGhlTSg2LrcPMDsrityCm2R/QjGmljIKFHFomL6DNf7UGBZeiHdFRTzMTt9vfVhkskK8qxPqDt1RFrwaM2f5QrcO+vj5GOz0o2dTBCmI+QiMv+qBmc8NST954lUEdUBKUZjp6GIweQYyRXUvHUHvT8bDqONNaXHEoFkhoIk9Fr+lG4AKXc8GUM0GWd5fAzPLTTkWv5sRiN/l7BpZeYzP2c2luYswELyP5nOfLFdIoG4D7+u+g9y23p8ZPqXaWTV+Ou5J/mJmcitY0iL49ri1XS4kaaypEVmMdT5niP67cQk+010ilv4rf5oIwPTJ9wbqb8tuKboPhSaZw66Ru6zu/1FZqk7cHv5fYxgM5X+RxLNsa6jqpRiEPICxMbG79I6rQqtRMljETEO6aLwRM7SmA7AFz3HBvISamX4YBpheQG9wkhPelQZ19zsj5Ig9k6BFdTa3/E4G/Aorit6AHnmCp7TYFfbdmxu3kylFtC4Pkfu4LhgkGIcEp7iiTi8NATpYJS/zfRQ1T1BT5OhoZ7pM/I96fEeqS9JRHQi4U/KJS8KyATEav857Grdr4y5k2hYaMnDkpz5uCJ/KR1xMqzkhvKdsYiUa3f7JmaCP8e6wluxPPcG/Oj4wzgf6KDhg/yyg45eiJmut1Nv/Rt+OEkxUpzseoRG24B5Of9C55rLevcahbSXO3Ie29t+Rd0uZbLwXp7trxcxrOZwLZOnX2GmczEjxz2jAhvdfV/6hJo2Q98lMd7GmHsAc7NuoqfLONTIyu+RdCFP4aT7ETSFthFRQirWl9NbSq3riUKL2dBZ6jCTVHZGjpAHnCaCzO7mAGm0MDA7rAvsVz9LrHMVFOebC3mujulxE1NpKil7RfqhFGmUHmXLFBar3qzIvXAtEfEu+Vz+pU1rbCLf6vFQ6acpsRXRiJZgce58okCpCrdHuk7gVfK1I10niaAdaYPmYWByMRrnFLR6ovYXMDBJmmxdgUernyXfbKX+4ii1taDM0kVedQPvmcWrR18DaVO7cRLD4jHq+wTbYBGdri+1SVFfWSxzFa85TzRanIH6SMdnFulFO2r8R0lhpqtp6yNpstuwZMptmKT1T0SDLNUZKpY6WhHPbgsfwCnP7xUilNmuJ/zeylC3gmglvcTC4aRxlBnQAPLY8GaGucPKm3LoScIBBEQkQ2yP1vB+Zxn2VvA6U7qReI8zPnKqcANyzU5VwUzeO1bDGasI0slz7QYbJjvKsaZgOeZlzSLyOhkum5hQHFNJwTH3KZU4SLllOEdlmTSyngYRXhhNRln/IA1yHynAawhFXczGT8Mb97LBEyi2tbIRfQzhxSi2Xi1PV98dixhJzEWH7eFD/HZckfm+90knUknSiIM0oDwVJQYajUQFB5OrU979jGpBctmFPDu8wyjDklhbS6Sq9u/ELOd6NpzcfHSQK9c1B3cwvD1Ba7dhmvMu5JoWqsqks5ZMypDhnBIqN0pCf5AkPaKMS0THBgCVXh86SKUWEa2mEP3ihPFsZqvM0hiOfORuC7Ov4DPeuKEeUX5PeM0if5uTNQOrSernuGbSKCwqE93TLi9d7FHZXSgZVsYkicFhZqe72vYzvG9jcvISXm8+SAphQTShocMUUGukD3qZ4+8nn4vz3Hres7D7yWMTaR/pPPUn5M2no2zb+TRyR/enIikirJPlPcUI0qB6AfQq4+wVqaOV9MMX76RzV6LYUqH6t+T8UKIMSx6+t/Nx2mAK87NvpYWOpsFk0DROHvAKlfU8OU0ujeo99OTRkUv5vsMwGeF4G9GJoU9jg8skKBSnUsX7K1V3w2SrzKNOqbAi/KXKW0UuFSBnM2GKY05G1LrckjYv6bnW0eNzqMM5WJ6/iNnYdNVfJtnlrrZ9ypi2t+7FgY7DJNWtzMCj1EKc17jJnybjlrIbGbp0TO9Pq552O43LzgYtsV7LNhl/BirRx6RzELX2U18RGld6rLBHZB686L0+eJjIlJsRWAS1xLhOew+q0D1D9cQPHdV093/p/g1NoRM44X1F9VuVWhfyIcOTMwl98VSICvgbmoJbaUzlmO56jwp7ozOqtAgpdhinwB8/h/bQAWWUVkMxjcjCtL9L9WsVW+fApsthNWUWZr4aPa/2nmPDNGGeaz5kHtBYnnmppcfIJHwXmPOwIm8JluYuwCwaWYE5F1cydF5XcjUzvDU81sKXOEUHqsc7J70bJZZSbGx4UTlUjonhUBdicrCcfHc273wxDiSolEuDOM8s/QRy1CS//q/r2Q251KmE7wbV1TPwfVKpk7zE6mO7yKyLYssklrFInc8k2gSzs5PezaoxK+wr+02bySTSNyVv2pzs+jVag3sJifPU0IJJm81HjK3y6QZwEnluZwg047TnMfhj5xVPE0iWp9X49ygDFBGDX5GzEsvzlhPWvfhb7V9V73m6i+PNJT1cTLJJp9GBxTnz8NZJN+HKwhUKyUqtJWzoTiKw9GrPZkNNxVn/WbRH2hm6omz4CHVqJQIuo5ZGzsJGFg2KLFeyXAlGmU38u7etpB0E0aY4VpBqnCdyHcoYtQS1Fmatpa0YsbdjE0P70FN0tO3Rc2gJn0Y5kUoyhOGbSKPI9jH3wwxH51BkXYupzjsZk9Nv4I5HpKLCASbb30YUDDAB+KMyXJexDPmmaRd643tEsrs1BVfRq0roECfxxDnpurgYb770IkYWZz3j5IrSHSJZpoTIQ51blSMvyl5LfhVFZddhXhtjOwRh1MSogxk0MMkEL15kSMhpnK66HDrCR8hTa9mavaFMnLbcuhROfRFO+7YysfDzbH9rED3LsM4s53JSlbNqgHoonqs97X1Nwd4U+yp5Om+V+Z+8jyej6ce7fqX6qibZbiKq3MjzaTTpe+1Y/0ml88wLVTYpiNUYfJ2KNan0NxhPjyHKCHtP6m7WWbC28Cr+TOFg5wE83/C08qZ0n9Ob/xAu5o3JG087iVxTGSlmE6nacNZ7hhlzkhxTEEvLxGUp63xxuu37T8yoxHqV+r2RFIa5Kc9Jy6b/yXjhFMdKuCMNdOhKZTSi876HRI85NCyr3k7db6GDRLvrJZ/1HrplH8naUEoeMz9rPZFHpkpIp+KAg3xBCOZRhr9owkO+cAdK7asJiZnn5Yzn0Os09NDJDHHH1fDPJPsVjOFlDIV7yTn8mJt9teofEhIvY4I5pjw2ToMixjX+WuQYsxmSZ6hhlox1eBMdMpf9UNfrOOreg2uL3km0no6tza+xHjWsl5+E3cv6TCLnvYH11ao6Z9LZWA8972U15KgRkXZmiDJrJctUqj6Tchl0eiJWDkNhJXXbgrkMeyaduV/Z5VqnwcV7eFHlr0QpudZk+wwmX7p+19GIU7TAdbTOHmuTWNr30DBza8fB9l+o9wFn0qiKbSto3/Qk/jdhB//JiHuJdSXDbTPaWHEbQ/NU+3KmwmeYJJxUHtTjETK+uSr/BhpjgpiZwhPkWzvbdqj7CCK8mQ9/3MvscCNDfQn51mJsbHwRW1u3Eq0SRKsw65YkH1qmul7SSJFBX+M6qDv+m+RYz/91zOY3p1GLn8khvQJC4qfaVzCDrKWhS6iToa8e2+g5tFiat45JlZ1ca7PqNpIMse+hLTBPRYFlmor7g0WLYKIN+zt+xmygCRWOm8irVqh4fClEeFq+RbheIc4xS00wNa5gOBQieca3s/uqtKhr6W3Sn2XW+3lNCo/XPKZ4ihDmobKVN1qkXMfc+9QY6MLsVdjZKq9/vcA6y7pcbqJZmNl1FvWwgHoe+S3ksYrozWmYROoxj1ngQXSGzygj6xHJvqc7rlB8+5j7FXKtIM/S6vqI1EFeaF6Su47oVoXjnn00wP5ZpO7Bh76wwWnItEiEFqFEO/a3/4Sw14ppjlsIeetVwS6lyDuHkYQbzeGDzDinEa7n8PezaGKCMcW+lEp3qoqJ0HmQaywk2TxILxdCbMCutp044TmmZpxKCJFebwmdgmRvBgnFA3jm/O/ZfHGGmUK83PQKM8I85Jo7RePkjSGU29fy7/QbxZdCRBM2QxEjgXTMdqkI1CvyQrFDkfdzgUMKwQR8BttHiu1TgJOe/WgK1WF+9koVRXpE9+WHHtow+EtQmZkYlT/WfNmMSkSUKcvvNAZ3KVJfzNAoxL3Ku4OhMYdkd47yqh4REikzHc76D2Npzgp6UhEaQw3Y37lPhcbKzsOQt4bTPKXbyMQiKW+EsZ32VuKV5idZURfrVEe9VmBedh4b+TgbWwszQ880520sq0y8uzSGJWLSuZRRNYX20ohnU4/5PNv9PKK/XZ+L2uAheKItDI3LlO4GisyVl+6qY569yCI4yTBbjzPoHvjyJ/ut6CciszkPdPwcHpLoqQx/Ffbr+jXmpRYDPUYIprwOXmJbqbymIXic5WnGTOeaC4YhImXPMRXgjPcwU+gO/PO0T2Jl3pVq/QejzkS0ayRJPkJyvEWFSSHIJ92y0l0DOcY5yKJvnqiHCUKUcN6bfituM8GGF+Mznm34I/liJ3mkljyqGLeV3aimKsmTTRqSdiJVmW31hQa6lGLWy4ss+4ii7ShVqNVbXxMjh8zNqwscptEUMzRnXojNacjGKe9BcuImLMxapUBARPfJL3/swjJGcuN4MoyDHb9EV+QsOdWNyrAup1GJSEorMVtQSyahldtWMbS5URs4iEKzTL+RYaPeSspLBzK+VeneoYjlbNcShtAiZjXzsDr/StZhOhwGJ9xErtpALU54j5MXHMMRTyX2te/FjrbtNLzNikRv4c/jnqNMt88rpJMO2HRINV20oTUGa/B0/VNsSAORuAh3T/0nGvoO1PmPwGUw8u4hotVbYWHmdulFOkVlgZZ2tIQrFe1Ij0emDVpqqlCL4VAceppD3ojuz6NELHobWyJJPraHGWYeI4qsfJOSpSJ7V/QTqex4RGVkFUx1pzpu5iWX16hE5JkWfR46widYqVp6E1GLlTzt3akI7TTHSl7T33uyWam64BlU+4+pKbQOpsTSMSlGWmAuUEa2NHcZ1hSsxTWF67G2UBaWW4TZzjlqhZRia6lCEFGSzD495j6K3eRru9t3YWvLFmxpeY1KlqnKbjUrQRBNptGIIY8kouhwPISn6v/KOtSSQ5XjfTQqLWLY1iYzShzMsOI0Xpmlei3vfZkG12k9Jl02WoL7yac7GR0k2++tj0lnJcqGqdN9NLwSyHpmg1CL97CT6Ff5DjOsnsOynGtYfqYDPYbFBJ4f7MNpz1PM/JZjput21cBvjEiXgZlViKMxIOtAlTBzXaDifV1QhkAW0bv6j3XJnGyJ80fdu9EWPo8FhGWpoIg0rPQaq85JrU4ReguVliseZitTPGe2ay4WZC/EFXmrsCLvCizPW4ElOcsU8lmIiIFEAPWBOuzr2Esj26xWItzTtptKr2ZYCyhUkyEOmQbTGe1Qg+UHOvcrg3z+/LN4teVV6raKxlvEcP3PDO8O8sBnGJJPqmkoXZEjfNZila1JbnZ5RJaHzCYqu0k7DhM1Jyld96KWTBSQdxGOoSN6nuVcrRx1oJiJWtFEhHWuVP2LpbZpPYYlyzH7yKt+QQVZSSbfp+ZHXS5RXsBGkflCcoAGIMZg05HUhg8wDLYzHK6FSW9lQ+5l47kx3blKXdMj8ruLSpB58se9e5kAODDJNnjlujRzSoe0HoOTQ7pbxDB6DFBmUshUnemOGQrppKd/We4KZXwV9inqHhIqq31nVTh9ufEl7Gnfg01NL+O5889gPw2wynuaRt7C2iUUX9FoQrit/DaUWafwfC22tz2qpgVlMQx6o/WY4Xo7jX7wiw+XVKgKC5Ol5uBe/qpFobX/XCsTyyNdDjX+AzT6ScgzTRKtdX+aFrk6y5RPND7EJKQBi7OvTBuWWGGV52mePMQGuw3ZJJCXzmuSiDEbCcXrSdCrqFBmbeFdaAu+gqbAU6zgP2jUVnr0DKKWkZmLF62RQ4RiWXN9LmTJaIHmLEMhia6su9nHcBie8kxFOOs7ymuOY4ZzITPJzIvujlYkNPbMQpVpzlnGbBr5JIbRJVhXeLUytikO4X05DAl2htRiLM5ZorjdNUXX0pDeTkMsYmjdpEYG1hTcTDNL4FDXC3SaMyS816M1tIcG5WCIlDDS26iXR2RKkpVOcoZhvhpl8o4no0WvaGh4DsW1kskY67qM5/o7q4i83OqLeehMh0g9yqD71Jfv2+CO1uBI1x/ooTNUz/qlCYGCih044/4ejedpdIa3Ef73wh89xfOdyuyNRCjxBm/0CLONeTCSZ8mLrM3BfUSiEInhFWxYmaZcyWzvjEKtgZPSZG6RDJQe6doFb6wL87NX8exgRVy8CN7JGgcWenkhZjlnKzSbl7WABjSTBlaiEC+U8OOJ2l8q47y17AOwG7PQETnPhOFRFJlnYIp9Ds74nqHTrFRpf+aySlLlp5O1UB925XADkfhiRGb4RoioHZHjdNgKOPu9r5hex6E1XENHOEvevYIZY+ZV/iQ7Fyrij3ug+8SXPrzhcOevEU56qJgPKM+5VA1R7/s9UeokC76ESl+NXOt6prHXU6HrkGNZQwRaw3hdTIPbzZjdSuSUXvUsGposXXianrAQMutB+k9OebfxmiCTDJle07+82TSsACt31LOHiGXv179yOaT3WRq81PQozniPYF3R25gkLFThdmfb40SIRvK5d6E9fJDcsZ5OcgvrJYuyDCynDKk1o9rzQ7QGn2disZdhVQzMQaey0SgGv9Y1diFqEaUamIXLsJlwvb73THNVDSMB9WnIRol6M35wRDPrhGuFcITGpa31v8bKyduw17DBJ2f8wsWLDh2hrfCE98NlWoEyxz/RkK5io09XKCVK0oJeSKXb9LOQZ1lPZR9Ce2gzq2MiPK9jgf04H9jO8qUYWpbQ05cwE9mJc4z9maZurC24VfXKv97yLLPcRp65vCFGynTCu4/86zVMc84j71hDzabo+dXkZXtVR6+gb0uoknWWt5gGc5e0UZ1HtfcHCq2yqDupf3voFZzq/AZOd/0HGvyPKlohupPRkvGI3NOsz2F5KohaVYgkvXxyr75kflyReRocxlzU+g4pzpVZZAbyKjqIGbpVH4lusBmKMSfrPeqDiRetUk69/3fkTjaU2u/u9rKhDFgylUnwx04TqQ5SmYtJpItpaNXMts4QfYQDmIhKpYpQtoaqMdW5nB43ICTqTAynOapvqyPchDmu5WzswRnNpRDhScI3/lb7c9X3c0vpB2A1OFUDCWEPxruwKu89TOVJ/v0bmUVdyUadzm/21b8GkW6kiqU8DK13oMB2i0J6q2EaeZGTBtBCHVXS0DbxeSeoVz11J1nd2J1IkjXpGG8NHUC+Wq2miGd7y2NkFu1nElUfPMrkYx6zxcxrxctIiIdJje6mT5RuWJB9Nz2mRClEMqKJPMSA6vy/pVfJ60V3scDToNVmeo7kJHokkpX8Sg0NJaRmM0biZ+DSNfN8G5HHTaPqQJF1CT1MFi8z4qx/N5+hYSiRnmPpt0rfT+qSay6ksgI46T2o+rUm22Wt+MGzDSb6kGc/d/6PqAkcx3XFd7Bs81luvRo92NP+BKbZV2JW1lp+/ipDWw0/v5nGkpvWAQ9BuziNrsb7cxqf9IrfybpIqi98SMuQk0euRmQ3r2W9mORoddRvLSmELFTnYltOV6/P9+p25EPKLNypkdmhjnots6/guV59Sr2Ev57x7lHOMs25jJ+ny9P3kOuyjXnQfeTzd26Y6bqJJ2Qqim5CDylIe3AzmoIbGd5WocR+gypUpmv1rFQq/jqCoW8iHtsCQ+oE4jQYbzwJTWI7bJoz6EoUIxCtQrk1DzYimUw7kZVQzvkPosI2n4ou471k7lF63pZkcTI3W8ay6oNnsCBrpepnkfOX8jjp2YeNjY+RyC/HzSXvVcqW8mxt+T29vot8625mkE5UdvyJYcOGednvuKAHQd5Eyosq94+IbDV0oregxHYj9TNAX6ynjHvKwh/ZpkWkGPOJFEfI3fbwXB7PzeI1dNV+3xn+kAH+tvBxBBl2pzpkIqVVPUMOmcIjL1NIv1sXqcXcrDXM3F0sb1rXPYdcm82QqZ2d9Vbl6RMtgj8BIk6t728wa6XT7K08N3hI4IIwawpHX+JPDWzW+2G1fROTHA9S8Q50JkrhtH2Z6fhq+JIGNAc28+5SCRNWFzCLpWdtaf2jms3YV4Q7yFjWTcXvZqhw47mGP5KrhLo/vTTSGWnFs+f/ABvLLc+VskmIkmGbGqbss5yryVemo53JSDDexmxwEQ0vHcaF18SY/Z3q+gl8DPvFtltJHUbXPhZdCVHwAaXrM57/RVNgE+82th58GUYrtEhnrSRLNSxQX84mQ0BWRTvcsRbUB47x/pmphZRWayVp681iJk6k1/ys5zdsbD/Kne8mzMquU0OJlqGuFtFYJQz6JTAYroNetwAm/TxyqreSY8TQEDpGwruSlbegIdyGROKUKnehuQLLc2+F7Fq1t/1phQ59RWo207VYzX067t6PPWpi2sTPcxKR+77U9DjT8kbcWPpu1acmw0pyfm/HMzAy85qfdR1kykxLsJI6kmUbe96KShvVafcv4ImwrjSoMtvbeH60bSMNX4CZYly6QhrXw0wMXuVdR0/oRZ8FalMnGxoD+/jN/t+Vukwmv5KQedS9hWcy82RJQrSXYtKeVKbBL6lxJaF8PYnpfFXooUQyv3DkORY8QqOSPQ/lWsmhYvzuXIa4lcwqZWGRMwx/M9ERdcMdPsZrtKrjcp7rapW17O14Fk2hs4OUKcZ2PblOrqmQaPJ7NR9K+rjSz5kYkfod7NqO3e2vqpC7KHu1cgjpfK4NVKLWX4mZzlUqXEfjkuHuhWys5FAbAMjMBw+R6qfoDO9FgXUdw58s4jHWthHjyiNyfYo/82lcv2GWt4/nR0fm5XmyRLeVxLwpdFj1bfX9rhiM05hPrpyeI9dCGjLQkXtk9OY8atGQR9SiIfAsY30pDeum7vNDSxKdiMX2QKcrgl4vOyb094Qi67XMSrLZGM+SOxQwGzKr+6df3acy6WGr8+9EjGFua8ufiBCDV6iTHvj3TLmPGc1UvN7yHH5x+hsKwSbCsSTbO+M7iudJ2IU73VD8LipWVCsdm1Ec7HxR/T43a516XmfkLHUkG4POJC+R1fiSOO9/jvzoEPJtV5ECvFtdPx6Re6WN65NsdLMKq53h/bzb6JpaHKHEKi+xtDMk1pCg9/+eJFjTHStpdAHFbYcyoQk3LHnTtsZLvpP0odxxO3nQ8B2uMpIfix/k9Y0wGq7k3/2HYNKGk0OOdjMzvDaEmFbLwGlT2I944jSvIGpJSCFvmZO9TsX+457XqdTB/CKPZP89Uz6p0MsdbcMfar6PJ+t+TdLb2X3F2MUdbcdTdb/Br6u+DXlD+C2l71HZqAwFSZZ0nplgje8gybB0CpewQhomHJVsaA0bUBa1S7JenfT+zTTKqUSDO/iJlH38aJo2riJMc32E99LhlJtIGDmgjGIkEX3nmWbTORMsE6MCw19/Ybi0TGEkKWS2vZPXhdW5gTKhhiVe0RR4SXGEAst6uAyyHsMIHa4pHSKR5/ldC8PgGp4YjCCiqGxmgTmmJfDHTpGg6uFhOJG3pyVNFpGfi7OuJ7+wM6X/Bz1O1iod6PVipBasLbgZ75/2OXK2qdjR9iIervomOcNe8sHRcy9JAvZ3bMHPT/0rE4dnaBDT8KFpX1SzKiQ898jBzo0sm46h8Xqln1DSzSy5kiGlnFkVwyDL3R7apoZs8i1rmZ1N1MxRhi3DLEzP+hgfocfJrh8xLO5XKD+ciKFnGSepEY6G4F62xuCEyCovGduXoj3M5Mx/jE48+J66T/Wb6Dd+EaUF4nVMk3+poLjC+U9MUUdau1RDEn4WwfAfGQJJCo1CVnukv3Ll/lZDGTrIQYR7hRIm+p8PJfbruq+QkXgbSb8Rp+hJ0sc11SHEOLNhu5gSz3etYOM6cdJziAawTaGP9MFIOFAhknYpoaBvL7R4sqyN+mT9r/Fq81MqDb+59L24qeQuNkiuUryIZIINwVMqNM9wrsC8rPVyM4alszjheQqTmOEKcY8m3IoLGXQuhsC7eMn4QmBmkfnr+azjVIa1g2gOyls5WtZ9Jmsk3CizAcv4rOzhLYiVb57N7/fvLBV9WOjAJ7zblK5mOAfPj5tAw9IwC3wEwUQ9pjjer4xgZDEgEv0ryewhWI1rYUgdhzb2DOs+iTccuJlBihW2K88OxE7Tk+z0fg8q7KvYuD1jbLLZZpEKP3XBo5hqX6K8K7PIIh5aIs0MkuoFyqgOde3AgY6t2N32inqTRiYOtoTqmV53kBP5GLL82Nr6LMPnr+j9LcxGr8btkz6KaY55vFv6X4+IgWxuegSeWCvWFf0TZHMGCUUn3E+r9yHnZr2TIT0HraEdzBA3odh2I1FMNkTK3NjjlzShdxnnEe1r0BZ6nfqr5d+zlT4zP090bSYava7GjktsC3mm/3UyLtgSqma2fpp1kdfu+4PIhBiWWH8XM8Ba31+RZZpPsn2D8vTMIucFOplxkCMFwo9Am/LBljoBTYLcI1nH77qR0st6UP29QHxZSLx4X5L1DPFjl8GEHLW+QTqMyZsidjbYCY+ElyiN5gp+lhm1RERhglpzXMtQYpmkQqWU3RvrZGZUR/Q7jAOdMl9+J41uG07z7xLrFLyr4mM0rGuoYAvv3v/+abQ6qdBqunO5ehlYniQv3h7t/IsaNZjpupnoF4VsVSJS4ZChrvGvKDO8yM5fLuppORIsQ0d4jwqLNn0ZLAPQSET0LMtQtUVO0gibFboO3DhTIoKsqXbGt1dxx1LLzH56mCDEYlbk/Q1TZjemOj/Ixh9qvQEdDacamujv1RGNPKFQx6LNhk6/BHHj7TSy87ymnYZ1JWso0zP6VlrmDjmIHM0MhfWIpqzQpEIosws36zVkGSxtD9UyaznMLHA2uYa8gTK8CNHOM5cQ1hdgYc4qLM1dq1YOFIOTl0rzTMU0CCtW5K/HTcV3qSkiAxukR2Q4ZWvrnyHbhawtkKWdsmmsOnRGGQbd/yDKXk0+NYcoeQz1/n+Qj17F+8lyTUM7wMWLdNQa4DLN63bOQyzfzvQ5A0OjCsG99ZHZDt7Yeci2KNKJayNZH9gWVjpIjwPPdkkb9MpFk3dBq/aQbBh+hA2zisofaqszBkuiUSr0NaSkzyoZQISP14Bk1fJJxE33kMcvQVLHzJCGpUme4lcy9ZGk2CiriXkaVp6ZCxsvFG/m+V5uomPIWZ73VjZTArvan+JHmQ1gsKSvk7AuSybmMazKYrrzs1bi6sK34h2T7sHK3PUMLWZeOYRRsczNwbM45dmJyQwhxZbeWayNAXmzWPa7nq1IcqN/o0KGXLPsknopjapH0ps6FVjWYWbWfUTnLFST3532/FJ1zqZ5V1pkxEC2ZxZtSJ+W1KuvSJ1sNCzZ/0ioR1ukjmd72+CiDUtWiKnz/41I4iRPeAtv3Xvz/qJDKnGY4a8dSdO7ETN9gpUhpOrnkCNVdF+TQAhZCDNTTMW204vkzRVZDEReOJWiykvgSXrKJGSZ5xKeQ4gmNaj3vUxo7p18JtfI7qBznGtV5+Qx9+sKkcYmojreiQYghF1+l74p6X2Wf0OJlPNQ10usWwSLcpgJdpNxeatbDEvmlOeYpsMXrVIdyNJ5bO+zwO/lEHmWgwR+dvZnkUcnbQ6+gqMd/04Odo6fduuJ9ZZNT2VwvJUkPqqmVvdvW0HhqY6lRKwYE6BtbI/e7oyLNqzGwAtM7esV+Rx+8VVmWYkTLA0tX78SkfheNloMBoPswC7IpYM3shnNgcdxPm5AXfgIGj1fR2foz/BHdrJi58kP0l0I0niFlqvpcQlWDjgfqkKSIXWgLM65kdfYsK/jGaJa/zlGl0Lk/jI79KRnu0ocSoR30CiFc7WFTjG0NKPYulgpvSW0hZqKsx7XKOO93CLLRwmtmOb8sOqQDcTP4ZzvUZYp3bksziOEvNi8CLKhk2xUmnbuXpG6lVhnqN74Ku9uJm7UcbfxXYRhyWYDrTSsjcxu5NWgtcqbh5Zuw9KWq+vSPe3lPKazMLKYWzU8oWeIXrn06rkwEbviiSp4w6+g1f8jNHi/hibft9Aa+Am936+eKct7yzpSndE6Vrw/FIti7ORaS3JuQnu4Hge6pD/pUhoW1U6n2dfxLFE0DNmwPb1wLIWGcz6wR36hIc2nIzbRsLazntNZxmk8L313NEuFqpfW+PuKwl6GPAGFPPOV6AwfoAMI7+pBHo0K24L2TYFD/GswNZGtBue7rkZ7tAEbG39B9JKdR3jlF7/2mQ19pz2M/tCjxvsYG/QkM5x/ZjY4leEoPcVi8GGAPtWCRPgx6AzLkKAnBCObYDPdCKtRPDiGFv9vaTBuFNr+CU5mlq7kdpUS26z3MnOZCgNTZhlLDEYPEHI1yLHIWzsOKkKm6uphJ++pcF6juhDkdXB1MJTmGEtQ7T+IxmAVifgqhp78C9NbJvYwoYPK3dT8G3rxTFxZIAvSyfJAJoSZ1Oxp+1/kmqYxMXgXCfvTJM+HMdV1F/U2XelIkwoilqgjR2X51HQVbR/9XdpD2tNhmMRscR/DYTVKbetZDqf6TN7nrA+IjjsxN+sWpVdpzws65pFvLifiuXHSuwO+WBtmOa6AVoY+xnpIqumL1aAp+CobfwbJpxiHrFky1D8jUvH99I84s7/ZCEdfpRc4YDHKctt6dIWeZ6g6Td50HWyGGdBoi0nJpkDHDFK28c8yX4V8210odX6KSDUdnsjr5DCN6j24Qou8hxdDffAIXTDExmTFmenIIfd2kCOsKXiXgundHU+zDLy1+mxiDwmD0ssuPf6LZASAIVjKIHPSav3b2TAeTOf5JALU22bWs/yC3hLJLtQy7Ne4v4ha99cY+ndRZ9QziX1fLV6qfyqF0hej3HEr/NFaNYFAyi0DzCb15pH089US+U+q8z36TetYB6vOhRuK72GSsw5H3FvwfMPPeM9xiHCjM54/K2ic7rqbNx/pHUQS4HglfxppWgaEY6dgMsyngZaSV+1FR3AjCexcZFuuZcNLWKAh6uYzgnhI+MUgpZj8n42Va3u76tZoD/ydkGtGmW0dLDotK97MsCzvxg2UFDOXpczQ5qOyazNDEsPxJZDG0Gkc7pKVY2aonnaZGybGFky4mSFuVMMk5dYVaAi8prZ/KSLiGslxEkSqes93EYgeJ1Geg1DsDKppXGc7v8RGPsw7SLOPPMZ3sSJ6L7FeRyeegxrfk0wsTqdDMz2xwnGl/MD+jj9R9zI22F/LElKFy95Q/BEsIgU45H557IYlxtQS3MkQtIfxVxbsmMyzw3ErWQrSgyT5kk5fTsXt498xhsCrEUk2o4VkXcuMLs96O1Gsp4MwQoK/hF9lWh87wJ/pYsr3bDRIG42wK7yZDVQF2TKl2DIL0VQYNf6d3Tylv8iEwKsK7mJJNHi56TfKCCdSOiNN+Ef9D5hghLGW6CiILiJcpZZlkolz05xXK++u9z1PFC4mr7laOWi99/vwRHcSke/AJNfXUJH1LWRZrmdoOYaqzgeJYv/KDHK/ul8mjjORIuWe4ryTThFBlfuP1Ldkw7IO2SzMyboVDcEDdM4nWK/B5VBkn+2VNq7rx25Y0t9x1vMYZN51heN2VnYUt0h2EnkaFXEPR4+RL5VDz1DXEniUvKIduda3w6iXSXGCVuoLJJW5jFllacLPJKFXdORXt7BRQmgn2Rcpt6+hx2jV6/exZECd6y8ygW0KVnRPCBQjcEf73nP84iGneLL+v9EcOosbSu7BNKbfapyRmokmQzjhfp5hJpe6WkfCvotIcJI88ipyEytDzi+Z9b7EJOQGGpZMlUky9JSi1H6vMjCXaRXD/k6c6WKI9P4nArET3fq+NARfyp3LhEgWLW4hiW8OblPGLGg2N+s2EvwZONL1JOt6lKXIZFyyIboB1zMsjsmwpD+pMfAqPewM4+7VDEHypsZIoiVaHVAWLS+mxhnGrMbl6oVVv7yYalpJHrSYRtV3ZoEozsBwuIi19fE4xb/TFZFKWo1zSPyXMGPcTrJ5jGFmOr2qlGjoRRXDTiaPElmV904127TafwBP1/8P0//27k/GJ6GEDy+QT9T4Disetzj7BjZOuutAdZSGKtlAxzDZsUrxkFrvs5D9m4usa5jdPol2/1/hMC5DoeMe1YBpv1e795A8V5BT3o/JWV+njlbAHXoVZ2lgbaGneJk849IYl3TcVjjeocJ0teevdNT0jrZGnQ2Lc96jsr697b9lW2begFwyfulcHrVhiad0hA6x4f7IUFSGcttblBJGFh3isb28gQGx7u6IGDlUe+hZFqCYvOomfpapGDQ03Sz1vVRMdmbv5RkS7nLMN5GfhNEefFadm2xntpkKqGWBJGRmEglFV+W/G0tz3oKzvgP4e913x21ckUSQRvVzktXXsCLvNobAd18wKhFJuU95XmId7Zhqv4YIeYohsxJ5lmWIxavQ6PsFUXoSShz3slzpyX79JW08Vob+UudnUO76EtHAgfPeH+Oc99ts8A5+mtmBLkakHLI8ern9Fpb5NM75RL8yU1emUc/DTOcNTD6OKOS60J0yQMQ4dQ8+9KkRxwrFqLzMAg+0fYsPiGBu9idoXKPZHYwWnfIhFv4tfzV0z5EyQjYwTySDKLK9j8od6j5UrCaHYVAGpuv59WX8u3fSoEGXjTjDqC+6B3bjQqbFk5m1HGLK3wqdNk7eJQPTgw1fjGuSbZ4aQD3meV2FsKkO2Zls8Na1Q4n0U21q/i32dDyDBVlXq4yob2gQtGoPn1FdDKW2peRXa1HZ/l1ySg8J8nwa1c9oJC6UZ32Rzy1hOXsoQCZJL5Fu0pfBYbpCIb6b/NLLzFF0Z9aNxHHHLoKbNmMZeewxdERkMZDljE6kJjyfQ0ohqw/WBXYzI58Dl7H/WmU9MgrD0qlpFvvavkEv8dCo7qPXLeXNhlNGj9DSk2cRDz+OpG4aQ4dsWh1i6h0mWb8VDrMsMDG4UGmhUbKBNEQhTfwgf53CoshamOnrxVsMugKGw61Ej2be7xZoU1Volq16g5VUeAL5Q9xfuiEm2efTSeI46t5KD6wiqixWU0FGEnkxYkvLH7Gt9XF+ZxFuK32A5eg/dij3P9T5OI32CNHsg2qgWfhKoXkes7+XFeKWOz+vuk5Gp0eRFB3GSuNaSX7mokMdQlfoBSQ1USZQ82jMMmNk4gxMNiw16XPQFNhCR+oiCMgIiZyXtTHKcY6RoS18ik4pSUmGiX7DGZYglWyidLDj3/mzBTNc7ye3uooNIqFmNDFeh2R8FxKx1xBlKPPFfQikTMjRy0AysyJNn3UvM4oYlxmaBEMpkU9j6DuVJsUEwEXUcsMX2aNSdacxB4Hwa/Alc1jpvST0KXrbUnXtQJH+o0k2WZk4od44EeSa4lgyrHHJtXvbn8XLzb9BBb/7trLPkAvZeL7XeAVdZBGzbS0/Rq55OhvBijOkD/LihAFVRGofQ9sDRNlFLFUPrxRd9j2GknQ9rIbZ6vvh+DlF/kPxM+Sdsu+ji9/uPzFxvCIoZNHLC78daAptUdm3bGaqdhIx5KknnPW+JjOtqUd59b+/Aw9jWBoqIYzD7f/JWHsCk53v4A1u5Q3SGc/oRIdY6H/JNevhJUp5EjpCPzmHzg8T3IjoZfGJ4XgCFanNomEdZ6s20LBW8dF9p+5qSSqL4Y1uJ2qdZ7b4VqTiT/E5NoYdK9P8XWx0mfK7fEhlT7LNZ41iOOrZqoyr1DqThuFl6JfJfV0XDl+8Ayc82/Fi4y+QZyrDOyY9SC4yeP8g6aA95n6GXG87DXUZGnxPqWkqecZsJisHUUyinu6vS0/5FZ6UJCL4Qo8gFNkIg34G0UcmJw7ncJJ9ZcOlhtES8ISZUEW2Qks0EX4pr95JbWUpAzWIr/6SY7h7DhZBVjGu9uBe6qMaxbaeTc6BHGMF2iNnGBJ3qk5qeZO+ry6GMCyZRZBAZcf3qeztzGrexjT6rn5fHJ3EyK9+yYbyoC3OAhF9Cqx3w6YNw5A4gaQ2DwmZLTpshYXYh6CJ7WNpy6DRy9pdPeFDZjo6iVpsmMhemGXoh6FBlzqNzkQ+jUrPxtyPOPlUAcNiZrKpoXHNVbzpGI1LjgMdL2Bvx3M40Pkijxewn8e+zmdV2JT1T2+f9AVkG2Tnq8H6kI7R3eRWbGGG5jo1FDLdeScbnxmhoZyG9TF+rwepDIjHz8Id+C+EI9tZj1o6cxMspnX8bCTnlR1pdQyNSxmyiuGPMDTSwNqDz5GDbVKjE9I/FoydYN1aqYsgdZWljG30BiazT7NUktQU2EqkKkK2aQ7PJpQDOQ3FzLBfR1e0lrTgat67V79DItbprt+gzv8MSmzXqrHAsQs9MXEEofBTaKRRyT7yeebV5ANLEdcUwpQ8BD05V5SoldKMtFM+wyFDquztB5kAqBo0rXj5X1DLF91BDtiJLPOVSDL0arVT0BVndsLKuiM0LiYLhZaVTEAHZ6Di0RLaZJ9EmWOUZy7ntRXqZdjJ9gWYQv41mZ9PY6hcWyAbfZaxBIO5kXSICqk92vU0ZAPLOJFaxgPN2ghD1ssosL0HFoYx1bDMGkORZ2lUP+R1LbBZ3sLQXohIbI+qk8k4HP/sK7JSzFQ4TVeScM9lmJzBkJhDAw0hwlAZiB2moW1jRv8CdbSPzlehDHG0xiW6ka3+2kktfLFzTEZkKpD0bcnM23yFmGe9m6lXDTPJlepvkX6GlfZoDao9j+O05w9U4BLMz/5ktyWOFUbpjQyD58Kn6cUGolSK0C0LtzqIVE7+TMAkGZ8mjphO3iUcSqQR7MwMzxKoaqAx0pv7DSEJ13LSIyNEBZmjTaWlGuAyWOBPTqW3JuCgsTSTJ8SYCOSZlygCOhBtZAS/3DqXxrOMx/ILxxQSdMkiBdXkp7wFlO4AHSjpuWI7237B0NlEftdBlGRGSMM67/0B62kkWn2Y+Gtgo7fBF/wlw99jdAA7XLb7iVK3MAzOJoIdRzi6i79Ph04naD4a46IDaW002goa1wLYTcvpYFeTGtyqevHFmU0MaYLqnUL41QjGbBU6R75/etauZPEtQZmRMQkuk0xeTH8v18QsMXyCZH4XnUl21JCd8JloPPjQAxvSvblJWnQ1qr2Pq3HALOMcLMp7kAU2qQvHJiSQ0T+gxvcXhFJmGlVQvdRgMqRXBCZ7I2qVwpA6C2OChseML6kZeqpvug+L5D1+iKGQhF8rr5X17VAV1CpkCr4D0ZSbz3KwOq3Isd2O+uAxevQkZoiTiMDydvZpKmcyzwkBlXr3PlNQSBQmUJ8+ZH3S9NqkPcdQuhDnawkfZxh9DEZNmIbtxMLcL8AfO4DWwOPItd5GVFlNRNoLt/+/+PMQzMYlcNo/w6xSskNZidkEPdEnEtvFYz/RZTkbP9NibJlESiZlFr2I4ac3b9eRk5oMZST7y3ksJtGvRmf4Jbb1YSL0JDXiMdL95c4WXtca2gnZVrDEdjXvnQ57sjJNnmkaHfcoTntf4ZUJotoS6D79lQ9uaAnuwkn3r3HK/Vsq/jjj6GzMy72fjcVKkSsIHI7+nxa62Iuo8/0EbpJ1J8OAkehkJvEWYpq+RvzbQOTKgil+GIZkM4m8LAWt6/584D9WXJNDw9rO34NELdkarW+5WEE1P16rjEt4lyZ5Hvk2WaM+F42hSlQ4blYT0mSH/RoS6i7WU5aithoKFV9QDXER/6SLYXf7r3jfGth0ESzI/TgNeArqPD+gcfpR4vgoggx9nsBP+awgHJY7GP7uoUFKuYUoyD9pqFweTqLWNqJEteJbMgjd85yx/xMRdIkRebJV4iCo6Y3uZniU9zll+5NZCkBE+n7zwj9eJDuyyk4WMtQjyyQ5jb3boMj6H2W2peggmT/jfVWd162859iGOt+z5Cd+WtoazM7+Z8xyfUBdnJ6r0zvvZqRDTbVIHMR5zzfREk/BQaNy6Y0MhQxH5vUMI0VsAFkuSZbXYYH5t1bTTkM8Rk91Qaufo85nWmJHr2MDJKuQjNfCbLqWXi5LJcrUkvQh88Ms5A7u8E7V72JgiDUTmic5P0zes1fB9VWFX2I4W6c6R2Xcrsb7d3hip8mBHMg1z2S4sNOZZNd885gOWbpSOg33tP2OSg1gQc7bMY1ZdCB2hGHwYeRZbiSKatDp/wmMNOQ8+5dgIy1I17W3DurQ6uiERPBUO4Ik3wbqymlaw58GHtTzRRwyzUV+ZpmX0zAWIBivQnvoRYTjpwgmSxSXkrlZMhg9uH1NirzL7qzSZTLVebOqt3wm44OyHc1k2wq0Us9iXLq7PzNvwzTnO7Eo91Mot93AeCw94Wk7v2Cxo/rH8JfyoaHrPjREfbBSQUX6MCIpQibTZ4fiV0x/eesLh4QabRlJfCW0iWokCf3yZo4s+NXvOnWtpM5hJOJ7aWQFNMKF/FumJvPZcgFFz+/qyTVk0Vz5TAs3ihwfgZ1hr9q3neSzmfW8i/VcjxL7GoYZIxoD26kIJhjB9MuXsraU7BgrqXZPn1DmfzKs40Ez0f5416N4rfmniJCHTLXPx+K8e6U0yqgC8dMoc91Dov4CSXoXilzfonGJt8cvlL3fof5pGb5mMxweo3Ftp8GX8ZhBJaTD28UegigmhjYxePmrM7wZrcEXyM9kB/4pPCP67v9PRDLESNKN88HNDH/zaPBp1LrwOflihX21ms2h+/F/PrMh37KMSpUFNiRGjyaeDxYJbe3Bp1BDi5ZldEqNTj7Ti2CSoZBE1GqUN1EGS4poIT3senIRDSJIynSZjGUgnEu2E5OxwCAMput5rj/PUgojCnqi+6kAH70zRF6xFDmmZQgk2hgCd0C2bCu0zOJ1uSi1XoVJdtls0sqs5zCqPH8jv3ySWc6TaArugGxPLJtVyQsjPVu8BeKNaAi8jgPtP+DxfdKHv6E6UKVC8fysa7Ak70P8XTZLP4dq93fYWLPJSd6JNt/DNJYpcJrfwVJmHsvsKxrez8hGDka207j2MCRKj3vm5RnHJ2IQWqLXapZxFrqY+HiiB5BvfQvLn/kNdnE2CYmNgS1EuxZFL/qK2I7Mk5/iWAPdFx765IbemQVpyxurSAH90eM40fVv9PYAw8AHiUAHEWbjRlNa8onriCTp5XoGiiBLUua+J2sYEo/SsMppbL27I/QVDRssxewwET8JHTMdTYb9ksVBjLp8tXZ8iiFPZqDaTGuZrUxDQ7BSIZeQ8jzzVIVHgnAlpABiYHZDCRVjp1F6SLrraEBbmR3/BVU0tLPeZ3DW9w+i02957q+q+0LGK70JeWFVh1X5H2QIlHVc02Vo9P+FYeZ1TMn6NBJMUHzhV2hU18HEENTbDzecyBBOLsuXg0B0KyLxE8z2BPV75qxNjAjZt+gnq6bvCG5i2S00tiv4SaZsMb0IbijeSqTeCdkN32kcuA2KvGGuE8O6f8ghndGJBrFEB453fo6xtwFTTQ5YTStIPl+mYdFoWFA70UWresyHEj2NKRv6+AFok01IGGQ90cHjTyo7ZDhIkHvoVDiUacmDDdDIz2TuUjjRDqO8UGZ5Jz3NSuifjrbwWZzxbSFCVatU2UJuJ2gk5DXPsgiTHddjputOTHHegjL7epLSq8i/5vOe0s0QU7/Py/kQOcadROczRLUWLGN4ne26iciWXk1Y5opVdX6T6X8hyh33wBt6CtF4HXL4u2RpmcqcWYgAdDSkIghEttGp/NTt6u7PJk7EMCxEUzX+SI6aa7mK+koPOg8UARGjLgfn/a/S2LswyXEDz4ph9QclhVjdv49LpFf2VOdX6OFHMckQQY5xGoHejHjsMAJ8noEprc0oyxNlmhrTI0zxWRFNqouodYS/FyGlk3cNB1ZMOIad4XAff22D3ihrb2WqvAwdFaAjvI2N7SG830o+5VT9LFMcq+gIQdQF95FbbSGiuRgqZQsT2WZFdqCIqUM6O81UoIMemWdegFL71TSmt6GcxiZ9Rq80fg/tkbNYnHMHFmS/jd+RIRrehTyyLfA8WgJPodj+LmaG09AR+AV5YRZ55m28pLd3enRCdmmYx0Y8Dn/kdT7bQnIv25KM1jhHJ0wZiMAuZoovE429yLNey7OZniELjRCpo7KhwD4UW68k6S/g2f7XDtfaI4r0stZ5fwZ3ZBeKTBUokF5tZnYxwnYsJW/yp2jdk0cH3yxXQk+k4rW65PHek/1EDCubSDWd4bAByYRMAMxUhaTqgXYaFzOsJamsvytjk8rLCwor8z+ANQX3MuxZSLp/jM1N30eIXEqMa7DIt1gT1lV++ohQL5z/JjqjtVic/U4syHo7jUqoRNpjxdFags/QkHLp+TcizKwzSj4ic/zTY4BjF+luyLd/jpyrHO3+X8IbeY5nx2qgw4uExCyiYbb5KobwV5hdMwGic2US0dN0tYlXjLz0Lxn1dlGG1Rj4E7nEn5kdXIFyhgipbEo7B6lEAw1LCmukd0lnZibLHygJJJghpqSXPSFz0uU7mTmf3iCDr8wQBbmGuEaSiSIbK8/btIV2E7n6r5c12b4CN5V+FRW2lUSurXi69kuoC8i4Yqxfh2jPIR2l/lgbjeob6KBRLc25CwtpWIJyvaJVL0V4IvtVOJHlBgLRbQqtLYbFfPxouFVm0ZIq5Dk+p3hXh++HDI2b5Gz6wwkSMZAS+90srwW13p8gmeHFCRFxNScz20LzCtQzWeuKnFLf7SvjCIXyIB0zwJdQ4/mOGkaYnvUV6COP8yOqWb8UEfKrYDKmlGAzrWdBM/GlTMLsMFlFntXC+8jEvkzL7AhqGZGMHyTnOM/sUMLhEAbIkBUiesrr7E6jbGLeN7xKBmMhp5K9EPNI7A/hpPcV/jxIsr4Zp72v9juqfK+hsuspeGJNWJH7HszPTq8N0VcEFavd32PGVI8prgdU31tn4BGWNwGX9QODlD9W0dG4TPppzBJ3MFl6nYY7m+FrqLUyxiPprFo6dNtJ5I36fOpNAGMwMKSTJBfq/S9TD3FFEfo62ZhrqkvK0kPPMS3/NzZMLqa5vgIjITMVP8u7MUQlz7FgET5Cw4LJFJCBK8YMJ8wRteW0Ty+0qRaiVya45734XK1+AUNhKw3sKM9lgmxZXcVCgr5K7koP/LnqNuhvhPQ9huzpjnW4ofTLKLMuUlOOZR2sQUfckw6jee/DvKy3KhTrK2I0gXg1acF2ZlXLYDXOQCRWRW50msYwR3Gvi5ckQ+pcItdn+bw4mr3fZHmreD6zY41HUuSKhVZZu6sUDd7fsO5NPDvYTMSIJJEpsC5j9ryFqNW/HGNCLG3iGDSee3HUt5HVMtIrPw+7cR408W1AdDNSpuvIKdKbCEVTepL21fQoWYBttIYlHQAR6GMHkaRnpmiomVNzjULBROx1/jRDZxwqPRYvNzDVf40cyg2bfgrLNINn+5dHoF3miU22r8RM5/qhD9e1alysf3qdFuFBjb4/MAzuYib4YTVMEohuUutOuCzvUEgzVBnHJjLoXqSOYHQLgpG9zLqvYqgcefbr6ES6OUhHqN/20Ev8O0He1X+Joh5JzxzVqt54QetS2zp+O91eo0Qsfi1Rg4T3M6iMdCFA4jLVXEiyt0h9lorJtiPykmkxErJuldKfdPANbsThhU2szeFhJR0RDxiKkySpyFJotIU0roP8U5aNziQyX6uImV0hSxkjbP8vM55Mr4elRfUxS4/7UMcQ6pIBWdmtrC30HMPTZLhMy+n5MTa8LFnp4DmZDTBUXcYnsgJilvVDiCWb0OJ7iFFiqI2Txi7Sr5lrvobceSlag/8gtzzBOg6OCpK0yFQkFzPf84HNRGtyLepJZETDUspMNqLN8xnsD7XDnXSgwpSPKTgGU+QxlsIPxGkEunw2GkNgwkc0Az20kPA/lj4bEWaSmlyGQBd5luzYJd/N9H2eY4YlnaTJRBvD70n+nTlLEs+zG1fBpIkgFKunop5R/GCiJG1UbTjZ+TmGy0YU2N5OT7Yz7DYjFD2jMjm9fiyoPXpxWt6OLKJhKHaKxrWBBhFie8mMUTEC0cd4Q6R0ctpI5N+rkpZaz0/5U4j8YJFNM2dm3c26t2NPm4RmD59KIBouFAq8yfJB1R2fRG3kvILfqc6PIc92JwyJPdDHdyKhK0Ay8jIBi6SdITCcOI0IdWgzXUOlpsedxiYW8rhTDLsk8DQIkCcNJRqS70RsJ39LwGCUjQf68560CNowwDJsRpNJpupHldFbDRJmL05EP0Fyy2Pt99OrT6LM8c8oIj/Rak0MUTvhCb/AMHgb9TCTV0+8YYkjWgxLyPfaGHK30JD3sA1qEKWjy3oQkuSoqcS0rzTayjHacpDPEe1jqS50hF4l55pKB5VJigPDuUz4K1MTKOv9r6ppNWX264Y2LClIV3gHjnd8Dj42cp5pLqa6HiR3mMpbmUiyS2CMb4EmtosPjyJlvJ4V24RwMogEjcFBw5IsZqwi1dZJR2lChndmMizKZLchwgiRIRWvZFrcAD1DgwxyZxJBzlBUFtoIIk40lBF9qZ/dOJffydxXM5LI92VTyuMdn2JDtmNy9qdRpNZol8bTwR36K53yHLJs7xmXHkYvKVjo1El5GZg0JBQ7gGBsHwI0NE/oHzyeYN130LFOMhLEVRYpDjEqoUNaiLhdkR1q/la+5UYimTj6QOPUINs4i0bdifP+V5Qj677y9c+qPaHTh6wmY1Dxus7zS1R7vy/mhRm25Shyfp6FElJH1fEA0UtLg4ozK0tqshDXT0ckuheBZJwemweX+QbCpHhM+vrRHjKtRmY96KSPShBPP5/nZPEy7aAjvRhJCLHoPhgN82DQV7C8cp+B15oYphsQozeXOj5GVO2gcb1M+G5AjmUVs1vZqVSu69HDSIcebcEXcKrra9RxAjNyNlDpN/B8+tnQJNHm+wkRPg85VhqW+s7AMk3kYVDDaA7q3GV5G6PFcqJyBZ9r5mc61r2LdT3CJOZlhOKyLVwWUbSYBm/t/v7A+vUcWra5OIWW9WVUooPnW2Vuu9Sz/3WSJOVbFqMzfBR1gVehqfc9fcH8ZBUX2bRSNv7uiuxDljaK2VY2mPM/+WmGeJ1sQcRzHw0sjog2Hx3Rc/Dzblm07GyGAD6++8KxiSblBvzfhsawkiz1EzzTfxH7XpH3Fs/B4/kPmEg0nY7/Is/ITM5DsSOo6/o8Ch33Uvl34oz7u4T4LaprYH7u95RnjoZgC6eq9fwvv0+n0zgwN+8bcPWZny5oEKKznWm7B9nWm1Hk+Awd9dLuODa0SL6bYFj08vCoFzrcspJgykfDm4EC+z1wmq9lPeRN7Ew0QoRunorgaPtnSSOOYWH+D1Bou5nnMs3QkA1Az2Fzw/3QXfPh5zY0BZ5BU+BpNSenk+EvImu168OYbZkBneObhNBMIYbZIPlXPPIi88EgvAk3PElJVfPJwd5NbxnpBYlhhM/T0MM0qU5mAcO9sSKdpS6iZiUS8TqYzdIZO9SUD+l2eJ0K7kQeSansQS0e2BbcyGxuEzOgeeQrsit+j4EoCqo8Mp3pyOxWEKW+hbOeH6l5S/Py/pMZZ/+BcC2f3xX8GzyRPShwfIDoMHC14Uwi95e7i2FLmJLwPFSdxyJSB0EkC8uRA5vxCqLaFayPno52HF3MYn2R13hdRBlaeuhtcFmFp9mNM6mnV1T4z7esJ8pnXhbUrJMVGcuge+jrX9lQbLsN+bbrCXM3okhvRAUOoNBYgYTj60hqM49yS6GTse1IMqUGYbc97lGrIOdYbyUxlle0LkKoCA0JqIbGAtPV/Hu4zkVp+BgiJMt6GoZel7mLQ5QbZSgMRA+oISiTrhwuopwsOdkafgXNdC5RlpOGIh4q3Qf+eDW85BadoW0MnVtQzeyoyf8kU/HVmJf/HdWJOPBZMne9xf8z2meUBvwBPne4DJQNT8OV3TkCgZ8iGPqNGrxPJqVTMkxK4eD3pYNZjGx86D9QdMrAlqkXW2RAOxCtJA97lQb2Mh1rOkOfONdA/cl0mULqx0EAeoZUQqZ9yyTBzB3TThnU7wqfSAlgqkoyjBi8H6NSziFm/x6Sal2AocJDAjH/v9K4ziCgLUYzY6uFSiqll6rhmIsRQazoi0DwaWgcD9Jwpb9sqHLQwJPNDIffgNE4h+Hwv3kuM0yHGaJqO7+olk0qcX2NtU7ybPr1sNOd/w7ZY9ksfXEMXWJcMmOh56ek37KTgyDd9KwH2eBi7AMbQKYQteJ06+2wm5ahLGuD+m5mkRWYfQiHn0GIRD9FrqhlopJi4pIkJZFpQ7JumLypYzQsVVvC6GnIWiKgjBYMfvZYRYxVdifrYHj8B7oCT9Ep8zEl71c0IHnBYmD2J09MoKrru2imc83I/hKmZt3H+g2xAMsXH/o4s0IpJGE/cRC64MNImm5EQo3BDcM5qJR48Fcsmx2diTBTeR/KiHbOZA1iJNwp1Yc1flEvUMg7dhJK1DSRoTgAr9U46fWneFQxHMqisuLlg0XHsBmIHVTDLFmWW2gcEjblDZQSZJtXMjzQoVgvmQHrMi4kT1yDYtvbUE5nKXXciWIe+RaZW5bJqKTMBvKQF9hQryjjtRjm8qqB5ZYGlcV998Hn/08i7avQkuRbrffAYnkHnWM1j2U0qGwaEA2cThON7uZ1L/HYjCQzdJ2uRKGZtFmmcoxeZNjLCllWSqaDe5lJirE7zbKUwWARhJK9qCUctoc2E/WWMTplfkWNhvWJ7u4GHQyBH/JZtYhbPso0f7jXsYgS8WNIhF9CnKjWHq2BkWlsPhvClDhCX4zQuIR7jJcn8LlaZjTRbcpANQZ5SXVow0orOEzFy0rMecwOM08AFA4h5F5eR7cYpsF8IWSnGAKyFBUocdyBQustNIy1yGKotDE8SL+XQScvl0hWPHSqLnysxfdLhtF28qsPq0brXw49DaVVhb1A6Hf09gDMNHCb5YPQkbOpV+P4Hcmq9frZTEjW0MhWwmBYwHqV0XGamQHvV4P8icQZ3t/B84IuEpIGN+7oRULddCL6SXjDO6mXqTykz2ogsMhIhk0lOm3Mqn2xY0pXmcZBuwM3k/RkLTREiBQ5SlKFnuEKSjyJy1qecbTH6vm4hJrOGjVcQU7GRmABtdIjP27DEiGayMJu0gNPLxpeEipUiBdHortUeTJLiiFqNZHLrgi2TL3pK2I04pXpvi0xoJ7yi3EM5WQ9kg6DYdlPkeRVXvjo1WEapcKRv8Ht+RRCRB4DUd1hfxAW851E5UzvDsp3GEWYDBkMi4jEt8HpfAg2+8dY1+l0ot3wer8Kj/ezrLOM6YnjCZ+TV8V6yt+TFIwkkgRZUWj/uEq6mn0/ZkLUwG9m4nUphe7ljn9mlnkIZ5kdZ5q90o1Yeugif4eWBC5OSE7JkkHDKjKKROhRhJJt6EoE1TRWeSET4nHkCvrESSIZDXRUGdFQwgwlcY7HWVIuZnBDhLceEcUk4ucZYvYrLzaoNR4yQLTMaY+TLMcOMEO6kkgkHn/xIlwtSKLfFngM2UxgbEaZXi0GLjNTgwj4v49A8Akaih0263tpKLeznKPVT/oaQVwdI4SBnEuMTc7G46fI07YwXL6GFEm/dBYnyJtk3dd0N42MPMjbUWYe8nv6SEuP0aUNXziWbAXoDr5KntlO/cgs0kyiYTY8hzo8jRZ5u0c/WWXVfQfnlWHJAK3e/2+83kDDupc/h4Z7JfKSRODXaEuGyUlSzJJWwNQTVmhcZpkrxYaOEc7HLSyLJtVG3Rzhrwuon5HWG9BS6QVU8AFEibzp8JGp4QQFjMyC5MVKKsh8lfp5cSIeb0BH4Lfpfatt76PBCmJJYtEIn0xvie6BybwKNtt9LJtMfkw35thFniV1yGWYlI3Zl7CeNhpYHWLxSj5nJ5HxNcXHJGSmudmL1MvrdLojdL4ziqcRP7sRW4xPXqMTtNYyBM5EVBYqCW8lBSiARc3HGuyg0gXhMM5FG42wI7wDedZrCDA5/CRdJxrW/Ru0sZ3Qhf+MpPEaknaZHD/c60n0wPgBeEOb4E4mYaT1uyw3sqLSQys3NcGQOMZQ2IaoXsb6MqWkoxMNMzLE9vIWRBVD5gr2iizX7eo2rt2Iswwm4zoqbGC/lvCEHBrWdoatKmTJfPiLTDREJ9IX1OT9Ie9lpmF9gH8bVELh821giKwmMX8LCfoHWB5B3uHqMVoRXUudHQyNc8jHVtDIltHYlvNYzHPlDLEOZb7JlJdlaaHxnaVxHaLBbaPxvUCke55OuIl/v0ad7VD3kizUrJ8Gf2Qbj/1EIlKHjF1OKZioR+GezYFn6VDVquO0Z4hMIZYu9CtmhFWI2z6rkGZ40SMe/gdaoscV0Dv1uenXwLvjbIqWrJeZnTI9mbxHXqMfn1D5LIuMRUrDaYySqQxH4NMiYVCjSVJZO4gWbTQuGZzu35BSVsm4/DQuIfCWnhVgxinivdI/1hb4k+ojczBzjMf3wev7lgpLNssdJOnv5JWi9PE/Z2gRFDPRMHJU/WUrGb1+BoyGxTSyK2A2rafhXc3fV/BYpBBTT+NRxkcdy+S+GDPlWGw3w+Fk9T0VEsOvIs7w6jSnV5gZKNJNZTVMYbLSxpBIZGSb51uvoz500H35K3ds0Pv/C0lmP3GzvBc3vDdpCJ1CfLviTTDTHWTE2yDzuS8IFUdFm4hqSVZUuNa4lcn7aHgfCb2KZylCOpKQ7ZDcCtyLcWm1BvItKV//esmurH6GDBlHyyKapFFnPKIhHwmg3v153suLIucDDDeHaVTfUQ1mt38YRpNwFbn/pTCqgSLP6HmO4JXoTBISMTwXj2Ia3hQalbTbfIVuJuOV/HuSohGyaIkkFhbjKhpbDbyR11U2LCvYZEqKJHy6TAsRTjTSuJ5T75fmW9dDq4s8z3L4GAZl9eKRwhaVSA/siJ5UxbUwDIoXDJS4thRJjZOodYbfYDgbt7AxpCeYz0Sy/56Ew4mGxN9me7/yvmDwj/REmVrT33Cke8RimEPjOsyQONTbPiOJIIURzd7/okJPItf+HtKEtFEJT3U4Pkv99My+vBxGNZT0GJsc4mBiIHLI72xzGp2BGb2dGWcqGWL4/jeGzhoUOD5JPZWgxf8bhKJHWddMowhCLVyYlfMQyuzvRUtgIw61foKoFd3ILDAXccM6XjPSYKkGXaFXESFpN2s0jK+59AAhyP1FFlJLsOH0iVryJOkqGJ1BDBbyOVkjSmZHMqyNvvGTVEIWOc37+bsefv//0CFq+Xvf72vJryQ8xYjAL7CEYy+j8Cl36CmS9mcYAlfBgk74A/9LNMxXXQkGGm7/9evfrCIGl2B5VzBjvZv8tB1e79eoEQ+N61/I0UJo9P43qUXmtd3l+7LEwDTXZzDJ9UG0h16jbmLHiVbrqSV5I2Z4iTKdl4XvjeQwZraRVlehIHawMGCqEBhTfVrjNix5XUoMV+I7s5m0AkYr0rc1k8b1XhpVF71QZj70fwVMXk0zkWfI1nUytDEWkRAQjlehyfM9Olge7Noo0fEZ9Uy740GG48n/R4yqV2Q6tZF8ymZ7L/VB4/J8Q02IzLG+ncnOAbT4fsCrhopqbBuqdqrrU6qPSxNosKaSrodVx+jQ5FgGWGKo6foqY+42ZGupVoZNp+1fGI9lJubABhe1+2AL/Bfi+gUIWz7Oc2Mxih5hSYl4cabrWsMk6B3f5N8jL6jRK2JECaLI7xEMPw2r+e1wOb7Icz1l0aLN/2s0+36Hqbn/jWyLLNgxmEcMFtFHFFXt/wR/+BCyDVnQ0TBNpmvgtH+UziZOOpr7vBlFBsb1CIb+AV/gdwp1RWcN3p/RuPaqtekLWEfhj5kl7bi6r355yQat9aM0FDEWQ8ZD9jHuDG1kWvkoHDQqi1anUlmH9U5FjsXI+h8y48AOXeI4tAxhKfWKvV2dH/vhQEp605OdMFrexr8zPW+4w8gsZw7h/YxaQFaW8LYYFqaHZ6gA2bumM/QMg6efpPOODN8ffMiwRqPn2+gIylqjZHQk7xZmXdmOzzAMyowEmaqS+bv/Fw7p7DXKOCdDYCi8hSGwFnnkX/7IEXjCm5mwzYGsd8paZvx+uo1Md1C9w6GJQP451Hl+BANRK1fvRDwlMwtl5/PhsjQt0Wou0huHn+8+Nz7RMmNJJd1ISU/8sM/MJDLQKut8fopGVYYu7/+gpete+ENPUnF+GtY0ptNXwBfep/q1er7T/6AJqueKIs1qHlMzkU6vSckqFTSqq5Dl/Bz1kR7U/r8v6Xo7bR+Ew3YHQpGDRP3foCzrszyvR3Xn59VgvppPn1Gocwi/GpCK9xUJgXXeH5HAdSBfz5RVZi6kUt3L8aRhL7PIomoSXg3QM40dv8KZeRlm8plBJGngmecAjSQkl7pi5Li+TnJ6s0qj28kfmjs/BG/wYaiVlumdHbIBElWS5hFiSFI/hj2ZysJnJ5P1dLLjONfFcJoKw86PrcqovsCrhGsO31Xzf0+0jEofoYHdTiqxE8nodlQw+4vLDnDtn0AkUafCZibRBKPHU6KkTCKpdKv/MVR3/QeydUkUmZfAm/CQDHuQ5/oaL8iUfvYVhk0/eRGVHrF/XWWLYxcdkaoaUc/XYLDeBb3tY6pRxyMKYYm2spZ6MPwigpEtiMWaSLSz0BV3i5ehwjybP6kSEu9kIr3zldo/sXvcrTGZQoBZk4PhzmleQ6P6Is++kUYlxi/O0JPxCrfrQdqJELl/Ep2ef0UkdgIF2T+CL3Ya9e7vQ7asmZX3B1IlG5/Wv/66rz704SHf0gnH63C26+swaPwoNtgRNyxAMLofNsvVjMGZp6b0F3Ikcixd4iQSvD7FNHw8Ii9wJNUOYjHoSJDHL1Je6YqQlYSXk8yvYxZXTP51XmWFkVREGV08xUPWj0i5kaRzJbUWlt0GL9Ntd6JTdQzLMj9ZJLWyUOzENeJYRAyJoYjJTSJ2AInwn6ij3awby6PNUnVMl+viyybzzPT6EjriZrWSUD65ZEoTRXvwReqpDdlW6WDu/xzdV7/2sSEMS0Ok+iZj6REU6mVNzIXwRg+qHmun9Z/S3j+ipK1dL9NxtOVIqgxyHEKkS8YrkYqfh95yG0+IUuXeFyMy8Gpm1jMPZi25UXQHAqkkgnQ8Pw8vDzcPmb3RFffy6EKQSClmlEfCbpWkhjrRqIFXQYmLLc8oRMKO8BqGbdmxNh5+DLHgI0iEnubfDTSsGiQir/LnTnLSRiZWYmAuqksihSDKeI2MbagrJF0IkszLHPkk8hyfJPCcYQLzHLkW4DLLoiC9vQqaYOTAoKdJVtDGwlZ3fgMuXZxoVQA/DSMSPYws+0doZIt5k9Gl5dJBag58hxaejajtK/w5HoJrRCL4a8QjW2HO+gF5tEzrmYh0nkGPCUHE+1X+HkfYeBP8iTYY1VqfWmUqUs+ELD4i24YQ1WTKsV2bjXjor2y4HBicG8g7Z/LCodLviRA6ElE7FTuq9n5MRnYgEZeRiCSfXQadzABlm8jiJ4noIXIhZtGsh3BFnXEqtMaV/HkVnaCU96FxKvWL/saiQ+GaQXR4yLEStch3/Q+dcjaq2u9HKH4KU3K+jQLb+2l86ZGWDIalQYTIcKztg9CSsE/SGZAyLiJaHVKolSV9GMr6RysGmEK/gC5+GhHbQ0iqlfrG8n0RHeF+G+KyeaTzAYbDt1I54+NZfUU6TKOer9Lb62Gwfxx64xpVN+2FBCFtWuKJ6c5OHkSMFOuUiLyCeOC3bG8by/R5aOVVtUtiXOSYDMuxwC/4zH3p8sk7gaQl8kwNs1pBXrX5J41PJR4q0TlDA9zGn0dZP2bmLKdW5szTAMURtDJjVU1FkrqKCUibDN8uEqVkJecON/muvgiF2b9DlFn/mY77SSE6MS33Z8ixyKthkYGhUKwyhmr3BlrhSRQx2pkZW70Jv3xCo7qXN8/U0z6MsK4aZpG62H6Gw2L1dnO6ImMTSfMTYdmzJYsGIK+EXUwjSgNE2Vg/ZNg4zPB6O431BpZKVu1jTTVC3tO7PKTP9ZRXGkGUH4dWN42NO1WhQzy8TW2krlFblIy9bkOLluWrRNT3n/x5DFrzShhtH2Z530EUWqV0kS6PIE+PcaTDkUZbBK1MeTFeSWOi8fF8MtHM8h4g8m+m0W2hoW7nVxr4TekUlcMh3+QxVB1kXFBGQkgZGBLlng7LO2GlHrpCm3hshGwzbNJN7kUsSRtlYfhm/yM47/0VsnRaFOr0CGknIxg7y5TzLtVfM3bFscDJLobD/6ZRzUVUeuFVNjm2+6ipHd6vsPHyYXJ+l2eG967hhSYT+iP5yaNU/CpmmrIuu46oWqkcIH2FrE4oyYYoulcSRAlZLzVdfoNChHjgJww9IRjJO3TmG/mRlG2sehooLKPM6A08zN/IcSx38t7klzT4sdVdyiGzG8wqlCPRRAM7Q756jOG0nrdy8xKivy5PIbbe8i5llMOOwqT8KkuMxmuIWj+CyXiF6jg91/VV0gM75uT9CZpQ5LDSgOwMVef5MX8egVVnQbEmQEsvhTfmI+wVMAP6tLrt+CQJE6Fcx+wr7Pj3Pg0zNon5f0CF1MKU/QsqSghp/0YfnRAF6K0R//eJOoUwyBCPWt2mCSbpGknJm9Q9obCv8FnyKpgsfmL9HP/uaVw6iXSHqLFIknsr0c/MxlHDOkM1znAiDZeg4f8Z8eATdCQ79NYPEaGW8/njuV9f6dGXcDbpegnyltU0NB5E3mSsWk0q0NvuJoJfx+sEADI9U0ujOsKQuAFGfQUKc37Nc3p0hZ+jcf07EawCus9/5dYN9d7/Qa37x0yvPSi03op8jY/8KoIgLDznh8vOio17wp6IXvE14VlJxnY1Y2HMhkXST8OU0KUzLiG3kJdFx4pastXdGRrBd+mVevKqT/OnzM7QQM2iZfYaM70DMfM7EGcI6XvEqGgt2njNaaLWEupWXoAQYRlklodxAcGEnzEsJuN7yGEkNA5+oXV4YTLBZCcW+B4zvRfo2FNZxs/wXjLt+2KNqq9ImQT5mKDoCuhgs2lIsnhbDvVzUlGOVPw4Pyvh0fNOQN96yBvvshZ/WIVEmWxgMa2FRS8rOJrRHnwButs+ummDN1qp1jCYkvUQ7OgiZ9iEmC4HIWYWNjO/oPqOLrJibAi9zIuSkCtbm4zZsJI0AhtJ/A423DEqQl6nH2pIIZMQCUgwY77/4KPbiQIfZcPJ3HMxTiJq6BH1M8qQIyFQdszoe6idXZkNGqKv8yfDir7vkthSthwa/Fr+NDEzO0jDkMVnPXzGbD65J90fKIIg0nUihxj9aUS936HzHCTqXUWjup91FIceS/Y2VpE6SNlIWch/daYr+SuRWXGxV5WutMKZtAN3fU3ByLpFY4cQju6EWd6F1BYwa2Z2qolD97HPFW4od30cJY5P0HaOIOz/KRHFgkCCzJ4o5bJ/WN3kYkWmPOvoBVqm7QmB9TH3wrPxxLs0WpJONi4dQBoyc4MNlpS868hQmmTKrpdwpbZMSXutlkiqj2xE0rAIcYNMXJT6DjxI4mk8eupIlzjLOpBv9jNsKQfvxcxZMrZUqjaNXtG9RAIXyy6Lr/X0/fEnQ5EM8iLVSXQ4RWPcSKT6FY2/VfEpvfV9vE4MT559uYT6kAFkciatYZbqC0uEdyqSLxxcq5e9c8QJ0iLvE6qOUyJcjMmFzSwdpTrIFna67/77oxsstLwks8Cw91u0TAKcpggxopUQdoOks3wg7fni/knBSB7VtiZsQFkGSTLNPleM4p+k2jLt+BwbYjd/L1HTkEX9fa8a+E+gOh74NZX0KvSEfL3l/TzbUyctjeplaIkWMTaooFLP9wb+UxufkyCk10hlKGBmmH523+fzvrpsNs46GpSFhnyIEUA8v4PlnaIMKRnbxbKQmIf/wQTiD/z9Jer/NK8nn7J9CHrzTcyke8o3fN0uzT9B4Dw6LuvAMJ9i2eJq8ZAEnUb2106XSYxez3Avr5rJ8Ji8/WSWbFXaNRw5lEolvfC6P8vKnUDKtJIk7HU1tSTP+QV+eXSIMLIwnvP+Wv/3kDIztbd9kuUaT3iVkHEOwa7P85ZGOHJ+ywqxcuJtGUQqGSYRDsl6VYY5sLpk7LK3k1amTmvcLAs9NWn/ancH7lBCdabc0Ho+h5R+HlLOb/BcTz9QBpGwH6tBOPBzxGOVysPVugtSb5lawyRGa5hC55hJz5/NRpzB0wy7Q9Tl8osYjywV1YmQ95usyxlYnJ+F2foe1iPdduKYMl29qevj/NmO8ty/M0QypD70tY9uCJB3xBheEoYy+GINRBcd8hyfYsWlkhMFxbwPuYu8daNJNAIyF3xcg9L0BnqRRudElFxQVmYxWDJvfSJKiUY2USnfI1oUUSkPdXOWnmtp7PJ6WeQlOtRVSKmXQoarr9TBAU2yjt87gJRxNf8ebkcuIqzwMvN6GhU5StJD454Ho+WtMNs+CqP1TvWihZ7JiJYpfnpO+UQ58kQJ9U1+qWdyEo/tU6MfWqKU7sILwenPdbos+Ena5a0eh/Wd0H3us9gQDv6FZN0KX9xDXmWiUd0LkxrXm+hK0rs1hFnG7JSEElW44RpyKEmyYswuZcEMZiXCBfXGpTzfFwH1avgj4JHwrofN+TV+R7LRvnViWUKP0tAbkLR+kGUb6f1CerBwBSKPJrKNf5MrDbMUeFrS4V7PCGAw38gwt56/y1pUsi9jD+dKN9ClFeFG6SDWe4gM/DuTyPuL2WradSyyk8dm1mEWnbV3FMVI3cbidQiQVsg2grr7P3l0Q0SbRDAVU7Mp852fh0nmP01oetsjVJ4mh6n9HqJ9PSDZJsPF+ITciFldPLabFd1Bw1rJyvegB+GbvCDg+Rp/DxKpvkSjktkYA0JMkmFNVsxhmpy68IrWSEJFagsZ1vfRjmX9LiYQNJKRRZ4tZRO9Xu5Qx3rFK/lYlldeSiGqiF5YifRPlURE+bcgZo+xDxRJngppXJOp7+2IMTs2CNKq7ggxLj1MxtkIhDciTDTXffyziQ1h1lderS5wfhnylvDoBpjHKUQXDZFGQ1KbEjKbcZev0Ul6/IueQnKcSpyBkaRXcbBkC43qy+Ri52FxPEAFyGo1g+ukib7CYxeSks3oxjLUJFkdqyKr4ZAjpb87HGq9saIJ/43c9n9Y19d5vEa03cxzLxGtn+XxNH9/Goi+qNoFkr0qCpRJGNolDMpbXbxXPLoHBjqWLCOlPpO2hRG+0HPQ3fsZMCtcpJAq3Zs9vkYevbCBdTms3Jb07yqUDOUlIwkhWldCIummcW2mBvXMWmaTqH9H8QGz7W4a29t4XQZHSUWIVr+lPoIMg0ztxzwG6lCdqpANpczSSz1cKHmjREL9X3n8TqFsynoXQACBUAJZjlycQjp6BXUY4hXfjG7l70QunWyOKd0pA+2BOlcD2CZFQ5KJEzSung06UwyD0xCS/SS/8KWrN+Ta77vwwWURjYu85kwatWRL33EO8fSIjsqKx47Qg7YhET9OqN5CPnMNDUveDhoCfWXZpvBjfP4SHtIfNkbD0GTTw6sZTU6lSb/qwb9M+hulaIK/Zh0fRUo3FSn7/YwQC/n7dJZ3fvexgHVfpZKQlGE1Q12J2hFEGVf8II2ngMqVeg2mCDq14EtYOXRKZocIlaAKZZKC0UBU++Y3vr2hd1mbyySKVxloWLIugxBg2S96/I0iSKszTEY09IpakExe+beqMcChkVATeYbPP4KU+XYWYTw7R8jmndnQMqwoz2YDDWnEl1XEQXjQqLThJ2lMs5C0P0DbkAF1KZ/Us+8hIZwHbUA2H1VvVMlbWLHDKmQiUZvWD+uavnePSEIyX70LEWfmnUq1q7AoIn1bui9/9Y6MM0gvrbAi0oseP0TPOJNulEFDBmMRyVryqRAn03o7LPbP8veeqbkZJOWHNvAwdWNE0vIu/hTSOg6RztSeOjDbG3M4vSQShcb/M+U4silpyvZplmu0uuU11ElK1t7SL4LseKuR9faj5GQyTUnCqBqY7hbqTU/US9D4YhFm58wUZV0IkcsMVX2FjSqzGmVrjh7kuihhfDffTLJOrqje6h6KTDP5jx3nxy1UIJFSzUIYp2jMKoxokp0MHzt4ou8Qz+UWokkqDI3/ByzL8yzKciRtD0I2Fh2dUfUIr5UMkSiVtN2vDo04behP0AR+xs9khmgPcpHDkWuabbI8kxmR4B94TjJNEnn1+RsiDCUGIhUzDMhwAVHk4kUU2HMMJTS4mBiBjs+XccG+8D5GoZIldKh1Usnv0m/zvBEiRuUnp/oxy8EMmbwxqRbQEwQdi1H1FXHMFHW0FAn7l/hzGfnaJhoYE55+90wwgZoMvXEZOW4lkUv0YHgjDYuFI4lPSoETNQwpB3lunCFp1CIN4KbyD9CumE1OxIxPTQ7RYSmN9Rh1XMMTb4RK40STHypOlKQhJC2ym8dEJWPCv6w01I/SuObTuJ6jcf2++7Me0cJoeTt/MmCGHuPP+BtpWCIstJpmTCWEN/LnpSa/TKmjsgO+W4WwCeFEmhRDOu8ljRt+lj8vt0r1fO5fWa/XWac1SFlJ1EfVYTsWERCwq9EJWZ9WE3qCz/w7z/d0bqeYJU4neV+j9jWKR3e90YZFqNVNoyfMJe85progLp5rDSPCHYTPkQ8kmV4rhV2syEsM8nY4MyEZP0SilScvIryOSaTvaT8b+jGV0SUtH+KjhxkUvyghCGgLybk+ziYqZNj93QDj0sMks0NgUqj1BhuWSBwpte4p+UnkFf68VI3C+ybb6NmVypDT6fNEiJSXziDDUzQqTfSF7nOXQSRpCPyYvxiIVB/jYy91VirGVUDkonFpc2hcv2V9t/O8gIF0nE6DwbwG0cguliv45KUw7zEKMwnf19kwXUDWw8xC5CVQIY8TKSTrzGwQfAQQr1Nrmk5U6KV/yvLXvm+w2A28/5doaG+heV3Mm0TDixp2E53JXkbkP/K8idfZUML6Js7y+d/m72w7x0NIb/KQVB3U3q4HiWbhjW+8YUnaHnwMycCvoHF8GlrLB3hyLOtgjUJSXbRbWffBB23Wz5keX+xKyQOFKBU/i6TvX5FKafkMkmnDHD5P0u+JFlKI0B+Q9P8EGuNaaNWaX5c7+GjILPayDN/lo23QOr/PKJweb/V7vvxmCIUUpu0aGRKQ8BR+STX+hAp5kDSCjO5rLbervpeJF/qnbM0maEj0Svq+SQcey/KWoxVtukGlg5e8TquGrS4hLx1SUjTqK8jpmYEmWF/v55V+heOZLO98kxiWNIo2l8C1jjZwjIrbxgJK18NEgCk9i9wtJWmydG2YJS2+VCEjQWVfDa2VJDp2msr+Fs8J8k4U5xKe2Ei/k41JtdDZPsUfvYv2X35JQCv1tX2Ev7bQmf6N5QvCYFzxZjEsiiYJrekWOp8LKf8PGVakT2i8c7V6hNWTV80DP+KvWWxwcpFL7t0xaOixGsttNOjtRMrvd5+fACGyJ33/jVSynnUhXTDIGmWXi1cNLRrzDepIRfdR1z/lCdObyLDkDV/dJIaq9/NXWr/EbsVPxuvt/J6sYeD/H+q+kw1Br9L37qB6aYUZkvUeotcCcqGnePyZ5y5W1UTe8F8UmmuM17IhZdWdN54ep8ugT6O0YSHr+hjL+fSbyLCUMJSYb05bf2QHCygNMt7xN3IRVjAV28f7XUcnkuUBLpd3U9kaC0PEA8xwC+nFv6BBMC1XU5PGIwbW4wDvw/ReVwat/XIg71hE6muDzv5x8lcXnfkHbzbDEtHQ+j9MvZWzgL8haB3mqbEO9ej4vVPdBLeEbXyPOnd5RRCY5NpOHiQD7r5vsUxVrMtY+ppoiExkUsGfIukhOU6RMMsUGM1E9cFNsOim0r7+hVX3vBkNi9bP7FDnoALJV1K+fyfQ9CzZOBqREOjn977DnxGixr3d/WJvRNigcRkWsQz3sQ7MnDxfIHLtZlFkySGZrZvuWBws5JY8nQr/Hcmue+hgv1LIp3V+nUnXSJtVvZGSpN+sJzDc8Sbpx8ooJir2d1Tq79kGb6NSmXGoKRsjCblI4CdMTn4HWN4NnXCrS9hROTrRIyncI/gX1iFA45jGBliX7kRVb0OJgXX3d8kb0jIB0f9jRYZVD7f1dmZfTGy0MgY4UZ26l0rEsaNvZsOipDQ0rG/QnvZA63qIXOl2nhuu41T6eHYSGWSDciYCzu8w5l/qGROjFZYtUUVn2cIqbKMdNfOUEVo9MzvzLTS061i3BB3iEV7zRLqeputpS3crtEqj1Ju3qfqLFv8PuBnMQY0+/lMAAAAASUVORK5CYII=' />\n",
    "\n",
    "# Topological phase data and analysis\n",
    "\n",
    "This Jupyter notebook performs all the analysis and generates the plots that appear in the paper _\"InAs-Al Hybrid Devices Passing the Topological Gap Protocol\"_ by Microsoft Quantum, available as a pre-print on [arXiv:2207.02472](https://arxiv.org/abs/2207.02472).\n",
    "See that paper for complete explanations and in-depth descriptions.\n",
    "\n",
    "The analysis here relies on the topological gap protocol as described in [arXiv:2103.12217](https://arxiv.org/abs/2103.12217).\n",
    "The protocol itself is a method to analyze a set of device conductance data and automatically determine whether a topological phase is present.\n",
    "It does this by testing for the presence of conductance peaks at zero applied voltage at both ends of the device and simultaneously testing for a gap closing and reopening in the bulk of the device.\n",
    "The topological phase is present when both these conditions are met, and the topological gap is the size of the gap opening in the bulk of the device.\n",
    "For further information on the structure of the devices themselves please refer to [the pre-print](https://arxiv.org/abs/2207.02472).\n",
    "\n",
    "The topological protocol is split in two stages.\n",
    "Stage 1 analyzes a broad set of _local_ conductance data to find zero bias peaks (ZBPs) and then identify regions of interest in the data containing zero bias peaks at each end of the device.\n",
    "These regions are then further analyzed in stage 2 by examining the _non-local_ conductance (current response at one end of the device to a voltage applied at the other end).\n",
    "When the non-local conductance shows a gap closing and re-opening (i.e., the peak in conductance vs.\n",
    "voltage bias applied converges at zero-bias and then moves away from zero bias), this is an indication that the device is in the topological phase.\n",
    "An annotated guide to stage 1 and stage 2 can be found in the notebooks \"stage-one-analysis\" and \"stage-two-analysis\".\n",
    "The exact same procedure is applied in this notebook.\n",
    "\n",
    "### Notes about input parameters\n",
    "\n",
    "This notebook contains a few one-off plots and two analysis/plot types that occur multiple times: stage 1 clusters and stage 2 phase diagrams.\n",
    "The analysis/plot parameters can be modified; main parameters are summarized in the following table.\n",
    "\n",
    "| Parameter name | Description | Example value |\n",
    "| --- | --- | --- |\n",
    "| `name` | Name of the device and stage, the saved plot's filename is based on this. | `\"deviceA_stage1\"` |\n",
    "| `selected_cutter` | The paired cutter gate index of which the gap value is used. | `0` |\n",
    "| `selected_field` | Magnetic field value at which a cut of the data is plotted. | `1.82` |\n",
    "| `selected_plunger` | Plunger gate value at which a cut of the data is plotted. | `-2.721` or `[-1.36575, -1.37275]` |\n",
    "| `T_mK` | Temperature in milli-Kelvin with which the simulation data is broadened. | `30.0` |\n",
    "| `thresholds` | Dictionary of stage 1 input parameters, see the `stage-one-analysis.ipynb` notebook for more details. | `{\"th_2w_p\": 0.5, \"th_3w\": 1e3, \"th_3w_tat\": 0.7}` |\n",
    "| `zbp_cluster_numbers` | The selected gapped ZBP cluster(s) that is (are) plotted. The clusters are sorted on size, where the biggest cluster is `1`.  | `[1]` |\n",
    "\n",
    "For the first step we load the analysis functions from [the installable Python package `tgp` available on PyPI and GitHub](https://github.com/microsoft/azure-quantum-tgp).\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "f6eb3d4a",
   "metadata": {},
   "outputs": [],
   "source": [
    "import shutil\n",
    "from pathlib import Path\n",
    "\n",
    "folder = Path(\"../data/\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "7e731ec9",
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "import xarray as xr\n",
    "from matplotlib.patches import Rectangle\n",
    "from paper_figures import (\n",
    "    analyze_two,\n",
    "    correct_two,\n",
    "    load_cached_broadened,\n",
    "    load_hmp_datasets,\n",
    "    load,\n",
    "    maybe_save_fig,\n",
    "    plot_1d_phase_diagram,\n",
    "    plot_clean_gap,\n",
    "    plot_conductance_cuts,\n",
    "    plot_conductance_traces,\n",
    "    plot_device_characterization,\n",
    "    plot_disorder_strength_vs_n2D,\n",
    "    plot_hall_bar_trans_long_resistance,\n",
    "    plot_Hall_bar,\n",
    "    plot_hybrid_mobility_protocol,\n",
    "    plot_local_and_non_local_conductance,\n",
    "    plot_pf,\n",
    "    plot_xi,\n",
    "    set_saving_figures,\n",
    ")\n",
    "\n",
    "import tgp\n",
    "import tgp.plot.paper\n",
    "\n",
    "# Change this flag to False to not save the figures to pdf\n",
    "# and set a folder where to save the figures\n",
    "set_saving_figures(True, folder=\"figures\")\n",
    "\n",
    "tgp.plot.set_mpl_rc()\n",
    "\n",
    "\n",
    "def text_kwargs(axs):\n",
    "    return dict(\n",
    "        c=\"k\",\n",
    "        ha=\"left\",\n",
    "        va=\"bottom\",\n",
    "        size=16,\n",
    "        transform=axs[0].transAxes,\n",
    "        zorder=2000,\n",
    "    )"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "accdd9cd",
   "metadata": {},
   "source": [
    "## Additional information about experimental devices\n",
    "The extraction of the effective parameters below are decribed in the main text.\n",
    "* `lever_arm`: see Sec. \"Electrostatic calibration from Hall bars\"\n",
    "* `g_factor` is defined as $|g_\\mathrm{eff}| \\equiv 2\\Delta_\\mathrm{ind} / \\mu_\\mathrm{B} B_\\mathrm{min}$, where $B_\\mathrm{min}$ is the lowest field at which the gap closes; see also Sec. \"TGP Stage 2: Measurement A1\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "3998646a",
   "metadata": {},
   "outputs": [],
   "source": [
    "device_info = dict(\n",
    "    deviceA1=dict(lever_arm=85, g_factor=5.6),\n",
    "    deviceA2=dict(lever_arm=85, g_factor=6.4),\n",
    "    deviceA3=dict(lever_arm=85, g_factor=6.4),\n",
    "    deviceB=dict(lever_arm=78, g_factor=3.7),\n",
    "    deviceC=dict(lever_arm=79, g_factor=4.4),\n",
    "    deviceD=dict(lever_arm=79, g_factor=4.5),\n",
    ")"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "806f885e",
   "metadata": {},
   "source": [
    "# Fig. 1: Energy spectrum & topological phase diagram\n",
    "\n",
    "The simplest, idealized model of Microsoft Quantum topological devices consists of a semiconductor nanowire coupled to a conventional superconductor, where the semiconductor nanowire has Rashba spin-orbit interaction.\n",
    "The effective model for this is shown in equation (1) of [the paper](https://arxiv.org/abs/2207.02472) and results in an energy spectrum (a) and topological phase diagram (b) as in figure 1 from the paper.\n",
    "The topological phase (in red) arises as Zeeman energy (from a magnetic field in-line with the nanowire) increases and creates a gap in the energy spectrum.\n",
    "The trivial superconducting phase is in blue."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0449a57c",
   "metadata": {},
   "outputs": [],
   "source": [
    "plot_1d_phase_diagram(\"energy_spectrum_topophase.pdf\")"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "6dc7625d",
   "metadata": {},
   "source": [
    "# Fig. 4: Phase diagram simulation of SLG-$\\beta$ and DLG-$\\varepsilon$ infinite clean device\n",
    "\n",
    "A realistic device is much more complicated than the above idealized device.\n",
    "Microsoft Quantum has modeled realistic device geometries including the actual semiconductor and superconductor materials as well as the 3-dimensional gate and junction shapes.\n",
    "The model here uses the single-layer junction device design on the $\\beta$ stack, but ignores any disorder or defects, and uses an infinite wire length.\n",
    "The simulation identifies multiple topological phases (as shown in red) at different \"plunger\" gate voltages (vertical axis) and magnetic fields (horizontal axis).\n",
    "As before, blue reflects the trivial superconducting phase.\n",
    "In real devices, the non-zero disorder and finite wire length will decrease the size of the topological phases.\n",
    "\n",
    "*(May take around 45 seconds to complete.)*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "e91194a6",
   "metadata": {},
   "outputs": [],
   "source": [
    "ds = xr.load_dataset(\n",
    "    folder / \"simulated\" / \"infinite_clean_gap_SLG.nc\", engine=\"h5netcdf\"\n",
    ")\n",
    "plot_clean_gap(\n",
    "    ds, plunger_min=-1.4, label=\"a\", fname=\"simulated_SLG_beta_infinite_clean_gap.pdf\"\n",
    ")"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "6cfb25a3",
   "metadata": {},
   "source": [
    "Like the previous figure, the following figure is a model of a realistic device, this time with a dual-layer junction device design on the $\\epsilon$ stack.\n",
    "Disorder and defects are again set to zero, and the wire is modeled with an infinite length.\n",
    "\n",
    "*(May take around 45 seconds to complete.)*"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "f31e9ea6",
   "metadata": {},
   "outputs": [],
   "source": [
    "ds = xr.load_dataset(folder / \"simulated\" / \"infinite_clean_gap_DLG.nc\")\n",
    "plot_clean_gap(\n",
    "    ds, plunger_min=-1.2, label=\"b\", fname=\"simulated_DLG_beta_infinite_clean_gap.pdf\"\n",
    ")"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "charitable-founder",
   "metadata": {},
   "source": [
    "# Fig. 6: Experimentally-measured non-local conductance across sections of lengths 1, 2, 3.5, 4, 6.5, 8 μm in the same wire"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "terminal-virtue",
   "metadata": {},
   "outputs": [],
   "source": [
    "hmp_datasets = load_hmp_datasets(folder / \"experimental\" / \"HMP\")\n",
    "fig, ds, ds_localization_length = plot_hybrid_mobility_protocol(\n",
    "    datasets=hmp_datasets,\n",
    "    ylims=(8e2 * 1e-3, 1e4 * 1e-3),\n",
    "    fname=\"HMP.pdf\",\n",
    ")"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "e39054a5",
   "metadata": {},
   "source": [
    "# Fig. 7: Phase diagram simulation of single-band disordered device\n",
    "\n",
    "The simulations now add a further degree of realism, including both an expected level of disorder as well as finite device length.\n",
    "The disorder is modeled as local variations in potential and is randomized to create different disorder realizations.\n",
    "The figures below are for simulations of two different disorder realizations corresponding to different disorder strengths.\n",
    "Here, red indicates points in the topological phase, and black indicates points in the topological case that additionally have a gap in excess of 10 µeV."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "de9914d9",
   "metadata": {},
   "outputs": [],
   "source": [
    "ds_pf = xr.load_dataset(folder / \"simulated\" / \"DLG_epsilon_stack_pfaffian.nc\")\n",
    "\n",
    "plot_pf(\n",
    "    ds_pf.isel(rms=0, seed=0),\n",
    "    ds_pf.isel(rms=1, seed=2),\n",
    "    ylim=(-1, 2),\n",
    "    fname=\"pd_pfaffian.pdf\",\n",
    ")"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "better-currency",
   "metadata": {},
   "source": [
    "# Fig. 9: Device A, local and non-local conductance\n",
    "\n",
    "Microsoft Quantum has measured a number of real devices, and has also simulated conductance measurements on models of real devices including the effects of disorder and finite wire length.\n",
    "The first real device, device \"A\" includes a full dataset of local and non-local conductance measurements.\n",
    "The following figure shows the data for the local conductance at the left end of the device (a), the right end (b), and the non-local conductance from one side to the other (c and d).\n",
    "These values are all measured simultaneously in the experiment.\n",
    "The horizontal axis is the plunger voltage, and the vertical axis is the bias applied to the junction at the end of the device (left or right side).\n",
    "The color scale indicates the magnitude of the conductance.\n",
    "From these conductance plots the induced gap (superconducting gap induced in the semiconducting nanowire) and the depletion point of the wire are determined."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "square-banana",
   "metadata": {},
   "outputs": [],
   "source": [
    "fname = folder / \"experimental\" / \"deviceA_conductance.nc\"\n",
    "ds_left, ds_right = (\n",
    "    xr.load_dataset(fname, engine=\"h5netcdf\", group=gr) for gr in [\"left\", \"right\"]\n",
    ")\n",
    "\n",
    "plot_local_and_non_local_conductance(ds_left, ds_right, fname=\"deviceA_conductance.pdf\")"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "53e740e6",
   "metadata": {},
   "source": [
    "# Fig. 10: Device A, measurement A1 (stage 1)\n",
    "\n",
    "The protocol has also been applied to experimental measurements of real devices.\n",
    "Device A is analyzed for stage 1 in the same way as the simulated data.\n",
    "The regions in the magnetic field and plunger gate voltage space with zero bias peaks on the left (a), and right of the device (b), followed by the joint fraction — having a simultaneous zero bias peak on both sides — in (c).\n",
    "\n",
    "We use a manually set threshold `th_3w = 0.3e7` to filter out noise in $3\\omega$ signal. \n",
    "Note that the positions and sizes of the regions with zero bias peaks are insensitive to the choice of `th_3w` threshold within an order of magnitude."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "e7107354",
   "metadata": {},
   "outputs": [],
   "source": [
    "name = \"deviceA_stage1\"\n",
    "thresholds = dict(\n",
    "    set_gapped={\"th_2w_p\": 0.7},\n",
    "    set_3w_th={\"th_3w\": 0.3e7},\n",
    "    set_3w_tat={\"th_3w_tat\": 0.7},\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "e1ce9951",
   "metadata": {},
   "outputs": [],
   "source": [
    "ds = load(folder / \"experimental\" / f\"{name}.nc\")\n",
    "tgp.one.analyze(ds, thresholds=thresholds)\n",
    "\n",
    "tgp.plot.paper.plot_stage1_clusters(\n",
    "    ds.L_3w_ta,\n",
    "    ds.R_3w_ta,\n",
    "    ds.L_3w_tat,\n",
    "    ds.R_3w_tat,\n",
    "    thp=ds.attrs[\"set_3w_tat.th_3w_tat\"],\n",
    "    xy_text=(0.04, 0.70),\n",
    "    text=\"Device A\",\n",
    "    ylim=(-1.18, -1.15),\n",
    ")\n",
    "maybe_save_fig(f\"{name}.pdf\")"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "8345dbbc",
   "metadata": {},
   "source": [
    "# Fig. 11: Device A, measurement A1 (stage 2)\n",
    "\n",
    "Device A is then analyzed in the region of interest using non-local conductance for stage 2.\n",
    "As in figure 13, a region of magnetic field and plunger voltage is identified that has a gapless boundary, a gap in the interior, and zero bias peaks (a), and therefore constitutes an SOI$_2$.\n",
    "The size of the gap is shown by the red scale in (b), with blue being the gap size in the trivial superconducting regime.\n",
    "Selected cuts of the conductance through the identified topological region are shown for local conductance (c, d) as well as non-local (e, f).\n",
    "As with the simulated data, the real device data here shows zero bias peaks in the local conductance and a closing and re-opening of a gap in the non-local conductance as magnetic field is increased.\n",
    "(g, h, i, and j) show a waterfall view of the conductance data with bias now on the x-axis and magnetic field on the y-axis.\n",
    "The presence of zero bias peaks are now visible as a peak replicated vertically (g, h) and finite conductance values converging and then re-opening (i, j).\n",
    "The dots in (i, j) show extracted gap values."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "8ca84e36",
   "metadata": {},
   "outputs": [],
   "source": [
    "name = \"deviceA1_stage2\"\n",
    "zbp_cluster_numbers = [1]\n",
    "selected_cutter = 0\n",
    "selected_plunger = -1.17175\n",
    "selected_field = 1.66"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "fa6ff27c",
   "metadata": {},
   "outputs": [],
   "source": [
    "fnames = (\n",
    "    folder / \"experimental\" / f\"{name}_left.nc\",\n",
    "    folder / \"experimental\" / f\"{name}_right.nc\",\n",
    ")\n",
    "ds_left, ds_right = (load(fn) for fn in fnames)\n",
    "\n",
    "(ds_left, ds_right), (ds_c_left, ds_c_right) = correct_two(\n",
    "    ds_left,\n",
    "    ds_right,\n",
    "    lockin_left=\"lockin_1\",\n",
    "    lockin_right=\"lockin_2\",\n",
    "    fridge_parameters=folder / \"fridge\" / \"deviceA1.yaml\",\n",
    "    drop_indices=[0, 1, 2],\n",
    "    max_bias_index=-2,\n",
    "    norm=1e3,\n",
    "    phase_shift_left=-3.3,\n",
    "    phase_shift_right=-5.57,\n",
    ")\n",
    "\n",
    "result = analyze_two(ds_c_left, ds_c_right)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "6c36f459",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, axs = tgp.plot.paper.plot_stage2_diagram(\n",
    "    result.zbp_ds,\n",
    "    cutter_value=selected_cutter,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    plunger_cut=selected_plunger,\n",
    "    field_cut=selected_field,\n",
    "    device_info=device_info[name.split(\"_\")[0]],\n",
    "    pct_boundary_shifts={1: [0.28, -0.001]},\n",
    ")\n",
    "plt.text(-0.476, 1 - 0.07, \"Device A,\", **text_kwargs(axs))\n",
    "plt.text(-0.476, 1 - 0.13, \"measurement A1\", **text_kwargs(axs))\n",
    "maybe_save_fig(f\"{name}_diagram.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "9b67687a",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, (ax1, ax2, ax3, ax4) = tgp.plot.paper.plot_conductance(\n",
    "    result.zbp_ds,\n",
    "    result.ds_left,\n",
    "    result.ds_right,\n",
    "    selected_cutter=selected_cutter,\n",
    "    selected_plunger=selected_plunger,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    field_ticks=[0.5, 1, 1.5, 2.0, 2.5],\n",
    "    g_local_max=0.8,\n",
    "    g_nonlocal_max=0.08,\n",
    ")\n",
    "maybe_save_fig(f\"{name}_conductance.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "a615df59",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, (ax1, ax2, ax3, ax4) = tgp.plot.paper.plot_conductance_waterfall(\n",
    "    result.zbp_ds,\n",
    "    result.ds_left,\n",
    "    result.ds_right,\n",
    "    selected_cutter=selected_cutter,\n",
    "    selected_plunger=selected_plunger,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    field_lim=[0.72, 2.73],\n",
    "    bias_max=0.06,\n",
    "    scale_g_local=0.15,\n",
    "    gbar_local=0.5,\n",
    "    scale_g_nonlocal=0.6,\n",
    "    gbar_nonlocal=0.1,\n",
    ")\n",
    "maybe_save_fig(f\"{name}_conductance_waterfall.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "0320c7f0",
   "metadata": {},
   "source": [
    "# Fig. 12: Device A, measurement A1 (stage 2), plunger waterfall plot\n",
    "\n",
    "The zero-bias peaks in local conductance (a,b) and gap closing and re-opening in non-local conductance (c, d) are also visible over a range of plunger voltage values (with magnetic field held constant), demonstrating that the identified topological phase is robust over both magnetic field and plunger voltage."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "17758ae9",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, axs = tgp.plot.paper.plot_conductance_waterfall_plunger(\n",
    "    result.zbp_ds,\n",
    "    result.ds_left,\n",
    "    result.ds_right,\n",
    "    selected_cutter,\n",
    "    selected_field,\n",
    "    zbp_cluster_numbers,\n",
    "    plunger_lim=(-1.1753, -1.1663),\n",
    "    bias_max=0.06,\n",
    "    bias_ticks=np.arange(-50, 51, 25),\n",
    "    scale_g_local=0.0007,\n",
    "    gbar_local=0.5,\n",
    "    scale_g_nonlocal=(0.006, 0.0015),\n",
    "    gbar_nonlocal=(0.05, 0.25),\n",
    "    g_bar_pos=-0.0003,\n",
    "    text=\"Device A, measurement A1\",\n",
    ")\n",
    "maybe_save_fig(f\"{name}_conductance_waterfall_plunger.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "53fbcd02",
   "metadata": {},
   "source": [
    "# Fig. 13: Device A, measurement A2 (stage 2)\n",
    "\n",
    "The following generates the same plots as for figures 13 and 15, but now for a different measurement of device A after the device had been warmed up and re-measured in a different experimental setup."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "324bd5d5",
   "metadata": {},
   "outputs": [],
   "source": [
    "name = \"deviceA2_stage2\"\n",
    "zbp_cluster_numbers = [1]\n",
    "selected_cutter = 1\n",
    "selected_plunger = -1.4045\n",
    "selected_field = 1.14666667"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "65b911d5",
   "metadata": {},
   "outputs": [],
   "source": [
    "fnames = (\n",
    "    folder / \"experimental\" / f\"{name}_left.nc\",\n",
    "    folder / \"experimental\" / f\"{name}_right.nc\",\n",
    ")\n",
    "ds_left, ds_right = (load(fn) for fn in fnames)\n",
    "\n",
    "(ds_left, ds_right), (ds_c_left, ds_c_right) = correct_two(\n",
    "    ds_left,\n",
    "    ds_right,\n",
    "    lockin_left=\"mfli_5510\",\n",
    "    lockin_right=\"mfli_5602\",\n",
    "    fridge_parameters=folder / \"fridge\" / \"deviceA2.yaml\",\n",
    "    drop_indices=[0, 1, 2],\n",
    "    max_bias_index=-2,\n",
    "    norm=1e3,\n",
    ")\n",
    "\n",
    "result = analyze_two(ds_c_left, ds_c_right)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "b5fd4f96",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, axs = tgp.plot.paper.plot_stage2_diagram(\n",
    "    result.zbp_ds,\n",
    "    cutter_value=selected_cutter,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    plunger_cut=selected_plunger,\n",
    "    device_info=device_info[name.split(\"_\")[0]],\n",
    "    pct_boundary_shifts={1: [0.24, -0.0016]},\n",
    ")\n",
    "plt.text(-0.467, 1 - 0.07, \"Device A,\", **text_kwargs(axs))\n",
    "plt.text(-0.467, 1 - 0.13, \"measurement A2\", **text_kwargs(axs))\n",
    "maybe_save_fig(f\"{name}_diagram.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "d3935c6d",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, (ax1, ax2, ax3, ax4) = tgp.plot.paper.plot_conductance(\n",
    "    result.zbp_ds,\n",
    "    result.ds_left,\n",
    "    result.ds_right,\n",
    "    selected_cutter=selected_cutter,\n",
    "    selected_plunger=selected_plunger,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    bias_max=145,\n",
    "    g_local_max=0.8,\n",
    "    g_nonlocal_max=0.025,\n",
    ")\n",
    "maybe_save_fig(f\"{name}_conductance.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "beda5210",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, (ax1, ax2, ax3, ax4) = tgp.plot.paper.plot_conductance_waterfall(\n",
    "    result.zbp_ds,\n",
    "    result.ds_left,\n",
    "    result.ds_right,\n",
    "    selected_cutter=selected_cutter,\n",
    "    selected_plunger=selected_plunger,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    field_lim=[0.37, 2.25],\n",
    "    bias_max=0.1,\n",
    "    scale_g_local=0.2,\n",
    "    gbar_local=0.5,\n",
    "    scale_g_nonlocal=2.0,\n",
    "    gbar_nonlocal=0.05,\n",
    ")\n",
    "maybe_save_fig(f\"{name}_conductance_waterfall.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "42ca9cb0",
   "metadata": {},
   "source": [
    "# Fig. 14: Device A, measurement A3 (stage 2)\n",
    "\n",
    "The same device was again re-measured after one week, yielding the following data and analysis (same figure order as figures 13 and 15)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0b4b20cc",
   "metadata": {},
   "outputs": [],
   "source": [
    "name = \"deviceA3_stage2\"\n",
    "zbp_cluster_numbers = [1]\n",
    "selected_cutter = 1\n",
    "selected_plunger = -1.4083\n",
    "selected_field = 1.1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "469f2faf",
   "metadata": {},
   "outputs": [],
   "source": [
    "fnames = (\n",
    "    folder / \"experimental\" / f\"{name}_left.nc\",\n",
    "    folder / \"experimental\" / f\"{name}_right.nc\",\n",
    ")\n",
    "ds_left, ds_right = (load(fn) for fn in fnames)\n",
    "\n",
    "(ds_left, ds_right), (ds_c_left, ds_c_right) = correct_two(\n",
    "    ds_left,\n",
    "    ds_right,\n",
    "    lockin_left=\"mfli_5510\",\n",
    "    lockin_right=\"mfli_5602\",\n",
    "    fridge_parameters=folder / \"fridge\" / \"deviceA2.yaml\",\n",
    "    drop_indices=[],\n",
    "    max_bias_index=-1,\n",
    "    norm=1e3,\n",
    ")\n",
    "\n",
    "result = analyze_two(ds_c_left, ds_c_right)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "8d93fa40",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, axs = tgp.plot.paper.plot_stage2_diagram(\n",
    "    result.zbp_ds,\n",
    "    cutter_value=selected_cutter,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    plunger_cut=selected_plunger,\n",
    "    device_info=device_info[name.split(\"_\")[0]],\n",
    "    pct_boundary_shifts={1: [0.34, -0.0002]},\n",
    ")\n",
    "plt.text(-0.48, 1 - 0.07, \"Device A,\", **text_kwargs(axs))\n",
    "plt.text(-0.48, 1 - 0.13, \"measurement A3\", **text_kwargs(axs))\n",
    "maybe_save_fig(f\"{name}_diagram.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "88b51ab4",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, (ax1, ax2, ax3, ax4) = tgp.plot.paper.plot_conductance(\n",
    "    result.zbp_ds,\n",
    "    result.ds_left,\n",
    "    result.ds_right,\n",
    "    selected_cutter=selected_cutter,\n",
    "    selected_plunger=selected_plunger,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    bias_max=100,\n",
    "    bias_ticks=np.arange(-90, 91, 30),\n",
    "    g_local_max=1,\n",
    "    g_nonlocal_max=0.08,\n",
    ")\n",
    "maybe_save_fig(f\"{name}_conductance.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "af5f924e",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, (ax1, ax2, ax3, ax4) = tgp.plot.paper.plot_conductance_waterfall(\n",
    "    result.zbp_ds,\n",
    "    result.ds_left,\n",
    "    result.ds_right,\n",
    "    selected_cutter=selected_cutter,\n",
    "    selected_plunger=selected_plunger,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    field_lim=[0.57, 1.8],\n",
    "    bias_max=0.1,\n",
    "    scale_g_local=0.1,\n",
    "    gbar_local=1.0,\n",
    "    scale_g_nonlocal=0.6,\n",
    "    gbar_nonlocal=0.1,\n",
    ")\n",
    "maybe_save_fig(f\"{name}_conductance_waterfall.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "569caae6",
   "metadata": {},
   "source": [
    "# Fig. 15: Device A, measurement A3 (stage 2) bias-field cuts\n",
    "\n",
    "Plots (a) through (d) show the stability of the zero bias peak for 4 different cutter gate settings.\n",
    "Plot (e) shows a composite of the 4 conductance measurements through the zero bias peak."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "ca263c73",
   "metadata": {},
   "outputs": [],
   "source": [
    "plot_conductance_cuts(\n",
    "    ds_c_left, selected_plunger, selected_field, f\"{name}_left_conductance.pdf\"\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "5cc96dc5",
   "metadata": {},
   "outputs": [],
   "source": [
    "plot_conductance_traces(\n",
    "    ds_c_left, selected_field, selected_plunger, f\"{name}_left_conductance_cut.pdf\"\n",
    ")"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "5869b64b",
   "metadata": {},
   "source": [
    "# Fig. 16: Device B (stage 2)\n",
    "\n",
    "The stage 2 analysis applied to device B."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "bf3071f3",
   "metadata": {},
   "outputs": [],
   "source": [
    "name = \"deviceB_stage2\"\n",
    "zbp_cluster_numbers = [1]\n",
    "selected_cutter = 1\n",
    "selected_plunger = -1.15775\n",
    "selected_field = 2.12"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "80d520fc",
   "metadata": {},
   "outputs": [],
   "source": [
    "fnames = (\n",
    "    folder / \"experimental\" / f\"{name}_left.nc\",\n",
    "    folder / \"experimental\" / f\"{name}_right.nc\",\n",
    ")\n",
    "ds_left, ds_right = (load(fn) for fn in fnames)\n",
    "\n",
    "(ds_left, ds_right), (ds_c_left, ds_c_right) = correct_two(\n",
    "    ds_left,\n",
    "    ds_right,\n",
    "    lockin_left=\"lockin_1\",\n",
    "    lockin_right=\"lockin_2\",\n",
    "    fridge_parameters=folder / \"fridge\" / \"deviceB.yaml\",\n",
    "    drop_indices=[0, 1, 2, 3],\n",
    "    max_bias_index=-4,\n",
    "    norm=1e3,\n",
    "    phase_shift_left=-3.78,\n",
    "    phase_shift_right=-6.89,\n",
    ")\n",
    "\n",
    "result = analyze_two(ds_c_left, ds_c_right, gap_threshold_high=0.1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "66571685",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, axs = tgp.plot.paper.plot_stage2_diagram(\n",
    "    result.zbp_ds,\n",
    "    cutter_value=selected_cutter,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    plunger_cut=selected_plunger,\n",
    "    field_cut=selected_field,\n",
    "    device_info=device_info[name.split(\"_\")[0]],\n",
    "    pct_boundary_shifts={1: [0.33, 0.00033]},\n",
    ")\n",
    "\n",
    "plt.text(-0.495, 1 - 0.07, \"Device B\", text_kwargs(axs))\n",
    "\n",
    "maybe_save_fig(f\"{name}_diagram.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "6f6ff8ac",
   "metadata": {},
   "outputs": [],
   "source": [
    "p_test_arr = (\n",
    "    [selected_plunger]\n",
    "    if isinstance(selected_plunger, (float, np.floating))\n",
    "    else selected_plunger\n",
    ")\n",
    "p = selected_plunger\n",
    "fig, (ax1, ax2, ax3, ax4) = tgp.plot.paper.plot_conductance(\n",
    "    result.zbp_ds,\n",
    "    ds_c_left,\n",
    "    ds_c_right,\n",
    "    selected_cutter=selected_cutter,\n",
    "    selected_plunger=p,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    bias_max=90,\n",
    "    bias_ticks=np.arange(-80, 81, 40),\n",
    "    g_local_ticks=None,\n",
    "    g_local_max=0.5,\n",
    "    g_nonlocal_ticks=None,\n",
    "    g_nonlocal_max=0.03,\n",
    ")\n",
    "\n",
    "maybe_save_fig(f\"{name}_conductance.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "3c6bada3",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, (ax1, ax2, ax3, ax4) = tgp.plot.paper.plot_conductance_waterfall(\n",
    "    result.zbp_ds,\n",
    "    ds_c_left,\n",
    "    ds_c_right,\n",
    "    selected_cutter=selected_cutter,\n",
    "    selected_plunger=selected_plunger,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    field_lim=[1.35, 3],\n",
    "    bias_max=0.09,\n",
    "    bias_ticks=np.arange(-80, 81, 40),\n",
    "    scale_g_local=0.19,\n",
    "    gbar_local=0.5,\n",
    "    scale_g_nonlocal=0.7,\n",
    "    gbar_nonlocal=0.1,\n",
    ")\n",
    "\n",
    "maybe_save_fig(f\"{name}_waterfall.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "cbd3090a",
   "metadata": {},
   "source": [
    "# Fig. 17: Device B (stage 2), plunger waterfall plot\n",
    "\n",
    "Device B waterfall plunger, compare with Fig. 20."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "3bb23902",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, axs = tgp.plot.paper.plot_conductance_waterfall_plunger(\n",
    "    result.zbp_ds,\n",
    "    result.ds_left,\n",
    "    result.ds_right,\n",
    "    selected_cutter,\n",
    "    selected_field,\n",
    "    zbp_cluster_numbers,\n",
    "    bias_max=0.09,\n",
    "    bias_ticks=np.arange(-80, 81, 40),\n",
    "    plunger_lim=(-1.1582, -1.1543),\n",
    "    scale_g_local=0.0007,\n",
    "    gbar_local=0.2,\n",
    "    scale_g_nonlocal=0.01,\n",
    "    gbar_nonlocal=0.02,\n",
    "    g_bar_pos=-0.00015,\n",
    "    text=\"Device B\",\n",
    ")\n",
    "maybe_save_fig(f\"{name}_waterfall_plunger.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "61024f6e",
   "metadata": {},
   "source": [
    "# Fig. 18: Device C (stage 2)\n",
    "\n",
    "The stage 2 analysis applied to the device C."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "4f5905cf",
   "metadata": {},
   "outputs": [],
   "source": [
    "name = \"deviceC_stage2\"\n",
    "zbp_cluster_numbers = [1]\n",
    "selected_cutter = 0\n",
    "selected_plunger = -2.3655\n",
    "selected_field = 0.98"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "988af58b",
   "metadata": {},
   "outputs": [],
   "source": [
    "fnames = (\n",
    "    folder / \"experimental\" / f\"{name}_left.nc\",\n",
    "    folder / \"experimental\" / f\"{name}_right.nc\",\n",
    ")\n",
    "ds_left, ds_right = (load(fn) for fn in fnames)\n",
    "\n",
    "(ds_left, ds_right), (ds_c_left, ds_c_right) = correct_two(\n",
    "    ds_left,\n",
    "    ds_right,\n",
    "    lockin_left=\"mfli_5583\",\n",
    "    lockin_right=\"mfli_5591\",\n",
    "    fridge_parameters=folder / \"fridge\" / \"deviceC.yaml\",\n",
    "    drop_indices=([0, 1, 2, 3, 4, 5], []),\n",
    "    max_bias_index=(-4, -1),\n",
    "    norm=1e3,\n",
    ")\n",
    "\n",
    "result = analyze_two(ds_c_left, ds_c_right)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "61637204",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, axs = tgp.plot.paper.plot_stage2_diagram(\n",
    "    result.zbp_ds,\n",
    "    cutter_value=selected_cutter,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    plunger_cut=selected_plunger,\n",
    "    device_info=device_info[name.split(\"_\")[0]],\n",
    "    pct_boundary_shifts={1: [0.20, 0.0001]},\n",
    ")\n",
    "plt.text(-0.474, 1 - 0.07, \"Device C\", **text_kwargs(axs))\n",
    "maybe_save_fig(f\"{name}_diagram.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "bab56ed0",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, (ax1, ax2, ax3, ax4) = tgp.plot.paper.plot_conductance(\n",
    "    result.zbp_ds,\n",
    "    result.ds_left,\n",
    "    result.ds_right,\n",
    "    selected_cutter=selected_cutter,\n",
    "    selected_plunger=selected_plunger,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    field_ticks=[0.5, 1, 1.5, 2.0, 2.5],\n",
    "    g_local_max=0.4,\n",
    "    g_nonlocal_max=0.08,\n",
    ")\n",
    "maybe_save_fig(f\"{name}_conductance.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "5c556caa",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, (ax1, ax2, ax3, ax4) = tgp.plot.paper.plot_conductance_waterfall(\n",
    "    result.zbp_ds,\n",
    "    result.ds_left,\n",
    "    result.ds_right,\n",
    "    selected_cutter=selected_cutter,\n",
    "    selected_plunger=selected_plunger,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    field_lim=[0.685, 1.66],\n",
    "    bias_max=0.06,\n",
    "    scale_g_local=0.1,\n",
    "    gbar_local=0.5,\n",
    "    scale_g_nonlocal=0.5,\n",
    "    gbar_nonlocal=0.1,\n",
    "    g_bar_pos=-0.04,\n",
    ")\n",
    "maybe_save_fig(f\"{name}_conductance_waterfall.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "e111268e",
   "metadata": {},
   "source": [
    "# Fig. 19: Device D (stage 2)\n",
    "\n",
    "The stage 2 analysis applied to the device D."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "ec53cac2",
   "metadata": {},
   "outputs": [],
   "source": [
    "name = \"deviceD_stage2\"\n",
    "zbp_cluster_numbers = [1, 2]\n",
    "selected_cutter = 0\n",
    "selected_plunger = -2.721\n",
    "selected_field = 1.82"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "4ed44308",
   "metadata": {},
   "outputs": [],
   "source": [
    "fnames = (\n",
    "    folder / \"experimental\" / f\"{name}_left.nc\",\n",
    "    folder / \"experimental\" / f\"{name}_right.nc\",\n",
    ")\n",
    "ds_left, ds_right = (load(fn) for fn in fnames)\n",
    "\n",
    "(ds_left, ds_right), (ds_c_left, ds_c_right) = correct_two(\n",
    "    ds_left,\n",
    "    ds_right,\n",
    "    lockin_left=\"mfli_4909\",\n",
    "    lockin_right=\"mfli_5654\",\n",
    "    fridge_parameters=folder / \"fridge\" / \"deviceD.yaml\",\n",
    "    drop_indices=[0, 1],\n",
    "    max_bias_index=-4,\n",
    "    norm=1e3,\n",
    ")\n",
    "\n",
    "result = analyze_two(ds_c_left, ds_c_right)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "4f864893",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, axs = tgp.plot.paper.plot_stage2_diagram(\n",
    "    result.zbp_ds,\n",
    "    cutter_value=selected_cutter,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    plunger_cut=selected_plunger,\n",
    "    device_info=device_info[name.split(\"_\")[0]],\n",
    "    pct_boundary_shifts={1: [-0.3, 0.0012], 2: [0.32, -0.0012]},\n",
    ")\n",
    "plt.text(-0.495, 1 - 0.07, \"Device D\", **text_kwargs(axs))\n",
    "maybe_save_fig(f\"{name}_diagram.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "19c3f2af",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, (ax1, ax2, ax3, ax4) = tgp.plot.paper.plot_conductance(\n",
    "    result.zbp_ds,\n",
    "    result.ds_left,\n",
    "    result.ds_right,\n",
    "    selected_cutter=selected_cutter,\n",
    "    selected_plunger=selected_plunger,\n",
    "    zbp_cluster_numbers=[1],\n",
    "    g_local_max=1,\n",
    "    g_nonlocal_max=0.045,\n",
    ")\n",
    "maybe_save_fig(f\"{name}_conductance.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "9589626a",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, (ax1, ax2, ax3, ax4) = tgp.plot.paper.plot_conductance_waterfall(\n",
    "    result.zbp_ds,\n",
    "    result.ds_left,\n",
    "    result.ds_right,\n",
    "    selected_cutter=selected_cutter,\n",
    "    selected_plunger=selected_plunger,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    field_lim=[0.69, 2.18],\n",
    "    bias_max=0.1,\n",
    "    scale_g_local=0.1,\n",
    "    gbar_local=0.5,\n",
    "    scale_g_nonlocal=0.8,\n",
    "    gbar_nonlocal=0.1,\n",
    "    g_bar_pos=-0.06,\n",
    ")\n",
    "maybe_save_fig(f\"{name}_conductance_waterfall.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "349ad770",
   "metadata": {},
   "source": [
    "# Fig. 20: Device E (stage 2)\n",
    "\n",
    "The stage 2 analysis applied to the device E.\n",
    "This device does _not_ pass the protocol. Although it has regions with zero bias peaks on both ends of the device (a), the gap value is zero (b).\n",
    "Plots (c)-(f) show the local and non-local conductance along the dashed line indicated in plot (b), demonstrating the presence of local zero bias peaks, but no gap re-opening."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "31729bcb",
   "metadata": {},
   "outputs": [],
   "source": [
    "name = \"deviceE_stage2\"\n",
    "zbp_cluster_numbers = []\n",
    "selected_cutter = 1\n",
    "selected_plunger = -1.3893"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0c3c72a0",
   "metadata": {},
   "outputs": [],
   "source": [
    "fnames = (\n",
    "    folder / \"experimental\" / f\"{name}_left.nc\",\n",
    "    folder / \"experimental\" / f\"{name}_right.nc\",\n",
    ")\n",
    "ds_left, ds_right = (load(fn) for fn in fnames)\n",
    "\n",
    "(ds_left, ds_right), (ds_c_left, ds_c_right) = correct_two(\n",
    "    ds_left,\n",
    "    ds_right,\n",
    "    lockin_left=\"mfli_5591\",\n",
    "    lockin_right=\"mfli_5583\",\n",
    "    fridge_parameters=folder / \"fridge\" / \"deviceE.yaml\",\n",
    "    drop_indices=[0],\n",
    "    max_bias_index=-1,\n",
    "    norm=1e3,\n",
    ")\n",
    "\n",
    "result = analyze_two(ds_c_left, ds_c_right)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "b5ae73d1",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, axs = tgp.plot.paper.plot_stage2_diagram(\n",
    "    result.zbp_ds,\n",
    "    cutter_value=selected_cutter,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    plunger_cut=selected_plunger,\n",
    "    pct_boundary_shifts={1: [0.34, -0.001], 2: [0.34, -0.001]},\n",
    ")\n",
    "plt.text(-0.495, 1 - 0.07, \"Device E\", **text_kwargs(axs))\n",
    "maybe_save_fig(f\"{name}_diagram.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "7c02adcc",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, (ax1, ax2, ax3, ax4) = tgp.plot.paper.plot_conductance(\n",
    "    result.zbp_ds,\n",
    "    result.ds_left,\n",
    "    result.ds_right,\n",
    "    selected_cutter=selected_cutter,\n",
    "    selected_plunger=selected_plunger,\n",
    "    zbp_cluster_numbers=[],\n",
    "    g_local_max=1,\n",
    "    g_nonlocal_max=0.01,\n",
    ")\n",
    "maybe_save_fig(f\"{name}_conductance.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "cbc90e38",
   "metadata": {},
   "source": [
    "# Fig. 21: Device F (stage 1)\n",
    "\n",
    "The device F does not pass stage 1 of the protocol.\n",
    "Presence of zero bias peaks on the left (a) and right (b) of the device, and the joint probability (c).\n",
    "Zero bias peaks exist on both sides, but not at the same time.\n",
    "\n",
    "Here we use threshold `th_3w = 1e7` (slightly higher compared to Measurement A1). The device F does not pass stage 1 of the topological gap protocol for `th_3w` one order of magnitude larger or one order of magnitude smaller than the used value."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "222aaf4f",
   "metadata": {},
   "outputs": [],
   "source": [
    "name = \"deviceF_stage1\"\n",
    "thresholds = dict(\n",
    "    set_gapped={\"th_2w_p\": 0.5},\n",
    "    set_3w_th={\"th_3w\": 1e7},\n",
    "    set_3w_tat={\"th_3w_tat\": 0.7},\n",
    "    set_clusters={\"min_samples\": 1},\n",
    ")\n",
    "ds = load(folder / \"experimental\" / f\"{name}.nc\")\n",
    "tgp.one.analyze(ds, thresholds=thresholds)\n",
    "tgp.plot.paper.plot_stage1_clusters(\n",
    "    ds.L_3w_ta,\n",
    "    ds.R_3w_ta,\n",
    "    ds.L_3w_tat,\n",
    "    ds.R_3w_tat,\n",
    "    thp=ds.attrs[\"set_3w_tat.th_3w_tat\"],\n",
    "    xy_text=(0.04, 0.70),\n",
    "    text=\"Device F\",\n",
    ")\n",
    "maybe_save_fig(f\"{name}.pdf\")"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "38d7ab86",
   "metadata": {},
   "source": [
    "# Fig. 23: Effective $\\delta V$ describing disorder strength\n",
    "\n",
    "Self-consistent Poisson-Schrödinger simulations are used to determine a mapping between different charge defect density (x-axis) and resulting disorder parameter (y-axis).\n",
    "The charge density occurs at the semiconductor/dielectric interface."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "48a93259",
   "metadata": {},
   "outputs": [],
   "source": [
    "df_β = pd.read_csv(folder / \"simulated\" / \"SLG_beta_stack.csv\")\n",
    "df_δ = pd.read_csv(folder / \"simulated\" / \"DLG_delta_stack.csv\")\n",
    "df_ε = pd.read_csv(folder / \"simulated\" / \"DLG_epsilon_stack.csv\")\n",
    "plot_disorder_strength_vs_n2D(df_β, df_δ, df_ε, \"disorder_strength_vs_n2D.pdf\")"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "81972eae",
   "metadata": {},
   "source": [
    "#  Fig. 24: Majorana length for DLG $\\varepsilon$-stack\n",
    "\n",
    "The zero-energy disordered coherence length $\\xi(0)$ evaluated numerically for the dual-layer junction device design on $\\epsilon$ stack.\n",
    "The colorscale highlights the divergence of $\\xi(0)$, indicating the phase transition from the low-field trivial phase to the high-field topological phase.\n",
    "The solid and dashed lines mark contours of constant $\\xi(0) = 1$ and 3 µm respectively.\n",
    "For our 3 µm devices, the region enclosed by dashed lines roughly represents the \"critical\" regime, while the area between dashed and solid lines is in the \"crossover\" regime where $L > \\xi(0)$, but we still expect significant finite-size effects and mesoscopic fluctuations.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "3d21ee74",
   "metadata": {},
   "outputs": [],
   "source": [
    "xi = xr.load_dataset(folder / \"simulated\" / \"DLG_epsilon_stack_xi.nc\")\n",
    "plot_xi(\n",
    "    xi.isel(rms=0),\n",
    "    xi.isel(rms=1),\n",
    "    xlims=(0.5, 3.6),\n",
    "    ylims=(-1, 2),\n",
    "    fname=\"xi0_epsilon_stack.pdf\",\n",
    ")"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "noticed-courage",
   "metadata": {},
   "source": [
    "# Fig. 25: Device A, transverse and longitudinal resistance in Hall bar"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "dynamic-opinion",
   "metadata": {},
   "outputs": [],
   "source": [
    "ds_hb = xr.load_dataset(\n",
    "    folder / \"experimental\" / \"Hall_bar_B_field.nc\", engine=\"h5netcdf\"\n",
    ")\n",
    "device_image_fname = folder / \"experimental\" / \"Hall_bar_device.png\"\n",
    "plot_hall_bar_trans_long_resistance(ds_hb, device_image_fname, \"Hall_bar_measurement.pdf\")"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "76c3fdd5",
   "metadata": {},
   "source": [
    "# Fig. 26: Density and mobility from Hall bar close to device A\n",
    "\n",
    "Measured density and mobility vs. gate voltage from a Hall bar nearby to device A.\n",
    "The orange points are used for the fit line. The solid lines are simulation results."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "49b55b19",
   "metadata": {},
   "outputs": [],
   "source": [
    "plot_Hall_bar(folder / \"experimental\" / \"Hall_bar.nc\")"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "020fc249",
   "metadata": {},
   "source": [
    "# Fig. 28: Simulated device DLG-$\\varepsilon$-R1 (stage 1)\n",
    "Microsoft Quantum has applied the topological gap protocol to _simulated_ conductance for realistic device models.\n",
    "Since simulations allow computation of topological invariants, they allow us to validate the protocol against for known topological phases.\n",
    "The first step of the protocol detects zero-bias peaks in the local conductance of the left end of the device (a), then the right end (b), and then generates the overlapping regions (c).\n",
    "The data is shown with magnetic field on the x-axis and plunger voltage on the y-axis.\n",
    "\n",
    "The plots below combine data for different cutter gate voltages, with the color scale indicating the fraction of cutter values (out of 5 in total) for which the ZBP detection algorithm returned a positive result. It therefore indicates the stability of the ZBPs with respect to changes in cutter gate voltage."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "7f11bdd0",
   "metadata": {},
   "outputs": [],
   "source": [
    "T_mK = 40.0\n",
    "\n",
    "name = \"simulated_DLG_epsilon_R1_stage1\"\n",
    "fname = folder / \"simulated\" / f\"{name}.nc\"\n",
    "\n",
    "ds = load_cached_broadened(fname, T_mK)\n",
    "\n",
    "ds = tgp.prepare.prepare_sim(ds, T_mK=0)  # data is already broadened\n",
    "\n",
    "thresholds = dict(\n",
    "    set_gapped={\"th_2w_p\": 0.5},\n",
    "    set_3w_th={\"th_3w\": 1e3},\n",
    "    set_3w_tat={\"th_3w_tat\": 0.7},\n",
    ")\n",
    "tgp.one.analyze(ds, thresholds=thresholds)\n",
    "\n",
    "fig, (ax1, ax2, ax3) = tgp.plot.paper.plot_stage1_clusters(\n",
    "    ds.L_3w_ta,\n",
    "    ds.R_3w_ta,\n",
    "    ds.L_3w_tat,\n",
    "    ds.R_3w_tat,\n",
    "    thp=0.7,\n",
    "    text=\"Simulated,\" + \"\\n\" + r\"DLG-$\\varepsilon$-R1\",\n",
    ")\n",
    "\n",
    "ax3.add_patch(\n",
    "    Rectangle(\n",
    "        xy=(1.75, -0.73),\n",
    "        width=3 - 1.75,\n",
    "        height=-0.707 + 0.73,\n",
    "        fill=False,\n",
    "        ec=\"k\",\n",
    "        ls=\"--\",\n",
    "    )\n",
    ")\n",
    "\n",
    "maybe_save_fig(f\"{name}.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "ad1f4d63",
   "metadata": {},
   "source": [
    "# Fig. 29: Simulated device DLG-$\\varepsilon$-R1 (stage 2)\n",
    "\n",
    "Once promising regions are identified from stage 1 as in the previous figure, the 2nd stage of the protocol can be run.\n",
    "For this simulated device conductance, the 2nd stage identifies a region of magnetic field and plunger voltage that has a gap in the non-local conductance, has a closing of the gap at the boundaries, and has zero-bias peaks in the local conductance throughout the region, and determines that region to be a subregion of interest (SOI$_2$) (a).\n",
    "The simulation independently verifies the same region is a known topological phase (b).\n",
    "The simulated local (c, d) and non-local (e, f) conductance is shown for a cut through the topological region, showing the zero bias peaks in local conductance and the closing and re-opening of the gap in non-local conductance.\n",
    "(g,h,i and j) show another view (waterfall plot) of the simulated conductance data.\n",
    "The dots in (i, j) show extracted gap values."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "01cb43c8",
   "metadata": {},
   "outputs": [],
   "source": [
    "T_mK = 40.0\n",
    "\n",
    "name = \"simulated_DLG_epsilon_R1_stage2\"\n",
    "\n",
    "fname = folder / \"simulated\" / f\"{name}.nc\"\n",
    "\n",
    "selected_cutter = 4\n",
    "zbp_cluster_numbers = [1]\n",
    "selected_plunger = -0.7205\n",
    "\n",
    "ds = load_cached_broadened(fname, T_mK)\n",
    "\n",
    "ds_left, ds_right = ds.rename({\"bias\": \"left_bias\"}), ds.rename({\"bias\": \"right_bias\"})\n",
    "\n",
    "result = analyze_two(ds_left, ds_right)\n",
    "\n",
    "fig, axs = tgp.plot.paper.plot_stage2_diagram(\n",
    "    ds=result.zbp_ds,\n",
    "    cutter_value=selected_cutter,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    plunger_cut=selected_plunger,\n",
    "    pct_boundary_shifts={1: [-0.55, 0.009]},\n",
    "    invariant=\"SI\",\n",
    ")\n",
    "axs[0].plot([2.18, 2.28], [-0.7115, -0.7203], lw=1, color=\"k\", linestyle=\"--\")\n",
    "plt.text(-0.462, 1 - 0.0, \"Simulated,\", **dict(text_kwargs(axs), va=\"top\"))\n",
    "plt.text(-0.462, 1 - 0.06, r\"DLG-$\\varepsilon$-R1\", **dict(text_kwargs(axs), va=\"top\"))\n",
    "maybe_save_fig(f\"{name}_diagram.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "606e55e1",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, (ax1, ax2, ax3, ax4) = tgp.plot.paper.plot_conductance(\n",
    "    result.zbp_ds,\n",
    "    result.ds_left,\n",
    "    result.ds_right,\n",
    "    selected_cutter=selected_cutter,\n",
    "    selected_plunger=selected_plunger,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    g_local_max=0.8,\n",
    "    g_nonlocal_max=0.02,\n",
    "    labels=\"cdef\",\n",
    ")\n",
    "maybe_save_fig(f\"{name}_conductance.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "fe237c4a",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, (ax1, ax2, ax3, ax4) = tgp.plot.paper.plot_conductance_waterfall(\n",
    "    result.zbp_ds,\n",
    "    result.ds_left,\n",
    "    result.ds_right,\n",
    "    selected_cutter=selected_cutter,\n",
    "    selected_plunger=selected_plunger,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    field_lim=[1.72, 3.3],\n",
    "    bias_max=0.07,\n",
    "    scale_g_local=0.15,\n",
    "    gbar_local=0.5,\n",
    "    scale_g_nonlocal=3.0,\n",
    "    gbar_nonlocal=0.025,\n",
    ")\n",
    "maybe_save_fig(f\"{name}_conductance_waterfall.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "9f30618c",
   "metadata": {},
   "source": [
    "# Fig. 30: Simulated device SLG-$\\beta$-R1 (stage 1)\n",
    "\n",
    "This figure is similar to Fig. 12, however, the simulations below have been performed for a single layer gate device."
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "c65d238a",
   "metadata": {},
   "source": [
    "Panel (a) and (b) show the fraction of cutter values where the ZBP detection algorithm returned a positive result, measuring at the left and right side, respectively. Panel (c) is the joint probability of the former panels."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "957b0a14",
   "metadata": {},
   "outputs": [],
   "source": [
    "name = \"simulated_SLG_beta_R1_stage2\"\n",
    "zbp_cluster_numbers = [1]\n",
    "selected_plunger = [-1.37275, -1.36575]\n",
    "selected_cutter = 2\n",
    "T_mK = 30.0"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "9c226624",
   "metadata": {},
   "outputs": [],
   "source": [
    "fname = folder / \"simulated\" / f\"{name}.nc\"\n",
    "\n",
    "ds = load_cached_broadened(fname, T_mK)\n",
    "ds_left, ds_right = ds.rename({\"bias\": \"left_bias\"}), ds.rename({\"bias\": \"right_bias\"})\n",
    "\n",
    "result = analyze_two(ds_left, ds_right)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0307e66f",
   "metadata": {},
   "outputs": [],
   "source": [
    "def plot_simulation_cutter_averaged(zbp_ds, thp=0.7):\n",
    "    L_3w_ta = zbp_ds.left\n",
    "    R_3w_ta = zbp_ds.right\n",
    "\n",
    "    L_3w_tat = (L_3w_ta > thp) * 1.0\n",
    "    R_3w_tat = (R_3w_ta > thp) * 1.0\n",
    "    return tgp.plot.paper.plot_stage1_clusters(\n",
    "        L_3w_ta,\n",
    "        R_3w_ta,\n",
    "        L_3w_tat,\n",
    "        R_3w_tat,\n",
    "        thp,\n",
    "        text=\"Simulated,\\nSLG-$\\\\beta$-R1\",\n",
    "    )\n",
    "\n",
    "\n",
    "plot_simulation_cutter_averaged(result.zbp_ds, thp=0.7)\n",
    "maybe_save_fig(f\"{name.replace('stage2', 'stage1')}.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "ccc1af0a",
   "metadata": {},
   "source": [
    "# Fig. 31: Simulated device SLG-$\\beta$-R1 (stage 2)\n",
    "\n",
    "Once promising regions are identified from stage 1 as in the previous figure, the 2nd stage of the protocol can be run.\n",
    "For this simulated device conductance, the 2nd stage identifies a region of magnetic field and plunger voltage that has a gap in the non-local conductance, has a closing of the gap at the boundaries, and has zero-bias peaks in the local conductance throughout the region, and determines that region to be topological (a).\n",
    "The simulation independently verifies the same region is a known topological phase (b).\n",
    "The simulated local (c, d) and non-local (e, f) conductance is shown from the simulated data through a section of the topological region, showing the zero bias peaks in local conductance and the closing and re-opening of the gap in non-local conductance.\n",
    "\n",
    "Plots (g-j) show the local and non-local conductance outside the topological region where there is no closing and re-opening of the transport gap.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "4f1ca806",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, axs = tgp.plot.paper.plot_stage2_diagram(\n",
    "    ds=result.zbp_ds,\n",
    "    cutter_value=selected_cutter,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    plunger_cut=selected_plunger,\n",
    "    pct_boundary_shifts={1: [-0.25, 0.0041]},\n",
    "    plunger_lim=[-1.385, -1.36],\n",
    "    invariant=\"SI\",\n",
    ")\n",
    "plt.text(-0.462, 1 - 0.0, \"Simulated,\", **dict(text_kwargs(axs), va=\"top\"))\n",
    "plt.text(-0.462, 1 - 0.06, \"SLG-$\\\\beta$-R1\", **dict(text_kwargs(axs), va=\"top\"))\n",
    "\n",
    "for i, labels in enumerate([\"(c-f)\", \"(g-j)\"]):\n",
    "    axs[1].text(\n",
    "        2.22,\n",
    "        selected_plunger[i] - 0.0004,\n",
    "        labels,\n",
    "        size=text_kwargs(axs)[\"size\"],\n",
    "        va=\"top\",\n",
    "    )\n",
    "\n",
    "maybe_save_fig(f\"{name}_diagram.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "7ecfbc89",
   "metadata": {},
   "outputs": [],
   "source": [
    "for pc, labels in zip(selected_plunger, [\"cdef\", \"ghij\"]):\n",
    "    fig, axs = tgp.plot.paper.plot_conductance(\n",
    "        result.zbp_ds,\n",
    "        result.ds_left,\n",
    "        result.ds_right,\n",
    "        selected_cutter=selected_cutter,\n",
    "        selected_plunger=pc,\n",
    "        zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "        g_local_max=0.8,\n",
    "        g_nonlocal_max=0.05,\n",
    "        labels=labels,\n",
    "    )\n",
    "    maybe_save_fig(f\"{name}_conductance_{labels}.pdf\")\n",
    "    plt.show()"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "04519bcc",
   "metadata": {},
   "source": [
    "# Fig. 32: Simulated device SLG-beta-R2 (stage 2)\n",
    "\n",
    "Stage 2 analysis for a simulated device that fails the protocol.\n",
    "The local and non-local conductances in (c)-(f) are taken along the dashed line in (b)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "4eac3112",
   "metadata": {},
   "outputs": [],
   "source": [
    "name = \"simulated_SLG_beta_R2_stage2\"\n",
    "zbp_cluster_numbers = []\n",
    "selected_cutter = 2\n",
    "selected_plunger = -1.36"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "5ff2ccd5",
   "metadata": {},
   "outputs": [],
   "source": [
    "ds = xr.load_dataset(folder / \"simulated\" / f\"{name}.nc\")\n",
    "ds = tgp.prepare.prepare(ds)\n",
    "\n",
    "T_mK = 30.0\n",
    "ds = tgp.prepare.broaden_with_temperature(ds, T_mK)\n",
    "\n",
    "ds_left, ds_right = ds.rename({\"bias\": \"left_bias\"}), ds.rename({\"bias\": \"right_bias\"})\n",
    "\n",
    "result = analyze_two(ds_left, ds_right)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "22116907",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, axs = tgp.plot.paper.plot_stage2_diagram(\n",
    "    result.zbp_ds,\n",
    "    cutter_value=selected_cutter,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    plunger_cut=selected_plunger,\n",
    "    pct_boundary_shifts={1: [-0.18, 0.0041]},\n",
    ")\n",
    "plt.text(-0.46, 1 - 0.0, \"Simulated,\", **dict(text_kwargs(axs), va=\"top\"))\n",
    "plt.text(-0.46, 1 - 0.06, \"SLG-$\\\\beta$-R2\", **dict(text_kwargs(axs), va=\"top\"))\n",
    "maybe_save_fig(f\"{name}_diagram.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "52a9d259",
   "metadata": {},
   "outputs": [],
   "source": [
    "fig, (ax1, ax2, ax3, ax4) = tgp.plot.paper.plot_conductance(\n",
    "    result.zbp_ds,\n",
    "    result.ds_left,\n",
    "    result.ds_right,\n",
    "    selected_cutter=selected_cutter,\n",
    "    selected_plunger=selected_plunger,\n",
    "    zbp_cluster_numbers=zbp_cluster_numbers,\n",
    "    bias_max=100,\n",
    "    g_local_max=1.5,\n",
    "    g_nonlocal_max=0.08,\n",
    ")\n",
    "maybe_save_fig(f\"{name}_conductance.pdf\")\n",
    "plt.show()"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "downtown-exchange",
   "metadata": {},
   "source": [
    "# Fig. 33 and 34: Device A, measurement A1 (stage 2) for different cutter values\n",
    "Figures 33 and 34 can be reproduced using the code for Fig. 11 above and setting `selected_cutter = 1` and `selected_cutter = 2`, accordingly."
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "loving-opinion",
   "metadata": {},
   "source": [
    "# Fig. 36: Device A, characterization"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "neural-utilization",
   "metadata": {},
   "outputs": [],
   "source": [
    "name = \"deviceA_characterization\"\n",
    "ds = xr.load_dataset(folder / \"experimental\" / f\"{name}.nc\")\n",
    "plot_device_characterization(ds, fname=f\"{name}.pdf\")"
   ]
  }
 ],
 "metadata": {
  "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.11.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}