From f6feff2f077866fa6d227734dee68787adf5bf30 Mon Sep 17 00:00:00 2001 From: Unknown <unknown> Date: Wed, 2 Oct 2013 10:00:25 +0000 Subject: [PATCH] --- Inputs/Reaction/10He.reaction | 6 +++--- NPAnalysis/10He_Riken/Analysis | Bin 57698 -> 58961 bytes .../10He_Riken/include/ObjectManager.hh | 9 +++++---- NPAnalysis/10He_Riken/src/Analysis.cc | 14 ++++++++++++-- NPLib/Tools/NPEnergyLoss.cxx | 14 +++++++++++++- NPLib/Tools/NPReaction.cxx | 4 ++-- 6 files changed, 35 insertions(+), 12 deletions(-) diff --git a/Inputs/Reaction/10He.reaction b/Inputs/Reaction/10He.reaction index 88131ef73..9b4a9fc28 100644 --- a/Inputs/Reaction/10He.reaction +++ b/Inputs/Reaction/10He.reaction @@ -7,13 +7,13 @@ TransfertToResonance Target= 2H Light= 3He Heavy= 10He - ExcitationEnergy= 5.0 + ExcitationEnergy= 1.0 BeamEnergy= 550 BeamEnergySpread= 0 BeamFWHMX= 6.232 BeamFWHMY= 9.069 - BeamEmmitanceTheta= 0.69 - BeamEmmitancePhi= 0 + BeamEmmitanceTheta= 0.01208 + BeamEmmitancePhi= 0.01681 ResonanceDecayZ= 2 ResonanceDecayA= 8 CrossSectionPath= 11Li(d,3He)10He.txt diff --git a/NPAnalysis/10He_Riken/Analysis b/NPAnalysis/10He_Riken/Analysis index e4e550e76ef53467904d0205641e15791e9579a0..f5b004c15c19c002bffe3614626ee9f880d82d44 100755 GIT binary patch delta 21990 zcmb81e_T{m{{PPm_u_yc;4e@R(9o<^L^H!OBv*>kloE^eCnnZdW>jRD=*4gh7RS)5 zt(6*WZA+i+2Ze2vR#TZ_xqZ;gutcfc9MX(T&3@>7UavFvI!AWDe|(3>;eDU)bI$v` z&yRcVxp(fIuip#Xy(_FHA+p$~xkd^h>Q=8Aa9h^IDj^tjXhf_K4hX61qq`;|#D)n` zt_cxt#XBM0Nbg3RyU&{UMn!I1B*ezknrq`})9F~$=ey_=EjPV$={5Uad@^ofe%^<_ z{j#NW4l<21*Env8?=NDrGa}<+M+=wG5=6OcoY2bqJI!;BWMdpsjYr@o2O$^X27>jO z2zevIB!oOmce5p@K;~P=Qz36f_!Gi(OGhyiVK%~T2t^3<5bi*5Bj97m|98p%p?`Pd z@IHhR1XNT!h_D!8DZ(-Y`%{jC6$l<nQE5pEG?%DW*o)&w5!P6`b&!v-{J-E}1H$77 zn-KWfjPR5^65<(4J`4Fb>zJYr;eQbRf$$>2%Lx3u5)|T9$Snx3A-s;zi0~!?KU;%B zY}0V?mbocxf~H?w5M%novR%LbHgN9D+rF=87@{q3+*I$fl0tMx*cf>sU)AM8uF@~$ z4V~hbnC<#7=V;B03{Q^xryO9x3^cye937tGjL^&*!qbO8OW&*yBM)`*<grdAEjOP> z<|fU2!j(MYzD|*6I=Qi<Q+!^e*|$q_+|`|GS8870CB?DcoZls7z}FGBH^G+psgpm> zF6RGqNgiI=NqtMFI#zZnqS^eV%VdYgyxx@(m)WTZ(_9Q~Z71z3=JT!;$KB>mSIUSl zJJqA6Q=|Dhm402g**7BD_(LZrJ)N2-rW0S@DgDAuW`}n}+kB<DdO&m+(c(Q!Dmz*= zS6A==@aJ(tj6Qgt4Y9wc5YwQjvl6uT6(ZBGxlGHzt|TFTPz=PQfZ|Wj8xT8UgeXYT zrCX_@kLVvQ#IH!-&$3t6M~KvZx{EVgoaigQ#t16?PIFC`@kX@Bg@+T{HCH(l^muim z5bqrd)_@)_h6DTBG#A^KyZ{b2ej7?3jskif(_EZ)7{9f@5UW4aTsaoMg#ydqFk6K2 zQ_(S%Z)?)w?wobn;)U3Y3^!W|{(%-qfXD1Wdi+5GCOmj-8~7#-6XO@57RbS*oS?Kw zD)PCk0tVptMYR3oP?L)ybXg%<Y!%`R7O(&b^5Fm~7W5DWtk|cyKEUKpkFP?>W0XGK zg0_9=`;bGorJ#mrp)4!^n`n`6*h5dq{Fh!N#30z?c&EK;@Ih1n@rG2iLg7S{53&+m z2?xs20-0bsScoQ_1Os;bGMon+-m;nyX<mZm8P*tx6+Z@2Q0*GDNE$kw1#;s15G`0> z+5aI`h&$0`d7O2b;B^eDLZ!*a^hTee2IvVHAC4{($Y6KDgD{Y(jDfFF;CpBRJH7$= z-$Z(Q%tWEdOWz4O9EYmcq6KyKKL_1#lvEGJG>d;m6ZZX5b8(c=!QB|;bI?NeX>@-w zY6u=~#eal&UfF@tULpz@bKK@`;#Q22G-Zrjh4PxVYA%+K3|+<R81=U(gX-ZZM4&}z zkzHUy2KIz%l!o}h%9t463w?St)aP&DtjR|HR4e~nv``+>AHk?*g%`8Sj%$WB%45W9 zo%;|Dy@@8ZTi`4P-FRgbAB5#3rOQk>v<ve_ISOO}@50_L!~?Og`7yx3nSn86&kKJX zD8zv_Bh+Va!9Y6t9D3R0`RL4JsFB^L6LB6FDFxnwJ}!Ylw7g7z9Zpjp`n<&AV~{Dx zZx6ZwDVUg0A$v^pLJKW0(1#fO(m;AD=7PN;2a-@gL9^z%0r9L+4f-$xeR!FbejEC* z@kpq~n{fV)hdsN78&LNa<t%v~W6FgJ*j;)J@*886iRCz|nu|eb531$p({v2NCab_f zn7e<*7(rJ`ho1zOqK2IEIY=Ld!F_PRK25*CARP=wPs;dLF({8IW9;fEoPR4}z;5#G zXiIx8kHuWR12ryCYEZ+eRvC0*XrUryRPV=WT7XfHF@`}KBQ}A3uy40WA-Zq?DrBEU zf9!|xU!=6a^B6Q~?_z-rN7&>yplYX);eV}(<V6_%8D{~TkTsr*8f`(B4Y7C}3hY21 z!;FmIfrCt?i|&P9gTwZi`m2#D#7Jc-g-sEw^wB-IIIZSHVHrGzL9|D?x_t}><8fBl zUGfg%-O3<-5hcBV7O<!GY*g$vbU~SAulEUfI1L`#2G-)>_d^ZIaXqxTv!_j(aM`51 ziNiB;7S6lxj>Qx2x&QtV!*dqRS@6J|l6luJoL@43&Vpc6PEk>jK>E(f|2RBn%Kh`_ z&;M_6<6){Hqu1_~-Y8&F&dt~TY2xH5Q%2{`yDWdseRs?&F<%*etGVsa81vD~Pni{? zdSE5u#Gf!9ALVrjbKn&N&D<-t7_(6Ww30weirD${(STjVvm8<QIPiDbf5eAh5I;_d z&l|KDV?!4&ozL$S_gL{6uy=pTd|Y~5;-83TnJ5n{2oa_jbc(s?EJpJ4>C1P4H1VF5 zq2!Xa7&}h!FEBqDiMlHseZuRpJ_A!GxC;i`cuLeYmof|s6iN*1JW5#1r^H32Akn}9 zT4n|nFkMIqAB!mAxtkJoT0n^zu9y-omr!Exlv2WhGD>(|PKon+1tm_GN=gi(Rg@S- zHI$ea)=@r(%do)%oQE4JG5l*OF%>>b>Bh|oCC=LyC^13QQ({KgLiq${FiM=D+bA(% z_$a3dv7Hj5xrs6rb2(*yA(|;MNcU1=>NHw-fCa+=N;GR5Wxf!HDKQoKDKEy%Lpd1d zC*>tVoTR)|h%=Obz^JF>iHxr3G8COgnU3D39415*B?dt(<#0?ylp`>MCGg-f<fFVF z1CeqhCKAfaF$GeN5@HDD73g!yOq}<WqcK5HUWw~9<y9C1lw&XkD6hsCpv2{CBIQ_& z0m^KQ0m^F(A*NRSdrHMI%@5zj^yTjt&c^+w`FkJ`@a!qk{pX<d$6?I)lSW{`OyBu_ z;POn=H85}t{xeTtU^b4>S3<P6AzNVJcrxSKn<eLxIi}lvk|&bU$)dg9kcvDiQmLqw zoKGH1u8~|oM&pV0a>+BuXamt+EV+=JO)ivNM9wAWNp_R-$k~z?kPFC0hEx<&QAkCq z<PtKg)E+Cjlw3>}lFP`Y<c@Pt3@j&?liMV(AXk!`C0CMb$Ue!d$Q!^$d%aZDP*KYS zwUXD7>&P{dH<0Ve<&rm&8_C6zYso%xq2y=DP2@bub>wDpw&WMcEnGpjXGldo6>UtA zDtQapPmYz`NbVpD$=k?h$Q@_d0zNY4OSFLG?PT6IwKq#{BFB<_l6R3492kG?^-|GH zMJf~2O5RHzOs<jKLe3zUOFlr(Bo|9=BWIHfB_Af|lJg|{$$8{#$;ZeA1{E1n(LqHa zIaTsWvYQ+$`3$+3EF}L<E+u#HY#GQi-(5~_lgt^wT}f`1%$dMlL-tALjNmplP*E=x zoEqG<<XXv`BHVT48p)g~-1X#g$(%9VjpSm<oH^V+a-n3-Anqn|o@CA>?q+hfWX>pt zyM>Ajso>P&ZX>5k<`m=hlVc@ws&RLag=9`S?la_$-`E12dfdFmwzo;<6y$c1n<aB9 za>tT=k})OSgh4<>y;N{+a;K7OC3B8)4<^@uG1%m1`1iA_zN$Ffk(Yn7o3}x3-um2h zi<?U?|0CwYv5CEoA)|QUn(AI1IA|=y<Hn&g`k0%?#`mEj&~yC^sDiPhzqs$)9s02( z^OLcCdKW?M@s)JvnI&)#@u8Mh)W<wKwy(j4uKKv*Dy<~_!|`dF8w*857Q&o=09u1O z-lJ7%Y5s9IL$=aF`{%eudObVA)PIC=Q~gJ+!lzXW?|P}tPjgivr`NOHO0yC&l%r14 zNBD2hz@KdcE*hY^>bwCT32Cl{JepkfY1NL9ncEoGQS}Amcc3M`9+?-Eix{uxV=M2a zw*4sDms;Dt;n`*-`WA_fCI$<z%{Z+B+CReLAY|_-l6r#~WhHl6CK`STnXt=#2WjbJ z`)p~#PFgA5-F40g)^3_oxSix$JLfs9P}mPYpl`gMW~;EDVS(D;rMAV=ekJY7G}1L& zAA!E~dTcuSD;!N=U5$|Pb13QQ(b!NS?=zv*D0aawTj|>GpkG>ACPsumlFIk1TG$m? zXojkpKCAPk8*P?*<{8-QRkg3`!>Xh-zvq{L5vV#-b)vl=OnE&r71CGd_-N<vNNg!V zqfrCf72DiE=w#!zVlf9-66i`fi?>`6SVL9K-0t<rVI*BSY`GF4T`_E5egjHZOqlX| zWGd;(cDiy3x22X6^fDURcEvXLBy`f1DagP7c<ZKmltam}PFZWo4VJ8d{BXRHrX~Hr zZuJ9MQW_eP-Fm6ic39eO&{}C^C2ajDHal`H!U<U&yZpT|IGd{@kw2CX;~}Fir=9h5 zhUXM5UkMHJ1a^Bp9p|j^cZiK42LzQtn;)}9tnpxb_d3&A<H6QFW2I}Kj+HeF2v+Ef zWhNipL2a-?TxZZ>juEftj4WT)$d0(3y@v$P;dJ@$HJTu6#AdZ?lozUzZ6-#kkgfae z=TL=g?~Ef*+BLG3pIQmqt8vc|s*%g`eFbcXYQ&Wk9bO2n*CW%&3fcO{*gS~ehavu- z6*6Eb>SkFZ4iLLWw}xtD+ZjtotuC^4e^6><dl#=%$kv6~724`jE0lowSNgUE@u3=} zS~c1Xt=EG@f#6(j>vzx^Y83uoYs3ouNmhs>-LBBgP=#zW_nbzB468=Aaypduppzw| zoa>m-8g#bqGAmvCQQYRy_28fzY`K0I(V+@uSlUmZ^?GC)*+tMBvL>f#4JD1jx%;1e z#2U?(HL6i+G&fWu+x}ChLS1C*s-d$hWPA4@(*-MJ1l>W|bFGB!iKvm+GuHCCO!_Qo zeDM68M^`+VR{R**Nj7EC2#no^Q;`+fi$z9tcUFMcZAcd;2-#l6LHCWoo{+7tP?kN! zY>WS5({L8Z4f*s7bo>4L+gdkOqCvc#>nvY(qjZZd@HN|tUqky=N6H#$X@>trtAO@_ zxO4Y<?6iwc25a&JPC?eBi8bNOCTrqd&ywxIWz~f<nw>qKMOkCZ)`eLm9nA{4Yz%|a z>zQnM7$<8W$*RFp+1j~E4ec7(wf^4jxOupptj^i)T`yDTr)lhoF7#hi<?Ql$Y$IjJ zVEbm+O1EW&&e`6qliu8FmCY9MdL&)Yx_MUoD4Mle*S2u8WwreyJSOpa?4*f51siv) zs}R@6oaRu>*C@pH%})IZyMjZYuT?y%^(%BA=IrIDYM+KjF~N0{?dkz5M|-|zewjPw zV!51b>j<8c6~U3`ul+H|_VUvAlbLn><$rn+ms+$U+pZg)j8y0B;^C*r+&Z(x>nTFN zS5MuAu9wZ52qU=5Z}o6I>_5JL=NIjV+BV1h#kay&(k-ZM&1pIKF7ldn*C+MWQ1hUX z;G9-wHQtO9=E3WG8|zS&)|qWykFCES<=X9G>t;h|pV9?Z{s~NHHNLI8!b;a(jgzi5 zr>%8HTF+KJ<CNyYQ54*Q(6DOZmpc9@en15yX0}kRXPz>r-Z0R3ggUv9I2zTeXEUQD zrjhFk8r5^5Kln><fnnv=To*LNrB;p6j`<j)t(!PAhjbIAj&s}51eVQWn8L7uQ8@F0 zb8Jp6`=b?gqe@=YbLy!22r4gPIkg5gXV?IU$>fQFbC`FL)!m`Nq!UcB12x2YkgNKj z;x#B_+eUkHZjdECt0a*M(@xS?or5?oSk~W;hdRR6tvPNsPPoqSV6{k}*JA+M_c%pJ zUOjcMJ)At>2a`N)jnN!LGmAY&=ULU3SDT!A7PuDfgj%AJDg!DnVr5H6&5O_&{!*!z z3NPKXDlQdeP}mihcJESBNX0ZeMK5_1X5E0GxNG2)Jx1-jnW1cYYw;%6zFOJ@e;*~G zrS$7c!tGWY{s=$qdo5en3LOoye#^f%b8o!Z@L@{Ec~zYw?+$-O)8W*q&WU1@W1s^! zwjEg7(u#%tBpKJyn$xlKeLR=D^%gvtn^CdU<%sbl!^rZbQM?U~@x%exmZG=L<$NPg zjI~G6zFOKFJbn&M(8zCZkjP~e@?I6G8O1W+U@3aR35@=>)_y){85<@Xw8gL8{8unP zzdK=mJ3C{T{UKyW!in&~|0YH%?;%?Cyx{4I6Qdw_6pOr3e#{e>(a1KAC6uzy8Ze{S zNiEM5vT!YJ0d3uL-d5o$vvyK%qtG(B-;W7Enk&aNN^7l_wUd^$*T0iJYuzrjAA@b& zEz#M3p;b^0zk2iL9HHhaxX(kFp+mpHgg2*+F%I~I@nKnc#!$>W!p6nj0;)p?qTs4@ zqzIC^VF=YIXWJ#g7t7$&D9RqZvV<Ba<Y>Fpa^@W<8u@)bYh^by&A^{=ekVMzCd;7j zYpneA_7uvg*gj+Wtl2X;Rw|AqusR$c@Pz*PUp#A$ot*6a>6kft@*-#Tv4+nlM?0KD z+RbnCzIHZ#YkqXo4rlYX4Qp?%(4C2g8{!LkJDuNrWe%G*$l3KP^VVq}I&W?@Q>XWe z-+BmcWA5?)l?T?~-f_^JIDLw9cdPmI^!3gg4w#uU<~Sd0HJ_YuQ}1WL<ZKaD)$H&& zn%`{;qrc4uOy|tYoD&Y1SIjJp{_zX8l8+JJeqnB#+2_iwt#Vko(dqt24<L5x$6U7r zCx>lJxcDIRe9Xhm^01k$R(#5gnRTV})`RARSxLtA%o)5C9*04@_2sp(6$;t<m@}Jc zWfxNf7e7G-ePD|A`RHh`=UMAK9}j1(^K~&&uvFY12<h7SMDL)^UgJEjfg`xKE?}oT ziaMcF<kkHR*rLv!Dfc5C`9&qodX|Ugv##B@1k2RkLlU~ZsFR7uV50pRIVIpWd6#17 z#tRK~vkZM_8A|VDXc7%gg&{B274Rw40(*lgqkg1uRNor^SS$q-!>j&C;xP};j_vo8 zl?ZF=;4BqN^y7Xru5hSvPVsTk{tFxAF|?!FKr6{w8om+{{%J^SKeE^iJ)WcF91eLo z9nt7<>jZ9XK>|KB$%MDnIftu7bpm^Fe46=Iq089%IY-J~))SFv*TAPl)MPm<+rvUC zidhtg#Zfq64TLdt>%R_!hv9`c=MXL3fi~eG(y>ZOr#aWOFlMOZOEYKAfIEIc3VV>s z2tO3|AeA=a>5X*-l!|c6o1^1{HP`vY0&V^aSt~jmD||7Y_OmGX{8&tDF+Z6z*qFA@ zGMyNTx)wTbjy%04D6UAwm5M9!Oua-&Ar%RBirySCaG?tn-gzf0whK8g1s5zb=KKA) zcClhiw&wGO<+pWrCY_a~|9oNlbTGyHpPMfirCz%Ron%*D>dT>!&C}hR^D`5^3`NEE z)akQkz8@@2Y9F_fF~!-v4ZY?iIx?F$J3slEa%Hq9T*ffvOM8Pu&$_JcWQH5|nC{#E z?sR-=#@&%PEPpqOZ>gRdYoD88pTfbua>&W<>jC@rwdeIc=JY!*8TNZ~r#wqK<=Kj2 z?L4#hm<@Nho$Ehs=;xm8aIR}Mm(9PeYfb259-5yVxp=4$6~~hNfq$7@?z|}~ACIjo z&IC$&Y|Zau7TtMqmt#Y4vP|4#uD!D%Vt7XGyqo9WHE*GL<6XD*n|sf~+vneL?feDv zM813e!YT8E-@uC-r%uTqVQ#r=ar9-m^Jd?;c)^_e(o0H9%pnU#nhWk5-%z`tPK&kQ z3u{dF9KES=_1K2IdkY-~wT{HFt1>WW&cjfw{3sB}!_2S^=RqZ8IL1Ud=7O=9U7m%U z4w;9vb0HT%x*;ne%NP$?3F(7e2YCRp3Gxi&UdW!X*9OU#I!{9SG=7DOi@{+i#zW%k zcX20V5@abPm(*(^xt?!?EQZ_#xeD?GWF2HCu4>yM{|tE;avh|=P4=shNsuj&nUH?S ze8`u+2?R<Y--N6Q$5&W8aPR^Y;|>P`O^|tzhasm!;(Y?K3^ECFC1fV#TF89JY5#^H z$V-m|0_z}0LDoY~f!qZ-3(^m{1d?xKZGudNtb-f_*#KE!;NWc>ltA9@4+Qw)&ML?) zkadvFkRL!EgKUO$;YRx;WE$k4_CR1f<du+xkcE(CkZ(isb<sA+ZIH=WYqUTPf;<U% zC1fnRG8@tuf`fb<<U-DZEP^b8JOEh=c?NP7WKTRas)fve+yXffawAsD#gIu@U9TcT z){$|$x}6N!2AOdJEr2H)^^mEM8Mr-Ng=e`dehL_YIvmtO;X^@NAjcp<6J!D80mu@_ zKj2P#17ro{HtHeUAisnB9WoXV021&hH4Snk<T%JY$UMjekVTM{kfo5%LVg12gG@es z1!{_eG&~<U1Gxb*2@bsgnF)CiG9U5;WHF@nGtxu$G{1VFPmc{)V<$#wTBU=J<WmvC z%;-gl`VY8v>1N!b9%l9;&Cz1cd9a6>x=7PYkfPb#4LykyzfV%lB9~r|JeSEl^~lq7 zKGk8QY6J~RWoaT<99|%RIt?LX2WpIWIJ8FHjDK*bKIF4N;0Bt}OK{yBiJ}t+4uNT| zDYd7%=&ozJxyE^%(fXe%w3Rq5>G&4J4Kbg7Fwrr_j4O?gI0n8Jr+>6`r{h6Q&jj_A zR@1P)7==`cj%;&iX`+l7DAN^|YW>C`T@EY?yOirH!tkI&tV6B^=JL`M$7b_s<XvUH zg}kegdXV`p8m9w!H-J!)P+lA|ZyLN@WyUQ|af~pBFHUp}!2kG&GHAPFK~yI5moL`z z`#{oTeGB4>%%_odiuu;!MCY_tv)iKHX2)X9aV;|Uh>&J6<)9)DrSw0DjC{%Kywka% zrYS<I)#gnPB|4UwcjLc-<9M_7Ax&S1>7<`*lD$ZCkNGZQA7HE=g-NJLri)#o!G)V1 z>S1P-X)(D-Hx9Y#P?Xl##eAo{mzlRj(|158n;~KaayEk!X*ax9lDT?GU&n4~&bLJC z=R@>GFmT>k9I`qmk+{kHYDtRY3v>6<`0xx&V0#Xk$z^>99zd3vNTMB3y8rqLeWmtD z7`wk1Nqd_!%M$gZwm=|Dc4keP=FDniv+3IreX03M8C$;p(jI18nadGnI?CeLF4c4m zIU<=OVkBm|6DUbH`<L~&gcI$u{|y92qvJKdQf^+TT+W@QIcI4i&XBm}@%jtkayoC; zFV!5SJkk$CSB#@Xha3NUtX<|BScpRSCaIR@3XSH;A+y<-DbviM%M$fUOw{GFI)%$L z{g_SqI;6-mmoH1~x&>3X$);&=u{-p=GVi;~5*?N1yUXGuqA*()V|MN#t;8Cd?v{iq z#^fSxA!su5H@VEKj6Pkf!3E^KE_1fgXV5z~A69%-n!X<PPaZKh8-1PM9x>m3u(#Qa zGCGi}r)*RWbD|$42U>V4=({6kxAMOF0Vsn`WtVIEUYjCHkYpEVDzob~pm1|3-fy`Q zmkil6W_`KVH4VhMNQ;g@wfi_q*#2e*EG<CHaG5s^^Y+8gCpsQ7hb~XlH}a|{o1t*I z=G=SOynXoyN2>Yw@<gL2GU2VQ(3qSQ><y{{0ms*2L-EqsV%_HuZ#eW?r`YenV13-F zf8!Jfoce41#3Nz)#xSuyjK@v#|G_Z5K@+FL^lDA4(0KfU#@Nr~|F1N?LKnZw<2rq) z{+2HOsWbYJ{C`x}z2RbII5n?^>mP;-pFDoD3p8fo!@XUIIO-xSd$_j&*c!>SyCd~O zT}69j*6UqGeODg0bfvyJivORE(qE1iZ$xE15-rN3d0ZDw{g?9pkI{Odo2crRg_j1J zyYcu;H^#mnqo0b!E8HX3$BNalJU$+)zt&wm9yjVxoM?&DUo*Okm%8iCJw#h~{nzf| zCwaWHhrXkyIN3vA*Hd_V^7yG<e|Es$(cb#<L{Sm1e;F^HjMp3D#S8LyZ#?zu6ZSYD zPW93E^c9cvO{(uJp6|=!kNYzE*S>UVTjC=8_a*7o{luQ6$N}w1;%f;j`Z4;YetK&^ zv9F)*Nfzh&p{F({n>7`aM~xb`Scoqj`nyhXoE`Y4Q-35(>~-R>JWPK~6HkQU@Pjb@ z^DxmY73eb2aJFKYBm7IeH=FvMIl}X{aW_s~PPR)edD4=6&6WDTmK<Tp$(EdJNxss_ zbbNJ_lCOYL@>NZ&1mv6TxRA>?;7^?k+_P$V@D1AHHJS^P(%quQ6MRQ^|3HCv--EAf zVtmREE*xSc-VDO5B{J}@CnR2o6nM8$Y<NR+@wgHaFGtGPJn?cQ-Ut=qTp(~yQ={ez zK*CRfb<jOGypA8Te6@0z*|AbL_?tQC<)_;=&6SLhj*x{g31Jq(LWE@qs}VLMyoB&J z!Y2p^5so9AL+FNeKr%u)LKea#gjom+5tbpWM%aw-62jXEpS*=v?GNJMIKnxEZdmds zBcvl_AxuJ;g|HA|8NzCW%?K|cyp8Y)!a;=N2<H&G;Z2!jgmi=~gh>dq5EddVLs*Tl z8Q~>_w-G*hCk+O0a2(+rLN^}_AfzK?AxuJ;<!cyJ86M`;G&8+AIfA`ZgQjdVr&gy# zu=9%1`1{SJAhzrcc$u=y{9AQOINQA`)BLD<KsdXma<qA(dO!qwYXe@e%rN_{LYl2e zbCo%6)qsdII6Q+lEsM>&SEa~bL>6ab&Kzn!wkkQCKbsh{(0mi=_%3<DLfn>^|Aqni zE5}7<q<27g^DmmKd9gXdJ0PM_$NNbhEN{%2-js0uW^!Vc`7qM(1*I8P=KtW3FC`UK znR{?3e>hoX{))p446=-BbHM5W5gczLt9i}3adk>~0bWTeSY_Ua2tII{v5E_)O{*_* z&RAu>w|bECFJAMz)r0hMT(#=WUXNZBo>4DE|50YvqqiGuGmfw@-bUfWR-+8#?GQd} zH_9;H)!@UHqYRT{SjdqbWL%c;VH*Y+=N3L}#US&%!P~RvnI|4@<?LUfgz;5&eC%JL zgb6H}@nM@%hT;8ue%Pv%Vb(j-Y}z36M8mol?HgpCaCi&WI*Y;tPdvQ%e4crt=407) zUJ*Q_?b9wy;4Ni**zA;H@{A6i>p|v12o9pjgN(P9@hMC)PuHa1W|!+gx%jf~e7O!x z5QD*Xu@2Eo;H6(ti}4I0oM?v$kbI$ybCd(Qd7<QMY?$^WyCKE9=E^n6BML0<S<G^1 zw^-G8U=Df|tQPYmIQF=?Z%wkn?+<M2%~11;0-FzkePFu7(tn5iD}GnPbs4=If`t4k zAs8PcvcM4?cp@Ov$B0}UL@GQ590d+~A0y^C(Cx8`Vj&XrRCozEL1BJLn56LI;8cbG z4o*{;_a=iC_JM~eyo)cPW+;j;kYJ?3N5Gj1p9E)v?HcQ7iCl$yf%6m|2rf|aUjd$> z@HGY$g^GgTtQ0BCpI^BZo&#Q>@Ir90!i&Kr3a<c{D!dk4rm(RIigGAWG4yeacp1E0 z;rGA~E4&X}q3}`g3b5U$=fIT;cg0n4mBPv3nh+Zz9g1}!g~$SL0NYJ`3%FL{d%$%H z^YKx=!nNQ=h2H@C6y5`FI?qOoI1EMec}0x)8Qh|9B&yb?a4OiZa3;7zVLtylqwrj? z_%T$WC197$2AX^g6tT7f{2Vwz;YM((!aKo(6+Q&cQ1}!$Q{hNV=-CRVfO9RzOVBYQ z9f~|9!PVdbh56-op~5r3ZiN?sixpl1E>*Z1T(0nw;L0Eya?m{wMU9f65xha+55ctx z?*rE<{4Ka%;nUzog>_uQd<rLkn<V4>XA2C6qFG5W9^9hvG;o{3cY^&2F9CNb{AchP zg`WTmWu|-;>>{Il>A(k2#3~66f)f<}8Jwyxe_cIT;r`$Zg|7f-Dtse2Tj3etTri&? zN(b(QB2P)M7+k3EpTTa0H-n26-UcpJcptc2;gjG>g}bBRniE)C+74vkV1tri61Z03 zLU5hJ#o&5{S(8SE9|ijqehS>A@XO%l6GrH~d<TjaB>`{e+Z6sc*st(!;0}c&vF<se za01ssKZWdH40b8ZU+Tviieds32@1~zrz(6uc(B5J^D;x>&EQOhUk7I^+zifD*bg@H z6vbI63KZ^+m0O|07lGXh^KPnG;T&+O!ZX3;3O@j@RJb}4*Z&$t@iY=_P`CkHtMCWl zI)y(6*DL%zxKZH%*r#xJEIpePzNjnCzh*@-8VOnyo&auBcm~+7@V(#;g&zT*QFs$r zoD7YDSHLcX-#4I$RTQ6s6BIrIPF0w9se=_x0B0zCF*sA<(co-_CxUYoHfBJPrzq|P z7bxrj7b^TD*sX9qxLDzx;8KOZ1eYs(8eAD-LquV<UlUS@OZgyO;cLOQ3iFq<bqX&6 z*DJgd+^FzlV4uSO05_dyT>n3SqWQccMjQdRC>)8~k2Zyef&B{Sf;$v01fNm33@q@a zH>L|Y)2;`*Y`z{BpI4!XwH4r9-~@$_fKwHY#QhJSiSWrSm*IoJgHM^yZAv!yZDw#K z*H^rSXnu`pGgt4suw4c2bYwQJ)(7Ao%3k^9;@~2Mxmr&HV>KcD#%f&*R+x+RAqr!$ zZioy;!PWXmg}GYKRG6#vF$#0Fo~^JCJYL~l;9P~dSf8ly5w7O*6a`o71z@}7xmqt& zn5%WS!d$HvEBU!vFHxB9S(GZw#d?{-w{kUKt|+)#U!gEp>y-*~wZ2MWuGVW5=4ySN z!gr#lQ}L^FwdXe|@y2Q-h_vpLoz7_SH?V#GDV#mTKRfYGFdreHi-p6PAP#};+g^vW zxA24c6RKc*PjQND*>^Z&gfkXh8v2>~bs`-K^;h9_V0DL(cCfmoNFnak)!)s}0}uYn zTv?m!TzS-dt~T8%`kDJ`(+8+Oo;PtrukP5ujZ%tz=Aft3o#~&Mw>+I}sK1227q|QN zvxM$=mhd20ZlQrn@u&KafYp6pc7oNP<|p9sfV%m>Qn0#@%wNFj#x?alaUQ6@{qKjz z357Tjc_MK~jN$eJL+JqarUI)!#ec6;{1;&Lr}wk)IKqC`z|9Ws0;}7&EbWc+V21TD zf;(7LA%VIF2%lf58`oS6R`&{c2&`_sa}aEQ_$=jACmtSr{vnev{#v7xVi8!~|KTWD z-T&cYJTOuA_7NMw>W&%pV0BB5*TCv_4Bh!$#rmW|8vUVAx6Zj19NN1lT08~j6S7rE z%r5z-Vt_R`3|2R;ISW>|^tdGviL4LPKXf4*;z=mfjbrwJ)$J5|@D#Q_G>`{YchL9* ztnL*Q0IOScT-EQw&5xb~hjvlH)&k!^p>Dx355JICcf42zwm)=WZ6|&PtnMN5K3Lt@ zrj_w{984R`|Gk=@cQt)d7!~|5!7Q-4>&ezm@t=U@m()nh^gn>rop8pbLaFSevka^b zau4IJ555(_mkFsh1`QtR!~w9n<;JxBq4Z3;5Zt2F&{X1S{~55lN5~<twaZFybE4Bw zsGEwM?PMTv02TsxBF<uQC@$we>%#(az<WOnZSV8~SluUvyCbRFJKY5ievDv^mqMX# zF!d8S)mjm7)N}I#`@;eUg4JD8vO2|Y?8L{x>b@`mJ09iBKI}8_!u_0X0jrx6)qvH_ zi)z94hYoE4tNYBH0;`)5ahC?ZLCZW$p91EG<x9m?P^cTr6oJ(pvK|1dJK_A*GGOQ5 z23GeJJJBgV>S8P=tq0lM45>T#%{J?<b^tCe%b~D8^spAJ?h&>VtnQw4@nAc3umYLj z1&2ag{cHrQdjstTt9y*aUBdjb{y_s5Nkz~o9k?2tjWLjlV~&A5u)0gyhn>=Y*@^pH z8nng#v^Nf{ZiLgoc&mJw;1m?<o?lsiz-h8Cw6)H?V0C+@BVhZ(0vtnZ3Kkd#RyVLJ z?G*nS*ibhQIt#_GZK0iTt{#eO12_wI*+ox)?aL|mf%_YnZ|JgEI`BQ%{*b$;V~vE< zl4F(m$AZ<Jb(W{2{_1v2A0dJLp@G&;{2N%^Axj%(rw%&M2dwTHHXf{Q{Zrm4{oBLr z_J<u7xCaT;9k0Ft2X`?G)*wIQ!UCTGtGlM11FIX;-8LLu@Tqowm$c`=>K0v}8BqL= zQ9BcM*+c_JTv)(0V0G8CDzLiG+jroD*lZ};%0K2Z+_JzP&pz6F99)h;oC)S2-U6P1 zWvM1PMi}oyVSl>g;7}(%4OX{P3Lj~w4z}1Ju)2}mmQL~CgVmk%ZoK@0{5`~TP#nM@ z>j6DI?lP+3)qjlZmGJ-C=E!^{ZK|0yAlmFmiEo(y?h5CS|DO$%n7bisdscMM|J+Dv a{P2d--OHl;Oq=xIH%rQCIC*fP^Zx*JGaOj} delta 20725 zcmZ{s3tUxI`uEp4?2Si6Ar%BsAwp#uQ#4I7^N%HtlU7IbLZ*e+(9EdFyl=c51(8IT zmX$VSX3UHgg=5TLD>C!)j?G}1p@~wN9@31QEGozI{(kH1?OC3A-_2*8XMdmdtovoJ zz1e#mJ{0lQo`|x9DDV4*r>_uV?y9w^lLw{ygrGB_5wSwJAe1i7^rXWX5kmAegow2L z-H=0IKLmb5;@0lDC2GeLLhQI`cy?TrZr7Z|f6m{t>v-?G?)~uoA-fKIv~B#Blc#)> z5GYHob&c(SG_o?Hu8X}@c!ZH43O#oSqtGSf6;~H43s!f)@wgLiDBLiDecS~(9Bu^M zJ+{tk%aOLE7zH;PZY*35oO6uF$ppBE;FzPya8ux>!r`#g@5Aah^fLpev*GgLP*CwC z+*5G#;O4_Q$3mPG!Y#2C%WX-4>Jlq8F2;EYT&b-qgIve-B~I4EZGd|oj>ijdFRC*k z{$$H?$d~PNimh<l;9i6KGu&U{c)T7I;tj}3xVPY{;C8_6gyZo}P>5;+Cx4flBZkZK z5n0B6e(O4=s;@E2HFAW<wt#30ml<^}f3CVlP)Gl+sw;)KZ;NZ~JIEi*Ue`b4x>h?X zES70<WMr}|j7;gB#f;ep<6jV#pFUyfb(V)CldY*?p&G;DZxa@3eON9hhxwmx>Et0; z3(+AIDL6lNt-o8emyfqdcJ-EPTO?bHuk{pb!iw-xSmtJMl-UQ1oEBz(F)VXwVNEjY z+Vo^E&tTW}@^Me{bw|Qt>Lp)>cHp)4;%<4^lk93KfAAz*Plu(|BdlU4gysF~F#F7~ zhTb1m(btFZU$2eN5{JVKqOLU%AGgMoP;YqJp`zM~*5bif$m52m6*v|A<aHQ<8b8$@ zg9Th+`};ZyvHY~*;oxKVyiP*gcGcc$Evk^f?tg_MoE(EGGQm_enIaO!yXXr;5umdj zVL~EK5={>aex2wj7NZIGo-sW4+Ww~ze&aWWrx45pZ%h)R&52+EnBh2NAQc5-MW$L% zc#vb`>rj9#2}0!h4bK=m!g1FNas061$+r1^BzOQBW{ohyg_w3zJ~2F_Z2kjE+>8RS zC7HoVh;I`z=-5w96XMyU!2(;1ARGO0s#XIj2q@k)R1Ku!^pB|P7q3<sstZk&g#=c^ zj|s3X%8@~oLuKgb4MNl+Lm%Pd!uUT&hKsZ&^r9;FHL46*%HT>Dln7Nc$c|tvs(2J4 zL{iFsA9C>oBIGz=gx!(Ajh~17=Yj9}!tivn{V%6rVuiiK9WemPb_WZDyokFHQ4iI` zn3ZG24}o2%;uENWG&Cv`_#2Gcoi#kXzK}~0;Z8J(*XH}$pa#$foGPx1Lz5xHsCi|- zB_5p?_D+?zL-^%ig_?LV3e<)*!mAP!xEBE;Q6;KVfGcp4`(7x=ucH$-*#7KPAESnv zz~d2~EWy)JBWyFq_b58y>!=ZDl*FKhPk?*d{?p;ll>*aO5xj#07Jlm_*hY*)pUBtx z#OFw$;eQ?fR^mSN>MeFBVuD?mK~yO-z&^p@HXPxd0@WZV1zI0yKxetpZ&+8`cnpR` z3aT{4jv#<0*@OURQO7J{9NKIXG8nA#7%}o{8`@2Xx{rjtQv=VVla*<mc0MA$tc{9Z z$Z*0pR;W!TBEkvC0Mb`|f}L?8I*~K(_hQ(j+a+L$??(jl(dJ!YM{YvT%GBE2kBMoB z)<-U)kJKPO)V#9)8pAf)LK_y?it~_<Xr1mowAoN3=v3ibsPbyGAv(K?;Qu6}7i$@) zMgk9fY<TX1KWnTEO}G$EaI<ZngBmh06BgS3x1tHGGDPT<urJD9s|}NEw7Fkv(|T}` zHj%7C_g;?zIGw69GB6*V@PM7*9!%ADg3;MkhWn#N_M(6oJc_$?#HuJ2<K7vkb5P;} zWC#US2KQhB$V3EAFFuBoPY{7q<+-R!XR40HRJ{)+_G$&_feh!`O~xAe6Eaw-b?T=& zp#N{wI@PUg6SSdIBbSlm*;)opp_5i>HIR)yW7HU)amXkWI0Jb{o#FYHJ&(MOg4~G& zolY4*pAl%X9(DmwS{U~Okwe6zN<0}S%e8C4<IsB%fisZ|O2Il5RqpU6MA)P?Q7csW zarirZrU-qci+1sP7<E4d>@>;H2Vk&6>qLJ<%+07`XD&a10B>NNFSIlG3^KS4P3A<< z6DPw?R87CGol#XESDP&Nr;U_<Ieohvd&{Lbw}=piaH~Ir$ANP>t$(R(=pQAo^zSG; zX1;CxhS?d}7AbJCGiIO)TZpb~VjM2~{p-K{Bd_s4srh@m%|jnCMX#&j-J%t;%=k0l zo=lzrYmev;f9jAPmH^0RC_%RfV2+_Dk5Aux1a)bMJM91kzgvgKbc=iJ1fn~bo=C*N z>!EqS5S_p+(E2Hq=#FWWxDKXMBHkX9h%ti_*U`QmESz8slt~4m9!QBiW>F&NLnu+E zVU%d?5tK;KONquAO^I=kLy25ZphUOLrNoppl@g7bN4Z6aS(Iqqd`gUm0*fc;hI1*W zV#=VzI4GpV(qTCzM?U4VLae652rZ+;cvw$44&#Ut17{N@268zide>IU6d@`oF{y8- zM6ao	U_W;0Y$zYD!e=Zc0oEHI&^kWl*MLg+Tc`Ar4XAh^dP5_d*<@#E7b+yb0Gi z%AU9!QT7tzG$k6}PuW|D^OTQa+)(zx6!arcZbrYS?2F3;Wj{1N+O)+jXno55=r)v@ zm?0={MIWIYAVho0+h9&P5Pg91cJu*CTt?C;|3`>1#oyew?7ZP;zj|QFhiMaVpC~U| zfk43bNrCA<3AO(i`iOrBa-H+Q!4Cud+9NB0t_k>Mpg`A7ILqT{V*`Q(x~7on*I27K zjm*B=Sfe<dj7AiVm5O_iQ^@6tGdh@-Coi3fGNtHC15}o1EL5CHMwvxpzT$!8EOM^m zEb=h2SMd<Cmz<?|7&(WWp?Cy2mz<*5%l@CoA~nV;#b_GjlZE0O@?3J$6)3t+AQzGw z6z7tQ$hC^6lFP_7iu1@D$(4#{k;}>Diu2k3^D3w)Q;Gs|CAm=XTyixzU-3e64LMhF zA^8B=t9UuNmYk)yh+Id`P`sMlKu%F?l~LiRB33EZlbgsw@kVkpx#>TwflXx0i>Lv` z<z(KDG}bEKN{%JhD6SwUkSi5$C#R6h$yV1&D$=PaQwBT88RSC6)nwl2HRdbcP0k|c zDy|_9BYPF^C40$PiVu);$Qg<ck#or@U@OqImWn(oVwJ%Waz0rot|QMSH(h28oFEsH z8x%K?i^#Q#Pm{~YHH!V@jpRzj=gH;d^2<m+(6xz*3M$H!!H?uha-rg8ay2<$@nv!i zIae`fio644uVT&^d9~y$#hf|v>c|<2IfLXikW<w7@5-qp&re0HGT@Yw*F+YIIkn_9 zlbe2H4RDId!(|9HpqNulo`+nkm{U$(EV%~U`@08=zg%{@$vb*f9`7RZc<V57Y#x^` ze%Jf*nXII^LkJ<B8e9@ri<9bsc>Z@&W}JL4t3x6cf%Y3Rp$humdi;qq`^-ZL@<LW( zd=}Kcnu4|*Ap$2Y_Zw+hakBH^juvaY_#ew|GYWd{AChKNV4+1%MCTuXQc=fyoMI!* ze;0cCF2>OKE3R>+zGGnOKSBGK+-WSUF-k<~2BjUH<|#(VQr{8V<{8LPh)PAr@oO=H zbxs5xMnHA()d*@yu(=2FT$kcci}!`1nMuE<;^Xu`i7G7hsjwg~d`f*^*kSuP@kcX$ zrFG)Be1~kK3otsD6imR0Cdy8r@jG1HL-G9^rlrAvs*oq_D9-&7io!{IKWv%D#-%C> z$Fy8!_rw)+Z03Hlmbhbbj}!8Ggk(GPI%jxFeGPVE0mMM<K&AED+C~+xvQb%c^w*&+ zOMMPS!0KEA%W8#`-~SrS>ABcYB40q~REm@Ee%r2b6;?jzE1BpK{!}VIEUsf=rlT5) zC)U`Fud>l(XYU`)h%c`AXz~8yPHFzlzXYs6adYv-#w<it>Qh#*FUhW9Jf~o6D?x)% z0w*g@bT>h#Dpv%vlEAENWb$@aX5n_Kc;entUoA7Bvht&ymBA`2mXnuu48^SM#JaB3 zr>s;~_A)Dz5Ra_{y~FCWyNVOtE6}N|%t82L7wR|HARWqTTUOe#!j|Qb`-fO*#sjJW zn}JNJ2P%@yIzwrBd8IZU+ESmgQH5~ygIMkG{|a5CBzDQ3ff&stQ3xL^7O74XUBGxY zFf$xejCme32opG5>JzAc75-;dhUySh26bOU#}ctWIJx^5?d*~}x?Q$i<6=*+p9Tx$ zv7?!QiVYQr>j!3dD72+Mk4j&a$nluXc7Q?K|8I#_sS;IdC0Y|IkrPdKtw4^homL<x zcUIu%P>CGnIp|bdzJWDA^Bycwf}Qu(h$mE{6k9t7+ESmgQ3Z1JFG0@&-GMIp-%Vs8 zqUg1%MC>3=8$BB;krU52X4GyXM|Y>4Q8tlE8R_)W3gqaLoC2N1-At%J>AbKrZ@b|i zDp7`AqHWNY`e0;_97q38#)eQ&pr`+LiCCbGszB`NPJuRs3gkqy@+TC?vP<MBpM=us zbSh<}^AHW~PUq-`+jfm@t|?Kbo%QqZ4izZN)*gkn)TeAz6G3mOl3ZbINNEdt@c+9+ zFR2o7%sM675-O1s|2vmLP2}iafzBzAle;yv3l_)<W(R36unijr;m)(vcdwo2<topL zh6Km&R%XRF%=RCG#wf|&Te?9hMyDAp$8YGjC89X{P{=MCc24(k!Z!cJ@^gp}3MH7P zT78@yeWXfWk)8ZsRPtWie;@pvwsErj6>F21yr*!VRO)j~cl`)8GIujFSDalN^cbT$ z5+kuYf@A%l9m%<lr9LM!tDxJTeMBumOx&*rFV>FlEGOAOgj_a!V0X2KDZ3c#2!9K< zmOty?L3Y+O-b?b#@P5|IxK~CkGBa-6MwetC;^~tJTtBg{)VC1rUNYtonqF1!F6eQO zc+6ivA+0^;_z7vz>KE&s9J^ksQP^O(;8ka#t&dfDWu>}cP(LNjKLzEfpV)vfwq<`L z;#7*Gy8$|<GdU^6(9Rwmj_$V~>=qc0>p^{XUA^i3K@PjStMvptf|sUq(e-A|4~8d5 zv{6^8`s{{!GZ*@kztnH$jAVy4JjV_mZ#>oTQp|HT%WDd5w9AG{l4#`m%^d3?U5V0h zw46&|+DwKPmIJ|ypFOc=m$R~~QI4xaHoJnVSD`9}FKa5O`IHp^pQn)k2P*F`O4>pb z_{6q&3njz~gX_d#kaD{&&b2Ph&Qd9@(M*(Lv15X1lugYPUT6wctD{@dB>#5LAFSh- zgE`%gKIN>H?}vHGm?KUPS@nI;<Wsxb4S+Wb=XATf?82%GOLnCy5;74ge>GK8p(=&1 zD!<aKgU0etRC=XY%52)jRf@-<P+bJYRq^I(7D{oiV-cqo^Y%g($#p}foF3<_irccn z?L~lELm6qq{a28_V`wY=-$LgYI*b0Jw6oXjj&6@_$MVhjUS@l{TP<-5gke^at?n|u zN7Z2<lw?QK<V(<ryBkg{s;kF>@lNH}RG-~+@IySP88;S>YDO$u<Z;FL9!5`CvM4$b z5F*CMZqK?D@g`TO!|;UP|Jrw`UnA|0JTE~NRP!l1j67C3?|EQNFQ)ktQ!(dWK%Z}@ z@9c$^z7fL3xcKi6e<z0LLu`h3g42h=A3|^#jzFpW5y%MFqV;AT2f>6i42S&SSuDav z`!VV*Msk6D38{F_2k*6#*hqD+@c&ncYZl|;&08?8xxpBEs0k<y!BPt0hb-G^D7}iG zdPX0EXc>JsJNn({RST+h_<%J3KQACPEPk2&J*X$7%>Rm@jqR2C>zcvKFj%&C!~+jL zG*Kqc?2tMYBStj`nrtGQMg1&0pI|rWn(t-1`{ua2e_vH{Uu&1U?wovK^jGfXjdJmr zeeQ*gRYS)vGu`i;srttQ@ou+uT7EsToBOq|Wvd7G_w3TZi`$PJ0Cr?gphErQzCso; z!u^AIg5Jq|Kk?<dhwgJPYLH_d+Te~qCC@)J$$dkEyf=4beEvxe&gkM=SB<OogN6v2 zFFYl8=HBf7`J_CTJGb?_^&DyHuH)@`nLjCUz&Iq0nlD6^`foc0-!cE-l|MKZcQb~# zFB#?^JRP7;YgyApKgbU!4RHUvK?WvuntTZ<mj*8;s}P|xQq|RE9u&?<tz|5Ev<O}X zvz3B*phe>hbdyrwR(om^$gDj);u_?l$qCjYOy;WpO~@|h=vi#=*N}e{LtP7}Aslti zB=i*Q$RBAogBdrf%D&QFv`vn~gjq}ti=q$w8wXqHo{Z?U7SS&!u8HU`c0}%r!HAB- zThGh4C_&>G#8Zl^Fro^TrZAXfbW4Op2oE|yaBM@;7`f{6F=F%8#H`1qDaN-WIQXMm z*#WTb4iB*Tg#2jAO%KP@!mgauunhG<hV%@kl8rPxc^3Xi*gAKZ1E9xkZ}#a>j$6Px zbdy{Bg}i@iO4iEb><LGLUHLR}9r%=p+;>I{PO~dTG^1p*pF?7HL)XcLlM}o5M9ihx zrx+>%)xuNM2(#H2MJdg$qtC&5`ODPQej8xm^f#qn4u#X-R1{U3g;H#(m)>a|l4d}^ ztjWc?Z2XVjv1*UYrPI=<bvtH9loTp+Tj)x&)wnntEWWw|E5%V(7q20>B&i|xk!GP3 zI~<F6yG36)A!DW|v?;drEY|JEtNKlE>$aX(DgAtH!JBCD5ULo{MCq6&r~4+JV<(I) zjxEhr7B|7dsS~B>1cl04+tTcx8OLxamhDY>;K7L>2G<Zu+s`(GMTGy)pI05sOLArA z@Ji6XPP>vgR|5Jlvb#PH_96Q+bUy>ce<nN4e9d_Mvp_&5Ka$uh<p^?6$4uHo?oIfd z>!>DP)ogDf5@&9n@|lcz<afOed>R&}!_jNQjMKuj`b<9jNS?d*(W=9bJmhj`eJcAr zI$4%K+Bqr?|A$$2zLWp+TDkwxk<l4=NVu#yP|$8yMx0EW-MvK-{$Jxie?*R&UF8`) z?a};MlM1HEHjj;qy?e}kqx*;{kIkMw<Kdw*W=)ezADbsz=HFbEm5=|+I&Y;^Wj>Yb zvPiC^h}$wSg-t_$_I?})c(HnW4OLnM`2%DjhUk-M>aCD#Aic1C4ssUcCdeYlO8P_Y zhOB`+0C@uPGGsHP2c5P(;%g7N6*3*N#=t+0Mqo6ZhGH%x{^=-Qg6ssj9WoPgALMAr z^N{(F0yEWW$Zn8)!($fYUdZ<#PeUGn<jXZ@AUi=uW0KB<>;yR)@~f`{fda@2kY$mU z5WnJNI}~$H2LcBmmqVV0Tmy;!qKda6J3;<0WG3W3$kC9`e1nJ}A3Pff@YT4Pkd=@t zArC>WgY-kb3CVXp4nw9uo`4((c@{Fq!pTLP6hOY<4+J(s)<EunJONn?`3vNENTV?j z@Zg?49Wo7aJme6_*^s%A>me6HUW6=%Z2xT_up4p|WF6#q$R8nRL&l;h3m~l?I4Q!( zFvxX~Qz17(o`BpA*$i0&*&fe(jzDHW`XNU^)?y`<59!5XXf+wKf{ZofUNU3@<WwwD z4Lp)Q4VeNt6*qaS@o(4af8+kP0w+hHs6j%0$bm4p44DJTyRibuhafjXz5}_NdfbaQ zKz4z=3^@cc0sqg<f$R%82ht0<9C8-qM#v(_?T}j`-MAI4fgJVI!zd|Ea-e92JOtSZ z8ESycglzjW+8VMOWIkjMNd7l}AY?NnA87J$&G(KHu7#dR^@>zdv|$vv=$8W59}}WS z_MYFaCH}1`ZUC92$Q)4gG1+r|N7I0^6_qk~zG41|IZKhb9VVyclKDxl%X0JlB=fxT z`Ep@9Svx<*UA9|(Ilo=2<^zEM-<>dOT{1BLCU=LUvgd-1W*%(YNtcB<7Z|N_L3}aG zkP-5+1&L-ZQtTlI!?%2aVdjInlFak4nJ+(BkmQ;r>lY-MGaVlpyU;KTK<$}Na{`tE zvt;juNv=8aZv2+Z84HumPv9S=%25vM8c=6teIC|(<Oi@jB<t~8BY%PQbW8!kQfDkO z%&`uc3o%^`l_M7=xwgqi7bTe+9UoZ^yOp5SV63ovMt%Uh`LZ6rGvzOflHF}j%8rY> zxc11$7k7}}#fJGdCW2r#Hes5nf<DQ$MV2m3GGBp@ib~cjHe7#Z2s0W}$VT}kY*vEX z$ym$NY7Xow&=kh@2)SlyoS6rkU*+vql4)R4Nm2#Qw+uxprIx#4VZtxijCG(N<?9xk z@t~FLY9kL9c34+vn2|8)#x{!a!Xy?nk;&9X$ln(xns<SNC0&i_$SWrmCb^!Fi}3p( zBaoGdVJvJ^cQP--CJTC8`Q&$nNoF_Y!}d>EVzio38wf0A`=5xAy_O`J`wyd|tBmI^ zG2GFO=!Ws;cFa~g<dP*xuI>2kV4enV<(aHqVz@TpjP+r}q(HwOb;8%m40NgAmL$4z zWbdUN%n6vLcH=C`o$`&mXK6=sIVQ<M)e8%k8m;^w$r7$HP?3pd&Kdf$-`6ZP%-s%| z(U_&S$}g8DA<-|FBGF}@`_^NYyq96O8m4D(8_yj+cWZO_GGhgLK??M<ko5bUQfX7E z+5B;txx!c(L7Otz{M)gqG2|o563ukXu5FlymWPmPN01r`k$w=Kp1ilLquB{c)vmH> znPGMTwPpZwAf~KX`2qe}or5zE&srAXOu%AHuyU6h<^Yga*>8v6Y`J84QmY*33YgJG zRz#e69zGQa&w*5r-<^z0#`+A`8?bLDvoI}8;EXDJ`;y!q)UB_NYnWW(OL7gAn|)ne z@5$Z1F74|Q#mWeCQ-s(MVdDG%f8_;Vyao8opz2G*TxN>P>b$}<-!;WflloKY_nc{# zMv4`Y)Vv*O?vE5T>ioqP)V$TgeBUDuwixuB2Q|s_Y7h12)bEm(=9(z+e9J*+TZ#sC zUKB;&UHH8{Q5=pkPqh+_QG?!TB`RC-ysj0!OQQMva<utoYw`E!K`UE}!qz;mXifbW z{FOB;R0%$eF%QIwx)}3fj5w#xKZrHoyH1>sHCJCJeAn^(xwflZ2={(Fb883jc{}sX z_Tu|?=F0Y>uszQ!+Eah2{hwSATjI??yx0-nsXksDRp-s|^xoWoNuTO42ftq@m~SMC zO9?kbtxgoaM8a!{^gf(u7IqZBB$^c+#mgNrdOzqWe<<oK|Gu(Mm*igag!sZ`e&7}t zT=4o|xA~P@9C73HlH07hzNnWgvO8W#n{`HxEq>2xiD87HDH3g2W=k%UsNZSJf7<e# zE%}NS{rPSeCEwwq<hx*$d}j+6NA=1nt`F+f)Js<akF8!3e77xqt$_;`{xw(R3%)!z zUT$91VMHwcFOT_JB;cn#B;JD(c$G;Us5U%2uZH9cz_D<6Ckiiu32`M5c&zyy!xMl+ zDFqe+k1gGa7sg#BH(TmM7%0`zZI|J>1#T$ZSh(qMbKq9MJqx!DZYSJ6xTA1q;C_OO z#M7L3xNdN_zzu~P3pX8Z4%`a3XW_QN?S$J0cNFf-ZoHoQ6HX$r(2a-d26qeGP`I&h z)8Xd8t$=$LZX4WAxP5R(;m*MQ1Q&__vd6=9gS!Q8DBM`M>2P!4R=_<Aw+(J5+&;LY zaAy+n^AlVoUKNgq>jrlV+|c){#;%Txa9@s)53cFbl8sP=x+|ByHOVd6HtDGC@$z*L zt8yt`qAHM|tx1k#MK|}CzphD*WLs2c%G9#dmTaB_nP@h7cUf{IUo&nPAfJFoENU+n zZ&qc=&1K2z(<U$8qw4W{xv#8CBwtrvI8%NLJHErYaV8eZGHz{hr21Ux2{~YGYNY=X zp0PY7C$CLyS?e-AwM(#GlBKZW$4}`#`8I6$Vne1+9>*zPX2|kM_qx<b^#znq-iXs^ z^s4BUa@0Cx2&b_tnbSq<5D8voNH3CG;K8RLJ&KrvkJsJc?olK!uIuK0xL9_S-AwNq zylyA6<qeTFf5ltV{p4&p-C~`wS4Hr~0S9Z1GJ-b%I9P9#5nRXPV9ilRs7@Q~%0b4E z!NEEVGKLNg=Vk2(fw>n4tCBK8jVs)|@b&MYf&(l##Da{)Ar97QkYg0HR)fsZg9`<! zHpm=VnBP&qLFNcUfl!x0=19W}j#rr@b2J+DssgVe;$V$aMhLu(h=bKm8KH(SN^q6A zJVgSia>@w2`O70KL3Z2FbF!1J3nLu=TECjE3o}AG&Fo}cB2Mrnze2kn+<1R?IG8WO za0+rE^ms_V{({*`aUNvD`*P=oE`9o5ji#OYD(sM5d>1AiwD&d1ya@JwD<hxnV)1#g z6FXlL<AY&`Ptgp_QZe<*kgwu}ZZ3Y9yB?5ybR6`L5re=t;i+pde2f@|lb#yi2hIQo za~~rnxzOx=HN|WgWNN$sJW%5jaF)i;gNJDR8hDt-d^|Qn;~KD6<3qS-7o#=BaTw%i zd=@-G;~&AfV5h_;=Jh;{@rJR;*SITqt`@#Oc%jDXablsS;FI*_8uPP<B8?}3S8F^Q zT&D3n@Oq7xgEwlt4!lWYYcmw(P@rIF;~4QK_+^dvfVXIT6#R<D=fGRRPMcl<S7_V{ zOUmsUcL7(1*b+UV*b!2QLEvhzQ^jM!H5xw#KA<rlnAB=q4zAPq@8AZFKLPu%vK1pv zL(z0q5hH#EH)|ZlqJ0<2Pzu<iaV9ucV?Oao(0B?sMdJnFbcZce`C2G490hnAI8$SO zDw(D6LGUn*Pl3G}UjpZ79EAxzSL0-Gp3N59yeAa-n!)Yhxf+iE7iv5KT%_?VaGAyn zz#BC#0heq1BDf;Rmg;nWhN4n4s0LSSydPYn@lo&rjn9B<HU0@)r?H7kScAq1V83FF zf7U>6D4H~bA>d|>$Ag77)6E2XG+qFX)p!*+LE{&|DH^{GPA4OMm4T0-$j}V<@=B)0 zKZCP0j=@EEn8y69(W`NPaE``zgL5^W0L}yRnVHJKOepd-gL&XWjaPw-G=2$Ortxm@ zMvaex%QgNHT%mDWBwTp`-&;8u$iPXpW^fO<M&n%Y0gdy)wHmV|bsF=d*9MLM1omtE zCb;Q>6&lyL2NTVj!4a_dKIH!m*rV}p;8=~Lu<l9FI02lZad&XK#+l#@OH&MoB2(ij z;4F=w1P{}=1nkxLC2)?$?|^ePt_A05%)4?cUsL=Cin$uM#mcQv;~T(58s7{q)0kf$ zZPb`wu$6246u3g;k|<pND>cQ-FsRnJ3S6V{N8ke*e-5tI_&ac&#sP4H#(dGvukj77 zF#eh}#jP-C)_6GAs!gFjz~=@YjUNZcYP=Ghpz&sKipFn&(>4Ckf+9mxd<xFg_$)X} zV=k75X`BG|YTO;1qw%faT#ZM7^E9?5K#{L09tY3W*at4u_(gD$#+BeQjSqr1YWxMb zT;rd>6(P1nG*<hSA%(aRT&?j?aE-=Ozy~y*1FqG01-MS*XTc2`{}t@N%DDc21Vz(T zMT|HLZq_)8w;%X#FQ)@uD|&%F8V>`<YMcvB(0Cy@MdJ<Nbcg?l2;YVx!%=__fipEe z3(nFw3im&J<k1Nk;xfD&c-Rl}=*wL!{x>eTlItkmg*X2V=P+0A3HVkJcRDH<SL*|C z59O@<hT-Hkjk#JM2*zqcl^d&dk)<&g>q9igV%-wMGzC}dBQ)k}-K#NI>!UU1YCT8e z8t?>-4}o(v=3;%S#%H;j&(joKt>=TCmg8!DuEt!g7i!GadXW~MtM%0y4+fWM%*Fb8 zjmPojfQ_1htMyG9bG2TsF<0wbHRfu)LSwGhw`)8TEj<fOsJDEj=5IX%gS+heWVgGu zcopp2e+qXyu``Shf}Pt~mpeh60z0?8E_b}}gZb$Nx`1%G+lxzJ{^vKyF~W^oZ|(EY zKZu@C=$~U&fb~!D8^QXg`?<JR*FRXD1}><VJGXXme{xnH-P+ULkRT&p?U|~7wJLE# zpYToS!{K+q-uLC$S9`kWAC@a#?P4|I9+{csm*0=$c3=OJ{AsZM>H12rzMsNDJ{~w6 z`X)93j|cRReiwoD{TZGE>l;E;w#PWoKf2{-$b7kp1Bn^x4c5PT-4=Y@z&gc8yZG-b zSpUTJ4m^Hvo-OdJVt$*iZ;vnu+{1qMz|8|5jmJ39Kej&&1APaGR(J%df0X+mSl_~8 z4<Aw52QzdmjN^jOJCp_eZv^Y#RsYF?LjSP(3$VT?h8GV=^iR>Jg7r_`XMpw3oc9GE zsHg-}@cS89-*Mtsu)gJmbsrwI@Cn#z7_-V}^1+LJFvlCg`VJ6pf%R|YqwtQbbI`tL z7~c!l_hVQA*1l4<#Qv}V28OKuY5G#IzQ;xtSpW3>d$7K1$oS4`dxERG6FdjjKW$IM z2d<$n$?*wsA1IuIB^VmU<G}jI`A>uO9V6C)^-tGd4h#Qr7`IKvW%?6)Hv#S<@)#7( z!3-<~m*Lqlb>vsT`VJF6f%Ppa`lMXjt8WE62gA=xaVQwyiyB)rsDeWOl)cUMp#Zc^ z1@pZ~YH9xm&7aI|PV{f@C0O6S;4fhPTm82};Vl*6zo@{|Zv2W=xW<ry^8W>_?~{@X zqiB2xSD@@$i<4k|pOoYq7*p_I2780`eOX?kzq;UC!Ct>BEWj}as1JUu&)oqohWUHb z?ARfgz!tDyn@Ov}{67Zk+oD8uwKx4h`Y3U0F$4;ISCG+Q=U@VPV0{CJKZEs6K6ZuK ze-767XAs>m?2thNadDVm2j-3vRv{ED;UFlSg9+w<_04yl3G=T2I}fVc;z9KrV0{yh z&fT5Zf;BV{tnZ~W4UAn#@PBUDu?9X23vdLi@0pRAe(fF`kAkmm#)I&$gPnr~ct4Df zfm6`tDez+roCn`_9Cx)&2B!SZ2_DSx(_#F#Fg^*^ca*sI#%sf`1>1Y&r~)1b>${S) z{ym1drR{UF3<~F92{woE>tKB+nD#v!>tKdQg7rN(o(Jo@kQ@W+d-dFM6YgrA%O&?q z8Vm(r17)<#Krz@k$Zvr49Vfm3chR<}=-d+%5!M6ji)_;AVCP`?d13rqPh9`?t!B!@ z4E_ezH{WRh>w9Q)>g7Zq%<%Xyo(<M_j<LY{W<X87IR5N|IsBiDYquYI0jzJ}bOo&M zHgrdCOgxw;oULvi2M;@8c=$pEYv2R0b1?pTu-ZKc=k)Jt^|^L`q!nO&f22CFzCTjG zn~@-vmBDcwBZh#Tg9+q@@zY>^7m}4>{-a)lLf`SJUtcG>V3p1WJ6ldM$NvoT@7@nf zE_<e=e^I}x)XF>J<e-$NsxE!B+%4zbI8v@jX<gOgVDHvd-lOwdSAFyKEcgEf55@gr diff --git a/NPAnalysis/10He_Riken/include/ObjectManager.hh b/NPAnalysis/10He_Riken/include/ObjectManager.hh index e0a66004b..f42981842 100644 --- a/NPAnalysis/10He_Riken/include/ObjectManager.hh +++ b/NPAnalysis/10He_Riken/include/ObjectManager.hh @@ -104,10 +104,11 @@ namespace ENERGYLOSS { // Declare your Energy loss here : - /* EnergyLoss ProtonTarget = EnergyLoss ( "CD2.txt" , - 100 , - 1 ); - */ + EnergyLoss He3Target = EnergyLoss ( "3He_Mylar.txt" , + 100 , + 1 , + 3 ); + } using namespace ENERGYLOSS ; diff --git a/NPAnalysis/10He_Riken/src/Analysis.cc b/NPAnalysis/10He_Riken/src/Analysis.cc index a86e9e19e..e23f787f3 100644 --- a/NPAnalysis/10He_Riken/src/Analysis.cc +++ b/NPAnalysis/10He_Riken/src/Analysis.cc @@ -39,8 +39,9 @@ int main(int argc,char** argv) RootOutput::getInstance()->GetTree()->Branch("Y",&Y,"Y/D") ; // Open the ThinSi Branch - RootInput::getInstance() -> GetTree()->SetBranchStatus(ThinSi,true) ; - RootInput::getInstance() -> GetTree()->SetBranchAddress(ThinSiEnergy,true) ; + RootInput::getInstance() -> GetChain()->SetBranchStatus("ThinSiEnergy",true) ; + double ThinSi=-1 ; + RootInput::getInstance() -> GetChain()->SetBranchAddress("ThinSiEnergy",&ThinSi) ; // Get Must2 Pointer: MUST2Array* M2 = (MUST2Array*) myDetector -> m_Detector["MUST2"] ; @@ -58,7 +59,15 @@ int main(int argc,char** argv) double E = M2 -> GetEnergyDeposit(); TVector3 A = M2 -> GetPositionOfInteraction(); + + if(ThinSi > 0) E = E + ThinSi ; + double Theta = ThetaCalculation ( A , TVector3(0,0,1) ) ; + + E= He3Target.EvaluateInitialEnergy( E , // Energy of the detected particle + 15*micrometer , // Target Thickness at 0 degree + Theta ); + if(E>-1000) Ex = myReaction -> ReconstructRelativistic( E , Theta ) ; else Ex = -100 ; EE = E ; TT = Theta/deg ; @@ -68,6 +77,7 @@ int main(int argc,char** argv) else{X = -1000 ; Y = -1000;} RootOutput::getInstance()->GetTree()->Fill() ; + ThinSi = -1 ; } cout << "A total of " << i << " event has been annalysed " << endl ; diff --git a/NPLib/Tools/NPEnergyLoss.cxx b/NPLib/Tools/NPEnergyLoss.cxx index db84d18f9..d478a388f 100644 --- a/NPLib/Tools/NPEnergyLoss.cxx +++ b/NPLib/Tools/NPEnergyLoss.cxx @@ -31,8 +31,16 @@ EnergyLoss::~EnergyLoss() //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... EnergyLoss::EnergyLoss(string Path , int NumberOfSlice=100 , int LiseColumn=0 , int NumberOfMass=1) { + fNumberOfSlice = NumberOfSlice ; fNumberOfMass = NumberOfMass ; + + string globalPath = getenv("NPTOOL"); + Path = globalPath + "/Inputs/dEdX/" + Path; + + cout << "///////////////////////////////// " << endl ; + cout << "Initialising an EnergyLoss object " << endl ; + //If LiseColumn is set to 0 File type is expected to be from SRIM if (LiseColumn == 0) { @@ -91,7 +99,9 @@ EnergyLoss::EnergyLoss(string Path , int NumberOfSlice=100 , int LiseColumn=0 , { cout << "Reading Energy Loss File: " << Path << endl ; // Reading Data double energy=0, energyloss=0; - string dummy; + string dummy; + // skipping comment first line + getline(TableFile,dummy); while ( TableFile >> energy ) { @@ -110,6 +120,8 @@ EnergyLoss::EnergyLoss(string Path , int NumberOfSlice=100 , int LiseColumn=0 , } + + cout << "///////////////////////////////// " << endl ; } //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... diff --git a/NPLib/Tools/NPReaction.cxx b/NPLib/Tools/NPReaction.cxx index 872dae270..8517a5a33 100644 --- a/NPLib/Tools/NPReaction.cxx +++ b/NPLib/Tools/NPReaction.cxx @@ -265,8 +265,8 @@ void Reaction::ReadConfigurationFile(string Path) ////////////////////////////////////////////////////////////////////////////////////////// ifstream ReactionFile; - string GlobalPath = getenv("NPTOOL"); - Path = GlobalPath + "/Inputs/Reaction/" + Path; + string GlobalPath = getenv("NPTOOL"); + Path = GlobalPath + "/Inputs/Reaction/" + Path; ReactionFile.open(Path.c_str()); if (ReactionFile.is_open()) {cout << "Reading Reaction File " << Path << endl ;} -- GitLab