From edef94eb38cad96a33b9ff99a356ba2d2021ca77 Mon Sep 17 00:00:00 2001 From: matta <matta@npt> Date: Mon, 19 Oct 2009 07:47:43 +0000 Subject: [PATCH] * Fixing minor bug in NPReaction EnergyLabToThetaCM - Default value for second argument is now in .hh file instead of .cxx file so last argument could be ommited --- Inputs/EventGenerator/10He.reaction | 2 +- NPAnalysis/10He_Riken/Analysis | Bin 70134 -> 70174 bytes .../10He_Riken/include/ObjectManager.hh | 16 ++++----- NPAnalysis/10He_Riken/src/Analysis.cc | 31 +++++++++--------- NPLib/Tools/NPReaction.cxx | 2 +- NPLib/Tools/NPReaction.h | 2 +- NPSimulation/src/PhysicsList.cc | 2 +- 7 files changed, 28 insertions(+), 27 deletions(-) diff --git a/Inputs/EventGenerator/10He.reaction b/Inputs/EventGenerator/10He.reaction index 2c1407b0f..2c89aa85c 100644 --- a/Inputs/EventGenerator/10He.reaction +++ b/Inputs/EventGenerator/10He.reaction @@ -19,7 +19,7 @@ TransfertToResonance ResonanceDecayA= 8 CrossSectionPath= 11Li(d,3He)10He.txt ShootLight= 1 - ShootHeavy= 1 + ShootHeavy= 0 ShootDecayProduct= 0 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% diff --git a/NPAnalysis/10He_Riken/Analysis b/NPAnalysis/10He_Riken/Analysis index e45aba2337d941b4290b9fffcde8bb1d1a00e1f0..59f527f9d00c3ebb0de4ab83fe2eb345fcc3f2f7 100755 GIT binary patch delta 18921 zcmbW9eOy%4*2nk3IXWmWgCMW!fby1sKzb8R5xP)hs92~}*u|H+)MJH8MX7gCK~W&F zO-71378<1$1}T-+m~U9VprTP}Vp0xihGj;j&huM)&XKL%=bwkqhgti(_CD*Zz4ki$ z>^XDJx$=f#ZMC61Sh=twWcsMgcNwF9mEq6Cm;xbm392y@|Eb*=>kEB%t=<JP5AE|% zm$zc)FxRWEGIrHZHD2Y4^1RQ#o+)_m{XgdC82&o6dd!6@(Mx*xsCA`Ey{l!>HG{Di zHXgv-*)TI>ZVrRdTjRb182kFV!5F9UU@yit|6wrB)3~+=W7Z_qSPJ%Hfy}KRV>A9V z7>!!}9qx>+y=*Z0XguvU#=@C;!~UT_HjUM|t44z`;K@$GK<3>#2eL&TjJ0+cj1a=X z57FbPa2S3G9)u`l{bDfo)_4xKU*_7igYmwMt>htYW=q~BgE3TVRfXm!-BcsRmbyJ6 zoUvM12}O{J7s32&bUR3^|BN<FR*fSesGILV=D;hxn;ZNXi|VBs(Oa?otiFt8qZc%0 z)c&4+%)%JD`7ji{v4c5i`I;2>AHjKAT&TVb0h@_9>f^FEfHB)YyBx?F2#3)j442s9 zv*^$&qiQVB>?c6J_?y9qXo&i0ZrK0g?p+O15zs_<NIf7AHz2N7Defr|jOD^WPtCvu z^e9n^zywdmuK#K<Mrrj{@ST#w(}Ec*fISFd|K*{9jCmkLx<>~wFx@4Gy1}95o~kiS zGq@2BMDk1SW=k;wt;eAk`n?zVN&TaN+Jb?}lHC3T0ojUyLo-qTAQCwVZt3ZA9eaHV z^|}``{28l69CiB>`ZLxc#kHnCyMq}pF#1TxBTk8b8H|UuOlil!n2<U8KxCs}<vTsh zA@3vY_0DvnGilf@(oDGZ7ZPh4S3J#Li+iX>4lND(j~MhS9_MMc7^U2#q-v5}nh%## z+YQDs+E&}ZpJ1dg0K)dh5Oh&Wtr5u0YH*EK|0xdPJ0yEk(7qA-(+~LxZpc+SrqLer z&Lzi11gksTEC7?6FTl-eDduaD8mVACm5M0&rCfdwoQ`1V2S*#mtO1#j4tq3)p71;w z4(em_66I<&23a4pFHlH^LH%gmgjC6wqI&^`o<izmY6dBuxkwtl{TUZbTHqmC{m(d9 zOj3lNg+mLa2%NP<qNQ}SeTsDID|t8(^@kAhJ>8r=7>3*)0SEPiWh92z0|xW79j=4} z6;c*2=!5o%fFANS7{W)<A(#;n;b~qzQ!iozJ%1XJZ3W0OeLgylmi^(t1KI{-(1XQl z*AV>^jy2=x#BA#fu^S>F&9YO#d3>jr+2Vst&;##>+_hmKM`A0Q6I?McrPvCHf%jm) z!R<d0XN<qF`Vrh&E6ugDpnn!#>)jfIy~v3)2|Fr~!uo75GHtjt5EIeg2jCDgQtZHi zf$q=_1<JfeSl$Z<pceK2!SPguxc_KyW@bH7yGqO4U>3lxq25i)M9Q>_*f3qnRI_P1 zj!Ufp?cfv)6kwKXb9F{wF5(!AxJPL_%s~W-rS=zs8GB!fz;`&-@-f{sYV8$_%OIZQ z?Hyf+*wP-pgyqyd1}&!O1}p@_qMobQ;lU=pz}xKAf`mxX>;z)G6ZsZzv!#!e=kK6< z_ux30qP1FvxoDh}R(E2udm7%uqzG3H($FfUp$pQZ0w;5Q@>4OlRU`a5&&6<5BFpqy z%d-y~!uF!O<28ftOYt@#-t`#ZU^t|2@DSo1Xslc4vpC9<oTwVRWBe!~R2=#l;2ucp zOp<UcizLFGO%mfdgCwGoL$W*0h9t4t$R&yJ<&i`%^GV{sn@19!7m!5a6_Ui1Q%uqW ziB?JlWOO-6^t6H`I#5XxGeQ+fa*!lVaU?O`H6#(q%_I-uv`Z4>Tu%}QZUad~emlu` za0HRWc<m;MBczdJGp4iyRKQ%=L=ur|CW-S+3rVEMNs>ia?vTVG&_;3qV|J3Ij9nl( z5Jv<_B-Uk;7@|&+2>Df#IMFafxJM*o29g+JBT35#I25Qb2y-9FrHuKIyqz(Bl1TBs zB%?8eBnLAVN-_qgVUn>};gF2O5kc|}#)gs{!dN`Xp^S|r`440&$vZJ;lN^S8hqUyF z$Du}w;W!3KjzA8OJc1k`c^7hkBxav%lB19VByp<GAvqd3Kr#_IKr#tAK=N+n0LgzM z2S_F(2S}zM2T0z797t088++V~*^Csa$N`e~AqPm}L|Vz$|I^n}F^{nle%0Y{6dlb| z?Fq<E`w<l%rD&S@$o^f9f21MCj<9rmBF@YamVr_Z728{}nIkNdnCjY_1ZNQ=*;u<x za5gcLhPBrUo<ST;Y;{<m3CkfxBq_>8gPFvLJ!>xzoJ)*Qvi4lTdBlmt*@E+lQ;E|B z&m&GJP83`~oJkxjxR5xTIM^Z-#iXDvwVMQ&66X>#!R5sH#2we62&*72AZ`_0NnA|a zB)E#WoY*GVN?b`?E4Z52N^G$TMGYxxNKr0$GjT0(f#6!=2I5@7^~5&fY{3o0jl}7K zw-YxJCknO^w-Cn)-c8&JwzLNeMI$Nf)W9V80C5K~6Wm1HN!)RbBG61sU(ecG1-B5> z(!ae)@JV75u}yF*u|ILG;5Oo51yfYJRVeJFh@=MPf-exq5*G;WAdV-_6?~aEkvLm$ zCvhrqy5OtC>BNbG8L~bnlQ>qefjC>i6xALq6h>0yPy>@-H{x7kCfG!rPuxLAOPCLF z0dcEfnox3ziJJt|#FA4^Y!e(zTuEFjIF#6GA%#^aB1ut0TrN0@xR$s;a4c~HajxK@ z#5Uq=!89Y~G!mx^rkN?Hi8xVkB5@0Gtl(r~ODidYg(8&{c4Cv@G~y0oCODnAleptA zia-W2!)dI&RWQwHIY#0p!8Ehwn22qHvx)tQYX#2$TXKR)VHJuTY7j|WE_fz!EOCM0 zT;h1*T)}z7iNx80X-dpVB~BMiQ)5m#aiZV?;!NUL!G%fcI3xinf`y`(8sre01eX%$ z5;MW&#QDS>e~Nj4xPZ7-a3%jMHKYfv-g1~bA9!E5PZU;O?AfG}KCvh?%*17i=o~)f zK7W5=htDfH&^RlHe?9w$y=oLMyDzX$F0@6qyxw%gGZb9+8A5aU_WK4I8?ZD5Soy{K z`YRFKeauCd!2bN&m?Y&rKKlN6*Ks(<?n@3ejOGRR_fgXL%KP(`ZQL_0O?i@MrQMB9 z)}{F=8~L`hp#fjxL|s-CkL+HGElP@V(50n@Q2Q&~Fm{mnRn(HDnI5XKtY{>)LLEat zN&6xnJ2qMg=X1vfco$KNS2~a*OEW@^RJ)qLId-yn9~J~%7Sc2ev47|m`i+ZrnKY2k z8|&}AS(u;$O*gTXKR<4=`FE@_x=dtgCbBzq6F<Sk;sL_M$HGLb<i&n|@A%2)h?`8z z)J&BBuAA5~KH6n#h%j+fm?)M^oZ^1zlg*R;WiRGwCTv%96J_bqE<Xed6F&$OHp#?K z{2WXy!Q!c_Bc+;&j^A_>Qzk^a_&QDeEKJl&Ca&;p6DFItVLjGmqDnIne8WWa1JN$` zn8^#DPGO=+GNF3%;s+*~&tYxXWn;5uBmJ`O#x>Y@uD`J1?n$F+m2CKW^2`S(nf-3E zv0bxK@T+d);DgaF8~cf9gbN$xl8sm%kuh1hms>K-N+MsEF;F?mcV$dd?&skXBb0nT zZeqBS$={q9hkV7Se=aPP6wP#IrsF75zaEE#eIf6dI1(F0PYO|<=i?_u`e**FVJtN~ z%a>0Ysmw(~|DNd`NUzSNzqGgTi<8VsD)*Q?Kv~4&CWk4L_~glV8%tqB99-)s$0@(^ zW0Qv(>vff9<|t(dADelXSE{dSJXN)a#<(P9^C>lzzm{3vy&6-2{W~7>&}9E~H28M^ zsWl8ODA*@-8*3lxuYAvUJk;MsG4WFm-R`T`(6$BVovpi>c>jkdg};t9Y1bjiuwpGK zD%1?ptUw2&j~`$4a9|IA<kqP*`F!2OeIol(S<Q5MFs0z{+XOiq@Xm+1@7?$cQ8MK~ zNy-5+0i;xt(J_4eBd3(_xFu^0;`VvgLrMt`oibiI!t<s~g!`XQ8K*??fJY}{{GNF< z%0=nT*E||L@Hu>afuU1{d*~=eRN%r_`@e2K-MZHERa>FWv+yt8`e=aR49*|?^mJc- z?q8lg67j{@e(-yqIK!V`oZiDXheDeb>V|&<27<mf+beHWbN{KH*lF3n!hG{@Zb)C? z?SrYI*|*^4D!Qa$|Dz36e*M*Asjj4GhBn{1+R5QeJ~(@r%ceGddWNqH2Emf;%jaj` zfo!bJ?vG?SoE@T+^Gn$gzLk;{cqw+i`K;av^3vI{=l_fDZqk|mtN-qzs?UG(zlAH) zZtrp7`z}9U`i}hM3#Z+YDF*uwWMnJSquqwH5siatB`>S?(6lVGGi|%-X`Yem$%vD= z+w=$}j1QeY$mP#7JbSt?g7EzG5T%B{Ha)`GOI$Uu8_vi*?@rfC?T7j38Nq#S9V8lw zRiftBKGHUip~jDw%<$>)!r2@BRD$@oGwyrxrY@dOm;L?!_s<!hCiE94L3~p3H}Vya zc`NmN!((B}2L8okdq+7>f5}>^EN#_hIhtUL(lqTHXv>O38yxOk`qli<<Drk;dD`JP z@4xie(!E_Z1E5>lvGnWnCT%<H+e2!n@9)oVyV_-GKKep;KFDcL*SYY?|L$kmjtPz? zoSr7<JVw(>$0>~G#44Ot4b(v|AT{wr6XGKp+xHm+Cutodrm2ujL~FXC+BR!$*~VR~ zx9CQ$p(kBunGA<BwtkxSMG?C!O{+j#RwUYp*y;M)ptoN}SIdfY)d}oZ52UX93_3lK z{q#M(D>c-0tN8L9f6EzI)lJYcpF(gDI^Beiru$H8r|ULI?fkWNE2v$6=UC{v`C2=& zz4Ih&x0F5Bv1m7Y=#fWx_!E)+M!>sMQ(8+>+^9tq=;YYelk;0n^$oSZ#Q*1s;g%WW z5cTC>BZ!93$EZ9IacD5^#gHV28f$3@h*(h#S?Eso)Ld`Ox2>o8<v>fd-I&5?3Z<E` z9xWP*aRE-=oN8$EeQ*&C(=cz?<8jSt>8_lmEvM94=OpKUN^NkKOxW1Xj*C@fWDaE* zjhB`Zb=bqwmN%@h30pJC)=#*IuyhxhvN;VzePzeI9*WI^3vdF1o|ilsV5!hdZPHA& zeI<tI#L|6B5s*m~K%BG5)1KdmozP|YvXpJ4sX;9j;SqWBK1`IQ)Kbk_ID~=Lte7M# zyD6^8p(gZ%x?*wD3P$Q1Z4E={gP7ItgD;xF1Txr1jI?tP(bQ!+gvV>n%!Hz$;g&%n zGvv>B_+FOsDHShkMW@pigH3qdKvh4&Pwmilb}&HO301n^?D#EyH)wu~^WDoWv^QaL z37ORV7K-O67MkDi#}aD)0Ghu!PIhRomDp=p%04QN5{~cFqy44u+eS4*B)>%;r)fK( z%1H>K6!2U8w%9bs6Y2Eu&6l({VY3N)(;XLzPodBq7yev_;%|PF8QSX^F?5HixL9=I zu<rMF!tX|^nfHyuawFiv<|b_?RJ#sQ{J$M<)ch_Wze9!J!e#(0=za@@Q8P)2BK)}* ziW`1U*AHj1Lwg-09RHGvhqa=X-aiS)o2X_Fx}+zHuz48@3Jz7;;SBkWl@?{VTQ{ye z-N#}hKc6~5@f9YfkqIp^gkpkbg8US_FNflWpSqi5hITSXq}XXH`U^Kt>+!uV+-#v5 z588+GRAq769Hs4qYU6J8)9J5vAhu}nrJtbfTaT09!saIIP4`<URzacXh4AMqDE{X6 zBV>p6x?K#_IVw&SeK@B(?(0R<M=RA#rkyw)PZ2g3YdfLJ7ZJpp<NM4reJo<(<7g&7 zg^9~BuKOtzXQ2>&K2nw<c0b_1`AKGIC#o|z1{GK9X&$F2%z}l6mRn{Rs=)?P_GY-H z`z-A5(soGZwjqEEcH@FNd*so7P2m%k`tf^mgO%QO|IQ6}arx_b-SXM)iZYU~nKQh5 z76#CMke{4$NEydB<?T~e*3EggP*o<>b-w85q6Bc${0Jq4$Iow!w(p_U;$<4Dk`&?h zY-F8%OC!391&_U$3R+rj=F$JjP+sK=|FcqQ<pGvy$`n4&@<5*vyJ>~uwzNsHDNQ?D zyOFya`B}?QC64<S%yoNZ*Nr)08DCcrIQ;Nl#NQCwh|`LF4jd?(_yt`xAScRFw$tW2 z_R%I^Q1KU0Y{GaMoA@sUcPX3r;067rM$>lq-h&o*poKnALJ<Om7$|q`P%%*p=L*hw z>^R^7wP>G7=ZvE2>p)oYho{<UeFxvXAi(kly4F>_0sAbp<;{dIhM(a)`XJU<W!Q{# zH<nLmN9>OKSoAhekZA--TgJUjAwZJ>XhM4)bBl?OV4^+Qg7Yx7ijoY4+-zu?W@wCN z$mJG8v&j%$GAm2jEdqR_vklIcZWG8jqN~j|_!{VJ46ou;v?;QAR?zpqqHSn_hHejP z8{o^Q@J`y`IKRAb@Z*n4?vDBF=Gc6Mh*FpISiMbb2*nRE%sv=xX~jVSISM^WDQ#Wh zAB}2^v$iDI$%_|-;sV(_ivl7(Lt}k>n`!GUv^9<IdDzuPbpCGs&I14N>9A7foD=fV zJ1T0^^E|k%%~4PD+Y1L;D6#d#5>DNwnLUrh60=H<)IumG>n)O<eyinhfIRX*#i9;{ z;`qZgIjPj0lyf2p=-RPX?1vHK2MbLz&VS<U$@!qJKzs?EcA%7^3%ZNq$e9a;p8Gc1 z@sm;up?FAd;Y%rj7O6a`D8PNNrY6nox_=k-b_w#d(SClvIX)fKK;;!nLQ;;QUf;P; ze+Grvd2iYb?G$uw(eB_!9-Xqwx$+X)^I9`%ak}o$C4H6nTv{dXK)?RBN~RiWXx-sV z9G3u<=F?fU!G(|bxTS9z&VKB0@CnNT2LJFO_TN%6(WD<{-=pSyphy`pfD<X2+G_MQ zp15q#;FoUQWd8P>H@SxW>YF6-%4JU|+v;3OvK8eyKE3oQexNkObrgP1Ty&wYeL25U z`arL2T>D+r>B#f`IGaCM7TLWPKRVj0>K2rRbyMd&GxLcV{H^jV{>ZCwy!_R~y8IRR zZJUlDmT0ln<*dwCER9<n4g(hQX-6E6k&r7PZEpAh5~qVU$eEB$n78sEFF?KkNw=bv zLWV+CLMB31LuNtNLgquZK~_Ni0$B@5r~LzvcS71BX`N?K@rzDA3c*<YS3!=1tbxpg ztcT2lY=$g{ybQS+(g$mTM#!O%nQ%B8G7oY#WI1FZ<Yq`KWFurfWE*4))=@_MI_`qS zW+>!f$iyD_MJEG=EGTFdnh!Y-vI6pTNV<7%E93#llaO}E0a#^I6!RcMA^!uJ2w4T0 z1xa5r^C6EyRzOysaX9KBt09{pYauT{Mtui|EUviE1%*(^6v#x#9LOxlT*!RLd5{&5 z$!B2*(&u}JqZu+7vIBAyq#NePRLDrk#~_m-mq2Diz6?1Jat&l9WDR6JWcm*dhouFD z=TNu|`7)%b7h~^1MnQfAnF{$W<P1pDIh>LpLm{n@V<5Lf)<CvGehzsR(uhlB{>U8@ z<WR^^$aKg-kTW3@Ad4YW?GB5h8ijNy`nu!SY{+=Xe8>#QQpj9LD`Y9;vY#*{kUcMA zNQfcr#6M$5a3{qw$Y99+xU`WBxfU{;_!8QikRy<G$d@2-rwM!gk_DkbgRM|hLVgZe z4;k=_!_fjc3-T)IA%k&eL<{6d$m@`qkiBsoI}b7zvK%r2ax>%ukd2U!L$*O~hon1H z-7ce}ki8-4KDot_C}g3q3Njz^0AvMZD`YL?Wyro*)EW4K)we4Z{GHW7YC*HZq4F)O zy?N_ugZc!>n@G*Wah=VTsvspDpML6QI?k!Y)2j?>EZPkb?b@JE<+IRk1TU)!Qj<~B zPuR0n8I*hYgf)I@2F?y+_!(jA8a^|Ho!B)7T-GDoyHOvHmO)qrQv0jmL@IIX8Us&X z<HH-)ME2ai2?>MyHVm2k$2Egh;};Ic7_p<+wFb5M7_v+_v>WH2Ha=l(kn%NuYHg6( zgqk4H&Wd(K-h9Cu-kt{Rq!C2-pc>3k{EM{#Y9pM#O<2jp>~w-F>w=UY_^fq)uA9Lf zsC5^X($^W(??8c~Z5!%Z(JDy!k(aFt5;b9>oo$_=-vzWwL6BIIhw<(r7kBmEMcr^$ z3d_K}<;%~kqhA9QYmjm$?`I8C<FQx{=OOF7dAilm?`{wVylYz&Y1=lM`}i!YSvkp9 zVA})u^mC;NZ#!nld1Cukt3mCHX-;^o(&XpCmGwa`*IIc0^?@#iulT#GeR=wNgVKSm zy<LS}v_enTo0Vd|0-y1G%laTCk{`q;;<DcF<<||WKjuC$!fFL(!T=ufdXQqKO7#HP zgGxLX?XF?2>?K0o40kVpYT&vdK8l-b`tVw`$b}`*Sv46`;fwqTEWOBo!slFJCHM_P zzbDaVAhx>^J`~B-dTcS9C%s`-3i;GGf_VMwUcBs$USW&TehTUhSy8G{5BYtI)WwD+ z-SAr)i^XJVgU^1hX(;VQ>jW`W9d8(1?B7sK)IwA=bFXTO>w;=;o>*;EwqhrKJiFSU zzKczIigCY;c2#H*6m|$pq>s>fLyol52fKDU6cg(Vz7rM~^I0{1>P&EyuxQ#~P(rDs zRzug9$888we8Ap3dxJ5|4f~>-?F_NeU7Z=<WtU-|J<7{A1gUqnIvfSgoZDbfSLx)M zffm;VA?H2l%rXsU3{C<1l=rXkRV##vkd2gpqc#Sq{c&CpGXT%sXi$9x(WcGV)B`l0 zHf`psHU{<_4$LHM<)3T}40}q)4#`<5^0}DIU+3302Dwz6<^9+B^4K>G>RN2%Bciz* z)hqdgH-nT){?wa6>N3<g9kHU_LY-7P2fWC4qFp^dg?4t-sFY^hRAcCM6hu>(Az(1l zEI=7R@#g6@1~stF;fND^qLoN2^px3UH9_ha#3#cUyBdSCkngD(qBQY8P`{h^d27_u zMG0(4H?^u8d!?I-a-;bCwwqdKU{|`SB?eY#pz>w|)gBk0XAJ6Mm0?W+&AY0)OJ&<t zDt|9NFRALv?reE?(tOfgJ=UEa5M>8_a@$+JMpPU&5`6oXuLZEKC$;^cr@E^bJJ>T} zbuU)Yi^|)3k^X1#S>&d2cechY;k+Ab6J^{?irP=b=W%!StOxtaJ>g>yR_{UOlO9yN z+><^xc&c?Kw!<@_#Ka0sRNibN{b_uTN?^a4)J45nnODL$UaZ-R%2&Op`n}%txu>^! z%$uF+ov_K9ReMwUfH&zK^vP@A_8ks5*GK)u7q_&kmA-6+FO@I)?o?2*)nEO?pZ(_V zg0d~Z1E0qO)XG3s5eN-;yzSQms%4?7mG|4^XSq0xHHWGAybz`~hO@$O)f&!T3#anV zaME8FpDQEOixKSI2z74+`&5+A;d2<XM$+fUkt)B99g9@YN3u3iF1?NPdv3ez261kf zYQ2kH8BU+2Bh=0jta?P>rV(ua2sIIKX#`cgb{BQ=^j$AeAFhv7TSl>!qk7hkVsDM2 z@{v(geSMVrP6FGQpdL?PhY~QkUab3dQ<96ChL!Wvb;I9{QIt~t#CyYV3RwT1MaioR zsr^({=J7Khm|f{k;L<Xh6DHJq^S)c%T<N)u8fY?j+*Y&deAHy()`4nmkZP>u^HIvY z4R>Yp4OGHXK9e8ZI?$DF^T>G@>nLHzZJV1)w|M08e%l7B;@K6Rf>JW>SxV+lp_GbJ zDt`r~D*OUd#dn}Y*9xrs`)v`5mAf}YsNcL#6TMf1H&1IY;V7BX5TR7_qK3aL`u)q@ zaQ`y>k<TQ%F>xLf56Czf>k<5s?@q?L1b@^6C*u}TV)*K0+$xGc>W!0q1XEv}j9UX3 zyIVYt<5b|D1^iLRoQz+e@kiZrGMx=Dfun0q#?2b|v+|^mVnX!YD$E-<W_RsYVLmJp zoXb~yWcH#vY_NL3&mDp6VaQm%^P`x!8JbJ9?E+{sv@R*Q$$L3%d!xP&TMf?UAs?H? z4==i{CTQs=7M)Ln>Bkjvl=gKM@?|`dK&KoREIp&pkbZ=5)_XF#lQI-n6rJ09vQ!kp zBz_Pa3C5WP_B`2Cx*0A?Qap_Yu@cV(50$tC953-};E@u)2~L!_9-J((4V)_R0lFD3 zO;UV~2I&%?17}El8Jr2$d#qy7mo0H0aE`=b;9P0@Vc<N8ljwR?zNDbvq~=LH8C)Ro zRB)lhPlJmkegRx6aUr-|;+Me{5?j_lQ3(Y)CdTAF@M{u(23{rc3Gix(&x5PL`j}n^ zTP5y=MQXLg{lPU|Y+-|;*xaRH3E*0=9^y&h28o{m+axXqH%e>;H%a^<xJBZl;MN;# z@np1&w%<^AvfsfS68A*cIwcOKuJ-P7Fdl4_I2~+~cskf$;<?~poh>*FRzMM{E5IAU zu@X0c<0ak?PL%jGI91}`!08hA#0el%Vly~fW894C$p%A_BQ+QW&XqU~oG)<}xIp5W z;9`mAg3Bc?0ar@A8f<m4MI_x@P}E2b8o;#@?*lhTd;)Bf_y=&K#8<#g604Z7S|s)d zw+hDbPZ5ZN!Y(yP26ssOFt}6VC&5gb`R0O+5-$UrBwhvfm-u~fFfsNo9M}Uzq}1SB zaID0?gX1Ok#Mh8Si9^7t5)T8XOFR~wDKT9u$_88X!|X{Ya-;??fb%6@1}>0zEx1_X z?cj2WPk<{Wz6`cX+#5Tt@m4K2j&u5gR4X+Y2X2r!8*Gy}7u+Z@^`uGS<=_^H*MVCl zeh+N-wsalGpF+_gHK1=Eof3Z!W<Fi*{{kB&?ujo_CW-yQ{t`!kgC&j!M_MGs{ZPb8 zJRKY_@f>iX#3kTViPwVDCEf<kl(-3;EwLSJ$&nP-pvaZDH@>6gOFR%<An_1zvBW9h za*3yaD<yswY?ZjAC(i#hl43m?)Jj|jZjg8n*e3Cp;6{mm0yjzQ0Jlio8{hR>B_7xd z$DdtN3`c_wiSGw@N}L5|eY$dB7T75960k|)HDG^<w}68s-erLzQc@fP$4Yz-94|4R zb7P4T(-q28i6g-25)TJwN}LAHme`U7MUJGH1<sYY2%InRYH)$Xwcuij_k+tN{svqr z@fEPOi!ICzmoRI(6l@T<R^ofW4H8cW+a!Jt+$iy4aFfI<!7UPR0=M2^od5SgVZWj9 zWaq#g68FS6icX1R=o<!p%cSXm&K>lKwNYZaR%?=&F3tK&{0cZ&=TsPYABsp_0X_hZ zmG~SuUgDnkN<-h|=yDV-!^6Re{L)^tg?{;QF608)4%E_5K04FlJ<wG(W@uZ{VtpZg zztb0fsVFRxm=^1WU@Rs?H?de}#S+t6y;NeX)h(=CQqW?(LSkC1S4vEabzFqKu>)GH zTP3!Et0g`F#x<-P?P;yPS>kiFn6H%-v{-Kd>jOoLb(_SrSZ|b=7VAyY_Ow`UmN*IA zA~CJiPf9$Q7W1u=f)?v-64PScE-@|EFGx&_^$v+?v3^<NCo#}ZV+iGe@098-FQLK1 z+BaVp7k6A<!O^dM?PD(9Z0jw&AFO{dRb2epX)t{^$B;9{#h2N^^Z<vGeb{ecx+})X zp3DVbrlq$H+{p$*AwSYx1C}3ZKMI!L6XA_-->>nV`@P-u2WOW<lgH2OH!ICt`8-BR z<Z++J^p;=3umj)F1C4z4=Yw%i#d<W!=SMz|Ns=G<y^b&K`V|3sqS=5i?eYW0K4AG# z^hx0BJn4&=B>CM2L-5sJeol5QSbmE5aj^W3fMZ}Q-|~f7naYoRVYbLG(D1=k0R1aJ zJwh7-mY>ca4W`>#^e2+D!SV|cW`pT!2t}SA*8dN96W%m{(}Q^C_cyToVu-%}xC$V@ z`r;W}b}(Jz3lEH&AipS~6;~hRN3Hwfs)PKZkQHG0B>)@1d6-pnKi&q<<RJ&mAys$+ znXZ-41K4-rK^pl1;tH^pPd#X!AU{?83$*g%@%JOI<>%BNr##pG=&H;Tu>4%~C*VK% zk%KW7`RV7icpPXY9!HA97UWA`JbEj?fMP6I|0AnUfaO<3YyiveQ81Zre&@$nu>9=z zVKXM-O1|JwjF<enhUj2C>(6%{iV2Zl+VBUo@}t?kLMYq44#&jFuY&jhTKNI*OJI5M zJHa+S>+onV`IQu34WQG&&ckLe`BC=Bfe13cba-@#{8kQnLyG*W1T(mXPiPt)BEN3n zIJENn4IT``(aX0qnM3546zqa_8XiZVh+tB{hlgW{!7nwLm2B>HBnH<Qla82^<Tpc9 zMZi-e`$%j-FCf?smft3G1uVZsrF$d}4!-_~*-L)qOFpy}{K%2fN%FJkozTiJmT|r9 z=C@kh3D*B;_|w6Q`GlipFZq2S$DqB)7aWc8l3xn4Xb{DHC)&#I&N&aQm0vm<6C%Hn z<nZk`zpTO#g|j|SIu;W`cOeI3&(!lp(8_N{_!FGPR~#Fyyw8suGl$48>{vGV=6Aj9 z0Kde&nn&XY#|h2mB>6=)nX#DQSl8<YD#7ymRDJ{N(-FNvf#GPW<y)F#;^fzv?1ol; zA4c~(a6*(`H!&DIh<hE6QF3|G@fa`p)f9n4u)N?4j>p8w?``=IS|*MBH(>c4A!$SL zvok+)+?*u8kf85BaN?C-lTZMz#dJeIp;E@}2Fve<_!WE)PdX8kB)=7+0gtN7ueu2s zhGiFC;(+Czu<}1(`Q<H*;LUu?i5TTse&&Q(ImDGO%^?gogX(!3JRF;9f0Rr2fGc>? Omu4lrZt9n>jQC&bKa%DE delta 18831 zcmbW9dt6l2+Q-+1*)kx?AczRa3<`>-f@O({3Lfd8p`xLpK`&lT?`YJaqVwYMGE+i< zNK(s;3Kg&6rJ^)L!=k+5E%6$b=v0(LS_exDOP%-ktl1lu*601><?~^F`}?f5pY^QE zUVF`+z4z)3hF{kkN-fI1mm_A6h<%4K`d1nLOpGZIVH>F$WAUHt$yg-p9kg^0XbSSD zAT8zPPltG)evPry-Bjagt|<Q#xV!Gq&tnT7KQ=q$mwgd?W(>NwA%FMwWy%h3+kf9O z7&~I;o~#oaVrH!6w!s*naad2rJ_Zlcc%VOHr8f-5g&J3RF=kIxjmyCPEQI;?U@ZA} zgVCs^-`k0?=dT%zff~>1%a{*S+x&kX!e+6;PO8x$9Qd(DIEWjqa|m0~nX!E>1|viS zxC=E-LcoYia6hy{!p{a{7mc66{!e<h7ceo1v48UjU$ZU!s=*kk<tj(^T3^*jZA;Z2 z)|;^v@DhneCSC&fQ&H_uE&Vg>==hszj0aOS--f0jD!rN;yP=CZsz%gS?Eh34qCzdG z&nW+cJ(!I-P~Z$0x}bn(kaLmL?mvLjwRWNO5;RyU+EMS91Kk;`kOHvwLcpjHx=R%J zEGo3vs2c5>|B=u~uN#bL4Us;}*T&eq4(%C|(V!MYNHw4cH=<qbQoAShVJr;}d^87F zP$SEj_7<4r$JjT&7>xb3^a^l@6!0tyV_EPA5&mBp9KzVoK5AP;$Ivl9sO<rHBA{xd zM`;c>Apl=~wUgPFjfU3Sp(E<O#bq#3{ivhfL&wBRVSkJUc>^7XY$E+(4CE+;r4N_e zDD`Qi>+Lcpn6ZUuN8SIVa2zvIyH<v?dzljj>OCEYc8dJNU_7pklxB3yHED+|1dz`^ z?QD)%iD9o-rWKWmMX?xWBCJ0!upZ=!pV@z&mulp&Qm6lbPPg+xerDShj3s?+(xBQQ zm9PjQ_ir{B$7y?g2wsPtLI()nn<7v}X=n|@*enEBYU!Wg5Dt?3Wg`DJ3{CxzpX7_N znulW=`7!TYb^V0K`av2x^ATpg)aI)(G@`-!P$~eAk;d}7;K67N{orUqpOvG6saimO zh&&nr>V5JOjnzW5x!!3<un`T1dW&quP|1*5_X-?M!q7?A98!CxVbJLLzxKeS1-?&9 z{|^q9Yd70F{aFO`nA8H_+M;o|qek)AfewyOXgEm`PC@$nX!8A@?j|3Cu^o;8^@AlI z-TO0Yn5z}I5&`5%qxiY5$d7cr$ycEZN25Y;Bf5~M_y^8Bi5>LuQ-{$u2lJCYADux? z4+Jn-+hH7PFi$#|{)WJ+adcv~6%PG;8CxdJvYFsC{;9v&){0W}hVO>4TOn238y%C4 zz2F4nvopJnurHvU(f`8hZiMr^G}k@_`*irzt2GX#$OCCA3d)iuv1ct9X+5N&ItBH8 z07DZaQWW4qN7rl9G#wyy@Vo^9z%0@i<9M>8-G4B+51<_w+PT`uwXp8&7SdZVko16l z!j7@hNHv>g!ypaSMmWgAEZ5}i?)U|0$3V1uf2}|(S|D4>f5pPsN~s0D!LgQs>84K0 zub^K#^TdFF0ammvRpASGPTp_O+7#7*2Ln9nWA!#7Sk9jdF#FeFKqP5?LeSq9{$7CD zb_1p9<NR$@uOF(HspVRVx#(VLSUrHrZZe{WOVM1F7>4=MF!aFi$ivB8pZrwJZG~ul zofn`x7Gjj?vzA|1b|2e<048b<-<R6kg!bNn4mKbleTV62@7rG6pAUSYzb!05HFiKp z&}u=&q2C|ug<%~_iogbwqPgQp(Vy|8Xq5!g4wz9$vD!!`Mf0VQqA^oRap0wqBJy-n z47?0dOgUMkoiWhn&;~}dl@vA2CPf8uNHHVik|Mx7QVj8hr08!uDOxh0^mDWjDf)RC zDGuCHQndUE(zh8aBSn9eli~=eAU%yKZ3As!F03R)OKl^?`KF2#!(#_&0S;VJ90E0@ zJ+V+BU4~^4X)ng=NHMUElA??1Nzvq|NpYfSAdP11JSn=^L2CPeu_oH+$C#6}h_NfA z_uvpA#Sp(nI)JfO(t(WKB#l8|qszQv(Y2(57&DUIi_<gdeK<*w4#tr}dOu^qqz^C_ zMmhu|7i#m0!-<m&Lou0?4#OBAJ%KSmIvitw6q9Zo=?IJgQk?1&NJn7|kS1UZkS1ab zkdDR}ApIN00O=Tv0n#Ll0n!IC1`^e9qu1Xthms)~V}Nul#sDc!q&a-;--2w}3mGfs zH(f4Q!O2|JITEAOc|yfUHnOHaa%ivX{zx>5D=HeFXcp#*>W?ieZ8kfwn=2}on9`c- z1rH|1U}Md-g5!uWXjpTV;CNySafQnUOH=|GqRCJyGK?Wc+p}i7;ACPnC2P(RoI;#H zoFzDwIGH$Aa2jzcaf0A<;&kFz!5PF^#1@+{WRZcY)NB$whd76r3APgF5x3liAu5~L zPV5kzLtIE)FF2RDl(<%K9&s6QmEeWM6~wj*VX%{-k_@GS^NFj7?Sc!5tBG?2FC(rc z&JtWoTt}QLcm;7iaf0A7;s)Ya!R5pbu&vo53>9Q>QU;UY4a6<POmHP}D{;#$YJqLU z41=ZFA-IZ|mj2E4f_D&`h-(E`69*Gl39ccwD43#}D}<qz4AGRKRPa9HSYo^2I^sCu z9KlD46Ns|}*ApibrwTqzoJyP^xPdsGI9Bj^;w%MIRI^1G9AvOk29w|>;v8Zo*h!p6 z+(Ji7)D>bou|seRaUpTN;A_OC#I=H3iOYzq1m7gCu#urc7-+_^RuY#ArkTfDMQj&L zGmy2KI7cweMAlm3EWtD*S?h>X1=GxAttU<pOf!_VfjCw$%~V#KgA5j7pef7hBsK}A znakQj%mmX6W^E;I`GZ=3W-=?oX{^~Hm}WGqk+@zk&1_Z^ajjsQ;jF>LRf1`zvso=< zs1OF4@~qLsrGjbdv&ItJ1=AF0jU&zxOjDsXfjCPrO^McI;#9#jHCj`N69m%~X-y}N z6^yBJ0tNvYEW&^}6JvnbB$%d7YYs6JOjD>ekGSP`F%J;ii5-G-_#ep;UbK3%I(&Jr zvAqL*vGQWiCKh)!VWWBsE>jFh;F)8CgNa>%uO`6aP8|B>)88LZeR;{)kgmzF7S!f; zp(CDQ!~38iGMU$m?Pn~*(iC9lKaCAnI`U5Ae)728gWnpLsBk{&p*Zh*agIGWCemQ$ z_J_JE{rSp=@{~8ZUrLHHfoG(Q#x848x+zQfhbe=*e~c4#NkJS&cM<j|F0i6XMTSV{ zd~O)u&%6Ms<SC`IYAh*;r(8&5s3%!x@bTjZC_a3_`0fG!pd7EZV2l){MH(q}0e@@! zH1k$02--cQXdYs3=pK4Z7~nCi7hgC&IAED@K?j=dVkQ5d3DeA9VU5x5B13bL)vCMr z5iaKS6fWKrE*w%6oB7`-PBVAB!$r2{qVzZ2#h!@+Jl=>9F5VX|3MChNc(>GP=HbC| z6bm&MwZG~vN>T@S9JUA-2ZW1S$;B~#2`(PT;;FqN%QP1)*L4?}lLmPFZWb<%3l~+A zi?jU0Nz=@4Vm;RGqFi%fX>&1P@&J#1!bM<%a8WP0_^vZAoIKTh1Z%r?AKNq^sn_%{ zZo$Wla0(;vN8zJF^6_hDp8oJu^B-6zw)?2jeAs`{eH?pufXC7vqADHyXk?U1K1{q% z+BD@JZc8&O7QQB}m$H-ZO`8Qv?<swh3_f8>ZzYbuHD!=8jy{7^;i0%7+dVQ}Co%MQ z;E-@W##^SuW2XUABa}z@#HrE2aevAfOBoOH<x}I4Wx&*6pVSr%uhycUn``(_Q_V^= z?>w!iGMf*Y7Nrd1)259!X2XX#xYkY^q<q0oO&e@np_}~DM=0I+`1Ila(Lt)Qv3x)E zadFbNMm3tho?hA^A5($z0FRkIEjSh#E*xrH#gKy<YNzgF_4IIMKi@Mw+~aFM-Z=f9 zAU%clwO@AU{??C&&zRc#1*}Qi56LF1V2ca#HHS1S(81Um#NU_^;&l@Pqj6OpUo)d? zG{ZP+RO8$cWZM6{kBc&^alCZ~59;j$ubD@Sla7iBAgPR;4&-YeX;cnyTgEu$zxd&d z>B>KOWadQW1D=~XMLEO|XHHN|y!)e5lt}*1NBetx=g(I?YUwotUti#;F@HZD#b_0T zaPH8T%@-Z3{a&lhul38n$Q_S%H|)pxgNM%y3QC}6%!u^GzX6SI!FOcmvG4fHbAou* zOz)uNJCo@vw6m-&*_3)`vW1uYS5#2?oksdv?;Pmf#|4*JcP7)t4(E*~l-SR!lapeO z;}d5MF}~Y`?efMJUX(Q&WASKKIEKpQtO&)*JI(CFSI_kG5}xb6Z4379Xee0ze<wfc zPNum2D(t@P$$u60wQoe&#_!w0R<`gLXN^U?CuRjD=21<{_tTV2eWGgY9<uFBbBvts zXk#8dyN{B{Q)l<<9`udNmH(Arei%-3cXTl4vm+Ec|6q2Xu6NfTR}-8Tes#9xFH<h_ z(mCBYYTotknp3;&Y)ien;*|OvQt$D4_1pG(58=Ph8EcC}jW66lov99&Q{`Iy|69S` ziuLe_(<DCW@MyXS|Di>viVn2f>M!!o{t>Bs$}j!nzzFyGa*TE^-RRJ!MVd_uQZ(yw zSW5~-9?TEz_7C|Fk44%hQxU=LBJ@2*!dBE$^yOufwjV<2C*{-23zG6>Xd&Nwp*>{X z77iQvYaWN=aY<6SQk+yiJgLI3yf)uf>sQc#sgkCyAGE29CMaL*Qaq&`CujpT#DGS3 zmfw(t_(zSM2MvNVwVD&tq)INvYPLX=d!Jo*?bDaiZ75fF(p}rZ-Kc3UcOTq`N>X$) zM5AVC)*i5y6o@>cQFVI~>`s$N)J;F5T)l;L+XdJ}3)|4ddP#?+jJj>7p0VpkG!d5+ z(Dwyufm<j;_o~~PVbl99P|J5l%BR~tZO_-QAXv-zHf0QV4+P!zx|Wd&Y~+LIKG0(- z;%Uru6eszTEg5uvyyN_$hQ_c+CnkltLv3RwpzXJQiAFO-CeijNbY->i06Jz&q_LEi zlxQ(ZArC#tpPCwg`P<RhBN<jo^<)ae?pZWx)+0xCAuj02yHg77{_S~DP70<AXB@6! z741zeYG_mo-IH%{qgv|TGGS+5CoX=GlZP=NsFWwPfpHuq6g8~3!zT*OCSO0|qD0YN za#hBc&+BF@(;S|2;tV33+0h&1ELU?@tT|(-Ir`;n(ZM1#$<v6BRsj@i@2^BLbYZ_F z=|i&YLn>{;4aM^cT$ChLQOa(F)b6E!nB>MmO*<yigvw9=z6wu#ysJ%|B>r)Z{0Rp$ z$U#@pn{LBcsw5rL6Scw%VW_UYt6RtomH#+mElK)>HccYPPxK%wM3B{#(t*lw&rUXX zh3?l1q13*2$Eh_K^+HKfH7x?r+71-xe5Zj*6F%3YG(B!%c!jF2#Vz7F0mGl-euVr` zsXI_=Nzy^uoG$8cP!IeA5qK@7Or%2GfxGLtSSy4QpGG6z*}m^)ceRNX!*5@ZbK&9& zoa=E4!&lV1TAZTzkiUwP+)yF6L{A;3&2=K;<9e92BFs8UDWfv<_7y%4YK2hhM`*M= z;=H4M>$LV2OPK@_xA1wt6t^%$N^y&LX2bBO_O+28D)m7T_!qQ!S=-b`?-3DrJ*9k& zD(USje43=V4O-kd6}ZRm9dXxdaaYj!x9BXjH@UPq7n4tIK#M}2fI)9?5zjUl{uDU* zp^%n~A$E~A2aC8b>J5Hg#N9wC{iL{s%jsGfMK+|4M=Mf;|8?LETHp=}{Nx!5T=+bQ z()756VJ8fF+#;UqF#IWQazi0Ch_1V2qwPha4wv-6uZh4Nl=2J}q77)_bDdTQC9Xs( z{#D#_(C{Tm4lQo_>Dd`7;ubDV<Wg&Ihp4{`#f!inDM`9Sjz|7g;N*uABiyZJQ1Jr* zP37di1Rkmz?wWEag?4B`n0g(A`-@t6l&GBIBXh<m1Gg{9>Fwds?d<lgPj^z3biV7E zp&jDU4$fxY`poCbJbp0upz{9q<<I7;$`jiI7j*Mb?&tj%^-)Ifv_*9T3<qcpdX1`F zoFqzp8olW}R)?Blk>%V#8#sum$T~iDahkH7uU@=TxykRh%~Iy`GTY>?Q}@xz#J8wk zsa5Jfb$C+V&$!DrSQ*a;+ZXu0x3_JQ*u+1whYY=pLeT;hIE6Tu??>X4BXr$?f+$I< zq213Nq+O2CW@kizA&>qtR`Slz4Ob5GanJRzjivqY%?UXsAhX^v!Y~2`(J!6c$BMH* z<!~?Q#KGeX)p9gv)7)FI;Wl=Fr@@F+JCUEG7*VDXV?qShmgy5zx5U<lpCXX_+T3hp zHXKBaP%*LC+6-T0Yg<dX9U8uq@)S^cXgO_7cR<VXJF-wXkKN_sEQ)F#gR}9$F}6s9 z0r%PHg^jq=(PGU}Ph=rS;deQDi5$_DxRRuOcUJDSJEw17a*pb2hc>=lx--LB1hN$H zrsqvz_h^~$J(eaeInz~cS~Ae~fE40mpWWFWXVHvu`y3!PMB+Cc<{XE-wAi57e}Wxb zNloqD^azsC;o6F#@l#}=>k7jWY;jUkF<Vk=EDlCL@;ysB+N`xSP)@r?coRIjP7%WM z1$fp+g)khWqA=K-FCz$T9K23p{pC2YA&QbDF}bWjlh7veQI*nAZEdn!%Kx6<YcL%l z`Un-3O43TYj1ek4<V!h(;TgU<KRBcx%u8Ane19+a{uat;8+8%9KELlAuTR_oUk+_g z`Q>AZMI1$y2um+k7+!-xKQP2O`UNS6Fg&a02%@2ZMw-Jj3%Yk2r<uv}()J|<T|7eX z5$*chopo-Z3=zEQ#fYT;AYCt9m@mK}3hz>q^c&^tfsR|U&+^ElnS0&KGGX<i?$l<S zKeh)IhADBWv}itudi`n9Oey55W{-O?d30B4Zk<Fs7(e3kirzA~K61JE6HB`f^w<r@ z4aHMT`r)S3Ai@waj6|nAg%b3o=?6SxX}^J6cHOnh>v!$aUD_p+?^rrlIki2!I7?Ap z;kIQ@@Rnsgyzj@al}oOKIp5_$C6hbG;nMPwR#$GoM{)d_lIRYV_)XJUvwdwzl&9x^ z&-^z(QkuaFUVCPH*~=}4?WHU86kE%CE|&od`ScSmS3Gnbw3dEs#>t=wng^}NEM|vZ zffhjN*0wSz9g-WM3D7EN2DAp62W^3_fZm2yK_hYQItm>PbwX!B4Jv*i%EyKUtN%)9 zJah*%9a;;`g*u?6P+I42gIb_<&;)2Y0?vl!LTSBH3N3|hgVO4x4q6Lsf;M4EWW=xS zov}=egbsryc;OeKENo=J@FX-3S_oYMeFsW61AYcQ3UxxA(1BQU)7=hsXe4wgGyz%( z&47Lm&4bcc&K1xNU%Om8pjFU%XbtoVH0~P&Wb?*tQ`m@vPJ$*tbD<f~h0r``A#?>a z^;<ZCTE26+PDA@cTcBg1zL+G_pwZB$pktsMngv}AT?pL_ErafW?to@|?{e81u<<-L zu0dBrO}OFUBWQo<0cbMx2k0DV*d?5cps~;j=u~J8bO+P{Jqx`F4RpF(!5BMX(817H zXe#spXf|{lv=Ex+blF@Ru#o{nSSS2w4UL25L(`yT&>ZMC=rZWaAJHYypr6ns#84;k zf6yhkGh-#x0_}$j9AltcpjpIMk>7+d0(C;)fa1;-_Rdusng$s@gP{z17P<o(^|Q;> z09^#VNp`3OcUm+-<Dovd9-a;jhvq`#p{3Aq&~4C1p>@!wp-s>_DBbZId<_+ahC}JT zyH}#Ik%5g(&^%}ZbOrPZv<i9?8io~}FW>y;J!-{i3>ofU9-z{FT+8^?H$#<Yd8hJF zC7aJI@5U3#4Y<mVEp;2xQ?M1P#Pdbvp=uiT>`Z&|D&%>PpIFmP4Z{g7m7hSK$^5-l z-MrJlQ;{!1<g=_Ys2QM0QQk77&EbiwLX|Y(e;;@Z?Zr!18I)w&Qg6bRf~`;`6`ujz zu`1f9X{*acHwqbId6(4#l|Y`jI#hA+%+;Z4%PE(uizupewL$$J<S&Z4iG{;?{vPt3 z=f_ruMx90q-NkD7?H<*b*x5L`z@w8oeTm@(%trX>hNVMK^G<6*)o&1PfByBl0G_bM zp#BVsB2u>@$4x$SO{nq{Uxa*D^b}r&e7}MMMaf~9mwrOigetB4M|}3=-W6sw4Qq5C z5obb$K^+b1A#!g+?lF93MJRrJTZGTi{B`7ZV!rDra@QkwgH9?<ZFT%A@<sDbYeUr# zOl$6Z32P0iH>juEEApB6%(bCd124ko75@5Kv$BHk!si(7zpfi^LH@y*%-m6AU~0P` z_E21~AFwV|jY5j>#Ix2J@EZ#DP>&){2!9TFyula(>x>;+F)dc2^$e!|JZ4i@M74H( zApdb)NA(4G2o*Ijtv3jAm$VRb--~?E`cUOxd^|oE^SSHIN-2M8LpNT5oQ2r0H$SmD zp!0du;T2FO`5(*=uMbg=>S3uVnEFm~WkaZPjrYK(LpVs?fW|@|RWw2k@_Yt+sIr@v zYzS2=MG9Q=+6@MuFR^ObjVc)u_~#o!)SM$YL}<5HHX3kZoLSjT&Bm-4$0Ig|DzW5U z-2jdj<>Vlb1zQ-Ms|}sEp-RE<O?TTFjQov_A?gKOfwPMs!Z0I_<|p7}6FzZxfPFA` z+o_J+vdQ4{B#2WTU*iKeVONnz$5_@T3>3`TZbwJ4V?X}fCK^!BZ9?;JGENS`Tt1O@ z-8P_oP(S1z)E3xV(E>@xlTLZ6r95Mos4o~^^rSptsMUNuPaXGpD@1(<C=PRV8Rp>% zp7>TMemcl}D^y*F6w!~o^euzB0n|y<lIH7`{5|A*3Bh&qruq5>q)rhnT4_*g0jVM@ zP2JacL}jSY1=w;ava6_gkh%>C1JUrIN^iahpYhZQyb3!HZqmDyruV7*1X5q(KUUtS zH1LScp~_J{V)KZ=<`L{gPqo~Wz3QoAyN*8jOPhmifHMX%eT{*rGAtOd{f?^cRauQn z+uw=LtE#%P16$sKEFX7JPjz5N#kNa)zUrl}H?p_AM*iW&u8Hk(Bjr3UKEE@nzj(7- z#*s(7**<UDzT!=(YkcrIf_>nl?(N8q`HXzCBU{mtwre_4@_)o<fv?Itu~oh!FZ;44 zv5i|&k@^XJ@~kaEM!<4Eg2h{cG5{A%lz!c$F6qKb{6~J}&rbW(_Dz4XzuSdA_jgfG z1+d00Bew>y4FR-$G=S_b@%d_?y1pxWJ8&d^sl6t)%ex|#w`}P~?X@#ly%Ef=2YX<< zwtHuMp6RZZg|HPNWVsQd{uat!4pp~?G8|Xfc7$H^Lc-EWHGeqU6sfL{Vn-s?ZzI|H zNVPbM<wvO}hp^M)^GcLj*PG?{Rx5h5*L&0Ur@bllw)kAxN98YU?cxPPMYMX7ui4t| zF>LSe+u#d1GKBnpIg~zIhN|a>vBF_tJBG2hhtc+#VPwBGjB2@e_)7}pn|QTm1p7VS zXT=CsGJ>{07(vP3j8I=5$(D{(ca3E4kHkcCawNa?c8~2txApQ+V{t3!<n1%xiBWJu zeEHoWIMfflYg1CTkFEMd#hKB2r&*=jf3x|3oxQy2A&tTkI%iFX-Thq2qMg0e!cf@x zMr_6Q#e*sQ81mBNDfFm{iBH-Qz<odTRq5W21m5GrUfy&+M>5h3A}s@H;(-<Z#D~4q zFx&?e#$Uyj1zQ%r2l;dHV@d9Ndj6ZpPv7YCc&F+<Y98#v`ONA7o>FbX#gNSEK1w++ zz~@5#cJ&`N{pMv)+`LSGR0mQ|hE)dss1~H2wA#QQ)rHg(cYV_z)riy+x9-v()r-`V z;XYCPQLRWl8SVqXAJvW2li}V1%!lGx9Jc{KHRF#e=w|#jjX$cUo9W!#e)F;?Zl<L_ zJD<BdCPFV(L9w`xq^(#5#YTgZ`L5k&f4aQ}Cw}}k5yECbO}u$`%%FHJB--~mSo>>L zQc%U^v~OE_SN0}2fsg&jEPf@?ebvKCzlZ32k@f;p;8fO4=)dqh0-a<$uuklcjP#3% zJKc}bZImAV`u=__8SE+X!(am#XBGJKV>9VaIHP2E5*d6YUH~>pTnrAB_;qlw#BYJa zB;EnGNL&k!l=vv!2^TFHzC?!p5?=zxN_-6*2i9w>VojDHaaVA%#8Kc>Y5yVMG>H>! zFr-Tc`T;3J;%VS4iD!c6Nc<$&D)D@9w#5119Eo26=Spl_1w$STsF>)Jcfs=|{tWzY ziO+%;NPHQ*5Uls<ZLnS9j#y*nOB@a^Y-bxA2*a{=0~-l01?x>b6<j9qKfx6e7lJD# zt^ij_ybD|{@kwxP8{7O?6AX2220!*2xL#r(RINc`3)mrX9M~ywD!4`B+2B@*7l2t8 zt%5cj1~0>4)D7TGV3Wkv;9!Xlfh`hW1V>AJ9ULpM56%X05}UyZ8skn(KQ<7CWGTZ4 zaH_;9;B<*Iz*!P!gRK%T0Ov?t49=7IO|ad~HZkZn!%!$?s0Noxd=OkF@mX+%#NUG} zCH@s$C9#SLt6Jh<aIIh*|I`A5V5pNai~-k6JOkVy@jS3Y;ss!*#7n^~62Ad%mH2%y zqnQ)sivae+V3elP3t*GPzk!1#_QML@B5?#bTH+z#Sc%7j<0PgFI|*Q$ewfXJAz8{W zADk}nQgD{UtHD-@Yrr`Yp9SYhd<|@uxC;s{3{Y)t$9XIarBa3o;4+D`z!ehbfGZ`Y znp8==99%8&8gQ+|?}F<BZ0!@xCot4Y8BT*6B>oQUkoXU<Q(_-{?P!rW7~CpxA217S zFE0*kv`L1CU@%EM8yqb0GhmCv#o%a(SA%0E{tz4|aXmOeVkg*^EE#UWkScK(d_hZ> zxEDA};`_i>iIc!N5@&+*Bz_ibm$=vm=l?><uof9gCEgA$lXyS4LgFvLl@k95u9DaV zu9mn9zSPx9+^ZvwzdFe<6dCFzehAzkaR%5S@l#-@#4mzdBwht>mH0g{>)JjB_S#@D zN`_Nllf;+6!4flkrL;&)mlLBU?gNgMcqlkd;uLU##I_6=k|o1a;8ckV!08gd3C@zZ z3T&155I9HTufTZ{{|dIZvyJ)U3vgk(f%OBIO8g+WOyb$#3W@&(u9Wx%aFxU>!POFP z1=qGQ&j0&isB1I$u}k23iGA>mqCw&qutQ>czS=1<U07|An69I?O8hFA;ipDSA52Nc z$>)6-jJg4Q6l{|C5;$05AAF^uZ**aBK+EvnU<>bl&}^gMW!ww75Vi-Y^n;Afw0OUb z=f=~ty=bvM1K*1Dg<mo@9+8+9>lt7yCYZKZXIT={T78bhSgYHZRWi_GJzHX0tmjBf zi}hTIX|bLsaV>bE#7Dt)iD|8#FYzT>%oj=qTCA6XMMI-qX|Y}=F)h|BB&NkWE|Rs? zm=^2XBu)fZNla_?9THEY#eB77pv8KP#I#tim6#Um`y{5tdY#0ySU)Q9JaqJAZ+vmq zx{x})UP`yU1j8`xo3DpQC$<r+f2Cv|0c_`8d<d+6F;zT**+npYM~8o=cmy#gn4aEn zb0E78rW@zn?8iLtWm<aGz*II62KjmJO0fJy`$@3;l869&`<7phuoS%Skin=wL%STD z&b<$tm1^Gia15?nOg|jcMSdT{9(-TF?Zw|XJW$!k4<Snizi~JwQGWLKHon5^*9Yj~ zW&^&$%g-7Ig5{^ur-HxXbC1L%%C9-N4`1}<M`p)^<p+u%1Iw=qI0d%zb4Scd62Ec8 zY?I%m5r}I7bhU%pk)EcF0Lu^Rj{?)}DEh<6Sz!6y2v39a+V;QdS`1!}mk!X`h2E8L z9W1{iA}ko!1mw41{1X=;uHEEYj!u{)zcaysYZ3BO*I~E@A-^-^Ww87{fQ{fZ%rWkI zTy<>$r}DAK%n^1xj7(QfD4OATrbd32cm>$bi;kIl$`4ro46FR4{6iS!tu6fUv2cI+ zN$8JZzQAuBi|Hvp@4Olh0R0;eAPv%r48y~=Kd5}%=i!O+J0Lc~EWaSZWJderF>1{C zBGU0-`JwLP;Dvni=P~~BYZwMt=mL22=P?oT%NTCJDnD}FF@lEf(Bm<K<X1d=0IU4q z_EoUllda$i{`&Dz{_?9P&h?}jzWKP>Uw%YAx))lEcds86A-`>d-h3jz0>KO}<nxhB zevQBxSmpN?JRF4slb@?MN60T9*bD1}XrC!s`wZ=k)dTN-f_`!tdLjltj?6t_PL$vH zP~HcTV&ukS4|<<M4Oo6l%&%bi?I<0hG1B>=6J~$;RW5n3=J6XRM)j1RF>i%ceqoGv zUwT^wpK&rA_gl?BY4$%Yz3bu>5>N8YCu98OcX}-8N6p)eyz;AWF2kD7yPt}QkY6Zr z{GL1CHDTyamnP>TR~D|np|?c^m%u8&*x+|?0^fCNl(LfFIAxBI-?y=J;GM5(*#mx> z4?R6f`Iyf?ZBCTmLz5nhNr|<;Euai6zYgU(Sf4`Z1qcj7e+55xI%bgkHj;gGKuE8% z=z#kk<hMl(1b5~`&&1&8%eiM_{N*=GgxrUBW$?{saAgtorDXth=q_04o<O~WzXHpz z_(&Oy-?O>*S#zTNK7p|NaVC`Ba$pC`Z{gVomS6Ai3%DPjdp0IfeyKzC1DI~`mIAEm z#o^k02rd`m9Sr*6`#4yB_evdj89#S6Mw!mNzc4HB^S)o0JshL?xG%z$JU$mDpY26o IygKZE0e#e_k^lez diff --git a/NPAnalysis/10He_Riken/include/ObjectManager.hh b/NPAnalysis/10He_Riken/include/ObjectManager.hh index bc9b0684a..9f84a86bd 100644 --- a/NPAnalysis/10He_Riken/include/ObjectManager.hh +++ b/NPAnalysis/10He_Riken/include/ObjectManager.hh @@ -109,39 +109,39 @@ namespace ENERGYLOSS // 3He Energy Loss EnergyLoss He3TargetWind = EnergyLoss ( "3He_Mylar.txt" , 1000 , - 1 , + 3 , 3 ); - EnergyLoss He3TargetGaz = EnergyLoss ( "3He_D2gaz_1b_26K.txt" , + EnergyLoss He3TargetGaz = EnergyLoss ( "3He_D2solid_1b_26K.txt" , 1000 , - 1 , + 3 , 3 ); EnergyLoss He3StripAl = EnergyLoss ( "3He_Al.txt" , 1000 , - 1 , + 3 , 3 ); EnergyLoss He3StripSi = EnergyLoss ( "3He_Si.txt" , 1000 , - 1 , + 3 , 3 ); // proton Energy Loss EnergyLoss protonTargetWind = EnergyLoss ( "proton_Mylar.txt" , 1000 , - 1 , + 3 , 1 ); EnergyLoss protonTargetGaz = EnergyLoss ( "proton_D2gaz_1b_26K.txt" , 1000 , - 1 , + 3 , 1 ); EnergyLoss protonStripAl = EnergyLoss ( "proton_Al.txt" , 100 , - 1 , + 3 , 1 ); diff --git a/NPAnalysis/10He_Riken/src/Analysis.cc b/NPAnalysis/10He_Riken/src/Analysis.cc index ca47ae099..e6b193ec1 100644 --- a/NPAnalysis/10He_Riken/src/Analysis.cc +++ b/NPAnalysis/10He_Riken/src/Analysis.cc @@ -141,30 +141,31 @@ int main(int argc,char** argv) { if(ELab[hit]>-1000 && ThinSi_E>0 ) { - ELab[hit]= He3StripAl.EvaluateInitialEnergy( ELab[hit] , // Energy of the detected particle - 2*0.4*micrometer , // Target Thickness at 0 degree - ThetaMM2Surface ); + ELab[hit]= He3StripAl.EvaluateInitialEnergy( ELab[hit] , // Energy of the detected particle + 2*0.4*micrometer , // Target Thickness at 0 degree + ThetaMM2Surface ); - ELab[hit]= He3StripSi.EvaluateInitialEnergy( ELab[hit] , // Energy of the detected particle - 20*micrometer , // Target Thickness at 0 degree - ThetaMM2Surface ); + ELab[hit]= He3StripSi.EvaluateInitialEnergy( ELab[hit] , // Energy of the detected particle + 20*micrometer , // Target Thickness at 0 degree + ThetaMM2Surface ); - ELab[hit]= He3StripAl.EvaluateInitialEnergy( ELab[hit] , // Energy of the detected particle - 0.4*micrometer , // Target Thickness at 0 degree - ThetaMM2Surface ); + ELab[hit]= He3StripAl.EvaluateInitialEnergy( ELab[hit] , // Energy of the detected particle + 0.4*micrometer , // Target Thickness at 0 degree + ThetaMM2Surface ); ELab[hit]= He3TargetWind.EvaluateInitialEnergy( ELab[hit] , // Energy of the detected particle - 15*micrometer , // Target Thickness at 0 degree - ThetaN ); + 15*micrometer , // Target Thickness at 0 degree + ThetaN ); ELab[hit]= He3TargetGaz.EvaluateInitialEnergy( ELab[hit] , // Energy of the detected particle - 1.5*mm , // Target Thickness at 0 degree - ThetaN ); + 1.5*mm , // Target Thickness at 0 degree + ThetaN ); - ThetaCM[hit] = myReaction -> EnergyLabToThetaCM( ELab[hit] , 1 ) /deg ; + ThetaCM[hit] = myReaction -> EnergyLabToThetaCM( ELab[hit] ) /deg ; ExcitationEnergy[hit] = myReaction -> ReconstructRelativistic( ELab[hit] , ThetaLab[hit] ) ; X[hit] = HitDirection . X(); + X[hit] = HitDirection . X(); Y[hit] = HitDirection . Y(); ThetaLab[hit] = ThetaLab[hit] / deg ; } @@ -244,7 +245,7 @@ int main(int argc,char** argv) // Plastic for(int yy = 0 ; yy < Plastic->GetEnergySize() ; yy++) { - if(Plastic->GetPlasticNumber(yy)==1) EPl1[yy]=Plastic->GetEnergy(yy); + if(Plastic->GetPlasticNumber(yy)==1) EPl1[yy]=Plastic->GetEnergy(yy); else if(Plastic->GetPlasticNumber(yy)==2) EPl2[yy]=Plastic->GetEnergy(yy); } diff --git a/NPLib/Tools/NPReaction.cxx b/NPLib/Tools/NPReaction.cxx index b2cc194c0..1eaec349f 100644 --- a/NPLib/Tools/NPReaction.cxx +++ b/NPLib/Tools/NPReaction.cxx @@ -231,7 +231,7 @@ double Reaction::ReconstructRelativistic(double EnergyLab, double ThetaLab) cons //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... //Return ThetaCM -double Reaction::EnergyLabToThetaCM( double EnergyLab , double ExcitationEnergy = -500) const +double Reaction::EnergyLabToThetaCM( double EnergyLab , double ExcitationEnergy ) const { if(ExcitationEnergy == -500) ExcitationEnergy = fExcitation; diff --git a/NPLib/Tools/NPReaction.h b/NPLib/Tools/NPReaction.h index 3ab556da1..d8d638aea 100644 --- a/NPLib/Tools/NPReaction.h +++ b/NPLib/Tools/NPReaction.h @@ -94,7 +94,7 @@ namespace NPL double ThetaLab) const; // Return ThetaCM double EnergyLabToThetaCM( double EnergyLab , // Energy Measured in lab ref - double ExcitationEnergy ) // Excitation Energy previously calculated. If no argument given, fExcitation is used + double ExcitationEnergy = -500) // Excitation Energy previously calculated. If no argument given, fExcitation is used const; // Other // diff --git a/NPSimulation/src/PhysicsList.cc b/NPSimulation/src/PhysicsList.cc index 5b8679f84..7c7c186d4 100644 --- a/NPSimulation/src/PhysicsList.cc +++ b/NPSimulation/src/PhysicsList.cc @@ -212,7 +212,7 @@ void PhysicsList::ConstructEM() pmanager->AddProcess(new G4MultipleScattering(), -1, 1, 1) ; G4ionIonisation* iI = new G4ionIonisation ; // mod by Nicolas [07/05/09] -// iI->ActivateNuclearStopping(true) ; + iI->ActivateNuclearStopping(true) ; iI->ActivateStoppingData(true) ; pmanager->AddProcess(iI , -1, 2, 2) ; -- GitLab