From 18c6550f2ae852cb3d91f2f57fd21ff3c4e05a7f Mon Sep 17 00:00:00 2001 From: Lisa Date: Tue, 20 Aug 2019 16:45:13 +0200 Subject: [PATCH] stimuli figures --- contrast_sweep.dat.pdf | Bin 0 -> 22012 bytes ephys_analysis.py | 48 ++++++++++++++++++++++------------------- phase_sweep.dat.pdf | Bin 0 -> 23876 bytes stimuli_figures.py | 33 ++++++++++++++++++++++++++++ zap.dat.pdf | Bin 0 -> 13585 bytes 5 files changed, 59 insertions(+), 22 deletions(-) create mode 100644 contrast_sweep.dat.pdf create mode 100644 phase_sweep.dat.pdf create mode 100644 stimuli_figures.py create mode 100644 zap.dat.pdf diff --git a/contrast_sweep.dat.pdf b/contrast_sweep.dat.pdf new file mode 100644 index 0000000000000000000000000000000000000000..aec55b51e506c9e55a42ffb835c2255e32edf33a GIT binary patch literal 22012 zcma&OWmKF?ur3@35InfMySoN=cMb0DPH=a3*WhkJ2baMmxO-sGK@vXpk?)>;?p}B8 zUo+jSyQ}J{>h7wU_jxIm#U&Y8nKXRvWugsnzkQk6K zt63SFSvy#gaQ)|Ml8z2;B<%mO)c;Y3mF2zTzYP#^aBy^UeP925va+#@`TO|)RG(SZ z+||+D#q`~tzx03CMasp|-HGJ?SpA=EW-&DqW^r>5Yg2Posdt^;Tk0;x4z5mrEoSQd z7Yxh_=4RH$qK;ngV*ZlN$;HLQ%f`#f&d$Th!%o7^!NSDK$<4vT!p_S1mmdEK>R(vC zL;G&QyHc(sy#JxZ|FPn)_J128X>Ioo8VR$c-8+)v=BAEj?=$7i9W33fNLc@J*ME4! z)y>7+*dE?1x7J^--a*e@|Lrw;Qy$N5*$=qJ;Szk1PngO<4z&VhW@qa3`l!ZPJ))1_ zvAUufTIZ(NIP-jD+Vl0h-|Lay6L@L#7KH!$=tz9HUje?-+9rCvY~LQqEvU$!=lk>g z`Rm)y=UbhTGyLmYslXTE?cU$N+lk*W4R3s3O9}#B+VX*dmZz8r{z5*8^-MrZ`{n?N zw_W<@Dd7C14d!)*mS5O@MZjay8}r;l?=Nd}zdH}1b)WoJt-v!WgZt)1XOH@lfY(V+ z;uB?T6`W_?6EGeO48>f7YfNMJDnKcm`6DWAJds|H`rX(A zDSbRGy{y-;yZ5|W{?E~#j=A^9UV+QeNastnJBOmcT~4Ttd{y5IRBBaV<$8A(DCGDl z*S?nh?r5t)a}YWr&im;wqZUh!SFKFfL)#1+g_GBD6%Mc(Gv zViL-X_R3zk1g2XGu2U}^smc6|pyq7Q|KWRQ$uv-y!FqppyNEh!Ks1{QrZP0N<=+Go z9efrKp>ei3f zRI{10r?g(tc-;!)??2VGYiaEH4UHFkf+aN0JIGzf3=Oq|z889Z;U!6dD5))#us=#4 zDXPBU`H356Tj2b3flnb;r!m~M>K^Bq${1v1~$!TD%=(J|OkFGkEt0MnPRSg<(BaX9*00%KgR!dwB9Y@jDV=DejSs z7pf;^j6u7Iv{k;rW#JxV`N?dFONfEql12rc-YHSQT?V%JXBSby&a_+#O1eFD0J8K6 z4)bWY0f|1`3W_qy++fjUK2|V7y`%C-7vs+6d>YTMHhe?0x`8LPPy;d1`7gVZj|)YI zvBU^kOcp|i05be~{Xvc5?kS!L7|EmjZ0uVGIH(ey;RI*PXdAqOyKt3{(LmN;$?CoQM!Y$hEX= zAu`WBT!b8wprRlbyT##YE_%@&%Q5ktGx}B1fON@6$Q-cuNvgG1k0gneZ~m$V^{#%k zuu{g@X=%<(8T5uPf?6_|Ze{Vq^le0>m~m|S&);hemu;YX3T|bgywHS0UZ;jYre|9tbm5V@C8rXjkyWS-PoQtY9(BGP1!L#bD#-%CT#*Ei!XEj)QsptV`IR z33B=m;0{?-kE363rUOHJw>qXzzofk~@{D|hmAHFrm?E)lq5BgTjn@H&=QrmdZ2c4R z2W?{kT;cjOzoP~#JPHBn-QuMo@JUp=YaWDj9y#6zWX77)-I zl8;bVS-G!uSUED4x|AzxH1X{lup5$hhPkE}eXEaxa&oB5k=&fD#S8(dtwcp1SV7AXJ5F5ZekRbi!rSZeVuzS_6MPaoo&X_bFQTQo;x8&0c#(s=E zr~G50=5z-oGxqil9Doe>*V)5H5V3Tlc~9RTm5TN16^^r(@{&+18b18~O<6Rr5cl#N zNFhe%(?}>*M(t4G{Sp6o5q`nKF9B=cwFSst!Xt=Yr=F%d?9e@S>y43HcTQBn>i~Oa zWgr%24V;865Y9H}Q3t8r1g5kNf|_AJq-njTfuRWE%KJxPW|0JHqQra(Dc;>J#v7eA za7!3V*ft}gkt}#(AVOpMacw_HieKBbc4JVOlM9EfEo!dcT3%VB!bgS ze=$H&sP8-utKWwA=gNqJofJYtK%^!d*3LQWCzIdc7o8EMv?#nOmUx0Th2&lf)J2Pw zQ+NW4`eNA94JyJ0`lVxlet+d8xFz@iRNr#$855{KoKu>c_TB-7 z%Q6oKN`WDrka)Bo;6hi#b6X6uaPdZS&`dNrVto6+ zV^?aVm>;kyR}P{=4QJamR#kBuDp@nLcr_v*78))t2#YD#Vh9hKKOaIEk?||@j7ZaJ z@Y$m_wV5?{$7DuTbskqdjL!#oMdj=k{Q^TSW>7o}r8Y-@-Ps{lK4ICaCGK?NQE38j zDQP^jo|G#v-m(~l^yi>r6^}Dpu5Bs zRpY$v!)Lfz5OPi5fV#oRIt#JM@%JpHZy_2LX&1DH$lG|IQf(&3ANx)Q+?e6EVXCjndomk?*xHTpd8{Hu47VFR9q{v8+>j6uDPF@^p zC;?rP6quAt#s7kUwqU7O*2f_JOf$=e#zV246-9hLh9Xr9?lO|;6{3hQ9q-{Q75${1 zNgibW{gxCn`?HMKp;w+$Ljr?_P#^6|Z2d}9sK}MkanVY?hhKrwn~j&QjLQ>>Et?>& zV7ssCm`{W2hoO>hxITH(`p>cJl*n^M((Fo4vC?RZ5nP@Q@JzvWf`g2LC!EJ8IjV(XtJvZLy!A`&->;zEaVWWJ(Be8vvMH+;bLPM^Y5 z-ye-mo~&q_wI#u?r41=_ei#CK`H#;?7c1JcXmi<(ey*k+npN2Y3=Y)@$&ZaN>2Yd9 z7Sn=KdJoXCIs(!#PWtzeH#xIGrV=u*y!I!Y2+R2fM_0FdNPSu!B zyyr<|#6}SZD)bmiMQCN~h$X8WesEBg9LTDs4HV~&kZure=d~=O3E|$FWu8h;2L1UG za_VlMnPr@Y_UjT4=p-h3f^K*!Gv5`SRX6;(jz>U}iD=K3bAO==%?p}x#lcTcE1cs@ zk3*XZJuxfe;u_+CtCtXNSVu7u(;lAAXDf4ND_$3+URHkCup*c^)Uqo%7FqNGo*H_B z?GSYc>T(_2QhjiC?Cp1la(60ga?hg~6*kzgH0?JCmljRl%kH7;bx3W-aDgj{Q_UC^ z^06jL>$|K&yJ_AqSlTifH;XkPpBiq$`rz4lJefa@!b))igAt%tbXz18Cfw@nw0R~dxq!e+wVeL+8fM(Vedq2k| zcYqH`fvMJ1qFF7(?>@^2PmvxS%A4yVH!^w=zxW(JDg@j1Jlwt5t^P#LGFH#|DJK?N z+;|zkIQ$<$tNc=v1hF%do@R(W+AJ(p^A2|~o|UlBAolk!Gv@a`J34wAXc4$A@u27D zR6@NwU#b-3tgU`Vp{9h)&@Aqn5aHM)uk|Bw&10TcxgN{KGz$TB6QMrCuju>&y>)~d z>@(W?t9ITjNVk>G60Ak68K4wOtQZgPJYC_Bq4Pzg5VCR|74J`^rQL*^=)U;hQW4h2W;$tLuMij{b zZaK`UW|yy!<)W)az^K_Gh~OcB&Qw4C!zI@tuNVVM3$^0NBuR29;kWkK8O`t%V%#rN z`#k})i`%+S`GSsra^V$TE_O?14qiJ4ugH4!&l_<*>TX#cmy0chGX){fb|U)h|D-{0 z<@~BNPn8gPOAKrJ!;;OP33>EPgpG{S!Zm*3(u*i}0sZTJ`MUGG2F1ABOz$qK_G^XK zK6{LqW++9zUx3c1(O#mk!^naa1A8OpD}=!Jpc8G@4n_qtaoa%!)ee^a{A~LP@SN5r z%1W5nLO2WETa88xSU(@l7v16zdd!)L2+shR6^~|39^oqmZ`UfbfgBE}?%A{EV-pxmc6&2FaP}I$rTd7D&Q`^-?W(52Q^41alPnnD ztL*X1>7AejN0oz8NO$@Z5ScAv48<%}6!3@+m$6yksGNnR_mS_}R6g$a;5%x9`d5BY ztKI=B+y#K)RSAATZ(DNh@2E702*+Y#()|7MetHp0ERd8jpD<%$kilTSUMCcAkHn`6 z?lkRh655T4AnNcC%I+X7YzP;SQ=sG?0;UMkEqhcbd{J8Wefqou8VVF{_Lze^>oJ?xFRVh%RuxS&R=I$7B*<|; zjlxB^8%5$Q>s=`os=aQ0T}=4s^)D&q#Kq2*)H$gg&`U6|P`As84oR?;|f)}I{G=_^gqXYpoTp%S(@)%kvYNa1bq|W z$8+0^=zHY-TRDZ}EhV!c8~wx(Lly562~8M|;WtD<0TXV+glWjZkMb#wSr7fR5{2KQ zclWJ0OK=nQ$(+ioa`PyQ=1zZ(v`<6sT$`uLrh^C44LRx`T;-vEl^ZDNs&{KJw?}*Vq`%7ZoQ{vkL|r)tIn} zLQt~ZC&OmMYPJTe(=mN*pPPO8Mq$Jcky3{B!th1oHv`p3NWnFMzM9YLNh2ENPAK7xg&l z9!nj;xY4Q|0nm@F1N7mfEwP`THAsuzBz$PZZP=NpmJdI03s$gj;4pxrS!zGoa=|5e z66>}x^Q!_k)Kv@t8jOA8#?dMy#jo5Pev58@N>5z-d{$dt4O%R@8 z8`!zx&fv@a1h7WMLpD0m{sVU`XwS1uLF?tK-OXHut`+mOb_5p_pR#NY&7b%t>SN$4 z6M6iS*Q8Ay=15=Zjj_WB2u}8VvsIBRz>FZ+8Wh{jW|;vB$Fs%*;a{7(nn}z(8%AB# z+*HSs9Qoq82V|1R*GIl(9nf=oxADdJ#u8IE(k6u4V$%D{*~L43RTfS5?6=hn1QAr% zAor6Gp{sK#NO%C0`gl3DkrMzIu7NMn*y*BEfQp!m&P|A%7QC;Nv;O_%G3IV^s|ljad;{ASTVZ2y{~d+U zHWqr~OzMQWq?2}$Ks+yjg|uSz+ce8Lhfs4cBDU`rxTn{)OSF3yYUiIL+^kZghKi%c z>7OLuAjwIY-Id`)lggU|usrASUsDy($9qyW0E?r{l89d29vm@YpElx15r{1m& zw@XuXg>A$>ya<|Mjee=LH!1@-Hw`Q-&)67yNQNhSO^WiUKB(xbre&m!dyNoA$aH zpnmk`)??kOcYIU^Pzp4c#@~y}!{Kb>n(JI_;u+Kv+r3ROQI4phZQ*k4_I76n7j~?k zQ9>*XdPHsrfSzI0`h{%{$c5OfR8Kq>%O4ZT^_OKO)@}tm=DbVFQI?xfqY@a}^ycD@ zC*yMDOtd#uT&(E1<}Hqz1_EvmlPT*+iat?l8iy~1BRxXTuXVXKg~nWGtAE+Mk8A!C zT3r_dCTRCRD0;%-VXr@yMblIbMN;TV!&fl>T#F%Ro|71J`z@pzAex35!~F-6NE4mr zeM4ZuU~2whlJg}{GA?^<*!k4k%*SlzlLTud`w-XM|0M@aFin2f3x4vn-DUzWXH?&f zWLW$4&j^uj6yoYOil}OA+TvBUpf09Zmi6B4h)yAOwdTSa`o;3{vH`4UM?orf%%#@m zb<@)zp4Mx@UXAR$C+$~+?#mO-kT^$N|DY?*b=Iv%h0H+|Vq1LSd-zzP0?Tv7N=gFr zxDPh8n_gfJK!5+bcAyxWe@8+j*oVKnAJ6=pa2JzqWUvOwE*)UG4CZXg$3S=e;dm0( zK$N0u5b4OhsCpGvA_AV}N>~1-0ITJ%vCk>@rd|56Y%TVY8_Y&q91i-; zwIfmn%|u_yWaRLY3Ni$0$QTD&SDP`V6PFshB0@_qEm z@#09WjAUW7Wb1eZ-`AKZ5aG%$SUE8a1W$0w(%W>I0Y77o3`nR&ts_o|9L?gGfgc>t zJ%9_$@%oMCXKd4pQ{b56^zm3`xlydXGdbhzl{U$pWXlBu;$b*>#H7d4oK!8k@Wv15 z7sBQ$JoFemy_gdh7FQ;Z0&?cWZT>ZHHig@*AFS}HqAfq+ih+j$nYC))4DDvbjKn}a zP`0OU1lg5Mo&K>(iPLv`gvBo)wK- zT(6b9aG${ni&FTMOnFwuWOP*waT!4*@q7I?+8hBF{Un6VYtoG}qbe`xGKXud^Zonh zIIz!p%dSYI_o^aTUk3uAJaUcCQkO@WKT8Z}o}ga}OIWcSN!>-l?*~Ns&jg=`b+5tl zX)eAEU0T14n+yd8m6cRTHeDZ+CyC8<)%@fqQuIhK z*||hV9*0FZJ!E7&L!{ZodomoRzq8dPp(bidgr|qL{mV&L6T2Y6HicFvI%5lGJfk!+ zUPt0jR6$K@>AE1t?;=+6OD7rs*Jz^qYNBJ8@_9LwnTDGBUL!FHPW|rv28y4$6uU6T zY#koSN3Uy~U)i^C;Kkf~yHQQ82@t-)%8XDN-s-icPaPBCtho7 zU-hZq;Q`&Y)j#GHg@{NxO|uUp>9xcG1MY{-%tcsxo1jwK4IU?sf%qa2`|*(;+!pnd z({Bz(uU(+oGOh+|fsyB|iWy}%HQAS8xn+wH-m!DD=jXUW3|MgbyH8$i^0Kq$Ec`8T z8Wq8n^?6=E@z3pY5tCrQI!}R&rk|mnFv`V2S`y(6q&&W74{y!_V(h&J7-&?Agnx7t zDH_3fvfV8kYsHp+KS1amNP5W7DVvw&F^|eP>DPHP93{I`T4ToI@hJ@mu5Kx?I>XRK zX1&k)J%;Q!6nBv!WDGkzGWxYPAUxMzWc!-K^`lS>Da+0|2ZIWW5Mx5QuTRz^mf?y{ zDG}o~`^T4z7p+hFGyDwElE4z1EUbv6EvbF!1Pid)3MmsNo^pLuX^X8Cu62le61avt zOyH-3NuJeAi{N7MD`k^8;A)lDWB;eVCU>H&taZs0(Z}^4+Qqq~*C6D@V$ucz6cfV@55)HS zjt`mO33zpv8>jy8r<)5KGnnu-@%0O!TaqNQk5!M+_4h#3ul8on39yZchZ_?mb^s_Lh(}N!+&-?L^ z2an|^_?=?YM5E;+UC$O{P-o_Fla(<86g~b3u%}3;KXd5qdWK4_RM_8M-FBtR&HNjqY)J7gXsQ?yoE zB4_#BO07slW^-CQ`9bBxQl{LI+i>}Dnrt|xrAxjCz*M8-C(feAI0=ZqGHEo;&*8bnibWvE~&WQ7)_ z9)CXaQ$D6JvGSXzR!2Lm-W9)_r@XYfDN6^9M6MZ_DH_B>gFLcB1*|TlOlXVR!Drte zrn6KCsKk_k57BJXOmLdom`!opAS<ZTH zL{8Xpl9%w*$Ie%G$(^qOTi#n)b{!joy}0`?giz{F`QO(11BT^7RoXyVRmAb8KI+=C ze&!kspiR%-V>7$*u-jtr2PvA~%{&x&1umhSuk6cgSuV+f9%NWdr-Mhn+@2`)IvA5Q zmU(c;{q5JV;Xb`iQdj-ln%$lnNN9cN+uD7HBJh9+M#5k4O)9T}&v{T0IfuEK)?h!R zX#pvWVGf@_OFU76$i_^jZ7gPX6A#@Vzc98WlI#PWl60u#41>&BDefc5#_HvVB)fgjx~+q6 z8-gv0$W$(;2qh4DX;J{o_yO}GOSLiv5WiDG_#FEy5i(Vi*(rwAWCicGh+iIrQW;TX z0&`OCAhPqmW#14*ySQlJGJk=iOMbysdRR5`FXU(~g>ce**Oyg40YwFH4=8U~$cB&# z-{^`4OS-ca@yZz*A-?;>nk`x%Ow^gJiIh<@bIJO&d^bH(kH3v;*CIT!vVIetLs;~Y zZ{gr&VwsKeF!j16X^74m78(_!BGkg0X}zP`rY7=Oxql1=8&kop9Oz0z!jw*4S`OmJ zNR((+uz;KiJB5tOn87HJm)zrJXoyfG7*(?(-$eQ;&CpNn@2usL3IN z38>aH<9&=RnOmXPlX^eHY}OkC65eO0_so&JB-mf~Kz>Y{z6L01{z@Lpc}3C|sku`m zM3hvSqXy|XLT1;RngmkKx-twO(8!f}L}}skA%^UrMunk?@oLP)yqEbDhMC%dR)s9T$A2yY3?R z=l2U;Gy~8XUvQjDmC@877pOtVhJNzfB8R#{6>*=R-pzV>z90M-I<3w>j1^ed?u%s4~>t_viqsS|6 zH!*{AJ$U*%#n?Te3J||Q5^2dAh3UCZ1vPS<0@_Tnw^UdaHlsPw z{|_s7_!1-H)fBmjy^3un?n%)c zyto$D@qC@>TpP%AGAw6LJX@QxV&1~pfgm_rFLH(3jUDF`|HyGCyRz{N$C_4cFFB4_ zQKq}DzZA;@fdzSw<*idiVAYw&-YwFm6Tvt6;N_Eh^AxA3jV2I(c9JS)pw7sZm1KO> z$pk0rcAFJWm0 zW9iNm@wH!1+Am)I>dOqa&JV^{<{SwE-s2(a68(1cwgx+7E~|By-F%*3COF$VIrs~@ z(L6n2VNM^pKWf^s3Y`%){nU2WjN67wP^UJ1ntkBkO}d1j)B+k4cymqB#PsuOJR)td zAlA@&I2@8n_pG&gs(8neQA75%| zm++XO{lX{DV2a<@z35D91yj7Z9Nls4+;+`bc5;5b{$uhWFW!AZ2ARNeYe`lRzjkt% zj##>Y<@1r9UwCReTmifPb_deT!p5>Cz;7eJr}(8GJTb5T(CPJ)AC$s=yQ{+aPzgO# zEW2_v|M7NAZI3g??4(y>RRwIrowUizq|h!Z>1@w58ij(a$&sPZ zPsc#j(;I59sj;GGi4LB2a?5|<5U78;cEa=H^!udlG7I>H){~2ryE#=;i(kEXl8f{j zZtt|Zd>8^Zd!fs8ArrR}i$;vpXt#UqI=(N$<*b~!?N{I)lEn^yJayK~EX}?_qD6Kv6$4 zOvePD4?1R3TYZiA`cGK*Q(zO-c)a}#jstyGxxPBDAnbWhp+HMVxBKt{ z;ZjS&fOGtwT|}p|N2NF1-aP2DVwAF<#VECeEN3q}GI|l`=vrn4kXhA9d&xuF zD@jiOPEO$nw!)j#tI+>}`V@sDkF5s{Ed)|8i%6f2r9DFr|T0eX}mA$=V zo_3AC?YLOp6!!%f@lBa5T0vHfoI0wW03cfAw&*4Kdce&f3AJXvW|L_q zpm!Is#@t)YN#CtDp~531>(`fg9P`Q`?a-?)s_E{I=-d$^mY*LYfn@4Ko{F@F%`%-q z89RIlk@S38qLfDwho2KpC?M?xsRDT0rj@Y^6ZA^de|C7{;m<=3{V-x z59#p-1uSPumblF_+8BUd0DOrY*kfN%E=Wo{+=TY%viX2YwmQ6p^JRk_b=cT|VA9U- zU@D#HzQ|B1A<7x`h4M-4m*zgdTqq*gr2KJeg+`RpVHXY#D|{dm61b$B#Oq{KXm9M5* zV$*Pjf-Hr-vFjWIuTf9>oyTElMIav7Feliz%1VYvk9N$B!I*2|f6ztr?BSS3tadag z6NiR1r(W6VnHN&iOwiHHs{U~#a-}9eivaOxmb>}-GTrh=$_I_jVF}BscmfIvv)98ItFR?g6lOxNysGIb&cnP%__+` zML<(ocK(uz)Ys#UbHtysu~R*(yOI0+!VZl@Uq@ZHvT6i{=-biz*mZ`~_-;O|&E*sS zF^U-kosT0i9hV@=vj@300U3|?0dtY~@zjzKa&&i*CgC;j-ned$&_TUm$p>FZPVaqS zI>sT+5{`EMb9PS|&~%*vLa$!ww{1zmMkS$c0DI5)xE?aM^J7?12=xLioyf=wgtLVC zKvAhQkT^=dmK`J5qesJ6h$}?f?TEC}e&$zh#{k(=lD*?NAcITKRtbhxK>4>MpuJxL zJys?YzpilS;M~`EsJ#>8eb0__SOgQseNNlp`!R=+(RDQ>M{)D}knN5ZN;cfI?{s(M zz4{0E=^G3S*W<}Jl7LtLf&tCDDF$Z};Ow3F9zVUm%x{HOqlO$IHCPuPxj8K(*bIg{ z^NVtf`uyf?2ZnI>P`Ilp`aGIKTfP=6dq(+gWANB+2l3Mr$3<~`8 zxSUY%{+K^%+aNB%%a)K$Eb}yw`RyWO`wxO%-rM;oFyI!n9pJGt#Oe3@d29vRFD!5s zez9b9%O|}n@C|L`kts;bHvG)HmuZ*&btt{RNW0}D^0C#!m4WEfA2qP*^PM}~*InW^ zC@iHWr($ubmq+H^7p2{ovEziuvk1m-HkgtSqb)UkPB8cR0|mKTu><)+f*z-F^>5do z`~F;dZg=1GeHr`>oD(*9y6VMb8At$UIdVm~lqUyMOhpR0P+p(+fiDXZ@nqn(7Vy9; z`gpnIbr<-W58aLgty=s0yH^F*;lyvit3QXw!q1oWZ=kP%Z=jX8yGutWR3-ImwT=?~zx$k9RWk3~ft48VW3Ur>lqYy7S z2lGoc0HpK=jtG2ucj|Rt;G5&N|0A%V@AZlJMkDa$@(ltdJe2t_%FciBME;As^B??? zzlbKRY&>kd|4r#({fl7pZ|EmnJby7@{(tHTJ6pRli*f@81X?&275tF}3_NYP{+HgW zs7_HPY-vQm1Vl69w-(CW;7C1D4w2MgZIm(#Dhm`Do{wZ1FdxQSFc1Gh!@syh|3()! zj=w9zsKT54Q3Q7K7}H18b>XUUGU@T{6n+9CBwUPykI`9qLqmReIU|UC+e~}PH9fsNNgy%0;#HW(HT`0yGSBb$#<2K zID+AT(@F9?V-0c}#20S4ju1|7BHCn~gpxtaICg8S!BuOBPR^sxQ`!2e0m`hRkX z*jfMPO^Jj_Bq&luu`Bgmv>RNZs)(ih)j##;EZHKWU=cEY5DRS0?G4xL@bz5r5ceOh zrSGO`48&bNs+ylq6C-+VFqN%YQ~+T*@`Pf>Gno@(>1^%zyKM>~HFjz0wILn?FrV8h z%)i41RM{J!F@O*UH4HmOk%SZFYCV5w^ahHmv1I&%E&rAm#mV+JQ`k9K-UD;CRehF| zjZ1@t6RHVmS{_psv2*7_9~PaWUONt`c6Uwq z4?S1&r>FW`#hKA3X@N!bfDMVKrmVQGCcfe7A(lM7yUN^{aD|lYV>&t9c7*())~Mrj zf52yMj>=E&adLx5n@d(Y@Y7E2Z~Q-`);u`{uyfSU7O@ksRWOg9NR^5sPG$G}1ZYV0 zeOV;2&<@n*zQi}{*WMV4J|MT_H4Tc`$Q(CQ3#_{PC^0!;PosP`>r@9Ts$5*oviDxd z!U>|a?P~w>=g)3D#mnKj2)L45o`|ipA`KYRqCTfI+EAmCKr>Tp({_ki%KueFQ|=n> zic4~~79Fm7WTN#w5Ru=Jfq9K#EOiQjO|n8EpT}z~ATImPJMGx(H-z1-EB`+L|34{- z|1aR#-noDO^El#EoixD>@oF?iYQ&a)f3Lc85vp0 z{VYasY3I&ipDQ6p43oPtq&O^qU1{@b2WQ0cZ8XGJ{U#?asSaBl0JkDzMs;3zf@n>a{P z#7tX3iS;ldhq(Jm?{n`#l7?6ykud88JA>f(!J_a{=eRmXFJ^U4rZTq4?r$ro9xOLz zGWv;6+X?1~t}~B2B`6t2R2)V%jxto(OFisngHc8CLn1lqIGzd^40V{(>|>Xj+)aGV zBz9hu=zAj{#G%%(6Q_+zVjI&Ehh2y&;gg|tjaE+aRbqnaQ}A3V6l zm~VC}yS)4n&*d5hPs)55{OUxBr|#tWJ!Qy$T$iJxT6~69&B{wk;}bOTZ~yXEU``Ky z83*3}M!=-er~3!|{$eQq8~k|wo^)`=cBA&Qz>4`Q4Zy<0Bs`9n$v!xST+V(g^?VEt z;S;&(z^Bb+a#+*y!B0rZr~1%%CogunZ6vW2ql>}#)xmVugsxZknSvzU{l<{zE)i{u zco>_|_aF58w;38M7yIA+o5+v+*Jb2c5R45=a<&ijR&}$v>eF;VD-0RNMT6Z39YpB zs_hH407K52W~)kCgYQMkFEW|hpCtgVWP}J7_77c5kXxy!F13 z&1K`1Fk}zuxGq<{sY7LiubK~&r~G*0x%$^@2g5`!c?27@-7+sSm(yXJaV&btw5Dgo zi+Ey1N>$08N1`1#BxR&`VhOWDk?ZA*66aHJKoOk$jY6S&_?P`7S|1m<2at33?U#o3 z9=K9F*n=^YSY(lI#kxS}^Vh96zM;J1pNLxVo*;gC+wOg4*74}hZxX*S7HUuf*q%pOvF6BL4x$f13#M{Qa(wi0wvtk7;q=awDh!DvJV| zwB|s364?zWr7qtNJ~TG9E#vSAsC-0*a5WH2R*4a3hoov;-bG1+!&mp9BOkTz_m7|G zlsr6~1*a!*@BTrrf1A{?v2*>+zY5B74&tPU>#ouIJ2Dic5Hn&+8pnUu8;8^{>B4aW zr9YoGX%L;OnYc$?+>|L~Ot5)lRvs#?zo{TXI=sz}iL|uG=9wI!Ee{V;2pn>d<5RCN z5GU{};ah}_0OI!Zm~f2>upl_~(x!uW?fJOBbS2`oS2u5&DI=1diB5k?ValzT%ZTL7 zuQZbNq-ssX*B4qU`sl3t1HRKVCdCmzpi|UuJS*?K{1sxka=b}rj6^1!bjH8LHx)jv za>AP)bMSLY5$b*TQXx!QXIOF7Ve%{CPe>N%wpjO5MA@L~4+GnD75;q?6mbS>|2Jsw z_xZm&px)m-@HfyT5<5}flOl0f4G}oB2H3O4(K*kzD(olO7RL{a5%rX6)x^|f1qzg!469UaDR?gUSR%5lnO=&l-%HwyvLY7GJ1VO?+%3`;tS2Xi18=Rg*i!@t zInSuC6;8kf0xeZ1dWH1}ss3MDdGchGS^euxf=Y5d0?Bf5c7$-HvS`VjO&s7a9&tZd z^lEjC>_cxtaC`rQmH#^R}&U%5L@?EM-OgvN$7`pw{*?d1gO!y zwb1I(o&1zQZ;(C+N>Sk0We^riHkm3+Re5A-jlXklN|JJ)9be1kWz1UpRdFfr_gR}N zOXZL@@0d-9y6RS6(f?O7;oL!SnZ`@*@2#Zt8?x5*qNP(UN|mg*rHaH3fw$NK(4UX- zg*yN7K90W^_w4^VJ7oKNjBqOPH+>{UJm(E#O2mz|Hv%?isWU%17ML|5q5YA*xf7S7 zGCUfCtyu%LVc|yNg0a*<^@5i}p=W{}5wM4QqM61Nxt!-6ITicEeFlf@%LCm=y(a0{ zP;D$J_i4jflO6-wc+}2{P57oO;&!G6=cgFUsPZ%)N%JmQD8~peY5ki@VtM%#52B7u zs)w2Su8jQ$%2|sC3ACVW@q2h99L)G#aswyJCjZz|xMpaZ%DqFsC$7KF-2}{G*OA9ygcLo5eUXtat+u z*y*bA2^ZwJvYqJAy9zgk*q@Q?-9BnY5H;$N%$6Fhw2`=se#Su@$U=`HLzGp3d8>7K zRNc4|{Qp|H8W=08E4&DnI4!Y>P_vCO*8-xjyL0D%2BmSAx9cKG<8Fi+BxQI%yW_r} z?YvnR(O819>XOjfwA(^r)6m5RX=qw#ZH*es&$?vMCe^l35}UYS)wPk5wyZwKcOPM!PP_-~qDTDNLf=Eaw%hWftx z+)63a-|^sup};-Qn%npDhJ$Z!o3eWE3>vq8c-_+z`o0@hU;h0sUz?o! zKw#>0RQysx=LuEk-P|f4C?5CS2zwrDcTf0Y&?q70hV&~;A9kJQgTckunY zmc4NAqQMVNJaFh&FFm_P-FUHm==FtfOdi?&_rdR}2VTGJ;oaRc8I=s+!;Q$oXc`U-$RZEiI$pp4`;F9;*}`n>OU z17@hXt5R~ylH6Vf>JK3w)$xfUmLzAqLLpCZL;y3|$aIs*{dwQ3gOg~^pN zAcKOz3Nth&7rck^%rTfeP|oH3YPRCz7S8E1$6<1o$vGzHC55Ag#}p1BMJAV+T!tLT z2Ii2?WZ)51AreAxk7k4Lba^cq0B6 z-JUALngR4vHA|>bVo}Sj;SH^8pknoR=DnU?fL&oPCU*z<(t5_F8oX1gd-=1d!pQ&# z(54MRj?F9!_n2e@om)hBX{p0cR@)gz;R6_Jx6Xb@#l$!_9WA> z9h4!DN--GClBgslXi!8Z9Ff`@7;y>f^RA+1{n6MqhCXX zVH|3wGXtd!V^JwcrFy0?0+rOim^3mVbbvvHu?(~kMv0W8l6+8Tqtgv16Or09G%7)2 zKB&8pBm+YfEy=$13n&a3$Hx4TUm#Y-De{S^B%@&7B~j9KjKo}r&#BHdkzWYEL%C@>IK=l? z2=%nhk)z-*LQa;cwF*Uea9E+1&%}NZGRh1-C1fB-*HjWoO%(zqIS(fdD-0qJ zg;xVi4aZjxj){;D2Mg7yKaVH!RyCe`r*RGCtjiQ6gI4@f`fF7_rW(Ckp_;Y$9GihU zGq(GzGvix5+8#8Qyz&6XpB_{c@asE;5VjFmg9AupPzXK9zuMYXPYo%s-6un$28BuUwug)?o3Q7eWrEcIEtk8<6Ong zDp?PZku)iZwu+PMIBDNdY|TvThyOO?z83hdSIqkYrM27JhJc4H o>6&r@K^47HwC<`@XtmQd7I-yq{Xi@wRZ$F8TDr9TOC8dG0On&EhX4Qo literal 0 HcmV?d00001 diff --git a/ephys_analysis.py b/ephys_analysis.py index d627a3f..13d6aa4 100644 --- a/ephys_analysis.py +++ b/ephys_analysis.py @@ -2,6 +2,8 @@ import os import numpy as np import matplotlib.pyplot as plt import nixio as nix +import glob +import os from IPython import embed def analyze_sams(filename): @@ -13,37 +15,39 @@ def analyze_sams(filename): for g in b.groups: if 'sam' in g.name.lower(): # go through loop, until 'sam' is found break - rtag_data = g.tags[0] - rtag_data.metadata.pprint(max_depth=-1) - # print(40*'*') - stim_tag = g.multi_tags[0] - stim_pos = stim_tag.positions[:] # beginnings of stimulations - stim_extent = stim_tag.extents[:] # duration of stimulations - - for r in rtag_data.references: - print(r.name, r.type) - # embed() - # exit() - voltage_trace = rtag_data.references['V-1'] + rtag_data = g.tags[0] + rtag_data.metadata.pprint(max_depth=-1) + # print(40*'*') + stim_tag = g.multi_tags[0] + stim_pos = stim_tag.positions[:] # beginnings of stimulations + stim_extent = stim_tag.extents[:] # duration of stimulations + # embed() + # exit() + for r in rtag_data.references: + print(r.name, r.type) + # embed() + # exit() + voltage_trace = rtag_data.references['V-1'] - spike_data = [] + spike_data = [] - for idx in range(len(stim_pos)): - spike_data.append(stim_tag.retrieve_data(idx, 'Spikes-1')[:]) + for idx in range(len(stim_pos)): + spike_data.append(stim_tag.retrieve_data(idx, 'Spikes-1')[:]) - # embed() - # exit() - dims = voltage_trace.dimensions[0].axis(len(voltage_trace)) + # embed() + # exit() + dims = voltage_trace.dimensions[0].axis(len(voltage_trace)) - # # f.close() - # embed() + # # f.close() + # embed() # exit() return stim_pos, stim_extent, spike_data, voltage_trace, dims if __name__ == '__main__': data_dir = '/home/lisa/data' - data_set = '2019-06-24-aa-invivo-1' + os.chdir(data_dir) + data_set = glob.glob('2019-*') stim_pos, stim_extent, spike_data, voltage_trace, dims = analyze_sams(os.path.join(data_dir, data_set, data_set + '.nix')) print(len(spike_data)) @@ -59,5 +63,5 @@ if __name__ == '__main__': for k in range(len(spike_data[j])): print('current spike time is: ', spike_data[j][k]) plt.plot(spike_data[j][k], -40, 'o', color='k') - plt.ylim([-100, -30]) + # plt.ylim([-45, 0]) plt.show() \ No newline at end of file diff --git a/phase_sweep.dat.pdf b/phase_sweep.dat.pdf new file mode 100644 index 0000000000000000000000000000000000000000..7fbd83c25ef7de6ddf5a8dbbacd860f2dad61d7b GIT binary patch literal 23876 zcmafaV{m0r*JgLr>Dabyvt!%ljm;Z(?A+Msm>t`;I<}p3(9w-2?>FDnOifMI{93in zUVGs=XV>1V_IW6k#U&Y8nKMASV$aAZQuk1NSM{VoB$-uVkT}Tc8-=L z%*rN~09O*WPXkpFAt5+`gZY09dH#^~ig{A@UXD6a*`N)rfT8>_zeC(tNm{z{sRCBvns&V(cQ)D6N>+o zuKEekXO^F~t|Yww(cze-TpZnn$#C(_)iI9{IWFfKZ?!Ha4T|A|?`Vm+A!^!8q z`=!TI18_8cTut2bq$%uwuOj4+{_%DUdcL^@U)6M{fJ?XivuypgJZyw`J?{?Pu~$ICshyUhOY)Wkinn~ztj`yj7nV48uj zPe#M=+xk0$XJePw*)#8Ym1jUBFpB@mAxtZ;cpHQ@;c>)k2DfFFjm~}ga$DUki|Me=KO|bct}7gL?ojIL4zkp}5eRw#5ZZL) zSuU%T%+8;;5#Majxb9Is&F(i0Db)avz_F}*0ikjDSNy_zTzRpD&d@db|9;xOy7d&p z3;m9n!q5G=yrghsc_|Gtvdg2sR7u;xbX0}BbWYMn-1PC^@}3sFS>ViVpA;ZDZ{Xl6 zY;KhhIL#DYH?G{rx^><5s9?_Z8TmIvRXu|$=lSoOCqe9Mmi@-e{XN@D1IQi1;3mow zM8)!8zDZl3iyJ~PM} zs19{RPC?szlDl7t?1I{XbH1>*mV!JdI=ja+6q!x>+LF2i_2{D&>9dVlgHOK!y{%n> zH`jD&$#9VIlt>Cu7YXp?@VA)EZ;9MqrGDVv(f-?@v>N2IK_2pJNcExHlD(7SHak>% zAtlRc_*h)hvm0P1GGm5yt!RzFbMShxS-PWEENTIAKiRLpvh+_C+CPzP*kD^j6#%73 zz={imV`0(8YZ`m$b}ASWp3kc_HQ~Q3>`@@~)g_RYI=d638%XB%{v0{W@&y@gRA@diSPj(T=ph?n=69C9nN; zz6_ixoz{HKQ(i1QgHAv?jnb~2UJSdE=o*0^o?l)m^i6uvLFeUP+LER(S6@R+Sx%7c zE8JWK9uN84P4kbKI@j=kNU^`+o{E}vsVRlT=d~Z-4At2Le(N+%@S^Rb{;h3O*QANp zsHaw-jFI%%lWz-s#o;QjncI}hbT^GVPBBRkfY6ldba6ed_lSn98>RU9>2r}`)zu~8 z5Ds({8-DjAz;|{qsbPX3Q9wPTLHHq8?8$p77^Cc9HyR)Fqjx`S(_BJvqiD$*xV%7L zr)xFsb-QKa5`HIpa=<%6C*}0|vY|^vr^WV1+ZwAQZDkeR(UGS*mRk{c?_$|KHMrq9 z9qoP?k}cBSsHij{}`3=7Mm zRi!dfQB+4{78#+R-RAJSul2y#V^>O zO0aWW0z;yE7el<;G8(VzePP-q&Nv@Ax_*x(meuko8(9%)p9nlCie{@5>!2Y@{9U)^P#I2w zHNkhI?lZ*8mXaGp3zmHVOH5tb1w=yNuL)Nh9g$u?pAay&6XG~c7FYk1xU zfT4?IO(rCC2Yo-iB(qrGUl*w9qi~g@JJVemhnBaFtJgGlm6gg{d#Cg7X za11$t9c=kVe#ceE43?5L%O>7Lh7>i8r*hcJg7(qMEeYb|uAHG-sOa#6FgT)ofn=r* zUMdR*oCCy_K?zD3X;UJg%Ymce7VL!jE4wpaMG>b`wwgjv8^+71CJPYL?}stmH+l8Q zT$(<=e80|}R@kG@)nda3qk*i!G`llJhF^P38LUcKI{u#UCM3I=?+B9OKWgqvWkY*= z22UysaIebrI(LiSdU;89`Xr|6H}Tl6rS$IqEY#I)^~UZ|0tHtUCC#TUG!qMNz^jlE z#O{3c!M{Y;Vn4OVIbuUla>l-Jw+hWCjy21X;vGHe!#&-ir;x}5gqV*|H9Po5b`P?J zm7iVhCM26{&4O9ly3T|)jLFCBZR++$B-E94Opvnjm=9EZ1iuOVMS6RSe{i!&S~Ws2 z*DGU~=x5%fa8PBdHhwUz zJ7rIy`3ktXIZdh(`wp4?ddu&X%FG6*0G-P*xjaJ)?Z7c^UP04PeKg!}q8r_0R0DW! zRo2`@YYhu2hSZu5-RB*x@^+py6da(N0p5yJm}wsMdIQ9#c+E&*nESeB>(jX6#q!hA z+^1{c0Dk|Fj$!8*oByDdcBt}_-WqKxU3UAU0^IgCm zLHH=9Z5vjVocT8abzNd)M^k>3-o~!5<3i3l;kr&$dd8xH(wg9R35o4+WK{e~r`GEN z^c~WHEOh*{lELTmL6}Ptczi-b<)rpV&nv-2*{*gaeJf|3hRK*`>&6{R&x^k_c^NTmE+iR;n$DrX6 zpQ-oRd?)o@lYR;l1n0y%B!Bo z>_dV*`d=C1nD9wPkq3;B$Q-x}f_ zl`=hn``M3gfv!UPZbpT27cipud?ef7+=bODdmBV$W;SEO;bSeuF^n&Dw+{Ko6oBS= zPLN|dy;7tdzJ$}in69!bD7a1HNY@G!#brF6D+QKdA&ViG=F^=5cxLYFV4L)Y?C0>b z!}#k_FXx|OrAW41fF|D(AUIG|puvvp1S8$ysZoWfSCFjoz6EFV{Cs0+Sq}~5crKt{ z%h<$yV7Bcx<)qdRR{u#)(Xns%j%antYEzw5@}@j_a+wU*jD}*rYO31XcUkJd(tYM) zzvOfG6Dm)T%q#M!5;qdk)P+qn&CoL+h0Tr3sWn+AQ_5YDu?IJxzhOfIV;f3*ijy$l zZ12N~lhRl!1>ItD3HZWyOJR=x>?3fG8HC9twL=BBc8vSSv{D=h6P!d! zZD3}6aoA&oPr^F`;B4V)y-LfW;m0%C;oQt{0Kj5UKOEbgd*6JY#Oi9Uq-@bFRuE9h zLeBBFE!cs2jGZ2Db`#e$jnK}Bmzvc#v}gtTI~VfhF8q?XbfG$)+R`xbKm1S38ZR?nuvgRqbGv{)fygl-2V^ z7&0P+8h%)Lf-~SHEir|$DP~>S4a17uY41-@-Po0ZSy??P4rQ=jJemhQRnGd%DV889 z2G*oTHVBHVyxJE;cIJdfWbRf>YU##ammT0<(z}Y~lm7Mf2 z;=4bQ@1)Y+Ak=bPQW{YE>kaxl$4vYnf;dT6Zrt_s{;hk&zL*}$mG((9qGfGKJqxLa zgt>+Bcw!^UI)aw29+a~oZW5&_Dxb?Jk7DiwkrEKIFTc?f9anlTT#A5x2`~tmcSHro zaXuLysOhKM4f~cq$%6D}!e8~<;M2fWJZ~mWAO#i8BJjb61H(6)Z_5?hb>0WYcSsgS zlX;=#GSE?0*T_ZLd79Z{SSTV zwoghNKn6b{^db(0Av*QtiW4Xy84($#khMh&1-X=Yr31_he>^Ibi41ZGv7Z$y@_yK4Frx?jGObIB4IVkX(VtmNEySQFtNwv_Jr$GqP2_Ie&8~~%`*^W#8(&twq2k% z#3gygok;-USZH6ridz50uc*S$Eg*ZPMfiOwGOd0_)m7kkxK%PkbD!YV@@wzvnUV~` z`FmcOpFvLw4BT;qTk_io%cO*Q2qr&(p$84MS)FIBDM)s~(~sa=vh*u$b|on-vbSFH zzevm#3Zm*=Wnd*z^jtj%Cil2Epx3VzIs$)*HJf;19P8(-iQYT_a53KOFaj9o-AwFX=&NX-tl2mAF$ zl+eZyd5TQr>`O6~B+*AIp(RlD#!tGE6@hKez%(+58&5qt6@r){5u&iqN^f99>DQgE zUDzi2K#}};p7&Vpq%%0VAq2z9fK8oWl0kGa*6iJy-ke1sQBN?Fc%NEIw*1$Z6Nezq zl_L3*zo*oZY3z=qEJT;T$@~sDl8I#nJqCo%bmF5IC^}IBKZIgDziIXlF*$p8#YBXu ziRUNGV0FO(d3Qu7nPE4V548v=G2J_4xeBIC3+I2+LQ$B-R^0q{Ls`xGE01%XfjiU% zmMN9K-)qd6LkVbAQszb@--PUmhAj3JWt;P}7Mwb{83+j{y%bHRw%X7Sdzn&a|I@GR z>|r949nsb34SU1^t%MZHF#a9Uft|pY65s#eYNwY0SbJj+p{D)!+*o}_gnpD8JvQz| ze!FyhzRSx_6kA?JbO;DK~n(!|7 z=Unm}oG(R&HP~mDdeZaBJuITOm)3E0rK9mAXYxl`pcyc83I&cY#4xrd!P;so$~W8& zX_ZnflShz{F!B1`Zk&*ljkVQ&0DM`l3oO@2eolRJp2b4)+wK^ar11<{)Sjxdesz>G zF#Vn#gV9Z3RF1o{Zu&$qp;45$B{j0hxy3G%tmuow%_^vRWWm3n+N3klF~zms=Loy4 zKuxz`$oe~wGgOWl`;d(G^h-h8Pl$7Q^c$CudVtMlSes z%s%5*iCUijXrR0iRcnU_8Wf>}!=e{;khcrC7|$t0b1QH54^T_3E?N-sDa9S6xk7i= z3aQ#I`O;8Y@lA!c?;lsh_G-g5P>F(6FG!~jnE2``c&B`*1Z+m3@(ZSXr0Zb(iTnp0 zYs`AC!P>23BUAos?NR#P;Y6-KWp>Rn^nbaE4uzigP>Nj7bJm$aMqci@} zJx&-Z9TIxtPftOaB1bd_Otl3@?#9D8pgdap(t9n8{uxn7t|QAPOB_i3CZy)1C!*4> zSiPeiw=lL-x?-|M5hf?;rsov0r;iFv>rMak;%TZ$n^!Q*FNwa(VvPDsi&tHy`*D$g&9+uUQ1(=+|u-bM6Z z4@w=KYT?d);yHrQlwti&2~9qVcSL*uU9=p@#VGw($MLu8c2biOYN`CKE;JBzDm%DGTX7io5rZ?W!dJJ}=wFNS52;Zz5-v=hnE(OW7Dn0CeA7su=FDf<-1aeKO?wt-Ms*C@QXsN~o^w>h+f_+@(lU+P5Rl!zD=WZ&(_; z;MZ_W$E!1PhafmgqrVekY8k;K+~l3{X1>wE#@#TiH#V>BsaX4l8c?H<*ed8dEIXXR z{TT=toPp06<0+Xpr*!Lz6tp~6?{;D+>PqlKCYIqhq(~4$Xi+B7#<0KvMY0 zA(1)0VMBm?3+J*2NS!-{<&3jiUp?Aa&CLI^#S`)wm8MJl;-3O-bE<~*RIZ-M>5Ot_ zBSuO?0N<_O7nA`@MV4t~k%a1Sp zR*?cbSFLkgCw|#2O_?`;``0q&9QkHl3e~1>L`6~FI$Ckspv))tJ~k|X5wRA19_fiT z;3%J4ef|S>M>LKyzLg*%8aa3SoxhY(u#$sgJ;ULxGDimki@8@q{CTL(NN_F{Of|{i zQipue1B8^sU?s{BrxLly96z4%r===e_&l9h!krQ4bf!bHr;AaLn~G+cW~fjvDB5wp zK16gHa*qE`z=E9S91@)g^Tn5i_&pf>@an9+Uo%qX#>X%QD)DVnF}&eN;vR z(=>JbDF#ZYVv}a=~}2*i;YmT zxcHlBfv$c!>i0>2Z^q=t1xh+u8eEaNj_@K1W7Q0cP{RwIQB?c@N2i1ZSKnZKp9Pr4 zdC+?ZU2&75^R9>u_e&I`on-hv17d);WjxcwuX;_&9H6OM(06*8YB>ONTc1g1C4eR; zUKg|Dvu_X-r@K^7#(j<0R+hrD)KeuGv8llEEPfKU#oZcn99^FzF}d0+l+{1Gp`k@Z zcv18aD|*lMD&0NANT|0Pl8ZJU{a*0Hv%0y7DO#1A7NgZ9d(yg-*FD+co^kJ3N1N6g zJ&ZeU;z5mxK^&`oA?d1U0yS=ddiH0i3QOmpLl6KX)a*-(s#KAfn6PvTl&^$6$)Jca zf?72rG-!Kq%mhbx@g14VAvyryJN`aXRHEiZb#2ohOw1N8YuMNE9oeg7CI9vsU>eM{ zt`7S%eDMmd0L5-&3o-GifRlVbVB8Lpd2#!nwvrxA$0I`(u~qARORY7CjYXDzg^8QjX`u z66d&PiirgcHbZpR6h|1_Xv=pgNasZ&aKhGwPh6`nMSRiQ5F3uJu1M8|(!jJS5p0bd zHXfDAJ~C_VsDypTlLg@AlCkP6vT0Cd`bCf9Cio#Ct7~ZPn5~Jy?PDEjds4ez9?Z5s zN!eTOyl`UO$n91fWQB$6Oxo`W#cmj~Z%vt$c_FMWc39X5kTaxr?I6@+>!4LLkxhEu zqDtUwNDPJ{rYL+-LQSt?KcfqEJGZDvDC$t)rMk|KS{G;sOU;AX0!d}rK@h$>FZb99 zcFgIlMjw@o3iFG_rn3X?=_Lk1&XB2mTy>RsNGJM1Vc)8xw!d8-* z4)|%FC>qUxBKvm^hQ?;D*}!9E@J&{SXRpI+nB3^0a0q%7-qIpT753kr4(fWUFAJEn z4_2oh>rd6ZP0*Km_Hw;o8a%QLN51%>Kl7;%*Uv>hjFC_u4k8DRx~lU*5=AGj1u_@| zo0qR;MElQ`hqz)5i~VZ%ho)9sQq)r7?{;I51|}q$27l~1oukTx8kL8=?zjDOW+9fn zGOw>resF(tnn3|}mCTRz&pnbx6NOZDPh3@2NK`7_rzaF<6j0E!qgzmz52S47!>}g@ z7h&7d7l_trUqnp`M(MjZv3Ia9IEj%JMwJ}Mah%J~Z50r5u2Z!X(q*y#Wm{d5gwrp* zIjl$ln~@@r2oQX;^q5Q7(v&IcGN3*$C@}#>`_vMfZQoZ5$zC|IwcCEfPw-ShCZZRR z2&J=v)zS_3Hiz4qVezk2|H3-LELo2nYIiz3>moI2kK^b7*M7>oGDMEgq}cB_t-?%s z?P0)3>DOgot?-!zap~OM+#846T?zI;-<8}wIb#QGQZk9);PkdaGn{SynTh{F%)A7X ztSGe|cr}bZ^p?L%#A>k_Zs+Spg}nAt6)7vp{@*i{$iyIbTh(^~5RpmL7Ta%4X43wZ(KP!Z zoF8pfmT#C%=39V@iPkmd%t{QmlxkT3@n-7vEPsg=R1>QLrY_V{Wd_V2+q<2CjyZN7 z`p-^~=}wX2t4hCHs^>k|3=%J6?sEKJuOnRY^_eg&i$6kX_T!9l1UWS%474R8wLuQ# z%wM5jaQny>T*kcd918GoDZ1C`f0Aj&eKRo75EVn|_P}meiL2Z4S5}~{5P`G;dmio% zFisA>$dG(4h)@;SWd+Ra_{geiCDIKX4w+1J%7VSMn8`by?^n;i$=2L(%cFIMK}iBT zCU@otR+!E>p>Ywhp`1d4Tyy-F{TViUk!j)9M>edf&?aVD~<$1n2i`~lmt78y!a|lxJ zlEm~q&+MV(xh#!-UG^{3=7TxopLN+Ew z)<*@s-KK_<71tiKK>1l@bZb4dsOUvhds6b?l6z3Hiy|2g_lkQl$n-n=bbN$_6w~|{?~0qcn!Et-1-Dyaho8g zoJ%a-IfRh7+H7fXZkCyE_tKF(fFccQ>(VZa+TfWN&G953+L$!95Y|<2idV}6DFay<1-kGJtDsd|D(yD z=hm(DsEX*5>fP@xTwB5rmx z?_O3Za}|eP8Vw(x+e((C)l2UzzfC0HiO}C`o<5AHeHE4^`pIk=FbhlZg{XUzw6RR+ zHyrl8XkoN30UVRHy5Mj}c%3@goPBk;QQ9h+k)X5btrw0TWmwWk^rFCU2b{f{&9E{r zEO_B>28S3z0B*PLL8@z9MVz#R^1&B?Y#qM-z~xMSsY0aK^&4M!eeuOcxKRry`yEyG zoO)DFd~1?uV?kI4Re7l{i=%59E_n&e-8(+&b4-|Phai59OO9V0- z!}@Qqhx}-xyQ zAVv$3QmVJ3K`bP$+N1bTJpA`k@(x(SqpJ8=!X|=3pGD(?Dv@llKN)Rp7ts-WU3F!d z&>d3v(y`Yd#;-u+<+efM{hA3%s2)1~-N9i})kgkT{{FWNnDKOz zjdwiiYgy$KTC4$W&+klrd?~Bk!=)WjWirN^DY-egnq#}=0XbT{4Y_9XMY;}(0>BfG zaz=Lo&vMh277o3GNbA3F`;9~bf^pDG_=>+cpk=f)*5;_DCvm#oiw_tlUcPUICNYlL zA$Yx1Gy;r#$s!XeN5)i@Y_ghf8Vt6S{%$F+)TzmzI5UWQikdfF1S)B*S?D~e=o0wS z@-Nh9cP$}MJurBslF4Sy+W(X(?07Cd_##MOHr>i?jaB9*5RBm34(*~kZ!<|YcEz~< zB4diMwS0qreFnF|TApRqExaTp;T{7R2$xIW!qQk^yJt-fzTjDPsmqRu96Gnk;1^R> z7m%{cwG!j>*MDQ1pGwMF8!PZ3k{ia*2>#Y-V6|?IRO5ouTceHiK3hU=xI9Bp@T1sX zsLs19FiTspxh~h9UaX-QO=8uIxSm#1WH+d^P8o9WehTPW%&6iNWz~L@FON$4AS-;< z;olYL+K7xrdB^AF2C_D3*46&1gB!$wL|hOmdM^y@m`$xHkQ^z!A^!&{v;wgZ>E#V?}m&1%`9Ub5IJn3-Iu2#Jx>Wzv1d^Vlg z&vy|_DOu}sTvlJ?*n$J3RRi$5g~RN4gbXthNCBek>X}W&yTv=l&RA>gqX4Sue5HU0 zxa3p_v(9cA7jT^cxoL0Kxg#r`tEV%~-Oo|^rNOxN*Ji613TkT34+1 zO^39rN{XocX6-2^f8VUtT{b0UpL@;Sdu2*dUkR(<^BdlT7L9{*!UI9`X>1jq>$qaS27?reIFuY zT7!B@{S#+ooHI`04D~L`rYUhRi%8TRQ}fXbXUl`u>>4}I3chWTktEA-hmU)=@4%HR zn~Eu|3@m4T^UQ!~TN@Dkl{y@h&xqomgIMrEst%Hduy%UF>HDplH(X+)CM|bW?j4`| zgseVTFmFt+33>Ll2?HWEf7*EGGVUnHFe#4wyR|b~#dkd|aX2ttzbYHd;Q;3{zKS0E zfJjE-&0z>~=}Gk-WUW9E7dobhV5jDot7;w#2V<^=A&>n76=n4B4bS)fnJj7k`fv@i z`=1?IX=CUd(8WJXX<%^2QJ`rpbO#8T0^gSka|?y8vSj39e?%;!kdF<<`lV%8>9{RX zTN-4L;Juuwl)`{H2S%xIr)o|=clA(JHtyd{SWd_#Q+GGbBjR*=fJR>QHaZh{{S>$Y z5^MTE$h98tv;mB=L>yqWK8uU1&)yH}t~ZOW`-?qH$Wo~Or(D7Lb}O%4K`s%fH`DUR zju{V60_TZvmRt5Pq?XV6WDC*P-xWbVe zl-DOocC)=(g?xASn^{~Eyvp13)gN{tyKOyuP>;U(?I};|<|qU@1r+VMSecVho0CM$ zV-6mJMK-@iZCY=?t_)ax?`<~rSwnI;wb{)$F-~FpYAb0U#}|3)7pv$=zMtu+r~Ys8 zG&w&0&V?OuJhb${7nE7kD;URRDH;;;;Q4lLQ;=Y8^^^3P>*LoP?VTNJ9P0QAJ)o2$ zg=ipCTCOqHc#6||vqelnM(v@(^2KWo0j8oP`CJvzh})7HE>E?VDTr zxZtu#gB?6>^+mtt*K9|o~vr-BZeKZPhGmd^nHMm@7#0%7>g}}MHdPHVQ|DDGC zm`oMq>Zw&)4fby@C8hH7TUMNB9&^(6K0$g9h?hQL&NHP$#vxIvPtw%Ke}{?rT2H41 zhAV9Wr`{k-RQ6v68?nOl#}rm*+QPgtpr5vJ>ButWVV`R-U({R&w3{MPhoA&9qYe|1 zRfMFiqytk;Q=zfZL#fxD^+rS8tBj6JjN_4T6@XUS^`sPz>X|cMO!AEJDOL(SN7<9*`+9P*Z<%mgGC@8x zjPq1le|({@xs#jgndR)4ck~#4?D(|LVZXMDw6P)nqtn{X_$}0ZL_A2vSZ7kdDZ5qM zmCxf2bMk-{ZUXh?Fn8W5sc#xDsL*C3H(3#*#JqFC|NA{toDGX@TwHE=Pevo}qqELw z62zmjE?Mlc6%GUVS@xd6TZ7m)dPFL*5p-MVNx4{MhkRHyrXACh-#A(OBj=D6Hc=|* zH+!e;}&GX!|$0ga=-S0 zN0?ysTO|+i-hl1<3BNI@8f%G~;hBFM!LK%9=-;5&1O!9N6d9ISUQ6auW4;hy(e{&a zin%H(6w-HQjIs2vPuEYXraD`D!?`x)z~B@G38Z+K!`F-#_in4NxfO$8)`OKi83u0p zE0wYkyN3Y4b&0~pOC_<~=$#XIXJ~40maJT>7XmNij**%*`#NWm$WJ3r2J~7z#c%}4 z-&xalcRUzgU-OtO+UVHjC-|hzddeHQ&0cJA1Jta`OX=qldp8Lu=wmmw=D{d$V{2BD zj)ePi;^B&kPUV=*1a7uhH79fsCnLn4-Z67EBbmW%)iX0Qo-~)7(y>`{$`v1H<-I*6 zi$q7wZXcJgnX=JWZfgU}%j&v8GbUJFjPz1g2R|X{ZftMfJ~!!_59K5x1|NY9V*2sA zFXZoGmp$9&y83~A`W+^(u_eX2R68~2EzcR>1Vj=%W^Hwo_P@>O5mFJ+{SiC;lePHu z9cCVlOjt|PKKTV^o~VAt;M8d6`8i|Nm3Ndp`V_;%ovfK)UvDsMl;-9xlLr9M6paZq~l_t4!6 zS~NR{2PZeUfU*n3?H6qg^pRdcH!ga_mg9z>J|}idJfWbidqW;|B4-xb->UKOon2di zcN^m4-Fb!a5RG65OvYPpGfP(sg3#A>;A6{AO(b|NCyb3?%SCF{`(mi;Q=+`$;hWp) z?aLbz98}la`;Hfd+^GT*@A}G{b1DD?qz~PekQLB{6Yd^;T}_}jtx^R~RMO9K?@=Ld z>M`neE8H2BH)IVblVKR(b)v>J49AKy1$IedTWubY>ggS`&K>iK011l2Ro#ngLmef`WCG}$0Gf*G39rBh|It4TP( z9*>GLHjZx+&ui|~Pkp~Qmixoc6}^vKs=U^KW>)Oq3An1y0LNKkSqlZ#%DK# zoMuX=PKlu4o6@t6hV1yM>Isocl43k5BohHFF>Kz;M%1=0E;~yMd-+{cU$#RJNoHZl(l3{ur6*F@OVQcF+z>1O2hU=Tr^d1UV zlw8plZ=VLshP`x(56-QCR|u5+$cJiAxUbo6%fMF4Ig1?KYsen#1cxRZ}}Uf;oDoSF1XF9SB%EjIctC)JRp?&e!XP*f1?zSpAzryFsq7oF?RgGm{9yWCb^_%pr=LjB% ziRDM2t10ee>OiE)mrpg};_?aEi0F&_F6xw=E*VE#x%h<}VzHl9-{=r@4Kodm)$$OETfd*}@2*2r>QHlXwtCI3 z|GhK!0F_Tet*BZ35iVNFGEvH(diSlDMeZ4g`P0xhGwR0vW0Mr?V&X!AOfci>30B6a zb}8gMOKhUlXRo*<0#3&kdvq^h%NKdqU`&;+@2U5QLqMN9i}KfUNw7uR`4dcx8vR3v z1sey0Y#h=m-Cz7R`z7Yq)jEttu(KZnyUJ6E)Fsc8<9(^J2!D8ccc5TFDpsYaWglC$ zNh7aA0l_-)c;ukr$#|B0FY87M?)X@Qf*DS@H&Y|jY(T&P<|`j-NZ|~Iq=9gc#h*P9 zdrqR*?EWGvIOx>fwL3Fkdgf#WP3MD=j#$eDRoJ;P>_c9X@rwl7Q;J;1X%XG^0W1DG zPrpLCkR8Or;v7ctne{1x=9sly`Jhw$0P;es^2mrFOTVsNolW1+^Jg6X_HAPD7^>Ag zDX~RKk4_`6!VgNVo%NJDO3O}N)}v)yaVslm^4wps?ZVS#l;s_Fqy|ni1;x=Usq$f7 zDcn)1;=BKjA*B4pD3Zih%lJ4?ky@jtet5b5YTw;oKyek@>Re6OZ2adBTi3kDq8V0) zDYdKjbCY-yB2N(!5H2&?&+>jW|1V*}7+JOb^H`TEx_VGzu-1A!d&GmKOO4=zVE>SZ zoC1M_*Sm<`xU!tzX4A_V4PbxVS_fKTU-9CRQ}{`sZJ4Q?^k(U~>+NLz$V!r_e8vO)(L(bs2!PMYYkCtG6T0;^9lnC zu%1Os6cX__G{NfZ?XQ#oDo3r9Qv|gA+$Tu^V=wFne8FgP@jzQIEl=vuQR=C>nT{=T z>qAOSjgc6J70uspqzCI4Bneq+H{T0}wYRKS3N7o3?PuOUFC!345(h@EO-Y%`y;ykb zs#R@y2PsdH%nsL`Xy{)EInY8IE$_cLI)QSa!(VBUGV>{DalGxoce$k4WOPL7jS{5{NIvw33-!DWOpHkAyp0+24l)4I;sXb5K z-TTKLGg?!k7!rxO)5hwuhJ3=!Tj!~TqxcU=#L^T1>P9xuSxN6|u?Iza)z8r)-uj%i zF)A6!wKz*vI?sQ_fa4v&wijynt>Mft5O7kRMC>k%%Q?U_12^{yUS)@f6W($k9Rn`w zxKpY_hTF3HnPl^|HW&ZrFG~N+6)D`CjB&Hh%7*r%g#x4ok#~M00b#{1KO`kDbPL+B zxN!6YO4`e6VhT#|aj>8#*9f*jEG*H%sv1ee_rdqeaPM|If4=0JYsUTK%Cz=D6ZP77 zkXc7dIP0)m7vlaDfH9ds8OJG>KWealyX3iEBvnzOK7P!vSVS?P!HT!`P7!fn;7y3S zM~;N>n*@q$U;v_2g%#doc&c?Z-a8#}>$g8gQk?1G|D2Eu?(ee6EyHFn=KE~UmQW~DnY|*C7=ubf zGs3cyJ?~K+caw1Y+$lqf293l(Fd}rAV&{?g-aOrfHk$lkXEEd;SqD+jG!55`=&zw_ zS`7rDlr!E8^FMCp#>k~PysPUWVI|I}X(D+W%k9dyxlc!TuFZ^>V-akP&=tMPQ=a$t z+-Z6iY{^cbP+*cZdzb*ugvCaZbwf{!{JaVbL+i+~IsV1YjV#Sa9vUY34eg8IYh zcN9%zI?4iBO??1eOs$EyO*j3dyat`TFCE@Hyl%?t7^~5uAtYkynB_lLI+K}_KY*?gtiOfckE2g7`p9a zX-=Fc|I%j@I~WusCqplO-(v&2B}QnPgOxbD4dz8vBwY>OeyME86l8a5zwx}pR!R_K zEfYcd=rcO&!rXBVXyF#*G<@N3aND>ITglCq+(9^3AQowAVHY&dBZVsT^Ye2l8~fVI z%%9Qtt@$iC(7Y_P?Pz&)u^FcOa@i@Qzuol<)suZfLK1s6hX^jX%2n>q{!75?sPy%5QX{X*ulsh_ zvEOcn^pBXHOE=7s1PY=b!w)VF-Jl2VH zqo}E{8*bCo(b=4+y+9>&jA5Qog_Vg`shKnQaqhOK1B{FabHN{v`^T_^^eA^l4P75t z?|*yzUT#4F;7^uk3&`m0ve55!Rq73V4Ep!~qKthiei|B?d8dCq*~boe zv!^nMVE=P>s@)#TIco0zx>wWvMqPsbNqRMUJ=*Vi-8=UJhF$e6-o8QlB7c}D#Shi( zwp&r^{>1`+mS~4+N}`vuITul};!Ad0!x3$TJMRf&+C2gk0?kU!a!h zAizlFmv;+(po6aGE3ef?{yrk@5_Ww8b5lb%c&iTxo_e$Koh!P}BcY}=KkXCsKsoHs?2nS0g+`=>M=C03PXiE!=*8zIuPW5_-Q2_=Fss7qHLh z!jpWl`dD*A`1-A91VD2zi*qNuBe*B?9?(6)blLmbwG0;hcykpe^*eurF3}vO2&1E6 zoc`#u{$*S@zW?<8wLAEOIQPqX=haqq7TBTbz3UcSNcG+Y`Q!uljO-VbEj;INyOMA< zz1f3CiPTT*Ofg#W3dH+2Fm1A1>C)B%{@bz_`!TWo{`Q$%-ut7`+cmQzHiH150iYzw zQg?9 z{});MKLl!4PA(SC|3y`2{SSBi|Dsj1v2*?ZIn@fvat`982RoNE>JK!Yn&fTK`Z%fN=Pj9TjP8iODfNLR%W@ zrw};gAjhX(W*`RfE8$y&4iCoe=P=u^ODimrLwVlolz2* zFw$wiLZ4)~*s^giddz<0q*sQo{n{112k{{oGjh2#GZnqVS-;s6VRm`_d+_Gk<^>EPC!I_2^GuMpdq=1-Z( zFEq^V8TBOsB-T{lKmP3Qn(`lduH;Tl_BM$#qfgMzFxZ=MFSQ9;f;ZB6D+;eRGeM>qp#Nw9e@aRar$io3TG zlLPh?3bJ{dI#^NV;_obb_l4|tLA2&w?MpxY%(_#&ES?MhE6Jttm})Ch|4}XKb4ue4 zH7W@-bH!$Dhv>!J%Y2#=*DzOHlC#yQFx4Ydt;zrden$r8Rff^zNq9ENQiWU|;HZCW z=ABo{G4Snwv~nfjRBd0JD7x{;n8t%B<30D>*^pUe9-Gv2an;Rmxz|)uqM1qqrO;p~ zDX)kWX%Ly3Q>piFK*^M7px5BP_d#9U@74GJ{lO-}4QH zC|(d*VHgB}M*g+p9g+g?HabK`?4ob?nwMqiQ1&O;e1K|7e`df zo?S`Mo58+u{>Uk(=KeXaE^dp7{hV;|$R6#@g}D=|Vy%xDhJD;qm>Rvop<+>pTvk?h zvFRKm-9wF@MJ~x_qC5A#RkJ<)B=i2HM9B>+q`J;Xohd!!G+ng0C@1I7gIaBe9+PJs zIRDP(kh$O1CAEE_?FuWR2SiLd7-S)96#{sGQK`RPtJQA1Jz|KYWnV^|dhFL0qtt^Q zyn!yt0ef1+{uAFeE>O>r(`{8)jee`{z3ErnRx~{?eaE=%Y6|r&+8-8eXh_#;sxPuE z`gj)O{oOVrzDc?@Kh3?M#HGsMq=6zMr?vK62Z^6a~K4ktu3wc${N|evM z=^44@i)GUBy}6l95lMml@q4SvYa0|x9BE(LHKs_!tvfo=zSrSM_-51Q*IH|M6F>EY zmcBPVbG~e1R-j}{wQS-aW~Gk#>q9%fJ{vHoNN;^2wQy%-kD+})X)}>jv;9J>)mw!J z?~RpDcSz0Xu1WRf-H)+o^*U4cAG4g>UuSjn$vAKK-#S*0m>~Makk63G*st}d=gQnU znJrfO$lCwNwReBz)o<#y zd0K6-HK5wct0v{z3QOG<(siTjaotb*Rc3xZRbNvmQAKR*Q}OJ~omgjHl|VjMp}n(0 zyZWqmhuX>LD7UPN+}V$t69(p0@|EM1%GFn1A%1SnQ;#aovU+0d|Mc8anH0|aWQCOY zF3zpG7P}7Jg*Tc%BaP1#4b4|Aj6H@%R;||Gi3PDO16V8Kv#KgTy?(qaTYTA@Ds~oQ z+~-v7Yhnuy)(0L#zE`)Y9>-*J^4k{7Iq#wv*&Za;YMLEXFO`}3?`` za^d=E$b!E^oGwO}7|3C6&sdwch-XOYZg%ULkecaDEwOgpk=~hkCiu_r4X@X@<$dgv z&~}Q52sOT0r}9CFToBx3RAk2}^j{0+;8Gf7Ns{V#uHf}Tq$6N;#w6EepSTt^l1bN) zqeQLa8cLt4?4({wER`v5+3jJmex?bjldNl8)AGRN^@SRDg>uD{hYmO^`@HbWqUESG-BC?((D; z7!)FbQ5y11RU($+jmk+ky)zTfCMQ+$W}YDI8NbB)QQI_^ggBjj4=08ldwBg~{foTT zPZ4wrV(*I+7o%dNI_%C|dy^flyEfl^R{aNi!$@28OYd^zC;Q8ccTG*-oorz26La2l z;4oXK|4Ys7_unQd&h^q3qCbdCGfICH=C9R9rk!&}oly!ZvB!F|Y0B$ReVQnBYW(&Dx*HK&;Ei(4d{?XT%CYs(N*+do01=A7ayA+kUK zo>8)3*l#!y)yc3xOjW9bxk*1sB{9hRc44YJv5$M$uRvPvyG`#0%lVqiUsp?>YV0Xi zQrW7~ujJPAmvL9c#&+Y9Et0RAz<+t2okr~_Idr_a^767*j36l5kc4IEpxwfdZ_)W6)# zin*FYI(BXP!tqXb{G3yayETu@Q2y||VaIb@O2YDS*lurD`s~-Y7CZJg1kAtTlAuPK zlcmyvO}>|w??X5`vC_bhpshH5e>=bSKl27wE?ep#L|4MjpOoJ)oSqAFkC>TcU8~|6 zbN>DAvvs@s2HFUSx+`t+thk<1wiEk%$47Qm|B17vj&~%csL3zPwEft% zo$TKBvcV@m^OzUWKuxZGPjjBB`dWj68x~8Srmi#iYkj9`xBP&%c55(8RR82#)6cH6 zdLHa*bZI|$beFbHz|2)oE=!iG?JJJ)a%87njFvceZH3OlFEL^#S`A{yu9Yn0?1BL} zuU{fwQJw3gEl853D+o#_6vl1v%kaDuJ`fOe$B;F7zW0#&9@VcfI&EAb?~1NEUmck=JHgw`@4krWp9F#sIfb4m zeq)1LIhXFR04ECL@ye3k-*lrhB5jhNUN&A}9kNMviOHTNatF3vv6DKxkJ~ig+RS*; zW#vsxkVMFt4>;6Rm!5<{iR&5%Dmhy4H=KW|J8vWG&AovQ#0&AB~QKYrQNA~ zclTHY{X4B8&o^0j*2Oiho-%DoTgrMBojIt>^-YWIXIQXq9jL8N!Ol;6S?myhDlpmo zbbN%-1exOF>^#YJxk0RU-5#lH9?MH^NSh?tosPWZQNG^X4c#xq6NDU5fG-O5LP7L8 zYGUbRaOSm$cDgD#Az@cvP4KpcLp}oV0&&mz!wBLEd&%`C?8A*CIAoMI2#F~2QU$FA;&LlC0OndBqCJW zzXY7&-~_u7h~_Ll8`l0sCGfQ+oA1lxgoJT<2x$=ZQ2^LcaCgAYzfiFm-xmTn!srLF z5?AmMcn@A=!{)458HP|PbOJ!O3>kW)gXDvwP5$PJehLZ<1cZj1@Q6!vU`64I0t=u3 zezpKRljCa^ydsbd-NG%59pnir0?Qf>nM)BGu+RU22>)4_N+ zCzHVhDIY3{)7q62kKoJa_sB|V?GU-5G&`u*^ zV1a{tDS$5ki~u~~I0`}qB?Bl6mN;Rs9}QOmEFpi~Xi$Mb2B276SyV6vSIV$(u*Cfe za-zTo6sib{4VMN>K)78XyFnsQ!4mP2P`^JbK)}H%FtB3y21*4ahY~?I7&g~%2~0*9 z76g}yLI<5dIDq0qKNuk@Y{Z=j%mhdMm<0-jq6GZi&;(2x6B#Dx2!lQV{a_^Fx(0p$ zeL{e6e}FzQOt`NM`2_R_L4tDO&tMMd8^eUJ9Q;!O4MYOm2cRE}A|?aa5Bdd&k%Bbr z6T>A`3a~o?FJWDPa9a<)hkbq!^b2TJST2M43f`&l)G~az1_2z0(ZvJv5cDiS#~qB5 zFo46mf-4N(_*XzsgLwg18P85bNdwp#@5Dhcm;pg2KLpqvufYt0!7D0g4Ft|TSX&9y z3Tp(xgn|w6$CC{glq=*fPq9a>=)fgpqoGh7WjbilP?s1Nbx@#1#h1)qw586g2O>x z!Bs3U<8Vm-na&5HlZkd5KbXXEjAB;^N^HRm2blv3^g|&6;y}D)zi?lG_B%I-6&4c6 z4GZM>AR&JKh!4v*0BBGPk(`dr?M;UOGuuq6a) zo@;40Yz~Uf(Grj;ipTaRfVe^Max964 zj^#@t(?^fRNF)&PMvlc$I@~ny_d^y8qrmMmq794`!j}Sf_o#k!Ch-4}ZFCgwp;2vQ z>X>!V=^!zV?#F}yW#m`}2CO)`jSQ0O=r#&-Y#V)S8OhJDM 0 and '#' not in l: + l = l.strip().split() + time.append(float(l[0])) + amplitude.append(float(l[1])) + plt.figure(figsize=(20/2.54, 12/2.54), facecolor='white', edgecolor='white') + ax = plt.subplot(111) + ax.spines["top"].set_visible(False) + # ax.spines["bottom"].set_visible(False) + ax.spines["right"].set_visible(False) + # ax.spines["left"].set_visible(False) + ax.get_xaxis().tick_bottom() + ax.get_yaxis().tick_left() + if filenames[i] == 'zap.dat': + ax.set_xlim([0, 1]) + else: + ax.set_xlim([0, .2]) + ax.set_title(filenames[i]) + ax.set_xlabel('time') + ax.set_ylabel('amplitude') + ax.plot(time, amplitude) + plt.savefig(filenames[i]+'.pdf') + plt.show() \ No newline at end of file diff --git a/zap.dat.pdf b/zap.dat.pdf new file mode 100644 index 0000000000000000000000000000000000000000..1661966e3bab9a0b9a7cde872d120a167110a13a GIT binary patch literal 13585 zcmb_@1zc3y6EGo&3rI*vhzpVu8|>1Z(kW6ROG?AiAd-TJbf+R+(jX-WlG3RnqNH?* zs33fILElrK-}}Gc_xqOL;oLiS&YU@O=FFTsXD*Anj2t(N2LWO!egiJ52f@Hlu!D&; zNJIn-(e`vS2ScP$E+|_EOE5$oWohmVh64ziU@;tsH0q5%+dB>B<2~#3_#M}-oeEgpbs=wN1@FD{1f&e zn&!?9u4q$WJQ)4~UgXgZu8!cpM}2Y+k=6o3WX#>HOwBdr0iJ+6Z8XZ>*%3p`)Dt5H zh>E$H6-vs%1Ly|RH9o+$XVF}LIXqOYynBin43D70WFoy?JZqyf?=5HI-U(@7qmIb4&;$m zlA=<|(CKt%Z|_u&WH`xXCM)L7ca;t}tL}Pxa%dok&*P&i z5Dsdtu496Eg={ybMdw#iIj>J|EpA+UDz-SluNLateZ7*q zcPZWX+rY@k#Rp>Fh{F!LBE@{YQkINN=X5+)RkzM7ExIHM^(d|EK}xLN)XnP=E>J5q zBKEI}8%7>ow()5&EUddxvfj`iw~)deU(}0U9Nx3eNn9x{Ej@aA=;HG=mC5qO$i}O$ z?gwA?bJiADNFLmHaQJEW==T2DiukR=SG!}8kiD2MMc2$BT(47(cx;RsCKyY%vS$@* z>k{K*`A0T}ypwUApN`m5N^$cH#bBJj=fCpTkG>SXX$*Yrta-pvkB*ZcN)b9t#!`Jq6q^BS!0Nut8z^% zH}8pk(hQZ;SiLvXu5_WaH>FtkixZ=ps~x=O)1863aUaUoA>G{lYUGq@Rf5?Df9I3e zHP=Hmu2~waHs;m0U7WwUr`VbZmkc8Wh1wV6nm(2oW^DDM+&?BqnZq;MQ!@@gpH?Z1o{?3ovY7yg1SGMjPAn{wU`x-q%wnc=#u^0b6v#ohdIL_lw|4eL--2_Hpa@|>84 z#TT-bUT^l+Rj+WW3@rxfJCSuxe1XhyO?eGei~jwgxI9LgmTxLT%w4x1)>62|1bomQ zM%%}`>e7Yn0~Si}+RTC7MO&$%OR4Nv;jNP!%<^n*pIcWrC=``GNxBpn$7VBMD^R9Z zWj!}?DdZWs(e%OVAy?YhT5k&v+x;gaEN3a)H6=IJp#k)>tt^?nT^F zW7+5;!eT;VQ$iG)G`DMXaD?M-Gzp{fE-=c}RVh9uTM&v-zq0VmsfVm(mxsbPjrHN6 z#iU}r&z@xcbEYuLOF3-DiIULNy}oC{Ia?_N!E9+diabz7Gv-eTk)CD|w|aW5!YS0W zgmE=rV}G6X!&-nn8(YtFW{}Hn^gnF!sOZAi`7A-d9^= zeQd|?T**N?D|1lXJoiEU2_6;2mDffRJm>2D;v{0^YXaH2bJ)S{kiqFsxeucQgPHuA zBXw^?v)E%>(SJqTcd0VPV-5Ft%Qd>OU&6`j#!+MxU$uKMu}$DdkLpmx8lp0?v2W z=Tn~ws`IJhNh|{%l&S$X3hD`upoESjr_O97AmqXxvE{kW z&8^;{#WPE&E2~l|__D-vFNk9G(%lES#`v1poN=M3sYws1y2mm!uVT+q6w5mWZjNCY z7>A5A=#B{$Z$UR=>CNc0$q%=_>foLh)M9$l0UOx$gPtu9_b{Ir@)v?T8{0*gu2)wO zCv&ucKjG||V5v?4J3`g86l{ch3K~Kr*6Lk|;-?*6(qiHPI;91xVm2U4tvRYpHN&| zZwXvsbdS-aj{>v#Oi*%(`|nz_M3Wh`T}mzz`H3r3$VWetI5Srsx#f2s=DlUaa`9o| zp+vC$rP;EHch4orLyJu>J|mUOV0DbTaHyKrtD9LKh3}NFCoh_%?x!aR~7!Rc=Jjn{38T2pt;IXw9fFC91jMg{Gm#Ytkj+4(Pm2B zi7;}Q>;->^LCI7L$59AQPbla9OI0qypwQ({%d)T5$O77)?Jo+rdC~bRD6@Vr=qZm4 z$&yvNBXg#Jd;e_j$XA(?i@tHCd5cy(Yqv5&VhUaIS5mn27K~cMDDL8iyo;i>qu
>ngb2%N49So2Qk}4kA6t13)fjFv5Gbbi-OY9#1Sm#jcuA2 znfO#B&>+taYEGS((Bl*6Wy;mxVEgzbM^IiNGqgzmjqRG`x)V-&3)ZcemXB9(y97^`o69bBd^%#CwIGT+1MxA-4(Bzc`3v;p%=1d90d#|mp2)(GjB$+&iO^eWQh{+4HbALl0yVffl0&klM|{3L3rdN(rWzZ`%lb5 zpm=-gV_(^9o#l ziCrg!7Y6TDu-nUQk3*V7v?dJ1uXSORlFWFi`EmY z&L~d5C*(8Zw3-s)Su5*<%p;X&+C})3t@k{aEEoIsx*{!uyhzWfy;vcO76mnw^2`uP z>Jhny05@1sH`U~3!o)PjshbXyEJ-dXekcc_V|g-PqUP!}G=;tK{F0)vn>)V9 zNH)O`>MYkTIL0qbbU}~J;4bYeZahR4AIrG~7W`oO#IafW_iuF!O4<^Gri8~{GP5c1 z*p%gqC2ucf6zjmWl%0`xgRn#5g(y3PU`O;?rtfGvXElnQRk6|&bjRIP?)aoIdTJU| zSZJy)%Q5-C<`7-I+#kCv^ypQa&#JJTA_KuE3d{0r0<{~LHR`2vy8Ah^Swy&A4ykvE z^nbh&lc?-#gugWRn0cv1Gi=ccM}i;>Y?V+w+}2nxk7qwW6l`>7#%d(ievC|JLYm)p z>|yYJV!H#Zu*ZmnrUaKjbjmWc{$+q=P<)47P~hzJmI#klmCI$?g3TAO)^u)b=57dH z>SnTH$|7gNH~SkFH|3U9dh==)$~xB@%z!#)6-fJVb1Et9ZN5;gF?`SS zv4(N!g&mz_cQ01v4|biLIRjbWv>CCBSCVYJt|&j(5D#-)Z$m4_2U)n=b#ZC-@htHy z@!GAY*j&~%k(L>BQm37MO?Y^I*ItJ}o6!c+H5`A3^z5xIZUPS1 zES{XJ*h<9Nnn^j|xYBx`=ujsWrG31bhDDT=-YK*reTF*?zJqn3twtXrbYz&(SB2O| zTp3RJmSlCavu=mZr-M+TcG1h9N_wS-XhKA&Iz$1dF7ABZ=vu^$g%TClFLH9#SL*bq zLTu>!r}N`z)c6Z1w(4>f*CNrQg?1_Q_HK5n)BA)A8tmYxwENIc}a&JMba2;q&1K{Cr% zb+t_}x?qoU<_QK+qnQ&B3(282?AVmKw^oZfNjW<*VR!H-4a@9v^Fi#x2iy&E*IU_# zhTg#9qOPyG5i@qFxxemwn&LojnEmcqgsVE+%9T>??$H^n@>1H`Py^=(ZCrS2N~rU< ztVv}ozSV3e_rM}DDlUdcLA-buxDiMkZd=%=v36spb6Ex?Q_(?9^@<)web)%!Pt;4Q zVw8mf5xaOQBg{zQ1W`FrC%>C0;vHvY*1ZCm;(d0=p>X-#d$FXCBje ze|U$UWJkM`Y(5}W!ysu-+7n6qK5g}`e9-jAGHK8Yw2TfxOIAqE8RL6g@$iU3maIWi z*y99G{0B>K!x7F(Av^9xk$%kcfM!rn*aor6$1Y&G$?&$X%ZU=^k#6VuKG9h#K`YgA zUC7v zng|8n8AVr3MXm*BkmOTbN~q;8x=VaRhR;`ov=z*(qr95hxgHWx-{Bzo*5;~%)k^ur z5_MmZ(&k9dc2cXnE9bkv*>062>oQUYGQ97hh%X_#@P#rwZMCxFI4Vs4er& z@!C?Bj`Pss$MMKHp3Q@f*%lZ9N-tI=9d~|SMnZI8Ci9+tTP~eeGRtDWpMK@N(KdHK zo-c}u(jI7wIn!`Lj@?&X!68l&SVTTHnO=emQzGL~fxx)4GnZSzZ3HkGX#JTOJ@pE1g6+(7)}zQ; zl)S0@$E#5gH*3DC2;B#&z6`=onC@?x>-b-w@3BLZeo!3S%%FHza=-(d&+FP7=8jz3 zv6 z*g1o-cyLH|O4P2ps$_~^$MQ^f6oJDIKR6<>*M)@Ss@BS+tQ%3QvTW8UbOaQcB>tsF zY*1sN6K`I-P}ZXB);gPk+tW`_lBx3qO0HpY=;D{*Jge!{2;KbBIj?9>Huud_my0^u zh|WG2k&+-5c(?ER^usa(iCGrqtC!-&zF!COAdfF*GRc~1cUg7TBG>&>CX*N@XL%_JCV~(M#a0YOu0X8eSB54GIBI%zQ4O1sb$6 z@PfWw@qf)~uDZP8u?f&)$lty5=rZk_!2LasvbQVG_G9+a;|`zfx$F0oGhZ@XX};xq zJ2hm3A%nHgMTCyHhh!{yH?ettX_CKmdgXKJsM>zc?qbBzW?Vsw+DE5cT^i8K#UeEg zGw#K;8%8#b&{uV_r{>}AJ{ zRJ*twqNEVI8LAl|!sAdWdBaxqgJ7rS^U}AU8*I5h&ZjotI$Bfn-T%5fkcfL^n5zCs zrZQswb+ORkK=xr03i~}&;KJ5 z3&X@k|4YlN6)S`+LN+lcc5oxUi zU4H_nv{C0sic4&{3N<8rz5b}@5U!w{9z5b*iyK(6Z-)3#;gx7ycB7}|LP*91=#u7Guv(BR2;ed$;$wAJQO^H?m!AXI{}y8v_*Gi6_Zx^?pae4BxyGk% zvsyf0kFWHl11r99RIB%XB18_?dW#AUIh8@g6|A*|pPXfkd?KeAp1zt@W`D!gf5A(u zXm{`QWi>Z9C(*taG#~%K*ME)u!V$j;iCaT323L~c#UiDbl=H|&`HSeiqa{!ec3`kH zScuz6y+e&Of%ENl*VU@`8_DONxj)OC?s27~5+tR(GJIZJCI5OOL(^$DdJCMLY$R78 z%?hr#DTv?8$wRs0)tV~}@|5?4jhclzoRiwA`#GML>&vd>HnrUQglmQa#9eh{}42C3-b~*UG?-i1zimMDO zS|+ZctJ^s$i9K*n%KUk?$s-LS9Q&ieF3GCe`{^bNGr>7qk=Nm#B}Ttu3s4Yo<=bv@F$kI+iCbZz%DuU`D>NqC0DK z6rG{iSSGHMeiNw9G%VKPIlW$H3u;i=Y5WHEcz$<$>g1S(8zz7Nb7Le4g9`|NfzuZ7bTn+QHhO`drX!R-P& z)H2r=QR`KsC!NINQA#E5O;?}#N@+n8|1cKJQQUtS3ojJ$%UDEXgk##EMAF{r0n}ah zKgZ7Pn6bt09j=MNyLJ2&Lb2E(uE`a-B48`lGe<-3-usr)ft-86F_6jm z{Eh6FisbP63ZeGGHfXxR$NaRgV3pXE#mh=GwS*aM)uD?CKF>+`5&5TG!kA$1#U%4Yis8h9%dhSQYc7~vEAk~0 zcHn}%;_8a)C4|f6sbmOvblnb5`REzH=T!G6- zv?6$iCyM62o$)|rvNmMT;elBNG5O(@?T<3@tVRo61jVnmYN%-L(tM`U+A0-YYYSQy zeG)=2KO>cdVopP@JW0))ZEwCl&pZCCf%vIO0`;jHXzXr$d-4mnXS%d(PO24jrV}!@ zrq@@ND8|W@*>M@u*Ptw2^pqSv6^E6yDxsxe1e48Lj{f~g%|00}71ry`B|&VQW!w@{ zJP^HMrQ4${lR0g3*Dmg{luL<4uBG&9w(n6@QVejH$`6p$a)VB3oD=;FXQN~q9T@!} zJToN8c;Fdh6YS!^^dDf6=wr7)E%aiR(bob$Wc@d%gO ztennYywL>|)>+*IlwjSum#4EiIT)DFiH2UOKsWL0+~ zd{t0#k`&&8uHTqTW&KxrW!>4bpsmX&sQ!=11l)*i?3}X@O9Xs^zYE zdI}ej-A$GTuZTM>gK3pHsE&z@iJimKvIgE+f0%{;x?6(-hweXbD3R(fh=EyXumn-% zQ{ecUW#}}y9@*Fp3Zux5q;$LD2s{hA`#F}B#DtVdYIvHez!{tynaLI9 zLoa8H9cfvO{HXYiKUjs+B3<4G1Kh?RJVxy4(1Sbv>1e zZk5P`C{^s4quXuk+;f>RJ3Se6U+B`-Hiawg-?Le}F%6+9;iHv^&hrH4WEn&}Wpf>8 zD#zW&5agp4gsKVHk5$F*pT#m}GSc%~eJmrWG#ruhz@f2ok3<4XtLhJ;`LCNt{@;w~ z{b$5}P$Jx;bf!y)P!%}o4Bo!46=9o^14sXL+lUbSWv+Pns(`IE6ahUQBNr@@5FILYkCFNc z*7M%^fHMnzdiWBwj|5IXW)0ya#_IS$%=*ucdH?Hni&qf-OV6U3O7?ibS{>kcB9XJH zg3c&-yFTS>klL$V2!i0&z~T~#3RuLKWA5@Q-`+)EpV8QEvbDBcU|i?n%l)j2zEL0i z50)x1I01{}X%TT~QvR!lAwlsmg%^gm?q_cr zJi8u>Qa3P|Jmq0ANc!q+bT}FB+7OMC9+{1}=u(yjyjp{XM8)#f))ZWGX9yrW+%?rB zZ8;#tX3V^xkhJJ{>3zneBq8jhQ)z;l?OqaxPV^?8a1S)V;RD|9xb7QHt}3%_ z4q4^VpDTIFVDc?c+*T>ou#+=HmLk8Ao#6I_QRZfr48v#ynJ9Zj`+OQ@S>8c$UpajS39Gk5Rb*UmZ+;NvsGC=Q z&4B5JQTL4fuLsyXSeZ0`n9%>cl|}x#ftHO{iopdF$n0eVT=F}Z6G=|FiauIHBq{T` zCFXw%QoNK)&b3J43H{c!-Dg&_R83!$9*%R5iK2Jz@<(0&v2>1^-fEfZ-Euhk&C-Ry z8REsqY1$P}G&+O|g(wAb7kj=It{*(%>Unx81deArE_jyDalQt1|6Ep|j(fFhZp+Hr!A+39f`8q#xvQr z1an=a_s>c;s6O^}S~+-oB-EC^_?qZig!`7TvQ5J*1lE18&QEroVsuB8x@3ju!bW$p zDSE_5deN}##gt_3z`@0pJrebM7I*(3$Den%|F&y{BYt&aDjE%TGBBc*X$r%4SSvQ3 z8HAeGKRHNu6GXLwGa(##EL5L5l!)G+HXwOA7kyRQmq8r1%F89%JvS=8P}i+Z@szSK zHA$V?7&$iJbfQ-er*{8nYsdiyahLsmB_5&00ee& zP;fn_?BnfP)8VfWss!75u$-Iez>+!jiaM&#IO~3Sd{jG{As-Zjr*Ay5M5hrJz!^7w z#^urY^ZMR*MT?)@_@yD+@AB(Ce27+*%343qq4BO!lkCo+;e0WBzT(1 z>Cy2*q9#7I5=RfrIX~}r5A1v+q+~a|{0I5{*IhUKw--5C^@eb1Kz^cl)C9)gxV$sm zuf$I_iz>=F-eJ@ikI^Sjy4#{eR37IrBQGZ@O4v+&XNLFr0QDWo{-A7$v|Sp~WRKw2 za`qqG9QzyFc<&41P`B$l=8hzPGJH9gcAAZ}-!5+JzEr(}Hl>sH11IIJ!is_L3ICjG z+scvXRmSijB?T7N#M@ezkE%EfE`N9}d#>!toPVqe;ystRbc{*w!#IsSXm!L#r;2EK z*TL>rX@cCzuh#R{ly8yfvL;zm5wJ5I*?HA2MuU-$ViKHgu$-;BQH^ zRz_6Tkg7(HJzJfddn1KoDfK1q#iZKJE2&auclCSroR%ppy^pZOR|I(fsVT5u#tP)r zsi0iYR>ui>JWwdcoB@A-T>&|GQYdHh+tzf%H#fItu#D+>#AOa>s3KWG3(@PQ$Y=4dMiGcXX6fq0pt9l#I=dmu~C#oYmn zgo7a#K&~P1-^vYuIRlI#ffnXgmN#9%g23A<3Jiln!4NYL#2mOnf+1Ew1rV|YL+pSI zKcKxMPytPV-*ys)2$5%XWh+}RXu<>=yo2E&fC zGJkpx$KW2T*^l_;FKgQtW$6s&1%`zofl1Olez%86p4?zw?|~##$15ii;|Dez5G@yT zI~_0-^b-?5bq~?~j*A2=D?|ZEU$ruow70Z12ZRDqL3tcY1z-{Z2mP-xp9tmO<@R^l z=CPbH?|24aK441!KNQgaD-B+l02l%2K7tPl7T^^C^YQY51d&iMKM>A@^CQ8$KoA&H zkN}^U=NNomKG;bO#83eIe=DFL&v^k{Ou^9jQIJ3w1W;)}06?6EPXHh!C;;Z?2kwDM zMgpDlV+yb^;DSgHP{4e^fgYf*(Kv^cTZEMi2g?kHNwKkCh3agt`BV3Wg~t6mNb=ZXV_XL8I2|100FGnX0EQzNnBbG=U|@!UD-hrt z0O=mXP+<7?1bZ{U5uCsgI|v|o z-xY9G0JBO~W|(cj@s|ghK#Kq0Z2=SVcgzp)%x12pKqf!C9m>Vg*1^Tr$^`6aW&t)q znc4vCbATYyKn}l)m4m$u=9>w286h|nCI}S({sRjRCddhevq7P3KkxuQ`!gjN50VGQ z1-|Lx;wS`x{M`kQ1KN_~C!Rde<`y8puJC~X3;q2E7+@s8s|ENk7)CWQKd}9G7z~c& z2SoN8jF%6X)8Ao&|7?o@7WD6U$UpS~c#>c7fIGh5+d_eG#qWJUp>SZ&^jljb68Sp} z25hZhTXXdB zr%PH^Ugm&~1G_(a3kSgGpI9$V2M5fG#n?1pkvhAe&@RWz0R|X*B#4DYR!t7{e*l+P BaPR;C literal 0 HcmV?d00001