From e3d544cb53a9ec3badb49bbcf12ec8c294769aef Mon Sep 17 00:00:00 2001 From: Unknown <unknown> Date: Wed, 2 Oct 2013 10:00:25 +0000 Subject: [PATCH] --- NPAnalysis/10He_Riken/Analysis | Bin 65696 -> 65917 bytes .../10He_Riken/include/ObjectManager.hh | 22 ++++++++++--- NPAnalysis/10He_Riken/src/Analysis.cc | 24 +++++++++----- NPLib/MUST2/Must2Array.cxx | 30 ++++++++++++++++++ NPLib/MUST2/Must2Array.h | 2 +- 5 files changed, 64 insertions(+), 14 deletions(-) diff --git a/NPAnalysis/10He_Riken/Analysis b/NPAnalysis/10He_Riken/Analysis index 26e6e27779acdfb7f4c0a90264d8b8cda462c3a6..cf79ecf074b7abdf28d3a2bb5b17cac87ad2e359 100755 GIT binary patch delta 21856 zcma)k3tUyj7WW)D8xDkm_yR>aA{G`tv%*p$^OWz<tgN()rbT9ore>*aOizJyQrUXd z0>#HQG&A$`pswD?Ctkg3WJV7RU)O_LQMZTr;`^^TGn~cg-tYT1zd5u2Yu2p!&zi@c zz1e$8=Gm{lWM3Q?k+awC=p}?$u%yG|VM~jJSX#VfVVkhDP$5P}FWEFKa&w*#n;nK@ zv%?Cp&5SSk_psh!>+Ns%UgYha^<uZ1!(Xf$F(3H{TT5(XY~8INZ5>Rfb7@4#xQOPt z_TgEBkv0@z7{YLb5eRhLk1!J90fdJTMj_Dgus@Q&kKlK-!ebDRRp%2_`~>1B)wzzI z7Q$K`Vwf6FRwq*s=cx0kDt;F6OoZnYrf3$z3kY)w;K-G~LQowSAh;FQ(aSh_1p!`? zM?QWR3SkNlPF_P;gs=p`Lij5J9n1ZZSgB%)UPmZGKof`(gf$3j5!NC64PgU<K1%h8 ziZ>#D3*l{qEeP+Bd$;05AMdJ@_Yi-GP^R!lh(AX7JHmE^9SENwlq2j!*oCkgfsVb_ z)zG%>4&m%D!V!ee5GoP$anx#O|IuE2tgq9GGHM*R)=c>3Z>_8QTXjZuze^EXONdyU zp1ZCwb~_s5&%+wylMd?zM|;~I>!Yy5kPq$Frik{#I~t93w%QtF5oEJ@5Qnw)#`^!U z28Jj08yqBkFQ}Bxk&ScJ1TpOqYF!F%ZypFL?A4IQ*{XufE)OERJ;>~9LHetL()%?o zFIL2aH!g6M)!UKStRTqmJFTZ2iRP&wY5O1_qz0w;4f4T|Ao8CERkpG*IZ@J>TXS(N zsDL)tS#2V^n)V<bvx4j%Zfr+<6WX}WF9q?e4#KYn72Bb4ep8$a%5XWTqF)3N`Y=fU zT2PmKKBz*Sf-3MrqgIA;q{F%t*50<ndN?xCoDt;4!$E8w4#H-TepygyIn6Os{b+ZP z7m`F15#L@2$4_<##1da~3URE??&#sirYMUMqV&AoF-IwEfx(Q9hT~P>W}=OVO%NjM zg55zp$l&QHAzu8q-O)nfD>p#?bD)6ocp+BQ+a2VUW+GagLjH@e=f`ct-6(Gi2qww= zW)pD+{+xK)?w~1z3V0qF7F-S(e5Iuji?1~>77rpl^;5g_ORTs~P{w0OCpT06qsh>R zy=2%WJfVXSlTqWgDxGH3Hm6jDOej>h65{%+c1IJXu$(lL4F|O|d0+?%tV02G?TUqd zCnzpsE-nJjQ68Z5Kj7<ZcuY6JAXIe~3QSb?Pc*|N`I6m{3v5z`bezmQYj?C$3N7H` z+>3q%O232xYT!buG7;(q7iNH<cIkJWgqVc<)KRI%-!v8C_!+w+RpA?>h3Lr$^#>Ct z@;m9GMh5qx&Dz6dn&ARq2Rxz)8Ss!`gx&|CLs#q$lq(C&xE>Be0T4^v4>hD!2%z*> zF^1Ib@PPKn_h^!x=#=RwUuLj63UP&Xk`naVZD=!yW%?&*<Gv`s-)WkNTR~_c(lx=N z)<V3@3U~+!$C=9)k|LT&cjBG!NT2WRj>D$vB;UjdagAl@2LeSH6`Ju$47;LTx<*1w z2a`gK0_}l9cwhhsA(DhUpbw3P2Wc%tO|ly0OxhjD-vNYj|FAnu3kuZ64`MvcK@;e4 z>V*8DDv)YD1O!7^FaH>Q<R6Ui`|v=g<95d&W$#DekI^SkQAr@r36HT(I}IKz`O$7V zN|nN$P}qqAG=WeQkcK8fvWO8s)&GJ~Aiy{m5fMHE#r~`l?naG$tdqV)9$=&53*heX znC>I%(YD*Jns!Gf43djhqSt1Cp>ER*5GX+>(QW=VI!y))>T!Mx^hdEKxBxp(qEn7i z_Q<nI@Q~K`g<z=)e6vb_6n!d+Pn8dW(H!Q1xtQ5XP{TgRK%-<0n8Y!HKccb!2^a4P z^;e)XM$N4tsE6UWHW)>)H(eF@4hY->0=iFqk0v*Jp_6Ere}OhjM~#qAdgK`xzOpmW zX)4f(t|1+sl<AjIU>gt^p>QfHko$X}iCoD49_E40Dt#YHE@D%2bQH$l$Zml)z7xJ0 z!g}c+Xp1c70Shjl0|T&!KG0O$gQ@u!xLjWgYQTI2o4Th$?*`OZSMUl5@IEy<0ps7< zDd6IVFazDfdTBlitbq%u4XM-hf~9g42(k1)TqhI=m!HK{Na^XAczUT$>lE=~Ino`f z6H*_Wg^YEqlWsIIJ_aZSGB6MZvOwrG#yhq7T(~eD!!lJBcobD#!1ND+@Mc!yHfX}> zn5Zh0{wegCNY+VzhX<1P*&R0MQF*2v1`=S9YDf$xg20=YNOUK85SQc8d$h|jt>Wrd z27@l8AB8?p!<wi9y>G~_fd2c?i)N$ZUU;M*M%5U9ep9T*#4><&!Xotg#cbwV1P>Lm zUfU7`mZH<?nQk^bQVtK?rwH|7O%Q`7_!#p>0t(dn*P#hosH<PB@ZQthvb9eeIb!VO z@l&3;`R*RBiQ+Gh3~Lbd@8Ozhbv~Zp{i(%MwwT+LUaxz{&Un@}e!R8mVyZR%Y>c<C zRr`?UY4=W;K61jO3C})0C1=90Dbt=AJIPxd_r`SzJ$jC~YgqsLhPsB`J8+l_NEYUr z{Nz)UCs@v&k9a-Fm94C-J5n(>VT|Gkq1ak>$6%YVzPn={e*e;EzE#twy%p2PZoZAi z!6jHY5n`Tt3c13?m=pjUHvBDS>6zC*4-_VzWa)ERx<9|c@`s2WD!n`GJ)8Ixx}rl& zL4GO|`h!tmXrl~^!u*Lu;PZE$Mqjdv*($^I&V>%D(94Rz=?;b?3~fYfJmWs=P~UwK z=t85C3`h8z)~fW@!_i*ZgfR9eQT$JII*M^coK7)1>@13%m<uSzw46(EE=C5$n1u2u zMn8F#;#aJ&JG+@TVwFkgCX5$~x8QD;VhpVJC|-bWO>r08wNl&_y^7*)xLc(dzNnzM z2QG^g_r(1w#hA^HQrt_3DvECr;#-Pu#jv3GHX%+UHpA0!c_Kk?AudpSJ7xlk_hAa7 z_ztXaD8_7mh2p;OC0rDq4rfw)Ct88x3|uuRz6&m=xF5PR#dqW4KyiO`cZxAUofO}L z0pHQM*AafN5Gf=WfaNR2193g1co3eL<e#~J{sp@aMT$q~?b@0pM1ggszP{dFIo<HJ z$F%Qz8;$C_*w_0=#xeWj^_^hSk)4H3TI;+Fy-s8~JDYHf#8g6d4q+iNm6T1lTD5hT z>g%bT9^XHfe`5aWy3COyvT2IVrX}w12eWB*@}0kA-RCMb2VnNgKQR9e`}D5+`X}2< zG4=1WCsWVy?S)yA(cF}8Pxf7Zt-k*6WVd!bRu&%jnydBoBwvNN;BNc;a(jXBER)ia z$&P&F^tg+Z&O+o3<k&3HTBMm|;4y8$K?X>i-(X;cGB5~P&bQD1Jb!P%%oC91*X4gr zHWq@c$1U>$4Uq+`Q+cn!&yrvi*_Tpn-*hikO6L%tZ|g5Wn+a0|)Go(W(BszjeuS#W zpHWtFy)x1IkAMlSyNlAmu{Kqj(5gcu-St;ZzYEC&x^S&BNas9&>6hA}n=l+6cd07u zfNGsSQo2=1-<Ni!j-;mL)l#m}JLn+i<EXBt$E}5LKuAo_#{>#_g%p)jbivOmz1m*H z%btu;CwyB__55A=Rn!zaFXDX6V%7MPMwz1b`|q&VGXGfqzWjt_U*|vS&HDU{`9Ig5 z!$sfYmRis+a4jY~ZDB_Vf%>8mH5F~HGi1`an~)z=aGevrWJYBGX2kq4OFeGcjU<)r zipn#RimAD*CelPDOHz?qlFCw|(ie6-Zk^Q&XKE^1wi+@?#fki1`M#>DF3v+qU;R*B zoU7tlD$Y@HHexGtVBA1C$Ye8?Lqz2rM;%ByQA#V6v;tC(Tk6QFXn7sg7`a}wJ{uTs z&Y(QYs8AY3WaR}&kY4|p$E~yEqq_bc5RlCv7pM-PiHs#Z)d93@sM4!#isf~H`w>6O zG+)z3<f|0<i^%11OC?EO%bkk+7HjI@c(Xu~Cr@eeiv#3!mXE0EWcT&Bwd`+@Y4VDR zZwcus@>(`m>DAIYGC;mYlTQcDfaiURe2OG5l_YsB&m>n+ds%l4i8l)+`5Z=md4Rmm zQu|wgyq0|dnI^B8_&y~)MPAF^R(iFQu}%n(uhZlwqjCZA7Zv$Ykb2xwNs`y{Tyh1- zH<jd9O7hu^d~txh&JxeaYnj2wD<;05ehrY<vTDfW2w8>YW`Mk)ITo{sZx!qW$UFRE zUmm0$w^WkkwR|g;3i9bzv!U_kYDqqekzW@eud@s!QmWf%*)56`O?rxnuLC2mWlc2s zLs;7e$j4~%hmam1AE(HdLF#c!B}rb(YsebNPqt2g{NE(`JVt&?fV|G~{4W9WS~dYP z-E9;T-*D1X-A2p$D!tlkSXKteJ2m-hNDq)tQRID)dfZY;lGpMWMShj_?yz|CLrFf9 zk^eM6UT4|(p8$C+TLzgXubBAek)9&2WzQ(R+OCbrch}^*HX@&<$j3?YQc056@(jqy z^M|bAnepa+Nq!U~Ulky)vwZb)fV`INflQNEOnmQ>o+7Vh#Y(UCF+Az;xbIc$_M%Qe zQwa9=&(UXyiaSFod?4v-G>fc&^d(GL`L3dXnh7=CqqXwWR70A@`vl5c4q1V#$XYYJ zU1}W)!a3~KS8BSN*F!3A?F!tUs0`6pSq3FDqiq2rxbD73*WLFqC+Q`?usV>a%6AWe zYmo=DdIc7y1u^Jla0IzRc6U|xMwA*>2j_l*4Z0$KgZtgaMoevm>}wBlp1ACeDWRWI z`q@f9Wz|@S8;<HR<sSEa%FP#1ghGl+OH|C`9-z{FaD0Jlsb`rKj7+vuE8NyhJRu6u z5ETK%t<?5vYSEdCRNJ3#6DY6(vVE?D<bv9t9PpQ1&uTiGRF2y7xF0MSdeFM>{?7Lt z6w_}ma2=%Ov|gBZ{lPz8+kz}HYWkxe>xKGcRLk$M-MFtPKyx3o{&9b6Q#Sf%s1RlP zq!W2no;q<qSKz83X$P2wiwb1R(VRe3DWMG}?XP^*qm!c$hzUJH@g+(1l0U0o_g~aa zhwne%lX+^C(a9~UU8T~VtAtE?+?mSjuc7mJ+|z*y9;m=FKzcm_CE@<b<4!@7dzMLE zEN?u^NKY#ry`Yq&5Ujl@Pq)i-Q`uh)p&TJ8fieq7@n(N>PzhA!=d_+`Ce80JN?sd- z7W{M7d*TXFZN#9j9p7JV#8IE~d(CWc+YjK0Zu#fDt^ZfvhOAzhwXE8;wHhnvmT39v z;v-lFSG&qcu=G3$s-!?_%qM}H1jo(&)x|Xd)vu}TihTjqhoSndM+)51B(;%S=1|S< zaod$|zDH(-<ovoq`SunmmmZYz?MmJbaudsNf{hY};o(>Hm?KzRD@{8psGCg7>cC7- zx4NWnA=C4hmQ^Ud+IF~=mc6xw>}pZ+*F^8BQKfVX;%Ubnwr+E!C24ATDIUL6kEueH zmEvPmyc+RZHi*&G1>G8;X*QIeq<KDv@r5o8o2;vkL~27<_j_xID)lQk$FJX59x3V7 zc-&gC{9J&ZmfgP@LhVK>lPbJ|6xCHr%N8iT+OP0%pxRaLpE<wU|8?zmHLIeRmoF%f zUN9FkOSL=CDu2jnI&~k1s(*5z$CbXqT6C~6<<+hlNLy4JDSuL<=fb|HgVehJNljYo zuhpvj*3uf(OkFBH$M&@g(#vSj(`qH8_M1<VQ&XCsxVrDW2~Ji6_>byhTG~{*<!^Bv z2FsrZ?c>kj0`=|NlvrI%i-~I25jo@=S=~Sd%6mco97fL1#5`YJoQ*62*#Ie{rEKju z^0D9j)?E*`GSCnM!mZH{$D5gds+6qz)w3vrdHjF4^4R~M+B;h9|KahxAie+L@v&O% zzc^<Cr3E=B8nbJ#bHb%8%Y=8%Cs!Mz$PF@^J||!5bMhg%uXvU`ITC2o90&sI;2bp+ zI-utSmcg@>^g2jAZmA<@1d^N7csC<~;!iPy6~xTjnM>E&SV-zCJ#L-*2UJW?6I%8) zWO^9tviCAQE!(E`+60~_c;pN$%U}QRfN4D~FNRE;)|at4q^CzXDo&iUm0oQ;NY8g2 zw9yJ<SN`5luO<2y4)jt^K@84ieLeyby^qal()*Zc@4ai9RrhG5`GdM&>w;B?UYXnh zy9Ey%)Ys!aXZ#-RtrjW?NS1zAPZZ~1Of<BEE0>t9M5P-V!Iej%mm%^XmE7-_ya$Qt zGo@UrETB5Jl$US6#<Qw|=bKPY7#YVjq+qz#Xepm*s^@vooCQ+4$@Ig)@sMf1uJHSC zkS~)dDxt3rQ!G->P9Z(@=u*pCC_SP!%X(x?N@`17vT%hdaLFazPw-szn3@7t6lt6S z`gUl|w{xrL8t*z|Z5b19PRF%l-rOj9)En)75`f$$S{d{ptUvza=j(k6uCyn+=zIri zTSoughLWR%9`Zt&QpkKc6=rm|8xB0SYBvB<N)5%oIgs5+v*j?#PXDJ2WH*!TltlUG zQ8U|<hhT~Civbb33#c}-{Bu-XwUO;VldFFpE-7T>Y3em7#jU!=P->2<br~=TQ<KOP z{nOQND#vdiRnFGQ$mYzY*9xi)l@KQrr2Sm<RI>l7vcFE*pL|lbs(SKRcN(SP-h$Zg zhMQ0_{iECCI!uy&s1KdO3{l(!6CVHaI1Mb(H^@pew8SQ6=EYOi3u6<^3hJh)ft-&& zKS72n$YvF6>K-Zu`yh~A>ONKKZxEExzh;TPZXz;6*&|Sr$Mq$hzNStUyE4h{2&F7I zQ=g!Ry1yptR3VhQ@1!JnL=h~3WnY3MC^L6p1f?K`5tIer`QHeV8L~bCUi7%Wqth?d zsUr9<Nl-o+{?}<d@Ktk>zcLrdgScCzd@*Pu)IA2xPE8m^lMeX;S6AV2RVm&z#Ct*w zl_qUYgaz$bDR@|!JP*f8K6wx{IM&pjB|B8=8z|M|I!mWcN&KvC_bN%ehEiHlA^thk zZ}S#a2qiuYC;o5VHHx=f>umYXq|&5KFG|zIrC==tvPV3s_<ROI198op>`<wv{iWIs zdXA}X$GTGndIyEvrZdou8fEYz!}+%UC9A)_(^P?!`d90^aUE<ASzX2@+KzgMjvHyS zE%L5?ys6DrXT3YVci2cc!#B^W8-K|5v-Qb@Qz1>i@veVjzG2(%O_|g>1h2YfOuoUk z(#oD(Vau_4Olj3>;89cu_qRSLov1eMZ;hEUG9<dv@=jT1%eDsQjI;e}ZO$3da_3>X zazy1Hv6b77d|DGqoVr@gr}naKv<6I_-J~7TJ+5*R#~!whOpU+&-by*^%);E@yY4fj zK5&4Z8I!v`uI;3F%jcBm0G&P}Pmhpaa>K2T({2xG{<-y+X$j-@93h>9q_YVIb%RL3 zN(f|w;36W2L>_4<4@$vvN~88EbXt#ly}G-|J%q~mr|w%St?FrU<{eZ_LwY;VT~I!K z66n}hgFo6&{>QB)DoHYZBR>Orgm*BVpU5^GrS6Mk!GKgMp+c3YL^3zg6DDf!GBKZ$ z^jl_VV`D@ARApUKh86}HI)@t6J^(u&*KX<KhT3HMSx4FMzBvk!FV#_D7`9Rwta=o~ zn!7xDDcw7`-^0j;LR6^AKm$A|!!uU$^sY}dB?aZ1%6|n_416e)_Ei#Ldoun>EPNk8 zS3imx1UVi>xlRY%|0R;q(CS*y+bWyh;jS+aP(L_IH7jtPh9~PkCkXC+VL&&!6r_?x z^yk|1C`&c-#lzIh*I$EWHS7EbHKtQIglMLsDr&~6wMNb8*z0|GP4_lAZ>@(w_cl3i zea<wbV5jxYjJP)W5YMl(kt1ud!U$BW#5ysf^Z1$nRA$-+QpQ21#zZ(XDhJ7L;Hit1 z!|MzLnwS(E)f%m+>EV%G*5sM-O{xD-Mv^o+=v_85HpI-AUi`grrRzw;JV-;r;d!7= z_uS8Lms#*YOn!{VB{haZLzgQBeISs`Vm+?k$;Q(V%-`*N^szC!{7XG4eL$%}BgFT- zH|6;R+ijV&=D7{E{?nR=Qpncr`~BTVEtrN=hMCsv7fS7O_u=1i=Ir=xIeSok6;V`U z%|y!iHnQKy7CQ_ZdKGc=K5NSC&Ryez^MqG4&hs40^NqE4_ERA}D!e`B{KXa$v&Sl& z`$XZ3?IONU!lL{_g0G`>{fqAKEFfRhKJT^{Cx<)r`-13sc)w#&dGF#8Yw!7P$({aq z&7Iz~m-B2LNXynX^p07Wl_o=R)<5d&ld-;h{)_th(TI0>PrPzpXmT&C;GKAC+aGZs z#Psh|I&3dwju5svjxfAt<Q-Y?ip?osH)<PYx94H#(hy5V81hklJ-y^)kBYZ0zuLuU ziE%N=T2s*6%J7U0S=ZLP&T}p_WP4lhq9tKATaESC!b_3Aqxy?5GuSU!uUiK~CM8%; zE-ehXoM0_l($b1qHa8@uowaCLQb<BO%Zn(ro%hRSx7$J<fD{TT%XizNta;0u7p}l8 zV@0fJZ!5G?R<t)?>LOkVHHt#TvQPu(74r9ZsNuDX-$IQ7yI5eS^UZcj{X+hpup9FY zaYde&8paMod|*)Ock=hVVR*vC%V8wh7G~@V6Xo*!{cw^T2sgfRh*RO|?>fXLhw;@- z4pHf#1a}1ez7b(;j1=C8^m&oul007<N%Ak~*Q#65$^r1VC<0Aawj%PsHlg%_rpC%> zv94+QPff)+d0rSz@;&nRn`ompM*Jr_{nHrnaSWZGjG@$J&FJ^7X2wUcVpp^DC9z^* zES+zQCHWcodpXv~Zy~Ja=_i_tubb2POD!n%UoDKoW=ru+3*&MN@e4ppaiHaV8&bY# zWfZm+ueA=r`RDTYjX2{<oY)#?{4-AMlILH<|4OP~w>7>@5cO@1<q2YOg7Ha$C`&L7 z-6-a_^VtwxYj140QT#R0_#shjNi_B*icjSE*+i1BcbXJ!cN$AOh@(#9f>WH4=bjEE z|3v<N-ND$FB!2E-EKCyflj!`dBuZV=@fbyUosC~Qi;~U>2Rn;B^8EYGBwu|a(b#t* z)inR6nfUdj8da%cb!y}Xsp1`Z{#j~kO1_+GZ0aJ`bulW@x?M0t&Udk<7PYf}eSKN$ z3l8x~I7H{djVs~e*Kh-Z3dfz^HE%2rEiAEzT|0?aYc5-nCFy2!TpwxK+Fr#ssrXhE z-=$)@RU*AHDxRp~SPWQ_FIVw;6~Ck6f2jD7iu+(Fk^V`4Y>J1#lLSwy*p3H-gyU3v zqlzC>@oW{->t&=z@2ycx@3&D*@4-<_@2k;&Zs@Ii{2xxfBKqsK`YEr@^S{?vV%afW zKq>QO|J$*J*6=ln{c}+-yd{G7EQP>KAvXWR?x6D`#CUl?zV3$iCFCbRSUXMmc8A@8 z=POgdUeqD(x!vxtJ}c^F`PM|4(@^$H2!#l15w;=hM5shKiSRQ*2<}m15t0ykBJ@KT zfiMnX8p2Bmg$QdAwju09s6;r4@H0XPX3kiIB!r#_{SZbVjN66xsHfrNC4@qRwFuh~ zb|O?FoJ9B;Aq1=MScD{mo(TOAMj(tsn1=8YLLtIhglz~r5h@W*BK(XHvIqGQk`Q_# z^g|edFm6u<4B+G?ghGV12-^^LB2*%rMEDsYWG@ULBq8)f=!Y-@VccGC(z>uv+it7# zhQtVJrcAVXrZpUhTCW7{P-s05L{9ICcjPu%Yc?br)Fj_pdp2|oqed%Aw|>}=6hZB{ z8L!RdTJdithEb<Ty~|2}GbwU3PDkH`hBd5-Z?>1;DdyaRu465EGckhRE1rY<_5{ne zp}BSN&8P_aiev>OEvyr;BR^Qov%*S~!cw9PN6M>K*V3ejsCdH>wGwX}S`U~0Xd7Y; z+?Z(43&6dsXExptmbMPLqpa6ACfU-g9UE^jo?S0QQ|r@>&8;gNW5Pyn62jTuirthH zHm3qzxwm!urYB6QJ$e{kv80ZMSHMw!ibL@hCmrMgibL_%CmrMiibK`6Y2*n%#=}}1 z<PJZkrEK8al~5WHILIp$hsrU8&;0b46QOt=m=68MbtvArri1)KaVWm7p@Td_ai|<I zc*_v(@nafI@Nc{|?XAhSUe<}Xri7I2wMK028Y1>tFKzB>I^aH1v0)G{LLBGRxTSal zI8zlB0xIv4qB`vOAYcd1=}{8R05-fTeFQN*G}8DxV0v6c1yK4G#8Zyj9dv~!d#N@& zrlIFbeuZd}j-01(C+Rm3Ee7D^G{g4;`+$)j2BXDT;0p{-2CidxHt;VD7XV*mc%==U z>k1RR2?Y@nD1dH7><pI!I~YC)9L4b0z%dM;1#ZFcFTioYy21wCc{!P&C2$JEoq)SD z1GfS9VwkQNX$;eqqYuLm0cS8g7Pvpd=41#4Fv0V{Ll|BFoXKz@@JNQ~!=6#VsF-XL z`o!U3hIax#!tgQRM+4Xt=OGvkfo{`lz*!7mhiln*hT8*YGu#z;62s}hIl$V*4+76( zcnWZC1IF9`c@X3^2%<#^@M4BP0$#yzC2$GDr-4fu{vCKL!;z?28N*KCa*a*&;xq^< zv;a60_z1&~16MIT8@PtyCBQz0Hvrc$Tn2oR;r+lO)bAmD{S+;}gTSE#gntE&VK@@g zU>w7Vz)ps{0;e$C7q~maLx9s5rhAMGKQ?8jnFPTArZ5{gli>p3Q4GHUoW<});B1D= zfO8n$2Rw`6W5Br*WBgN%{|!MNQ}_*dF~bHX{uK<<dZ~or9>Apx_Xpm}@FT!w3{M0u zCyeqX0a|EPFolJ{M;Lw`xQgLzz%>l-0roL`9Jr3*-+?bO98Fg)%)F)^CLJJfu&MYq z;24I70>?2t8raG3MBo&Lp9fB3cp-2G!>fP?*bTE`oNk98lPMen9>wr4z*!8(f_XN> z-GOr$9tJ#%;Vj@>hI4@P>}Fu1AXhAA3X6bOFuWSLgyDC9OBwzIcq_w4fy)>^16<B9 zeRW)6GC?F(3r84^1FmBD2H+Zo(|~;p4+gGd_;KKi3_lAjj6j78fK3Mztb!nh;SYi1 z7(N2*WY`Cs!thn#?hMCbxs%3lN8k*GZ;r(Ee*hD>ppeP%7~oM1&j8M1xBxht;WfZH z410lRF}w#jm*FoXG5+$H;Cm=6X1E@B1;epeu$3^}5xA7$-oRTK9t>Q@@Mz$2hNqek zR4~CSz(*J^2CicGUEmsqcLDnt{v5cD;j_RO8I~)|us|O$<8a&NV1gTfV;JrY9LMkw zU?;<)fm0Zs3f!GxH*gxmZvbZmuqn1eFd!fh`++kV{swpy!@mJ%F&u>@ayG*^0_QM% z5AZC89|z8Dz_|X;gdnd$5G`H?Ud-@cfmbm67H|o}dw@$BJ`22+;Yh66%NXtkT&^)) z{|7)&p#{KMz(*LK16;-MGT<79w*mVYJ_uaL@LAxC3`gRofgVcWGcSSlXA^Nf-9_NJ zKCKJrYWR@<t^x?x2(0sIT~7~2{cF57azhb%BC0X1>s#Ohz>$c_9xda0;+90O2?yZh zW`=1Q-wPP!`LF-9j!*LoBmr8-_hFdU@fi%$I=(-{D}e_vOw0Ho4AU|`li_mUkrJ!x zKds|OF$G%3k7k(G@mUPhIzAg1Q@A8-V9}bxFfHR}F-*(&T$Z1f@e2s+>p!jI^Oypy z<6mW%*71uOrgi*MhG`kUf?-<57corB_!5R`8NVL5t8Vj&X!C5eA$1DfCrVj{*N{G_ z-rD)`WAb05&=R(gXc2<v!hEBGuE6}0?NVU=k$NRC-#el?ZczD73JZag+1Ig)fcee} zrOgCA5b_`Ni0CvF4qdUf{r$0|hAkb?1$Q7V@QrB`aTYjduhnyVddR4MThq38HT69; zeZ{;LHu*=!<-q(y^bkCxI({bb4e)S!Xr&G!F&>z2udpAOe=&Osn15q@rG*H^?gLiR zj;>~4CyXYdDL$Oj2id(InD5N730QyJOkdW12u%OK)!!;t0rOoJ&H?kUt$zdNpB6vU z8Us*oL_r@TFU3<hzLUmT;3)Pr_F(-#0FWbs<AU(2AnXmor-1n`4x{2R0Ql$JOYtQt z{~|vE-;?t_2HFAh@70rl`A!=qus)nPc@3Cvs^A6Y-zrC$ZPC;n1Dl<61l9-Da5ymk z+<PK0|6+evkp9^qOaqF4eSIe||6tzCr4gqNDsV>-{tTFJesUka%H=z8<Puf~<u3x} z8)H-fW1EMDT^{bnSG;Vm6I0~UXT9u$^^nBIJ8859)&~{PEeK}-^L;0>fcYK<F97pB z3>G%ZZwm22BY_Z?f%(Q2qnx^Uf2Vso2rmKVAG~h^<{M=g9WeT?2KFGh3%Fran7F<T zfIuHq!)JrAHwb?L%(n^&OTzS#9N5F4Bd|Wm-W|Yv3xNlK^)?Y^5{BJYmcfY(yMg&$ zBeg7@4E_PkcW1e!W8;l2h5+*&c<u+*2j!m%to~JK3bI`gr0@nX-wx(5FyCk5)?^U4 z7TDP0Jz&1=!>_>lpaR1?HQu)(8Q8x$3fiAM@HzwnF0>({T^F_xRtFh4NEz6!3ul7T z@4jAZ`}@c;V7?d0hM@FvV7V6w%9jMLLcljUy6pxmqEN#}kWMq>y}<NtO&?$)I38FZ zL~s@`-?Qi=V7}GI0bstd#wnGb%9m4VJp}Y_j>a)5jdxAx4SW$36<x52z#w408_RQq z)j@=o0Mq_57S5>x`-0L>03Sbve_LR*Ne{H>9C$KKyDGH-<~!!(181>bzLC<^K?PO- z^L?g%0xn|PNZfa$*7jF83wT?d8Q9`yHw1k9mOp^`#u&YCYHVO|5S{?+1VJa%$V0P$ z`Box70`qOA<dz=N|Hwc!NJzzxItV`!gy#eEtv-GUO24&B;|+97V7{x<HekL3PfAx> z0J8%7L7)#JFgys41!h}!iHe|fqg&&Ru7(2h4Rbue{#}>+O@6K$%GV1D+AZZ*D6Bda z*aj-GJBA(pYeb}o(38OWAp9yY-vsF+V7@`lWnjMlOvfILH$$1!qoMx(0!%3Ag9<1H z=3Dq|2ulBZ5IzmeH^7PRiR(XxF&1>PlMe>g2ieazg9OFEd>^VUz}s-~$b>RYREan1 z?0&+Nf%#@A`+@oHEM0ox@_pLguou}>V7_0CSposujz!b~>w_5Du%CrK2)6^~8`?|& z=39U51?D?(#odalAKRB|25=GnspFqVqD`>?0)0>cB|-Q@V7`UVzM%B8z<d*=8MkR` zeh+N{E@AV=zkmbV%XDlax}-JE@G>ypPwuR@NyYwF?f!3jxVXHP$7_v0?+!CO#p|7T xrZ&X7{Zi9bJ^st~Fs_;2;l63%5e?hTSh?Mv@ut+xj%v~%qAh8>tFPGX{|DsoOhEtu delta 21519 zcma)E3tUx2x1R%N!vRt70pbHWAo-RsEwz-$PyG1u5t%8K8d{llO|i@@t&NXk5*}Y$ z%~A`ql+21cl~h`&sF+z~Wmuv3K4ck|nQ!iY&6(jWj`#lV-8yIXf6baT|5@|cvp0Ls z)!FufS@y-TVJSQ9jvhjYsY?<b3i)ca5MQlc@?y)74#7f9G?r`}5?=a(5T&7pqcqeC zvQ2C8PVdlj&&G{=ykORu+^2TjG41|8n`h5K{tRo0ZIms_`oorB3g@!0ps`_1bM1F$ zXCUoPgn<Y{5bi>tBNO3ngrNw-5bi~wW4JGpza#K_pThSe9<9#Ds`w$q536$>jTgc) zgQ9Pq2*Kk@I7!7hh@V85sxU=QAv}%nECK7=pr+<5I%CIdoKp?wDwzxMJcK-i7ZLIi z+z5*hmLkw$`697g#T2bTSc&j50{kr2AiRq38p2wH*AZ}-^7jq->;HKZryCI7LfDA# zHo_(ZeY~SiiV%AcK2-Q4#2+JUMkq%31mRPJ&k(jEe2(x10v+2BN=0`ecH(3g!fu2z z1UmL0e2K6Z;VXoF2>STOinIS=MHyB0o>h<hY+W$Mk9jprCkoL7rwbd@#~mH@@zWvo z@w|Za&H?!+SldG4%~|&PxhjL|<KF|)7uxElzY$s=Zw@GRVnF__^|AF}Xk#<FzD#5p z^=<A8@WkAJnw|{MpB$k7SU{a-*4G~@ZVMot9Z*2O0DLk)|HFXve%2UAeB&Df96ruk z;)pk22#|U)P3VJ~xjbAWU$Y+xs7y_O!(R)i+_V7w{`K{nif01|d>(G~3rjY?4&X67 z!0wav?T9x6^xFmSo)v&o0*bAzZ{HLP12Vi8;ABtz4B~+R{a*rl=Dz}Z)t-Px>R+!H zG!^3AfEwOxEeVe|TLk#?$pBhc19G(u&>tC){udhgr89O1xg$|D6szNeaMsuzARzH9 zrx0iUusgc>uqj4H2~qN^-7!-s428k8wuWN?aAVO@EN>;m^wV|+@gjrsBQQ~(usfP5 zY_=ET*0cTs{?!8W*fqO@+|yV@igB%lxBz=T+)|v1#FP$#i88<0Q0#zXXCJdWsCQBU zS5V;UOMZjl(Lxkntz#?>B0Z(VF8$O*G!XfXQ6r?2zbXHSB<RCl66_N8I&n2Xjgd`a znuu>Ysw!kcAv;Eh)GKyJL#5D?G?NSmH8^?T8{jC^h%RtVpuYzmc!#+-4H@Sv4^a9{ z_&OIJ(@n510W&-bj92y_2H{N??T%bvlQN{@WZ_A>15G6v79;1%e|!p*{wfNnf(xn2 zMCd!XFbxE?OUK@TYajAckEI%$jf6OV-0tY2@Yle1GD5rCpidybldffC@J_VZEpVA; z7z*r!M>L`NXo6Tq=splSdD-qjxw60=IO&4|v<FY3Ptmdk$<kvpFt&!m1KJ}w@W@Z- zl&N$rCB`vrg}A~xNn0@d5^V;tOz(%zaW4w+b()6a6A*eG>6&1Cb6lxe0Y9PBoo6m@ z2*p@5k?zFf;K5<P*&W}Qs*{Y46(WXtU?(zsh*6;#&%?0G-L7jS#6PIPMvMaOf#`NZ zWP=bQN%)PHSj51C_i7A6D=*j`mh#9!^pRyK&(!16iG(6m0M+sv)Tkxv)uYkr4zL>a zLi(giyW>t}@Co4k=mV&jB+$SK53x@AH#}AhT&U8wqpg3sYB%d9h7G7{2Q&fFMU*&O z`!pEGf^jY|b;9vr+Mjj0e`8d5S*J_}p2|i+7H|eUr29mBwB?uZNG9x2llMiZ%sONG z+weFD6r)q<Hco+yvtZCr84kiy4s-c$keld)St^}edn-Jo^>?AuJqFxer60nuOJW}R z2?Z@P;bPqc0y9|wYS;%EXjCLYAq5$9Cwl^odlW9-iAxjJU=#Yx5D?VE?gWh**rUZO zrMp33GzjQEl!GaA5ImsE%R-Zx!%-vj6`7$Y4A*S)ce3^9oG~a6o|Ng&qkx-$hbsIY z2(CQuZ=y)#{|e*2qe}lTl)Qya#U9{^H~5=4%)t0d+2Jp6CaRjvJkS$uv>FCr4}GAK z*n=8Gpay#Ge+bNr*i=0Qdbgv-x`qWPkoT#>sEHE<w1<CcgYiF@wMkbLC{Q4^A$7XV z821%uV-k}~S7B~1g3C`~>ZA0dm}rtzr*(=Jq9w+oQ*}bpZ;gyKtdkDL@W@u_l>gr* zTs9ozP*0U3P_@-e;cg6*QdYw?xK=E{L{z2-jYAV8vQD`QH9o%6?sySxO9i}uiL3|g z=}!4O2!DZj0&22C<|Ql}QcAUph4`y2)@U%8t`rublU1=MnTXz&vfZ!$5%k)#QSc=^ zbU#MXD3w1E`LkK4i$^Cf^iOOfM0dDkBkQFrL0}U)nV!k6p)o7qfxDE!D(3PixV$gs zi5@6W>%WI#+U25WMblL_&ppwzY}T_EJ9wg+CkI(~9dGWL(PByi&r^qo+pMgUnbw3u z4Lv38rbb)c`*aP%Xu%OAf~?2;475J&lWX<qJKHMh+uF+Q8)CcLI?}gk!z!?+!-ju7 zEL&>RhQBgNf|US?UnfY6vF=E18uur3z*<C2m;iFYf@^fQrldA)cTH(b?YI=3JxGLg zkj04O35FvC^9C)5j6v4M)ZOMf%u0!dgD#OHDaH`UrWjrOVTyMOF`i<$Wg^8GQaKc( zBTS_j9cnto&j~S;V$5>66r&@}rFae&6BMHv7Ep}#Tud=~_A<m)&Fu;1MzkTxI-#j3 z?u=1EF?#Wb6nDi;N^vsAKgBm<iAQk?Iy1%HaMwujO_(nz?v6f3aStKBr}$<xImNdK z@e9R0(FZB+h0cxG4Bd{op9H;y_>JOQG0Z5wO^9<8V}M_vxGyg26sN+SaAD}}Lf9!z z6T(699hgrj?uW*q_+LUaqqsk=JrrXUw5E6fhBd|i#xfx9MEiT@{ATx33q3Gv$EGp3 z!MCo|*4DcAPc^)4FlTxv82D+H>KT5=0egJ6-$&;inRC1*bJ)=FG`EbWwcOpq$I~$N z-gwzsG@#HNAB?93d*|F{pPIb8f0Dfn8oTXD)au?#muqWDMuRlZp5(ptT5auTRCaY& zjB&SX*45ftl7E0W|6lewCH8#bUM{7>k{o%+>2}RjIt!84pQA{k*+?_Vz!TblgA9;3 zug<`9Wnc)hoNkl%UEVIgna3c@tI7L`Y`hDyZkNmplqn0Cqw+S<_9Mu?lxq8?YpPN@ zgZOkSUjf=oxGJF9iOJ9H()Rv>s@s=QR&uT~@xVpD39Y-E(!jC$BP>JQF0Gm)>CU}s z`dm1c6|Pl==$yllliJ~~TD!yTny(5wtXgNJlrC1%J<_h!k<_&Otdwi?-6ar>rMjAK zmll2jAu&B2<uBxIQdCaS1;3>9s<+~bfxd#FA-qRX^}HQ<71R_VbUrFqHNK>=QqjBZ z0_;WS9mw0A*E-2N>SC=~n|C4a&*~>(%I%U`(9ciLB|8s7SqXveKu6b9w7IU3N#}O= zQNfHayjhG&0cNwjQFGlc*^MNX)r!jd{}7e?G?(vCn50sTS>NrFT9V3KqOu%zln`h$ zoT;g3bL$|JRNg}VgWiKx*2X*AOgMO^a!m#bq&Q8*y;Pinn3{E_Y`_8_D(O$KC+Uk) zTBxMIK<ai$9a#x2w^Lgn*8!{Xz!v6m%CkJ1@}Pf{l~*7^dhNGvm(KDIHGyn<zwBkI zylQ$)#6^0l>9uT%(yM+6s}n!>A|K0M-o$$3ixv3@NnR>R@><>lav~oh$-gbhXEE|c ze)2j?B{iBNuVs56)8rKs??<Gk$ZOeJrB}Thmo-265>0+6D(ClnnIhi@QnyPgN%C6$ zFu4NcZ?TTR^B+j^nT-5qKY5*{AtSG47tj02D<<9(kZJN-_N~&Z{s)(HKlyS^el~3T z$yX@y6CriGq>?1B<;%zwAV1D}e^3ka6G=XUk>Bbkue0<eQp)pM){&7{OuW&Iyq4KC z`S$h5S84JeAl*;itH`g0)a{Z=lDw9ek~NTDVck8Xh53afpT@}V@RQeB9{<}<Ud!%- zOm`c_#G6ifs@rH;H>FoS0?S3e=W8_iqe%CYzo5u}1*zL5l_YsBKS$O;ewUSS*Z)CY zXIXR3PhQLNAk*X(6Yo=`r^svB7^PSJ0&WNV<ON;VFmrpG)FbckiM>~l%k7d%lDw8D zDe~v7vmn1ydVW0f{MUZZ>n!{J@{`xHPaxCe6%+59q^HPh*%GB!y%}p)Klvz4z8`4% zJs+#ccah|!k|eL?!^jozd^<`0M@c?~kw4-mud`hJ(@$Q@sv*<l6%+3dq^HPh**2wD zeJX;JcWUxapmKil9TfSIkh)z`Ns`y{x#S9vA8u{SY++VQ^1T@O3x4uCOIJo-%UUq< ziiy`?<hAVVSzjfqn}fXD)l;zxm+X2@CfMJ1fiZ-EHA#w6_*E)s^c7hF>DicX^U@di z)jUwsbBR`djcQ1jpg#Wc_Cl7QzQAfcG_Fe&3c@+;)$Z56*1Ymn-s-z?bEz^!#<2`a zX5PqTSh;*i%Z0v}ll02vgFir`BG1(Yu0@`FVU*$CEI$ej1xJu8WOr9}uScnFQIkxy zrZvqn+^^O*Vrna7Uwe@A#4=&jV)|K0KLzwtNIymNQ|xx#uKau(%1}sgxmmVDFO}X3 zZJM7x*S%Z{h9%jl89tAC6ok+cC2p5iJO)F$CA91tssOb_OMiivA={l^N<OGw;(){C zI$YDKR5_~G;aOq+;8N@JdpZs%6;r$Cr<c-gfWDzv`@65t-v+Y8Xv!HM);AdM*J!tw z;$9{n?Okeh9oF2Gt^O@4M41NQc2J#ez@1}$dI?FF!F1)QvV7TeG_Hs$C2YgQSGlGT z3&rcxTu6gwM<z9^bdd(mA`o!flBOJa6&=EDx*i71zqbUdD(T6UC=B;uZkGe?=w2?B zvEXtqCq1q7<Zntz3c;iwa%o>{Su$kO*B$(2wj{moJ`Ym?c_uQ4Xhr3EP4Em#Ap_Ot zLg9<+YJUPVXQk00#BdzigAR4Q!tr_vg_Y?Al~|Y;qhTu76p>5P3oA`X9z9KxVkwa- z!%5&G0W7a6QL4~<PUbE1Yj%O=FK#JtNt4pwGDn@ik3*eWYcVprT@uOpZ$3Gc>R*BS z`~2;@uiT^>s6@)8n#`c(<KT|UQRQxz7VHKK?R_otKqkFUcRb4AYow><awQA$=35CP zoMyRF*&o2=oRU++VYNj_a?$b<Jm{(%RpEANK?38VWeph@UC_l}{anhgsawSp5Xu&R z8hzaDipA|e>zb#KTGuswH2Oe;d^;0e#iw6io+#;6xm{Xu4i&4GdRjL5*V@cFUs0K4 z@d3t9%l@tD9mIo$%Ji~IBL%(h;GUnV4_Cb!`BKT;lE}HwVRToz##@WBoMtbi%5eu( z3_VO0^vuuuNHx0GsItoRDoC4E8fn+5(bH6KYJghCb!yUDbFEh4vzA$>X6jPu;jg#G zt4dXRS#^3^t$@___|y)18F6*b`ZbzS4Y-SyYiNB@>5{)|%F)HIAAb)G!&Ty!cPO!P z4PDVI)61{7`V1AgCYxp(Ip>h`dt$Ej-CX)*H%b{ThO5VtkA3d9qVA6|&=A=nR_FU$ znB%Xjl&pLCDU`uH{=c~Lt^c3eY_0Zxc>I9?z5m7IJ+#{Y;GABR7T}!oH2meAx8FHu zq%6&ZcTOi&8ad<!nN6Q(ko%RKd`RwFlO<2e{cV~F0bO8xx9!<#D3qdw{L8uVO8PFO zZkN=NGXlv?YP<tTptz7`?WkGXbLo251d`hQSYxQ%Yf&*h4z<h;nVzzA+0QUNEqg@k zwWMbem1}6-;+t4y`H!CTo7QvhJ&<YB^11_lx09Zl4YjPZ(yQ(U(sR;FZFH62k+*Bx zwRqohr-&*Og>zY-Pk==412Y@;K498=?|8!cVq~~EMBT0q#!5i14z8V`nXFV_KO+Ij zqxY%B!eNE=Vj&(b^tVHf@>y@x+ush;m9#CSZkN=N?O<(vz+n!A&>JIDv|b117A<SW z+@ebnOi#;x13&4$X!(p(Js*S*|3VSE8TG@)!>4NRuXFZvv<-EJy-ZQdHbbTxML$bj zM|$c}sFp2KdSs`m)p&FV^EfUKxa#Dm%ca%=u&W$Zm7g9#8qWj$GPvrO=PKxGUv-?) zYAVxf_Uynj<`Kj3ka^;q=Nz_37rp8+>$!;9S|K7`w*ZiDMT|ku{YKIReC{axZBH6S z=Na%yG5x6+N)EGxo+Lw?lBmpq#K6c5#>JqjGA<iZN)5)JGmzg&)9OKJ6w_Y;DZkE6 zDU>||iBKGb^|<#iJ+sEc#!6$lZ!U|iG;)1sa;Z*l{+MJ03wT+|o2n00z(`)jx-{5? zt%+o7S)Hx<J_}vsw4Q|Q&Rq1|N<+zGrCi!4Iho(2%<oX<JHi=gY58<>c#?Mp{7QEM zME}8KvSj)*vOE17lI%k&ox%=m*}sGhcX}D6c)(JdN%G#zth6|0bsHOF7R<8=(nU5@ zL(bj{kB~K){T|95E!U4eK|iWBE$Hu3MNfgCnEv`nmC%kRGgSIAFm<Q@NT)3%#UC}r z3Q4g%FA+(#NKzAnuR3j1A(U_IQAV>|(JX*{?~hevyquI)lLbXn3i7F<il!|7QwXl3 zNoFW}C4B8pKP-!rj1Oy?ze}1Gc?r<UU^J!8M^zz|?`<n(Y;%)3i#xqSv8*DNPRUX} z%RT|ynxzzcPYtP9%F=`XCzfP}ifAWUo|2`hQ`O!VB+Dx4%q3KaZ<h7ds7Mt;`OZ~I zXZmPXDVlQeGvP3~po)}|8BHmO(=<D|(`E64Ah^z%WQMX&M(@!Cq;?xdU@W{b$NFb% zoNc}p{ZPCu*K_kj!)&(OJTE@n$Y$GOz5Gb;kg2GRcd%9Q$k(<*)`oG%f-e2!S@7r_ z!{+itJ>ERXHrVPgxxMXPD|2#L*OA|1o=-nVHBpn(6Szcpjql-D+*W#j`_AWKA{J$( z=L`$F{H?VrXSr>d)osdH+cs<2l%dh@enS_Fh`e%JiLHEFRWK=sTbHNwu(_<1CuU@a zA>Hk}HN5t9ZSAh~DzPj52*Oc>V`%Wblp}b*oKzNL<a>WCM~-{<(vx$d=1$*28b_d^ zltaJ^eV4|}$`Gr;led~*Ad&h)^(pirT6L#?<*SvPWS5XZNML0+6m>6>f^?PhbVD^7 zdXR>1&Fvc=ZME$Q)IVKssC$fwsEnG#kC<>IvBCYgkd*7RLaFFFbP#B6$<&E7&1^<h zF=xoX@12D$lIbgIH=jNG4@_qtl?OlTX9!B=SE%HTM9tm}uyF=uRo{V$l?<O_hFtXx zO;d)N;O%Yl(6|6Q>#f(OwiwnLtq4OAEue*1zL!RZQDHP;DP^31jINT%WprcV`Vbw7 z^jj%?tgU@lEWdtlE8?l-aqUP!x#8V?_1$pBS7@DE4iUB|^?(iUQt0Z3MIopv9u}n^ z_j`IBlF`oU&s)(}!BYuoJN8mH_}<qsj)QXTcLYJb8isY7OF;z$vOk?hS*n%au<&Ag zyDwo_&35mh!gQ+Uk^ZO&O<%CL+)A6)cJLUos^(QG&Qw;ZugMw9WExU1No(}?RmE70 z^2y@OR?)Or+eB;6w2tQ3FMZ@rxANC_2xMxYctA@RD@V@j5XgxV#g>q<B}_vKa<xW` zYRkzGxF0z;r&K&J`oU2<e3vlU%J+VaH3(TFLkIVIo_o4UP^)z5*?;PLqJ}g&TN|HA zNWX}5O-YKq5J*Z*-0A<2Uh}UowUs&_7_AGFQX8bYtaMM*v#o8nji>d^Uwbfx>h+Dj zY*GrDI$!4NMQU+0fijHURa?8qx;nPG_52LC?LSt;%+@!K+zI0q`S(WYAvJI}ShkdX zMmE<Bl%iJ+wY#hlGdm{R1M}<%$kUeP*<igpb3)MXyFD$Qd(am2OR4qj^W%aSFJ5f@ z{Cr~AzP2<`wf2TuA*L%d4amD`muIOtDKu8Ug4cV^(ANhXUV`@(aHK6N@i<?YZ%d$5 zTdUyPreijmgaqg3wY5oj*I?rhwY4J=pY!aQe^+o=AN*~`i3h5Ktb_$Ef>O47`Yd?C z=9F*kw2H9X=VN@*I7&h|S6o|5?+V$IVy$EDuC@VI<Au|M4z~2<FZ?w)=v+(B(@R5a zwpi<V>r!|J^s~_z_4Xud_OiV}Z?v-RTE5iQ-8#D53d(40J-s3^Xjp4&6{3e*d%jq4 zt1Zag+PeJGL4y|WhPBq+)S$moJFK%S+t`Y%sF&NA`N`siU}JT#SRQQPyiEQc3N}1; zaUSP(G1pG#o9vYOgZw>WH|7}PiaakgjID<F*r3$I^7pi1xI@Ie5R!ZlV(bnPCGz~E zP?GEoH4Zw&vC!1_9b%)yIA}V=eg`GE!svHZnDJJ)@Pwt}<?c)J{Pl2>{~&))hZ`3g zh&c^XzilACY(VFKHlWluBk1?D2xCt}@m)mf+J@qlhIGEYA<6%izxj=fmm|fxMyY2T ziC^XU(nzExiJcUQpCgUxDDhWh>b5BHX%w9wjUxH-#`v`|Ux{%5oM=K|(kn6KwU3%n z`rc;7!Dw-;S?c@IVq-L&?~f+=ni%7?=Hl%bV{eQomFKIP&#@ul%~<1TtT+@Kg!6SR z==Zl4#=e$fR!d`5OR=maov&}Thg54@8&~4Q8*#>#IPp=OaXd~`#u;C97KLrh4{eD5 zh&Nv9EM9OLhn?azr?J^7isX5PljO@2==Yrj!<8s@CK#0o;($EACgpfp2&sD$jdg9s zu|(r?BE7YP^VMxhzBK6oMHf37Cpw8Gomy}0BtGs$=Z89xd~u)7#G|+~)$~&5XYf0> zt8u8SSlu=J<F4X8dH!QpN}Zc*c#_4YWaH;#@qIE4kYuZ1Rcq_hRV$l+<`7?mLiB5> zaV1pz9cn=Et;6t4Sp8D4RkWsTX@NZ?2JfZiUb4aqQq8j%H25#BsF+?CBiu~I@ha}5 z;@ec5q2d}9kHOHTe3MlCjEd*0c%_O%FwjVElP@+!0=SZ(n~DonyhX)dt2hV`5=k#! z#obk$fWOL+yoZYEuQY_mAf`X7&@1hD3@zWc{QFw%<OQ>QZw+Qxc1H+aa~B!&e6N~D zTX(OG@1KqO;B^nYs44{J4^h6=?x6D$#CR=1zVC*2Fa#}&uxgqtAjfln$zUmJ5Rd*0 zFC(vA)5G$vjWDO7?0E<)5Y{7njIayg0K!Rx^9UihC5=I7hj0tR0EFQP;}E7H%tKg# zupZ%Kgk1;+5KbbTM+m{(8H3Oc;TD7e2*VM^ZO2Q<({M5mVFki^gpUz+Asj$BiEthv z1h;B22<;GVK^TBA9AO;7G=zBwD-hNre2lOQ;Q+!(g!2d?c!4Vhp&i042m=s?BaACe zfdQP%Ls)^Z9^qqzT?hvdP9mH~2-yh(2<>)yZdw--j75L%_3>fUEE#BHVT~b#*2;k2 zqpdt3a`I@r6IWopzdqieCJ@#S>)VD<BbD^EuB}fDqjoFDYjW9E#|`lz)Cp43tRWi` z!?SUkorZ?BThDH2BfkPnNyj^d*4hp6Vf0e&>gVxtoYi=JQ|p%v5n=SnMKL51)>+t* zp98*N#S|unq%|-cY57)aVPaTPtl>ynDMX<)q3{n|AM3ui;thIxHpzPKt@a@u*5Sp? zAZz1WiM9^bzPH*N54|BoBkS9@npzPXqe2F41eZ7~absdg>oOtc_OkBWIL@TnqwaXA zk~$WBWT-K@1n)n?0~7}fJmaH-e1IcZc1&M)^<g};#X;`yVOq8NMuiV!1mGaA_^|v_ z<UMtt83f~HTRQaX%E5SDnGW&`#liR#g%0uz#ldof;MxxN_%MwgA#fd@{q`hVl6Cg& z$wApWt+AVugGzQ<OE)E(PPmU$Y#4Nl5XU+-juxwcb5v15pz=N`s>6=A>9^vX9u?7u zx8dFLa>VpZNaMr6^fZVHp!6$<i!0H4F#V9dE;hWXO^=Oy3XvieIfHH5K%~e3wljPW zumc$RVK7pR0ghmJ5^xm5Gk}{hoDUq!@XIz6-Hi!0K*7l{-E|}}TmsyI;je(ZF#Hp6 z3d5&>doX+sxEHXluwi1xOJjm);0%Uu0M29vdIAq)nC{227^a74BN-k6oXzkU;D;GD zCqXct37!F-$nad?9EO(yPi1%w@N{5QOg71Tz|Sze9r#&>4*<{bV^f@lU?v2*O|Jpx zGTZ>yuDJ}i0iMrrGVlV1Q-K!)YZng(Udixe;DS1AMhdz}7S;(OMFDUT!=C^bGrS+T zgyG}BWeooVT+VPfYF5Fp6Szua6TP?>1YRuw&IGPu_#xm649@@-LH-&p0d_FF9yp5O zV&GVY_W(N;#-~D&;xGgqn8M${DGY~W8tla|W_^*yFfD>I814(4$?zcHEQaZ>BHM>e z*=Z(1FrF#Q0M21JA9y;$tAKMEehYX$!^Oah8Qu-NlHmiu1rlTYlLzQ)*g~do9=M2M z0~3ES!?D0640i)AW4J$XIm7n@S1|k-a1~*cFA2<mz{?a~1g>HD72pdDe*i4l%()ZT z!SEsAD2D$5j%7F!S1u<oZdQB*5+LZn6nX-uFgzHz7sK?4W*Wnf0cSA$3~(01F9K&X z{3`JHV8g5%r&}P%VG0L;r!!1%LF6*r1kC3%oC3U<;UT~)8O{bSU^oZ3Fxd1DBXUI% zQ&<FC%rJdhTf*>rz-0`70bI`T_rMhlp8&36_%g89WP)(47HSxd1-`&=ds<D{{e*h~ zI~X1a9L4ZMz_AQJ3G8GzAK2``1oZWM3d5U$dof%NoW`&hID_G<z?lp;!Ez^y;kLlp z40jL5^?y7Qq(dQx;ZeZT8GZ^lm*ITi`3$cGUd*rucqPL-feRS^F&yKskO_W=LJ`BY zz{LzV!Gf)X;kLkK4EF{uXLulR1;ZnOs~DbQLf~Zr`qI6I;WfY)7=9mE82&!61K7dv zcfe5$p8}3$nC{A)44bhKbYOz^z$py(2JXf1AmB8HM*?RsJOwzDVHa=~!>fR^{n!*m z5RCT=<cG*P4F8Pu=?tF-&Sf|POXT?scLH9_@BrYI3_k>1P=|5-p9Vo;ogh-o11@5C z1#mILZv&Swyc4*L;Zwln42NUQUcvB<z*QR4^*;jwuNDAj1J^J-6Zitd%Yg-+zhnN8 z*NzW>9Sna39L4Y{;8=#kannE#A@DJm!1}YHXa|843b^b`gY-JIK4_g!>w0<^>RaQr zlskA#V3^kRUGO#CFvMh!mhpGumPD@!X&FC|VOquy0!Deh>p!jIGkpR{fY$ND7^Zc6 z7Q?iTAIb2`z}XDbGX7zPX&FDB;S%79607S!t>bf;0<GhxGED3E=?v34J{MRwfq_Nq ze1>TmznEcK#;;`gX&Jwou)hA&I=+A@&^rDNhG`vN$S|$rH!)1h_#%dB8UG=}w2Uui zn3nNdfRlBbKZZ8XMH^D5(0!tWWmt&x$F5oXKOHUq{0}W*3yKs$cp}Wd>Q4sd-#Zrq z^Dn&j1M`g*n&JkPZ#(cJF#pta5is9|ps=x^2SC1q9tj<X!pY0l_RmHq*6n|QF1QtG zUDS32^vUWe;FO(Kzb&ajIVY^yTarzEPfecy7m-ce=wj6;KX5Gp&SD=)2jQU<-+SP0 zdPt=XBJnUV-;01gDdwM_9s}lI=w4|ig7K-ob<@^l)4y#*L(vG|qUnR|wgcwdSZoB= zUk1}>s+)oN76f~M`KQT00P_!X2jQWb-lT!N@?>)iK)%hyE-3KN$ve?!SLz@_cLd<) z0`SrRyd9WtanKVVnDS4bpKO5vz<2iWLV<sfe+8I--EG5{seBuX(ZKp};^c8)+Jb-@ zG!Zpi3d}$BJq=8MAzgrU`W!yUY=tv*Pz_Un`M1*l0p{PgzaF5!D*&GY=ASON!GPjl zmp@1YOdV97xf%k#6G+QA%+_pMih;oTpbU=!^Iam=0rNc&8nvmv@y9@5{&Dv{V79A- zDNaJ54=UhN0B#Upe^Z9e!2ARMbYT90|GmKaAbYa{@GfBgrViLErL9w|)9Ayd9l-o+ z@Dad#*NBb4{Oj+Nz`WQ04y+Hd7n`7w&wduiCeogTfbWy>HZZ;B(FbYNNy>Nz$f|(( z_w|>6`JM^)Ci)9>;&cKq-)CbXOUGU*Xb%YZ)+Drx1K;K0EHK}G<^nT7`Dv#S{YO4Z zN&(gfrPE&nZ^5XbmY`8F6<F=LU`oNy5b#|d?n{C}wnM-x!1^GAn}GQy4ts$4&LFqn z03~?PhKM|P4`Fps{wct|eKp{Jnu%TvQ1AxeTiaocgmFt3TrxNiSRZ8Y5n#SC$Lqj+ zgP&i3`F1zI0n;loeW?A(U|4%?&euuefcb_!qk-$~djZUMnV27t{|#WigTOvuzGKiK zVA_eq%tg86fd(Dw6Sx7GZ?5t{K>Aa_=dp^@ZCVJ-cMy6PnC}h~*Ab(Jb?UyrIH>-< z8jpm4Z>r=1=DWol^JS>hKL^aWv&rfNYi!d7H!$BNVjr+RsL2jX>?>c6`->1b!O)3B zYLbY~_4m=31I+hSS`DlZGFTLV&jjG+UFvUpF&|jf-{*lZA;9)Abr+AOUF&aQGys@y zKeGi`A4K3l06r0b{|08;REd$v{_>><UWS10>U0E{?-rAAquyf3H;$hN*6RV<DI^d0 zW7a1;!1M|bktBkbf%QSynS%QB9e)O=AOl<YECS}6uWSZRJYjdx=%k56bkpX11vCx7 z-GKS_HMa(&-_=dEKbav$)1kn3e_02-j5W#O00VJ1)hEymnD4pN514Pz;|b6YGP`3o zV>^S40_IzlYy;*SZy7z%(XjW4-mIkuFyESFF0ekD;G_VUZ{FesuD7k#f3Pn_AGXcZ zDqz0p(golxSex2mjJ6)5{VViA1+)#ow*vF+YVHh3e+-yy_%l>o2vF$O6Rg=pkqOK; zR}n?P`XB;fz3T6k_L%4KvOO`L5r<D0*6)|^ww}Kn>-prwjlrHq)j6To+4wlmrkWiQ Kp4_WfgZ~?ms%D1( diff --git a/NPAnalysis/10He_Riken/include/ObjectManager.hh b/NPAnalysis/10He_Riken/include/ObjectManager.hh index d8120bdb7..3e2c4435f 100644 --- a/NPAnalysis/10He_Riken/include/ObjectManager.hh +++ b/NPAnalysis/10He_Riken/include/ObjectManager.hh @@ -106,22 +106,22 @@ namespace ENERGYLOSS // Declare your Energy loss here : EnergyLoss He3TargetWind = EnergyLoss ( "3He_Mylar.txt" , - 100000 , + 1000 , 1 , 3 ); EnergyLoss He3TargetGaz = EnergyLoss ( "3He_D2gaz_1b_26K.txt" , - 100000 , + 1000 , 1 , 3 ); EnergyLoss He3StripAl = EnergyLoss ( "3He_Al.txt" , - 1000 , + 10 , 1 , 3 ); EnergyLoss He3StripSi = EnergyLoss ( "3He_Si.txt" , - 1000 , + 10 , 1 , 3 ); @@ -129,6 +129,18 @@ namespace ENERGYLOSS using namespace ENERGYLOSS ; // ---------------------------------------------------------------------------------------------- -///////////////////////////////////////////////////////////////////////////////////////////////// +///////////////////////////////////////////////////////////////////////////////////////////////// +// -----------------------------------Random Engine---------------------------------------------- +#include "TRandom3.h" +namespace RANDOMENGINE + { + + TRandom3 RandomEngine = TRandom3(); + + } + +using namespace RANDOMENGINE ; +// ---------------------------------------------------------------------------------------------- +///////////////////////////////////////////////////////////////////////////////////////////////// diff --git a/NPAnalysis/10He_Riken/src/Analysis.cc b/NPAnalysis/10He_Riken/src/Analysis.cc index f1003bc7c..d7562d0db 100644 --- a/NPAnalysis/10He_Riken/src/Analysis.cc +++ b/NPAnalysis/10He_Riken/src/Analysis.cc @@ -38,7 +38,6 @@ int main(int argc,char** argv) RootOutput::getInstance()->GetTree()->Branch("A",&TT,"TT/D") ; RootOutput::getInstance()->GetTree()->Branch("X",&X,"X/D") ; RootOutput::getInstance()->GetTree()->Branch("Y",&Y,"Y/D") ; - RootOutput::getInstance()->GetTree()->Branch("Y",&Y,"Y/D") ; RootOutput::getInstance()->GetTree()->Branch("ThinSi_E",&ThinSi,"ThinSi/D") ; RootOutput::getInstance()->GetTree()->Branch("ThetaCM",&ThetaCM,"ThetaCM/D") ; RootOutput::getInstance()->GetTree()->Branch("ResolThetaCM",&ResolThetaCM,"ResolThetaCM/D") ; @@ -55,8 +54,7 @@ int main(int argc,char** argv) Chain->SetBranchAddress("InitialConditions" ,&Init ); - double TargetX=0 ; double TargetY=0; double BeamTheta = 0 ; double BeamPhi = 0 ; -double TrueE=0 ; double TrueTheta=0 ; + double XTarget=0 ; double YTarget=0; double BeamTheta = 0 ; double BeamPhi = 0 ; double E=-1000; // Get Must2 Pointer: MUST2Array* M2 = (MUST2Array*) myDetector -> m_Detector["MUST2"] ; @@ -69,15 +67,25 @@ double TrueE=0 ; double TrueTheta=0 ; myDetector -> ClearEventPhysics() ; myDetector -> BuildPhysicalEvent() ; - double E = M2 -> GetEnergyDeposit(); - TVector3 HitDirection = M2 -> GetPositionOfInteraction() - TVector3(Init->GetICPositionX(0),Init->GetICPositionY(0),0); + E = M2 -> GetEnergyDeposit(); + + XTarget = RandomEngine.Gaus(Init->GetICPositionX(0),1); + YTarget = RandomEngine.Gaus(Init->GetICPositionY(0),1); + + TVector3 HitDirection = M2 -> GetPositionOfInteraction() - TVector3(XTarget,YTarget,0); - BeamTheta = Init->GetICIncidentAngleTheta(0)*deg ; BeamPhi = Init->GetICIncidentAnglePhi(0)*deg ; + BeamTheta = RandomEngine.Gaus( Init->GetICIncidentAngleTheta(0)*deg , 2*deg ) ; + BeamPhi = RandomEngine.Gaus( Init->GetICIncidentAnglePhi(0)*deg , 2*deg ) ; + +// BeamTheta = Init->GetICIncidentAngleTheta(0)*deg ; BeamPhi = Init->GetICIncidentAnglePhi(0)*deg ; TVector3 BeamDirection = TVector3(cos(BeamPhi)*sin(BeamTheta) , sin(BeamPhi)*sin(BeamTheta) , cos(BeamTheta)) ; + // Angle between beam and particle double Theta = ThetaCalculation ( HitDirection , BeamDirection ) ; + // Angle between particule and z axis (target Normal) double ThetaN = ThetaCalculation ( HitDirection , TVector3(0,0,1) ) ; - double ThetaMM2Surface = ThetaCalculation ( HitDirection , M2 -> GetPositionOfInteraction() ); + // ANgle between particule and Must2 Si surface + double ThetaMM2Surface = ThetaCalculation ( HitDirection , M2 -> GetTelescopeNormal() ); if(E>-1000 && ThinSi>0 ) { @@ -89,7 +97,7 @@ double TrueE=0 ; double TrueTheta=0 ; 20*micrometer , // Target Thickness at 0 degree ThetaMM2Surface ); - //E = E + ThinSi ; +// E = E + ThinSi ; E= He3StripAl.EvaluateInitialEnergy( E , // Energy of the detected particle 0.4*micrometer , // Target Thickness at 0 degree diff --git a/NPLib/MUST2/Must2Array.cxx b/NPLib/MUST2/Must2Array.cxx index 8041b3ae2..ff4a1f8a3 100644 --- a/NPLib/MUST2/Must2Array.cxx +++ b/NPLib/MUST2/Must2Array.cxx @@ -563,6 +563,36 @@ TVector3 MUST2Array::GetPositionOfInteraction() } +TVector3 MUST2Array::GetTelescopeNormal() + { + TVector3 Normal = TVector3(-1000,-1000,-1000); + + if(EventPhysics->TelescopeNumber.size()==1) + { + TVector3 U = TVector3 ( GetStripPositionX( EventPhysics->TelescopeNumber[0] , 128 , 1 ) , + GetStripPositionY( EventPhysics->TelescopeNumber[0] , 128 , 1 ) , + GetStripPositionZ( EventPhysics->TelescopeNumber[0] , 128 , 1 ) ) + + - TVector3 ( GetStripPositionX( EventPhysics->TelescopeNumber[0] , 1 , 1 ) , + GetStripPositionY( EventPhysics->TelescopeNumber[0] , 1 , 1 ) , + GetStripPositionZ( EventPhysics->TelescopeNumber[0] , 1 , 1 ) ); + + TVector3 V = TVector3 ( GetStripPositionX( EventPhysics->TelescopeNumber[0] , 128 , 128 ) , + GetStripPositionY( EventPhysics->TelescopeNumber[0] , 128 , 128 ) , + GetStripPositionZ( EventPhysics->TelescopeNumber[0] , 128 , 128 ) ) + + - TVector3 ( GetStripPositionX( EventPhysics->TelescopeNumber[0] , 128 , 1 ) , + GetStripPositionY( EventPhysics->TelescopeNumber[0] , 128 , 1 ) , + GetStripPositionZ( EventPhysics->TelescopeNumber[0] , 128 , 1 ) ); + + Normal = U.Cross(V); + } + + return(Normal.Unit()) ; + + } + + void MUST2Array::Print() { cout << "Number of telescope : " << NumberOfTelescope << endl ; diff --git a/NPLib/MUST2/Must2Array.h b/NPLib/MUST2/Must2Array.h index 538fb699c..4f6349f26 100644 --- a/NPLib/MUST2/Must2Array.h +++ b/NPLib/MUST2/Must2Array.h @@ -74,7 +74,7 @@ class MUST2Array : public NPA::VDetector double GetEnergyDeposit() ; TVector3 GetPositionOfInteraction(); - + TVector3 GetTelescopeNormal() ; void Print() ; private: // Root Input and Output tree classes -- GitLab