From a66d706bf2bc1ca86a4fabb479f4ef7819b0d83e Mon Sep 17 00:00:00 2001 From: himyss Date: Sun, 29 Oct 2023 18:39:42 +0300 Subject: [PATCH] data added --- figures/mm_spec.png | Bin 0 -> 18413 bytes root/mm.C | 281 ++++++++++++++++++++++++++++++++++++++++++++ tetra.tex | 12 +- 3 files changed, 288 insertions(+), 5 deletions(-) create mode 100644 figures/mm_spec.png create mode 100644 root/mm.C diff --git a/figures/mm_spec.png b/figures/mm_spec.png new file mode 100644 index 0000000000000000000000000000000000000000..6dbfbba59ab20b0b4abcdbbe4b297b45947bd991 GIT binary patch literal 18413 zcmb7sc|28X^!IblF&4^@h)Q)$88T%aj#SDV2^k_&l6frQ)RiVhg<~En6q(C#jL|hu zHBH&ffcZo;7^eTHm#{>jrw7%p10CKnO7(J#x?x zA%Zvjb!VW1PaJVV+3-IauLDPqGr*sV40hMxXKwF9rrt(w4&Ht?p7zMm)y>6T%*)o( z-rm*gjGOlip-L4YK6Lb;`f>k+NvgT;gN6dsG5jRy?@wxHqAz|D4Lw^V`Rj_`UKgoqw1gXfq@ZaORT>{JFkMs(6&S+ydYB<^Lq5`$>w=I^}B?RQq&hy-M0 zAMZ&aknRQsS_JDx&~FY^0EwV0Y#lh_W4rgLChN9bmf9hMx3gqIiHA7Ga`VT6U$ix3|x&s?AARuC!!} z_S#o32pO&IG(YSXj>dcS9qVk-Sv9Zs}!ULiU*Q?h)fqp?rn4V9wc2taK z7&azbr4L>#b1d3P7raokr32^Q>BN#)piOG!<{^xGD%`P&&UWdUPhvnG)paBIj0saqd{mWnomRSPR>Nzr34ZS$t+RIgL5lAgFaC%)sav-=t-OFfWr~5 z^CK=;p?dX)?+6bgi}CdDyRnU^A!AnV2iiJtrjStG`sLMw;W@ z1vKPRI}=JNIaOjiSkOJtSEMyQs?%Ad8;Caw6@M_9kS>HeUygWc^yZ<4L=ubj zN$_NV-BdWyuogQ93gQDiVq=&_CF+fs-=_kW@R`)=^h$uYuf80SvZeUL#zBa1iZx(> zW%xA7{L$M|k^7$*1Px5oDl4KKhBXjM8wQPeo1GR`qOWGASUBix-3J$h z#Q8LUmR>$sgI@xK$744xWn%ogvnNg;tS6sgO_U5ru4MGS zfzpV2fM(V-AhL!vyx6sl++Ya}joF{I;PdH@|7G)A7qmotz)V_FbD~(9bo$<2fwXAAR$6&O~=1&1k4yiS!YbU zTfrz1;?M$Y1n9H`)4RD9m3EGRMXxIY_D$MPb-sp2`hcA1jjzZf>^{MU-}gD+y-paw zMEN8$^H|T45k5}_)5d<^FzKt}UaqF#^(vqfYcWL&{Hy?+$@CroOSyXJuOMc6Q->G1 z{YRB>^S#txq3+G>d@88c4R(dzcyqssHKDI&_#HN?w*@Un-y#gS{J)e5$I`~OdmS(l zIsXm+lKbbtzrOju#i@HNs6$3qOE|4^MsxLC|H4@P2R|Q6xo$5)3bFREP!kCH-reSZ z0z{h1RW(2Do+5o(Uo9;pJKpC>dxg#q4xOGh?0F#9t}|%EsG-(}xjCbBIJAE}Y5I1~ z#LGcl~LJZ{;2B;YBX({YyW7TznkHPO3>%wBU@pL*Uz|CBC{pVdCre@8{>2 z9j(#qGgSD2ec9h*8@&-a)7RbjV7dv&|f=i(4UfFVbB+P9Cm|5RGI?PB_Q zrsObU)X@O96IBO)ci>BjZ*ad(&}K#MBV98y;e%bcfkT-6Fyj@F02aV(9GjD+ylZE| z&B}QZvW5tA$AqA$&5LM0&3{^klRdh9(uNAsIInOSL}qjKk!DYLcUv;O^`2>34|cNE z;SxXjy3Y&ry>eeKM5e}f}4?p`Hc+wV(+X`TR;kDdzo%R#{`@AOtx?BDA>x6u-Q=Hk4E>Kq}G zOEWW*`c6Cs3BTB#(wxA@j8KWDQDj(b3e&&V#lkDByPn^UcHpXCVJN!mx3gnQ|8JXJ zzX?zaZYo{BNg^%aIWKA$5Ap3t44cugu_edUzV2GOTxsIsZzkq2?R6!_j zcWg=PHm`k2A5+yCXGH;;9f?3nkrD%FzYh?ewcykwVJuK@_$S}? zgM2%W4a4;z4rF_%#1B%%3Wj=yG;C=I9FQ!R4^3!*6iz$<)6yA~dL*68LZbZ#-e9 z_AX6`y8ClUFly6#0>VgiJoi~Ln1d5niEtT>HUGby;Mby@BBkfz1bi&XitNq54Y@Vpa0RI7adt2>PN&A=V zIPWDw?+4C-r>vd(5T~WMnwlEm>{x*gdds0}fEr#jb5j7*$HvAwd%AoQDLLA@o#h2( z5j*VF-99yLtKJ|bQV|y6bSE`_?X;pNJkP&e2>eJ9-_CM1BlxiU0`H$9Yw}YpOSbSaTSUa;oP9lC3RE?^uF zwPKupSEQ@k8NxEnl!IcMp4QrHN&d6>XdE>ypV|1n_@CBrxJZ(#@vmdC;;0b=)7WEi zQX*Ah|2*}^#_`ZPIzMiec{45x+?ojx*Ci#Sn8E{xm=9caR|wHbUbKT1!?u7VQ`XOU zP9TNnY}}}WwMn&6@XhC2_%z6KV^Tu@*)w3>YV4=d^0b+yS2DVpnP*?!A3=GKa>JK& z4Ux4&aP~4ysFhhCD3C8x6t9>ljuWc0zn(V9i$qSkuq5(LJAUk5+f*&@x_!ivbjl2B z?ZpZ83Rkng#pT^?#8dLsEwm)Z{sO?$8ObJ1zvjzv+QuJ+A~gZ{q~<)y%^0S=DKF_7 z{|TOypXOfh!BhC!QlY_x{srMQ?>3#oLAn6?A9RvZrT0xI6J4>85@gvm#y7*79F4OZ zK5^bYwABFo#qOu_ZFB2ab@DrlTslAb>(mcPUX)mCmb|rav5mP`{)`8w4T?FSqb)e3tIiwq0#!R1!frjxy@>-a-hqqnDEgKbh8`6m)85@lKRiaEL&uoN2$*lXt#h z&(+#pQ@&-#jg7yynuMlBhU1?N=8pwI93)n9ZXDOU)T^YbrdHT?hWXFJ+QC_|GWiy* zBqjn8M3%pD^0eRr=dR$t>{BA+hd%jF`<`zS?>mk2C0*`48+3n|jl|TFlv+kUrj%{t zS1Q|+J%u^c`#)11rEihvm@+P+SvtszIpxLzh79|q+qZRqsb)t7YXx)m1j(j2foXOu zO!Ahe85_?}#q<=ZsEKdg5toIMFV8t+%fw#{9+RWXa|l$ytA_|y4&2$XP-Gj zC*CaXC3#V=op)kZUT@?fvO6+BN37$tG`-kd=g=G4!!QbxGHK!LXx;KK(je(ui&u+e zGXw11KA(wtJdxLnIV7}x!RW!g=B@*wEPS?W8j84NMxZ?AiLr2~Oh5Zn>Pos(<#sn? zV>xLb-Oh^}zqyJfKuR|o1{XEAbe6?FqPog`n2CR9kBWPya!$lmS$HcK+ahZs5G^gO z+d3?xId`yJEsM?4T&!FzI^(CcYPagi6pE#mI#iVz(E^t013}A`s|fT}bnba$P*qja z6EF%dqVGN-CE7sCE~;o(ivrg+)rIDPXA=}r!ZDz@lwW`HS8l7gS2>e=nc_FOu(VnI zn6a_bCACe6hi$Mg|BK`W-Ca9bb5`@0EdEH*cA}(+Ye|srKHp;PvoP9WZ{f*3*vC4^ zOnnKc*|T|Yac5ckU(sS`;sk&fAuAabxy0YDw~euelMV>pr1cgAwr1#43`X)i2j<49 zsy_pMoKH{y`R)hAZuAGY!O(BdfyURrp4naml^Me7aqzYc`ybCo9k#H=n6{_%k9O9QQ1K%23}Rd$42$%D$gW6`)mPnZ=u*kwj-X z3uny@8HT~NMrrTg{mJC$v$>SC>WWl4{jIV)2p;Q(n3tEMzVkzbns*=C%F{vjA=W&@#e;9)L^(PRvsO3 z8ob-=7$)4@GWl3v?`ORImB$%)Zw#EthevDTytLzr;eSH zRwa{vaJ;}(zu&+(imnOXyxUwkW7H}0I5YEU&(sb&FCWW{<)9g(KTx`$J0V3PE6U<_ zLgxLhAm{Wcojc(-x=oz$6-5G5qGb{k4W;1ghkW_Uo!|NLS}~q-9U)2LErT3Td215qFPQnp&usj|IXMN%2oRjb!?m0rFCn^_spmSJ7bs9@fnK9}{ zQ4c3^A|yaFMBGpNrhC{qjgX1pN!9+Qj{UlGZQ)hfGfTz`n~0b6V&(10=v=iqHzhMO z{i65-m*snZK-6y;IN>xBdbx~Y;Wk5Z$G)Fl6b!{$wdrU(&frloo&xCet|;+iDXRV< z8+nUBYK*)8B{12yjQYB(rzp`Shx*Rz+vu;_@89N{zm!nkJTg*FSq`fDUhOPSXR=oF z{M!t2Rj)#N^f(xj{-KH5iwa2A{xd0Of`WZg?#9)Sh*ts(xe~9OVq^uyU>;L5C)!+_ z+mDGyA-MIE_5?-mYOT%fQ_N_F^%X8}*FC>doe#`om}esQ4(uxy#H)K z8XMA=r}gJ)Z962Cc+&$5eg+5s|KEMfj13shQqo|or+=~}F9 zY@7xG&!+^?)>_pq%HnHaGNd4GK~<12%IJQiOec^H6ZzCRZTth z3$(Jt%sYIj)B#Ay-W6@>y=CW4>)%-sD#XJw30d6k@>KI*SDHQ1=<>0|sCbsE@%?T+ zvio&y`>311bKTfFTnt0xynh*w+Re4BUJ*R^g>2uTR2tWR_8;tvqZ+K`s%6g;Z^!ln zH?S6pYIrEWe6-4++$^bQDG^?8`)+38an6KWf0raXsnOJcmxMT5JE#wAM_p9K2>u$? z>zu;DZy`ayVZqY#dsy=2vdH(|Eixte^6_K7u0xkBF$4c+k5%7D(=Q*S@Aa{Cqlk76 zY?uu-SEttfSsHozzo(9eg29#6k5*;*Z#XIY;^O8GmxKR$PH-oV+VR*n9eni%>wjbP zCiad)pgw4PGx zTBDS*k2(Ggh)oI>W$Q@9muhOWKK!CD0oltAS1);98XwfVFyZ&sQyaaH4P>$G@{uEy znwJGemX>kYh$+rbt_MjQme${K{wxbse`Clg*MP@sgk5J)mA;Q<@KRr4m7hySDM`2U zT&0L0ILZ?Z+53o(YEeTBZIQ`u8V-2z(|r0xj1-H5V9FA%q@l}X@ze7*2k_Ah%Xn6h5>$xe)w7Ef? zbE6Q^Lysy(>yBtXcNr>eWABLCHq4Z~Ai?hU;pv;c5B*3-$AF%=84{)NFyX`&1+>j! zENFuk6cjD(bbOmz4G_cGBGRin&=Ulg&clki%&eDAB zMT0Bg=bvleb661MW#{-n@6qhxyZ>;)JW@edLBji6>w)6wKWXNF7 zBN~Z`n)kvgRSmvc#arJ6aeNzoiX1R~mr!!r>t8J_zjl^%$e({a@?{Zs?(c0o`wnsS zkb^8Al&7=azW3S8^vAnNOYy5>i#wfqd|&kz_EHY^?JX_el9K*q`7p%#TAc^aw+vju z$E9%9XN>Kq{bqY7U8sKOlZ}GCW!Hly!5LaDXUav@@T4ILt_qWA#_A)db^e=2N=U zVR~ByhoK;;uiyL3A8aD+WNH?*A4LpjAm7i@V)U+4Z8yBVF@+1xkK-dj(OWqqxKB7? z$gk3wD@A3lSM(rcTjkDg>^FvLb3S!ZYNU7o(?E6gTUkv|v}7(bS@}|))O}*34M|D7 zSwyZOVh$D$Xat949ay@Cs0j zi%abRw$|L+#zVbc9_lT+O4S=1KZgxP)|6@ik7xy__pItqUd<d2o2a5(t8(Q&Rv5d z<@St#VSnZHBA47sSJmy+p+6Rvp@pRE@t0qw{@v%eDa|_X zjx+1X8~}X$R)5-SJpHV-r#x@ls%$b#ia}fos{RA)Oq`J8SLjThco(ZrpNE5{ynpu- zh=!=BXUdA@BPFM-=`;0Y3G|crZUE~)(VKj$IwSRRrR#G(~A0_Ybu%{myvihGgC^|C?u`@b^eS;`uhu_{oBSq zz8q_P8z|1!ktQ>Bg=ggSm&0h?uXUvW9jwfsG>^diOqouBi}}{w#1As+I8vi%bSe+! z^}z6*A{Px=;-{l*c?t`cOK!DgV8&~<2{54|e|g(lYB<(L4(tzhj$xC$#v?V0;`gNd zDSG~%4gjBA(J7^rAyvYrJf~W0pdstCKrG;c+qr!y&E!k;|K7xsbKRJyc`5h-#epVv zO8c=E`l=qq%MPB9DE?&F%nnbDm`~wod+nShbBJ4!);ptxs#4&KcyjGN9baYOcO~=HG(fDZd-w$a2!jX=O^>K`Bha3$+VP%yex!}0*{FfqDt4%4J zpKFWZtcCZ{LmpRG`1~L+s?pLNLL-vu+wvGs?F6@udHx>jrc~N@Lj{ime$$%sP!y0~ z=s6D(J#54Qfe#yqi@D80MtK0IP{MMwKw}Br{7poVoPI2d@zA4<&+aKB4YX5pbaVFQt^aXM-1d>Bc@nYOwKa+~H`}x8mx&NK}v9e8t;RFb1d)fqBnR zOpOc6A^P?`;AdSy8&D@DK@A%uyBVO65>Ujr;e3l6O$>Fddr^HFn)2tkzP<6O=wk<1 zl!Z9k6(ddgBOFs5S@CuTja&)l9&8kwhChbkE7VL5AQKt4L~HT!bMRbOM_5`JLsZ zk?F$E61*D!9Qq{RGeg%`_{LLE`}RiGXm80j!@K;sN4x5j-*A~Mi**U6#CJ&W$QeIQ zyGir@*;^B<9|HUK=3d;T;TT4U=CEloTjVONG}%I_GWjYD z9h14~)+iVWK__fjs0ZH_EqlawdeR98o_bk*Q7MJ1Do8&XxU9a6VNZ3-uEalqKv}ir z+nOqaO0q-P!V9q+llXOam_ie4bAP|c(s8h^0J||{jg4t0Ip**pJ8YLpUK>Wm|5U`9 z?zh%W3r4O?LmJDfUApyZlJaWPBw0)iLnJk8BR(%IG^x^hbKjJhMa5C>1Ywt8;yKp3 z1csp!+wydbTV9*=4Ej2cP5;;(fmUr&2KKK*p@67AzndsN2P+=|a_i7|9Zm0v118xo zztX?pv}w66aa*a*`8Rbp8v}+Pfz6!xW6l$?9kw`KMvi4~!4Wy3ZCv8E)=wo(Lc2Gg zff_^VDfo~%KQC!6+r~vNDqIuufMEB1DZoA%CNfHE1(7#SR^g! zsV9m8!FZHNU%fLNFQ*Rw)&vF(rfG&&mlu+QP2E_~ja1{oGw2&*#JX%>;SlB0&PeYo7rFLU{)FJi`+(93=#8M*z zK|n3+^N7cUHB`7yVR+XD2PB<_#`OJ?*TV+|U_ls!xi-*uvGh`wI2jSYKG7iOGfbqR zV%?x+HAzkJns=&1i(gU~(4gP^+#WgJQb2EMxG7|U5h^<6L();f=Re_iOj+{oSJ z1*T#P0AlX5&z~K#BnKFpH4S-*28ztcVkJfwxKqPuuZh&eu-owZg~=6d?HCcMsV4NE zs9}EibQzcbd)d3%Br(GijA$2}u?U-!ds|-hCJ2vO&b@%pOEg=b0KjJLcI))fy2gK5 zRI?u!6=k9MWfD9+8^=Fyu^GU*iVV^jzFQ%^{k_vh{_5y6A$XLEJ^P3?J=*nqljNLL zO!;#bFC62 zBVR^H{4vQx8c`fId)?Y;iE#SnudBEfNn$WsnIxJ|gPLIA(E{d}+{)(l6RmY^bN&)8 z_fb(Rjrb;}y(2?%JAKaXJ+|7wUa!+&-=I|+? zluJnKSj+sQFKh$p440~3YBBg~^?>Vgbl8j55K)@D3O-Wo0AJoyw#$hxzk)NaN)LJA ztX*MZmec{mk3xx3+K55B(Y|i?!;h&PW>(>GsYqGuz~jD`T0)wlF2T8zMveT$$J4Yf zP?9mgBonbQ^m+VuZ?yULU2w=FuQh}%!jb(=gUsTmY)!1%*y4a22KRF}aVsDgPBDk<{4@=py8qUK>ax|%8xWoUs{?}$3K3^zjnxyeUeuo%nqE8&iGkUjT)tv!0od(1- zs=2pjg}x~qIhv$`Lk*Ln*&&n0PME_21n_e^ypJ@+6aA$CZe|-wn(bYk+mJc zfYTu^hM!yg_h}&bu8-PA6gGeNY7)hfa9|sEF9x2ui$@9I>g)56e2(j+el_$VxM8;a z)z_6siDU+96Dkrt;P)aessmT>6hjhza;(yhaD6q>tvj~}Fj2Z@fs;gG&Z-@5v|dow zmx~^7aCG9g%Uv^!$W(&Rm9Xu|U7biDv8P)e4KPwH@WWb=7!1%uTT%4D)c(ktc?N?b z2Cv-c6E!~y*31SHl;P0{4OtfCe}l=gZui9kLYnF1h4QTpu}1W>(G0PPX;Bdsd)Mq` z9lC+!9=40Qj=9AT;xGGJQ=HIi^OKAt*Vg2`7wZWMMS@T*n6yRrtK(1-h`Qkph}!=$ ze9C)7%#($fUwsp7 zz2TC85USot<1Gzi=|7GL154hviB#5-Oh8F?EL5iTVtP}tQ;#eVjqdRb&$9;6@Mu$YV%|y zcHLv)rVwZvbG85vhV_dz;f9%=ELV3HdJ~D@QCK^i9ftLRK#zPmbhoKlw68);}X}mX< zFb%^%y&g5h~+`B_CZGAJd>^vJcs*X+Ge{Q-LFwus)a$b7m@r;dPv9~ngt;EuauKDvWiq^KSQ936zSdeMV-DXA@wt@i_t)mk3 zHbP_?bmAe5zt$Qy3tOGA_nN!5{FSlaS=^SWn&u*|37yNwImqnB)<8L7Rp^Z{k#cfB`TldsVY#ZbvkM5@>`D?{!K(;*l#7005!X+O8Fz_C0jAlsJ z{u)W~>0@ST;7_uTLzu9!YJDu@m~w-{gS17sl*T_AQYE4THni5=tGv5BhxN9Hla&#= zrU{{vv!M=$V(+Fl$Ane?{7$m?vclXb;)CgLlGxpR0+NCuJ^$e@H-LLSb8r;Fa^-Hj2gwzZ;@2Qca$h=s(BEOfY*^@Y( zagJ9`m+{C>he%dX-}XeOhy^voz@v7-7E}1*?%V0GhyGmibfv1{;(U+u284_1;5Y;? zN4W8gu@3s&au3sAKT59UITCJ`J7&p1vh5ij(F$}mL& zugE70Ii-OFZu{nXr|i<^l-AEQF!KIHgZ{;=&ZsVkb-l}*DYq`xi;X<73c14Q1q?V6 z^;*aTBzu@he~t_cVKt0RlUy!5t?Lzj#NLZTS72h@c6733-ENYlY<%0`*q>rF=<^;X z-5HKT2YD!n|g3zixGd)`$ShKuu*j+8hIoc4x^`EQkR@ zedO5F%7p!J{Qzc7*SoNT$E($C!I!pHrK9RCpf0aGUK}E7LH-~C**YYg|O4=PIxzhb63EYL9X+) z5;6B46A(TkmC!Ud@1ex*x8XhHDhP(|A!>fzTO;|I1{?z!YDsxlbH>Qgqgf8z_Pk1S zZUsQ3Sj+a9u5ts`8&J%64H=73e$C6k~ z`zs(@LH7yu@8BS7PrRbFYPgmdtsg>em{G8fk`Kw(NrXP{$?fl2-2{6D z_Z#&&N%JE^!M&ER?bq#dr7bOtAE9E0H}$6D6${h9upw(hyyAf4!f5TSWq76S-lFrA zAmVt>@vwp|1-sM%tA%!ZlfvA?{^IHU0k6dshXn?UVwtP8V2vc_*u;>BG&_zLRiC5r z=3(3Q_PTr6A=VNcdeWrm4ZVGB5iv=GB2=JI!CK@?Y^#xgZ4{!i5ZjhGcnO1e!fu?@^6%%yY0whs5E$jr2ymj+j|Hlp#4C^XmAO>+NHwYu^>*@$y`4PG&a zRL_~F^{h>4!4~X*>uzc`ZwzaExxAQ|xw@><>V-#q&UnR!=h-Y#xe?b1h+zaqiL0Z5 z9|?LGX$WlY1jRVHQy)0X`XP_049Tpt@K0xhc_#Hr7o92311X%v`Rxq#C3L~{pBIWl z8FHNfx=%9Jn{M|C4yA=oKa91N^+UUXcUa{#pU}l*IFI?MV2>bRY z=ko?=leZEvSrdv-9?#iAFQm#(p0{9M)$q|bz64~=kEbv)Hw(E%8O!RN-V}> zC1d@hJ%m-rkpGed8O8-lw$)YE@k-3}1p7vuumoDGlWexkImGYLTi?oN)@mT>Jl_{yQWBGV=W7Ct3^oR&A zM6$=0ncD-6T%U02)P<^ruPWU$tBc>}DnNA~Cl-QM-`1C@ovKkSu2cg7E-Qcl+a=$f z5?=BS>v=nNvR*!vBu;3j`LmX`DRB9FZq;0wq1KpAuek%Te(7J>AK zY>3oiT%W&-PY!YyjKRj<-&h)2^sM#HR0Ec2cEc~)P(VsFKA}u~1B|U|;VF$w1PTnR zy7$;rJAQ}Z_5qlpW`qG&O*6{i4@&arJGi%Z4@PilF9;d%ILrcVhL9D!Y4;+OD&po- z_T$46P+L2Nm)U0<-+KT7{b2N5Tt9c$huQ}wZMY#r3)^-cTSJTDy6c@o)PCFY?jpao z?Y3=Iw{^}|_3t#&>fhmGxpRV2x%f_BZQ)FS{El#wvDfS+S#q(f7X-=1e=PLXa=D$39qPTVI?6D9 zLRR_lOkx;WP&I+8q2y`*Ot?!=-}09^?%L!EKUbAP=8>59<=lN)7r&JR*2?vUuk@(} z%3J^1!qs#k5$Cwk-g(ce{{BR?IO+T1?}fm3Q@+A!Joa9gRer}0SJJB6xs)i^N zD$U{3a?G=Iw_8rEmSq>rjHW{?Jsi5XE_yYC?(S`$_UzpD=|*#Y8kmLc%h}d9V`CYT zPPNeS`5Id8(z8>0lrL=zler(Ba=>LiWcWM`K2UOCI@PcJDp3BkKD9VDYpo2)4}bT@ zRL=N_^p@mI7y=%I9np$@GpJzOj_K=7Dj4mvkhq@~(8QZKnTy=^l;Z`vHkb)5qefvqF>g8mB*B^SV+REAmW* z%2YDpv&dSnufeT8o2_){Q%Uzr+4REU(!zPn0Cu|- zm}!LQXZllPJe3+BDt#T8c2bhVhNF8ncg`yNzj@&0T-CqoKYgq+ zf2?&j;nBq6$0hREEsz0uAzCurJWS3*=RAfnK#f)hRjI8m&CN-}fEL&~%$REncllIM z^*Onw3hb<*1Cg0osS+AX9XwB`#+|e;!)roDoLOz^T+fNuNH}d~tYy%>WLvm+7US)v z(E}9YX@~MbZt~b|$+1m1OS-u9(hG%w^1PK5-2G=;_suBm2v>_vPM~O(GD_7+4?*Wl z)$aY|8wIasC!PBiMshI{eJ>Sv)cQ$1{UqMk+fcwp{l9-NLyYd)@RHdM{;~xMQ zaP`{m+`mv*c^j1ii7mn7_QVao)+uF}h8IGg&7|r2JSo1N;S_f?NRG$|RB^=h0QI+~ zMs%S@Fn80ls_*jz<{36E5c1$PaMaMi!Tr>tz+I#Ga?7gQ%xzDe*Uu$LiZ7`q*7ns+ zDI#-$k6!M|(byI9N@ncH`SGv=Iu1J)Kdk00*>V&g@w7kfGy*QMiD(~>XSHOOZ_asY#uOUOl5;RNTsEmmIrcZloX?3LK0*nWalo31v&Db80Pb?#RQ$o&|N zlCc$tik24y)c$HWzs0I5Rb>7Ia#KY;iAUchgMj}Q3_kYS+Y|@w{&6|e5qFD((^ji)s4@n znC6=%>|bm0v1N3oxs^D1GTF~^Yn=zvkbG#)L;+9E6K0<@XR?5AwG#W!r|^G_E^X(j zmwdsc(Z7)D6bEkwk1wg}baK>rSoqvZ@V@2jO#U$M+}p@L5Vf*21x;_&P#Ymny?0_3 zKMr+G%!?#sXRE+s>8WRV>;H0weV-2VGpYy=56^8q9{nY->U(!-@3z&Cs|3y@tb1lyO9SeMxP}< zrR_iW?5jfkm)PevTm_dbscd}uAkYoj-BSg}sZz{0`DJs=M2bf(xQ-XoG#J6_h#Rs& z)`~hm7|dmjk7l@(=B2wbMmQdzM_W!Bl__VsmX0)B)+{Pl%eG7()TOcwNMdS!b{yNa zeW#3Hms{LuT2x)~Su$8lICkE&mjhC(-y0|+(pN9fz=g1ozYdu3)5Xe7WO$1qYt%jO zspOGHiam4x^9Qw!Ac|;4cq1dfch^wbCU*V))fREGlR;lu=w?;I6AFK|R6G1Cjgz{a zPgJih{%w~L`r*}G@7jpoQ6c_mo0amusW|5>a$5Qs%dRUnT+=yOu8etW%S9ikm(e|? z{%C@tdoNK&ru@Rah_q)d1bv1P-{O`jnbJ}BS*N(l><`^DnKMge&S^F|!+J?5_zE<| znwe@1ylZbgke%@#dH&A%VP|2fA$t~#BM0jB6y3D(-w_vX+L@KJ(cJ;R9qltVH&y)8 z`+}^=1;6Gm--~UyY9yHa96NP8zs&Ya)|EIjJ571&r4YWsUs99G1}V+sn{C6ob#+x$ zRfCv@hUBMLemwO({{Ry^bPnY&^&Z}%sr%$~SUmpKnhl(uXryLAyWumKAB(iE(KX2( zUlp6C7Bv&kM+fb5zHBgEW-(Xk*vDSuNTSe|Y4*c7vkhSetOptStat(0); + c_angn4->Range(-4.558442,-1.575,30.50649,14.175); + c_angn4->SetFillColor(0); + c_angn4->SetBorderMode(0); + c_angn4->SetBorderSize(2); + c_angn4->SetLeftMargin(0.13); + c_angn4->SetFrameBorderMode(0); + c_angn4->SetFrameBorderMode(0); + + TH1F *h1__1 = new TH1F("h1__1","lv_4n_mm.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && lv_4n_mm.Mag()-4*939.565>-1 && 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)}",27,0,27); + h1__1->SetBinContent(3,1); + h1__1->SetBinContent(4,2); + h1__1->SetBinContent(5,2); + h1__1->SetBinContent(6,1); + h1__1->SetBinContent(7,1); + h1__1->SetBinContent(8,1); + h1__1->SetBinContent(9,3); + h1__1->SetBinContent(10,4); + h1__1->SetBinContent(11,3); + h1__1->SetBinContent(12,4); + h1__1->SetBinContent(13,4); + h1__1->SetBinContent(14,6); + h1__1->SetBinContent(15,10); + h1__1->SetBinContent(16,12); + h1__1->SetBinContent(17,6); + h1__1->SetBinContent(18,9); + h1__1->SetBinContent(19,10); + h1__1->SetBinContent(20,4); + h1__1->SetBinContent(21,5); + h1__1->SetBinContent(22,2); + h1__1->SetBinContent(23,3); + h1__1->SetBinContent(24,4); + h1__1->SetBinContent(25,2); + h1__1->SetBinContent(26,8); + h1__1->SetBinContent(27,1); + h1__1->SetEntries(108); + h1__1->SetStats(0); + h1__1->GetXaxis()->SetTitle("E_{T}, MeV"); + h1__1->GetXaxis()->CenterTitle(true); + h1__1->GetXaxis()->SetLabelFont(42); + h1__1->GetXaxis()->SetLabelSize(0.035); + h1__1->GetXaxis()->SetTitleSize(0.05); + h1__1->GetXaxis()->SetTitleOffset(0.8); + h1__1->GetXaxis()->SetTitleFont(42); + h1__1->GetYaxis()->SetTitle("Events/1 MeV"); + h1__1->GetYaxis()->CenterTitle(true); + h1__1->GetYaxis()->SetLabelFont(42); + h1__1->GetYaxis()->SetLabelSize(0.035); + h1__1->GetYaxis()->SetTitleSize(0.05); + h1__1->GetYaxis()->SetTitleOffset(0.9); + h1__1->GetYaxis()->SetTitleFont(42); + h1__1->GetZaxis()->SetLabelFont(42); + h1__1->GetZaxis()->SetLabelSize(0.035); + h1__1->GetZaxis()->SetTitleSize(0.035); + h1__1->GetZaxis()->SetTitleFont(42); + h1__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.Mag()-4*939.565 {(lv_6li_side.Mag()>0) && (ND_time-targetTime > 15) && lv_4n_mm.Mag()-4*939.565>-1 && 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)}"); + pt->Draw(); + + Double_t Graph0_fx1[84] = { + -1.75, + -1.25, + -0.75, + -0.25, + 0.25, + 0.75, + 1.25, + 1.75, + 2.25, + 2.75, + 3.25, + 3.75, + 4.25, + 4.75, + 5.25, + 5.75, + 6.25, + 6.75, + 7.25, + 7.75, + 8.25, + 8.75, + 9.25, + 9.75, + 10.25, + 10.75, + 11.25, + 11.75, + 12.25, + 12.75, + 13.25, + 13.75, + 14.25, + 14.75, + 15.25, + 15.75, + 16.25, + 16.75, + 17.25, + 17.75, + 18.25, + 18.75, + 19.25, + 19.75, + 20.25, + 20.75, + 21.25, + 21.75, + 22.25, + 22.75, + 23.25, + 23.75, + 24.25, + 24.75, + 25.25, + 25.75, + 26.25, + 26.75, + 27.25, + 27.75, + 28.25, + 28.75, + 29.25, + 29.75, + 30.25, + 30.75, + 31.25, + 31.75, + 32.25, + 32.75, + 33.25, + 33.75, + 34.25, + 34.75, + 35.25, + 35.75, + 36.25, + 36.75, + 37.25, + 37.75, + 38.25, + 38.75, + 39.25, + 39.75}; + Double_t Graph0_fy1[84] = { + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0.0181, + 0.10863, + 0.14484, + 0.23536, + 0.30778, + 0.59745, + 0.61556, + 0.92334, + 1.14059, + 1.33974, + 1.93719, + 1.9553, + 2.7338, + 3.34936, + 3.56661, + 4.54426, + 5.37707, + 6.08315, + 6.93407, + 8.52728, + 10.08428, + 11.58696, + 14.37507, + 15.2803, + 16.85541, + 19.44437, + 23.78948, + 26.50517, + 30.19851, + 33.2582, + 36.60755, + 42.72692, + 46.42026, + 51.41714, + 58.315, + 63.8188, + 72.09262, + 76.90845, + 84.91069, + 92.27927, + 102.0015, + 110.6374, + 117.1731, + 125.5013, + 133.3949, + 142.1394, + 149.852, + 154.9756, + 161.9821, + 168.4092, + 172.9173, + 176.502, + 175.7959, + 174.782, + 167.504, + 159.6647, + 145.0905, + 132.8698, + 115.598, + 96.44334, + 76.92655, + 52.4491, + 32.57022, + 15.96828, + 6.75303, + 1.46647, + 0.16294, + 0.0181, + 0, + 0, + 0, + 0, + 0, + 0}; + TGraph *graph = new TGraph(84,Graph0_fx1,Graph0_fy1); + graph->SetName("Graph0"); + graph->SetTitle("Graph"); + graph->SetFillColor(1); + + Int_t ci; // for color index setting + TColor *color; // for color definition with alpha + ci = TColor::GetColor("#ff9933"); + graph->SetLineColor(ci); + graph->SetLineStyle(2); + graph->SetLineWidth(6); + + TH1F *Graph_Graph1 = new TH1F("Graph_Graph1","Graph",100,-5.9,43.9); + Graph_Graph1->SetMinimum(0); + Graph_Graph1->SetMaximum(194.1522); + 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(""); + c_angn4->Modified(); + c_angn4->cd(); + c_angn4->SetSelected(c_angn4); +} diff --git a/tetra.tex b/tetra.tex index 1fa0dbb..ba87e81 100644 --- a/tetra.tex +++ b/tetra.tex @@ -91,19 +91,21 @@ The missing mass (MM) spectrum of 4n reconstructed after kinematic selection is \begin{tabular}{cc} \includegraphics[width=0.5\linewidth]{figures/triangle} & - \includegraphics[width=0.5\linewidth]{figures/mm} \\ + \includegraphics[width=0.5\linewidth]{figures/mm_spec} \\ (a) & (b) \end{tabular} \caption{ - (a) The dE-TAC correlation providing clear separation of gammas from neutrons. - (b) The ToF distribution. - The set of three peaks on the left corresponds to the gammas produced in the diaphragm installed 20 cm upstream the target plane, target frame, and CsI(Tl) array. - Green ad red-line histograms formed by the events, identified as gamma/neutron by the dE-TAC method. + (a) Correlation between the neutron energy in the 4n frame and the 4n decay energy. + The kinematical border showing the $E_{n} < 3/4 E_T (4$n) relation is shown with the red line. + (b) The 4n MM spectrum projected from red events in (a) by using the kinematical condition. + The orange dotted curve illustrates the 4-body phase volume $\sim E_T^{7/2}$. } \label{fig:mm_4n} \end{figure} %------------------------------------------------------------------------------- +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. + \section{Conclusion} \bibliographystyle{unsrt} -- 2.18.1