From d9427c1a254a14f79a41f9c542bc1c2457dc0b8a Mon Sep 17 00:00:00 2001 From: himyss Date: Mon, 30 Oct 2023 12:59:06 +0300 Subject: [PATCH] data added --- add.tex | 47 + figures/angular.png | Bin 0 -> 42140 bytes figures/single_angular.png | Bin 0 -> 16089 bytes root/angular.C | 2019 ++++++++++++++++++++++++++++++++++++ root/mm.C | 4 +- root/single_angular.C | 616 +++++++++++ tetra.tex | 25 +- 7 files changed, 2708 insertions(+), 3 deletions(-) create mode 100644 add.tex create mode 100644 figures/angular.png create mode 100644 figures/single_angular.png create mode 100644 root/angular.C create mode 100644 root/single_angular.C diff --git a/add.tex b/add.tex new file mode 100644 index 0000000..71861e1 --- /dev/null +++ b/add.tex @@ -0,0 +1,47 @@ +\documentclass[12pt,a4paper]{article} + +\usepackage[utf8]{inputenc} % размер + +\usepackage[T2A]{fontenc} % указывает внутреннюю кодировку TeX + +\usepackage[ngerman,english]{babel} + +\usepackage{epsfig} +\usepackage{xcolor} + +\usepackage[labelformat=simple]{subcaption} +\usepackage{array,graphicx,caption} +\usepackage{xcolor,color} + +\begin{document} + +%------------------------------------------------------------------------------- +\begin{figure}[t] + \begin{center} + \includegraphics[width=0.7\textwidth]{figures/single_angular.png} + \end{center} + % + \caption{ + Angular correlation. Red circles - events satisfied the kinematical triangle + } + % + \label{fig:s_angular} +\end{figure} +%------------------------------------------------------------------------------- + +%------------------------------------------------------------------------------- +\begin{figure}[t] + \begin{center} + \includegraphics[width=1.\textwidth]{figures/angular.png} + \end{center} + % + \caption{ + Distributions for different angular cuts. Full statistics, $<40$ deg, $<30$ deg. + } + % + \label{fig:angular_diff} +\end{figure} +%------------------------------------------------------------------------------ + + +\end{document} \ No newline at end of file diff --git a/figures/angular.png b/figures/angular.png new file mode 100644 index 0000000000000000000000000000000000000000..a76e96b673619f958891e9b91e1edd7241bd8ab1 GIT binary patch literal 42140 zcmdqJXFyZgyEnSAAWBtHDFQN%ohrQqRCFu|s7R9<5fJGmv=E4jVWc{v0#XMRP)ZVd z4@d-5&?pfS2oMk zO#AN%@U44(LeN<9b@20Rr`G-dYv%v*vj4nvqP|J0P|aPX-p!W6CfW5fJ*BTNi-5PF z{kIaKZxN=L;KmlSKKwB|3U1A&wPmR-&adg6G}Nv(^Z3Jd^74nv7-X>%1rZcIw5AwC z|9cAChOA+ZGYaD81nztkveMhBd1}M~XTJ7x(M~`&88I5sX_`K&KC>d_Yc5Z-&~vQ+=VMXSqXc zzZ!PTT$Cw6khf~4%Od#Nxs#GTjV9=JcCC8nO4iD*F+NSxr;&*Yk!7XZyqwZX+nyc{ zBplS}$zt?*oy&~B*@VJbX6g)FcFG7MMG7v~Cq^iBt_glFlI~7kD9N-g^BWu@%kR`m z8S(DEa5vS_+`{X8{j<6T*mIiCHJq0dhCCf^m*`5nxe^&Wm#7l=%@sr7D5O@bIrofm zMxALlX}Gr4q2-^CXmMgV#13Yw(vI59)-sL0=-1OV%pyPMy|yYSejI&mWmJX2p|5i5 z=X*5U650f;AOnU@#L9G~I*K$MzH(_!8xdW05#gF>zP3y<=+X0D3yRktb{=r1#mBNA zgcr1Mny>;r)Z9Rw+0_AOuWw~#gsQrdNcM`#(eAdNJ#8sRiCLukOARL8IX-3`T98@N zA#uhc)26iqg0fCH5TjPwA|8D1-?v+_pbvw`GM*C0=k-a%@nAi{k2OtqLQAWL{oMclJ(i2tOlc!kDK(MX7Xo93Us zt|Pb9tm_lQ1hI(xBAi8*BGv}xF+@IOu;A2d<3elL`OE_b!dHh<^nr>Xn`|LTKX}R z??Bg{Ze|mhYWxZwSKWCy!!uf?qB~rfMaqeAiW6`sA0}+CVV(JKV-nY|d8Sp6jV0AZ zIMMhw9h{nrM~av_WokXarpomEptX^>pr&$Gz>$?6_tj92&dgC{W+7^3fR5_^%Kf>< zUR!0~qECaPFjhkLddsh$!=BS#@*Nb>yIS4#4l^ryQOUf6e)RZqLDLW9%1WVcH{~_Y zzMs#oe?D1Sx;hz@$+#V)Bm$Wx7j}(BJhOG5xa=QnsrXxSJ~P{=dAt(K3^v?58q*lA z_q0~Aqxs$6t9h%l*>26_6@uXhSj@YcT;FKE!swT^Ar2)@L*b~i`yefqi;{fj%Aru+ zJkj6Q&EbCkzQekNoQPSiDXg}%vm1{!bFWBTV%G_G_h)zf&$Eo<#WZ_5cceCe8-SZ$ zBY%I+jC=vt%P2gidyc6`SnJ8u91m>_$9C7aO&l-gcl7l1^dreUX%>?=5dXPW`3pa2 zx)f(vft|3(4RpwRT~J>qAS^3OCC-y=%u!JUSOs^*hVcnuI`c52u(hB*v$Zv#w!VPN z&CKwe4LPLYKshea?;J}&?LF5sWX7$boOOqVFTW3p(`isDX-Z?&b8j!{vD(O?qsI2n zyB7$a7T0j7MM7;`G2+Q;9syhCM(f$R#>;H`dEh6%Z4Jrz{(SL5y6*zJkTbw(?5hnL zxwf{pWV05&81`JStREq$&JUsh4XwD-ynN!kmtE_qx&L9Am~b;!;^)2Ts=3yb!Yg*I zTf#LpDLBm}ne3^R#;TUN_q*B1-385v%FoJfQ=!p$K@N>WllO4}nOh;WOb7BAoU-Vf zwcx}Otu7S(x-Mb|`vbM+{N!vDUk#<=SfVWv7VeaA{v-Fczvwj|^d_o0qc47fD1oK{ zTlO!+yGkapRxWHz;F=Ceer0_xk!)d|_!e$_56G-tJgfcYO(c{5M1vF*x0v$*Dh|qt zx@r=BjtQ&aFD+}ZSXkU;;{-|*Y_F69g1Rg8^>p#%dAcH!F*5X_IfG{}IE!MA+oQtY z>m;@I4$+1`)~UE7*>wd2Ore4n z*hB7GF?PAWs7Eecq>VqqHhZMLx|V<9H`i{}7U{ylND29C`CPW`?8a&SkXI@R!J*y? z+G8=f_Zlw!kseI>IeCQ@6i<^+;IIr_Cs{-GeW)5H+l;sOY+rP01qE61J~pWsCCYi8 z^z=vY@vZZZu0IGjYg<+BKpK8Fv+mM*As8Z#M+6V(6%e9}#_d;PbK@1RbTWp#&TYgR zrv4bhM%UiOVM8S!*9zj$N+E#~YK3L+RUI7ibEc+yM zabAZdU%XW33BUU01xx=m+vvKgsi^YNop7py#t#uf)1k3te>wXMm5t8~abUCVR%4er zn*uB=I~Xt7tf($T zWv!+Jl;V|j8XE2Nm|?#6#pRnT2W@fTeap6Yb(apJ>T_?H$Q%5`^16;AP?m=X=7fVl zC$|s?N;@_{c*Q$-(x<}&HBPG%sX-$@~KBDyyT*f?e<5wsGcV`ZfA<6|(! zHsPEr(I{=YXjcze&&3mbMy15o4%bQ>BJ0H&3+MPJ0VuIxT|y}+P+HX#{bM?OX2tki zOiiT2>=Qd#Ho3v)TW3ScbZcAb3~XPrSDF_pIy&C9Ny62_YJGvpJ24``u+v{2=t3(g zRk4!E8kv-jrH4m}#d6+C>FJxVdo2_8c26*C+}%b?kd}6>1oKh;F@m|eGd(t~o5E9) zN>7Xy64L%Wsn|EWR=jfanai8(MDX}1^QxtMm+8d<*6ji1w93PC;-xF^c;A>%L`S1G zm|(OWgPOdfvTvt$%B?MK3;FIPWQK)4y4T16K*oKCFy}Yz(*&D)6&^#v$tB#|C$nx@(-9M#;{*hcS6vpdNgikumk_}aL0IXLkt$nZ1gtRYOr++ zK;vAkLzFxS!v940eCiE8=te=jU>aYw)8jn%v6a1-yaR&7KBU4mcrNp_Eh5mH!9Ge} zbpl*9;Znxr6!K1y%lt;D{TQQ5OPaqMf?BSv3`b?acQBizGJ7V8^`p5O@yHs13~le8nt_ggv}6yi=gYw z#+tPW{Id_{T9?$**ZB?dI8V2Ry{b~PAnO*Y^#HOyKP>3KgsXexA-(=a4;AlbAIBCO zDlmC2)mbiN-Dk?+hAFPJ^pu@lP&@^pl=L(T2`Ikq#M?iJBcPg8(#D%C_KYNG86yoXEIr7~nG|Uo0 z!Lt`%0?=btK*w5$2x(Q9C}Xau=N|0|ua5>s^Y^A(Go z0|){493yjebX7(8E(nUw`I!$tb2))Ly~K^&n;$&)2%IRO2pjWV*J|gxDhY6JefN`t z#TLKP|_`@taFajlBQ$;#GPFbjq1#aOb&Q|9bFoggnou zU`Z1TKe$L&3V_7D+Q-0asHv5_GX2x$K!|)Ygkn(8&<)={L61oU%bqVEE>nP@3Q1U&XCFzTPFYt%PI14+^cX|Q|s670bOSjUZx{4AKoFf5MX3n)fL?@t(wwy zKSV6>h6HqCH0f#Z2XP&IkVuNHtU0E5nV7312FFe;tJ5G4_Fu| zHaolK2KOvQ`Cps*)hk0wO)63D?lW$%FsOicu(Lxod`NW{9U<;HDFTfR_6`9y^04Ab z+?9YlA9=|Gt8Qi*2^0fCZYFLJ%uoj!eBsMKKwC7rZ1RNKp8^Fa-qA-F8LOUYk}4rj z&onMI)vH;;fd=yLtq-$m@Xfv?Gz>$RwnJKlML}9k2B=q=F&jO$E)WJ$Q$OxtNCx{1 ztbTZkb*lK}-k$Kr9Uz>0v3h$~ul~jB5*=UIbBdz;r+YU0r;l1gc+`S;Au7B0$xo37 zMa56pue1=-H*VHx%|@Ai(DU{9LxBd*!Niy>Sfe03(4WB^uS@%%36)MSV1m6Cu0dMV zJr>C#NQ!Ri2?=y3&{1jKo`3;V6ZQeVL#H0eesZm0!f+bdDArIie(~@LEQY(01}kR_ z_&gFUW_!xV=4*Dedj@j0Rf7Y#D$7^85=8QnhOd5Sx5lC0s-Xf+E#(Fx)A5G>uXi^d z9v}4ePYl>OKU`8q_`Eax9PjZB6#!2cg{N1;Xukobo|F?+JD9BHpYHMq-4#JS&OfAe z-}#wpEG*y90258EhyEA}^QdA}B?1C>S1&a#7=SG;Exjb$NejL|ry>faGfj*>{v$ek zRhFGzSKJr}j`ak_dn*x~2Ew&X?;}_vIK-)}U0SNCEPQ|!IW`K4eSmSV-ymeVs14N| z2|v2-k6=Upm3(;nu~Zg>w}^=cta}0%#K^h{ockzKjLfM97mdPo9ci}dMp!Q|v@t>E zlwqEAkcH0fGux#SEyPsCzCzHCntnE4=liI8J6B?4i|(YQF9d1j0sGy#lLule9G=Vd zH@lEkui^4NRT4niv(|vq3@O{mb~;-`tng3iQ264iwt@9~gaS3=uMUC@KsxR0@_2{*c*=e`o>_|Y zm5mU!Z4hG-&ejAEX1VIaz@Ek-pWmhcAk^eey(Fg}G0Q~)Krl;x&Oog1S^LwfX?c4) zJFr_0^r3+4gIE8D!?ffulW($?OilO@R7U0tm9hO}E5I$9z%~N}Vt~9ic~Elm1bc%> z{(TghwBz{(a9ZI~nbhcVxoz!xorUCH2IB2#XLmU13+u97fB)LReECUg7DLs5D`K`k zE}}Mmi&!DjktNf?`IuObF&FCg-Fne-8Qk`p0Qry@mnfWDs&Y`ufpeSZzKjahJzo5b z&bZtq?WGxos$eG3DxfvibL1BOJzO6CC!j)ZD9+#HgxI(*#=LiOSgvj=Uz1ARa=XM{ zmJh4XZk`Q`cM6+}p)WRh>77}aC?^MHoA(-{nuPXm8AiZ4Y?rhF_l$D$O-jMw7gPdVHP6wsS?4KT%A2xtrws z*7rES>nrhhdppj@HONDJllud0-YP&N=#rW@Fd`6)Q({66-2I#*qxJ_!iz zyQ`D`til(o0^Zq9&Zz-Pez*4Os8q*HQ+{X!0FWD7+t4Yp^=ZEDdT9YbM4p@{0wy3X z$L5Fr4J98=I(2trd-fmtQi(DN?(JA3#wG|oX%$Dp!=#hHDaOpl_6f|CzU)G87HbDz zKySXNjA1JAw^_LA4gv;;21XM~beB5~*pIHzk#Q}>r-OCj>N>>fRZ^li#rjKbiIv{AC`$C zA$E4G_K=LH=43o_TsndGVMEhVOBf%8->C|;;Bb05-zCc3(HEx$o9B-J4FC`_#=h&L)K_KTz zPsKn$w`WoqZTm%g0z5@uxBPbquMOs@0q+nPLUt=)D>PzPs>T7;O~l{RCf#m9_{%Zz z+N!@Cq4slUD(Q=q#vTHrX)Anf(nl#Q(4g*ZKQJxjK5~&E@kb%5%)07=0(W3S3|{Sl z%1^JK5xQ9jk?MXiCE8{?KwbX#*Vz|(xC8tsw*rKgd_qW1X<$mWKKksn+*V)3NsG}pmzTyp^QL|_q0;{WzVR7+~i^M>k4%fT&_7&Do8$lPm&0J+b=rw z&&qXSIU@`!O>Jm*;kqETPYG$6`}YS{7=I`RYVk$s%nbYK;ZbFTE6!fOxSs;>?7{!y zB~$T&SKy=VX*YIZ04uB`k9e~Jp>GCuE@SfhX3yQ2lbby3H{jFzR=$-f3Yj_e?NQ*W z2%RZuQa?iJUDaps%5Rw%cks~49nIT9KL{UU6C(6~EI?Mk;xLyZmE%c;wPsHB>sm2l zY}ruC7r=dL{z8i0-mBDv3KW}aS0)Q6jX4)tYFYcLW z_StaIW(1G~nPi~WNxCdM9&DJo}$JtJ9X#MQUgNu z00j1->!=+!{SCTtnvFxeZ~));NX$PdaaVmI^-{)|GJJLvM#Fm^6%qzZfaCY^0Ci}s|n62$X0;rYiX*Rs}nGnHEzii$XNCs$V zfO?-E1|Ejr9su%hega_jBd`L>SA__4CE(mrmD=X(fFra!@Tavu%VY)94?z*^0iS%5 zF(g-nWtKq~WdPi#?U%3=YlPlp1M$&H&dkh|IVHF2&s4#VDVR&ohUA=a#|STWcdZfA zhmF{X-vO6{zFq(;aC|)-z6~f6nZQKXJID4ZrumF~7H>a@5=P(9TmRc8w3jrV6{6z; zEeVLLCfOFa*#m0-`mQRL?<9K$U@GQrpe28_gR~#sU{0_Z$`QukApacjc-@CzA~yyH zR*U*uYa4b~XPAb8a6g2>Gqn0gewp(%%>|xM#k#Yn{2RDf`%R!%j*DB`A+)PwQef#J zOy^~HJ%uqSDq+-pHg!xGw0ncrQUi9`BQz!+=KXcg)vWvDM?}0QEVICNnSvEwJkQ4Y z0@&=Yjw#uc_4vQe4o&)HW+8oHUNc<={T>I{VI?npGzz@*1O^0`v@#(7=oy639U*WfAF0+q7JdaV#+p|P3P#2hkEA7@>i90&qAiZ=I z5+H5>J9YI;WqPQP(b`3U3Ol8NzBt*(KjxbDgYWk41DC^*s!j!V2r>bjwa*607W5C) zo5tMNK0w>iKFc-z?fm7R4H4V*D|}b+5!Z?BUuxQj)a!)RrneT7>)0kWs(2zua$%)SU+t>LesL zN0II`J}72m4abFUBRUc2!8rid=i1&k`&w9tr2WNgnr_Ynfi(yUwvB)TEazZu5f#v4 zCUin%W}4gxm)hrwcoBcLiFevR9ThGFZL+REuV#2fMMp2s4+B?K{8a+Hy>%rp)QmgD zKvq(}hrcg5+%>Yud~z30@QR-@UZ1zUYMamvSa-^RT>^0p)fQ1t5Hj*RbGcK7uJywK-L84x}Fe9PIB$89xZOzcA9+bITIfa&iL?P~l3-6!M zuWI(8J^(?dnEQy3HI;uVe6wdaRt~0A)GtwKx+Cly=8s{AuBRWF!mK@X{4b#FfH#I{ z!xw5v?DMM_^T~J+lvF4!xq00YzId3Q4+nBHu^{)$lOUp-!7 zmF2skh2+oo0$u58j3UN-J0PmKIH(#`JSS4uIbmU~I0a7&Lw@iLK%%t0U;BWHBT%Ksmh7({ORTA!a-T&j_U?Ida z&Gn!Q?x9l^W|T5|=F!7+ahc)=J?qaJZAOEK=Rsk>{l<-@x~G%c*m|-XkB`f3GOb2# zY<1ai0^E)|4x%<_pHGyWTj{OkD;rE-Df^$k>|#EYJta74wlYw{>D?OWEp z5RypBVzvxL3iVy;VvZY+hiQ z08YL*aWS}a_j`GJ_3d@XeQFs35$Q_3LZ^wtp ziIsg5AC=S3@<#ngt;=XBZK$mDsEfuP;sJndM5pEJ9uP+1gk54TpZMY)PCR92rR3Wg zrKO{B@zK#jxYRBRNc;WA2D6V351EJIY-tjZ{{@9O07lD#LEw!lO8mV_e6H~idkn_x z#@f{C7vgA}BfOce&)hO*3q-dWaPx^RFyao+jeB4%84s}!e;}Jzlb_ZtJ>5^Y_Zyztq3*#bLRftSh(UFTNzT!2TZaW-ux^tYu=tOyyd;jUJSYtybrufEJp3lL0%v$zo^Cgqz|U%-9Di+lz=N2qq;h;b+FB_3ZavN-m;ZQea_YTv&FfY zM;lsncTFRzbmXqYm{ds)NVk9a6~sCuaW~JAuk7nbink*vpndnX6&TYrnNsL`rtgh( zfr)x!;I{`_S#tf|Y@&1GBMkUvtI5Zbm}i~IE#{#a{O|8i<{_A6nd<$NjCWNvu=4d! z+WtWm1`Nf8DSTUKKBgM z7VLQs)iZ=4J|RFpkl^y*Nf9!>Xfe6X1gn*3boY5H*F&$!LyK*pHe%eHjT|?)-=_1k zPu4LN{x`TB9_2~Iy;QJgODcsnmMD+{W>^iaZB(25%zZ+~)A5Ox1hp8Zmuv8z@UCZ7 ztrBHbtHk3#RR!}-9Xk7Vx)ReP2S2}&D8~Y{7cd0i>0q(aDY5^Gn!W1fzsqIdxncVz zqc&CpuultSZj6742~Xwte{#F?9JT}vU2)2Pl`C=#SkuO+&%=fN<=mzbcI|L`oY7pn z|GIa#^H>t^tn=*o9V!crjWFYM0?+6^+$J!fyFgNZVA!V5(~2IFIE{)Jk=k4iXomJr z%ZsY{Tx8>5H01JN9gr>Cpvd~UWC}%xF^c4?4_f$1!v?yBJf8?T)*b*xxn?DAtr~Y& zm8SZ3F|}RB?X3e;(QeW6o{6Eu>jPddK?4(;X_<+#oK$CKsczZ4eOb2_~sCGVZKKYtD z3ili6Y`uL!>J?epUp-jp6juSLq`V*U)Xb1D2=yC_f9~_)lJvg0OwZPlYMVf94`Ie3 z;;fu#$|2}(T#T%wbJ&~a%G>XD-#uipInZyd(i{M(b`9{gILqNR@u{U*)ehOvt)5+b zaiLBF>NG#W)Bsax;q!ZeJa_~GV+(s%<5cIhD_~x%Dx< zvk1e+*2j@Sbqa_;?TCgXq_VBnY0-iVDO8l-4~sXS zvI2hQdj>RVcGANX-UZ$=>%QqH-Q3(SxKTDKgkf`6bZvD108xNeiz^a>Aaz$&>cG1< z7~mU?m^BAc>)l2agJHHK4>AvpFzTZMC5)(ig59s!^-w4Ec_Y)Cf>^uw>QcDZO|1{>!MX!O}$0PYk&zfH37h|Nj>XZDBayU@T9^V!{S1=Erxxjf}9+!%LL zFhACd#bi}}4I&}J^Ud>*lit-J;`Z`HOlcvH|-hqdnI`dmm-6HO{zznoi zPfXV9n>HsM#iEz8-yQf!OQ^I^@y3ssJ}wlso9eH}Q4{A3>OOo{hn^hkIS!I?fHph3 z3fH|5^v?lsta#wyy!tselAiD=eQ;&rLXh&@xFNLkC+}sN`z_f&fhppO5%o)bvpELvl;>Onq0g6!V@ zA-&QnD{^~8vzth$4*OLc!B_L-2tK3(FG;hQ?=RAPv;Qwowq|P2gzu#TTTmADfq6(X z5v@u?wRjN60%2aWf-Uam420X|pN6!WJ(Ud7m!~q0B3}bZ1-4S1a0Be(K#n!nCygHh8!(t`pDf}6IOIq_uH>j ze)9CqORl4`t$V_@z0TKec;_`ziS?lcp>u>0R|!AK$PXP zpB08;yBz@~R1vsbiYla7*WcwOnLhDVdG6xiL*IAfyEFB@4mUhTt$qrKWHOVY%WK!B z;!=R-Kr6{$PO6AdQFT`jDJXZ-!WtIU>Y6kLG#&J2L%rC9^n_Sl&KmuXycP}Vi= z^GeFUXh}ch_wEt}Mu>{AWo?{$BD81Cze`Ls4wKOb4B=Lm36Z3MfB1tdE>+uOXD;#V z2YC@K$hv?`lJox8Re2j=Nq==GdIuGdfC>Pv;rCsNshd4W`x7zU2NI)(uS;)%{rIlM zi*Z5HYV*wMD`8fYQS009cdX|lLClB#vFEfy)`P0m@3lNkw~z9d5KUfm&h>^)=cyq~ zrUpNQGVv_boiJ?GDDB7wq~VWs;1eGBb$971JGar6Bu*)N1Ge1vu%>BC^H(E+?`elCIZIE8wz8qxV`TuUdKtM`vsCK19fNGva6Mgv+6iY z)-bJmjmV#a>wARg(BgU?R-^LgqexTLFtxzez;xjhqMwwYu1FB^A*La62x(HxdRd+> z*BpE2o#4}i5LD%N2w_Ofu1*e!!H%jd1bxp46i3lIGXg4GEF@o9Qs2dy_`pn$v(y8R zS=3c2GeZ5AxpDqNRrV%VE5&Q%*1o-pag(F-+0|zi?^C1#%sl-qqUjL<$aphL+#+oD zEb^IkJ;@kMK$pO?PG@{FI8CjG~_(sN%f}p)8JJy^l!*Y$H=vNG(PyWm$dXnNN@;Zq70^r|Vs`|4Y@VUajyj)}exm zL@ttiX|jFUy{^>;0;edd~>=Fc*0# zC4~K5I^eA3XWBE7R}a+7*|khgMiKYnSth*Kc=VQ%&1I~QdD=bv2&X5iZo*oaH11TQ z{E0S?g4Y*GVF1VMML70=pp!9fW-jbip?EBl# z)z7q*ehtAL{G$(6&aBYir0E)&{+UQ^5A#Eqp1IVV7=YaB2};lY!DebH?w7b5JsU^= z9=#a~R$U4FRoFEHBN2L z6wD4(P7k04i(7&eRajmo!qf6TATbfbi&cy}R`>C6UC|0PrKMBG85p#jKAQOf%il#6E-{Fav8Bm zmne2Ut62u(=V5Rhfi#Q%i)T@32M(c2g=gITFt}@CmAUHgr*YS^t0R4 zE8X|^64Q_H?(JnJPFl4p4ITtv*eBkB_u+t~mM-?YQ{$|-C2dB47`pBFlRucY6BFs& zwR`T#TZZYxf{>cyze^wwHgZ=t)Ol6#nJED^fo|%V^a{+QPeUA?L_PMCtZ&v zqx1ia2&2g-{s8W(n#P5t)fsMY&j-)C-7>SCqlr7KKit+x19kzz`w1P^;Z#-t_rClu zKlPPv%)s<4Gu!Jq{|YJ2%T=1%QR&@r{6;l4%qdOYs~X$*wn@}O%FO-yjn=82ruPpy z*_zCZBj4of&c?)7crtlIzU1oZIA?`FKCJL_UnW`(n{b2+ZU+%e)@@{bMUbWHTXm#3 zt09B9LVG7bY&rq%LMLP}0hy!jrq^u>M+W>=b>I350*YEtla4y&&VDFsr7LqG%|hTf zcBIjVYHu3F#vSOgd=FPH$X?xj>TjIBtBvB5ZX!9p>JLGQ8%5Pc#Wu4Ecbz2ow!5O_ zNJg`#Q7AP(13!)W@^(twfU}8I0*Z}d1sBViF}0zy%hU%%#@E~_wY5=Z9Ell0LV)Zg zDS(Zr=1r4)|NSLQGSg}A`x4-O3W{+sr~QJd*`ie)6GNC1Q!~h-b(QpwtX+pvHH?6_SIK&r7fvR)r(?`)=^c_DcZW}g>St52M63RR@QBvF=1*%SteU~ zfOC9NZfDV4VSb=`xAH*e5ZV2@vXL5b2k5Gh&&b}QRdF!EsSVQOEzv0E4{hQZ-D@zL zc#>BEas{cuA12eo1{wLQ`8@>QG`n#&?B82q)<_Z~KKT*qWCtFXy(`9RdP`k=j!Pm8 zmZR@vZw73YuIzR8ONU^M9udfIaG5bluml|pIv_Ie#)=;oX3{nmXDnPQ_L0|e1Rca4 zJ#vJ>$a=!1upqdPaMsQX>^=4%=& zjmo{rOhWAYn7rd{Qb_ygX<@yV*UJ0&v>wy6*-~UckiD% zcK;FSb7l;nl$!<+VyYl}TQ3;i%2_QcAX|}JD#G}PxD{}N>*ryFm0p~R(Ne+kzfT%# znq#NGPbz2t!doSe#hDd=k}nC>@@mpH;z4@~r(a7xv&U;>gLY2C-3SyiB^G4cK*p*+ z3?Qgm=8<94gTzHc`2)F!SYl$0TTmXnW z0sq15n8aP1B$I&4vdrBHc0cxga0j^-mEi%-YNuy}>0fhAxMYE`v;}i0PO|-~i?rDm zM&fdaS%dp8N`uy`YN*Qxo)&QjgM11$-Z(mPp5HX`ux*J(_zhApR%D6lSaJ@$XrwVDf-w>(B@EbhtlYcSkf%jFRwU~V z>ETbO?&Ye$#dTdbQcvL2e~sA%;#9)%K%ED2A!EMbBXqAMB%vWlSQ+LI`f4B{h@NYi zW>I;?69gPcY5+pbH#vY}y>YK(qEX6LB1kpl*zDc*TVBv**`V~pX*f68&cvK^-Hhg` zyjAl&Hwu5nO_)AtlmQ=`$7`!zGP#+=mR}sjQn;MC3+^VmtfhEYRYy9X=5@YXA6I%y z>r%3Tcp3;WgY7x`z!)eW-)!K8jTZAk_W|gP0{jXpa7!;`c>O9v_^{liXxV-C+T0I9 zxI7Ifgw45PLO(<=oL54hu0 ztz(cco)Q;uO$gQ=Yy?>pDDpo1%H9qfX(6g*TWO%=fv6z)r^vfGplH*g!I7}E+%wr{ zl#W$)V!1h1GVF7GVy!5#o==82vWcft3a{d~f%rGe5##2#}lY#dh zf)C&gK+P~TrlDPae&b6G?WpsKOSC&6sU7pAsL)oyzKB!;zhZJ@COcwir<$n4bX8=( zg@Rm{R;lKrThd`vjk6A#HHLI7T!v-=3tX?rZoew5#v9Y9tS2wW5|*0^+cK{JPEqzp zopfv7<2b}G733GTmkOGeeCU2(Ooh!n(am>io*p)9(H9NTF86Q4srA$&3u-C!3Op!b zg}CK6F^q<#+?$K17=;Cv% z#f#DqoO%z1(#8ACxp*)8k{V-jN5}a5s*4l1O<#kNHJujCYW)%})KlpGZ}OrKl1pPG zle$bnMxVZuSY*_JZ_4*>LMWZO+>)^z8jaG>=w6^I>334PCijMt{MeW6aFfikyH&BS z*dH?4pgF4CZMw&GI+k&m2qTv>o9azjYn=BJI*ZMEtV&$YQ|0*8WrmH-{)BX|m+STC zj$H@pJ00^Ki(fKe-MY~Z@@m#)j>{cW+I^(qLAw71K1A+H;Q+#6&R zd3~W>ZDi56OYtlWX^FkmYH|Ay05gEFP=-sptJGAa+%O8ax<7jww-p`-RhXbG2IPhF zBpZBP>(>G^^)e2iTyrBS-%b1EP{w(4aW=UQ!^-cf_O_zvJ18N!XFvBt1=zc1>yrsh zvqvd%u82>g?#f%GzEla8N%6{dADqeQWMoZw{i#H4uX8+xsw*wF$vj*D8=!jHIJhFF zuN8;4u~tmXB70gzda0zpYc7V!3%6|iYExX(%=eOlxEQ3XyAiqRKL1io?0HJSHR)-R zGl!i2$I5Di)9Na(EpctRU~Oys8lR3Vd6>~0M5yQ93vX6vHpd5nlyBX!d?=r!a<{B1 zR+@vkQ&RYF=q2**84F3(N`uky`4!~eVr?nGx7Xs(m}>p0Bw;HfUg2(G&E;}t<2PW= z-!IN{Vl@KhP4-HSgYL8-T(V?}g(PlzX`+Re`VmDwuP92csxP7Gk+3p>FI-S$|3p95 zAAN+F6jCiy8x`NA(YO5E3v^43Pcoz99F$Te6@eNp?OQ8CeD42KF}{xS>5?aNlKE|q zuhEyPJnIr*%dmBKBSuPSA+c7muGb=20fZp;X}?_S`l4{Bp`HP~ba0jd`}4Zn1hlq$ z@joY2##PpSbSH+UZavjE*;aqDxiuA`H=HQbA?zS+94Yc@C@E)AxB(6?v8aGT5j4@`#JIJ8pBsv zO&599O9B>3BS9vPPoulFYQoHlNk!xoQb~7q9cXP_lX54-E?;tB1u#ZE#rU9*4D{xB zzx&S7UV2qPz%n~3OGguRQ{XmqKEHT7JX5MUlTYNTn4zu zyT8+yt&AbaO;QhCuZT-jI8N?*1?roo1_M)4QnSo|&J$^2u3j^KHI|Cv0eOQtJM~KrAoj zt#8)#_kRq@&DPde7kd2D+RD`nq?gu5Ohh4D>lmME)1T9<+@Fwng6}<=<7$Cc!yH15a_Iwc-znDW$DXqo|G*PPqp|y%DG+I39H|%Pe((!^7kQ{I}Yja1L zA6L>>mrBgE%TK5X7@Up{TPO*uW zVs;mGTKavxUpCN1Hg$xZpQAdWnQIU?Fr(C)cTe*kYYLDz zkt#!(P0tIjtPGp4Et)5;3FvDq7c4h|TDwa2wW+MAqm%9Cjy$KaJ}%L^*R*>+jh4(b zt;*b~Eiu&isvXH$UF%p_b&@vP|vLCzW??pz@Wcs z0SK{F>=3!tPG#2ildk2kp*J6d@Jlzmv6jeMNei;$D7zGiSW9#du~W&4i$kDsETWI` zQtl(#o(<0Dyr zG!p)(sQh1A6(umret)LkMor!_=;g`to6gn}Ul7p6OjP_P=8c zT8jTE z-_`+j`AI+9dIaER0S|Yf)9vz^g@x^j=lAXxJwJ_X?kb;?pJwFJiOMh>q`yNsq?{f>|9)cSUoZb z{DB1B#2DK=Yn6>J3g?evyh03#g_B8(sOUy1u?bz!=GWG}#BJL3_qyv9$&Ub76R(4M z^|9jP;RGJur#JY{LcX)6Ch@kTL>^|!HAAIDNB8~WvOpQWpzO{}MdTuyoMhJnSs4P$~6>>dmC(^NV{vh`s?xQ`60FG0ZIoT{IH}O)a?6^s#fu-KPopGXT+iZs|E_&0c`;?pT$#=`+1=5A zd|WiSb-c3u2upo8bL?_Rv?uKKKL3lqdToY_Zf0v61ip0_#;KtcS)|Xsd4(O4BWKmP zUoN76mXrBgEq%g|t!O9z*pgGar7oWJ?Fej{B#F->} zV{{hQdw4)R=;a|;&eOv#j{uMk4wMLMWZa>fi9shmfB<>!0bx5+m6`#4ryA(*TwlO? zCkHgv|EoX8k%9mX0xAT~)N^441c1p0-Vf;b2gOVHR2-*?dc@ z`GRut@jGgr)f+dS{ChhM-F9OOF7vwK?#U?kSo5VvtFk$tUpkq};>2ZEe*f!0e)iW2 z_gzqXZd`s0$SB{d3jhrg>q{@_Z!N9UcqMFw1pb!(4MiQb)ScXgg}?lEKlGzLoQF|ydaxE#9^Jt9t=bJc0%Lt|&Or$BFy@xwG>n*}&9(1zt=F3H##0wz6D zbL9;nl+4V`<|1%7CJO!d4Z?IcMx?y&>&Wl=R$(j!@COxSJS7mW?%gs%XS475Hgz}< z$a0%HZBXCeBd1@+;a(_Vd?^cM0Kq<%9kYlXMbc*qOd;_krm#HmoYQnp4lal)4a{^X54fn zs1A%|}1MUZ$RDoNRS7?{^70)E9^6qHT#Gy+^8C?&AZG{J} zOQ=#73>h^d+m7=ZLo&39O(#FQxb6fSg9>L*1cWsMuud@E=`Ux&wY!9v*SutIYgkF_ zdR=>{?OSqPz>CtOujN239|(kXmdrlaaLi&(&2dnm zuZsiGBU5DuR3Dp+TY3q4Be|fH_bZ&f!Ex(jIFw9_^`e5UNc;3ynCv$zD9PV-T{l!wcvDzV9?3S2#{+90h0r@;KO43i^AO-z$R}1U;S_ zI75$J2`-=vy3K#&_bxf1&|sN4baCOsU}Ra=X_W`iKI&SDji1%Zi2L<6w9yT5;wdo@ zHp_v^v3_8c>!bHzPZr=bz>}adl9&ly{P%AMwYyyRTN;(=@Nq~ZRlR|r!%iumj5Kd} zfE~IoX92kWgHZtbPA@7jw}UaKR6!Z^nKDTT zS4P%>^?Ok4UlqswZ|36l9Pa&k`ci&D*u|RlkSC6AV6_wYqiIgyR=c;UiCe7(DAe*W zb~;Kt^&{jhtqo&_|Jb#oL{)VjXi0r&2NdA{i+Y{u`< z#yns~<${jg%}@BlhL&EAGYXe?{TtM=e0SX2aoekv$N$|8IrRuEM?EG&Wo#W z#kE|&4xlC0_n|2Xl!G_^6+w5@%bFBz-`m$P=oxxdoxdj9eq-vJD{)QTLu9b&{Kb=w z|A)LckB9RA9=|UYNt9NksCP-q+Cqz^@-9V|NR${%l5Ay-u~nj?v{_1#cb2R(l##KF zD9T!dVT`F{9W%%>%$S+`T;u)u{J!_&zJHI$@1Ogxd;TfcTyrh2<-E@8obx;ftcy+r zHPZT9(|5HT*i9@pQWc!OGV*~NHQ>)uYh1c)eIAx}&{ z#Iv=lqC6kPtiJkoQ%56z*{NTry~W9!%wNbGq8M|$bku#w&j2|$dScz*@|7OVq4jS@)@lbHlBPij~Y z67o%otRU(LrxF%-?X2P45>s3eZ5<~k7UJ`E6QX@$8r98WIS1`Ppz#6;SXU?FirO7i}qq ztQqgnIphgCm+V|9J(ki5M*(@%o!y9z?7V~#RL&(`0MVadiF&EFfA`nhuYZmloTp)= zqoY5)6ad>UlHuMNv>`6qoN}>k+uLcRGk$R63Z1^&un4Cx%Xsmyyi;|CgI+FN%~~to zX<>bs%%*{k0Aw@(qr$cctkHa>TkC>2UdtiF1szX+igdzaZ$PLP2qN2cVox*VW3$F? zJC`4QC^1#+d5B{F>qp@)2V{VP3fs5!t)>8OquRrF|Io=2CGl!UgW6h{guTm;ofMRD z=!+vYd-4Ose-arrMl?m?Z=UN)GR``_^OSEM^o4|Kh0>Bp3b2k2&ZFq7Qz9v5l%;1R zDxyfL7jS`A?%s;yFO1S|^IqUso(9BCn~Qs zR_mE6GvFoYNL(57IqcMB!9#zJWHOiJ`$wzB&q!R=%|ln_`p1C9vE3&Tu|~tsqT?^} zbh7$1<}#%&K4K(2JaieWw|!*|;uy0eyQ=*)_J({HJP|?VNt_(0POF?_F$#8?-&y&U zGL*N4=2XsFb>wxd02W2$)AM(`Aqviv?$PT*y`>7FgVuz>xTz#X%iCLtoCK z+~X{qJ~(;cjCH)bHpIz46aMDMja0xJL5t%DSpDN^cLaoU+r5dp*LY8FR8gy9nVy$R z0HIH3-vYfQ^94lY;a+i;?>G_ScZzGupum`Z0=X|B5*voTit;j~NvS)RK2wZhflCC6m_ZAu5_1JfKWL$x^3GLnGu>Uk4K5Lc+~)}L7u0p0f38;=_{8HZIPrD2%9$+UHjZ+_aX|q`O%f#xT zdjC^jCp)}G4@xdSLS}i;ZQr=##2lY4w1-iEV$_jkkmaE#gR^{dGHhnRr!QR2VP!DU zfxzb>G!FhLNWv+WLY@NlHlCX`Y<_lCcSqSgt`s%hY48x6^2Hp?&kFb4Hf z{)Qn3R|L%qMzz&xpr)Xc4#Gw;B5>g;Z1up3P_+x14ZiWHH#gC7ydqS=y~7@8IVhU8 zsdi{1j<22PCx+aqig6X!XF@(R4-A2)rxVShCvH|!tL32Zu;46H%DRnEt%_-BZ zO+cers_NISi`r+p27HEzboe(X#}b+@mVu&(_}c%Vh(2&REbKG_Qm<=BGXk>qw?hjQ zy4K(;#2k@R&BfagdaFR1{(WVq{^7JR2a-sse+}N~hAjg!(zxLKe7}`eY3jI$Z0{qj ziSV8zFZto}$tKFbPTYhBjKiUSI1H|;JUonJ%{@Fi<^G4Nj=b_)=z^g;OtaQBP|&rG zACW^m+AEMo0ed*sy7TnLg`L8WlXePaI|HdDzayxd=4)A4p`N>(-!?b=!&US;12^ls%gcr^SqIU)AyXm}b=)1yKk42ifpINUU_! zhgj#qS}%^n4zaIW5KsaUN9StIq>~#CC9Aplfb zwXkV~5kf=5ca30Zql2yY`(9Is^XOIzwzdXM3bB=QvNb9X+)Ail*`>6*vY1blP)-I&dK6e*K1jYB7II z$sHzE}11Ot|T5Pf=$R^!VAM>LBA$ShCVYD8oUKfj@e;APwGm2OIx1 z8SO(Y5f+=_p>$0UZG#4X{hO$T=#>RLkt}vzxVfXC45ME*-#IQNTjDJb9e(m)K60d= z3I8xV>p>%f4+`DC5fUryHQf{bI--g3 zfcdae=APO|AEB|(VZ5zXD#`MHY;E=}f%e94kpDg3^`WK2naJL=H&r%Oe$oWRU_nWJ z^F}jI7|P)l76=)etgW|1wm8Br^eSktQq?57Ytpy4s$HRiQ{vi*#qZEIM3T;&cA{*J zt0wrZxn94N;F&R2_tv^1zjVZXN4oHYd#8JnoY;F!ee+k#JI^b=!h2GGEiUv(vG<~3 zy{Rk6>YWoodH1zO9>Uf(-?J-6|;EfF#bg8ck(}snlnlrs96afGJCcj zlIx&hPG~feG@y|UwK+A~R;JyaI8LQ!R(3EhlUw7OQWCk3xD46R`*skMEW`(;b^t!e z;0z4HevwMdsos58i*y1jqh0oLa7b*tSL%hDT2`UDoXM+XBg8hpDa8gjvZ!@3QTJCx zn&xudYNsk<4pLIGa(gS~BbxNqQylJ}5kbdyRSFZZrGr>n@SD}Q;^(|%bggV6@+W33 zEKel9mVg-tuLb_J9qqX^ME?-q#?3esqp6&>3yu>&COZ5Xb!@3R{&Y9Lv-g1Z+oz97 z10nmqu0i0^RhPtT4yH&2S5nvgNGMygQ63|1n~iIC%J5tHnPDBFdD~KW?k*JsG^2=@ zFwcIkcthG^$qCC%leG}9+w9olAT?#Y41sy<=o)sOB!bgAYg+tiXKo>4PFT8YRUZ>H zyKvndbrjCKlrK)V1~m|0%CBkA>PBzWYxe#%dox~ZMrMm`XJb}m!2yh?po6aQ(B9hG zLu{NIJR%m zE@jDNp~U~;&wFS%-AVN8wBdx+Pd~+R#y{Kczj$a#z#Cz{Zt{qGi$w6Lg2b4yFIpJ+nwwLXEtkUfeJDiJWUy~m7Ic23Xk8{2YIyM zq-~PPu*A$BRWc9q5rTwy#ya!pHtHv=IfNCYlNZ;-mL5ON@(=wy)6f4=(_sK0Ax4L> z^oQN^%yT3VgjGnq)20?N^_}GdZT8B;T>hUpt@Fxx)Yn*(K(7DR=a&8@%&7`2JJ{9$ zhRtCBf;TOS_m^9)Y&AQ5%0etg<=QA-wvHM$=8gMKtf5w&Y1qM@%F-qp2^TFCIQQ*LXl>1S(R zai;tjd*^Xh#U;}*C$Ah2FVErYm396l84H)&?j)Q0eN60;VnRWfT`=H(g`-)w>?1U|V&|mG2wJwvnBPFVZLpef+ zT(YD3%xnL*XiB7rn*ccGRV=}n*kIK31Re?%_?I4-mG|!}KTT+Y5Wj`nF+mL&c+Ay!Y&t}>Nm;hZe+rcF?eLout&fLR0kYkQmK8*zDbY zR>sLaM4qbXlQ{5O^0AM1-HgzB2x$^-N!bR+zRMUi7-_#(9N5kE7}D%WC_YOjnEHxA z^6dYC5h6~q9v@pjw6#u>Cyd2syvnTc1nryKh>cq7;7+y=H6?Aypekm2Dx9Um>&eTHY- z_tczno;M-7l`s{HOmQph9JCi^Uaj24UK^j38rD?HoB_g!6+pxga`Ebl=vpBXG;%s5n`g!uTjc&Dn6lx%)?QGQ9PzkKeB z8Oc)}kF-u6N%o7JUF`9sUMO3~q^| z$EO!RkGtv2(z6>&ZA3mYj4ycMN1rGy$%AvMe!)4DOB^-Lj(rQY8~WpQpGmFuY49}V z#_{Cfdu>P<5St9j$!E7Qzz1j6){YKWsJrAgl)TF^hwmxRU=AGi(yC>=&dJG0Ud)i` zi!;Jdd-WZe1KnmN!MDPVTU2iqPXD~_gE`sY8RJ=L1MFZFYu`u(EcKc%f9A zUaATPwWlR>)r#)2VD)!UB#_Vke3 zv^_I2xIa7IaLf{2aT^RS4~6RIBz0X0^s#dv>&V-8b)%)vXLQqZ$TCC+&_Q?hojB6- z>%_8qHVRf>kBQMcn4Iy`>}ZRBT!|Al4~l7_Dxnxln*giH@Yu?*k}`Lic}F*Bb;u8; zLfUCXDXrRX>MAUrr;6bgN_Y3?d)fS(PIQ{is9maq)ccut9y=P$*mL_* zCVnM`+mhsRpW?bag1g4Dbwpozl*YN*={tA@b$JApdiV4U)*iEoYWJ_`ny9YL`9Qp& zMC>Pe-yG4n-gxcu!3RpC3*?|Sen<<=kN~Ri{*9=uFRRUMcstWlXNgI}(j7c-EDCD` z(7or{9_ALcndi#XDO(IoUY;iyMY<8pX;?q@enR`T33wAn*Ly3#04`|{3jYInqG-k$ z3mQ(*ZcO2#Rvb~P!OyP8f|hXfK~t4Z?53q5Z%{S*Z2pkf>L+TVPFtRy-e+z6szjnC zAk?-h;O#E9*PF+kluBpCVA|V;#v;*IS78e(iv{}xsUXn{9D84_LR@4$9*k&w9&6{F zI5C$}tcX7|IWAEr;`sN5tb?I%PH%?!!q_Tlmre9snwv|s*q6C;R7&YSuYceVcIFXt zj%AC*y&LF@P&W>W?nAB`0f+Z!FnCx{7S;)ivLNy%k#N*7Bjq zsXx+b7xDHF_s_Nujg0QLW4zQ2TBNLvH9Pzn-Ex!5BJO5*%_V1|7E|!rvF8rpCxwDW zW(x$?$%>vuGJKo%Hy1fAj&xFW&OJged3uJ|J~*)3?hh+K<>UMZ{T7G&o0fX82{`c3 z)e!RUqxc{WHWAl}&)pJ_G{+Eon))sp=kc8ndcA4(zBh9028DPo@YIKE8;0BVwDFV(?t3D zMM%Utz`XTwSGE7T^8FsaJ-_V9m~hMa!YMZbV}y?fX^Nksr0D zk(FOEM+eCyZYqcpU<<9$^hcq4PX` z_xp*@4^?(P8&FU-d9%Fj*%iUQOX$wWpI|EL_dGG8zU(m~Eiwtr0o(`cEF!eE}*dF<_5ge=KYl-Mg^`q5MJ z_u2;^5 z=!l5+G1-nSzq<)Iow)o?f2!fvx0ie)55;+?Eor1>=0}ly@pYvow0zGj$+D6HpS*^B z8?oqAaFUhJ;J9;rtZhA}Yx^GmG&-;CH4tpjT6?m;$lV~_?bD+}B zX<{hIu16sO-k^w4vd{W(|AiHJGn_l*-8QbNOT)-z8}pc|hubl($kTB}!|(>+sn zD|ErCDT}R{XYCN(y5X0O{yA} z-B$1Ngi%`ZN0Vr6`o!q{KZVThnMz0#F{w6cWQToLkm`>bBg@i~{~pllIF?jMo^f@j zkwZIk?XNbKHsM|LRP~yvXBwQ|aBGTH`?RbAq)MlP8$G@23iJp<-ccvg>32#0y_5s8 ziNCd|#tQDFRr9BG{7Pq@jHyTS-OJKi12bO&2I!NQKJ~mc&V|P=Nm#%D%OGjfQ8|DyM}HRd~7Shc;bs-6H_AYm~i25 zXl=UQvaz8p=mTM$R|9f&D_d|?u-xVd%%uL9$L-fvU8O%=Cudh2&z^gk%(Hucx^h7N zgD*?@Kw`~?#EG@kj6RCHf|^W&=YfPnns1vh=3W)NQd@O|U|K4O#`AiAJb*omcKGR{ z^z1T`Ta`+SZdde-Y7AF|7&X_E)7&&lH+cf#7tb60Rych6RV90b)I%)Vgbs^y%WAIE zE9A4WI-CSzRI;=Cc?G+-QzTp+=w?L6FTCGfeILtiST6f%)I9z6Q~#chN~0Sbhed3q z-`x-LP&*i98<2EC+6Y6ycaYUIvL$3D=i0sb?AY!9yLOn7KrtWKZy7uHC6T%1}-Alq-NzF>NtM+T$C$&GB9opoA^X+BZm03Hb2;J{%T0q zia~UrN(idftPJ+M)qFp7tj%9RcjhJUV)M=iu32{HXSJrUd1y2Yc9AC)8c^g@cC%z( z+uC-F`@4Z;XFG=i+?IFqm1Aek$|t*Ra`z7$GHCl*`giIhPe`PV+}kWk z&8;{SURqio7+<8(=$YTgkQ8;s>6#7odMdcf?s+oSwE*`?Y?GuuvCJh$_@DC$dFkqP zUPEq)+NQ9)ed(tX9P|&U9=J5%5dVIY$g#WTQKR{NAhS*#m+d&+bL@*l4q4%vmE=qx zZRdjnacfLA;SwZOFol^!g1xWbPd7W!-gnO!C&AkcD%I;dF{|4#H8bb z0-ZZ_G;4)?p*9oLKMRs;Zd0Y!Eh?xj{Og+cANfBIwciuRs@z+IIoBDO_tSQdll;AN zof+kEvG9-XMf=4MZwWM9iW~q>O8642mgwKG_6&FjNANCrp?NR$<4y5)RysZ%bapw_ z-Z@)UXq^c-B}){o<)80;e)O>saeeo-+{o#YgbDwSR#+xsSw3A&D7NH%Iz(LCo!_(} z=aTXHaCf({^WEQyzjIz3+}QIg0QOqI!^uGw(B6$)_Xm7Slcc(mOmfu2-#WPtX1ku% zQQPQUraXH*BO^A}NQdDw!q*52)ctaUeY$#PZ}jS!yOt1P-oJ!v1s()4&)nnrr4I@7 zML){PS&@P#2l_q0WyDM9WF{mwB~VkS9f?gHS*dOR{*^#gs3w*qmoBSGnDUq{ew#d% z`!>EGL;2tqm`5(6<#MQxAI*6E-mwfhmp-2BNrh+!@UX=YL8{Ff>)iF-JtK4AzW$bxf$P+a$TBWrwX34icDl{RcN@dz^Q%o2%#SYO9cjMEf5y%2y(`^|by%D_&}&s{NrzRV1nJ_9t1$rs+UD z0(~n}O^|FoW)b^MdCBs`4&T^`;8mxBgz0hV^XmD&j0#nrzkTO)Qy4!V8?0g-E?OSP z5)89AL`#Lb)hXDu)b%lVO^MT zUm@0)5#L4hgjgYrS^5qX>k|eCZMw7E+S(dCY_osRkS2S;`sAYi#rhawtjJ(4NsvV) zRull%7gW1|_30O^Y#0b8Y%Ji`_aDUT77W+-ZzzDOL71|2;EE9Db{g7*m6GtWWR{AxXDYqK-ktq0v ztKRlhyhzeL#)r_WVYUi`x2XNk3oaD=tVaX16yk0v{`K5Fm3M%<@s$PbrC;7?${&9N zrb!Sm-K{D<5Q}j3do%gAa~lG ze*HdDHkpO$8Hp5E>nPJay5hV=wz>igf^!LRxI-lnDlZL7H+9-eyRYWj}|m=0c;Z70{sek|Hw zs%q`Lfl8`B0Jp?b9qy%hR8M)kaNhzO7{9{Lr`?y}e`Q*GD;w5gD`n9_5J(4#EL9nr zwn?rC=8n3T@hUYkgyB(SA=AfXCy;!vZ%4DB9Zdbc?(e`#*fMTC$;+*iB}@MGGR$D_ z7>i!kCLc_r_`;IHkv^zG5PQa<&0&hn$WB2uK1rK{Sv7!<%0Pve$rfP$4_rqM~jf}R} zRMr%Ug&Q>hPYhy=I{1FL@xBgT*X<1rBXd+LGIMgWcGhxQm}GTpvy2@46gy4gU#1Vj zZP0K*G`=F-tj+j_Rj?T+l+E)?eUhS+2_aN-ksOQ@r? zxb^bL@NS1sq7F0Ahm9msEq8RRyiT&b+$I%E3OB)&{cNWA9i7d4JQ9}r!YnJfDxbd2 ziT?AdERPC-Wc*6`_Zf6{Sxk8+i{%?ddNx|AT@jro@pdf>hM7C=L0#WvYun&Bw=TuF zI`=mT33fUZhCxXg$R6a?60|`b9#^M-zfDi-o4CH_jNXXxc)RvCyqJ!^eO%Pl>`)pV zqn`L7NYeV^bS}ZN|Lk+u;O4I(R>_h&a);cXn<{Q%UeGO@I6lS+?D|B*q>pl<-tm%p14Wkw_9aJgPd-&a11%B1BlG5KK4 zu2`uS?d6zM{ie_&hu-8w-ehQf*{#x$wt1#3_(DT4*pMT#-=N*1l-)fXqHhq=Tah}r z8ar?KL3{o0>D($(yA7|iXL%l{g%P$1XZE}f#=hH_#l!X3`9D2tKV}e$1SH7nk4>HK zq`wwYeMD6vn?r;!$CtxPFVIr0_IKnrj|47IeXO;3u&YGmqd|UEC9$H1@k+P; zZpo%S#r=JKSi)ak_cnWY^jhf#-s^<-1-CIg9r}4?D|BWF|Nr^Fh4czh4=?&z`(Kvt z{}+eQvF*MJdjYi#6E=8K)oH@mhe*x58d0zS&%Q!sh|4$Azl4p0d_q(Hq8|&@|4Ntt zmqVN&2xruMl({UeAQk6j9jdQWk?FR3-<6};FB455ol+d#x_#ESjT>@_)S`<{X@r8t z5@oz$w7eR6R##Uy=x+J`=~jp;am%Fc-STm^YoMWjGH^I_V_c1sZLn8xD!cqok85{V zsAREPcFthxa-_<>>8K)Wq*=#E-r?CLD5=667I=Cn(4Wy?yBToHRTBEpnn>Qj^Wi<+ z&=952lKZ}1G1+?ArdBP7>h?wtoPHii#Gmuie5<_4$!J z1#|Or+=l+JW1SUGlCsU#bD6p?^>aJb6I)*8;&Z~zyw?i`lhsN+-@csSKQ%^%fNs+Y zeJgPxzj^{yFPzX}lbELr0=-WJ2s+3UUA$6mH}1x9VT1<5S3#`~@(iwUVA0-$UOEsk z7>Jiz2pEBAfge&mtRj}6$J_mi2|LKG*Klm`W%4&j+S%Wk9+bg^D8{UkF{9RlW|JqP zys82_Ka^2FQHnNR>D&{Qa4bALNO}K>F{S~(&miT>;OpV6;RupuQ89JgfGoXXxPb$5{-0iHhnO+ zPLl}pv?@4bP1WR_I~c7Dj#l8T?b&8m$aYc*=>bIM;JjsL%y2c48g=G=;Qo-X)RbgW znBOl^WIVxW%Q1r!8%TxEZ>F(K3)ARI5Sx0irN(>EmnI(uj=gyf1K&9aA!BT>@_?+= z0UCLk9*{QoKX_5OZJVACT8IbY?Xhj6(eAwDbZmnbuYkO<3_n&7!wMmJqO4zE$>e;R z%XhL<_?Q?U7P2|JKshPeM^g-!G1S}C(coQy@uNQKr9u*4X#GovhBUaWfEOuGQabGe zc%ES%7`d_oG+F{_EJM&6HSmpR7fu~&?j!RV`)-qcHjPHiV0VYN3?I-4d%OaGZnjqh+uMO&_Nc4k z!2X%Cw4fR9%6ZnH2w=!F8*u|7 zxo|q%yJC9KfS8{-7rTkUZE*F_rwtpN7o*Xq+3jh4wH54V4pa=v05&kxD{xA24}5vN z{miAb(D2cUmBnN(2gkkkoKKmd*k8i7oeVdrqfUj-^&Siwc;7NJU+xFNLN_w!PAkPP zO7EN4vtuOX@4Zc(1uNUG+YL^}jc3_ii?1{fmmiVPG<{7C?X$hS_xp^NKd-Bo#0N}x zik>a){#(e*5| z%iEU|%CQSJWe;3sAHT}3bzo;QcR}npM998dY@che#$)={6P9UXn&z+5+^7U|w%9fH z@hieglhJ^1tl3R-$FLb+rG{8iKfLz)jb$<)vivvf8M~2B$yvGFrq$6p<_A?TNDqj zpYCkS)%nf!G-U*Ogjbe?ac_L0%;?Wd_YyVt!irXn7EKLVR^7@S3EN~RprO3r47HPY zFoP8=%QbkMB2zssh}1%B8+wI#J>L#0?}xs=*KQdCDca_uvimEF%WxhHMXNTp+d;#+ z+bs7s2FK}g;IoXusglB}N$?zy`5bqEdVaFW0$!6lS3q~59j8Q5R&0_A+)=xKHcrH` z&2oHiKnO$Z2u(9rKG+|VoXKwV%^lf80j)_<5o1W4>S6gQ7V}bA?%4TnN8I7vQOH4SH_6^Un#(32D= z;oEldxh}l=dp0Ssd-BwYPS4KQ0p+`H2ZPo7{1)QVQtuUm4n(Ujp_0#9&nCEuGcLO~ z_t7`>Vp^qYO}>U248NZq)zS-^p7V&-9Mw^uz44cxe3YI<45@{=^c;=sUz1k|ks+1M zyF`ru=r4>=sZ|{7|5pBdU@|EC3ghbRe9%z(z9x;~%B;Qgw0>G>zQCxe8fwy=8uSk5 z+Rc`Pb4xOk3xjaHbPz;*_Be>1u;ELHM|HQpUo~m|op7k9J&6qK3X$-R(%;mUzq(BM zwi9o7VE$v1MPLlx?}M#%>JeK1AGwrbrl-tZ!^uN0!&J$A6B9mi&o#h}D;;Ty#p`H` zIw6o$sZ=}^JpGQklwsp`Gec}x#XlIAecSKb&y4Ohl*NA^(N~VF!$$g-9o5}fc{Y${ zWsNC&sH@)q3ng~J(_Z$_vk8HAAH#Yp13V$HE-_XVWQTQ%Z&9c6s_h*EWA|;F$ipYWhtm~#ZAtMH@bk<2HbLxY*()e`g(bMZ} zINS<1|69ebq!YTS;b$^G&CYQO3HfgNsgoWyaXc2f2|yz9K<(8Xbhz_$XC=$pHG?c= zeeu5NNtOL19etY;DJ+S$USFcKbuxKi&pQ)qm%6^IL(`+z-M&zK%1c7f$B!sMAHYl^ z6r$p9y;$3uK0#g)J-8@4E{^)+4QMy9(~C2vd~?1UPnz1xWIDRPN&ZynWpfw5#LS=f zdRJY|Pu~=tg{t)<|6QLgKiQH)4I8%HS4$WWh83OKPDJ(5JNxcviN!41VJ3mCec$vb`^N z#THNaZ75h|;yE0Asb^n=+}fNZ zba>*G(OJR%7RrCyBAWC_Sm+RTp4|J(!Sg@+y9X5Gh{wQTpv{0pI72%KKReJUGW7F* zdWZ+|w_{uFdq31Ou;Zz_8C)IWLAwUx)g)t39N_9?S*Dv$lh`E|4=2Anb2%t~c!A91Bsc$Yv<=bO9pxk#RSaoxlP@RLt3l2Cu=VK2Zqo zfM0UXa$sQ1)rF&tDkeI2?b6G@Z;rkj8iI*bW&K4uIhw;bpLiEelGfMYqz1u%+Mbo~ z9_>BJeQqE!hz3`We_5cYU^!04OCAE64MMNbn%6wkjRc>8smbqa08)V6J;--!Yg4_s zycFiWqXg*Y!w@zG0uJK8Rz`U-&qJdcCwev#_L@=Z>Wm7`ohtC8W}%mnu8wo3`&C@m zRY(-Z7JmATjq;4Iq(XJEC%)M0rmb7*jjSIoz{qTK~G0+TPmI&4hKKL@h3+(V;pQyoJXCPSXAu0B=%Iyo; z$#)<$m_aAe_MNqkFGgKpqb6{(0tUkd_U`s@eM2MsU7t=pdc-^st-V<@-x|1J4<1jsGnC*#X{a5FrNDF9#jVoNHtEPPh@G;cH=tF;EJt&;+XUXcV3BC&Ah{)O)^m_8FwKB0M45(|+s7cJMq! zR0z|!&O0+s=|+yvXH6cxIswTSyPRX2gSIS2Nvn^tM+^F{HM=cddvQU?KyONKs7|+%<3E;f*RrBW_@T6RdqlVMfS=Y5c#O!k}D~{LAg;BeOo06 z=F3!mRP)m@&jV!Lv;PesxHFwW$2rL1Ip;`|OObmCKj(*|e6#I8Tgt7kj(alkUCo(u z&L3n6k}=GFxOb|9;(BEIR zAvzymNYB*k+o56SZ_h+!6l^(jW1%w$=^#!sr1LOx+c$ze>D69{N@yE?=eGx%Y}Gl^ z$)0~~#I%SgJNc}jXR{DV8}gANeu8U(FqQoIHus|}7HHS+{=Dypt4;Ro*Mg@!fc6Od z+xBI3R0d1h>gx9~b=3Pe^LZg0{oF2CToQgi0fz+uRKP{oq7^Mry%4JWNW&_d{|F@E ztBr<@Z2Q#DV8~k*duQ0Zs_eszt?f zm0uCGQVyVc+!6_BS#Tx3I1fyRU~xrEg@ryT7@-yYlf56T;T@LJv3tUQIj39Lnb zd>`D+L|cBn#p??9j zw1$0@^<0lf_lCA>MYiM+{%zTc>wLTf<-!i=900o3S zA14sCqM? zUFJRFoqSD>*h+i?>HfnPm4k%yLNvqC2be88SaYEJq5j2>*^9i9!phm9wi560Y_Rbi zm*XxE80=!0bv{NzyqtFe?O3W$5*8#BhDRsV$|wJUfuaL}=>??;u(r&tvqH%(7=b#e zMebB>pa=d_@$#w;*jXdKEI?$()>+GGPn|?9yJ<-i)lfcj{d)7&;JUW90StIn>zBeD`6v$~}so69bQ4kx;F;84~Em}vX)b}5OdbFu+qGU9BdLTYo3x{pZYLfr!MKKJrNn-C^>slDgP*g=n5*ho8uKTzZJWk2LYPYF4*Z_Zm5i^hG4G}d zCq~(ljDeHkY&ujheDUh$O-SYnw;ydCg*bcQ^@O!c9gpOYBcz@moUUB%2jfrwWY$Tj zb-h6PX)8uPD_DfQLUPELeo|yxzm_1GTBg75oqFLg3-6_{MgXP- z&sM4FjxHw9i@ufHAckvv=lIB)FjRzjC-G984w%QtDyKwVtHB}>a^h7Hn|2FxYyd*v z-{`p5Ax3Lt7ZOtI7mYl~0LVwQTGR0{e_bA}MeGk!KBbCG>Y?{gb!0mwt(kwi&APpV zaz;JpBQkqzf$`HSTT%1kN7mY-_Q+WP-QZjb!9Lo_;^y}>)_3`k#e_YJP!)Ug8i4p- z{D^kzCPEnZtRW&cUYuRmaR~%`#N4R1r0>nt)s!fa6ZJeYa!L6&M2Yn?wIy$@R$Vvp zpd#*u?%^8?K+)|kf7uz;6W7;kibE00f8jj1{%tDM4y{|`+CG1Wq+XO zJ_E4016a+7tc5EaLaDeNz-U4-@8A}UOgX3(A*i22$$}^PP!y$_R8gww;jRVbO@miF z15KcLI+h)mM&tz{ge2KMv><4A3@mk$sC$R98%|Azq?EwjyKBIl(L4$lASgjWv4`>M zOMrS{kzv9GB!Ld$=xLRuzcy0DGaW_-6<&31c7Z#4IY=F}eh(6{77mSqe1?gFu96OTaQ)?1oPfuQfUTLMddnis zpE8FFg&<-NlmP~%ECv_FtLtG3FL!!JUG2jtsQY}nnLYj`sLSwbize}lUd=t2Gllgi z`UivY$MT8x+iL}@o}*35KaUd`0C5PSn>6Aj`jlCJ+H6G87L zzC;(j0jo1YeJHhs92Q=B!E3RE@B78b6Y;7r20pIPAKDqp$N|kK^s`Sf0%`YStsz~) zYtGL^WVEF4hskDmUc}=9z%WVROKwxvU;dHPrqT?mfwsq$+2T$Bv0&{(1a>$POGY#X zzesDcXz}5n4)$+;+5XUynT2l=xm}gJ^TOi@UZlSfN24?0p)Dklr2FTN!we4s@DE68 zXW=@=m_~!0je+KkklM3d6h`W&0A-_J_?&9zSTQB!VXmu4(inSA2>Jb|(PHNmTnEZW zx1?i)vgSl7jk4Y5Er=VhYJ!cVF|v7U4tcZ6f6i`J&hmAb0k+2DiDupj6GeF$1Jl1e z9v#uFS~Aw7@kt*RebB3*T%U#DIEB`7R4YP31l8-1J!tK~U!RvbSUXo0JYz{myp)rd zk(ApCy$nNK51x(?;Yu&5#ER^AF#)my26_k9;^C28!LIhdvL+awuUfrlWV4?Sp&&MmkxV2so7+a1q=;X3h8 z0pp7k<$H`Kz3{uOc5V2RN&xs?juu;K;Kbb$6=9q9Uvh6143kW^?V`@5n8;^4$G_<1 zop+uSL1ud}?qI9ZeEnasTG=5AA=uH3=_}CX2xT)EogSjNEFp(f;0;k9uBaa30Vn|C z)7J}!<#^+e;JD{D7TQxAZs*8=7C5{fAQBw9Iamr(gm5f}%=+}crA6xSwO_e07x-B(1WE(fU6FIT4Dj=$6 zvK?EYRA4|3S%1&OBDCv~>!^d)*OIp{j<^o$m%uwMhbIj|i9O9n!||YlEb>tpI5}>RedSx{HUT>m(**C%a zpIV0p(1W|QgbT>BnqWz7L;{uh6!AV(UrB*fR=XEQ5pom$Z@I8OAqoHuKa9J}rT`xg zau8Zd8ovIargfYhWcyq=HIdn)LaZ)OAoa853j;R3!u%}#G7_I1Mn!kS7%%+$RG>Jr zR@JQhYhB@X7+-_@a#>_%`*bQ(z+h7lZe3kF9}J(f8AQo64O>CUgV<4Tb9 zXnPKW69DOxeU{brU1OppcE$K5GmFv#_HXps0_h&Wa#1X^?A3*qquW0Gg;>1$DLZ*} zf8N{8FudOUs#(D7%B4=ESVxKIwi1~&!>g_D-|k#@#?HYAW@CgLh@^r})DmuEOBRXi(k|wm*Q~|t3KU!s>%HWui#0oeDT+#mXESh z7rHK=-0yNNgE(AI31Jew} z{kLDhcS1UZ4hIlEju5+h4QsTD^E>+D z%D=>64l!3pO`if9@1=z1;P4VE4E-a`Gv}cwY$h2q8b;Juvikg$9xxcpgV;z{@R(Yv z2TIKT!#>bd3i4EgoYjk0?-oPTi$USJ-2u~>2iDGjmX(IC*s+kS1Bn0-RB=_vhgul7 z*W=t$Vf;Hnobd9iacb!F8%k=H(r-gExqz%#=0sUko4Ibu^_US24u@*Ue!F1b2m_|y zBYs4qHWti!AAjh<9XP8Ve*L~&$PryCXdpv>S*{0C90xNrN|D#62fp zB&T#X{8Q3cQR1b%gQvKSbu%mAcNUUPvgbDB{<2c4y<%Lz9|Ro6j0jG$3UfbmmBr$=Gr0v z%7$G>?=n$S+$&pf`KLV|6HGrKu9GGx!lmIQzykca)z&|0jy*>KNDTcwBB3m~)$XQ7 zY&&MWGK@SN)?tHQ@iQdm``r2QN^`Qg+0a<=!Kly;3@gZ#1RpuBJ5v<7ZWr5j zivZeN2uES8ToK~F&;U|21;raiF!K|fDF%B(D%)&4{txd=oyG1O8c0aG+jHsamUmpV z$0-U>SI*R|v8Jf|EO~LnO8*;_wOU`;OaU$I3YEO3oWI+7UZ=}lqtL-pJ-m5+=?l$o z*G5xa%-=Irugx=WnrZK7=NCx{wgu6+I$z$IchN!XqTI2;)SJ9~BLmGkj;2u(4_{O$ z41f#b%xnL=`_GG-q!)(e(tj}gbmw1B+EmOd2g3ePW)tok6I)^hLsrl}qb~qdYk8|t zZ-n`3vcx$D3kgk_WJX-3glI%G>>F=8+3$zTqXstHX2?w8XdlOMB98?agi|r1?x}^n zh<29Sb1E#@ls z1)U7q_0}X<^dJVI16#Fdk1cfO-O4NBkXG!FEj8wGr`@tt5sKbAClCsu)tWzWzBb`U zSjGoFjPZlg|7n$?XzhW0HQJq3u1g%p?c9`S>oO&Ik&etgN(Cj4chQjTz=6XFs2Uj& z`*C(s_6jw1>FBq&Yvw6BA@#x3^+L3(D1Eo2Y{^I1D0%|P)+U(yyTHlP%{mo&p0EeB zwokI&aXC314jUKCCX9{m4wM{|Z`90-!TG-S-atxHk`hh(b?PM>jr}xoeDag>Nr-4 zmCyOo#)((z-#*SPy7fKtR7I(%oxS>I_y4S0cbN|!6k%gu-@4~2AVTN{r6b{y6wLp) zVfjhTPdjmiZK1X6yGo6+{tvhp3%d76)23JPhvg4dcC_#{2h|+Ygk}1}nXM zdc9rAR#%5^k=Gm&#U^&y4h#dVQzdR2p-bVKNz4s0zO-9dnb z2h&U90}z6tRh_ZA0%~Cy{Y&yz-W&eSsa~ehBE_mOB zzI0y;_KoM5qFokG<@AbBZ#68V_TcNiL^rtM=Y)EKRKu3p;Pgm0g(tfTIk|#oq*UMe zjaA$!fAop8PXU1IBgLt-{8K@$)wheC3)dcc-zod=^TI0h{?SATtZV_rEa0PVR$I{l6hFo;w`6TXYC63IbmC)D{Hi1*--rZEtGIlBFFbX9eqC z!IZ5UiFDFDj?5@9mb6O0`w&omNR5DO3dn|dQUNF%3{><*G~Jotqg$?(j^PQFIu7BO zy+q6+$b%&Sx!?dcSMU%O%fQ8CE#zt|%YAn@F4O2pzv76i&0K7HWZ$wz_b+_vD#sf2 zeWsH|f2nS#7`$7)N%7s!ROw#^G37F%*XBtDXX>_f$9w-Sx6jmXZYlN7xS*5W852Mg z6_QVno>b87j#I3nTB(IoEp@(DV^2ztBTFwjW&hautS3=%=I?Nq7@uBs-<20FTMPPp zygIAJYl1T!(h_+;GW_t}r@4jJxB{_yr6n(dSI^5IP~6tv`XF^DZ!+GinV(C8hTpsk zR*Wq+wJ9%kI|``8=vUE$KFdX{G^`5Rx0M$6l^F{{a1(j?mhbvkiyVX|9S4R>H=AFd zDiMP_5T%%5Ea%lY(OoT*WHU0hm51uaXJzdAEw)4&SZzDMPCkp1yOBQP7fHz`8c!A27tf{J3>k-Mg$ z;Hx@4OZh4LDp`S@khV4C#?+LSqtbttWtUj%sW<#?Mdi6vYM(YaKzx7a%dL6C zhqwhsc(uL3%^ar{#I26EE}$t%U_&%f>U!%vzseI5Y_2!E=jj&u`x29gJ>Wj`B%6EJ zyk>rN=K51z`bv3<9f?;Cbyr_o@=#A`v>Ag9K}6%fD@ENnJI;?BK0W3;q%lePqw(p7 zjDy|lRMo=N>)OSRW|qa~DvW6V49Wwd-$r!qmYcfB;!X1r+(7%JWy)*-MxKJt=FPY9 zKE55%5#%6(}a?umdh8FIjzO4O@++ z25kx|BC{RnpcfZz`>URNHm$~8<`~Oc{4nBQIkCVSd7hi(AG9-Y^St+c@c(M;TEm*S z(lD-NQ7YnEdBAdATPy{{3sqJTVnK@#kt$N45CS5iVi3u43y@e4+$w@qAP8u5F(9`D zgn$G??6QCY3d$vrkVGC!SQtqF8xn|<>^HENKmF19lX)h2&YYZ?`R1JW{oe1bv6h7| zd$j!#CPy8{bX2Z99pDA;u*!Qlf^Ph300FI{1ZNFcvxN#DnJ=woBfuM30+3l3~X(DgB$Jf+p(8k2S^h zQR?u~MkVzOg7y#VJ$R$hpZ|Mbp?7N`UbM@gAYv~Ffi~HfNSb9lR@3ZV+j==5QTDME zAf2CG+Tmy*2c6K82(M}&eQs^r$j~6L-v6&hVhD5Tn&0`#M_;;x*fl#oX$rm&u8=)WF?_!+xv>) zHBPC(_hpZjheWf)T-!eU)$R|%QE!U-iIJCQlNyPmcLUIgQ>RkW@vG{h3F@IH5(c8P zkN`S7KI#`={fLiOY(M^KtKTOBwJB_e83Vl=Kwd11g>WV(1xSKe9|ANi7L*U{HAq-e z@;+41ny%I|k97EZY+ev5<{Tg0=AMtCW&5^iGX?m$&9XE6q%%mAB+@6pxSc{#?u|cl zjYFjUsInV(cAT@aP0MCVRu=`Dd%XVm-f=$l2)cAzp~K$(bhim)Qc|IP@?_@pbwOXs zD4Di=(ESS)>N_py;}HD$+ATe8HW4ETtsMo&7{;AGMqD#WSK!Va_myqo$t(?*dceIUkpphnd**O7g!D`(ZE%Z z8$kVB&o#)e!4V_9KoXT~eM{l7S`U0akkIcHn0G7}QsRK}(1?{39x}PUu6y=jZEN)t z4msG0!=w_*y=EQ~SL1iNJJRREgL5?9QKf-fBf0UL+Y6$TYyJvz(dDMi(khY0ot+7$ zVEMJhi5L*jygin5habxjTOP)Nks)A zP+jQBi;D<;p<-39p5E$t?}J~Tcno#4#LD;M1rEPdfOLo#rYPTaPV^n$sXr(L0nfA? z^uBXox%j898ue+OvgczrFn5@+>XH=dciYrYLV-oz$G)BIk zIUO%riCZ#K4O)2q{|cda0M7Li-@ZZ#vd5FxG@TH8o+e@`P05vKi?6uw?kqKnx!Q(c!@fxRaST;4HsId9g`&xa7KgN68$8& zoB~AYbb|{K2kz~CcVgM8z>EQrO%jQuf&|NDfG4|WVtzKg)f}Icr;Sgd33NH_lMr$a zcwQR+OJ2{1`NRoaZ=8(1 zOJHE3=o1u5BNPq#j=RX}pMf$#;P*L%e8J zsh2{uZ6MAq!GBBQ{qJ3vCsku!Ln){so5k}`TP`T+4wrHF$ExG7w3GPAVgQaMud{Ze z2c~~=N}lqYYSxGQw5oUMzCW&JTGwX1)yT8li!Ek9= z&9shYnyajC$o98*xhNsKCFNqZ*e-=)GN3mOk*oS$f9@Y?q{dv{lP75iDFkR0BtyA2 z=il)rljRmPY;M@HTtggo4%B07bJ&5+#sq`PvF&$iVSYm%^S?tiI)h&-lVbWe%?hv< Ny$>Bd$npr!_&0JJzK{R_ literal 0 HcmV?d00001 diff --git a/figures/single_angular.png b/figures/single_angular.png new file mode 100644 index 0000000000000000000000000000000000000000..a6d8620b85e703fddc1096e47c811433077f1037 GIT binary patch literal 16089 zcmeIZXIPV6yDj?ARH`T-C`D142qMzEf(W5W3%yG3AP}Skf>LZ0;iV%8(h_@? z{d`=+1D)=;xcCOT`vtAgbgDwoIS6*+y6N3q5>Y>;yccm~=&jNHR+S+!<;QZgL#&fx z%8Z7C%b_=420URYgreWxqW01hsapPY%SbS0yCK9?DA*Q91zvT2^=!Bw_zg+Nn1Y{$ z{6*9d#LCRi0YNDDV_>Av4GnEN@VYbVqTs_PR%7tggScccpKJf+BWmG(T}Fbix-B8_ zD9=q>TbH*YSf&WnFWcibbgaR~`m?AnaK;#by>{`~ z`k&h5|HJztHz!eka5%@K!pxw~f+)Ys{4vQc<3;l{5abXFv&C&pa<3&6igfQsMum5F z*at4vx>&}E^&QEVfTc{%(*j!_&D{&y*73ZcKm5sXe3YfkBRj58L>&$vjA|qLs*6BS zR}YIltiMi7nJ+1Mey2i~C;(GGR_xD|8^<$0*j>eU+ z3xn&d>X{=qnZUPq!1F$syuiqWlHKf*pVKiqD|zOE;FLOK>%xr*j|z`jrnf&U2r<|Q z7I_QSs^*nO{0IW2mm)B+vbB0wxqD{q*cMqu0QPI+&Q;19bRi7bj5W5U(i_D3lIKIU zhoxgy-E4juz6+{I!7KBDMJDCC(q|*00@+k^TwM8*w8}7^U%u;@NY9wzoF^`U?GTKw zDmtZkVxg6f=?=+kN(PU}$#0(n8#ArL!fRaWjz^Hy7cBgJgv`N!9F$6rqxB#ADFqFj z_G=A#d<7iQH=tQLCXCQSn3&|$D{X)nDGrCA5TxBs^!$hWTY&kU4U1{$70=`(y z0&`Lx)HoKEp7nE>C5N;~>cX2}^vSCZjRapeJsmlkuin(oq?cVXI5f^ue0oFqjN%!oEj}fF$`cSOD6@wjKrZ`5uzz@l-GJiKw*!elcCJelx|K)#ds=&qjB+XH)A9tEjVg*h}zzO3P=@|nm=-j`$0c|}oqF{x+ z2Tw|~(t?ojAXW!F3Vkw!?RQUPIl>pghaP&gJ*y&ke|j;I^B-^RCK)Jl^4b%MU*)k? z%9X`2>r%7m06Up}ABQQ$V7ji;|I4S#i;fyk#{jpFdU17`U7r7c`+O}~Q+>D#9ZT;? zX$aC1PNo(vOaVjvPYxfEY=L8eps*}?df5?{(*6y+GZ@KgQ3?b1X-pmlLk)G=tBbDY zonwZ&ZZ?p+kFZ#a`}l@$Z1z;!LUZRSh zw1<$$5cKhiuo|yok?iKN$LgY`l+CM2R@m{_Kuo!6v^vV@58L~+0n=bre>t@xu_fO3 zYgUy=s22-YJ!PS7yLA*!#Q~u>D(zs*UuB*f_ zB~#RNf2LHor{OHlT-+zU1Z?B5xqp8phV2LjiM}-YMbF3{%RS)sZ6{LCVR!>G?)SG6<(WLRvvq@UoDc^m98!#B) z%K8+`h~v?Z3-6we$9w}D=^p`83Dz8yVAuRaSSCEV&|kNt@Zs|it+ta_$y{hhsP?mZ zE^*LVTK|VzIWZ%warYP4iTd90BhvaImKWf_fKR0g@XJ?22HpmdcyB#VaJZ@#E$dh3 zLRGgP7kQlNt(y>hP7R5jhaK3o;;1fr*Xwwi19zzV;&foPi{)=Nel*w-`_qg?1CprJ zhjqj)Z?O5(m()f5`Hh-@@va^=Y)~`pjxQb!At39;>Yp${{>;JRWz2=4A%)0wF6F%qA1_5u#R7P5aD$1&&?4irP8_9T|75!o9kziin)GCtA;2gbm5j)LX z8AToj-%IVznQ4$q_fJ8n%E}H}eX@WYFr(KEE$u5!%QJm-PlPJ0?!`%+W_+<32tgej zq?6M5PsIDz_sX#}Gp4%ye3MB~%^hu@_+IqPtu}l;nE_p+oMt?hFZmT(#kukaKa?dZ zNYT)Ao#0|)fvUm&>n5GuZ#cvFgd0(1osYh!)Z;lB_;_n-^(=cTQ3 z6}Y4ak0dj%y)EQU1$GOWP%=_2Hzra4^Y-y0uKL5`z}@~H8eq(Q@VA599uKD}*>*<4 zA(SrKtDa0gTV@4@lf%N|BIk6*x z(r`;?x*Jb~^#?n0OiP>jY}pnv;3IR3-J$7kadCrOW;Fh1%9CI|5ms~wK0Flnzc>y; z^VZ>tDSYFkFP$0m^2(?&O4p)2C3}8neHR3HSc;zKzWAf_z*S{Lk%N8-LesQQl;~c@ zjN)L-{7jcSKa1(6x{OP`M3k2%-1=f+n;L9P<#5z#)HJ^j`jnlTFH4W+-$UPvRJ<(H zggSE%Wch7B88ND;G^wh6*nawVgess z@32g_e9QqgIg?BvMT7BWYxrnj!l?4#(fe*j>|Y~0w2qcmg(F(vdo z(c>uaNg9^v67|P#g@RB7r3berD+ZfQ3q1$chQWk2?KTa!f3b(gx`yNyF=bFj5GEOK z``J}nDPoiEQ%RxX2m=!&_7PtS?RA>Oj+AI_*;WTiWf|9B_R3%Yi{7OQ(xe>Osm+(4 z`Slxvb-Y&F#-qW`TMS%O2X2pI7}ENj)%O_C(1|*4`nRvcPguMLSp{%*wJ!o`W-SDR z%Q(1Dyn-L5Qon5mJ6KOX?t1hm;gp#%?5msDa}(y~Nv?dyuzDDm*aEwq@~r&jt`~3L zirPjySKlUeeBh5{JYT$G@Xstk`edjk8X@9YW&cxZ+r9JfQxoTd0Y?qb@RlsGzMIy_ z975~$Cd;#pgMSr!Lt_ujx))o+=v11f4vFECO(Z&&iwnuSgKCE{hdcR)Kif(6hr77- zu;E{X`N&Q+Y=vyvptQbl{E&)h0~9pU_?$s<)s9w{HCUT>eI@EpX>ZZ~kl~Qrc}NHk zSBhcLzbEZdCe^Jq7HV1_sCO+d_0_ZJMSgJIa*(mBcP>UhwYe+oRR)Sfa06Ca##>$#w^};d ze7Y~Y#6#rKH<$iV)*iPN{gk&$9Z>C=@;Q0L*+F5B-%<+rm8b?1ovfs=nx-wx)+F79 ze{kKubLtK7oiaNq5Cin|<*wVWPmzbqo1RohToru;gA)C0A1wq`pu|TA+b&~Ve z(ByE>O3jZMISz8F6T{2lc0%2XM^>&sp8*yMPJ*x)1|!O@>Da;vJB`mvOVoo<*S7qY?Pka1k0_poBcs-1o zt%A>(03ULl`1w~?mdm)VAySzSA4w&jaZ&gnP_(k9!LQ?-f^pye$ydUTk~s;$uY6?5 z0-;!WP%Tm!x?U#nQqB_G`aP+vsS*D6V*9SOw8LJ)$P1<$IHjP z^%hDiTdg9f9`t?@F6nVn&wX{&Q6>j&toH0#vx?B@RcUC+b?}3csnz&>32oG~3Okt3 z-~hoc(e+WGh+R)X6O9)W1v8*Jchw}RUr6gKCgS?lWzXxAl1T|*+HJ3|!cq*`-?_sZ zLE0n$1F7L)lvg~_VkB5)k@}i|*l8-USq9KLH9>IW$WJr|GY%hVlGcC0os%u_aGKu# zhLf;461x#1?4-(wlFAAH7()eW3?VjQ7t76n=OFa7q8Pej<>_a%O<$te#=H?ztOpI) zbQcAWi&=l%*rx;yy6LXmu0S)qq~-nn2hY4eFDn(HK_aq|+p;Ahwe7h@Ctr5U8n?^W zjGe!hQnvg{M@h)5E~=K2+~ShMoYh6UuJdlhTuM*Yi%4v}709$Haq1*7(-eBpb0XtO z4K0ht5={l^0alM?#chQdLiH=Od=p!+GIT*EW_f!*O>Q0iepASZdioA6D|C5oQ4&2> z6|JbH=X5OrM&}^a6SI~m*>~;x+;?H%g#RuLeg$jy*Li<#Ersrjp{Zu#dU(U8Ab<7r z=R3&=dka*mF9^jow9Kq(x}Fd;^e-@}Mr9NfHb*K{+x-5-*N6=s(!;Dev2Z)o{km;s z0F#;J%U`ryyCFTe4)(-zP_5@JY*lbk^`uh00ZrF!P%VRMkoCaSvK(Wz;r9E4u_Kl1 zmf+1zUC6@5l)sTaxtTzkof?)#HZ92>`YeZcPT4P8r_cT%h10_MIi_Pwx2YUBj6j`d zSB^HZ4gz&3eG{=iB7}^OF-5VR1c8yms4^9wH|;*a^%*4fkq?(|!2TF|{$VF>%oH-{ zcD939WA(8`k-uKxIfm=RqGCtqcA9;6_MIy#doNx4Zc6Folc_PGpw@OH_XV#G8B zUz@ls2#1|PqabHZT`i{O>0PKNBXCm81ySjfRfQsvaC<@yzMZIuhPnzT>A2^qu(vmT z#{iu1yZ%>Gj98y~c8P>fi@v*0&r{jx#w>Sqb!_uKV{k^)PSBhKgt{JFNqw~D{^%ui zLPzfpQ(~`|Qju@Fs%JPinELd0)2V^tm5twu-9m+y!i|+8j||giD=hD>D<2QG-D)~g z>T%|m@)V^|;B>Ss43rF)zxz@-mvl`tmX|fw$G0xl0dyR5&{(Uc~n-^)f1e zxI?2p7MY?1onxVt3jVxEs^Fz51^Q|Uh!LQZ(>h;{rNDk|g%WU7(=_E0*N0`LS%L>< zs2=2adKBTSbpw$PHx~|Q-TCvoxir0k|+KB}z3 z3g#0(QYi9e&j{*DF7RA17M?R$rLq{pq?x9fKo4%%9(LG_pSaW|sc4c@fut_1wbWLb zIa3Nt`Pd2!>BNLuhCmdJCwcc1cn<@t{ny7+f6Aj1>hd^&Gq9B0g)J&5(^)K8QI7h^ zM^M#Oq~eDTfDLj6rt2u>SgL^7>Gm^=RsF6qAAet+Nz(x}l$%#`-Kn?3|4Wm)g2Dn41A~FCQ&1RWuLy%K(-@{WOb8-gFY(YKu?zy91Zn^{S8E-1%LBi zXv3bua-S3J#s>F?l0cgqGC~}vNDcdMGKZ6Srzga|UXvCZ?N(B)-u^s@S9>njSB4Mz z+H3kf(MhZ;jM-8@<=7cLC$*M5h1c~}bnuQkfDTwH7=70+Y51$I(hpN;r*uj#_{r(= zJHFR{H;{kPkvh7D+tQ6Uzux%Eor#x`=Ys5u#(a~{rv^B|HqwDe{3Vjbr7XexoQq?8 zS&^&dDNyiOrjyz)`ipcVBoT(_%h*h0XSk?7!f~AvH8mQh#%GhTlQ~elDX#JPnU>n; zo??>aQY{~o0+*q>H4*F0{8X z!2Z0SdX^#VGjC}o5*xPK=34hJs|I2-U)%Kv?bLO4qboo~35NsF3%RfCk$tmq5<|}n z-Lwcp)mus)6NH^yU{xpE&hWg5pJCv!GSwXlqVcp$2URyXQhlAxgz?@wCXZ}YKM-Uv zCp~(t`GaR3=ij4mmTV+rin>EX zqoRHiKNNL}Ytjg?N>nU!9uhjH zolij9-LCcIj`Djfjq9g0*EUt5sWS{k_e5<*@GD4c`c7O2q}0plWEYqTK!AnPEuVM- z!h!DGZ7sE|n8q#mO6(wTX}U{dIW+Z^XN*gq@cZCHGb{>^OtWE%FDZQcAkdZ4lip8t z#OKShQnB%F;yK2O*8;GvipJX+impH8p&yNwf*Wghe;*lDAiYASN|4Z57A^~P&~{gL ze!+^BByRcGP=LSP!2YyTxw-MWa+F;s(HEF_Mh%_?-G4{QheU9)x+fHaYqTx3LEU9; z*J=-O)Z(r%Qe2#JWXV6fT1zROzHZ&Eme!A=8%11N8Im8({G!vCcpj>7e=H+e+4cdi zI*v_wdw{*DS#RL^Lagr-q0Ofw2|Htc4S)&RU^{qnJug~kN-d9xb*(nQK)e(K5W_F! z73Lhv^n@eh%7c|p86fl8cOP1GIIKZQWwEI6YOIJDl2|B`7cJdaH!moQM^ZUx;*l*U z$=~iY=xTWxL_WG)#1{)nv`hGfBILklshYTuQDgFq^kmHa04V9m$1J)`bIedF!0jL?oK0Kph3)QNz=8im z^zh)G;)~%(;ZMuJVJxlWjNf%)D)J%u<86|VIUzyW$yT6#JE;y3=-5ILvv3|z{G#fQ z%Ikjrm_A}+Zkk|m)mV0LJ>p6b0E+=Mv*PmyKED|?M}H6m+UD18;_7p~?7xQOIv7MM z?q!Z3uJu#O^wWQA8o+>*L`>AfS_u@|09s-S%~SIOo4DvJ5Aj5=nV@F6=J!adcYf%2 zD!Fg&pZVkXK81&~`I1E|{{ZaM0APQzs8su_B}6&Dfz9)JQSD~Hv@J@hFS0DQo7#vd`qvO0ZAD zNh?yah@hp=p6*Jp#cXA*@8(lsR(ea+i)b(6MJa$bfy~6^$$VRlQc(CXwX)!Z8C*t* zmZr;5zcM8kDQobuPQ3hwIlaFd8EMu8{F)Wj2652aT!X-N9{I}#@;_~2P}D$!UCfa< z{k%5#Ic9MxBXen=Z=5{8;;^CPDF>FL)yk!=T0g}5Vpc6ld-%6y039SBQ-CpbDx;2k zo$7GcYbEh%d(=UqEN|+rFX@vEKn6x#zw%n=2VZn$J>n4PL<;HtxN}$kXJDG_1b=wa z@FeKmdlaCW7oU*0fsld8bOo@=ihXvgsX?WjqMVoFm zp^4CNz|#2Jmzd4OlbOylEn?c`*xEiN5q`M$tP&J2(SAQyakSF~!4~n6h&DyQ5ZK;E z!l%EdeD}NNV}CSp=Fy5ZhkvH;%2wZ`T9g*7{WW_X@Zk5Ivr7s@9=*f66G$B-!||$# zxzK}4=eB%IKq8gIO&%72-QKi)iK%6?#>NJ=y&pL&Z62vRvw74ae{0+W#8iOwmT)w4 zLmHWK(MIH!)Qh z=bxfRwozNzcE~7suhC`qxXm~9uT;&W@=vIbxqmB?0-*e3C35rFUT}=x@4+u|9dl&e zwv~v)kbstM)w{7%Yxh9nJ(}^HO1_x!7UQj$fri5b_``JnwZ@}vr>*UD8@NNK<5fW0YASv7^15sO}jeJK2+K z8A|2ArX_EwpmA5R_ruoHgOWhrS()j*rm*~H$R{h;eN=@&5e~>KA!4&wpF2&tBpRxr zcKWm_BUcorzZvjbqUd(rhU9$rSGitEw5Cyv08HrU2yh=n-%oLX<}fkRRAE^^FXq&T z>`>_m_kJbEACn+TX-M9m47UcX(V`QTY!jv!V2Kqs(%jWV(Kh|3(oA{J^T)iWMEQM& zve3Nv)+!BK++Zyw>K|5$^<}J{Q8zTALzNo}Djun3tVsc>AxRQYfSU&C(nZ$M3Ftd$ z;8=50bT(fyCS6pnZ2#hF_iLVqf^mjv>sv{D6jsCV-T_Imrjb~qhy4lQcrV|$J8yzV zY8*S!NNy@GA#M`U_<|q<4g}-h37`Kpr@p_S$JW6Uk9Dmeh|iv)1$Z`sDBBCt^B@NV zx*nOlM*&2yEHEuv_<4`&U_K;5ki7J$ak_*0OG{a6;f-N4qM%hA6*PUJo7O4fVeKP0 z0h!@G>XQV6X94O0fdhO_u)VvOEATXJy9yz{o`Vv)+eb`_Ozr)y$O}5b#c)>rgG`pT zS1>Sq-x@tYH)vAm27wQKXV6y2=t~}h48C#x7EiZMnVdE`mSc&02Ff+sQN;9;C-n~t zU4L(>pH{*g7JU}yf+wP5GdjI<&cDi+B!*_Wemw=nZ~Zvn2^fQI2V{F`r1WAlW=Q3~_t ze~r_+D@-K0$`a0@J$}R1UMhxZe{wPeGO`dTYQm|L5-*atLt5VsQX#H3p&Hqgr$uVpa$RHZ;7i5ZS_uuVXIyD<*0s?~qjNHDsd$kvyL(l+Qz=e?#h(H5h)I%L#9yJIhkFcNGOO(WmCEhGK`>bnxj zITkJGu3ALQ=u=d0TVWuKnd=8Lp`Pu0qo-u@)Qb=ItvAv{+o>{E zfF}PTXrk^`q%AIp)MwXtwE)CDd@Nbzg5)8k6r+ZC8Tx0LTI zQ`Jc17n5!Bn@{5SWpY*T4va%Oz1JhZ<>>V^^ReVD?%ctrgq@^TqV;ZFS{u^ti%rNf z9k#)}Z3RDN15^R1$B&C$Hi-(irkdW?k!p9ptw*gN zYJ4V9?xvgH)~GB^mta5hQ}y{f4dlH0Sy`t6-&omUAPRNyM_dg^N?n^RlRldxgR&d~ z1dCW7tNIWo{2b^NWgEfexzM_-5@n=Ge4imDrGW#E(v3pl8hlS34~tg&Xwy$?G`LLA zLrF`U_!+VD!%_f+hwiDNhE4>YXU_i=jb#8#A;@wBU=&G0)>_DWfc-ZU#06mIE0!I8 zShK)^lqWBjyY3!H@OtaKjj!h>B- zFdrfy)aWP$UYL0m$X9P)U(^2u#_V`t1solPLf<>A6MK;*xV=Ony=@Y86f1=A4^!t7 z#X&avJ@jC0h6-{>sdX_$>4Q9;P*BG!JMZ31_lMD+Z?e?3?UqZGpn|SRj2b`a#qbCJ zTzf(M@E0_e#ilgP^h7DP%v)AyBeFBZp>(5YcM8IH3#c}xgpc-PFp-rBXZn6gK@}gxa%z&wJFyk1j7AoJ z{7y{Ya*-C`l*RfY=>`EWV13`o>VgdIVi7gx#xgAeJ@ajBFkmTm*77X7#W?WHFl`XJ z%M8-yLxUtnTCr>1rJE$mZl4)&$zrNZkUreZL=?v*28Z2;j_;og1J zB1K@3pC7P_tQ_l?BAfy9X$^;i#A347vyEB4WOK6WmcfGv!tvZ5d}Z_n>{#uK>KU<* z4RA&^o2Qp!IRU1TmvBl7mBvi_wN`?jH{&afL6IZxq_qSHy!mO(bdc)s`t_1ZPL#%K z=tf$9Pf_HBzB+(bj*vQ`id)6(W)4V0i@dEta~31meh^4xAt8Tt<{QrSf70~xLDZ}0 zFb5lozQ73W<;5E1o|?Sde+EsQr-Gu{|8)9v)NSi;!n-TUM!VtEvN&B$xUx6 z1JFEK2s;S6>o=OL=)&YNJm5t#g&aaK#ZVief*7Jx{G0cB`is9w)f?SwTn-^1C;-OX z&+yss06}i&B#*u=t`jheU1$Gy$q07v|6Vd4f#lHz633C~V>J)|E)gqAJhBfMFeI}{ z0OqE27 z+B99|JIf)jy4XVv4Y{QvrWEi_W6J>C*aFnLc)i%FlQPgfJ0_{|d zgMWNYKTQ0-Y2DY`5qvex1<0We`kFF*^##L?Akpg_1CBo))7i^k5rrMZp*Dk|>Aye( z`I9_>y05kB=q)`kl15)lll|^773}b?3lHvP`)4>)i(z)6T>jf@=%0T?p6 z?38Bh(h(}2Q4Y~3Zo9V7Ia&=1A&Tlk!y~)otcsRW<^XfcwTUNkZX@qH9Ot08O&D49 z$~%q@2fOze-Ju4JSwKqw;0HHIwNy>f@!;p(Rvmkca%-QFygQFyi!rM#%3eD5;|xXS zGT>jF>YP^Y^m*XOYHxslI{XsaC%cU+%0-1r{i7t=sBh2x$PUid6AhU*37(vQQpWPdDbBG}i9#rD=U!;t zC;UUY5iE%6}wq<#cg^)!B49p0@7MOLo z;qd;d!Y05E;YD*@#-~4aIwes^;+3sQLT?=p^UJ_fFJh}Vp8e*%d${OD-iW-H%EgGD zoQE=iF;z}3b9d*nTBgrxB|T#^1X(KB#5%G3G#qy36SpjZLMvrpW3F_{oHh&WmT6pQ zJqV^UP*>uH?D0h@cv-V{Z{1;Oc1=mEes6PkvS_|jvs_wk$CsH}0)8vtcE6CRAV^GK zf&^|o8UKZk5ubTAvEdWTD~-;ZCLTWK&QmJBz3U8U?{Bt^4Sr{(eSs|JJq5o6$h@1u zInr>9cUgJ$3(XAo3c|KTl6_XlH6{~bvU{yGfC?c9ElWmYXL$AL|VpXx_Q_ruxC zeO^ZzD$W>ZFx=HN`k3d_tm0vuUDAFXQK7E*MV%=J>{3SDaRihIY(6rvv6?2M>2BM~ z(A|OB@vt2Fzgi0r?xf+6))_+_CT_EGiQoIUKqQ z`Gdt8RQQto+lY4Q_dLv%GKQGnJL$d``&x500`d>?>v--p0Lwn)o0B4NS)*pv7e+9* zmC7!kHYBKZwoiP_tc5*Qt^C$~InLlq?iduq%3=&}H+C^P%6jGY(`3APrFTy>z&ycn zp#3L(nOJYzg{lV7 z`iG=R72D?Q-ToxM4oM)P@fjjbs+vXnsI3?}z+rT#?I&=P7l#vOmb8B)NXkVBoW#Cs zXr_X0S5TB$h|;6NZ5ae&DFhT#=M0eZN8@)s%F^@s8-sGVPnfv`X= z1M8vmG5$2r9nOfA3kUreT)p6@)TMrk8EsA8)F=aysCdm+$Pv2y9-ju)P&A0enIiP4 z>*iK%C!zB-+kjpM;qD^f!s9+>K0gjwP}I>H4eQihbN^_nPZle9P20tG(l0reWn+37 zpmK{K#1wXB4N)?I!MGAT`^8?(<5+I@Z)jA9G(%Br4^Ue?khc)JQ2H_9rkwh92vKIX~)})Z(O0 zjlp^O<3o#!h^g|=%WcY)m7O~|5k2I<^FHybn3c@7x8crWPr{BW8>1#rP=L(R?gP*^O-x?!5Od{NQHP z2FdM zulxrzZ)lBCAD+Ns-ygGvl(-n4b_-Hg?!gnp5=wr99^B(k~iueE&QLoH#`#zlpiVq*| zOY+$GwQuAQbbz*2P5+01*Gp3j_ag(f`dm2LHt){|8j~tA7c6 z0E7_K^=^NuovP~q|FxmT4!2>K1C+k8=%Vg-%OQp9Lj;PJmj=nC71nZA$LCJ#hE{ldg5x+JtZYx|gMkwuGviXCY2^5(gHPJqAj+em}7PaUjU zjN0f3@acbQtqd0zQ&Vtqn}Wxavt7pPjdy_p9%w_&3^i(AQRMV|DOP8~(V4E6vL1Cs z0x()jFaUa>YT)0{pcqt_Gzb(f1SI9~CjoE%pQX>B{eSCv%>My91hfA)uH!6;!cyAh zqTi2e-NmYQbCHj(j8Xnqrvl>DzcB}jRP92J5-+$e^q(Ey*10B|8`~clPvvi{yxj0Y zd_jKy7=^7rsHl9lvvPNC%1762Unws%YT|R&VWQi4i>;I4 z7ZOG0-uZzYxh*9Vc2Mz(w^Wb-gwkTS{;S(Q_2jx1ill%dCl_PatX-5n5s$$Lkr?|s z-UI_d=Z?i!u0R!&&EwdN#h*&2q%&HQ0+|#_TgVS)-p6^48RpFggR2_G>9dnL*nytS zq|I0_pvz+69|8DUd$rTY4Lzb7kA0P({_tzuU2ppUrk2E8nI!h-y2t(L=N`{13m_}b?H&JL|jn3x8~U5Z9)_ zN4YcNGG6VO;q*Q@qBTiA(c67K`D1L%Q@_$4-x&ss!%R9W0JYZJ?}VU<|l5T=$mLCU~3=+WMV`*DZg5mm2W=eP3tmS*D%w?2?~ zu)<_5Tb$kD8&Npi+MoZOYxis-XB_~Rt%W9WSXL$G$I0gov;}_`;j4cLhk@n-|aqM z_OnsR@Aw|dYWaZG(&=(|1Zs3Ze{X^v%{xn-gnh$#x=ZV`^nP#zdU1F5@s$B!6B)CQ zA`*vWg3kukS8e_FCk+Bc8u>!HGvJ z{Y8q7G|ygWsX{E5KbD8zipHABI~at zAoT&d2u$@;08ikT+8L|teb=&H--w|`8Ei%^n;teJa69OL8k@Yh9YmRWx+KNET*0jZ z9QUg%{YsStIb6zoD8PU24>va-_J5$ZT}hPG=ZtsyC?hlT*0FXUcXFbn!R6}8ojv-3s7P8YCy&NmY=En#wk& zmDwj5$Db&n_m!)g>(7y*c&aylYJvaWjmbTbL_jrh7|f~AJh49TLnbH+>(?Cob3SK% zJ$+&`(_RXi)Lm$;+L}BDSs4EjLq}0GB9~7gGR)v`z-~z#rc-3=-(Jiwl5@b+Jg0S2 zzRw=2sR(}fbLFf$|0qh{9f0Gh+d`QVCl=;=*CU?J1MPHi9vf;Od@=nSSn`hmEpP5h zjOQj|1unj9=J&WT)0TdD@^Dz^DzpZaRr2i?)(9NX;w?>u40hfEbeYbN{m6XYtJS~p z)HbmQAOe4iSX?v>Z`hS>EeP`cF|xquuL{Cx?zOh!P=g__sjcCrAAUjBsq@}#X7;E+ zAVip^#Z&nhSa`*E`>ET3(2c<);S&g7?;5~ER#41cilBg!k49+TG3fJCVvoFPWD3N( zuoSR~Aeu9ROaK@qfIhWToQ)m!8)lpL!}$VIE+x9;hXdT&!zVYWuLn!wHKqV0v zeB>Mq&_57S0IG-P!By8T*hin8{f5zI42Dv1ABgJ@gS&9=*Z0h)pI!Cx3;+D$V%qfP zy2l2cd`3~#V%SOogc1zKO)1Z4*~*E?b~>d9bRPI$kp#Q zA+I)V*=i%3C(FTpth0duI8V~GBQ6Dbr@i-eUyj3Xt~mLCq3wce6(C(J_Sj_#0uJ)c zkJ8Ex4K=5#iv=&+Nt#B=&6ed7T{56@SiK2{pH09{)PFW-e;0VWxZwE~X;WPCX5+1J zLcD0H6_RQ|SDOwcM*+>7rp5XY<(tfndy&HnIWFVp>xFLZDENRxxKN~Fm)a%V5O18!oAtv$ zTi>aZ1ePvIFqz!Rkqu??WnBZ75!!dY4mKwC8^3a3@xV=l?t!3gnZCcP)Oe|YgdLDH zl?2k}0evtvw#9P0%PQq7e_&FQ;6#Y91!gg9bMPnQ3(D-l?LQQowU>HnYn3qEd02X1 zZCCh;4xO-6y$28^Je$Bf#VmGcO$uZ?s~nE( zC_xI@#Vj+U(3$n>7&M><`el}ZYmI6_KT7}1knfG>X%6TIxcb(0v$u}jfNGi3dk(1_9$x#g_2jWL@Np1K)8Iyx Iy2HKy1<;aVlK=n! literal 0 HcmV?d00001 diff --git a/root/angular.C b/root/angular.C new file mode 100644 index 0000000..cd29cb3 --- /dev/null +++ b/root/angular.C @@ -0,0 +1,2019 @@ + void angular() +{ +//=========Macro generated from canvas: c_he7/title +//========= (Mon Oct 30 12:23:13 2023) by ROOT version6.10/08 + TCanvas *c_he7 = new TCanvas("c_he7", "title",67,55,1800,1000); + gStyle->SetOptStat(0); + c_he7->Range(0,0,1,1); + c_he7->SetFillColor(0); + c_he7->SetBorderMode(0); + c_he7->SetBorderSize(2); + c_he7->SetFrameBorderMode(0); + +// ------------>Primitives in pad: c_he7_1 + TPad *c_he7_1 = new TPad("c_he7_1", "c_he7_1",0.01,0.6766667,0.49,0.99); + c_he7_1->Draw(); + c_he7_1->cd(); + c_he7_1->Range(-4.483093,-3.336216,31.60923,30.37069); + c_he7_1->SetFillColor(0); + c_he7_1->SetBorderMode(0); + c_he7_1->SetBorderSize(2); + c_he7_1->SetFrameBorderMode(0); + c_he7_1->SetFrameBorderMode(0); + + TH2F *htemp__2 = new TH2F("htemp__2","ND_energy_CM4n:lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && ND_energy_CM4n>0 && lv_4n_mm.Mag()-4*939.565>-1 && ND_energy_CM4n<25}",40,-0.8738604,28,40,0.03447526,27); + htemp__2->SetDirectory(0); + htemp__2->SetStats(0); + + Int_t ci; // for color index setting + TColor *color; // for color definition with alpha + ci = TColor::GetColor("#000099"); + htemp__2->SetLineColor(ci); + htemp__2->SetMarkerStyle(20); + htemp__2->GetXaxis()->SetTitle("lv_4n_mm.Mag()-4*939.565"); + htemp__2->GetXaxis()->SetRange(1,40); + htemp__2->GetXaxis()->SetLabelFont(42); + htemp__2->GetXaxis()->SetLabelSize(0.035); + htemp__2->GetXaxis()->SetTitleSize(0.035); + htemp__2->GetXaxis()->SetTitleFont(42); + htemp__2->GetYaxis()->SetTitle("ND_energy_CM4n"); + htemp__2->GetYaxis()->SetRange(1,40); + htemp__2->GetYaxis()->SetLabelFont(42); + htemp__2->GetYaxis()->SetLabelSize(0.035); + htemp__2->GetYaxis()->SetTitleSize(0.035); + htemp__2->GetYaxis()->SetTitleOffset(0); + htemp__2->GetYaxis()->SetTitleFont(42); + htemp__2->GetZaxis()->SetLabelFont(42); + htemp__2->GetZaxis()->SetLabelSize(0.035); + htemp__2->GetZaxis()->SetTitleSize(0.035); + htemp__2->GetZaxis()->SetTitleFont(42); + htemp__2->Draw(""); + + TPaveText *pt = new TPaveText(0.15,0.9318933,0.85,0.995,"blNDC"); + pt->SetName("title"); + pt->SetBorderSize(0); + pt->SetFillColor(0); + pt->SetFillStyle(0); + pt->SetTextFont(42); + TText *pt_LaTex = pt->AddText("ND_energy_CM4n:lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && ND_energy_CM4n>0 && lv_4n_mm.Mag()-4*939.565>-1 && ND_energy_CM4n<25}"); + pt->Draw(); + + Double_t Graph0_fx3[136] = { + 2.293378, + 16.53346, + 2.775385, + 25.618, + 25.618, + 15.84522, + 23.25345, + 17.87187, + 11.87716, + 13.51793, + 9.709167, + 8.177823, + 15.57414, + 8.79528, + 15.49065, + 6.447484, + 5.468663, + 15.32576, + 15.44866, + 14.43255, + 6.329823, + 18.51257, + 4.628421, + 12.52081, + 19.88932, + 25.77312, + 16.85224, + 18.54831, + 14.66246, + 21.70794, + 6.069371, + 24.23959, + 15.92302, + 23.08188, + 13.04678, + 11.01585, + 9.533222, + 19.7312, + 2.250434, + 7.362024, + -0.8738604, + 14.52214, + 18.01523, + 25.71238, + 25.71238, + 19.00919, + 7.122373, + 14.77225, + 19.56621, + 11.88648, + 4.760872, + 25.25253, + 25.25253, + 13.78085, + 16.55374, + 16.32254, + 7.666889, + 17.17745, + 12.18927, + 12.18927, + 11.04283, + 15.78947, + 14.9445, + 12.80672, + 5.574539, + 1.294231, + 16.38851, + 12.18283, + 3.692083, + 20.43209, + 13.05403, + 3.35527, + 10.56624, + 11.00986, + 25.02393, + 21.30116, + 9.124693, + 18.07501, + 15.33917, + 19.30544, + 15.68121, + 14.6701, + 9.807625, + 13.23228, + 9.896872, + 8.678942, + 25.70097, + 17.79193, + 26.18868, + 22.8376, + 3.849799, + 18.20476, + 0.679571, + 14.87749, + 14.87749, + 2.729227, + 9.972701, + 18.46357, + 17.24962, + 16.66522, + 20.9598, + 18.17658, + 10.80134, + 4.07879, + 17.75641, + 3.633756, + 4.65884, + 17.17084, + 14.27111, + 23.64081, + 2.64591, + 10.31742, + 15.63272, + 13.63692, + 16.51991, + 15.26418, + 20.17847, + 9.149963, + 20.66901, + 24.80567, + 17.34251, + 0.2177003, + 3.090493, + 22.74068, + 22.66954, + 23.81961, + 14.55301, + 18.85983, + 17.55044, + 15.24921, + 7.34596, + 17.97639, + 18.21805, + 18.03173, + 20.16727, + 14.35985}; + Double_t Graph0_fy3[136] = { + 9.163973, + 6.449236, + 1.964816, + 1.359788, + 5.178768, + 0.03498398, + 8.769671, + 0.03447526, + 9.383176, + 2.666473, + 0.4635382, + 2.440106, + 5.390176, + 6.375852, + 1.925546, + 0.2982776, + 13.94083, + 4.734602, + 1.949477, + 0.9119229, + 17.27548, + 5.710012, + 24.2159, + 7.403358, + 3.93581, + 3.3992, + 9.162831, + 6.575353, + 4.866271, + 5.708108, + 7.070025, + 6.054096, + 1.499539, + 14.85669, + 0.9219248, + 2.782254, + 17.6368, + 5.868251, + 16.82547, + 7.0717, + 17.36037, + 5.493919, + 1.762917, + 8.676415, + 19.98228, + 11.87412, + 20.94647, + 9.686372, + 21.24376, + 14.36828, + 0.3806387, + 2.475904, + 8.730086, + 1.759002, + 21.01835, + 3.011658, + 3.843498, + 3.301521, + 3.249645, + 9.771441, + 8.041018, + 2.327877, + 19.21397, + 7.417117, + 1.852323, + 19.29377, + 2.539789, + 2.364044, + 15.25004, + 2.678995, + 0.9787855, + 7.585063, + 2.082057, + 0.2521605, + 3.357897, + 1.632345, + 3.1414, + 4.577259, + 4.179225, + 3.243581, + 8.190845, + 0.7892501, + 6.320226, + 9.764912, + 1.746821, + 2.202726, + 4.411285, + 5.531563, + 3.224871, + 11.51703, + 2.623029, + 4.777777, + 4.626465, + 1.353373, + 2.366684, + 4.382703, + 9.993218, + 7.656126, + 3.515664, + 8.187073, + 1.613118, + 2.693062, + 4.255403, + 1.414101, + 5.831221, + 7.573159, + 24.69945, + 0.8742778, + 1.442736, + 6.330581, + 4.266314, + 6.523858, + 7.44585, + 10.23886, + 10.57767, + 3.315528, + 5.560736, + 19.32783, + 2.242471, + 1.46184, + 2.058325, + 8.574683, + 1.994378, + 11.15542, + 7.585573, + 5.511144, + 3.669028, + 1.438863, + 6.547084, + 4.416309, + 17.42749, + 3.90872, + 6.114368, + 4.925547, + 7.415713, + 1.159745}; + TGraph *graph = new TGraph(136,Graph0_fx3,Graph0_fy3); + graph->SetName("Graph0"); + graph->SetTitle("Graph"); + + ci = TColor::GetColor("#000099"); + graph->SetLineColor(ci); + graph->SetMarkerStyle(20); + + TH1F *Graph_Graph3 = new TH1F("Graph_Graph3","Graph",136,-3.580114,28.89493); + Graph_Graph3->SetMinimum(0.03102773); + Graph_Graph3->SetMaximum(27.16594); + Graph_Graph3->SetDirectory(0); + Graph_Graph3->SetStats(0); + + ci = TColor::GetColor("#000099"); + Graph_Graph3->SetLineColor(ci); + Graph_Graph3->GetXaxis()->SetLabelFont(42); + Graph_Graph3->GetXaxis()->SetLabelSize(0.035); + Graph_Graph3->GetXaxis()->SetTitleSize(0.035); + Graph_Graph3->GetXaxis()->SetTitleFont(42); + Graph_Graph3->GetYaxis()->SetLabelFont(42); + Graph_Graph3->GetYaxis()->SetLabelSize(0.035); + Graph_Graph3->GetYaxis()->SetTitleSize(0.035); + Graph_Graph3->GetYaxis()->SetTitleOffset(0); + Graph_Graph3->GetYaxis()->SetTitleFont(42); + Graph_Graph3->GetZaxis()->SetLabelFont(42); + Graph_Graph3->GetZaxis()->SetLabelSize(0.035); + Graph_Graph3->GetZaxis()->SetTitleSize(0.035); + Graph_Graph3->GetZaxis()->SetTitleFont(42); + graph->SetHistogram(Graph_Graph3); + + graph->Draw("p"); + + Double_t Graph1_fx4[108] = { + 16.53346, + 2.775385, + 25.618, + 25.618, + 15.84522, + 23.25345, + 17.87187, + 11.87716, + 13.51793, + 9.709167, + 8.177823, + 15.57414, + 8.79528, + 15.49065, + 6.447484, + 15.32576, + 15.44866, + 14.43255, + 18.51257, + 12.52081, + 19.88932, + 25.77312, + 16.85224, + 18.54831, + 14.66246, + 21.70794, + 24.23959, + 15.92302, + 23.08188, + 13.04678, + 11.01585, + 19.7312, + 14.52214, + 18.01523, + 25.71238, + 19.00919, + 14.77225, + 4.760872, + 25.25253, + 25.25253, + 13.78085, + 16.32254, + 7.666889, + 17.17745, + 12.18927, + 11.04283, + 15.78947, + 12.80672, + 5.574539, + 16.38851, + 12.18283, + 20.43209, + 13.05403, + 10.56624, + 11.00986, + 25.02393, + 21.30116, + 9.124693, + 18.07501, + 15.33917, + 19.30544, + 15.68121, + 14.6701, + 9.807625, + 13.23228, + 9.896872, + 8.678942, + 25.70097, + 17.79193, + 26.18868, + 22.8376, + 3.849799, + 18.20476, + 14.87749, + 14.87749, + 18.46357, + 17.24962, + 16.66522, + 20.9598, + 18.17658, + 10.80134, + 4.07879, + 17.75641, + 17.17084, + 14.27111, + 23.64081, + 10.31742, + 15.63272, + 13.63692, + 16.51991, + 15.26418, + 20.17847, + 20.66901, + 24.80567, + 17.34251, + 3.090493, + 22.74068, + 22.66954, + 23.81961, + 14.55301, + 18.85983, + 17.55044, + 15.24921, + 17.97639, + 18.21805, + 18.03173, + 20.16727, + 14.35985}; + Double_t Graph1_fy4[108] = { + 6.449236, + 1.964816, + 1.359788, + 5.178768, + 0.03498398, + 8.769671, + 0.03447526, + 9.383176, + 2.666473, + 0.4635382, + 2.440106, + 5.390176, + 6.375852, + 1.925546, + 0.2982776, + 4.734602, + 1.949477, + 0.9119229, + 5.710012, + 7.403358, + 3.93581, + 3.3992, + 9.162831, + 6.575353, + 4.866271, + 5.708108, + 6.054096, + 1.499539, + 14.85669, + 0.9219248, + 2.782254, + 5.868251, + 5.493919, + 1.762917, + 8.676415, + 11.87412, + 9.686372, + 0.3806387, + 2.475904, + 8.730086, + 1.759002, + 3.011658, + 3.843498, + 3.301521, + 3.249645, + 8.041018, + 2.327877, + 7.417117, + 1.852323, + 2.539789, + 2.364044, + 2.678995, + 0.9787855, + 2.082057, + 0.2521605, + 3.357897, + 1.632345, + 3.1414, + 4.577259, + 4.179225, + 3.243581, + 8.190845, + 0.7892501, + 6.320226, + 9.764912, + 1.746821, + 2.202726, + 4.411285, + 5.531563, + 3.224871, + 11.51703, + 2.623029, + 4.777777, + 1.353373, + 2.366684, + 7.656126, + 3.515664, + 8.187073, + 1.613118, + 2.693062, + 4.255403, + 1.414101, + 5.831221, + 0.8742778, + 1.442736, + 6.330581, + 6.523858, + 7.44585, + 10.23886, + 10.57767, + 3.315528, + 5.560736, + 2.242471, + 1.46184, + 2.058325, + 1.994378, + 11.15542, + 7.585573, + 5.511144, + 3.669028, + 1.438863, + 6.547084, + 4.416309, + 3.90872, + 6.114368, + 4.925547, + 7.415713, + 1.159745}; + graph = new TGraph(108,Graph1_fx4,Graph1_fy4); + graph->SetName("Graph1"); + graph->SetTitle("Graph"); + + ci = TColor::GetColor("#000099"); + graph->SetLineColor(ci); + + ci = TColor::GetColor("#ff0000"); + graph->SetMarkerColor(ci); + graph->SetMarkerStyle(20); + + TH1F *Graph_Graph4 = new TH1F("Graph_Graph4","Graph",108,0.4340556,28.53001); + Graph_Graph4->SetMinimum(0.03102773); + Graph_Graph4->SetMaximum(16.33892); + Graph_Graph4->SetDirectory(0); + Graph_Graph4->SetStats(0); + + ci = TColor::GetColor("#000099"); + Graph_Graph4->SetLineColor(ci); + Graph_Graph4->GetXaxis()->SetLabelFont(42); + Graph_Graph4->GetXaxis()->SetLabelSize(0.035); + Graph_Graph4->GetXaxis()->SetTitleSize(0.035); + Graph_Graph4->GetXaxis()->SetTitleFont(42); + Graph_Graph4->GetYaxis()->SetLabelFont(42); + Graph_Graph4->GetYaxis()->SetLabelSize(0.035); + Graph_Graph4->GetYaxis()->SetTitleSize(0.035); + Graph_Graph4->GetYaxis()->SetTitleOffset(0); + Graph_Graph4->GetYaxis()->SetTitleFont(42); + Graph_Graph4->GetZaxis()->SetLabelFont(42); + Graph_Graph4->GetZaxis()->SetLabelSize(0.035); + Graph_Graph4->GetZaxis()->SetTitleSize(0.035); + Graph_Graph4->GetZaxis()->SetTitleFont(42); + graph->SetHistogram(Graph_Graph4); + + graph->Draw("p"); + + TF1 *fa0 = new TF1("fa","x*3/4",0,30); + fa0->SetFillColor(19); + fa0->SetFillStyle(0); + + ci = TColor::GetColor("#ff0000"); + fa0->SetLineColor(ci); + fa0->SetLineWidth(2); + fa0->GetXaxis()->SetLabelFont(42); + fa0->GetXaxis()->SetLabelSize(0.035); + fa0->GetXaxis()->SetTitleSize(0.035); + fa0->GetXaxis()->SetTitleFont(42); + fa0->GetYaxis()->SetLabelFont(42); + fa0->GetYaxis()->SetLabelSize(0.035); + fa0->GetYaxis()->SetTitleSize(0.035); + fa0->GetYaxis()->SetTitleOffset(0); + fa0->GetYaxis()->SetTitleFont(42); + fa0->Draw("same"); + + auto legend1 = new TLegend(0.1,0.9,0.9,1.); + legend1->SetHeader("All angles","C"); // option "C" allows to center the heade + legend1->SetTextColor(kRed); + legend1->SetTextSize(0.07); + legend1->Draw("same"); + + + c_he7_1->Modified(); + c_he7->cd(); + +// ------------>Primitives in pad: c_he7_2 + TPad *c_he7_2 = new TPad("c_he7_2", "c_he7_2",0.51,0.6766667,0.99,0.99); + c_he7_2->Draw(); + c_he7_2->cd(); + c_he7_2->Range(-3.75,-1.575,33.75,14.175); + c_he7_2->SetFillColor(0); + c_he7_2->SetBorderMode(0); + c_he7_2->SetBorderSize(2); + c_he7_2->SetFrameBorderMode(0); + c_he7_2->SetFrameBorderMode(0); + + TH1F *h1__3 = new TH1F("h1__3","lv_4n_mm.Mag()-4*939.565 {lv_4n_mm.Mag()-4*939.565>0 && lv_4n_mm.Mag()-4*939.565<30 && (lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0}",30,0,30); + h1__3->SetBinContent(1,2); + h1__3->SetBinContent(2,1); + h1__3->SetBinContent(3,5); + h1__3->SetBinContent(4,5); + h1__3->SetBinContent(5,5); + h1__3->SetBinContent(6,2); + h1__3->SetBinContent(7,3); + h1__3->SetBinContent(8,4); + h1__3->SetBinContent(9,3); + h1__3->SetBinContent(10,7); + h1__3->SetBinContent(11,3); + h1__3->SetBinContent(12,5); + h1__3->SetBinContent(13,5); + h1__3->SetBinContent(14,6); + h1__3->SetBinContent(15,11); + h1__3->SetBinContent(16,12); + h1__3->SetBinContent(17,7); + h1__3->SetBinContent(18,9); + h1__3->SetBinContent(19,10); + h1__3->SetBinContent(20,5); + h1__3->SetBinContent(21,5); + h1__3->SetBinContent(22,2); + h1__3->SetBinContent(23,3); + h1__3->SetBinContent(24,4); + h1__3->SetBinContent(25,2); + h1__3->SetBinContent(26,9); + h1__3->SetBinContent(27,1); + h1__3->SetEntries(136); + h1__3->SetStats(0); + + ci = TColor::GetColor("#ff0000"); + h1__3->SetMarkerColor(ci); + h1__3->SetMarkerStyle(20); + h1__3->GetXaxis()->SetLabelFont(42); + h1__3->GetXaxis()->SetLabelSize(0.035); + h1__3->GetXaxis()->SetTitleSize(0.035); + h1__3->GetXaxis()->SetTitleFont(42); + h1__3->GetYaxis()->SetLabelFont(42); + h1__3->GetYaxis()->SetLabelSize(0.035); + h1__3->GetYaxis()->SetTitleSize(0.035); + h1__3->GetYaxis()->SetTitleOffset(0); + h1__3->GetYaxis()->SetTitleFont(42); + h1__3->GetZaxis()->SetLabelFont(42); + h1__3->GetZaxis()->SetLabelSize(0.035); + h1__3->GetZaxis()->SetTitleSize(0.035); + h1__3->GetZaxis()->SetTitleFont(42); + h1__3->Draw(""); + + pt = new TPaveText(0.15,0.9318933,0.85,0.995,"blNDC"); + pt->SetName("title"); + pt->SetBorderSize(0); + pt->SetFillColor(0); + pt->SetFillStyle(0); + pt->SetTextFont(42); + pt_LaTex = pt->AddText("lv_4n_mm.Mag()-4*939.565 {lv_4n_mm.Mag()-4*939.565>0 && lv_4n_mm.Mag()-4*939.565<30 && (lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0}"); + pt->Draw(); + + TH1F *htemp__4 = new TH1F("htemp__4","lv_4n_mm.Mag()-4*939.565 {lv_4n_mm.Mag()-4*939.565>0 && lv_4n_mm.Mag()-4*939.565<30 && (lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && ND_energy_CM4n>0 && ( ND_energy_CM4n < 0.75*(lv_4n_mm.Mag()-4*939.565) + 0.5)}",30,0,30); + htemp__4->SetBinContent(3,1); + htemp__4->SetBinContent(4,2); + htemp__4->SetBinContent(5,2); + htemp__4->SetBinContent(6,1); + htemp__4->SetBinContent(7,1); + htemp__4->SetBinContent(8,1); + htemp__4->SetBinContent(9,3); + htemp__4->SetBinContent(10,4); + htemp__4->SetBinContent(11,3); + htemp__4->SetBinContent(12,4); + htemp__4->SetBinContent(13,4); + htemp__4->SetBinContent(14,6); + htemp__4->SetBinContent(15,10); + htemp__4->SetBinContent(16,12); + htemp__4->SetBinContent(17,6); + htemp__4->SetBinContent(18,9); + htemp__4->SetBinContent(19,10); + htemp__4->SetBinContent(20,4); + htemp__4->SetBinContent(21,5); + htemp__4->SetBinContent(22,2); + htemp__4->SetBinContent(23,3); + htemp__4->SetBinContent(24,4); + htemp__4->SetBinContent(25,2); + htemp__4->SetBinContent(26,8); + htemp__4->SetBinContent(27,1); + htemp__4->SetEntries(108); + htemp__4->SetDirectory(0); + htemp__4->SetStats(0); + + ci = TColor::GetColor("#ff0000"); + htemp__4->SetLineColor(ci); + + ci = TColor::GetColor("#ff0000"); + htemp__4->SetMarkerColor(ci); + htemp__4->SetMarkerStyle(20); + htemp__4->GetXaxis()->SetTitle("lv_4n_mm.Mag()-4*939.565"); + htemp__4->GetXaxis()->SetLabelFont(42); + htemp__4->GetXaxis()->SetLabelSize(0.035); + htemp__4->GetXaxis()->SetTitleSize(0.035); + htemp__4->GetXaxis()->SetTitleFont(42); + htemp__4->GetYaxis()->SetLabelFont(42); + htemp__4->GetYaxis()->SetLabelSize(0.035); + htemp__4->GetYaxis()->SetTitleSize(0.035); + htemp__4->GetYaxis()->SetTitleOffset(0); + htemp__4->GetYaxis()->SetTitleFont(42); + htemp__4->GetZaxis()->SetLabelFont(42); + htemp__4->GetZaxis()->SetLabelSize(0.035); + htemp__4->GetZaxis()->SetTitleSize(0.035); + htemp__4->GetZaxis()->SetTitleFont(42); + htemp__4->Draw("same"); + c_he7_2->Modified(); + c_he7->cd(); + +// ------------>Primitives in pad: c_he7_3 + TPad *c_he7_3 = new TPad("c_he7_3", "c_he7_3",0.01,0.3433333,0.49,0.6566667); + c_he7_3->Draw(); + c_he7_3->cd(); + c_he7_3->Range(-4.483093,-3.336216,31.60923,30.37069); + c_he7_3->SetFillColor(0); + c_he7_3->SetBorderMode(0); + c_he7_3->SetBorderSize(2); + c_he7_3->SetFrameBorderMode(0); + c_he7_3->SetFrameBorderMode(0); + + TH2F *htemp__5 = new TH2F("htemp__5","ND_energy_CM4n:lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && ND_energy_CM4n>0 && lv_4n_mm_CMR.Theta()*TMath::RadToDeg()<40 && lv_4n_mm.Mag()-4*939.565>-1 && ND_energy_CM4n<25}",40,-0.8738604,28,40,0.03447526,27); + htemp__5->SetDirectory(0); + htemp__5->SetStats(0); + + ci = TColor::GetColor("#ff0000"); + htemp__5->SetLineColor(ci); + htemp__5->SetMarkerStyle(20); + htemp__5->GetXaxis()->SetTitle("lv_4n_mm.Mag()-4*939.565"); + htemp__5->GetXaxis()->SetRange(1,40); + htemp__5->GetXaxis()->SetLabelFont(42); + htemp__5->GetXaxis()->SetLabelSize(0.035); + htemp__5->GetXaxis()->SetTitleSize(0.035); + htemp__5->GetXaxis()->SetTitleFont(42); + htemp__5->GetYaxis()->SetTitle("ND_energy_CM4n"); + htemp__5->GetYaxis()->SetRange(1,40); + htemp__5->GetYaxis()->SetLabelFont(42); + htemp__5->GetYaxis()->SetLabelSize(0.035); + htemp__5->GetYaxis()->SetTitleSize(0.035); + htemp__5->GetYaxis()->SetTitleOffset(0); + htemp__5->GetYaxis()->SetTitleFont(42); + htemp__5->GetZaxis()->SetLabelFont(42); + htemp__5->GetZaxis()->SetLabelSize(0.035); + htemp__5->GetZaxis()->SetTitleSize(0.035); + htemp__5->GetZaxis()->SetTitleFont(42); + htemp__5->Draw(""); + + pt = new TPaveText(0.15,0.9318933,0.85,0.995,"blNDC"); + pt->SetName("title"); + pt->SetBorderSize(0); + pt->SetFillColor(0); + pt->SetFillStyle(0); + pt->SetTextFont(42); + pt_LaTex = pt->AddText("ND_energy_CM4n:lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && ND_energy_CM4n>0 && lv_4n_mm_CMR.Theta()*TMath::RadToDeg()<40 && lv_4n_mm.Mag()-4*939.565>-1 && ND_energy_CM4n<25}"); + pt->Draw(); + + Double_t Graph0_fx5[124] = { + 2.293378, + 16.53346, + 2.775385, + 25.618, + 25.618, + 15.84522, + 23.25345, + 17.87187, + 11.87716, + 13.51793, + 9.709167, + 8.177823, + 15.57414, + 8.79528, + 15.49065, + 6.447484, + 5.468663, + 15.32576, + 15.44866, + 14.43255, + 6.329823, + 18.51257, + 4.628421, + 12.52081, + 19.88932, + 25.77312, + 16.85224, + 18.54831, + 14.66246, + 21.70794, + 6.069371, + 24.23959, + 15.92302, + 23.08188, + 13.04678, + 11.01585, + 7.362024, + -0.8738604, + 14.52214, + 18.01523, + 25.71238, + 25.71238, + 19.00919, + 7.122373, + 14.77225, + 19.56621, + 4.760872, + 25.25253, + 25.25253, + 13.78085, + 16.55374, + 16.32254, + 7.666889, + 12.18927, + 12.18927, + 11.04283, + 15.78947, + 14.9445, + 12.80672, + 5.574539, + 16.38851, + 12.18283, + 20.43209, + 13.05403, + 3.35527, + 10.56624, + 11.00986, + 25.02393, + 21.30116, + 9.124693, + 18.07501, + 15.33917, + 19.30544, + 15.68121, + 9.807625, + 13.23228, + 9.896872, + 8.678942, + 25.70097, + 17.79193, + 26.18868, + 22.8376, + 18.20476, + 0.679571, + 14.87749, + 14.87749, + 2.729227, + 9.972701, + 18.46357, + 16.66522, + 20.9598, + 18.17658, + 10.80134, + 4.07879, + 17.75641, + 3.633756, + 4.65884, + 17.17084, + 23.64081, + 2.64591, + 10.31742, + 15.63272, + 13.63692, + 15.26418, + 20.17847, + 9.149963, + 20.66901, + 24.80567, + 17.34251, + 0.2177003, + 3.090493, + 22.74068, + 22.66954, + 23.81961, + 14.55301, + 18.85983, + 17.55044, + 15.24921, + 7.34596, + 17.97639, + 18.21805, + 18.03173, + 20.16727, + 14.35985}; + Double_t Graph0_fy5[124] = { + 9.163973, + 6.449236, + 1.964816, + 1.359788, + 5.178768, + 0.03498398, + 8.769671, + 0.03447526, + 9.383176, + 2.666473, + 0.4635382, + 2.440106, + 5.390176, + 6.375852, + 1.925546, + 0.2982776, + 13.94083, + 4.734602, + 1.949477, + 0.9119229, + 17.27548, + 5.710012, + 24.2159, + 7.403358, + 3.93581, + 3.3992, + 9.162831, + 6.575353, + 4.866271, + 5.708108, + 7.070025, + 6.054096, + 1.499539, + 14.85669, + 0.9219248, + 2.782254, + 7.0717, + 17.36037, + 5.493919, + 1.762917, + 8.676415, + 19.98228, + 11.87412, + 20.94647, + 9.686372, + 21.24376, + 0.3806387, + 2.475904, + 8.730086, + 1.759002, + 21.01835, + 3.011658, + 3.843498, + 3.249645, + 9.771441, + 8.041018, + 2.327877, + 19.21397, + 7.417117, + 1.852323, + 2.539789, + 2.364044, + 2.678995, + 0.9787855, + 7.585063, + 2.082057, + 0.2521605, + 3.357897, + 1.632345, + 3.1414, + 4.577259, + 4.179225, + 3.243581, + 8.190845, + 6.320226, + 9.764912, + 1.746821, + 2.202726, + 4.411285, + 5.531563, + 3.224871, + 11.51703, + 4.777777, + 4.626465, + 1.353373, + 2.366684, + 4.382703, + 9.993218, + 7.656126, + 8.187073, + 1.613118, + 2.693062, + 4.255403, + 1.414101, + 5.831221, + 7.573159, + 24.69945, + 0.8742778, + 6.330581, + 4.266314, + 6.523858, + 7.44585, + 10.23886, + 3.315528, + 5.560736, + 19.32783, + 2.242471, + 1.46184, + 2.058325, + 8.574683, + 1.994378, + 11.15542, + 7.585573, + 5.511144, + 3.669028, + 1.438863, + 6.547084, + 4.416309, + 17.42749, + 3.90872, + 6.114368, + 4.925547, + 7.415713, + 1.159745}; + graph = new TGraph(124,Graph0_fx5,Graph0_fy5); + graph->SetName("Graph0"); + graph->SetTitle("Graph"); + + ci = TColor::GetColor("#ff0000"); + graph->SetLineColor(ci); + graph->SetMarkerStyle(20); + + TH1F *Graph_Graph5 = new TH1F("Graph_Graph5","Graph",124,-3.580114,28.89493); + Graph_Graph5->SetMinimum(0.03102773); + Graph_Graph5->SetMaximum(27.16594); + Graph_Graph5->SetDirectory(0); + Graph_Graph5->SetStats(0); + + ci = TColor::GetColor("#000099"); + Graph_Graph5->SetLineColor(ci); + Graph_Graph5->GetXaxis()->SetLabelFont(42); + Graph_Graph5->GetXaxis()->SetLabelSize(0.035); + Graph_Graph5->GetXaxis()->SetTitleSize(0.035); + Graph_Graph5->GetXaxis()->SetTitleFont(42); + Graph_Graph5->GetYaxis()->SetLabelFont(42); + Graph_Graph5->GetYaxis()->SetLabelSize(0.035); + Graph_Graph5->GetYaxis()->SetTitleSize(0.035); + Graph_Graph5->GetYaxis()->SetTitleOffset(0); + Graph_Graph5->GetYaxis()->SetTitleFont(42); + Graph_Graph5->GetZaxis()->SetLabelFont(42); + Graph_Graph5->GetZaxis()->SetLabelSize(0.035); + Graph_Graph5->GetZaxis()->SetTitleSize(0.035); + Graph_Graph5->GetZaxis()->SetTitleFont(42); + graph->SetHistogram(Graph_Graph5); + + graph->Draw("p"); + + Double_t Graph1_fx6[101] = { + 16.53346, + 2.775385, + 25.618, + 25.618, + 15.84522, + 23.25345, + 17.87187, + 11.87716, + 13.51793, + 9.709167, + 8.177823, + 15.57414, + 8.79528, + 15.49065, + 6.447484, + 15.32576, + 15.44866, + 14.43255, + 18.51257, + 12.52081, + 19.88932, + 25.77312, + 16.85224, + 18.54831, + 14.66246, + 21.70794, + 24.23959, + 15.92302, + 23.08188, + 13.04678, + 11.01585, + 14.52214, + 18.01523, + 25.71238, + 19.00919, + 14.77225, + 4.760872, + 25.25253, + 25.25253, + 13.78085, + 16.32254, + 7.666889, + 12.18927, + 11.04283, + 15.78947, + 12.80672, + 5.574539, + 16.38851, + 12.18283, + 20.43209, + 13.05403, + 10.56624, + 11.00986, + 25.02393, + 21.30116, + 9.124693, + 18.07501, + 15.33917, + 19.30544, + 15.68121, + 9.807625, + 13.23228, + 9.896872, + 8.678942, + 25.70097, + 17.79193, + 26.18868, + 22.8376, + 18.20476, + 14.87749, + 14.87749, + 18.46357, + 16.66522, + 20.9598, + 18.17658, + 10.80134, + 4.07879, + 17.75641, + 17.17084, + 23.64081, + 10.31742, + 15.63272, + 13.63692, + 15.26418, + 20.17847, + 20.66901, + 24.80567, + 17.34251, + 3.090493, + 22.74068, + 22.66954, + 23.81961, + 14.55301, + 18.85983, + 17.55044, + 15.24921, + 17.97639, + 18.21805, + 18.03173, + 20.16727, + 14.35985}; + Double_t Graph1_fy6[101] = { + 6.449236, + 1.964816, + 1.359788, + 5.178768, + 0.03498398, + 8.769671, + 0.03447526, + 9.383176, + 2.666473, + 0.4635382, + 2.440106, + 5.390176, + 6.375852, + 1.925546, + 0.2982776, + 4.734602, + 1.949477, + 0.9119229, + 5.710012, + 7.403358, + 3.93581, + 3.3992, + 9.162831, + 6.575353, + 4.866271, + 5.708108, + 6.054096, + 1.499539, + 14.85669, + 0.9219248, + 2.782254, + 5.493919, + 1.762917, + 8.676415, + 11.87412, + 9.686372, + 0.3806387, + 2.475904, + 8.730086, + 1.759002, + 3.011658, + 3.843498, + 3.249645, + 8.041018, + 2.327877, + 7.417117, + 1.852323, + 2.539789, + 2.364044, + 2.678995, + 0.9787855, + 2.082057, + 0.2521605, + 3.357897, + 1.632345, + 3.1414, + 4.577259, + 4.179225, + 3.243581, + 8.190845, + 6.320226, + 9.764912, + 1.746821, + 2.202726, + 4.411285, + 5.531563, + 3.224871, + 11.51703, + 4.777777, + 1.353373, + 2.366684, + 7.656126, + 8.187073, + 1.613118, + 2.693062, + 4.255403, + 1.414101, + 5.831221, + 0.8742778, + 6.330581, + 6.523858, + 7.44585, + 10.23886, + 3.315528, + 5.560736, + 2.242471, + 1.46184, + 2.058325, + 1.994378, + 11.15542, + 7.585573, + 5.511144, + 3.669028, + 1.438863, + 6.547084, + 4.416309, + 3.90872, + 6.114368, + 4.925547, + 7.415713, + 1.159745}; + graph = new TGraph(101,Graph1_fx6,Graph1_fy6); + graph->SetName("Graph1"); + graph->SetTitle("Graph"); + + ci = TColor::GetColor("#ff0000"); + graph->SetLineColor(ci); + + ci = TColor::GetColor("#ff0000"); + graph->SetMarkerColor(ci); + graph->SetMarkerStyle(20); + + TH1F *Graph_Graph6 = new TH1F("Graph_Graph6","Graph",101,0.4340556,28.53001); + Graph_Graph6->SetMinimum(0.03102773); + Graph_Graph6->SetMaximum(16.33892); + Graph_Graph6->SetDirectory(0); + Graph_Graph6->SetStats(0); + + ci = TColor::GetColor("#000099"); + Graph_Graph6->SetLineColor(ci); + Graph_Graph6->GetXaxis()->SetLabelFont(42); + Graph_Graph6->GetXaxis()->SetLabelSize(0.035); + Graph_Graph6->GetXaxis()->SetTitleSize(0.035); + Graph_Graph6->GetXaxis()->SetTitleFont(42); + Graph_Graph6->GetYaxis()->SetLabelFont(42); + Graph_Graph6->GetYaxis()->SetLabelSize(0.035); + Graph_Graph6->GetYaxis()->SetTitleSize(0.035); + Graph_Graph6->GetYaxis()->SetTitleOffset(0); + Graph_Graph6->GetYaxis()->SetTitleFont(42); + Graph_Graph6->GetZaxis()->SetLabelFont(42); + Graph_Graph6->GetZaxis()->SetLabelSize(0.035); + Graph_Graph6->GetZaxis()->SetTitleSize(0.035); + Graph_Graph6->GetZaxis()->SetTitleFont(42); + graph->SetHistogram(Graph_Graph6); + + graph->Draw("p"); + + TF1 *fa1 = new TF1("fa","x*3/4",0,30); + fa1->SetFillColor(19); + fa1->SetFillStyle(0); + + ci = TColor::GetColor("#ff0000"); + fa1->SetLineColor(ci); + fa1->SetLineWidth(2); + fa1->GetXaxis()->SetLabelFont(42); + fa1->GetXaxis()->SetLabelSize(0.035); + fa1->GetXaxis()->SetTitleSize(0.035); + fa1->GetXaxis()->SetTitleFont(42); + fa1->GetYaxis()->SetLabelFont(42); + fa1->GetYaxis()->SetLabelSize(0.035); + fa1->GetYaxis()->SetTitleSize(0.035); + fa1->GetYaxis()->SetTitleOffset(0); + fa1->GetYaxis()->SetTitleFont(42); + fa1->Draw("same"); + + auto legend2 = new TLegend(0.1,0.9,0.9,1.); + legend2->SetHeader("ThetaCM<40","C"); // option "C" allows to center the heade + legend2->SetTextColor(kRed); + legend2->SetTextSize(0.07); + legend2->Draw("same"); + + c_he7_3->Modified(); + c_he7->cd(); + +// ------------>Primitives in pad: c_he7_4 + TPad *c_he7_4 = new TPad("c_he7_4", "c_he7_4",0.51,0.3433333,0.99,0.6566667); + c_he7_4->Draw(); + c_he7_4->cd(); + c_he7_4->Range(-3.75,-1.575,33.75,14.175); + c_he7_4->SetFillColor(0); + c_he7_4->SetBorderMode(0); + c_he7_4->SetBorderSize(2); + c_he7_4->SetFrameBorderMode(0); + c_he7_4->SetFrameBorderMode(0); + + TH1F *h2__6 = new TH1F("h2__6","lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && lv_4n_mm_CMR.Theta()*TMath::RadToDeg()<40}",30,0,30); + h2__6->SetBinContent(0,108); + h2__6->SetBinContent(1,2); + h2__6->SetBinContent(3,4); + h2__6->SetBinContent(4,3); + h2__6->SetBinContent(5,5); + h2__6->SetBinContent(6,2); + h2__6->SetBinContent(7,3); + h2__6->SetBinContent(8,4); + h2__6->SetBinContent(9,3); + h2__6->SetBinContent(10,6); + h2__6->SetBinContent(11,3); + h2__6->SetBinContent(12,4); + h2__6->SetBinContent(13,5); + h2__6->SetBinContent(14,6); + h2__6->SetBinContent(15,9); + h2__6->SetBinContent(16,12); + h2__6->SetBinContent(17,6); + h2__6->SetBinContent(18,7); + h2__6->SetBinContent(19,10); + h2__6->SetBinContent(20,4); + h2__6->SetBinContent(21,5); + h2__6->SetBinContent(22,2); + h2__6->SetBinContent(23,3); + h2__6->SetBinContent(24,4); + h2__6->SetBinContent(25,2); + h2__6->SetBinContent(26,9); + h2__6->SetBinContent(27,1); + h2__6->SetEntries(232); + h2__6->SetStats(0); + + ci = TColor::GetColor("#ff0000"); + h2__6->SetMarkerColor(ci); + h2__6->SetMarkerStyle(20); + h2__6->GetXaxis()->SetLabelFont(42); + h2__6->GetXaxis()->SetLabelSize(0.035); + h2__6->GetXaxis()->SetTitleSize(0.035); + h2__6->GetXaxis()->SetTitleFont(42); + h2__6->GetYaxis()->SetLabelFont(42); + h2__6->GetYaxis()->SetLabelSize(0.035); + h2__6->GetYaxis()->SetTitleSize(0.035); + h2__6->GetYaxis()->SetTitleOffset(0); + h2__6->GetYaxis()->SetTitleFont(42); + h2__6->GetZaxis()->SetLabelFont(42); + h2__6->GetZaxis()->SetLabelSize(0.035); + h2__6->GetZaxis()->SetTitleSize(0.035); + h2__6->GetZaxis()->SetTitleFont(42); + h2__6->Draw(""); + + pt = new TPaveText(0.15,0.9318933,0.85,0.995,"blNDC"); + pt->SetName("title"); + pt->SetBorderSize(0); + pt->SetFillColor(0); + pt->SetFillStyle(0); + pt->SetTextFont(42); + pt_LaTex = pt->AddText("lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && lv_4n_mm_CMR.Theta()*TMath::RadToDeg()<40}"); + pt->Draw(); + + TH1F *htemp__7 = new TH1F("htemp__7","lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && lv_4n_mm_CMR.Theta()*TMath::RadToDeg()<40 && ND_energy_CM4n>0 && ( ND_energy_CM4n < 0.75*(lv_4n_mm.Mag()-4*939.565) + 0.5)}",30,0,30); + htemp__7->SetBinContent(3,1); + htemp__7->SetBinContent(4,1); + htemp__7->SetBinContent(5,2); + htemp__7->SetBinContent(6,1); + htemp__7->SetBinContent(7,1); + htemp__7->SetBinContent(8,1); + htemp__7->SetBinContent(9,3); + htemp__7->SetBinContent(10,4); + htemp__7->SetBinContent(11,3); + htemp__7->SetBinContent(12,4); + htemp__7->SetBinContent(13,4); + htemp__7->SetBinContent(14,6); + htemp__7->SetBinContent(15,8); + htemp__7->SetBinContent(16,12); + htemp__7->SetBinContent(17,5); + htemp__7->SetBinContent(18,7); + htemp__7->SetBinContent(19,10); + htemp__7->SetBinContent(20,3); + htemp__7->SetBinContent(21,5); + htemp__7->SetBinContent(22,2); + htemp__7->SetBinContent(23,3); + htemp__7->SetBinContent(24,4); + htemp__7->SetBinContent(25,2); + htemp__7->SetBinContent(26,8); + htemp__7->SetBinContent(27,1); + htemp__7->SetEntries(101); + htemp__7->SetDirectory(0); + htemp__7->SetStats(0); + + ci = TColor::GetColor("#ff0000"); + htemp__7->SetLineColor(ci); + + ci = TColor::GetColor("#ff0000"); + htemp__7->SetMarkerColor(ci); + htemp__7->SetMarkerStyle(20); + htemp__7->GetXaxis()->SetTitle("lv_4n_mm.Mag()-4*939.565"); + htemp__7->GetXaxis()->SetLabelFont(42); + htemp__7->GetXaxis()->SetLabelSize(0.035); + htemp__7->GetXaxis()->SetTitleSize(0.035); + htemp__7->GetXaxis()->SetTitleFont(42); + htemp__7->GetYaxis()->SetLabelFont(42); + htemp__7->GetYaxis()->SetLabelSize(0.035); + htemp__7->GetYaxis()->SetTitleSize(0.035); + htemp__7->GetYaxis()->SetTitleOffset(0); + htemp__7->GetYaxis()->SetTitleFont(42); + htemp__7->GetZaxis()->SetLabelFont(42); + htemp__7->GetZaxis()->SetLabelSize(0.035); + htemp__7->GetZaxis()->SetTitleSize(0.035); + htemp__7->GetZaxis()->SetTitleFont(42); + htemp__7->Draw("same"); + c_he7_4->Modified(); + c_he7->cd(); + +// ------------>Primitives in pad: c_he7_5 + TPad *c_he7_5 = new TPad("c_he7_5", "c_he7_5",0.01,0.01,0.49,0.3233333); + c_he7_5->Draw(); + c_he7_5->cd(); + c_he7_5->Range(-4.483093,-3.336216,31.60923,30.37069); + c_he7_5->SetFillColor(0); + c_he7_5->SetBorderMode(0); + c_he7_5->SetBorderSize(2); + c_he7_5->SetFrameBorderMode(0); + c_he7_5->SetFrameBorderMode(0); + + TH2F *htemp__8 = new TH2F("htemp__8","ND_energy_CM4n:lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && ND_energy_CM4n>0 && lv_4n_mm_CMR.Theta()*TMath::RadToDeg()<30 && lv_4n_mm.Mag()-4*939.565>-1 && ND_energy_CM4n<25}",40,-0.8738604,28,40,0.03447526,27); + htemp__8->SetDirectory(0); + htemp__8->SetStats(0); + + ci = TColor::GetColor("#ff0000"); + htemp__8->SetLineColor(ci); + htemp__8->SetMarkerStyle(20); + htemp__8->GetXaxis()->SetTitle("lv_4n_mm.Mag()-4*939.565"); + htemp__8->GetXaxis()->SetRange(1,40); + htemp__8->GetXaxis()->SetLabelFont(42); + htemp__8->GetXaxis()->SetLabelSize(0.035); + htemp__8->GetXaxis()->SetTitleSize(0.035); + htemp__8->GetXaxis()->SetTitleFont(42); + htemp__8->GetYaxis()->SetTitle("ND_energy_CM4n"); + htemp__8->GetYaxis()->SetRange(1,40); + htemp__8->GetYaxis()->SetLabelFont(42); + htemp__8->GetYaxis()->SetLabelSize(0.035); + htemp__8->GetYaxis()->SetTitleSize(0.035); + htemp__8->GetYaxis()->SetTitleOffset(0); + htemp__8->GetYaxis()->SetTitleFont(42); + htemp__8->GetZaxis()->SetLabelFont(42); + htemp__8->GetZaxis()->SetLabelSize(0.035); + htemp__8->GetZaxis()->SetTitleSize(0.035); + htemp__8->GetZaxis()->SetTitleFont(42); + htemp__8->Draw(""); + + pt = new TPaveText(0.15,0.9318933,0.85,0.995,"blNDC"); + pt->SetName("title"); + pt->SetBorderSize(0); + pt->SetFillColor(0); + pt->SetFillStyle(0); + pt->SetTextFont(42); + pt_LaTex = pt->AddText("ND_energy_CM4n:lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && ND_energy_CM4n>0 && lv_4n_mm_CMR.Theta()*TMath::RadToDeg()<30 && lv_4n_mm.Mag()-4*939.565>-1 && ND_energy_CM4n<25}"); + pt->Draw(); + + Double_t Graph0_fx7[72] = { + 2.293378, + 2.775385, + 25.618, + 25.618, + 15.84522, + 17.87187, + 13.51793, + 9.709167, + 8.177823, + 15.57414, + 8.79528, + 15.49065, + 6.447484, + 5.468663, + 15.44866, + 6.329823, + 18.51257, + 4.628421, + 12.52081, + 25.77312, + 14.66246, + 21.70794, + 15.92302, + 23.08188, + 11.01585, + -0.8738604, + 18.01523, + 19.00919, + 14.77225, + 19.56621, + 25.25253, + 25.25253, + 13.78085, + 16.32254, + 7.666889, + 15.78947, + 14.9445, + 5.574539, + 16.38851, + 20.43209, + 13.05403, + 11.00986, + 25.02393, + 18.07501, + 9.807625, + 8.678942, + 25.70097, + 26.18868, + 0.679571, + 14.87749, + 14.87749, + 2.729227, + 9.972701, + 16.66522, + 18.17658, + 4.07879, + 17.75641, + 3.633756, + 4.65884, + 2.64591, + 15.26418, + 24.80567, + 0.2177003, + 3.090493, + 23.81961, + 18.85983, + 17.55044, + 15.24921, + 7.34596, + 18.03173, + 20.16727, + 14.35985}; + Double_t Graph0_fy7[72] = { + 9.163973, + 1.964816, + 1.359788, + 5.178768, + 0.03498398, + 0.03447526, + 2.666473, + 0.4635382, + 2.440106, + 5.390176, + 6.375852, + 1.925546, + 0.2982776, + 13.94083, + 1.949477, + 17.27548, + 5.710012, + 24.2159, + 7.403358, + 3.3992, + 4.866271, + 5.708108, + 1.499539, + 14.85669, + 2.782254, + 17.36037, + 1.762917, + 11.87412, + 9.686372, + 21.24376, + 2.475904, + 8.730086, + 1.759002, + 3.011658, + 3.843498, + 2.327877, + 19.21397, + 1.852323, + 2.539789, + 2.678995, + 0.9787855, + 0.2521605, + 3.357897, + 4.577259, + 6.320226, + 2.202726, + 4.411285, + 3.224871, + 4.626465, + 1.353373, + 2.366684, + 4.382703, + 9.993218, + 8.187073, + 2.693062, + 1.414101, + 5.831221, + 7.573159, + 24.69945, + 4.266314, + 3.315528, + 1.46184, + 8.574683, + 1.994378, + 5.511144, + 1.438863, + 6.547084, + 4.416309, + 17.42749, + 4.925547, + 7.415713, + 1.159745}; + graph = new TGraph(72,Graph0_fx7,Graph0_fy7); + graph->SetName("Graph0"); + graph->SetTitle("Graph"); + + ci = TColor::GetColor("#ff0000"); + graph->SetLineColor(ci); + graph->SetMarkerStyle(20); + + TH1F *Graph_Graph7 = new TH1F("Graph_Graph7","Graph",100,-3.580114,28.89493); + Graph_Graph7->SetMinimum(0.03102773); + Graph_Graph7->SetMaximum(27.16594); + Graph_Graph7->SetDirectory(0); + Graph_Graph7->SetStats(0); + + ci = TColor::GetColor("#000099"); + Graph_Graph7->SetLineColor(ci); + Graph_Graph7->GetXaxis()->SetLabelFont(42); + Graph_Graph7->GetXaxis()->SetLabelSize(0.035); + Graph_Graph7->GetXaxis()->SetTitleSize(0.035); + Graph_Graph7->GetXaxis()->SetTitleFont(42); + Graph_Graph7->GetYaxis()->SetLabelFont(42); + Graph_Graph7->GetYaxis()->SetLabelSize(0.035); + Graph_Graph7->GetYaxis()->SetTitleSize(0.035); + Graph_Graph7->GetYaxis()->SetTitleOffset(0); + Graph_Graph7->GetYaxis()->SetTitleFont(42); + Graph_Graph7->GetZaxis()->SetLabelFont(42); + Graph_Graph7->GetZaxis()->SetLabelSize(0.035); + Graph_Graph7->GetZaxis()->SetTitleSize(0.035); + Graph_Graph7->GetZaxis()->SetTitleFont(42); + graph->SetHistogram(Graph_Graph7); + + graph->Draw("p"); + + Double_t Graph1_fx8[57] = { + 2.775385, + 25.618, + 25.618, + 15.84522, + 17.87187, + 13.51793, + 9.709167, + 8.177823, + 15.57414, + 8.79528, + 15.49065, + 6.447484, + 15.44866, + 18.51257, + 12.52081, + 25.77312, + 14.66246, + 21.70794, + 15.92302, + 23.08188, + 11.01585, + 18.01523, + 19.00919, + 14.77225, + 25.25253, + 25.25253, + 13.78085, + 16.32254, + 7.666889, + 15.78947, + 5.574539, + 16.38851, + 20.43209, + 13.05403, + 11.00986, + 25.02393, + 18.07501, + 9.807625, + 8.678942, + 25.70097, + 26.18868, + 14.87749, + 14.87749, + 16.66522, + 18.17658, + 4.07879, + 17.75641, + 15.26418, + 24.80567, + 3.090493, + 23.81961, + 18.85983, + 17.55044, + 15.24921, + 18.03173, + 20.16727, + 14.35985}; + Double_t Graph1_fy8[57] = { + 1.964816, + 1.359788, + 5.178768, + 0.03498398, + 0.03447526, + 2.666473, + 0.4635382, + 2.440106, + 5.390176, + 6.375852, + 1.925546, + 0.2982776, + 1.949477, + 5.710012, + 7.403358, + 3.3992, + 4.866271, + 5.708108, + 1.499539, + 14.85669, + 2.782254, + 1.762917, + 11.87412, + 9.686372, + 2.475904, + 8.730086, + 1.759002, + 3.011658, + 3.843498, + 2.327877, + 1.852323, + 2.539789, + 2.678995, + 0.9787855, + 0.2521605, + 3.357897, + 4.577259, + 6.320226, + 2.202726, + 4.411285, + 3.224871, + 1.353373, + 2.366684, + 8.187073, + 2.693062, + 1.414101, + 5.831221, + 3.315528, + 1.46184, + 1.994378, + 5.511144, + 1.438863, + 6.547084, + 4.416309, + 4.925547, + 7.415713, + 1.159745}; + graph = new TGraph(57,Graph1_fx8,Graph1_fy8); + graph->SetName("Graph1"); + graph->SetTitle("Graph"); + + ci = TColor::GetColor("#ff0000"); + graph->SetLineColor(ci); + + ci = TColor::GetColor("#ff0000"); + graph->SetMarkerColor(ci); + graph->SetMarkerStyle(20); + + TH1F *Graph_Graph8 = new TH1F("Graph_Graph8","Graph",100,0.4340556,28.53001); + Graph_Graph8->SetMinimum(0.03102773); + Graph_Graph8->SetMaximum(16.33892); + Graph_Graph8->SetDirectory(0); + Graph_Graph8->SetStats(0); + + ci = TColor::GetColor("#000099"); + Graph_Graph8->SetLineColor(ci); + Graph_Graph8->GetXaxis()->SetLabelFont(42); + Graph_Graph8->GetXaxis()->SetLabelSize(0.035); + Graph_Graph8->GetXaxis()->SetTitleSize(0.035); + Graph_Graph8->GetXaxis()->SetTitleFont(42); + Graph_Graph8->GetYaxis()->SetLabelFont(42); + Graph_Graph8->GetYaxis()->SetLabelSize(0.035); + Graph_Graph8->GetYaxis()->SetTitleSize(0.035); + Graph_Graph8->GetYaxis()->SetTitleOffset(0); + Graph_Graph8->GetYaxis()->SetTitleFont(42); + Graph_Graph8->GetZaxis()->SetLabelFont(42); + Graph_Graph8->GetZaxis()->SetLabelSize(0.035); + Graph_Graph8->GetZaxis()->SetTitleSize(0.035); + Graph_Graph8->GetZaxis()->SetTitleFont(42); + graph->SetHistogram(Graph_Graph8); + + graph->Draw("p"); + + TF1 *fa2 = new TF1("fa","x*3/4",0,30); + fa2->SetFillColor(19); + fa2->SetFillStyle(0); + + ci = TColor::GetColor("#ff0000"); + fa2->SetLineColor(ci); + fa2->SetLineWidth(2); + fa2->GetXaxis()->SetLabelFont(42); + fa2->GetXaxis()->SetLabelSize(0.035); + fa2->GetXaxis()->SetTitleSize(0.035); + fa2->GetXaxis()->SetTitleFont(42); + fa2->GetYaxis()->SetLabelFont(42); + fa2->GetYaxis()->SetLabelSize(0.035); + fa2->GetYaxis()->SetTitleSize(0.035); + fa2->GetYaxis()->SetTitleOffset(0); + fa2->GetYaxis()->SetTitleFont(42); + fa2->Draw("same"); + + auto legend3 = new TLegend(0.1,0.9,0.9,1.); + legend3->SetHeader("ThetaCM<30","C"); // option "C" allows to center the heade + legend3->SetTextColor(kRed); + legend3->SetTextSize(0.07); + legend3->Draw("same"); + + c_he7_5->Modified(); + c_he7->cd(); + + + +// ------------>Primitives in pad: c_he7_6 + TPad *c_he7_6 = new TPad("c_he7_6", "c_he7_6",0.51,0.01,0.99,0.3233333); + c_he7_6->Draw(); + c_he7_6->cd(); + c_he7_6->Range(-3.75,-1.05,33.75,9.45); + c_he7_6->SetFillColor(0); + c_he7_6->SetBorderMode(0); + c_he7_6->SetBorderSize(2); + c_he7_6->SetFrameBorderMode(0); + c_he7_6->SetFrameBorderMode(0); + + TH1F *h3__9 = new TH1F("h3__9","lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && lv_4n_mm_CMR.Theta()*TMath::RadToDeg()<30}",30,0,30); + h3__9->SetBinContent(0,99); + h3__9->SetBinContent(1,2); + h3__9->SetBinContent(3,4); + h3__9->SetBinContent(4,2); + h3__9->SetBinContent(5,4); + h3__9->SetBinContent(6,2); + h3__9->SetBinContent(7,2); + h3__9->SetBinContent(8,2); + h3__9->SetBinContent(9,3); + h3__9->SetBinContent(10,3); + h3__9->SetBinContent(12,2); + h3__9->SetBinContent(13,1); + h3__9->SetBinContent(14,3); + h3__9->SetBinContent(15,6); + h3__9->SetBinContent(16,8); + h3__9->SetBinContent(17,3); + h3__9->SetBinContent(18,3); + h3__9->SetBinContent(19,6); + h3__9->SetBinContent(20,2); + h3__9->SetBinContent(21,2); + h3__9->SetBinContent(22,1); + h3__9->SetBinContent(24,2); + h3__9->SetBinContent(25,1); + h3__9->SetBinContent(26,7); + h3__9->SetBinContent(27,1); + h3__9->SetEntries(171); + h3__9->SetStats(0); + + ci = TColor::GetColor("#ff0000"); + h3__9->SetMarkerColor(ci); + h3__9->SetMarkerStyle(20); + h3__9->GetXaxis()->SetLabelFont(42); + h3__9->GetXaxis()->SetLabelSize(0.035); + h3__9->GetXaxis()->SetTitleSize(0.035); + h3__9->GetXaxis()->SetTitleFont(42); + h3__9->GetYaxis()->SetLabelFont(42); + h3__9->GetYaxis()->SetLabelSize(0.035); + h3__9->GetYaxis()->SetTitleSize(0.035); + h3__9->GetYaxis()->SetTitleOffset(0); + h3__9->GetYaxis()->SetTitleFont(42); + h3__9->GetZaxis()->SetLabelFont(42); + h3__9->GetZaxis()->SetLabelSize(0.035); + h3__9->GetZaxis()->SetTitleSize(0.035); + h3__9->GetZaxis()->SetTitleFont(42); + h3__9->Draw(""); + + pt = new TPaveText(0.15,0.9318933,0.85,0.995,"blNDC"); + pt->SetName("title"); + pt->SetBorderSize(0); + pt->SetFillColor(0); + pt->SetFillStyle(0); + pt->SetTextFont(42); + pt_LaTex = pt->AddText("lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && lv_4n_mm_CMR.Theta()*TMath::RadToDeg()<30}"); + pt->Draw(); + + TH1F *htemp__10 = new TH1F("htemp__10","lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && lv_4n_mm_CMR.Theta()*TMath::RadToDeg()<30 && ND_energy_CM4n>0 && ( ND_energy_CM4n < 0.75*(lv_4n_mm.Mag()-4*939.565) + 0.5)}",30,0,30); + htemp__10->SetBinContent(3,1); + htemp__10->SetBinContent(4,1); + htemp__10->SetBinContent(5,1); + htemp__10->SetBinContent(6,1); + htemp__10->SetBinContent(7,1); + htemp__10->SetBinContent(8,1); + htemp__10->SetBinContent(9,3); + htemp__10->SetBinContent(10,2); + htemp__10->SetBinContent(12,2); + htemp__10->SetBinContent(13,1); + htemp__10->SetBinContent(14,3); + htemp__10->SetBinContent(15,5); + htemp__10->SetBinContent(16,8); + htemp__10->SetBinContent(17,3); + htemp__10->SetBinContent(18,3); + htemp__10->SetBinContent(19,6); + htemp__10->SetBinContent(20,1); + htemp__10->SetBinContent(21,2); + htemp__10->SetBinContent(22,1); + htemp__10->SetBinContent(24,2); + htemp__10->SetBinContent(25,1); + htemp__10->SetBinContent(26,7); + htemp__10->SetBinContent(27,1); + htemp__10->SetEntries(57); + htemp__10->SetDirectory(0); + htemp__10->SetStats(0); + + ci = TColor::GetColor("#ff0000"); + htemp__10->SetLineColor(ci); + + ci = TColor::GetColor("#ff0000"); + htemp__10->SetMarkerColor(ci); + htemp__10->SetMarkerStyle(20); + htemp__10->GetXaxis()->SetTitle("lv_4n_mm.Mag()-4*939.565"); + htemp__10->GetXaxis()->SetLabelFont(42); + htemp__10->GetXaxis()->SetLabelSize(0.035); + htemp__10->GetXaxis()->SetTitleSize(0.035); + htemp__10->GetXaxis()->SetTitleFont(42); + htemp__10->GetYaxis()->SetLabelFont(42); + htemp__10->GetYaxis()->SetLabelSize(0.035); + htemp__10->GetYaxis()->SetTitleSize(0.035); + htemp__10->GetYaxis()->SetTitleOffset(0); + htemp__10->GetYaxis()->SetTitleFont(42); + htemp__10->GetZaxis()->SetLabelFont(42); + htemp__10->GetZaxis()->SetLabelSize(0.035); + htemp__10->GetZaxis()->SetTitleSize(0.035); + htemp__10->GetZaxis()->SetTitleFont(42); + htemp__10->Draw("same"); + c_he7_6->Modified(); + c_he7->cd(); + c_he7->Modified(); + c_he7->cd(); + c_he7->SetSelected(c_he7); +} diff --git a/root/mm.C b/root/mm.C index b772c84..ad7e7a9 100644 --- a/root/mm.C +++ b/root/mm.C @@ -1,7 +1,7 @@ void mm() { //=========Macro generated from canvas: c_angn4/title -//========= (Sun Oct 29 18:31:05 2023) by ROOT version6.10/08 +//========= (Mon Oct 30 11:37:23 2023) by ROOT version6.10/08 TCanvas *c_angn4 = new TCanvas("c_angn4", "title",117,57,1000,1000); gStyle->SetOptStat(0); c_angn4->Range(-4.558442,-1.575,30.50649,14.175); @@ -274,7 +274,7 @@ void mm() Graph_Graph1->GetZaxis()->SetTitleFont(42); graph->SetHistogram(Graph_Graph1); - graph->Draw(""); + graph->Draw("same"); c_angn4->Modified(); c_angn4->cd(); c_angn4->SetSelected(c_angn4); diff --git a/root/single_angular.C b/root/single_angular.C new file mode 100644 index 0000000..3cd1c2b --- /dev/null +++ b/root/single_angular.C @@ -0,0 +1,616 @@ +void single_angular() +{ +//=========Macro generated from canvas: c_angn4/title +//========= (Mon Oct 30 12:28:11 2023) by ROOT version6.10/08 + TCanvas *c_angn4 = new TCanvas("c_angn4", "title",67,55,1800,1000); + gStyle->SetOptStat(0); + c_angn4->Range(-6.875,-5.000001,31.875,95); + c_angn4->SetFillColor(0); + c_angn4->SetBorderMode(0); + c_angn4->SetBorderSize(2); + c_angn4->SetFrameBorderMode(0); + c_angn4->SetFrameBorderMode(0); + + TH2F *htemp__1 = new TH2F("htemp__1","lv_4n_mm_CMR.Theta()*TMath::RadToDeg():lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && ND_energy_CM4n>0 && lv_4n_mm.Mag()-4*939.565>-1}",40,-3,28,40,5,85); + htemp__1->SetDirectory(0); + htemp__1->SetStats(0); + + Int_t ci; // for color index setting + TColor *color; // for color definition with alpha + ci = TColor::GetColor("#000099"); + htemp__1->SetLineColor(ci); + htemp__1->SetMarkerStyle(8); + htemp__1->GetXaxis()->SetTitle("lv_4n_mm.Mag()-4*939.565"); + htemp__1->GetXaxis()->SetRange(1,40); + htemp__1->GetXaxis()->SetLabelFont(42); + htemp__1->GetXaxis()->SetLabelSize(0.035); + htemp__1->GetXaxis()->SetTitleSize(0.035); + htemp__1->GetXaxis()->SetTitleFont(42); + htemp__1->GetYaxis()->SetTitle("lv_4n_mm_CMR.Theta()*TMath::RadToDeg()"); + htemp__1->GetYaxis()->SetRange(1,40); + htemp__1->GetYaxis()->SetLabelFont(42); + htemp__1->GetYaxis()->SetLabelSize(0.035); + htemp__1->GetYaxis()->SetTitleSize(0.035); + htemp__1->GetYaxis()->SetTitleOffset(0); + htemp__1->GetYaxis()->SetTitleFont(42); + htemp__1->GetZaxis()->SetLabelFont(42); + htemp__1->GetZaxis()->SetLabelSize(0.035); + htemp__1->GetZaxis()->SetTitleSize(0.035); + htemp__1->GetZaxis()->SetTitleFont(42); + htemp__1->Draw(""); + + TPaveText *pt = new TPaveText(0.15,0.9344439,0.85,0.995,"blNDC"); + pt->SetName("title"); + pt->SetBorderSize(0); + pt->SetFillColor(0); + pt->SetFillStyle(0); + pt->SetTextFont(42); + TText *pt_LaTex = pt->AddText("lv_4n_mm_CMR.Theta()*TMath::RadToDeg():lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && multv_1==0 && multv_2==0 && multv_3==0 && multv_4==0 && ND_energy_CM4n>0 && lv_4n_mm.Mag()-4*939.565>-1}"); + pt->Draw(); + + Double_t Graph0_fx1[137] = { + 2.293378, + 16.53346, + 2.775385, + 25.618, + 25.618, + 15.84522, + 23.25345, + 17.87187, + 11.87716, + 13.51793, + 9.709167, + 8.177823, + 15.57414, + 8.79528, + 15.49065, + 6.447484, + 5.468663, + 15.32576, + 15.44866, + 14.43255, + 6.329823, + 18.51257, + 4.628421, + 12.52081, + 19.88932, + 25.77312, + 16.85224, + 18.54831, + 4.118373, + 14.66246, + 21.70794, + 6.069371, + 24.23959, + 15.92302, + 23.08188, + 13.04678, + 11.01585, + 9.533222, + 19.7312, + 2.250434, + 7.362024, + -0.8738604, + 14.52214, + 18.01523, + 25.71238, + 25.71238, + 19.00919, + 7.122373, + 14.77225, + 19.56621, + 11.88648, + 4.760872, + 25.25253, + 25.25253, + 13.78085, + 16.55374, + 16.32254, + 7.666889, + 17.17745, + 12.18927, + 12.18927, + 11.04283, + 15.78947, + 14.9445, + 12.80672, + 5.574539, + 1.294231, + 16.38851, + 12.18283, + 3.692083, + 20.43209, + 13.05403, + 3.35527, + 10.56624, + 11.00986, + 25.02393, + 21.30116, + 9.124693, + 18.07501, + 15.33917, + 19.30544, + 15.68121, + 14.6701, + 9.807625, + 13.23228, + 9.896872, + 8.678942, + 25.70097, + 17.79193, + 26.18868, + 22.8376, + 3.849799, + 18.20476, + 0.679571, + 14.87749, + 14.87749, + 2.729227, + 9.972701, + 18.46357, + 17.24962, + 16.66522, + 20.9598, + 18.17658, + 10.80134, + 4.07879, + 17.75641, + 3.633756, + 4.65884, + 17.17084, + 14.27111, + 23.64081, + 2.64591, + 10.31742, + 15.63272, + 13.63692, + 16.51991, + 15.26418, + 20.17847, + 9.149963, + 20.66901, + 24.80567, + 17.34251, + 0.2177003, + 3.090493, + 22.74068, + 22.66954, + 23.81961, + 14.55301, + 18.85983, + 17.55044, + 15.24921, + 7.34596, + 17.97639, + 18.21805, + 18.03173, + 20.16727, + 14.35985}; + Double_t Graph0_fy1[137] = { + 22.96098, + 32.37453, + 20.42714, + 23.75573, + 23.75573, + 28.97319, + 33.52173, + 28.46252, + 38.64605, + 29.95049, + 24.18883, + 17.91793, + 22.82172, + 24.2163, + 25.36298, + 29.18495, + 17.86598, + 34.39316, + 17.95625, + 31.63668, + 26.6676, + 14.39871, + 21.97711, + 28.91384, + 35.50192, + 21.15304, + 39.67074, + 38.94612, + 15.47257, + 26.06359, + 28.86099, + 32.18519, + 31.19928, + 25.90003, + 25.30423, + 33.05285, + 18.97107, + 78.38055, + 43.17075, + 46.28093, + 33.2535, + 25.42674, + 39.67903, + 28.68732, + 30.61999, + 30.61999, + 25.82045, + 33.3565, + 19.97845, + 25.77351, + 43.49009, + 36.59207, + 27.43145, + 27.43145, + 20.90391, + 32.97262, + 18.1204, + 21.22906, + 47.32078, + 37.22993, + 37.22993, + 37.85647, + 21.61354, + 19.62796, + 34.49465, + 26.35171, + 48.47679, + 26.015, + 38.44679, + 43.79028, + 27.86231, + 26.99938, + 33.08223, + 37.60542, + 28.98927, + 29.76159, + 31.48392, + 32.51186, + 28.55604, + 31.75072, + 32.60382, + 35.39864, + 41.52173, + 21.38447, + 30.2175, + 39.88721, + 14.35446, + 28.19643, + 30.12882, + 29.8008, + 33.42579, + 46.05483, + 33.75885, + 21.86448, + 14.45356, + 14.45356, + 18.33214, + 23.83462, + 39.67193, + 42.60805, + 19.64076, + 31.97388, + 26.98648, + 36.02568, + 26.63101, + 23.97905, + 10.75056, + 20.97361, + 37.34012, + 40.78341, + 30.53927, + 17.6882, + 39.23266, + 32.33729, + 34.5398, + 44.46054, + 27.14791, + 34.37672, + 34.71008, + 38.29793, + 28.51064, + 30.15736, + 18.32953, + 20.052, + 37.68594, + 30.23097, + 29.02963, + 38.83401, + 19.82276, + 29.83649, + 29.06612, + 27.19644, + 38.30944, + 39.88242, + 26.82334, + 24.67148, + 28.65588}; + TGraph *graph = new TGraph(137,Graph0_fx1,Graph0_fy1); + graph->SetName("Graph0"); + graph->SetTitle("Graph"); + + ci = TColor::GetColor("#000099"); + graph->SetLineColor(ci); + graph->SetMarkerStyle(8); + + TH1F *Graph_Graph1 = new TH1F("Graph_Graph1","Graph",137,-3.580114,28.89493); + Graph_Graph1->SetMinimum(3.987557); + Graph_Graph1->SetMaximum(85.14355); + Graph_Graph1->SetDirectory(0); + Graph_Graph1->SetStats(0); + + ci = TColor::GetColor("#000099"); + Graph_Graph1->SetLineColor(ci); + Graph_Graph1->GetXaxis()->SetLabelFont(42); + Graph_Graph1->GetXaxis()->SetLabelSize(0.035); + Graph_Graph1->GetXaxis()->SetTitleSize(0.035); + Graph_Graph1->GetXaxis()->SetTitleFont(42); + Graph_Graph1->GetYaxis()->SetLabelFont(42); + Graph_Graph1->GetYaxis()->SetLabelSize(0.035); + Graph_Graph1->GetYaxis()->SetTitleSize(0.035); + Graph_Graph1->GetYaxis()->SetTitleOffset(0); + Graph_Graph1->GetYaxis()->SetTitleFont(42); + Graph_Graph1->GetZaxis()->SetLabelFont(42); + Graph_Graph1->GetZaxis()->SetLabelSize(0.035); + Graph_Graph1->GetZaxis()->SetTitleSize(0.035); + Graph_Graph1->GetZaxis()->SetTitleFont(42); + graph->SetHistogram(Graph_Graph1); + + graph->Draw("p"); + + Double_t Graph1_fx2[108] = { + 16.53346, + 2.775385, + 25.618, + 25.618, + 15.84522, + 23.25345, + 17.87187, + 11.87716, + 13.51793, + 9.709167, + 8.177823, + 15.57414, + 8.79528, + 15.49065, + 6.447484, + 15.32576, + 15.44866, + 14.43255, + 18.51257, + 12.52081, + 19.88932, + 25.77312, + 16.85224, + 18.54831, + 14.66246, + 21.70794, + 24.23959, + 15.92302, + 23.08188, + 13.04678, + 11.01585, + 19.7312, + 14.52214, + 18.01523, + 25.71238, + 19.00919, + 14.77225, + 4.760872, + 25.25253, + 25.25253, + 13.78085, + 16.32254, + 7.666889, + 17.17745, + 12.18927, + 11.04283, + 15.78947, + 12.80672, + 5.574539, + 16.38851, + 12.18283, + 20.43209, + 13.05403, + 10.56624, + 11.00986, + 25.02393, + 21.30116, + 9.124693, + 18.07501, + 15.33917, + 19.30544, + 15.68121, + 14.6701, + 9.807625, + 13.23228, + 9.896872, + 8.678942, + 25.70097, + 17.79193, + 26.18868, + 22.8376, + 3.849799, + 18.20476, + 14.87749, + 14.87749, + 18.46357, + 17.24962, + 16.66522, + 20.9598, + 18.17658, + 10.80134, + 4.07879, + 17.75641, + 17.17084, + 14.27111, + 23.64081, + 10.31742, + 15.63272, + 13.63692, + 16.51991, + 15.26418, + 20.17847, + 20.66901, + 24.80567, + 17.34251, + 3.090493, + 22.74068, + 22.66954, + 23.81961, + 14.55301, + 18.85983, + 17.55044, + 15.24921, + 17.97639, + 18.21805, + 18.03173, + 20.16727, + 14.35985}; + Double_t Graph1_fy2[108] = { + 32.37453, + 20.42714, + 23.75573, + 23.75573, + 28.97319, + 33.52173, + 28.46252, + 38.64605, + 29.95049, + 24.18883, + 17.91793, + 22.82172, + 24.2163, + 25.36298, + 29.18495, + 34.39316, + 17.95625, + 31.63668, + 14.39871, + 28.91384, + 35.50192, + 21.15304, + 39.67074, + 38.94612, + 26.06359, + 28.86099, + 31.19928, + 25.90003, + 25.30423, + 33.05285, + 18.97107, + 43.17075, + 39.67903, + 28.68732, + 30.61999, + 25.82045, + 19.97845, + 36.59207, + 27.43145, + 27.43145, + 20.90391, + 18.1204, + 21.22906, + 47.32078, + 37.22993, + 37.85647, + 21.61354, + 34.49465, + 26.35171, + 26.015, + 38.44679, + 27.86231, + 26.99938, + 37.60542, + 28.98927, + 29.76159, + 31.48392, + 32.51186, + 28.55604, + 31.75072, + 32.60382, + 35.39864, + 41.52173, + 21.38447, + 30.2175, + 39.88721, + 14.35446, + 28.19643, + 30.12882, + 29.8008, + 33.42579, + 46.05483, + 33.75885, + 14.45356, + 14.45356, + 39.67193, + 42.60805, + 19.64076, + 31.97388, + 26.98648, + 36.02568, + 26.63101, + 23.97905, + 37.34012, + 40.78341, + 30.53927, + 39.23266, + 32.33729, + 34.5398, + 44.46054, + 27.14791, + 34.37672, + 38.29793, + 28.51064, + 30.15736, + 20.052, + 37.68594, + 30.23097, + 29.02963, + 38.83401, + 19.82276, + 29.83649, + 29.06612, + 38.30944, + 39.88242, + 26.82334, + 24.67148, + 28.65588}; + graph = new TGraph(108,Graph1_fx2,Graph1_fy2); + graph->SetName("Graph1"); + graph->SetTitle("Graph"); + + ci = TColor::GetColor("#000099"); + graph->SetLineColor(ci); + + ci = TColor::GetColor("#ff0000"); + graph->SetMarkerColor(ci); + graph->SetMarkerStyle(8); + + TH1F *Graph_Graph2 = new TH1F("Graph_Graph2","Graph",108,0.4340556,28.53001); + Graph_Graph2->SetMinimum(11.05782); + Graph_Graph2->SetMaximum(50.61741); + Graph_Graph2->SetDirectory(0); + Graph_Graph2->SetStats(0); + + ci = TColor::GetColor("#000099"); + Graph_Graph2->SetLineColor(ci); + Graph_Graph2->GetXaxis()->SetLabelFont(42); + Graph_Graph2->GetXaxis()->SetLabelSize(0.035); + Graph_Graph2->GetXaxis()->SetTitleSize(0.035); + Graph_Graph2->GetXaxis()->SetTitleFont(42); + Graph_Graph2->GetYaxis()->SetLabelFont(42); + Graph_Graph2->GetYaxis()->SetLabelSize(0.035); + Graph_Graph2->GetYaxis()->SetTitleSize(0.035); + Graph_Graph2->GetYaxis()->SetTitleOffset(0); + Graph_Graph2->GetYaxis()->SetTitleFont(42); + Graph_Graph2->GetZaxis()->SetLabelFont(42); + Graph_Graph2->GetZaxis()->SetLabelSize(0.035); + Graph_Graph2->GetZaxis()->SetTitleSize(0.035); + Graph_Graph2->GetZaxis()->SetTitleFont(42); + graph->SetHistogram(Graph_Graph2); + + graph->Draw("p"); + c_angn4->Modified(); + c_angn4->cd(); + c_angn4->SetSelected(c_angn4); +} diff --git a/tetra.tex b/tetra.tex index ba87e81..b46c62b 100644 --- a/tetra.tex +++ b/tetra.tex @@ -1,7 +1,7 @@ \documentclass[12pt,a4paper]{article} \usepackage[utf8]{inputenc} % размер - +\usepackage[final]{pdfpages} \usepackage[T2A]{fontenc} % указывает внутреннюю кодировку TeX \usepackage[ngerman,english]{babel} @@ -105,12 +105,35 @@ The missing mass (MM) spectrum of 4n reconstructed after kinematic selection is %------------------------------------------------------------------------------- Obviously, this spectrum shape can not be described just as a contribution of the 4-body phase space volume, see the orange dotted curve in Fig.\ \ref{fig:mm_4n} (b), indicating the presence of some resonance states of 4n lying at 3-4\,MeV above the 4n decay threshold. +Its important to notice that the observed peak is consistent with the reported in Ref.\ \cite{Duer:2022} and previous Ref.\ \cite{Kisamori:2016} results. +Its important evidence for the tetraneutron resonance at low energies above the 4n decay threshold obtained at forward angles in the reaction of alpha core knockout from the $^8$He beam. + \section{Conclusion} +The $^2$H($^8$He,$^6$Li)4n reaction was studied at forward angles. +The observed peak at 3-4 MeV above the 4n threshold is in agreement with Ref.\ \cite{Duer:2022,Kisamori:2016}. +The obtained results are pushing the investigation of the neutron clusters. +The configuration of the detector setup specifically for such reaction and improved neutron detection efficiency is demanded in the future experiments. +The latter will not also increase statistics but also will allow to study neutron cluster decays in full kinematics reconstructed by the invariant mass methods. + + +\section{Acknowledgement} + +%Institute of Physics of Silesian University in Opava +We also acknowledge the principal support of this work by the Russian Science +Foundation grant No.\ 22-12-00054. The authors are grateful to +Profs.\ Yu.Ts.\ Oganessian and S.N.\ Dmitriev for the long-term support and +development of this activity. We acknowledge important contribution of Prof.\ +M.S.\ Golovkov to the development of the experimental methods and useful +discussions. Also, authors express their gratitude to the +acceleration team for the stable work of U-400M cyclotron during all runs. + \bibliographystyle{unsrt} \clearpage \bibliography{all1} +\includepdf[pages=-]{add.pdf} + \end{document} \ No newline at end of file -- 2.18.1