From d80e2c3fccefb2e3d184654d56c6b9a082d71b40 Mon Sep 17 00:00:00 2001 From: Radim Dvorak <dvorar10@fjfi.cvut.cz> Date: Wed, 29 Jan 2025 15:41:57 +0000 Subject: [PATCH] fsd v25 geometries macro allow arbitrary angles. Reduction of material budget. Provides binaries v25e, v25f, v25h, v25i. (cherry picked from commit 77882bafc2e89f8b1355b72a0de3345f8d3e1e04) --- fsd/fsd_v25e.geo.info | 43 ++++++++++++++++++ fsd/fsd_v25e.geo.root | Bin 0 -> 17479 bytes fsd/fsd_v25h.geo.info | 43 ++++++++++++++++++ fsd/fsd_v25h.geo.root | Bin 0 -> 17480 bytes fsd/fsd_v25i.geo.info | 43 ++++++++++++++++++ fsd/fsd_v25i.geo.root | Bin 0 -> 17480 bytes macro/fsd/fair/create_fsdgeo_firstRealistic.C | 35 +++++++------- 7 files changed, 148 insertions(+), 16 deletions(-) create mode 100644 fsd/fsd_v25e.geo.info create mode 100644 fsd/fsd_v25e.geo.root create mode 100644 fsd/fsd_v25h.geo.info create mode 100644 fsd/fsd_v25h.geo.root create mode 100644 fsd/fsd_v25i.geo.info create mode 100644 fsd/fsd_v25i.geo.root diff --git a/fsd/fsd_v25e.geo.info b/fsd/fsd_v25e.geo.info new file mode 100644 index 00000000..c07058e1 --- /dev/null +++ b/fsd/fsd_v25e.geo.info @@ -0,0 +1,43 @@ +FSD geometry v25e created with create_fsdgeo_firstRealistic.C + +Position for maximum deflection angle, e.g. E_kin_beam = 5A GeV and 100% magnetic field strength + +Number of SMALL modules per quadrant: 5 x 3 +Number of MEDIUM modules per quadrant: 4 x 3 +Number of LARGE modules per quadrant: 2 x 2 +FSD thickness: 5 cm +FSD front plane center coordinates: (31.6478, 0, 1000) cm +FSD rotation around y axis: 0.0329867 rad +Hole in the wall has a radii of 10 cm + +Parameters of module fsdmodule: +Size: 4 cm x 4 cm +Thickness: 5 cm +Channel Depth (ratio of thickness): 0.75 +Channel Width: 0.2 cm +Channel distance from the edge of scintillator: 0.5 cm +Channel radius at corners: 1 cm + +Parameters of module fsdmodule: +Size: 8 cm x 8 cm +Thickness: 5 cm +Channel Depth (ratio of thickness): 0.75 +Channel Width: 0.2 cm +Channel distance from the edge of scintillator: 0.5 cm +Channel radius at corners: 1 cm + +Parameters of module fsdmodule: +Size: 16 cm x 16 cm +Thickness: 5 cm +Channel Depth (ratio of thickness): 0.75 +Channel Width: 0.2 cm +Channel distance from the edge of scintillator: 0.5 cm +Channel radius at corners: 1 cm + + +FSD size is 192 cm x 160 cm x 5 cm +FSD contains 372 modules in total, +small sized modules: 156 +medium sized modules: 144 +large sized modules: 72 +FSD volume center coordinates: (31.7302, 0, 1002.5) cm diff --git a/fsd/fsd_v25e.geo.root b/fsd/fsd_v25e.geo.root new file mode 100644 index 0000000000000000000000000000000000000000..0b0531e6b8481084cebfcdbae2e80544dbcbd2a7 GIT binary patch literal 17479 zcmbWd1ymf(x9B@~a19pRgKKbi4{pKT-JRg>1b24{?oNUeT!IA&8r*pe-<SVA@2>aG zz3a?c(^FktTV|&BuXgQf2U}Yw5a=Qv1Ogd>KqBHG5N;9B20lUHz*`IwI2eOK4AmeI z#upG6D>)p($`tBJ8R|3`23+R}aP|M{6##lsB=6=*yai$bfdOs68VCfXEM{(HjB4U& zr0c@WW=wBtY)kL(;?e&{7YOP_8hBp<f-gRb0q<8WfX|@+lWRD*FCJ`o0P_1F_E!&f zUgXh~#f@!MZLOTGjRhSYjjau=+~53l2c#+UW(gbQu>o=e8_R^x_@v_AJ<*X)iz-Zn z%t+!-%p>xKf;(Y#Y%q;ivRi;G7K`c3Zj~&(MxU^zZG*Tn6LV5gRwDb4e>n2ew5k=2 za)XLkVq>r>jdkT$IM`6XIXLP$tWDQ0$4(YUlU?a!*YnG3$D5-TPt$SE^RAz_<Ii%( zD}9=dr(6lZ0e8kGinE98l#plX3YW$D{-&|>;9DjHw2t{iRyH5LPr8X>MbdhGY3k|~ z!UOb#hU5A%)<*50cslNsn5fI2Z;p@-S*$z7w+I$^+fP+1S}pW0Hm9s1Ebk}>p>2nM zpLF}kv6Qmyka-cB%|Y|CxI`hyT~y-#Cd|fOS!$v(@Ho`^u~;;ftKA3jeP{?KxCy}+ z6vJz!NJZG&$IayQCWO7Cwy&g37=JTgUM;F+EKMAl6<?hhEOt`%UOcC_qo#wDJN^n@ z#v+Q7yuFsmgreWs&1tzR*T3`6##3%?(P#5)OhSYU&P(_HY%GsO%duM74jHZhC)RQ# z+cl)mGmH|}T4ipK)+ucGV9FsOY}Vg%!SMvO<mY%EYb7h0&ig#y8SX^d;^U=a<Kh0! z#>TgFN7?sQN$^6|^u-vF&UK#sqLH8_A0U7AH6WyE+(JhlOtaR2RT{S*KgKama^o5^ z=K!Wr?jq(e$L-x`aV~Ppu$Y?6&1}67b7c4XmGcNz?^f6lw8^r&*{dh;%y^ERx0XK` z5r5;7^30C9ua8)Lf_S>!u|8R%vm|goKF*#wx3_%cwH!t6)kj=othGcieK0IS)M@F% zgED^Hn;y*ntmfr_6I#maxH{W}cO8Ls2)othxm)${>7<^00(Y^=vy#{M_nKiaJFnrq zw^d&;eTQR2hs)g6VeLat^^+AHtCf%X7ZMZd=s5FyEG<4*hN7^{gWnf-b`)0?rGzWG zZSJZk^@S5%*Hx4sC^~+2TxNwZSlXUE-jFt}em#}eE#m8Qp?*jv=YJr(@2j8tViW#7 zXc?MlgOWMl2nU*c81|coo)pTLLCw!92|mac7ls{<GH_k<lr>BVlr@G_iOJcOX^ok4 zWZC`f>98;Xj+_xXTF>j)Rx9X>3>d|_J**GrvoabsQMa8VinfGBI*)opm#2E`g&uP_ zo;E)dbic6g!gAZr@;a_L*r}WZPmKDWO>JE-+pf6oju!je#CjiZ#@%;0TyK?W%qK4- z5G9`0kgs5Cxh(h7fAKhbJm_>R*WT7uzfExM+v1fxnVH~syMJr?SmOFDNOsf#Wp<M4 zm{30VvDZGJHp-pwp&ccwX|STmYK3EtpBd>49}j;7tu8X359T|L{^72v>`_34h-Is; zmwA-C>03Hmn;w_Bg431U2_^i+^jl=`*pG`cBR7n*bt%cJV+!&;S08!!WIob(ba!SH zucT29WaSR=@Mt4JScoGpq?{{jo*UJph!Y`CC{|=RWAN(fj@LJ{#c+QYjB$5M=0qsL zf9aZ8Orxx~J-+$K=ecjFG_N^6Gn-b7&;9jwQ~0)`XJ$2RmkJj7S3UGZ{2b2=8BgK_ zE)Q>|yAH`lu2<G*ecNv@H<~-U7U1SD)Kh}U)=KBIzjkJQqA3sz6&E(P7s`%~64H-) z59-IhD_o#+N;Mr+5-SHWg=f)&<t3_u8Q^(GLWVkMLWRa%Ai$!Ke}zLa+&3f_P{%6r zYx??gtjN~Cz=YL119Dj26`yEuhJ{?fr^Lw*TVWy|G!>>|JdCY0kqn9pVfFhKS-LCX zH!Mso5Lp_N0(2?+VJl@5fsr7Wf|!bYQX$_Q1cZ3e->w=EP_q=dSk&IiRc#Z$B^vy} zkRlSwQvlxo<xT`Pb^<?>*ihyX4i?2-f?O!FbafywgHwb+5zIh_S00V{EEf7>5Yy!d z8qrxS;OcKBGQd@y?U3;k%yZy^vn2RGW4+5D|Iz8u`ETQn6*gAhBoog-=18n(hdbqV zcn|Q4PxwXWw^*R11c;_VCh{o{L?G@L*8_(rim!tiOK!FLE{_I%j6SnJ)XbIHSbtZw z+Ad&jvRhGx40={_R<+ohr0u<1xd^V{<9p}L`Q?+&00T?{V$tPw<A`F3T{Ro5&bI70 zGeFgb)1ARDTlv-3JYM;;xjpRoZMR%?^T_QMqlcJdlpB<;jlCi46Bf~`tLzLnDs}se z?|J3}YMv@{3|XP?SfM+8?Vsx^Y>T^RLcAKRZQVHA@RFP%`OmW~-F$u*&V9+uvWAED z^_Sdextqk9G1v<?l&$pjj798XdjIf1@L)92h33Qfi7_Xn!=ekECM2|(*Q%rVcU4*k zBC=VK<dHtQyS*Xv4Uw^kX^*^3-n_msUEU2|OkXl;2i(v8otX^_mAo4=!qD+#6zvKr zjO<?!Q=Did-$>biAR$Xx%Ea-6%b)uXXSU={@TAJR7YJiy??OzqxD{}QQN=igZlZe= zd3+5FJ(CEb;+uiV3wjatrG)wr)QkB&|3;ZGbOoz<vX~4QiYtGK7<||bF(xVEULf>i zgk#y8>a{K0&(TH_hMKwWi})@GJi6E~h6~5EtuKi0d^;}NP&Fjf-i4CDT#kfaQ%WSV zewA`NLxk_cOS6rfiG*-J5_jv69F2sldLM;p{5TFJ4zPu7+U5(sGfH|EhlILyI4@x! z!gqEXD4XlI^?p5{esJ;FE>}M6bzB!FGAkwOtps4@iWxR;g0TF>9W|kC$wWX)pIZ;a zHy(EpzUeSP&o4BRSA|yd?&H}1n5j<e5B1s>M@|<48&c{QkF1|aK$+Sl#U~kD+?5U; z7N-_jjjY$+MY-m2jPfb~RC^|oiSgxEzPJTXCXyu*2G|9Dw6i$5NWf|KkyP>Rd`Xdg z2o={`+<Cr&Es7ulX}bZIwqsbV=zPy#!FNM02p7)m)bFZvQp%5iUoN^qHTxLyeCn`K zUktEjS?$NFw*9sSmmS5RV_Vf4F>mQloCAtluBGs|+;RKkp+CpN*;ml`D*#n<adYxr zG1(8_mALe3V)ysLOs8h$AF8d)5wpEGEWsW{zB|+rthb>0y1Q}D4cF!jbc6^(Halp4 zyYtVIYGzV(`Nh=+OV5pYqw&Br50UKl%VRSCl5{C|o(1o{iqpoOnV!T3WxC9w=;C}I zGBiyL+`Qp06qKfrN>{bjpaY&J&sLtj=!rkrc?f$XhqF75=M)tckQmvU5R<sDD5ArU z2h&r6=EhCGZ{G}<-hj46z_oX9n8Rb?K6(Wp3#6>u??C3gX_<Jxz)X|#l_w;b^;^Gu zJ#EOG39pH^6iER%M)sO;;GCqNl=`ANJNgPzV&lia&3aFeLq14*Ck)^be$LG{f!~K0 zMQ6f=5E$9Z5R+oCv4qt4IfIU#!iu`oV0juv@xqa%nVJqrIPMZqp-h+W$WW>EjSamE zR&=H`2ETaJLx!59z5t6|%Fis>9CXQ8)F}flQZMr|_!k-xa$dYJZ)9n8ENF_ByY5Yk z^965hX&?9i4V1uFD0NP5S5aU_$vEQuX!RR(L2-!NI%-c}o`etcBbuN_7ew{J8T7la zD)&2&fOhoj?85-UNaq4|;IswPpM$2nbr*_s1P;IjYD*zL@?XJg?|^CU5aqaB_&vV% zO=7-igb<PBW#5eR5jdB`2ZXJYKnjZT&`U~VzQS=3ztu6yxlpVVP!v=TPzI`yL?F;e zq%86C)aU%o+JLau3<<5hM1n;nzuI4*!ccYbPdf{px>Nz8+zPnI=Rg0`5L@6T7fD>? z5{ay>PDW8MIY43hMi7spuYKcm!m@H=@{Qg=VsNn!jfDTeIkfrf3@Olh4Q6q2k!{i> zg#*0|@2J4g*8sWKg#9-hzx><d`AG-}Kw^UoUSi_j)^*Zsd7q}77%ptnllXswU6+y= z009vGojO9>=i;V6npJ+VON4q^y~OZOD&S&iq-Lu9l@=Lb?j<(aC4txk+^zT377BFk z5%;-xun{Co!ZCdg{_95NzrlX^Gx1+5YA?>dNOzGFpOJB;j#U;vS-&_-F7_5s`4a94 zVK&vR2jaFv_ZX^aKR*22oNDK&tIju<bXEs*h!`Q5?3wPXjw`S1{#L14@Kz|i5<k{v zz;z^S+jb=vDhde&xdcVA;>hpf1MIW~Ol}Z{CpS@tc7R9|=WRNE469UaGsN_XL@ZUx z!)mm6I7T`*(JF;t-X}$&@LBv=AGH9HmV!WxHb6%BrTMwPGSRyjGR$G*XK}b}Zp|P@ zT=(X}8nILxxkRG0l*AwLW4{LRSE2redly)=U-4*Yh9@0pWYxo%G}UmwfD@}QkK>uD zY;#SC$PQoKfNdBV_^{e*LW$VDWV^0`BKpA~=xE?vDs~u97}@bG5g?*7Ex0*1`cw-U z)?N_Uzu=kuZ{uDYjqEIb;=L7I5S>f(fzjck4J?|wBrKrMlT2cpoSbZM?Iell7SE#g zV`XRx7mnvoV!^YiDItTyN4}Sl=Ihv=UAlk)#g{IymoAtvEnZLHE4FXLpri}wN4S?p z)Lz4o|ErBKo+vIQF&5j?SN*RPDr~&`0Vd5{cu5^F67Yx+pAN9=OCm_Z@E#G*@Pt1J zkf9(wD&6UtAFHg{-wXeRZ}arb0z2hVv0{O$h?=*s>c3|LJ4GRzS0vitzc)cHf{{-i z9VtcA_B<G)>1uT>aDjzJ5*~-7KnpZefaa7t27c4m0R!OgMDbP@5@DJq4*!ZD3d~hW zH4f!p6`*!B8^Rm)3!}6u&Ck6RHqMI&rn5;zY8!+jpQO-8UqYp}9%4#YI+R`!pKWl_ zmy|fKWoHDB&?k|W*&?!@pB^|b66&XmY*!Ofc2F>{Z7mdtOlBPKlot|l9E-`?qT#ZT z#x{Mf3=G{#dzpIVB|^=rUIRuK7$F-%^;CfGbs-6hI`);6^z6EYAuJPr4%|yhr~yQi zp0*q~{cZeF1cRF!<q`i1-h1a3yz`LA!UkPRE0W5^h3L-n5bJJK_{XpR+Q<4bkoe5i z;YW>Zhu<8U@STFCCJ9a`-)DT`Z1>P)&hQt#n55yE_@M;dK0tZ_5X2Dz;VV5_{=%Yv zi5MWAH#$fBa-EOnKTU$fSH294BZ*`K)Otb9_)~jeDD%CO6$YG^RTa~?kA#p&Y9%n& z-V>7jL#6s>g_kZt<()aHrvzZI*2{E{MJA#?kFD8p<XR7Z<Y1bD#WcRgKzvrO^D-yD z!0HuI&l&)I&a__c+7G~Q{Ve(~j2=b*aw_F7Ri3H;#n4k-AdnWvc=6q?1^}2y@)G$g zWc!yYz9}fT4+^-zzxhl|d{(PKVvzCp_C=%V58#g~{!5|qzvG`@aC#a?$lIm$14uKv zDE^wVeGv)@)uOp|{_7C?Ux#9(LV3dEOj2-6`CrHP178db7-hz+;IAeAO*VqO@9_!< z7n+*7(p&DANsT-I+pmQsuzh>j_%427`1f-z0v<vE4VvsywM`Jy?Bo(pHkg%E1Di+^ z@=`Axs<fsz334noA21TQO9oIQ3Hx5vRlq~)KPcj*`v*|J@ME*$h!p5=J<)34V>C|U zr~i_&I`O-(e#MUBLmW4e_kURAOZEz=$c`==!L!}Qw60VPAiCjq03baBgh(w40)XX- zM#irvC9x7eVoEOI)7=2XIX@*3#{1jT27_xqQ~r(pXMw-5FZ>rb=l=^tJB0%P*vMX@ zpP1?DiU04A7;*`|kUd{Ee}J1%FaYA@$t3u=Uae)nz>W)Ieq`a^^hMDNGqa`t1)y#$ zVD<rIphRU=dvOdf4*svkh}W@-mY$Rni9tT*(a)EN0Q^e7cAyU>qK_0i%f7S|l$q@a zMB!grA6^70rAFgNjKTrHvErD^tF&sCtuz?KA9u(WNrlcC8Y|JX&?Tw5qy%WeDj$ch z;)gEPYVEO1cr9;ZCx#k0!{p_Y{R<~8TCv&|h$&WB)Bzw<$8;s~2n!4~a3T;LvbL_p z9#}&X5tk0765|r7ZG@PTCKn;Eu(H4clsBqJ8d$>Q?Jv+>q?%R%cPvSGA%19d9Lnt! zFeA@NQLvDS_%tMjZX&M`iw#KuR4Lt7VBfO3c4spsnEEvHGHWZL<)f*7U?Zey3L$KQ z&oFPpF`bA<`S+sMD}R9d5)Xh6+so`aBn{lOiUdsKe1gEWqhW(>)t1CIL(>^ezinf1 zk~?@PHY+8Dio^Duu(UZT=)9H%ri;|AcCqK@BizPPtVvDg`>+=(d2qTQ01Wx<fWg)T zTSzDNg0IZbM<D<{dvej;n?C+?q6ROW9pK!+S_6g75%Rwb4167k$q||+&DR7$-tZ5V z%m`koMr~-E|BAy*+7~c|DI6XOIEBA*<;60!WGGd=Gmv^@Qv;Xz@ZSN8p>I|ul}5|} zZT}$?@JARo0wlJY=j`X5I4lN0XOQ6?OeXf4Mb%Fv0?}SMtH>K!ulxQ5Nv76{$Xiyf zOh0$Nzs?6hKY1AH?Xc}55hl;f72`M`@I$^PN#4;f>=0De2mYKY$8?66e-fREr2c>N zt64OXit+}{vO3^L7Xn}~Y}yP%`(d6s0^5|~Z$jk?Bu`Hz6|!Zt1{p_4^E>I5+oYM; zOygBAE2)=M?5GoU9OvKVm7HYc%R*;VbEF?pWE`8gW7?ur3TErj=B1x6n@3NOkjYn# zS%c(XWc{x-+P<OjtF}j6$w;2CZ?9rjuG|u@t1ANl;9ui`bruEsTXi~>{m??$>t)Jg zZ3k#`s56?%<T9uJIV6dWc&gib7N<Edm9%7OXtGdbo}5n}H_Ney78Mn53mp!R{2K6n zl2f9sj&0gnJNmf^h&q(+#{eT!2sq<j0rpQKRoU<3NTNXho_OFllEkhgz_!F#16P=W z^M7)s+IEQecM%y7;w{@JlKuJq;F?%cCkty}=5MdyGk4-Z|D;avf4Rkw2%v7Ch)J3P z-F2N4Ci?j4l09K_N&u`fVq0-WHx;4f!s8H&Yb-(9r)a+dqKwhkn7KvbZ+e9bi?;e$ z&T1Fgp*^^R$MFwhq{$Tifs_u-zqHbMq1_9%bQplfZANd=$AthLgToHN*Cfm!S%%`> zKdI18{dFN4S++2jbOQAYP!b?V#Km6{xpyy#SE0Y*VTZ<FN{amjmhlU)2nAjyzE?^z z{3WU`@xOLE75`(%HlFCVVNN4WSs$?K{hC%4+*p1)o$zOWMI*$}LfcE~mKNDiOs!3) zVy?-AbKL%ugxDdw(n(cXh0D+tCJ#utU4!RHfMhnb5eOl%V|Xm9)VHcGrm0#As{4^+ z_@U4>x&Yt0oZi2Xsu^O+1>o7R_Se2-4gH7#sw)?wC4V4p<FQLo*ea6{MARfhrPMV{ zzbXf5(Qr}`pWWuyp_TMjCanN5m0$Q0!I4WqIJkC*TmAqD=5GdM4vu39`T)IiTjnlO z$RVa`{+)?EMTC@lzj2ipQfo)n1IrfH`G-Bu&ffe%j(gm{<eGowABw#t`OwSK-WafZ z+ZI343}jtP+t}LABs?FgaIshIRG!)X*4bCTxNk@$2(_${eas5%nqY)|fQf0WIE{Hf zXyekQ>B=^;@csg~%^R%E+zm(7r)UB@;V@wLclBal=2vM%ZJE2bn#jsy{$n+d&pepV z{FNnLj|d#ewtY5t(*2Avz@7-Om-2FSeBD`j-J1b+Q-Gb3|KGhC;3BX!^gr5A0YtA} zt#ty{PCy(l`!Mt^FV>J=_F;r=tqshT9rSG+5nuLT41v?vGhhI)BlY4X#((khWe@1T z@e)vtdjzO&ela2<?#h-3dVW4+O_L4*`vg-=Ih7~$d{I2%gXaxz<}7?CmsAxE<-fY^ zqh?M7Z{W9~S0V;I+KkWja3sy@ci^bpf0+x3jfyC1s9cP9Bqm*Lm6e^A_QCrWnm>AL z6t9}l?b_@vV)JGpNa+TgLopF7(%s{$JM<$sPO}ciB%##jwaqggR~oY-HG7M=T@N_{ z5kroX!NqfI10irfQKv%DAj9T~SpNfef5V~;gw`1yW#>jq0aX5vgp7<~y)1Nj(#=}E zRo)|~dTKI))BXBBtDkplH#*xAq@h793OZvb(q{at74YR{Mmhd8=WeuMA|j;6j%+8I zeW@JoO<7XHS@i|W{>h>lyfH^nntNmO8R7f0BoSt{qJ!YaNI1=%WOm4f8#eYs9aJl7 zXC&x{rK+Qhb@^guI@h1!P>=dk35n#NY&C^s?%;`cykTik&ZGC?$G@Wce1(g<c))#> zzg%!-;=A^E*CC*EXFw@CC5ua&-y<7QA_vh-4gMs<%z<j1ZWH6>d!KKf(o0qIJx(K_ z@yo9={9v<=-)>L3UC);L42BzoAeW%(#za3a$sH$oaLnPto_GRCY)s6_prCE3-9Q7C zEb-6jbI%cH9(08wA8<f!l7Z`X%&;Vi-#pY1aV_tvaOnI#T7S6AmBn70M`M?V&awi# zGQwPdvK?Tuj2ehXXXGgUP20;8h_XnlPBEL-%MKzo-M>9Zj&C`E^Bsl6y{R|7+r@j% z)3v7MGSuIrO)iu5_<hdb9-fOBu{V|+-!{*+7+6K(N{qYG5UN9Ary{=_1_?s7wuWeS zUHgNTNoAHM4L(`Hq0f2xJ8r_ncJ|M*8jSdi?jQ+!=zBzHy42})mDK4~SdsOl011ys z!Hb$n2sHMa+VRqwDLv|r0BCmmE_x*mjB%b!1eq{fE@s3hVXgK_r^oIF+2n|D7rY%} zACDgx!k*Zfa_ggbOZDVG2WF_%MShyS$LJ)XYNw}KFvke~l4hvjB!;1~&fuP_c2=hI z!|O*Nfv5D1(GnYxa!k9Z?SP4v!4m2xu~Jk{IP(F&)7ms-0<w|<zK7|MHNH%b0@mcJ z!86ew2hAEU#;9MXtoeJh8JS*^J#sPFzCMn7KSUl<PngL10+hdE@O}d~S0%#NfZA_R zt6Uqwni2D=MQYA|I30fH7msIdSC|aWPPsk~+YB?k4()CrpSAL5d&D`12tz|sU|kp? zDD@#I@;*J@Z7Ty=CFpu+!kO8wgSy+73t}bH$c9(X5xYNF@C!^i3W=AWLI4kSV|*wJ z;?~;v!nfQQNM&)_!1r{x8+Hb12%jjST$UrzMSgbF=*?%+5$pbmR6U477EI#f;j8F~ z^bJDb85_AtwKE&$Jy32R+?2Z#Wgu;2%2jkL!pmu_l`A<$IGDjzm4MDl#PubU9U{JN zP%G!?0aLT7-K?Wm%i$&d5u^g1IiHpHW!&*2EoDh5D$au`>kT>csj}|aAn|ox*P(K5 zNO^b4gIHy60z<ASJ)X<^gYZ8~*ejC+6le$*0bo3IPj(3Fcn*{UEj)Q+c6*Rn+EKhn z?Fih>WmsagMoWr&RR^u_XhK3D9)@EMe?K}zZGfwK=Z0})C`_4z)s8{b<rQZjHL>Wg zgnF*o%loxY9e3Mjc9$kChHPmF9e=V@CwVF;<k=?-UEY5alw?j}&80D`E9!`}th0Z? zS$prqw#g9hZ^F`B9yx&HFk7;`nAR8FrsA1Jp5G;bRf3g9C)i;p`>~4#^-WIvg5~re z_Amy;WEN>dt=vkgx=aKO?)AlzAuA51bz`g56RPE1kE%W7o3#GfcQ<hpSulr@gmzw8 z=z~Ee8qOGzvHS+>vLS!k)xt0_A`c#z_z_Ub?P>~3y1HO47ro=IWP?R7NE=K?Zm;Ej z!>^ux&M);%Md0COh(EVWzT+wS=xe&nc{LG`;TH(?Ucmbdj5SZbG0`Ct1=qO%s{c>? zRe6OLp@3!4WTzZlE-KRGLrZertp*N*hQ(W=Zr7SsBUg_(u}=*$F}|DoN_&C!K}xjO z_uNBK@4I#sBi?#)5Ip!_rAY;azfsuP!sckS-;3qtyE5&%j@Y-|T<EC@^%3+;dAnIG z(VW#^g77X{*~AGux&gD5@B`MbIJ%ZV<n~aziMyIpf1D3}!~1;$-z3BZ7X7;{VU1xr zQ<6-IteD-W&p~+*y5gRbPU}Z<X-HtkCXZ>g>=V+vK{jfKbOpwZ4tk#4h9?s60$d95 zBJ5Hi`fr=R{alXKC+>|WqHBa&jz!_oD5=w`@mD333xhXn250SNqN$m#G*3m!Txh8` z6xZ27u(c%xM$@$=IPX3{Nd_@2fAi@m&4=msKl~N`7G?YUjiZq|-XBbRZz|p=A`IhF z&7U(m{`b;yAHN<Py29A+gh_NJ!Jd&1C58X#k_aAvmkoZD!MzIQF4X7sk<7K{4C_as z_b~48nZceVS4Rlkj9ew<2(iE730F74>!}XTGSbr0|I<06cy+d%H+zrXF0dG9)FSZn z0dbc{3|k;XZOUKvSJ`Q>dyhp&Rq3IGJH7yTMNe(2c3F`l?asmtEJ59QD2vgb>)KMI zg;hn*7)nQ>80%x4=T`Jc$4T!dFOij~2`d^-@MLGqr9VbzHF#jjZ|2yTY2`v09b<~Y zZ0si%^qm43>=i@CYBeX49yN_T*-T2bElqPB!FKnv%-`?2D}G54DYA;+@?ugTZVl=7 zq6xWXm1{5(*iuJ+$Io}nptR55)fB{7B|50n5W<MPEWoI$6xMw+j>sy!S2WO~Ivg!m z&zn{>@#dYn2Fvk*T}IuZ5ty!X<WrxoJ<<nG16W4dyzg`avAdz!MS;ZE#b8Lyl%_7D zU*KJ4c{QCa&e49<K2?+oe}3Bq_NPDA9W;9UrQ<7T$7R$SuR9`51Z$+*OL=fhFHJ=F zSb2TS^J}Ggd|Ia$S>C;!Ls&pzg8SI@Xs$VABUb0MqUx_G8&z@>4s$--z4crOt+#DT zv5?`Hoo_IASXal<8H-rTMJ?6IWhkioNZP2BKKe*ASBv{fim91b1G{2nV!PA*s%(64 z81FPOr;#<k*?8$LINJ7lqIhy2Yy~D3<xxD|Aj@a><H2E0%}>JM<u!3}l6k}s#zLa1 z=#r>j{<ct};?#Uw%1718SKVu%@mM}Ni;U8G7T;tlySXm?+?K8=UU3iA_`&z5;`tp- zs#-jeiRTw^IS+)aS?;}_`a6eejn!kKw=>KLY?Zw+p0d(~V`-P3KS}geH*9Ac;CvT3 z>He61IXTM0{McLq<r5F!fb?t{Q#4?w89fPA`%n6wc#;Elv!jb1L$KC?mt!ol0j~kX zHPSe@eci`%UdM(^0d4U?cq4TY@*dt11;?MLI`NX>*0$P(?qkEd&Z%A@(0h%ILyJA( z*vP%`HJe=xBCQe6R0gpH78<Q8TtBp;as7o&Q*7n*oy_$AXzhXjp?{Jq#tW^l54A9} zLw8k6yn;#{Rv69NCBmwOVZ!wkwcgUWP*prp;=i*cAmTG1$FDO|AVyX4P>3U5Rdr>q z@+s{7b3=o`8LG6+LC|?VPR=O(yZC;5y$akmDQT=CU3{fQRJ1+>)oloc{$L#eb>kBr z=3p1tw76Xr4mCpQM}N{|{~w>GX~7k#j;lge7M6yiR=Ku)lt`~tnEif-ONg3`Fdy=r zxShi*+=#gLLdm^#@fu+|@YfYrjhQ(xDTRfKl{%-GI+ey9n}AO%^KhW0YnR|!j_bna z`++b0bpbv#dOGq>U8b^UiGR-Qu`}l~?V~My)SULQXfcmQn!^xwi=$jl@(STtNTG%} zE{qA2g__TYu&nG^)ynswApve7y&2Qa`ymMd$`!c76%2EavGdIBvD0<%@tPd-%3qW< z;G#mDl3XQJ+?NT!Oy~_GCA4s`lX8o-40eP0mdtluz&H#~`A%Ie=hzn4q{Oo78`Bl% z4|gD0{%q%P`-d-_;fHbVPXAy&1?4&}&f7?)BNzNyu^Vy;O$rx0CA$pMy3Ri`a*Uif zATm1{92ECg<e7TQT+u9$Y6XTG+)Rm&p_#v<)8FUylcOR$IEh}B{Fh@e+jdqHjcLYG zRd*GU;}FxSm}+rVaKcrUWRO5lYKf>{y|~3hyQ|*k*tyIp7uO|?z{=ja((c9J0>=lO z`Xgbu))9M?uUwm1@EMGdu^X*VlZ1h9W$XPAXmpyL%)UWcfgG38dhdJiKpJXFNJXEU zi?g_r;P#Cqx+5B<WGP(}It0TFC1P@%y!kP`n@c0<)u3C<zCF$Q`;Mt`63lrL#(-?t z4Q&g*y7u$-u+5}A7WS2*sQ&69FUNa+Ac(!m$a+9z_M47(%aAXK!*=;PB)OHNEmu6W zPuh$ro+oZ%;6;(pqO#1gVBKuA?jd^m&1j`baL?e|iEYjopN7qx#+X3QVA3Jjdm_6l zo<^B$!Le16)~P&kXnT3AQ<hFwEUCq>j255SSoGqg0+r|T5POD{Bjz<NoHs7jwaqLP zZq}#xhn%1%MLkzHsjImIpij`Nx$~3m@o+Gw>U&3AR8Y)tk5^pmNi@q&XIm{Qv9MQV zlm_jvkGBa_n)D_VEECJ4mrUcadgBK>#exUofwfg>gO8s7MvxDr0)qH*>}GWVXk<#= z)5aZ3K81A~9{kpdy*0`jYj=&H4Y}M$g~mrRSnXguG8@;P@}k|PHkk2oQ$xX*iCp4J zi;Lvh+22K&;@m1Na(Sf|tqk;`422j=-{514>>~ZZWCD5A6JomA7J)#?&o?~KMCC%2 zngST#?O`RST__jP#>Vx&us0|_4ovgjA(3K3TAbAE_hx?AiQvPY8oM=~``Cx~MJA;2 znn+hY5A&T0cU;9aN>#Ve;jM{u#dubA#H_h^%*WWwLvG{c#HAgP6TVtUJTgC|!woj@ zGh9#A4vbSlR=w{mF+iSwa);~ZHY@9in4?8crL?&|rrSR6(nTHPK(coipr*rk<k68f z<LIhvGXE7`)HScct3-RO!maAiuQ%D2-gOadJuURXM8<P*Yg<jHEf(H^1!j9u3*4O? zX1m3QxJ{F0M$nnUIEd%-;o5Sskjn?Gvy4Rab44r@g$sUbO#X}y%Ay#dDJPbnC6hN; zK)-{dF(5y@)vSNlu&4!PUHFdGW#VxhQS*JGjH#J#+2{E2Gcz?QuJ0vcE7V`cY#2?) zwBQ^I-;K5)_8m6bjjOkVxjNhms7=$4d|2Q+lqFCfv#BwOBiN#AN=L9Pd{@!(q@TO~ z{w&q!nJmF=*SIod)AC1qc@=-?jI-4h-Z!(8q3y=4ZM*nP+@CBVNXokNe<m0b^#9}) zg<6W?Kx<=MZJ?sug>G2b3TS-gZON7>#(p;_r;BuBYnI^yWcH<*==|xblVGkzPJ3wg zwcn|P41YxosC-)in@D8He-f(6yBUhoXDhF<aOtr~9((faFKc0{*uhYx(XL5bhOAI+ z)w@jIS;q?eT-M2)yrVajiL2{_!2F8^J_&YcwPFGB2U6@ei4z;R^13N==#0^hd|0dg zg?FNr?@O2n&vQR4vs<jXrAAoL1mLo;5L!5T=g&pu+FAs6QXm|;;e;Q-V|!TKPR$~y zj2|T1vta4S_RxGE%W|_uqH)~#O5$d$v4~zhwy~AMMyQL$r-SaupQSoVe2o#p!6M!x z1=EV-h}2s8F@jBbjX#4a$>xg2P>y|UG-A189onTTYWYyi{-7w$-98<CrF;6Teflz$ z&8ky<dfS6rT5;cTL5+EEM5dM9azq#N8o~v&L`>1*A_IHd(&CEx7@Vq}DMJM->HT2Y zJ3>c8HcEB>xQ?1>+QWC{mz*V6;>Ku?im@~kc#zNCon$|zqT|>-N5b)kdBPGLd3a;E z<JH!_!@G2Ap@!onxc*7h(7f!oNUjCH_s}1JcK)nhInk4)KE|Y~r^`^WkYR&Tu7W3` zT54BDzF5}#hhVYH+~+yKm2X}2N$qJ%jUVO|%Bc-vHMWu(fkeQRT?$&)Zg{ozNNd** zEPSR8?|p)Q@%K$;W@#?v^QDAkh+S0qAD@J-6zZa{6`m0Dx|S6ay~L>&t9JFF>eJ7C zsN83l8H#=Hf3E$?w1zg@%O{<2=uAkiyQ&SP2Ilf>EHjYW<xeKQFHDI%e6~hwJw04= zsu2XXD!KF_Yx@}wXeWl?N~Ky$4DKCCEqKL`E(Q}3?S?<ss{M~G4^DKF?q#}ieIYZ! z#|q1h)=Nn~N$FszGP)ZmW%uvs9X4&qWc$?rpyfnzQCXU0j&E{f2t3sV*@d%iu3b_k zdv(Fp7W*du#xCBt-G@1PzBHKn;A8CcjrV7X6jVIyx%f0LTf>#F$gQS~ugUZ~)Tfgn z`}T&Y=kg8qfI7#3=jZ-E3R!jIfiY=u_@mqkO}eLLkX@CYUKi6_;@QnSWkVc%xYrx< z#_ikAUKcYaV=YKL(R|+ET1uz550X}$S`O6N7CCazS=<}BeLGC&yW+BW&kz)bGMzw0 z&E~bv#{K%7^lLSc9_hcYvjLA@|Ie#|ysy^a{);uQ|K1v4i@xQhmIu84N2Npo5U3Rd z1{{fBkM98Q0i7Ta1O*86zf~Ot_*+lRyXjtSy$lhEQF{DvUXm;T0=Yf|ua}agqOFs@ zlew)82<S9;)rj$WMEZ|jvKLYQllurlC(`#X?UxsS&gK$~Kq1pVXTbUC-_}2$UH|*_ z|HuPR;RA(G$^O?usKATwFVE(`6hdXp9i0%AoE(hxt&JTdZA@&@L~Nab=h=y@%xx@< zjerkFC*X;8c??j$B=AK0uTPd%c=OJD+yyt2_xc)*##)9HjpU2rhHbJ+?;U~@x6`6U zsL|*Wl#MK1T!Zo4Tk6wN3!5V?Y?m?UaUtyS@}m<P?G-OjFE$GHNCgzSCu17S{W>3t z_p|fp!mo5cVwgqbnB8SMo{Wyux{he1l|{OQrY>r%c}}$KdQY^}V@7oer<VwDAat%b zFX#0(`Ggw&_QK$}3P6uk)6PCUnm1b@X+<0>6eg^KP(n9TDC)ZJE9MXGglcG|_oj?c z6Z_n4W7kFVBL!+|$c|B&7s6x|3IFNJ{deCb&hMd1|Ce&DoW61~#*vV8wNX}jMg}ub zu6@SqOXbvAeIAdyz<M;y-<98hUzbpTX586V*~Pvn9`3#j@O#q3y`V_hZIKpa>8MkR z;+b~@C&3+XhZ(OLB0@(40n0sEYP>#gpyo6Ii^qa!+MRip7hbOd5q{(E;pZ=jk!69V z-M@<FawHoJ(`En{Pl0rPAjj$pajkvyi5c`I?93rpUFg14Ojccm7y|KA>JTX`5ecqK zcBu!D)!dA5Unu`R%OGF<&~P}kc)6fi9lZ8w_I<hzTU{t-n$hUiF;FKy2fx<CJ*IFS zjp*$VWnk{fl3GL{<t_`^VdJMlpF+yG1v}e*&U@cqz>jLMme~!?)FhwjT9tK5^ER|# zS8czGdZ;t%Ps@Dz9#q=&?y`Bf&9ojYqK&!;+wVRV)YGB(o-niT%9h)nT7<G~*kb<L z;iDuC1GFlKj3Nhbhw7hiR&`~3ao~LkYiPV$N+Gx5nTxs6YXs;KlP`$6YZznbX~|NS zDKW~=@NidZp41ZOEvd=tMnK&WB#<HTk|=SCbIawE>NBomW(0v2wWToDy{DbY&791q zM({O5_1)v1-=j4u889;adN|*jg6$9o4x98am#V()<a@&@J(10xH0rH&V@NUCoxQ5W zzj1^0Am!bIgD7m8t8P9yHQ{pt7yN-XZ$4e>SjXHHy)}F#$jq*f#KXT^g*S5u4J6Xg z0<~VnIizcsm!nS)Cs2PFSlC@I%9q~pcBon#2V2q?Ct?Lu#EQadECr!vB;`tm;;#lr zvd;uP*(Zc7vC=KRT-+=6x%2`JccuF=O~zP|<6MN9k-T&N6ypZj`j61eCB)1J4f=0O z_FrIb_w|(!j)<MO{8=2$CX&6d<jeI?y(>#U-o)pemgVRbD_|R0M2k(CD=7PS?FH#f zmAffy>7#ob(rqS~tNDMVJHu=@PPYy-7q`P1@Ytp*cni@hUdHDZn-hO`T2JV8KJtUK zu-1t!c>*?BDd5+mIc-Q0wrKV5!*yvm7mxm%^MlXDgIg0`F~m5S1!bjs3TP56;QOa= z67n8N->kiutoyMn{q8rek-DI%M&h8=tf=kk!H+Ba+87XR{3x$(MyzA9KS96=#r=YW z<m><VWhDE1bC+Zo`}o~Q6*u39HZkZbHOQM76Tdd}(aL)(q|XbCZ_f4J3I!rFO2vu~ zsG`{y5wElnq;p=x`Ted5Vl0&G!;_PV^FIjewn?=Q*lXL}JxZ&8zfF@hWl{NuxMVtH zDCBF(9|mBy*TCUNM?!2h&8Sn;UFhlDgX;f4p$s>X2=N1R4t6$E-&?Z}ZlpUwaSOS4 zpJR`2U*)@8&CevQisp*%#%A(wob{6gc4v#Jiq2^<zQdMTh(JH*L3aTKUJIH<jy;{z zlnl5FJE)sH-RflbGtq|OpguWcZ6+P~&j(9=A_0ObjDE!fgwcE%pN|7MQbHPVkTSdL zX=;c3Q}kif{`7_h$!Ni?I?K8m{yMeuBfO2%Zc(dZ`qW)E#*PxgLcK)G)aWjde3>}2 z;X)Dl&Sir07&exD=`BYUI>#*9y&aGAa3(yGbUz_$66GVlJ-2mk4f9xmEM@uKce8Q* zW$7b6N|)>^u?J4(Nh?jaq*{*8Rk+SHT8<0$O$HOSMi-d-;Fm)A?@O?%;C?qzf`cie zdCD$fL!|a{EI9^#Sc>LGk1)5_*viuJj3n#J8#756<2R7yU!aW!2_`gqKs;~PAAoUh z!Qf@@2HNvkDC0qqZpUD|_385*B;~t9ZQ=TkB1803<C#(aJ~<HbcF`^HPS8j7COpz` zcxuVK6JmepoBGiu0-crhcEp`g+M7VBq$~dyo~NmxoIq~?Mh?L(2zQDyO{zuiP1Mpf z_@<lWy~<pr=VaRZhNes*vp@(|-=b18AF)B{WKF0yr!I}@hjJ4Op%KY@Q!)JZaG8u% zcHzH#bnly)eyM3&)EXcUeEvDFsZ{zAmH68yJcywWNceR_;?H~Rz7&tLc-+{SQ}UK4 zbuteau=5GnN&0<p-y)HE;#7sNg?9pY>c~@0{jpGwaR@Q_82iMqIeDAaarPOhmao-h zbFF^5)=iHl%2LRC5C}A!l$#-^U3DzV%?Rd<R_L-_?YQr^PgaI^1T~DRB3aoQAL5&X zNl;3~?}@>7Luq_EnjgnMaPQh)L{+7lZ#UD~6cIh2&3_kG_TwYoTr;lCk^=dl4t$ox zY$Fu;C%6YLB4W&>L>VqjHp9LmqS8}i>5*zSrwb5AY{#cs@Bl9y-(-Xc(lwt?gc33Q zAW@BqA1+O*>rZ*;%7?i3t6&-KB=Jowf}PvVicC+tBQ<WhBQ@DleFQF)Y!%`R)U5R6 zaoD}0)dlmmxcO<25&j=DYCYfVqKsFG^C#nk{LvMbuPdyqar<dO?Rkqh_1ha(C+E<v z?o&?nbDbj+=|yhe_r2~S2T`}acD=7UGs=yH?N_8=s=EuL^s3XGn2t_|kjRTZ{Wx_u z_Y2C*hjOJI;t8C&LzveMzGJY19DAL@d?d$%rj1E>VBfjmkKCBu)uTbukO$N_n3N)4 zT83so>@CGdj9jBB9Ho6b>Xo}vK5TUz)GR-$8tW!o{LY=^*^|zB-j7O?D18{Q@BC7N zvMQ<9nbhZrOPR);V_y*wMfmPKrR2u;=b{MqJc%Llk}!589T9eUsfeD}E_52PQ{Uz% z8YA)IMFP1Q(=Uh>6*$+8oN^j2b$-<xVoy6iIbqWAYJvU4slp<IVU~XRgx*-s8PNSo z!pw5fn0wnMK`F;9C+u%^U)9-M|2o|eH-ljwkV)g%$|q3$vdKL1I;U(NypbLs^_D%6 zY@VPi^XX3eZxdQcSzSIfDK3r{k)G+1_8rfiNmlS>cwGzX@-t1!ax|6<RzgIl(@Foz z@Kt&V@Bt59Ncj3Ro833qGy6$-Yz1f(Tfd0ov=iS`b&;pe<nuG+mp7;lF50N(JDMBR z;}SxJJJZ4MuNKw4`}y{2CFzu)f9=~H2Yq2_i6H;YP68!UNPy_Y+Y7Ley(kx?u(d5h ziGa{D2DEYy<HGlec$MMGCs^9qO8T`GPbz!AVcuxFk1T(suwTD@UK1@PoqGh=O7F5M zbAL=7rIAFb=am7Y<b|I2O?R=b;a$62VdESdZjKY2rn&*MT2we~61Ajl8_QUJnqr|D z#BZIQe8W0=ttt;s^=T5EUl70nn!9-ziJ;@kpR<~(o2ctNDNckH4elDO-Z8dL&S98Z zf?9u1_+4ubjgs3ePEEO-BGN5)X4<^G^{nlZmJ>#X5wd~C;-cwxlkM|XnNDE{C#}^j zD?`ypm<=fPXMR@E*86AGej>@)I;V1{k|vB<h>&sG)i3>flf6$(Q&5|$&j&{iMC0Qm zWqCb_E7}T6@*^VAdRII!Eup-jQwyXP?<Vrv;ytn-o_tS1e^6=mbI2p$!oEyLgd^;K zXhXNZq<I@Xq^my+mO(P~=k&;x-d{tJf7KJ$^{rcZg_X%F&*3;dRT^(-H+Qg|B#de3 z8`G^C#wS7|Z|vQ&+}<3xw;Y?t(qrUfBt$`*$C|@YRnaSZW5UO3tPZ*q5GC{x%0|hw z!pHT|?CgX=G0qg9)|K^!u<+|KZ`h13IXY}G-+cs?vf82`cI8YxROFPFQ4US_MikWJ z-@<P|WY#K2bj0pcjg^E9{GhIQ^p2Xdh1aZpb7?YTwtB2%Cxt#B>6#eFeE?}O@)+l( zp2k*SuZNkh`pfk5-mTJ785`uEk0`cEUzA$uAN3St!&M+sFN4wu&oqOv6g6E072$Ww z;1Z~%bN}R2ogtM`+kPkQ9yN|cHK&2dXTA|pI$T*7UpiWm1DPVM!93Z_wo)u>jSM@^ zxd{JG%goTD$ilN5Qb5&O%KG9+FG~l8jKY=TeU#YHLP@t9%-t$mp^OkfSt`H6MMy`P zZs?^0Bj`EL&Igyc8AhL=KLGX32tJgtpLaiJ5|{0Yh+yabPCtpQ!l{<}eKt#l;B|&J z-;uH_v<pruceMX}J(n|EMuV#?H6+})8yk2AuyeH20dWx6uoKzXsx%I6^|nzCY8;a+ zb8F&@HL~KQWHd|ptp2stdLXeM!t6!J8WF8!8kCXxpWZ@QrEx=uFEi5^YsfpG;j}ad zBK&@QUJvTG(<8BdSmycZa>o>&pD~7F2?Y1=bZGDFlyJ1G5FDlwAIl-0)2bbK3a<GG zpjYGu*Q6{@YCHohix_SqLKKNBzI~Qf<$SvK6xd|(IGY~<Nkweh^JR${49^l!MtyeS zJR&F#S=Yo1BS;PVV?g0p-e_90Y&X(<S+|}v(roKLIgt2MO`PdbVo}hV5cV2DJR1E< zaQ)E3Dvc{DtxeRaLlPNrJ&S(b)ayhpv94VI2e_q4H%Zc)(h4!)vEtGzUF)*hpm}5G z+)~kvNAV{BHNsd&&i86=tNN^N*>0=ml`EzrHKAFUa*bsH0yB1|P2uQMib2(L%9b77 zdq!!@!F_g=F6FIIp3xQ^%Jx6_653HP*{)QWX_E_lKE)JWLck+1T#$wua7+ha_>x0O zL=S2@es6PYEg{av8aCycvUV!xxSm_!RllUs{1d_4uIXy>@S}#iYWHgbTi`y{M7T%9 zd%-&0S=3tw*oFh<)Z!VVmCa~O1*jT@pf@dV&_A~teRFvqoJDJX@mv%3;TSB?BdXgF zR7d@z3`(RHFC~`D+BKeT>JTB=O;~zBmuju_gfT!c=L#vEex>E`>z+SqW`?c_B3r-z zh_adI^?4AvMXpXQ`iKYe1{`hwaW(tE@584u^-R@^H&BZ`ol1?erPgv6<A&UW0^oby zHaet4)wP}j(qM!F5RE*JW_A$Q(bZs~YO6K#Ag*5l%;ZRG4+Eao*QIAvlmS&^o`vo+ zRkPHbyslctF0yTtC(&3bNR4Hw#uG{S9~aKXe$jJ-KB33`Y%phV&v^J|iRB<@v+c*@ zs6@oD><9W29zj%vCa8P!VeyyXxtwGCaN`l$GFMp!2(zHvFl0Y%^p!xQq~Lv#?D| zFW&fu)vq;|6v`hlXl@VEddEai>liX>)v?;~vnm|<49*gvd454X>vUZFI9eGhJG@Cg z`Y6tMQSYZ?NJp1dMsB(c7R>62O!cLoaSai5h>b9w2{5Yp=oL5xLvG)Xp!+nBAuStS z)v4UVy^4PFkC3goYaHa>;Ii3N?;~$I;ivP4FfZLRwp4O(W*8L@0^5@ivue%{=fGBT zd&&L~TRs8=J>322n>XJ5c3(f_R4inzNEgfyz8iD_3bxhe89rJ66l}9pC;41k-A+mF zYo-6+1=}^D8lOKg2YnV|`m=j}vk2*r*P$wg-s)r270KLsY;9obbgtdM5yyHaU}03e zjakQfd|Vx@_ECouCuo%^@LcPy<`uXFs7X88-d7BZ6`R{JUq@5apc+<VeGSeT^coFW zl9?!R-ZZ8!ql32?9rY?1yNcAL4DPJLxeI4(I5MVcWYF6yEeug(R_3fu9~}GS#s*nU z%igy6m8oJhXdVbrV`H`U55Y0IKfJtDZroHxAJsLPU{;mzHf0FGXz=$;MmhJ2gy`lz zx<A`3DYq_`M{;JevM|`N1iV{l;j&>l7^qCJ*cq@^wS$mW8?tukXm*|+lRG{$G|xh6 zLGLXxr+2CAv8H>Ha$}jKA>~aoF|A%XG&2PXw500(nX5I^k%u_ta$ciRvlV(!VlI7B z^);pyb2~Xpgh`SY(t=g2-k=AF`}mLI++u5o<TNfs=Nxk5f1G$5DXZI4c8>g-NG=>@ z)_Rmo%xB-0;|#(b9t6LA_ruxX>j$7<!S-F<-+Q&jz_C@XLe&mMqp-KM$fz7{sLiOn z=;jfA?bZ+LR`)JN_XeRGYtov#l~_v!TflGJCK0IQ#I`wqZyE=Kz$iB6R3w3rPw6){ zN_HyQ`AX6BZk#l$q$Tt)-w6Mt=;57+Z*?ivt(MS;st^GUw!7>uq+&*vK!NAlBSp3S zmso~8EOaz2g!;h5pqR_&sMs%dqPtVyjn?-$H<lsHoBUFU{w$G8ZzHuc=fHiZN;$um z=bYZFgQX+zW6ym%m=<m?9%N7BYT6&Qf8_f(fkO0DJ^UOU7;PHpfnPFm_KYD{<x`s` z{TAVRTf#c_DXNHMOG{`_mG@mhd&Kz~)$N=n|M+Q?uSbZh3vNzVOmqdf<EnB_Bs4qi z)7vz@PqWp1UPN@LXLm5R+dl^*RJJiOdDhBA@sT=a_;2OgV~}ME<;+jEFKx1XGOpJA zd3Dd;p;tai?RB<3FXcNt@9UOuM9TR$Dyj{4D`^z)x@+t7Y=_DEfTwTO{)}W<y&JQr z-{V+Nu~L(w(pm4I>aEsW-{#YfdCEIcBRaphR>a?7#^CQ*zN;jSXgsBvNjPrt5T!nq z`SLSqu7tSB2X&8z>9I%0gRb{v4R0t0bKTYW!Q-QY-04X*9^uCO1F6YSoW-6Dbmu!{ zowu7!voSHBye@_fQ^YQ?hr~De9!djK@M2=TmMXaB{f6Qw^<*t?a6%Inyju4X`8C8d Zv3Xp3KZlGO?#M@f5qXACb()Cze*lR$l&k;% literal 0 HcmV?d00001 diff --git a/fsd/fsd_v25h.geo.info b/fsd/fsd_v25h.geo.info new file mode 100644 index 00000000..9fbd2807 --- /dev/null +++ b/fsd/fsd_v25h.geo.info @@ -0,0 +1,43 @@ +FSD geometry v25h created with create_fsdgeo_firstRealistic.C + +Position for Au beam at 12A GeV/c and 100% magnetic field strength + +Number of SMALL modules per quadrant: 5 x 3 +Number of MEDIUM modules per quadrant: 4 x 3 +Number of LARGE modules per quadrant: 2 x 2 +FSD thickness: 5 cm +FSD front plane center coordinates: (15.1229, 0, 1000) cm +FSD rotation around y axis: 0.0171042 rad +Hole in the wall has a radii of 10 cm + +Parameters of module fsdmodule: +Size: 4 cm x 4 cm +Thickness: 5 cm +Channel Depth (ratio of thickness): 0.75 +Channel Width: 0.2 cm +Channel distance from the edge of scintillator: 0.5 cm +Channel radius at corners: 1 cm + +Parameters of module fsdmodule: +Size: 8 cm x 8 cm +Thickness: 5 cm +Channel Depth (ratio of thickness): 0.75 +Channel Width: 0.2 cm +Channel distance from the edge of scintillator: 0.5 cm +Channel radius at corners: 1 cm + +Parameters of module fsdmodule: +Size: 16 cm x 16 cm +Thickness: 5 cm +Channel Depth (ratio of thickness): 0.75 +Channel Width: 0.2 cm +Channel distance from the edge of scintillator: 0.5 cm +Channel radius at corners: 1 cm + + +FSD size is 192 cm x 160 cm x 5 cm +FSD contains 372 modules in total, +small sized modules: 156 +medium sized modules: 144 +large sized modules: 72 +FSD volume center coordinates: (15.1657, 0, 1002.5) cm diff --git a/fsd/fsd_v25h.geo.root b/fsd/fsd_v25h.geo.root new file mode 100644 index 0000000000000000000000000000000000000000..c699b1f36ab94fc166f2ef329003bef518b7a800 GIT binary patch literal 17480 zcmbWd1ymi)v*^3=Ai;xsu;3DcyL)hVcXuaPaCZ;x?(QVGd$1tE-Qf<&_x;~_cfEJ+ zU1zT~vwOO`s(R0!{&jWL*xT4Rf<TwCAP~q91QHShfp7|evM6wX1BVDCa54gc=&L{= z^kfhi43Ev&$|TB3zU>qkI$XyI(E5M%@dv#sl67?<jsr1+z<@Gf4FrNx6fv_jLNRtQ z)Ny8FF{LvxvZ1qo9nt?s6$t888aNUF!2;kA0ggXffH8yppN58w`#ORRcR)TK#QJ9h zJFoJnilRm~$~KlxRz~~|4n|h`mTqtU9S%rc>di71$bAFk3O1SnpRTFo)-~RqMuQ?q zgv3DNN6aPkhMY5Qb#x$=TfB>pG#Zof+;){TtxAuus&#|7JOg7wURo^c%r6vaX-e6W zTCq+^B)%a?ncAwn9u79dcMgtf4s+9`)1iae!FX5l*yZBt#^HR~;btt;VaKQQBI6~a z-9dIP^~-$93~)+zwitmW_wzY5UA7phJK9>lub_6(Hd?1RTD+Jp>~)!z3av-={KQhb zoauM6&TfIY@<hV*NxQ~q_q@9L88>Iu&#Pec#&iXuHJGsNU1<dlR=ca?a#!l}*XnNE zntRTwtZ6Cp2zFQ<nAH02`EafGU(mcPg6{_@pc}7gLz4*}Fq|6998Kp>yUx$f{5EI> zvAYbA($K#LpPtGf;P-4u+Z;0dKK>khyEm<5U|qX22_X%WC<P*QoPI_rvZc2?8H&t% z0*_NSkb$MQ3{Qcz!zWm6urSTG@oe%jE!x*zD(r`W@RIQ&INJ@>zE*dzRG^QGB*Gsp z=?!b_UH^z6znn-#z`txHEee|5(34)aF9L~9%*C<D`dG?aj69w$4J5N~X~2k*Lk9=l zQ-}<z69~EM%#!K^%n=Ho%bRMk7|(6b#jVGAm9KB)*QU-Z_$rrSjm8dG4sK~0%HMWC z?cMv8qdmRnq+A<VkC7Qx-+XUbcuSES!7}XazRv=&4DIq(5ObK~wEut^21nLo^|s-y zT1ZXg)SAiG^Pzs@i?2_I+j6V(k5tjj)`#6o&lPo_ah!`UwyIxPr4Dg$hXgK?&Ayx? zI-tj2eUrqELiFB@`@tu5lyXqf@vBLvO;;A{=PIYmuU+@7cDK2qj#Z1JJB1qp=C*z8 z)?)))9;f{bjl=rd29u~?NZwZ_7CCwsy*=oR3U=XAxKCBT;!b-%Ry$~AXcLb0WwKm- zd9Yv||LxUqd3UpR)i;oL)7rRpacXpQ)K2KL-O!5Rz5ipzWI-nIq-Hy~($L3s4rf~0 zuGw}OMEWRX98c4m4av^0h9m_GKCNm|Wj=MnAK_i*5-Gm%a1Z&(_#Dz1A#m!vR%yCn zI(0klR-bu;&z8<)(&B&+?@NRS$wScQC3SKTZ+Ys6ALWccS=D|Y_0E(MBIFSeAUqaj zIt^ql-)v*<jJoxVX*i9iH{xl|wr((AWO-iRK3(J9tg>Ge6Y>~lwYqLUSFSuYna|ae zU6Z~c<FbTh9Y|dLR1IH{>G{%fvSN;TM_Bea+`4g=nRm0lOW5`lTI)U2`a%?Yy#`)) zGtxY)WI?&X*1a&$hSZ_`W?aKSS>cm4-T@&t%pGE6q#>-@fLzv>U*Qbb7qtb>{K~m3 ztE?S$18udTiL8yfoEGw~#}m87k%#?{xL{&MmIa2M$>xhuiOQpLvfbC8xp<^LQ@eL{ zq!+HFQuJr$403U4AwrmoA}u6cD5_r=)*_1%A&twIr8}W>>*|cvHnBu;{@{;tb4_GN zD8hfOnqEw$sI@u1{mkRBZ=f)*J~lm@T8PhCf43=kSJpkfnz~B~i}b4&dOUWHYnqfR zejJC3yWCBiWFyBjbELNQx2G%hy=}92Ib*)A0>p_%+M3;s)5m9wDWbk_IW=rKvV$Ty zRRicCeS|0ZLyVRY=B*Fklzg#-Ws^X15~U!<NFEVzA@*|cJ~8(o$dFidsF3Q*YE+^! zc!jRzr4!@nMm~AwR2~`d?ZQ?><gNV_R3ctQ4!-ar6LH_Yz7i(G@M055-yuO9KBIvd z2U5Om;&(&?GlQ>HRD9vZveBTJU`1ang-$8(0H<z2o}@=>Mx+F61<v{fPYRXuXgFkp z-${}rLIm<)`chuGkVCtXW6_KS-hq%ITzRO(12f0^0%Le2NW~z=f=r69$Zlc~MLo&) zhOj^hDunUs!|PB9&INZ}!rXv{7E-XkeR-6?|LWw{`KREO3OPhkD}%^b;QUSZBIAdr zp(8M^EaF$4Mngd910X65E<sqb*&lQ{y$GB{UO{Emw!qNpU9E@3YyP|RV+jt9lgcWE zCHxs$L!97N`yAf<3LCTZC9Ji(pfX`r1P`7RUFBx*mtpAXdxvGi4~k9msgU&N<R&RI zI1cKrb_OG7lY<@v9h1K&*3%+aJ@FS0V6<6`ox^tGPUAQ=HhRNMSVT*Y<1im77afz) zGcEX4Ulb=8aUi^KAhJ4}e$pH)e>lA)CbU|{^-kA{{LbND={C;7#d`(dPiDLU>Pu%= z-l?Y7Ubq>fV}EtQ9544!(9dM_4sWE6nq3`O9V|MmiC*n`9lr!Yp%sFLH3cV?(X}9$ zx}Lo2s;_NLjM$#hOu4n|)U1;SRW%rLo{@ukKSHzv>Gd!6ZR@GzJY$iBOe90-l*vNJ z&;9oBU{QX>;F<@+<g*lr5%d?{@okT-OqgVj6mU%uhmK$P?P+pNr3=N2atNJ$<w0&& z7Z7sG6GYG2hnN(2743`zPZQ7!`z7aDf+Tbdp}D6J3v`VryoVNaUiaG}&E=9O{%x3k zv=;wbo+8#y*NW*28N2KZiVy%BROk}T2(Mk;5Xf>n9~x&Y73^bcMT`1K5kp8zMl!H? zjBhVPLWuiXT#dzwiE<Q)u)9WuMah_V48bzF5epdtYyt9>zU&tvo_T=u1IueaPaZu5 z%5_(`+oCYL22Nm1oA2=f9hLnMqvAq5uY+ip64<^N{Ef^wcAS9hhWcuOLHpE~NtOVx zqTk2n(H!ts;j&XEv&i-zacVL~62iO5-h$6vXXs}rTZb2<LYL2C!6A0ncc>)iF>rFM zTbA?eeca=k?+Nz-)!JlVW5!Kif7AmwYz&1QiGI3Q&Z}Ou_rPxxw%<XyGewoV=bP2) z7^p|uZxYWLMAbYSsqp9t+}hJ_&UvoJhwxgB_#RX`|Er|c-RqM!SjA5)B$Y~B<p$kl zDb}<8gS48Y9+R!G`LxpG%`0`?H<Wx^&Sm23sqETq+Q}@kV+_Ss&AI`~4Q!VKA7HNt zUAA{D#I(PEM&Md18Ku@x@1-CdQ|IHItVu6}akcK#Ay|uKv8OS$tJFgVtI$lU@y;$p zwGL0vY96!p{}Bn*f#EZ9<x8gBOmJ!X6ym&hg&->NPOFIr2TT6;JGp|h-1XIIVE6ZL zJy)vikbOsb2dhdh8yw?DnRG}kTMk@7ln^a%JnHXO+F+@v!TJktL3IGRQ0zOkI5;TF z50$wIb8Zm`ebdFVe{czkpvZ%T%ASJlp>mC-^Pv;p@IYx8L)9hjPQ;3&+?pl+npi3u zERz%uQ1soL_6V8afc*?hjeB71zw9@K5&GHF8Hvvt9Z-aF-uw*$J!=JG0{T^vcafv- zrh^A)i4yEwFdt&Z)!Ihvw<H<g#~+IJA&m73l4Hd0|DekI3#kBEnPpj2{4zDi4=e?n zM9?IS0)V(D5m@9=rnlJy12akDzNU$<T<g#T-{x1eFA8S7dWTAsY9xw5+!)zG@ySUw z+ch6l)BvJ3*{#?Edx))VhY6L&_gxFn5=>$rGpeuI@?DyQxP}My7a;O^EwDJIz!d$k z<2wTSLiy?|z+|%;%ug5*5Vf4ZFVDc6YJn&}Ie~r;dguJ$;`@OIITJ}xrSqq5$sf?@ z=1W<;E5bF!YRXe=00@_(5=s0X=JP1Iv!S%-!w=I}g2}Z)CHZcSpSv*YXUZoaWL7SS z{(}SADhXLm`v94-VU3tLw~q%OH?ZeB-Dd>?#+D8V0il#~!M<Df&wmSVJPHFsm^6Yz z7_VmmRhIZ9B8!5kBxUd?$6i4Xt5S=?2*Ye1ssh_4?T9Fi&hA12+oF<RfjLt2`s6d2 z<A>K)#vlW;44A<>)WjrzW(q;vKQ&R3|IBF$QXy4L&+Uuh6O2nan5{#t*BOQgOu7S7 zp}*5{UNXe<c`qENG4vi9HYO?p6==jI0~!<LsF}yS!~U7T59F{YcYh{nrDRAI-R)3R z=1<zaZ%kW<KcoLK+NmV)mkLj5|5_>eSGti@C{u`{w)twNlNK0kabHmn+k-Kn(v)fP zNXg}C@5nW)M<Q-%M}sRBQ}o*DBW*DvHJ>%(A_Rx3`z?!o{X5m>+0Z%xvFdt=sq$MY z_L#>A&bsVLuuZUc0-^Qz(QD-V>PJL%(fT$pd5x0Md!HzLYPVuUK47@;-Q&Ze{)+8f zJ`;=4cYw)TCKZd9ln|?KftWgxh^9=6kJfjA$#W3{ijL7u2{kA2L~?X@4~o?v(C3m# zP*w4kIfcv`A-~#-!+!AofggSS+A}^H%i|{<f3Bvyfc{~-ivRYzMg3n&Ihy*18MJca zUQ1WZ((qraZ%!nmAAQY;8LyuqU{P05j~f?e`kN-YfNMcmR5x*0pbMdR^pGTn*h=am zZz565AH^qe{}ehU{e`s8*#o6>s`|(5GQ!(`((;qge*pA(M?rpV%M0jRcxXo@%@U-= zD0z1#PhoJvNh-#}p)LmWq<HQ57SG`1MiMX}0aWo40ag1%cBLl3PmD==1nDbKS~Ai* z0i((OSDPRz>3J+Fs=*1;f22tG&&^i#`_)mDoEfh>B>lI)eG3Ha%94t0ljZvv96l8S zGE;_m+_R52EzR#f;cHd5^^U&ykQ1v%C1}tJPzj*7aw4zN5-;V6&(wu{LVf}k1^v3d z7&$Yeo1eFG2a%)rA<qH=LX5yl*$))ufg(p#?p8yCYB=zB=0-_`P>OH}RX0FPa--&= zDqGioD(TV;Q=a@(R%chrV%32}<@sEvdan>bndlrcsUxocI#np0or&%(U?+u3z74QF z*pLc+oaai}>TQe_eqH?r$$)Fckfy+x8;KAJIx1S_M)Z(L)Pcnt@etFByGaL6VNomC zt;5a)<VnqGV5NJ4g+&=<d0qZ)B!M3L$o(gEU!r(E5a%{K8wx%1+>lk*_Nx2Jiicq6 z8YwuWeO6y=&I-{tfywRc_0AdRBJ4`S7aFw^FbhNJA43%KiFyJS>yQ?_|5BGV@Y+%6 zStb&vp%`;QU0;mXTw0P8V$!6_V2>GTl1Wq!(7U8BR9y=(39&$bnGzhr<YDcbGbgN% z&dG^9TLehQF9LT8Xg4AIPw~qizT_?}`t3RVMa{rs7C<WFtjGw!=NM%fSVR&Jp{#WR zQ1bHSFO}@Bd4aoL%M$W^MSa*O7U&FMZS2`U;LQP!U}@u%wR4Wv@!>geJ08e{yyvJ< z7{|N;Sp7v5?<;-0!(N}-*L&UhuY=J^`46X_WB=uX&MRqN@2XBBvYysy%s54-219Lg zeQeafLNC()Qso4dB#T$#_21+r72^3M0k|xkB(EC7akdU$;{TzL&wnVS^NLfR^KXDh zd=G$g{$Jv~{2>ssodsDZ-M<40`gcHbfQK<Td<bxD@U?#i7ULKX5ztGToXYk<tl6yE zu<0i-UCf}AKe6}(1B;rwiDEOJ1}70W9sl`_i79da5~7%_ShVti$cJlZzWV~m;h&?q z@A6ew!1l%eaHxK?zs#17r92ry>}M>1nh(VNzX1{P|DdRl>@T3&Af}j6Nxzf%$dY^t z?7?U;sMY*>gdg2O`v*IrKkZ@iPX1-l3%~;aPykO{$i$vklrLE0UT1eV0Hmi*R8rpQ zIsjNcvi{pJQm|h&nC{9Thmn#ofO)?AI*l_l)khW{9{!#Cgkt{?!-w);+-UzlAZn{i zuD*$S6(b`ZSbO=ynODM9Lk{mEKdGQCpu=g{u48$;2)&{rppQ1?vG^4)%(ei)e8~UN z_}u(2a-aYb4kHo0{zS<WnSr&==zr<=*Lr<}i@%WCUeeBOrOhkI%8xSG3ev@R$}Rqs z^J4ujDAiwU>;8~bPCrqKTmcwf{;-%w$<(a;uYz}Rjwe=3N&<BR)M@I&R>gHSt>`Nf z+4d@j$9yx9tH4CuStcLM7D6xOM8+p0&c&j&GgMsCXAA|vQ~=D3SNJqROiGdpp?|xe z<gm5~>!-YSCgK2`_*Z7t<A-Qd5mV_IDYOm&%KP=vH8_!_2d#n0>@d0;g(^pL#OD|u zFcg)Dh*GP}ltN5={w-e!y>kPRkwaNP$U&Mcss-Dm^`;c?zIRez7w42%tc0qX;ycx5 zZqHV&b*bf1Hjl)Be;OB2Nm4Cd{q<e3KOh1gP;G&2@%VG#8Hass6UBJn_&t9GpEGR^ zoRl67c~Xb6NH5J@&uEuLcBW3533;?XGG`obm*xJY(PE)Mo6lnSkm3&>el0KpP$~K$ zH?u+>-+(1|c*y5Sat-gbL!8UW4j%RZt&aKiwek|~5dD5cux<q9zjQ?X(@}k$$eiwO z1u=f2MJ)Uke{O$3RTQ2Of%->Z%`0FiL|z&7;Xhz1A|=ipnFBEE1$Ych4~~|=GelXN zQXf0)35#OX@On?CN&f&XyjE)$`G%`D>(zI%yT*>i`imwW1%UvkdS%HghB5E_fT3~s zjCyfyU{BNUnMNP=X!_w5Pzj{K7Nqm)i|VP=axdG1<-Z;c`>VSGEf=UpozgUB)xq`G znc`7FegD^2QUN@kiYfEr?*OOZUxAgiK~S}zSxN(JfW-cx(<<<<-hHe2R9qACx>riG zsrb7CY9;yl*KcZg5Mnb4|GU%TeBA;SBQsi%p|K*f8w6ZRl;%?b9$o97t&!)sFO*l` zqt>wD*MF$gD*BI=qk^b}zh0OHDdRaQ)nBpIBO-PK!1ekH+@^U6eIcGyDMj6S#r*e) zpH^+MnpB3d^WF_r)~T0Fp#Y$|lE~dY8J|#tQ8ZJh!;U86IrX*Q->?y^DXp$ij{jcW z2r*^vnUF5#$QssJKl<um4d)l|0zEa|w}o;w&l?-?Kge4i`_`d8jD8hr&I<sh4N(Tg zUns5OM~@}_h1xlym^F<(A|Bv+QJ_RY!<lMPwi|fB<S8xCJEn?)zxrU|?ex)qIs=Zy zXMaG~mp+;lsTIFnJhqmH`JYmuq7YDZbg`1sK1NL2QwP;Ldmm-#?otYB5sP{r6Na}* zU}FgxXaNdYkcv6pR89eol{0pl!}?!;hdxH=uaD$W{ts(BGX5>J1BGn|IAG5W?+!o{ z`qdkQ`Krd0V)^PC#YZRKtR+cAdm`@BNS823lE+D-k`m&-TC({=C+NuP8{>aHI9#%S zSeWz*Aa4-hI5s3FVK85JXadxHbr8J*|FQcZ@gGBN;{Myy>LzCLfq>OZ@n}!Sy6DZQ z551KTHh=s;z=TUDAoda4OjfA{`0wiUnO597g$1_PH1^99OB@G|5rs6!gpLu)g<rvC z6^ASrz5+wHz(iekIL6eNTQ-HA#)<q^w-X3M{B->A55$~dQTm6p>U$H_^!`$8$q(>E zVzG*mS;`X-gjB>sBvjQ*>J<YusM#t1JkqE|J>>~2AO1p%lQVn?*sr*>i(1TVvaa}< z0?)t0XdDnHu*_+myG$a3n5_Eeq4f|FQ0)260X(8w;kBZ#{?_9jdq>af`i=9iFEuX# zAguUxDF#;COd%t{?p<r_a1+3~rEM&&7a(L0MCMm*m0nmNehd!61T`Y8CcAA�fO7 zk$%n$_&Lr1`}jpvL-}?P-GH^#ro1)P5I_@ls~1?S*)3z)qf{0n`6AT#tdB90%N!kx zBjdh?Dfed<m`}aeFBhthKI(Elci6TxKZ8~;lT3A@VE!aiUXwEpe^M)dk~2U$1xStj z|0ZXEMj$uzKgy^8qSt>?K^=b*F>J4i7`m2MYlyFj7(p8=eKSRSJ!^-zuSpmK;PTHE z&;dwEy^a$7zl`!VQT5-95>SnE45)8@H6kSD!V(X9c{%);Dj5u>2~$WhnJe&eSvc;E z=LK)-BzP~AP!S2`x4P}EVnzh7@4KN}BmzCsgwOGKB>B<zz(KL^DhCn^<*l@VVj<p< zh-9T@W>#wIC$Bqb-pH*Hyh=jX8`Jx+&D(`Qg<Eho`FOB!H}`rs=qGUOCT;WyLWwVH zo9A3E)TRY0cIGj=?lOEr25ct-ix*h>0^q*Fj`_lY22JD9eg|%T1_c`kEz{bHP7M}( zD7>Eu85q9yFw^ErHfi)!cnzQGsz~uq_33%9e%Z0v=xB|Tga$FoX^$dHn)0re!Izd8 zX8TcJxYB?L3B5mdU^!9mO<{9u%#;w!tj$~YOB7D$jyj4^-y5Az58a<72{Wx09soZ^ z#BS;!wMEL`u(lg)r(97vCqX+bRvuxf`6gnjee(ql<)|-(kVy8~MqNPa9-er|3zi1? zB61&otRBs~9xmqc5$8$vYQcq(=f)kkoloIjpF(<48i(dvw{%#M3`7$Z__Guf8;Vt$ zb(E*i!#A^}9?Gg8F>3w|$-heQgG}3hyFTl5zF6$j8*C7QoC7Nx;(a~EcN}HGF^2NH zV+kO!Ffb+p1Ggo11N4<LMZcuYy@Z{+)8-3(!UnmD2dvvN!IH>-cUO6fV{u=BP3!yF z>eE$@G}hWYDyuAXrX|?5A;tofO+TYWM1L$A16$#5njWqI<V6})^4ZiLRuHkt{@p=h zZ1V}6&j=*WZLP`uF5XM7juj1uf!-cXVu`f-?+f0x(44nnd!vc5t@9j<0Tm=J#5gN; zA=>1&O0v6OAweisRuD}tYiF1ll%}bY;1gwRdhDmaW5$hbX3v&YU__^N21r;#KD>pd zO_^F(N|{=P6<Sa57jqBizpR>oKxMtH9xJYz)TL_ohi0|wq*G8sALGhEkos!F!Swc7 zP@`?a@u{m$Ix#HN8E=Q!+kFOI&;u(&W_<*2sg~?#|1_1V&`;A3=p7`KZFH0iX6Ql5 zsRnY6BIru%^lmvS=Ox-Ro-+Xi9+J0)ODsf+QEkFD{l*&lODLKm#VG7>X8pdW)u~7X zq(yl=k5j>GJQ?nJ9}_DE&V{?})vG)iB7U8I{I)lnp5ZCpEfa<1<L$6FBlMVZ!bsZd zuUL=H{T<v)nFwDEYQIjUd~FzWTEw#&u_^2EbO_fs7SGH!KM|ajVtov@31(^?+D%_J zbLDJ%*eRO`T}_;CT@WEK<uNe)AuZN*D;-HG@Mdt_iOIH|s>_E1VkN`Snp@Wat1pN@ z8KxAO#M4*KpNpy?HiY@@*4o9w_nat5MNyi7541QNw)&|EniNpZ%i(B3KijKx=QC)D zb<V<-4<eBG6L`3I%G$$y0ui`Ihi_ACO^3J-6x#+iWv+$k-#0Mk$h#KcWj9pI6dfZR zOyekvL1)I}_>f8u65rIRl(Kb$saRKU*3hYBa}&=5DuHLrXXYl4Ieey}C@Mz5el+=b zOU87ns53f1e3RRGsF)L6+LiPuQr;6spCe3%=ltOy^lS-hWrBbl6~WvejEnZ!7GWLF zo}#~*D|gg(4>D6Lf*Y|7fwQRuQ-sEFNq(>5paqvYI2ht_DC+R{lYPVnxUyHySGIJy zN#n28qYyQ@h3SZm%z7&!9&2{8zHO7oU3M8=#R-eSTWSKwnzpJW&w2S=`-CCO`)>jh z%*d@c)Mj;r9T1nb_b=J29~@aW>0|wjnR`mZ`?2k3i<TEtdm~$wJTl3?b&6pYVW!gZ zx7$j8?xaR}lO4NYF*SfSgibz@`M$1NW+g>cDvTQE=5op4BQ}OrLyP4zip71mvK{1` z)V^8V+nDi8n8R>FThC0ifxsd)C-m@WUj23H;IlTBuNdgz2ak-r2*{<jRry7ooiJC6 zUNP6wLBf~s>r964Ze)JLubzJSR_v34z{O1;dtsY+&sFf*$7GrPdfY$VHvsAbpVv9q z$6VQlc>4@w9H%^}zO&ft(lQM~K8u2h4jH%{6vT<g=EU4PHEen{^LIpDE>)|BF79(8 znsri9KAZardjWQV3N$wloP!Y`I(Ov5-g&SQJo;UyN(6?!k=xnAVr#J5i{|FJHtD<x z+qc<V=&lO!=J!Z?w^=CGl-XB=fE%f3?1&XvhtWd#3G-JBZF2xpTZrxWebuQS_NU&V z{oelXVxoMDew`MuMlc;o3C0DMOm0&bpj-$YQI847^&^>7L@*=cr_^fJamn34Yn4OV zJfjACT@Ox!6ES!`4!Kw%RtXT@_s!qF&c~|b4@TpWRRYb&!tkgRRH>Ett73}zL7P<r zv$oTbRE*c^r$QypG*lb%>#QKy>Y_Zusp=wZ+)q&Af%MDYz1xew!F2f@{tA7Ey#3?W z!B7?N48zWglKYtm-Kbdo=d`xpgQU#o`h!Ck7`vUXVjT&v=VXHkp=X_9LH+R3K~GXR z*CCwwdfeXPId<${`;h6}joQ7Zv1Z9s5dt>DSKqS*+g)>osv6^UR|aJoYG~-4bxg}& zpD*XmKA^SnEyfr&^Zk5$yUQhl#TTqH=_mcG<TS{w+q}J^_)yFZpAWpOyE;Xyq`-k@ zXW<r>pynck+3@V9y4Y}GRo)|t!a*R)>KOZ_1ufiR!mH6!XeDCYlA0Yn(FtSe%<#Ml z4=nM`919bTObCNRR3VtP-S~o@V*tIKeDG+s`gp>Vx}gV)agmmVNsa^9?tZ4(hg~=M z<Rqa2%h)YXMmge^;4V+<;G2&!b%uOfsz|uJJje73`@EfvfeaPG1KM@L3|Py249W^$ zyKcwceiYm*=x<gYij=A4PAwRJgR82>e0*S=UUO&&rsEX;-0Ne9_=#N~mVqYs2W@}! zZb()^0P#&B7-AEJiStM@yz?x#x|8_@>P+=>S+U@kcb#BoebH{9k>lj{deDyZh!tK} zSgH`_aF?gzz?N>Rkl?Z6`lv^JxoT``hbL+7gRT8n|NJ<&(VLMRGsp(ajwyNNUlG>I zWX5b}JUDynIS?A}S{0%pL$5mCVC;Nc9YbR%U@jH5P$iQhr|KnXrBwLrEy+|V>LV_q zVpa+43`<0Ir}~sxc;L`+)iI`!)W2JM>MS_e^mrh9a2{+0#24g}Kiwk9X7%C0VNA|X zz~JRJa<G%SM-fItqA2N*C|~_HSD<89e^<;y*}_xVqp$W<Ix&lc+;Se<Xd=D2F8R`$ zCNEm{0M+ox=coL|J#~soERnHCGPsO8Lgp;zUU%)ieWlv!G10qeCIpu9o+uA#$^6mO zE03QfddeF%vvqJji|n*#X2~Z<nHZm&ilDq>0UW$PpTZFK-)TZiK+!r&+Y?Q&$7*tL z)};^9IPi3cM$+fjhqyr;<Fu>!e8KHdm%*ncIsk8|Dn!=JJuK(&6Gb~#Jk-iYE8lH& zXxAylGZ=cWp<!^bI}{742fk{vvrecb?1fT4I?r6KMTuiZBNE3?&?L!5M$gey?@VJ4 z{EY5drVuZrwl>7v)E3P}CH@*JWk_x$bC(FS8io<aL)dCd?NV9(M1l9-hJc7ipA5go zP>vWy!Cfwfcvabjslq$I=j@gmfjvZFn~k93VvLMI@^|6G*m?!HO+vzGS(@ldv#@Y& zFpBFSGTp&C0?LLa9>zc?*p#Sk1vV8z@n^sH$9^-KQ#9c6l*bjpD+^0Q5vv^A-U{z; zl$d;fhl&Xs4>KL|oVZ@V%iRjO^gzkHbM_o&Jn+*IRgRiIFfN9Lik7$_pFEYs86AgD zEpfM}p=}f6SdQt$;+eshtY3gniJS_*SCuO7UgDiIed@@*O8sm@7cr-GEL_N?mTEtU z)$AaXow!0c8l0~tiUVWJXs+V@>1$@ztaAB>kYIn;;GXm;r~TkKf5kGKp)&fpr|5a6 zw&<xE_*iwedBtQ!HModi#{?HKCAVb)Fk?D{a4`*Rtc09G4gK99o+YzgXD~K{Q=U^7 zi#e9XH3^Z-+J-dw`NJJZ=Ckc=PQTEFbNsLDyHhhvr=T2%#d&M-G^D&=E4G8qAqk=U zr=(Y3HEzBg89IcI9}t<I3=D|+$#YG<V=8OnOR)q)32LIiM_2#0qutl*`ID_IG$;Yc zUjK3kV%g4Yq&7)ks_3d9au{Sh6;Up%2#UL|5D(<*PAL-htrayNZ*$T85<QnO>Flzk z7Es<ZSKPH2l;`k>U3(-5*D`EpT+gwY37^gY8NJc+JV6-nPP*0?fm*xC(eyi%CCFhZ zwdbK552U7|fLQRQsW6ix0dC(=tShW;QkudguAM*BKrAZT(Tf+utEo7gP6fKz^!xLy zpU<cYJHea>VHC)k)xajd-ldPbn`JubDZjTES@~BtSt;JjBSG|SdgkL>Cf})8*L2yu z7%b=dLGi6@Et$f>{r626qPe2R`kv%*%}UG63sy~r>+Zs*-wjt9`S<j{pV(w4d)IAd zH$(+^1ic@GeIT;E=4z12;vZclX_?Fwg|?H$JZ0{1!IW66XE6W5!mJx35uiAi`?h;f zF>GGl+-c)VRm;>|?sk2WchC`fLfB(<ld6)_ANmBXlJi@_10FWUWNpu|vl6l?&hd(~ z9f^9$>1>O6IVRSsl)`{5*6}uhQlsv;oJD+T<dR7&W>4%uhe%L=EU>pKuJhLQ+X(cA zltU0*j^3=y0}W5AdRV(*$|il>h6lg1WNnGC!rWaWXhkaZR-*P64^lZ83(vx_qquBy zt`1^&+EkPCVI&i~*5Dv{aq@E(B)_ms4PRcVMlAt-DnTa3)YJc*B)v%2KM_Y3@%%Q; zbc=wm=;s?QXrfYqa&<oRA9k?fQ_d8NsH0=L$*gsXPyJKe_lWPYAk9yz_IomZXovA& zO^)6f&3*30OO^_5xFOO}&Bee~;*2S~L9XZ$IJ`5KEE~(L44X9*jrts&amZ=39KW<9 zbiz~ZfJf?!c(}m=evadz+>U<A|55h`a}?m`pWNg4y3R^_yv^32qg2>jAJu7_ckZN$ zvM1TQ_gB$oIC5`Kopx|hG@h@A7k0_5^DNRDEpx3n^zBKsp>tUTTTcynG?wyM+}c*r zZjFYwXNK8c)BtxQgV}EOCT>-yp5}KVHwxtXa=5l!DB%1F^E^Es?Lr>YSniV73WGQO zlcF$sNYaVL7xBamX3+1TNOZ_g@6>B?>lQVjtnzU`I*&h%y;c7ZFJ)rtQ}QKt?A%mE zg5yV#$O={Rs5OJhs0N%vKJG~K+up+l+cDKPFc<qfK9woD;ZF-ZhtdS9qt;c%F$7z* zjcEuL`M71x&w4rQAI?*}Ur6Izca6$}H!Wt`N-KCnrkyOW@xGg$3~o1UZQI6X;QV9` zMpV?9KO3iy(>u#82(b{shSoyA-atXV57{uc;Zv*UZq5=b#KIkr(LubmF-`XdJp1Ad zG~P7j2{4xe$34`C>K~K>2EW4kmA<cljmOh}dlsn5y&a6uV=1jNckVV%9DVlaD`{pd z+d)^R)~ZTfhAdNV(Y;FCS;q|cQqsYcxT8CnfurM%!1Rk5J^^-cwQS+-3}W<mu@h^! z(wa##==72HZ?Kkq3%J7NABq?WFLFLDvzo8Eri7VO`{OV(6Pi1CeVdEOu`v(oAV)ZI z#ST4!$8tBno18^Z8aqg|W5(2$?xy}Rn(1nVNbRsuPvUB%wun|ay0MkSLa2kuqmAak zo2fiPe1jg$#w^+`0n>u*fY?&}IgCYdjW?Y!!TOroK!$a6By72C9oo4fV);<S?w}ym z%`OdXrE99*E^V38deyNut@Y70wXpX%ugWYaEW^@vIjoau4dIeXEUI8}k)E}6X>mn$ z6i!*!gud+K`-g!NTtWu}77A6rnD(kEnnT>uEB2ymQ6tnR`Dp5KJjj=>4$_&)$QV|S z;ZXb`uCH+pT-;Hdu_|jn;GMfPP(tzIT+ZUv)UW!?6RW`=-1YjQoxZ4+k9TLPjxs9i z>d==hq+27GD&Yw!7u%MQEtd405iFLNdB6C(@T?0zt2}S1@WPxzIkrNqMwe3|knowX zN<iz_4z0EvY3%xfg-+MteTefb{ISWzB*~$8u@tuqv5O)*qbYDLR}*<7_xv`ub6Gy# zQ<QSCVpk8UHtoWj(rtE`zR>64=i09fD`?ZbZ||q=JK_>+uB$_+fVKPv(-fp~^^=h& znK7P=$Hs82yPIQ9IgG$YDTgk2Z9n}H^~3;Np;%*y-mN{M8L#lkS${mN&EV%+rQfl| z!HIUlgH&gZ4`c@TXnv{TdNGNngf^x!gPXoWR^N{9VdI8WmUrz84LhQ<($XwbY@;iE zz^M+%HuU4>+7)G@XD3{Bp-<v(tip}EeVC({EB(n&-bRk!xqlW(K*hpdh)&_K)Lr`s z-Kk6Y7*FA%Jf95Owbezuly0#4*EsmUy!4&PW!8)ZM5V^yk8sL0>YSE9c9wg1UQTU^ zW;Jn@46^ay+-%4iwQW0jUQQd2HY0LH@_2!3D4gOvid%MQ*i&VhXUjlma&F}G?l4~L zic06cK#&_qbpR1H>pyWe&Oh<FKhZ!sr2jt726~OJ{ZFHT+<&aW{TFLq|GhQ97G2A0 zED!knk3xxjAW#bk3^)`2IpYGO2Xul!5ab}x{}y%Rqf?yfa=52>ZGQdmMJc@SiBxj% zfj}-Vz~?n&DR1Mb=V)eQ4FW3l|7b-2b4L7+UeZ@l-jnB1?q45lU(2r_-k2{D{5Z3h zf3JY+(|=2UzdQf;_W#HOui*oMQ1QM$fl$8JvA@2Y{~8FDGIMZ5P;j(2(z7zM7q>RH zK^3xb0$yk*vNW@{Ffs(b92|jH+GWu}ed53??Y}fFEb-=@csTQJCm!_F8jLgy$Q#HO zLk(J`6+YMp#c!uZ3Q?ib#wi+FIJ*Slxi#0OrsOw;o7*g-)8Rna;eCsYYp|2QM7i9^ z+k4L^*EJDUXXe}SSh$~+OB;Hv^BLVVEZg)x!{KCPjK*bHEwv=vIV5FKZOvo6dDm;a zxfUa$Q!uTFj}4(?y=ghOr_npa;I}6_+qFMhxQbTR>CwFD0!ho;(R@L|3J3)>Q@Mi9 zhu%Wopbn_I7CJAAFcpz6UDmdp)H6v?lY_Plirf&!BZ&CV*KWUiudsg)UirO-Yi0Bl z3(=1RBr6Rw)6&zKfN<>#UT+G!_Ug-6%q8ZN!M9!64fu62IcSEReWhKj%fg|qD}Uc- zU7Snuq}^soe&+TXg$S;B2XGRcem9u0ia{bYR1mP;lcvJ!^#Wo}<FI(lZ%w)~&T~U+ zl^{ZI{oH;1#L+X&QMLM3QJs&ZgJ4?q;bO@VFAiir`as-h9ceOwlE0qY2dN4?6pKi! zDiK4x)l3<D4@*RX<D6CO4tO=U!<?6jzt7XjRzKAp4lZ6Ts8<H9Y0iE~(`Km&!ALb6 z**XT|#24T<x;V$=E+b(*?ZWg--I)@L2=6&df_GSWDbXeo)9=8}w_kEU^yTrQ=&xmT zfipJBX1G*jo>IRH$=g-g@1z>+h&XGWPuqh^o!VVC3$>oog?($SD#Y@;R|(~GFt+=v z=?_JVZ4V7X=~gTeKdsOalDd8xr9%dxgLi|qFSn~YQa;%5K7>`&p3TLO+we?<oM=^i zbZ-+ciMpy7qUdNyla?vaOV9CeR;r#=;^!@>$ZCc`U121U!Lj1VG4gZEr4y>tE~BOd z0p``kFjhUM9f?isOs9tMRfDx%V;(;u)ynD7Gkm+*-<g2z5C;qycQX|$zw6+6!!9|V z#hNhUrEzOOKGBu6s=&K(i}@(w)s6jD&?HCIY+`cU`vflN6HV@Xn#8fTnF(4;=t`ie zZ7+$tUzZYh#vm$4sID1my@GvE$2K=xj}DHn_Ans7t5ldLt^M6#r4}}pxDR&N3aH>C z@<$^H2o*yK7gA(D6*!W8M(BxN0VJ`NF45(}9+5A_m#8=^T~Db}M*M6SLQD)~9s8%~ zw@6lJ0@GIz)1TDnzAxD&!`$ubDIgpXJ979jJD83qdSc3!>Y{j+7k|Eu%|0#3)-9C7 zGBl4AnKYAA^y}OU)SfJLmD|!ob3dfrj5Aa5`%HU|(Porp_0>$&7Q5eln=<bmM2~0* zk85;x?EPsiq36Z$%=`RmN0!8K*hB^YUr%N<!39{tmA?<yCE=Xi`))4|z7!5@jeABB zV`Jo%6z|EQiZO%lpTddBx+i?M@?^B?!?f^y*tkLLgr*#hfmX4kvaJO_F7s`re{1ba zaeX^%6_uq40Vfdi3lfs2?{o5S){mx6@vp37_n(zqeI8pypes}$Z=;NTThT_!A1o2S zEHJ#e(0wNmfW#mXE!wY)YF9wK(n^rVei`HYyDE?&U%VGjMk>bdAfU@S#m;}Pb$9nD zwf4g{b>^gb`5AH1RPbPMebN~{u-dC(^P(ZXZ81r&QPElG?%0Fs`$Vn?Hy#f$gE0p? zo1y2W-U~O}6(_%iRJhN!$Fr~WL#FCyf<{?W*$*RA**8vl34FV=g_H#sH0VEIOU#9! zUvi;4fdH>L^&;Dz_GwZ&+@&qlZLUsbqT9J}-B4h!jFA?jHvE@^rCuR_ekBIq!hXU? zp7bxr0c=UZb=Zg*UA5HJgMLYRFe+y~A%RjFaH~$zE(X6&ZG8#vVzio7Di}4pN=8|c zgPEz8Xc!yZ_!6(;r#GC*!*QL**^gnPS(o0iRiLrWqCVJiNe*SeBTDuWeoUZv!nfnJ z%Bf-+&6B1mz5ii4rnfA4#6#hnRUz`o&NN}E?wU}|_N4;HiCV*9!LCt%yxQ;*V;}rV z;M<2H%nG>QjTGQuil`paOIQ#oJ#0%30iTv4Inlz*?9{e0wLQW~dviyPlSX;<rFj=< zB0>CdP3{mc+qDN^oLew>S-SytJm!jcQ17>+uv~lfxDFD&xj}8=_>3Sy^ikoNQvE(T z5b$!=$@7ZSL-8U!QnP<<&bSv~ee9i_=@f#_%zQWO#vtiMpitEL?H8Vh3BL?qk3V`g z!5s)^k|I^2S>{c|(iHfntN4S`T)D?Y>W8|<3<1*sh>t!6#irgO1CojAP;XA18`2JC z#uq}u689#fc<taa7%FT-e|hUXG%@~C(K4^rN9zCbb6#Dc_%jOecTGHq!B2?zHG`rr zd#pa>PttgtSQwMC7AG}QkLa-Taaaj@y)oaz5xZlQ1#bj*{JCn#l1}|FQI4?*F?bky zMX=bpn^dv)87P-;RHSn(f4bC6jl@fn%eoWr)t!`@BBfrpFUm~wXOEQWuw3uB?YB*o zhqedSjVL2p+87<;n}LZ@NW|`mz;;2YX&%jw;UBnlZZD!JQ_i=UYHtb&U(9~P{aP~f z8E>u$$9hSQY(N`6Q*5>YicAykk%NdBBOzXj1B1n&w}7bl)JSr;lEv{7#1`GISq<** ziS3h!5cYn}Ta!>MiWelTRyO0@sJ#B1i>7#p^RNn*?nV;Z$jsld-6YTWygOXwnloIL zCDBXZOu<qnN>BBXjw}YNN4PR?-UjDeDrA`7=k#iiH@nDV6{5U}*uiHyg0eOF<yEdf z&8a-@-cJ7Zg4NDGw5`$1u6(I+Kzx6h)B9tutH55^wYN>TUVB=xA;0bVJ(%+D!U&!6 z6g!54;~^x{qIVy5&F#a2BGaKval2?7d(I%nO`Xpu>;T(dhaeBh@qkG~0v^~84)`Ni zCO6eckR;>*6*dNiP;&F&G>Emi@Q8t9B#EuKcSp5+SHhd6rk#p;rlO%{qS^PtQI<8~ zg6qSGB#FYOVY`mxBIH$Z-HwD_4;+eArfj>iun0oji=?7ko1cq9So0(XNK1lP;j~0p zrNu(Jp1aVg#E!k2n$(7(g^L6-(<aGp%gV5C8rWsjoNIh5*+iaqezL=);Z+0q#L4^u z{UPQ)*|?r)k7>}uO5F5v!Khp7CP6XVEIaIPRUhTq9KRZ!U|0PiE|77<=!zz&cG-9y zX`Nj<7v4~phic1?NII9_h3Rys?YA+FxU>!rsssmHvrzZcaNCZ@&csLXWq2KPs?u|H zic(bO^pAvZ9Zx6x%0pM_#K8OAbs*tuQ>}O3U`_8Q<g(<Ul5hPYj?s$$K-o!_GX0I0 z{#$9C%D|$v@;3)F{aPGCs8A<b`2E#_8r+}nu2&LH3HsK)-?P!>7Z>sK-tNRvFb4Yz zU%tBp3*L)xM*O<AMJVPUQbLbf>TZ<(As(+hRPh8;D@#GIy6jnL?>EdFEw|z2dUCt< zyO%ZL;`eh;;2LS2)+KIFDI?V4$hF*3U=-ZY<G*Py*VVkLm&>f3qC?HFgHn|@U{(wA zr;H<(v}~doN>7u_)dP90vJ!7Uj@&5A!c%Hafb;SLI6!qX`%1#^uyS@@RdE|}lPkfF zu%gCUh1oO8(!o9iQ%z9q=K;TK#imwtm&vX!lU?wB%Z-U9H+MaAyQukufqs~@zoD>T zs?B)&qD87h5W-Pob<5H~I2>jJO7(^Jqj1Z^i*g^4_-u`1sbf(i`Yc577|m*O-`+&e zbK@k`=IYDAkv-Aa7)eQP_uCaMxh2_Qp-A0pE|}&J?vTla_vX0cxvjD8S&z>?r=T+w z>iul8Fu1SDQ(>V9`=46T?5?QaMGord4S}VT44$1HxzPEk$@8vy;JCbVEiJP&UgbI* zqoYja4(Z|yvK5Ch33+3(RmJd3NaTgJTawe04fl?1^H_3}Y?Op3aPwGwD54^AWp7mQ zSmmR=4mm^-U6`U_B8}j2Z6qrzVPKRKx#qf}?jR<9EygX2;T2oEH3sfyQ1M3_<hPyK zlaFQD#U&JjQ$1mMwfJ}N8xR@Qiec^1`;?<a!TmE-Wlvrab2jkmm2a+$r%hLnwQVKP z`o&%1V>k~W&4-_2JXKR!^6YdmzA67Q`LcJXuvEeVdG;CEMj=_Dh3-jLJ~~thBIPPD zmGE3W2vc6&nO`1$w*)SZN;2mxyW$+Ngv#c}`>qkAa1=9Yh;K}{LJEf~>!M3XD>5Jx zgf*CF>)96aWsTt>hdF1#->Dhty5yO7wu5pg8cUhU4s_DAU`WUuNnS_sb<GsCyFr{S z(q)PW{uITsD;$Kh6ln&Y+A#bc^Q=5@@ta@i;`I8Vz8k`aF!XWnXHVd;ToV!OJlyLg zu#`DgQ+>!{F5|yR*Wx)+bb)rpPT`F7o3G_?Vo9%ak*0!#8*^m=PX{tbJM9n$0d+g! z4J`^|;Fj+iWT3_{NHey^ldX{CC&VL}OK0_NEY|~weGz6agV%^?EK;EiRnL0zr4_~u zAd;u2(ASW5Ktrji_C)x7_}uPP@1};MeKF0xrO6zVdwfA3j3(gUzt^U@w^hK_sz9)x zjDISHcuB3a=gPa`A%I>1UgMOoIH~gRw<w^$4GWehF8lsPQknhv#)EH@+5LQe7$gz4 zX~&Z(qCYfCKoRl9nf-{MFnC=Z?<+yd*E4-`htdX<s%6{ZuB)2$gyAL|zlr|%pDLn^ zPhyMwR)nxO2%?c_*Zk{;?v|+>5vi@hj_u+|Z`U*F)=fN5Wa4W|^=7~=jJrq@-V~RK z9Qa2SU+Y+v%m&UIIpq`!Z#;=U1E>+iJaT$aab4B>=$hraYF55tGF%msi6PTa;?Fm2 zYtk5sHYp!iIj3mR-nD0#$`sUVOW|DF0_72D-mYjj!xPtrjKOlP#6**r=dBr4a0LO6 zK!5o@M4xTSAKixxN-T0f-Qh>8YfBMv7Uqx%$E1~GDcjB50=McFwfb2YQ=7Vr@#9Pt zXT@%P9817H=6I-k*a!X^omrGSdf2)Hrj){I!<EfQ3^}MOxxhEgZ_vK97=CyD5R^${ zcKK5E_0us}fO|xj0jP#*rUXi;8ZRlD)XF86cJdG*$W>6XUx#w7_=LfqKl>Unjc%p+ zuzt@EB_mzO_$^DH->{;o$IV3`nR$+OHQKN{(gqw&-*F{t|L?=+64eal%QsMq-5m-I z(#2LXmtzK;1AO3nUDn#~i7Kl-`X#{#`5+p&987H?ZXzqeLR40(=0P04{F%rQ*B<*l ztZs_WDJlFbMm_S~rYmNt*tuOajGU!gCr%<UlMow9QjEqE@INn{kN%?L1Zkqh{H!yh zcT0c#Zh>jfZ@uly<)A=Bzw67jOl_!QVZ>21mRA;S620o-p)@qMBaj(qudn^PhbHv7 zaC9?xDyW1IYBRr8LpRpw_M>l0&U+}ow}EqekXAd!{2IrQ5vvZB4quev$fmKE2+h9b z)qb3ci5){NL1BeA{)RS!eNoV(c?{{`yvo2yo6d|;Ii8`q^fRU|tmbV4j7J=d@;9_H z?7TtOA4kx=>c@~44KAvbuA!a<KY52q*WA<&a&B>0tSk4CHXZTPxPzIN9vGU-+1S$! z3kQJoB*d(W)8hq@YHlmqA7uH406_=$aQfzrSD$VDr|hzY%oWMJX+qoqXCPo(X_oGt z`8QyjsXW2s(&Bnb@=z`L{|?x$2~>HXMIH2-i|Ebn`OYGwJ>7&T8F;CVQk5lg>aw(g zsnWW1{eC;zJq`<_<Yn}6wA<Uo-f|yhD1MwqksQyZ)^c8slaK0sd+Ud?A(28eTc(>x z@+uUAO3eD8>;cb_;3cW?BBxCwx)NG=^N|tHqS5Pcb&8;l3het(hPopo$_9Gfz2f{} z6(&XY%Cvz|O;;AkN*dPI&3eYNk-&LiiW(WIw0#PS(wXt}RJwIl8hKJxXM|Z*z}u7} z1f#~^Gag~zD-fWa`|S2&yQJ8%SQ^fr@sXL{n%N(Bp_#*)`Jlf%&U~leO4$}dQf1J} zxxLA0YE<U<+`udou^FwWz>LnhrrV14S;Cchg8Dsos<BDs%Au(V5TJdp>X)%vH5Gn{ zT`J=_5;0q*3nk**D^XizQZ~1fy+oKGeksXc@zD$P2=Ng6S(H;`?U0PxndpK|X6($7 zyMdyzEotY-w~^%1L29jA!Pso}T`Bee+~GmcJKPy3{rXQpz=Gw6s-M?tmA*rZOqsGR zvRZymae-ke++eF=X~FFi{Mwx_=AF)clFltc7v_W&XA7~05*DBDm~}i5$%$@t`q4NB z27z8^#4b+)A)C}^WSHn!v{O&sh&%Q^v#2@b@S7q2Nx>tokWXbX<(-DWu(ALFHI|$7 zE~I>VCSRV%+7o%DU2-&iE+!hP210E>d|=eoOGI?Ct?=&T55x6+_Kjr-vqs+}qO&D3 z$!)|orfj$$lt~v4vg}iPHL$b<zN|U#22w-qL<8-pT}=8S_K$p?#*vAhD~Dbp10qcV z-0_Qs&tK4GD!i*xCEp?3Y>Qb%KSvahY-tD#D0Acbw}oA-QQpm|^NyWH__znVIOAk@ zMn#r^JFF^ZheNZ{Jiklj(VVU9^(3N2IlqUo+5R~YrnHTL!L?Q*jE~qp&3h-?7KJ30 zFJpGHePx~LoqoOM$E|aYi&p+5vDeY^vh>aVWnZU=EnLR0L0)C3OF=D<+f7Tmd;6=j zH+b4s_0MqT)%#KN+C8=fB})|vO6~P_%AQKy^=%%lsOQ`h6{3sF8+rU4CUoBR<@@sY zVGXC$({abm?!r{3QprCP=8A|Ly;1h48K1hf-D!JH*6;?SFxFj+9^F6N%bcE6;t_6q zIFOhK!Cvf6M{~MI(tfwuI2#qE>3KP1kR)=6H7L5t^H>~^gclX%xm3n6?>iVnp(|~1 liyacT;MuYl&#NYyfyL$0^Cfu1U`IAGS?C2q*>ODTe*q%`*n9v0 literal 0 HcmV?d00001 diff --git a/fsd/fsd_v25i.geo.info b/fsd/fsd_v25i.geo.info new file mode 100644 index 00000000..1e0a3511 --- /dev/null +++ b/fsd/fsd_v25i.geo.info @@ -0,0 +1,43 @@ +FSD geometry v25i created with create_fsdgeo_firstRealistic.C + +Position for Au beam at 11 AGeV/c (10 AGeV) and 100% magnetic field strength + +Number of SMALL modules per quadrant: 5 x 3 +Number of MEDIUM modules per quadrant: 4 x 3 +Number of LARGE modules per quadrant: 2 x 2 +FSD thickness: 5 cm +FSD front plane center coordinates: (16.3936, 0, 1000) cm +FSD rotation around y axis: 0.018326 rad +Hole in the wall has a radii of 10 cm + +Parameters of module fsdmodule: +Size: 4 cm x 4 cm +Thickness: 5 cm +Channel Depth (ratio of thickness): 0.75 +Channel Width: 0.2 cm +Channel distance from the edge of scintillator: 0.5 cm +Channel radius at corners: 1 cm + +Parameters of module fsdmodule: +Size: 8 cm x 8 cm +Thickness: 5 cm +Channel Depth (ratio of thickness): 0.75 +Channel Width: 0.2 cm +Channel distance from the edge of scintillator: 0.5 cm +Channel radius at corners: 1 cm + +Parameters of module fsdmodule: +Size: 16 cm x 16 cm +Thickness: 5 cm +Channel Depth (ratio of thickness): 0.75 +Channel Width: 0.2 cm +Channel distance from the edge of scintillator: 0.5 cm +Channel radius at corners: 1 cm + + +FSD size is 192 cm x 160 cm x 5 cm +FSD contains 372 modules in total, +small sized modules: 156 +medium sized modules: 144 +large sized modules: 72 +FSD volume center coordinates: (16.4394, 0, 1002.5) cm diff --git a/fsd/fsd_v25i.geo.root b/fsd/fsd_v25i.geo.root new file mode 100644 index 0000000000000000000000000000000000000000..4d358295af5759aa8473bfdffdd00bcd993309e8 GIT binary patch literal 17480 zcmbWd1ymi)v*^3=Ai;yX2iM^4?jGD-gL|;xPJ#z_4+M8haCZyA-JLr-`M&=<@2>aG zz3a?cv*~H6su`w!U0prS4i2s$&}9M$1TqDI#H2tVydt1137n9?DFF@qn1Mh{)gTb& zXAlHlZ(!8&1lmcS+$01hLi-8O`hRr<fxv<kJ={s_L2MuhpbTgOfnd}mtZdED%w0?k z+*moR7%j{k7@ffu{eM(}V8Ek+a{w4v2%HkY`9}(1GsOSNG+cbJ1sh(#`0pUjKNf6* z$D^xBnmK4V*nY4x6LxVivop5!eD$w6ARW0^OE@5}4Uh-KNG4*&dv(vwv9@%2G*J>1 zX0kw1KCxF+yooC#{b~HtokA3G*evIcD-`L~M#R-E8>AJPSmP@4QrTyL5h#n38n$$5 z_39Ey4WSxzb`{?c;KKuE5ol+zH{Cm2+Sy#pcVv&<FRpG}&X-(nMzdVD{W~r)UozWV z6lc>u&V88%em=X|48v0e`kz`Z*$meo{aU)Oq;p5cyOs97e!Sc}U13i&u7J*#A31EW zSm*T^O@ngpcoiF+YP-<nV_W_rVd-v|ZFRwqchh;i0q4+HS_Uh@M6~j*y3pU)_B6Eo z8ToLBJyFL*p?faS-0qaaWNvU>Y;E!D^Vz4a3g1Yw!hv@2OglWnMGn4dA0Hn>kDp4g zr}I5jA?Kugtr((VjviC^yiv`()1(6Tt}~yP$~`ZIrUEH3K=3Q|{?N3`h)cY+<W?P> z`Lmx@yAknC%thXZNj4tg;qBrR9~0wQ^L)OxI_nb3R?U`t^!%FIHcPNH{WG6jO0sVM zbn}Rr>-ER!hJpymaX0NLG@%b3T2l=Kw4G@6SEeIV#Oc#`Y=W%~Eas!Pvt99@y1ssW z{xngO0vErk$uWfe^-40n6;nqu`U^+X6~(Qo0iB1(UitPFX3iTriMo82r86$<FPQBP z5Y2vu2VyUS$PtY!9{H-ZDYvhd-af;do5<H@>wArKHgxIrjEkSc`<&em2lh|tbGF-w zV&#_S)!!_97i^h75?s%ASiP=Xsojp>IezBmCvKcwtLuLxA(&l%jnUE%TRYBYzL?Y6 zLjYs;v^&|K_c8l+ixpMXd%L;Jf%_PnX&$5B{bp3}u4&c6A&aBd;YP>%dAe-?$-%3? z@?kzyQO|1!%loMOXu;&}qruaHj^+GgHM97z8Ugu~grMym#Ur!);d{-ZxiGdG3;Cz@ z=f#mauWSb2=Qel#8kW7M<BA(TA8l{%lL6wUqskD2d)6!52Gxxg`nwIf1>&r4{jRF- z8JFLI;05DxDStB42rDSfLa7FIvd`j`r(Jq0<q;obnyB2i3j{`g;R_@b{aXKMJ5zsD zPjyGpJ}<|Gn(<Tna`7H+!0=brVEb3DCxpTc#WM3B&Njn0@+n>q%vbwxP$gk-P=civ znKe7t{C6RqE?LzV8Y^r&3~-uP_P4{2)WrOL9Z&tf?xQ~YVJHhbj6_W@bnb$cX&>x9 z@xzP2qv&==*O)cR9|o0BuXP;rHZF=k!_oPl+k3Q9@A=-}!8vqnmJ?qc3x?;nHG+=# zUF@A_@MO3Ym>XD5K2~9OU(9}<Ns7iFxR>k~ZvU7^t}T+{G`(ZUO&<01SYDzRmuxA2 z6U|h?a#EiAdxT+rGz~M0_tLYJ+?(Y^0)$SR{uY<GQ0F?YK;30pwVEF@76oHEf~cBl z=WXM*bdAL>3yNZDPr%{?)b($^<XOH2`d+lbo;plA3qo;}`P2(~`;lLx&f_LTS!T#k z){@{V*M{w(!*pJOq%r;b`RQ*d7p~v7tp*hA+HaKvY))Mpx4C7yMeoYH>hd|p*x^xr z*TE_(+GM3ldC+82UmRBQ5U)lEyth}ek$XEo-T5$E=)UYr7~K)V{_4!Km|kIbsz#Zn zoCPKTN3;X8C4})}*yD%g8xo0$DN3QFNm`+(yx4e<vO*b%B_WFxA~#<H!kFH-9QxHk zH|A->b<4sfL)>&!VQeDkZm$xm*Tbt2dJrWU8I^or8=ZuwrW6U+WPbsj<U9^m=?lvh z0X(Xw6ihhF)fu=70o7z5gF<MF^bJY>G*cqTNbfXmjIu*_Sk$b4T7lu++E6lIK4kCb z``-z#;=|Mea8$EtAkLfe*i55J3`il!P)fwYS-!P=j$;cO6a_3T*^wp_Bbo%2G8YPH zLy8fpKp4yUMgXP;g*}entr!yl4UEzOfBPyXgZ|1KOZ}&iO$Bj2RGugrRkA`e@3X^7 zj7tdEWg!K~RKaBhN;1IE#DK}LQjZ$aCs(%`(zXdOsPLKI4Obf2;`GO+-P1;whDxlc z-@q@4JO&?ajeAsa3~+Fc%3=5hrLV@#3|Qo$z$N}HyWKFX?W%uXey5G_=J*>`S7hsu zHr2QO3@rU+92cL<v#c<OdLH)MLK98{$?N+-fJ!U5fbVXNYTTY~YL85mK5W4V%Df|_ z3fa*?ofUlJ8r3h9-QT&xbj<nl-VK4zcTjTaI+WNVHLo^i-DcJ0^_XJYG+y0Qon6;D z6Z2a>O=2Ds9V9)c@{%h>l-dnWs*tL+k5z|vn%&r;h_IK{nn#+3zfsW|!zymnR@Cbn zHy|pQOZaWcnC#_XiSwg&>2H6?U-Fcfi(_X0a++t&)!c}j_u+Y8A%mH{1vRmNo*SwD zX|v9Ughx_D%yTPraYNWD`fY;bnUw%)6C3B{7wWpVv995p7eEi}Saxc(-|@)GKDj7j zRVk>-M0==x7WYE+O5mX<`JIJI^TwW-+51rO4<PLOq%EN9uafqzI6{vc?oTwaYT-P% z7A>3w;IeBtj)f)#OD@|}l~JBj=NaayGW9Q8_wVBtS$yX%!~3W#qsP7EBhp)vvY*18 z8oa|Y5(ft)fTCL?g~S$lrm{&!o0y2$4hJDw{{lZc*?BxN+ArH1I>S*uGU7d8v|lv1 zU$){rCIVL$q`<y$#t|}L)2&okxYM=GnUD*5{2wEj5>f9SYk0(>>bxn%vcv!l@1n5z za)uKn0WCNzqxy`be+9Z_?nI&8{TVJJ9W^NFcUk=w(~(cyr<d!24<4;eaFWzG!?!&+ z+$UdTdW37~1G3W?+uWceShS|J@fO=6_@>&0M7E@pfSWg<0n6<L-&=f!r<g#RSdxI! zhtw*kBT*iEmE9U>O}@H=EXNJ5W=;EIf7<~e1NXezjr|Wz`<43H(V77Voc%pB2cgLf z^Y5nQfhnx;u;6szcJpexw8f(h3kMQm{5zoy`gijf>?&;Ykwi8*Yh6a1ugMI^ZB8;G zPXmNacX$b@te8--4b~SwQf5zw!s@Wb73=A#4<NGe&r-u0O`-`9b(U>#!~*E%3d)Ih zOM>xD%rq<O(${NE0;~FnUHA<@X;#AV>j>r&dhz$ftB7M)#=?eTw$~0N=_UE{`xsT9 z7d!H6AVWBL&wIJg=zSk`X({^o@%WY?mh?veW4FGY0H;Iu4;i7ew*sIFfLy`L4yiH7 zGITBKW0EuiJ#Zw{=)0_AamvAJ0h<^;pl0Y;HfmJy!8XON6m5u#-`w^2%2a-7X597f z4p<IU&`ICq++Grq)hj@%b^qvZ{<7|tQy)TxnY{uPhYz#Msu=8ZybJ21z5&%k0D|I( zX$f;~mx%)2$pkzq?SVM5mhcnqRp3!aSd$5i<&LHkTn7K*XN}~O5<DIL0r~Iih*1D8 z>6Ak3Kf#FWP=*w%lY=S_W4SUzCrL8{bB!q#<cXS0ux}*i8v-3}z7Q`ySdn+~NuCgk zM*l6Sl%ey-LKB@=XZx1zNuaSRTn?wv+_+hci+05!SB~3dz!RiUO$n3|vD<6CZHRqm z0wrFrPCtm8`DuvVk$Dc+FL2?<58dysAa$@o5%Poq_caE)tSIBGa1h_};12)RH5dXC ziAZGjo6uq}U#BVV1A{Hgpu-cqSjV%hpYbCIlm>FOW&$nKgTrfe!f4z+kl@9kD(f7e z_SGVbEm4Xk%>l1$g}YQZO^xD!#MkOt#OWl-sLHAYRAuu6RF-z4q`!sk)wiU%!?UEu zQNG`T!&!C1!j-S%f`Qql#XEb$k5knRO_X(%C<XOSyGpdTB9kZ{^=?9i%s7Le7I^bZ zqIBSX9okTcZIOdLcupq#ucr+BTl4E5pFV_#N1KHqzVkV6f$^R0;_^-fdzdKTbUg6? zxZ@Q%S+*cgUVlYyyRgWXNuC(cC;@IXoBk11w?_7l1Gh_yhzWq5^yvVdtQ=9`Ja5Mx z@xFMN>w)uMr4Grcf2r_N_>Yzru(4pCJTt(@Q%?-p+VfvI$O6XFvXnB9wrc|-)8t&e zXY>TUcvp+6C-WQbuwu9S;(ZC_WlMVy;&kP%1CUdQ{G#96!D!fjBdb7Fu{{cKrE8Z_ zpz!~N%wXvn;2{(eXKV+T-zXiY!TX9Ke|`yK?+Bfud1Pqq&?8?tqK`1nPclfXY3sy| zPbGhT8Dj4pJo@hhbd~bsfXfziiXW7e`SU9fdxYR3A%%UxAupv=lI*j9x-CKoM`1WL z*1$Eg@*;z9?t2xHh&5-peEvU!DEe0$W#ed4lj3l^Pm)PI#T&Fd)oD_3ye(v>U9#Pt z<fcT751$^;DFh95q@vnhgo4BxgC7{)ePCf^cLc7Lg2ePyz+EKKDb5p+(HqULV=1M! zC@BHK4;23tP7+AsQ+*VSRB3E$T;5ZgoVd%z4a&*@^47?O$(Y`fvRwOBK@aUHMTS_5 zw@weuKOn;=#KvHuo1EN8gL^^^%s)JV6`%uEC{jQbkJR%10&t64l!5^O-F1g$^`hDY zY5tM+90&cWH9xp_{-NjJQL$2SK01wBH5y{cKu^G;yIo5_?=uEGx~H_&I*#{C2{48o z*GiG&rOn<xL|#&p{>1U(C0jat*5tQFWLyMl-!Zb7yKMM7D|8Nq+}C3MjHJ0Z3{IwM zItx#ne$&k+E^xUWfD=Wzh8sHpMZiO%M)P*RS1@4Z%`XM+=An~C0Fz<)k&tU_aB*<5 zZbdXL^IWca{cAeAWo!BLlaC15f($9on4}Q}gvOVe2{>R^iO3`jf7TIkjo&W@-amYJ zv=I))?c37|(|BO^dvXF_zag_X+G<13O^B3<oMC`5t`By)_C+U4v+=tmFqlx*?dE!L zpcG@D15da9A|hfg^5E(33EYE6?2PwlJx@_4+_I7oH@;)cSNia<8f#q@Jr_o6VS1S= zNyF`!I(L2{aQCFzV3%zS?$iuU*(F_R(O9ge4r*e}q+{HID3|fS))fIfVnd1VNX!eh z+lQ%7!I_NJnL-~Dy2L%ixgQp-A65p-ZNO1JJQTVFT_gH#0|d#!H!S`TkJC4bEXJ|V zhJt3+fc&Rj8~#7VoIiY_adRWN&-jZPr|TVnR0~BWc|gqJ676&?ica>!7L1Y?lD|~S zc!~vAz0PHA)Hu(k06K$b;tmqh?&kjSeY=iyJ<C+ymR8_#_`%0#U&#QC!~m@RAnIuh zOrJGl@LLOh)=dAxi0NNWbx{1v1<ijklu9Y4h3&I-F?518Y-?Y!MgH#~|HOZ((h@wl z6FUk1Z=%4*$3~;0{W`3Y0+TPVpX))8@*fx){D(rR08WEqA9wFoas>cxar~>w{||x0 zmUV`zz@<MoM*5FIeriC3;mX<qaE%4rpIRcS4h5_v@e{%_kBYqxJ??QM;4V6W$?_v> zq`;)MimT2zFpDe|D;xLtgsUS}(*QMzfu29ptQ^rs7}x$*_}snDQv5BJyO5U(m_rVB zE-Z8zvhs9PP?L=m0BU5UYMP-YRsSCpS^WhR{Yk)O6!a(Op|!+t`LER>$Mf`>#@=xG z4u7$;*$A+V_g@yhkpBS`z!N?A1i{9c$B}%ncQXS(dI1PA#?An6Wi6594nkFAGeu5- zPWT+EFa+>(Zm=6ur(&^PKezwu`*8pIKC#R{xXDia4~STp7|ai!$ic(lKh@m4{NW5h zB!7D0b8fMh8R&N*-B<)%{i}lU04iumon)!}X#u?KA%Ovw{h#iP{{l=HAmKR__!j#n zTM}z-nEs_7(_ifsN0}tbSoa0o&j<Jb$W&~+sbk(!Cnz4I{VA82uJc7;fotu+W-BCp zrc4~d1kAHyoy;nA{ifDr)km@7oi3M*S~2jWP{Uc5w&IWxrV+h(`ArKi>Nk@{AZosA z7-FG8>)S^blavmpk&q?8O>WSNG6uk8^%p+$$3mCr@MwLzd=W8Ai^eR#|28}Z2;#xa zvV+S_l?j&t+@28-ZdZfTH;XJbRig`d{rii%JTXjJUQ|MI((}9m(cbc8IdiG-N6WWp zeGcD6qhgX#C*G?L<=OWy96wo)^5r%a0+X|MA%z5GxvjN;hCRxfJvP^r8{23WCFwtn z4Z-2tX+XSw$o&^YQFXDYsGjHrd_Vvl1hg&2_YGD!zjYrNE&K)=`Ux#|iH%@AU#PFO z;5mZ9QIU)!=2By_%QoN{{rqx-*HD5zuEYB935*9XZZ{+VD*Miu9F6b=3=){Uyr`Ms zPa@f#;kybkoty1}v*O0!$~H*e5P=I+1{di6(lPi?M=Y-JG+BX0D9U<9^gIB6mPCI* zMQdUPJP5EndN5!<ff?l*@gFeJ0#Ci}1Ay6VTA*2g5FdDle%P8<n6es(h<&L52Q2J= zfCem6=D2(L#}13Yu~UX`D22pdG>s*p8tGNcsPISC8N4MuuhX}%Q$YdNL@aCDKgk7k zz<@HK{DgKl=<q~>#W!`$JSw33PekJz0}fY&)r@9ks4893nPUGkVJx}#pIB)Ih<M8C zbxLai+=&2G7_|5fNB?e)HWJ5zDd8VF0TG|ib8)$RDLr_tG`F2@^LGVQDP{aOZaNY` z8^tC6v(h?L055@P#fEi;dS8o;;*hpYT9?QG5uH7_ULc##P?(s-Uz1h)-{9*w{vRzZ zdnRUoqcA+_D85fW{tUBq=am}80YWch;4ytm07_ICZ?qWH_reS0Z<eT!)tzCjV9x0( zQc7KW0Z_Fi(A+(+xz0kUr=`fkQiP-M<-GU0U5Y!jE-!ym;Cy%#P)}e$NsYcTvI(#X z2(ZMROZH-co~cCK@ieKA6?T#lnCZME6}^5gWXg<hNddQ!+L8Vf3mLn^<zIl~Cje@# zPzm^v8Bmfy=q25o7YIb%GAHw^5LR!l5i_^rK|l-!#3YD+2^iJJltiM;xJSYwONHUS z&J7oH{CvfkxH%yNQ4zVNGG*`usp-<|5SwQtQP!_;uN<nB+24${Nh)XyP7V>E5FR}G z%2OFQu+r<_CgA!T@5oXp|BaFK>HcAjKJ~u^CZKSa(G9G*gY=CL+rSK(6gEzgqr(3i zgl)Io6rfXNi}J`O(hesKnI=LHNlJn>;e+W^5NLAPru}aOC-E-}&A<R676JlChV6nr z;GV$L1P6!)r2pt`2d0336gdEqRQ;@Wnwk-y)fb(Dn2;g=*5t?OL1#MH1h9$~u$A%> zi>k?$DW&WcpOM^G9pX`&ZZvqx3TQKHTA(I0;Ny8kZT^63i9!sM{3Z~wVHkm~rEa)_ zrp&c}8(9o$p6T+-4yPVCf&w5`uO}kL;;ea8ZKg5*mtt;)9pgn560~CSHk42k>gc)X z8nfCA>Kn5%$TELkY3Nc%8%lSXzt9>*&U5<~Jg9dvf%on9wL1x30sycJQ}A44S;_Ha zHYCJO|IbtkELoDPYe#`s6#LhqGC1B6^oprE1?KN|iGO3MZ(xMcf~R5@V1b|v=rww- zG-U+>#HM%L5_8~`J&>8_d_Utq(s_P%7P!G6%4x3ueKKFEk`q3OdgF#o4)SGvBIV69 zOW*v+vF<@76)taxU%$OUr#+q@5}n|bMU=E#m&;Z5LQm#5+1Gk<_gFe8Y9oOo)lV;c zSTETc(*}O91V0;%IG-q`dt+by*-Qa%&ba*9TKTg%18k=NTO<F!n=?Qous8HS+ED?9 zUjNw&Y6rA+fw;gMF^tV%ZOGt_7*Pj1V=FahBYT(E;7u44;PTHE&;i(z0$YjsU#tXg zRQ)$r0;2Jb0r5>>C1O(U97&*;mqYe6*)WLra7EM;xgsx@MPq&hzKE6|MDG=nD`Q{+ zSGN2#tw<1!12zncC18iY6Y@MB$+8C=xTy7BeTBwBdo6FGRzz?lAzNjem7SLM&i4*h zFy_}VK^3vbjpcph=Iwll>MbOfN)kkrr`I=6*e6Kb@A{bI#4;aOH_!Rp=`0I1ovh<` zycC4QOt?<^7cOv&MIZyjT?@oROumoB1s-??niOszHBafQeQ2-|LK9RXW@i4>&Bl-` z`(3xY(s$_8P*YBLve(FO<>R))Mte)5EG&pkS$_mo)>3e#9I>p_G$)Ym!h;?{OpN^4 zh2uo0=L?r-W0s6)R$cy5V2XGKf9z4T&hE%uM#SC>S)^r+ct7MZGVb?w3P+TJ4ST17 zHkxJ4b25y>5{+Tz+B^wM{hN;nXh*$Yh)EQm9dtzG?h#41ec|a*FJkr(N55hCeM5-9 ze8hWFyqb4s5xDWfZxd3zH>Q@KkjJCX>ynQwR)G3W3;8U^%7tc^ZXfI8|Bz>u+D%ja zGhREW;q&iO!cfb$eUE2@ju)FfCX)?fkXuMqLsEc`^tP)aB-UU-R{{|<4i?sUNXVAV zPO!0hmgL9u*_X(3FNOlKceo%A>ELxoR(LX%A6}ZT@oerZaTx-X?A~2{mB(3~L+4b4 z&9a5KHpQBUap+^QiSA3lVCE{?r|;$qMqQxSqMAwT<^+*i?A;xtBs85M_zy$l-PT#$ z?-0D?8raeEm>BKSr<BTj?OzDCMtpr8xjT}Q&@#uf5L`*-PKvi&AFfa3sIIv42^xfE zX9xA&ef11Glg2Vl7IM6t%ZU4QKYq;IVdiW}6Hao<pr4F0{OxO4hA)%r>R%>T;KkOH zgQUEogfFYdq0l*RYeq||Ck$!Zf?zqFIv7>8F-Q3_k>ox(@UXsq7S(MXcYW%tmrsd| za3k0z_4As>6!pf*R9GJ-SgfO5>zkt05?ix;i`h;_)5=IQZ-p89In6}bRRU9eoyqg7 z=6R|9w9j-fk+<xv=^_V-T5PMhL!Y^>@gmxLi4rt!1gpM))0#9CB8uXCfyc?PRe?;e zeD;*e{&VpzXPs&v=IGz2?0LI08JRxPT?(-{{(dgI(_)WbPFN^<g4DiY^8bLe(jXz! zhS{svtXLhwo|5pXLH?focsht5kU(JNSdaqANxeP_{~c~}9oExWF>Cp3Yv@A`38uER z(7Gs6$d|{EsE70fk6#%m>LE7+V;@)@+h{xed7zdvP3`#&U2u9sg+IfUp_2IoC<pP; zHY9|zz5ca&G5_OhEVP;=eehcbybVX=G$7c9aa)SQ5L;`jHk`|3AT>CP(m04l5l$B1 z6DV(s@()4c8yUL&;%GU@f1uXdzo~F7&P3k8@>Rv7kRYd_Mxpo^>0k;^Lkc!43D2KG zet`6*UbBp=3qsSrX0w)2Gl!paIz$~Zb1p0Q^QemwJ#}#j8t$V7`z<BwshYt^Kj}?w z$D!KSu(Hn7M~RB=M5eFej0A3P4<gPMahAu4sL+wDgCO`Ao*j|a37o0>n)q@@9Cx9! z^rHEZTakFbmtsrMn=Y#CRvt9t(}jgWJr2em?msz4Z$N7Je*MIip*&&!sb&PKHn%7P zxslCiIox~INim>x;<(c(v$G_5A?%m7$nkqeEwbnQ0=_-s@TI+1A<0%$c0AfM2I4Nr zOZt14+%*rb9GgrDf#z)8Wl?>&&NIbJ3u!$uE$ZG`lzAOe*u~gs48m=W@=6_aXs>b- z=4~eXaRxD|#<R%lYZR8hXvsy=;oV#=ny}+y*)=rVKBL*(cWF35ze?+!!M}|k%Yr+M zB6jr2!srhv*8YGQ6(?xCE+2N*s`&{EGwR@xMGy(K%(1$lxT6E^YQZ=DT0T_#lDytx z=<Y^gA93aMV_u2>7bHG@ri2T}lzYBHC4Y-0?(4CjjDTR6w?e+>5bU{%4N1<KsCXap zVS3LJuFJ}GiG^$m$J-SUzM>(IKQ^W0-f819X<NS`>2$ANF?IKvm3Ut-7wf;dr@9;L z6rxIh^T0b0{kCIUCGw3o7tv$jb(%~_#4F|PUpQP1PP=jZ0@oHDH<5b|oAX`O;eNv2 zsc$xmq`qhM79-)usF}Os#MEOo6Tid$9na7djM5tJICfus8i@O@XK=5l?}wD6&_ZB` z4ZIm#dup<Ip)ITD<OL`f%0SY4+;#m(Aq^SA%={^>hI34IC&XU!kRji!!P(H8*W^SB zQHV!5L5x!d#Q0-#Kfvu+YwW>nET&qd=~x^QotieSif~0rtsr!>x_`!TDu$NjTIW=( z)Qz5YLuH*41Yc8}Z#r31jEnycMmmIP>4#riNgiBh;NkCxH>g`bZ(U5a2+pvad};Wf zNifYybk?Ty10Q4+l)fDty2Ckbf0AlXhCinqNRBw`kP7WXln;H9!@CaWEimHulm6<& z{izp~(aWsOZwhCIQVS_~GirsLE6nMdFG9<lpsOl0%T!m_=&XH8<@$Uncjf`3RcImJ zv`J{~@%0X$1ddRc=0u?U@6yvy&o1k>%92AVPeLKc@~)aMdZmRf^xN~d@I<v2;cTX7 zH#H@u^D8RevD7Xiv3AF}FU=THF5|w9K4Qz!W43hMkSQOq7SBx2s|g@dUd?i_(kp~B zyTlej*gK8Q8@UEEIjMw=)aZ;QKk1lybC?(F*;ssaf!Nv0vU<DYsq#5htk5>$mk*0F zX>(Yo4_(*|yF$IG&@U|%d_jR@Ce=N`j>Zt?O7VXE`Y>jkB_U=F)lZ$bqp#UTcMJQP zGzMc7>iE+N$6n!UX|o+4IA+uynnD<Sh<fhvcS3&0Z4A#$pZk-cFK#D1yD*sarU(N0 zJGF(|@MlD~8GfA)))(m0HP7WGq95OMK%Di)d4h(IKev4YZMzNI5p+hTiD3_Q`l$8) zGE5T_Jyu&E@%~n!m5|o%Ly`O7==>?DAklN=X85ZWbOUz#q>9GxXnPGxb1o|ZyxsM$ zP`YngRO6r{uG(K=ZL_b8VlWr7m5JMEQOZ%#_K>yEs4Dr%vQ|m@OG{{4RYldUluGPO z_G)klAYkI_U`?Xv{IK^in0ImL_D1#QJ@^%zRG3Tkbc>>x-AjOgH8D31N08gd!%g88 zOB@G{rfxu{akX!)O2e)5rbK|IS)i)hSo^7Ld<F%z`8=V~LVj~y_N66VMY8+>rs19c zn##pJ-51RS5_9j*kP2Q%Su?!5U3K@)RoW}ZByXlzkvJ;4W4+~N3r5neyw}K#G&US& z>Jj`GxEao@KA#+AVJUqthVe@Pa6o=Oi6tJi{T(A2P46syS2Ecd=evuWAycUCfsacZ ziZQ=2)D7||uT!nk1;0yurjVXwKccCY7-bj#kh04fntp<Ggq?$4f#=BJ&WA5PVX(Um z4Fd~Z5jZH_h}D}N^<vGDFEqw+`PSOa>O9lBF?fNZ7O4&jMy{4dXS%zPXN=DZMFin> zb>Y^Qju`HmN!Ku629<}ic1W;m;8^gy#qEA+UuvkFs0!XY5RnKNQxeviDwCqAdMU?~ zu4uTkR{9lmpWV_SafhpJaS^p&j8Zbo?iW3bu2(`jBqxuQr%Nt3iHp~Tp?M6TG9IiW zp>4b;!0PXSn3Qy^#HB?lQ3@nK4xD~JNe`()b6gp=Jij;?y~4BQr%HaK&Kj^EAti1; z#Cj-j;&Fkfd@JVO4WsbJ&1Z<^AkaWkBX;V*yaXO5PUeDY;#3xIWDGH_)XSNkp;d}! zDZT?oV46_&+dSfzn8~PnExC%WMZsCir}mtyG$jYd=vlpE@ghF$H0J@FCKrXAlx5<P zumWvKJUDX}YfZm*pR%%NG%DVPhXr|rb!SX|*b7SxQY*(BEN7a1iko9?jhn1ROwi$) zQ~Ru@jSwB?n(Qv6?zuz+Va{k0C8djtll-+v*LWvXV9{#F4T8($RN&OzW|m`NRYoGK zt|47z?r<BL?QAQDH!x!UobVI(&g3-fDd?-q!koQyI!gZUWyb-x@Z<>LQ;Mrkx;J@8 zrY=!q2PBp!{r!@GDtr@fSj)c)eX)f=3;j+_h^dpet>4?@v&K~(5t_`XLHXMylw&Ka zk<KDxv9hz0#ASfxR6?VuGBokJQaVJa>r1hCK%J!ZSgX6?$GF+d2{-pe?cj>;*^<tM z(0rFi+`1!Cgytb9^KU$xS%?|T&~Y2h&*Q|wZ{+I&km&TkyITH$u?4v-rgcAb5rDKc zRgnunelN=6Nk-T+mFkSFpOB|^Pizy8Fp-MQarG6%^8H>C#i$9}WclNHCeVLGlbdMP zn>ZF^&uQXN@XfuKzl&ok^{Jqz1XbgA7iAg2%Og?TZARAPYu13t1dj~G{CFI<Zv)c5 za`Y662KLCmV@c*pnj8C2B{r!qvCZ3kH(mD<KmB34+$g+j{NuzS=d)k^W==zFuy-i= z0Q>`q<27G{T(<DY3R&|+t|Y9JBK9d;yF0eb!Z&8?j~r}<@iM_`v$?Ok2Gk<wbgVyY zTxsc9S}WhKPY4dU!j6l3uWZs*@dm-3U{vwuB|i}0VolU_54owMTH+lqyE&2Rl%CEs zTUTJ?tjMYMJK`K~5veyCjw#zDmBlPtBw%+Z^tVfd_9XyotCD&@!@!LYKWJqn$)&i> zs(jGUgqF9xC$?fL`xYYPoh@f`v>o=&Dp3ncnV&kHpLD3^!Dv)Ao)h(Dt6NPd^V6ob zvOf!@)U_@T*~^DOH&Ln!+q9^q<r?%-(7RGpQfwn*rBwL^#=h}H%IN3U>6X8Ugo@W* z@xhXmiB#waVg7W2m!5Q^UO*oiHT=w3ulCe8$$yVbjstCdQoYxm`BOhq0B2(4&TLkx zhv2hZSi=p8fmSXSzB+Gw`3-7ir^w-*xor7pR#oJTm1L|^T;?IK*;3Ntw%Cb4jSB%q z0P^7m2jn@Pw?-T0sW7|YPqtVf&Of=w3-Fkc_kNwD%SfZTxjtghI_K6w8|zHAdmp5! z&wS+7mNw<$u4X>>4N=@Zx8A2%Z=~F#@-Uz~#eva%0b)HZ{Lx&_d*RoXrhZEtqB9%Z z)`BjiCnelglOJh|4&9XS2P(4=zK@5iOGP4X@37A^k}xh*u+5b(1?{i|Gv29*V}_@m z*nE^u*<b_hhsIz+zk8!ohhM*-3u9M+&+az%H2PZSZIYaYrGM$igwb<LO&Ok_#S+W3 zpGWMOEk<+^Tng}qn_l-EHaL!IwL-W%-wA0>G7i0)7dVtB(i*X^HjgLz#n70JWK)1& z-t=trb^Y!67rz&ZM2{V_im*+a>DIDJ!SJaMw$}tdEKdfu8h&j#CS>BRv4tV48O)uH zF(n$E<raq9NZ`WiVP0>bq27mYSUU)5f8%e;mMX%*?^iHDzICw7@B<?Il1vQ2bd7Nc z_d?fQ^oN?CG$JOyBm2~UEJKVXG37mrROj9fL>qCGRa?7tS*MIVd-s+$v6OFPYS8Ib zr!7I3Ycv~PrEITb2Y)PWXHD5QoXEs8@Izw#&4!o^Kd@3h|9ToZ?uXQgJwjRSgcWSY za9bX{ZSOq3c*Wad7UGMq@0K{NS3JH%TGIvLv9S?byZGkKMt^m%4sEAGI`Y7cI6}nn zvc8*`K~f(*NO5Ar)|c<1`#F;3VTVlTvhj_~!%TYtqiSU1S1Jdw0k(iXhKpd9#xUs( zW*8TnWS0zFGp-AAbBR(UhuW%O21~O2HJyn9=g4s6Qu#WpTV?dpp@h>xVVb8?I>vJ6 z<Tt1EB^vt`*ShqUM~}3kp5y##tI)_yTgRoy4%StqOIoSe!i5DU&X&c6Wvvke4MPj2 za(42E{!)Bm7ZVO@t-$!U>Ph-T{IV<V;%iAW^e2@#x-kOim(F&I>4}(lPVb=z!a=@I zi7tHnvAhYIt3MIlI(5+^2ol}TlC*WMdaYAxARoMp`d~kN)T$Wk%F-HP(J(Y%Dxc4= zM=et)5Ys4eETvp1?LH$~D7Er?333-$7k}1#{-r4hcM9X$0<{uXL5oBtWWgx|Yv4Gz z(tM=5695r0RZH+TF|g?8CM&BfkJ`mz;u6#jn&R|(k!$7Jm>cEi*SQ@_DoH+)Gz*nG zMlf~h7k)IJGfPZG{ts)bzccM%EqC+Cr<~gpQ);hk!fAoI{07?+q<OW*BJi0diBG`6 zbhWFCXI3MU$U*%pW7z6m#v}TP38ZR??jn<CTXGXY(UY6;SY)fo+G<tcvCYAWe)5A{ z$5(&oOvsUfGSl@EviCCj*c!~9#;V!9+lGgY8*<rxb<_0R$ZqP3Gpq@X9!$Zf1|Y`> z_RZBRniQW7gqkA%lzp6{jk`U#qn9h=iFbZxu0QzKie+FD;4dU6@i^+Q{l)Hd<owMi z@zI`72Ao>!qhHE4ID=|kf?i&F&y=%jM}uS2;t7X&l^YFCOQAa|ynQYwe@SM4=PMoH z62QCJP&8}Z`rvapWj@k`%oiiz3#qGmiuWjO+pg<On{Az=0Gq|T@wI20<zh!tKKBKR z%0#Xm$f()>$+PkP`B3^#HjokJzt6LQ!1=ZRX*Q7mk2Zw=qRsce*9PcfYzF7@fa8CZ zN)!Ttnn4i2FX^9Oe874@CI|#Y1p@tVSw|seM}NW6T6l0X_z;R!eIegn;uZpd++TnL zoU&ALa5ZwZa<B&hmBxQ0V*dF>{*PP=@KC{%=j_E>U0rY)d<Z7=^r<n=bN;;ou226h z{r%hVzqkKqJn$JlkO-CT{gVh40^1J$I3Jt{m9uhjMN)NjHZ!s_bC$L@cR&|&_yBy+ zPGW0iZ)0W({BdywK519P1ocV-pS1seZ(~a^_d$R+|91SrNV~yI*MzEpav{Q`MPBu- zb7<05T8tPiIzytGsg0X^D1m2FUD}s|?@`tcOPGv!P)-DSF^LUMDwk-N8~MBBLdu=v zvGrB~?T<x!*|`i6*9J<MmXSG@_n9sy!=v=>L)vMjQEuU17PMEr$C`G0$C~P}qB}&> zi-ow5+Sk7?<#spvg`4dAU~*jtVMJ-_WuG3+S<aI+zaA+NC9Z^0#jsQ^?0D!Y5)5sJ zsc&ZVrH<5;_}FRh*g-d)3Nta_$gIW>Wj>5d_<ZfT-*bh#KX4TYPS+|JsTE-!iO5!& zW~FCjumb7Y7lNKI-1;joqw$y6PbPUgiW`XQQp&K*+k5IeIG06(omW8t&xUxHRH-{n zvcha_wW`s4b1sl%czvF5qm=_B80a8ixhGFc(Bli_oW|e@*j`(7W}fFp)Tu*7+y;6D z1WIFOS)=RquAsXe$%n$V7$YQ5AzvIQu=_*Z=pDUh1%3W>?i{Km@=zinucb~3_4@sn z0djZ}GCa5J5-%XCxgFxYRNFt#pj>%Ze>kvkHLp_@y83?RZMr^3Z8%n%>F}>(AWwV& zd1HuoOyxct+1)12#M+f5vw%d-TN<{_AxMKUft+y%alZAE`?fb<5Y2csvlEh~Q8ClK zGV7G?O?du}=3WQwKzsCA(_H#4Oxonml2wHLq#^ukdoAFTiXL^e(}9GpPnJK`Y_`00 ziRD{xBm(s!hRN#t=+zII#SY#K)V<uU7|8kKBKi|o)A=-&KyM+k7V%<K3o*V<xg_bV zW{zc~r$}9*#w<H0z+0|<)=ZkSp{1-H0(C}`L5C$sqsFVuE|raIO}UR)5(Qh=l)%|_ zpSGub=Vm=MMXVmE>m2p|8KYgngqa!8#r?(tVw*I0(7cPaMB`1nz$<Rqv24!dVPD-_ z6RPpf>=jkPja%$T8Q(74*P<3*wXDV`#{5nYLf_Hn&ZWy7>swi1G)F9lSUUEQc?EW= z^JfmAgT(5aVAd<S2Mio@bBq`fgz6513p&fh1=8Ez3{>gi;z;}BMlOR2*-_ceWS}%n zW!x!H12qxI_E=!YdqhyAmOCYvin=8}mRzFaEq6Yp$(adrU5K$VQ?~D&V&0<Iorz3c zK~25WX8f_}^cn7M&qx*Nh}4xQkj=$%EX4;~vCI(7x1vPpHX-M<G{>+=8OPK*Mq<KB zSuL<*H$;D;%tQH?5r)?x!)BtDW}p(oIaaG#y4@!$Nk`m1uPvJVH&ES@r2-ytISKct zb;LdwL(}91HLe^fWAG`eLBF4@=)($e#H;oX*JTmhyn1gh4n7w3{~Gg&CB?<cFD=<s zMweoP+&e{(QuIpxVdukQ*Nbfv@UU@%+yP576c4LuOY2w%d0ZaQ!t~lcfcpA&$SyYf zJrsgS{BLMzfnKH0L)kyScSwKY9KBak_wavgk$|n#guab64`{&{u6VFT{y5M4>ca4i zNH7YsOq^t&2D(!r>2eEEI`?IKz<zZIbAfaZfr4Co;6ZSw{THX8-Ikr5qqMrWTXb0y z))i-@#gkzJVc$~En1I<{8&?nm`E|2JMy;m7d{_G}Oz%4?HH5JwsA;TO_?b*2U!5L= zq0U5=UnoU;T)P5$>OU2#*OGP1znA|svs8Tb!6;d1XQqgz@PZ!mCw!^380<?fYzL6w zwWeF(+SNZz%|N(xgt^T%s7moX7q1@-=}|D#W6?+ac(B+b79_0B98lCp93znN@i>?( zHLM;NIkU5lu4W)G)d)`WtUEkJP8VV2gS@-R?^DMB;=6dgCe2Eg_noC9oTy=Jw2SmC z4W2?NS4mSFZd6hDZe!fX@Nt}rZ@4NkxMt8F9QkAiGZB$xdx_bTsh<d)c<sJcvySA; zQ<vTUv>Y{Bl06cjcFV4mc;sdsx7G1TuHpJviT8m{*Ja+R(Ri%J^b%_i@=7G{Z83Hw z!hRz)B!n8exBMaw)R%6qMVH`ri!r<ykycLHzq0hbqbPcEN6b@41dZhd=jmfW!inF# zpkB7>4j_1c!4YKd1Um^>s}aDEZ^hzx^ce9SB<Fd;{KE4eMuF<3C9tI3KRFQbbu-BK zO*BIDB|g%2es0RV7vX&DnV9YngU!l%GvvuE>r13s+>!U2z}rGtL8v<jGl%F7gf~H* zCex(wDtd7ea??ZlL4CHudpzxJePgDGWiS-Ge_@HGpG3cGiVn=HQ@4inLxr*V@W_<i ziC94=giPj2$B5s41`pp^erxJk*BGPpeO#N<Q7utIBmME70BYbJGGXn2<jXFnKh={w z0WS{Lgrdz!t=uCf{9GbVvQba`k0|7>cn#4T(d{6<TFTVZKy0*QTw*K%<{k+gZvO9D zxO>brOE;SGUv1aiYbS@3<f#<Bh=l4-$}CaRuG<zAri62b%MCcLw>|e-$15V*Lh6S# zkZm2z4hgLwq^M;Qb|v6DVYJ^L&5aTscy??pplQ&|wOZ<Liiuy$<l%oRomL{4{f=k9 zs7%?fkC-Jj(*Q&H9^sLPgcK_|Nsb4L!=$H>q~z30cBqQO^%BGt*Y>^!GROzlKLsh0 zeAVwgu~e)eNL;&o+O1Jz{W%vy?GW!_1tP<fETNH2xP9xp3d{4(P_@U`q3Udz9wIkt zj&eyRT6RXtc${wWs{A<zyu38%$Uvox8t+#-sH2sVf+@IRX9l8*wFMQ`9&6UL-gmDj z_I=^?a}FJA-{(}l)Vd&(Uw-ZRx!YOjEbh_MYWPimO0A)w^_m<)V`qMtQDc%D%f<B& z8fC$+7q|BIVP1{(P@$wvGLif10M<>t{|J0P*KWJ00NHWBMME+H#7`c?BM(+jtr(Cj z^Z_j{7PZ*trhzFCXH(G;GtY1;S4q#dR>h8tA4hE)E!%WuL+yA|z=f+KXYvK#+hJKU z)ptWq?VpQLSELQwlY6}JsMA<;oXR7kiSaK|i*Fs)7Q}Go$V^ZcMRB4SNN~zZ#0-6Q zVADuldp6(GnMxKd5GhPqe12VCj(gL<t)T5z8&Jh1@w~mp4VO+(1MDYG6cic{vh^w^ zcE@>7fgY9<r<MvwJX<!2O1NgY;rF%tHD<mB)*6I)7!UG+%o|3Q--GIw%;!+nx#e>a zO$`NTe>suJ=L);Co^H49o6}3n8wjAw@NhMWbxjVnZhLQ!vqLT+8d%eoo$FASp|fSM z6Tfyn9S^LCSYecc?DH~!MyyM--+6^IwU?aBk&jOG>o;k<Uea5d4$3c6d4f!NW%Zi< z3-%g$E>^~Mc*HOf9~cn#RtjtJ*WO$&C!Z4auKu{^Vk{^r78bnSPNZfD3lhJ4a|sc) z8|{YtY4sPeR8V*+6MC7KS;5;Rf{F;W6KuU~Rim2nXZ78ExL10fLrdSNoYwDNR>e!m zXP+Q-(>v@-J)gb|(@CS&@ykI_^TUqqGhD7~`_?R#+kc3Qu)+;Z)7XGpDJ+;Yk6zSs zh+{50O|{ku5wy!rxn&=|(NILBc|Q&*C=B2L_zdk6nXt?9*?D#4ZS+m93^&rUHg7d{ z_XtNj_aIyiQB9yX;*K4ccJW;nw~j(iA^9&)R{GrB^{lPprW0nSA&S0+qQc2m^R0_! zxpq+~SKXCgwkG0Ha2qgMFM{mi%?~dcy(H2zwXS8Z#f_LVP+_C=E1!FJ$Ge{!Ctx;L zUJj0&Nk&J>N^`qjFY75UDh`Ro7+&+iHHGtsPt22B<B#RGBzR>%KKq}7&d})gawsDa zK7F2yj6m9Z*Mi}6MfWCVz`$q_B7<z;?DWW;F;H7YaK#(X{f$RixvlvM-{B}DO&WiA zCvT{uG@M2FD~n&%%+JImzBoIjU%PV<-f(Rm%Z^ZvkdcIJ9_tK7SH>*wj))#>vO614 zK@~Gbs+p$Hiyqg-aB>od#D1W9zpiFDfK6D3b<1IT#nonyg|7rEVRt}%-H|i#Se{c- zN<A>y9hqN8c!#(Fm06<}*%r4)Gg2JZH%(jq<QqNffT&aT>dJh|a^+axQ3j(=+C3?r z_W;^@=qcVuD~%)H$q*|~<G01f-8<FAQV!@dB~%C1&#KLgPlhUS5$aH1u0qm?&vinv zRdn2hRS<Vd5fW)-zn<k(o+FpiI{YN>95#zWv!a8_W4#qqJzQRwTs&G<09hcd!adv1 zG*d0<4h_1@x{2<mWo8&sWf3?ID5L2vW_@;HlxKiILFGyHJxZ!?qGs3$<!zQPS3?S- zE>T?OA!eXXH}TPj6ZW3t6hKJY{KS}O)Ccp!6fvB+mwzv39FOCgglPNW-YA))+_i@G zZ8lrE@J)uEz>%6etQ+nZ-k89-I-U<48TIb+w9p8n9vqMvz|PTj8`ME?{dQDCv+5|M z?VAP#m{BZ>%wJ=l?NC(4rDNF2W{hrZ*Mms|kY+B!R!Qh>(qK%r&bkZaRYy&rK2J?z zuA*#%2Gi1<NeFuh`Mqf0Ob*2bU|Z#-D;!gKf5aS!BNE=b*QdXCRK?Y+L~@=;dMbl@ zNvm?^%fAsIf?ZbVUzM>rsrC-CDP+2h3{xR3|M5{)gZuf$TWFKb>wInqBon#mB#<Rx zJUBx{9sSXb`-rG0Y+Z-o6VaDXXU0@6Wepb9OO8XGSGDWOL*E?&$NQ4jG$mP{q!xtj zh~aOLBx5kHh1U<gY}0t6(^|w`+oVxmuV*o?TlkzPB-NG~O+(t4cakN)Dk+x$z5y<| zHn1z537Iqd@U=vI<4N)vK#eH&(T4|3j};?!k8F<>tBPfdq3ZA~EQN;BAfYKoi^d3y z36+qlSv8xs&Rx?q*3ceDYPYgx81ESCHZ`Ydfy7o+ERJh+R{E5DzxT0)S5SyZOqb;0 z#$1y@nEsS7QZfBHE<am5nu|%Zu?H=9ChT0xxNc_W`L(X-bj~7KTXo#cAE&E%D|fyn zas=;Tk41Pzz7?)Dm_fT^g0DYd{Zce#y1W^Kr3_Q89P+B^6~@PA(;secL$m0uE?=rY zy*q{o_KNN_0oBq@m%@nE5TwRY*tsV#OdKMGdWg#Q8PKejoG=Fo=UgMFGcGqBe%lR1 z%git^f6dVwIHYFjeRB~)Y5i5d24lzzWdnh}_qd9)Z~ySQR4Y^C@)gWNSG#J1e2JaH z<){g7zYyeZr@cNoNmY$^pDYBi5L5%7i=`vfO-vO;xaLar9Ej(45Gy6}>SLd`-A&0k z4RuiEh<AbKROJjUH@~~CnVWpe_(=?QDsn^V7qhWsLZ$ihk>8BGp!XQ@YxP!4o*9on zY_Oe$?Y9E>TvSPzb^`d8=u9<j%y^1N^ULEb;#R!9)dxqnMY0l|jrI4t=_8(tMmEDH zLraNaHVayG4HL|6*#nxtlEVbP4w>DBw%axr);)%fUU8{%`KW<FIfc7OY?YT^$37XK zFp6G^#))X2hcS$MQP}<d7}~{cg_)NjgAJ={EK_T7Extao_H_fCcOsld9!5EC{(#5N zBiJ6DV`!TOcP$!^2%o|=!6AxOPwj)Rw|E@(ReLC#u7v6QVXTV}%uN+s+!>}t{lNAl z)Qskb#|vPqxwUw2fFlnHiV@-A^wlfhUdM0ma?0nkmSyv&i1GW~fP`(8Rfb>I--K<J z#<+lcv&SjfLyheJJ7K#jQtfvZd(dMoVKlQFFoTr-bQ7*_;;S`6Tb{ye$k76!#o*q# z|9Yfr3?5G1*NlCn%g^1}b`NbZX^dWtiom_jc21dBh?cyq<!$+(M3I#v>rD(*HJV8k z_P5ZSexKp6MY*x!51VF;r3{GH!^1wsBiB(n)S>N_xc3pv^+#qj4NQi+B?V!ctZLj< z>HQ<`JvgAN=s8<9zp<1Lhs*&kYG$U{`YtrqVA{t={nkT$_(@BL1#U%^U{j75f{t+4 ze3*N;P=sMt$@9f=QLTBQEQ&jmosG$!EeL<UiN~Jppsym)db`g~!x2hWbHL86?fZwx z5ryM(6RRxbCXDVvD@M24E<1*284tE`I&%IrbBn6wLrV)FK})U`n7L9t8Fh$Trr<Li zJyUK7BjMH~Q&(+KKD(W>NSrKvDJxvb?hAT^dPq=`<ds-Gq@;5rx!_V5J#*!6pss37 z-98FvB)fExTkTRcx0-oVhTD&Dco6yqfBJ**w|79og5#%Fpzlhxu}iZ;xrQUEc0qSZ zp=lYyK#OTv;q4RR>RkZ#oxy#o!7WlJ_P8ByGpVjRj!?j;eG-t#iEH`rvvCvx3bV+J zTZIftF}2ssG{v=e`x{jw{wR4?aZ~tVo+;r;;Um77e^m+1ovz4`h6oWIj;H(%v`R*n zP`>x-6IGSd=QyTZYz%Z=q`Kgwkl3r2=(x|0;yV*RP1pChH<qBR8Us>E&K4<Uw~$*| za}a*gq+UEIa!>Bo!ZQ#BaDIK$pBCXH8RA6eZqXaPcjW&xhD!2WHTV(}9AgpeMOZv^ z{(`Ac=~t5``v&P|OUf?pIl7SSm##>^20wmKYvjc$&E2ey;OJ?zzgL*M8(vOFY)m<% z%Zge~6f7tG^P4n*_cK*JJ|qli=l5_9TWkH1>RVV?e5<A6gvf1Ef_I9ou_$r{3RWju zSN2(c8P}_U{08Uv7!^-4yY0;{i+Rp3dj`c^Q3`<#Dw>0xs@nPdo_hLSTc70pAk%-< ztVOY{+>cn-?Q+el+iJ?t=&!fYbXOU!ZwctdKIfijl3ZNgs1R<mVhXk`-B*xDHk{H; zB_20<iPN6SeO^nREhcUBL))cedFs;lV(30uB^ZdsT6Z^l^ipzGI6bK%Al`U;ATu70 zyU>+^@!=jt|IKFOOl<6XpUXj$REbNR0m)5)$CBVwg4kG}#d4mxfPr{wLwTE9-0;MC gpXS{pL2b!Q96tB%k72_m+ln!t#a^H^T*qSn7gn9P2mk;8 literal 0 HcmV?d00001 diff --git a/macro/fsd/fair/create_fsdgeo_firstRealistic.C b/macro/fsd/fair/create_fsdgeo_firstRealistic.C index 2df8de2f..56790a4f 100644 --- a/macro/fsd/fair/create_fsdgeo_firstRealistic.C +++ b/macro/fsd/fair/create_fsdgeo_firstRealistic.C @@ -1,6 +1,6 @@ -/* Copyright (C) 2023 Physikalisches Institut, Eberhard Karls Universitaet Tuebingen, Tuebingen +/* Copyright (C) 2023-2025 Physikalisches Institut, Eberhard Karls Universitaet Tuebingen, Tuebingen SPDX-License-Identifier: GPL-3.0-only - Authors: Volker Friese, Lukas Chlad [committer] */ + Authors: Radim Dvorak, Volker Friese, Lukas Chlad [committer] */ /** @file create_fsdgeo_firstRealistic.C ** @author Lukas Chlad <l.chlad@gsi.de> @@ -8,8 +8,8 @@ ** @version 1.0 ** ** This macro creates a FSD geometry in ROOT format to be used as input - ** for the transport simulation. The labeling of versions is ment - ** to correspond with beampipe versions v21[e-h]. Letter 'z' is reserved + ** for the transport simulation. The labeling of versions is ment + ** to correspond with beampipe versions v24[e-h]. Letter 'z' is reserved ** for parking position. ** ** It allows to create module stacks of varying sizes. @@ -32,7 +32,7 @@ std::pair<Bool_t, std::pair<Double_t,Double_t> > calc_x_trans_and_y_rot(Double_t // ====== Main function ===== // ============================================================================ -void create_fsdgeo_firstRealistic(TString geoTag = "v23i") +void create_fsdgeo_firstRealistic(TString geoTag = "v25i") { // ----- Steering variables --------------------------------------------- Double_t fsdX; // x position (cm) of FSD in cave (front plane center) @@ -49,22 +49,25 @@ void create_fsdgeo_firstRealistic(TString geoTag = "v23i") fsdX = calcPosxRoty.second.first; fsdRotY = calcPosxRoty.second.second; - if (geoTag == "v23e") { + if (geoTag == "v25e") { comment = "Position for maximum deflection angle, e.g. E_kin_beam = 5A GeV and 100% magnetic field strength"; } - else if (geoTag == "v23f") { + else if (geoTag == "v25f") { comment = "Position for NO deflection angle"; } - else if (geoTag == "v23g") { - comment = "Position for Au beam at 3.3A GeV/c and 50% magnetic field strength"; + else if (geoTag == "v25g") { + comment = "to be add later"; } - else if (geoTag == "v23h") { + else if (geoTag == "v25h") { comment = "Position for Au beam at 12A GeV/c and 100% magnetic field strength"; } - else if (geoTag == "v23i") { + else if (geoTag == "v25i") { comment = "Position for Au beam at 11 AGeV/c (10 AGeV) and 100% magnetic field strength"; } - else if (geoTag == "v23z") { + else if (geoTag == "v25i") { + comment = "Position for Au beam at 11 AGeV/c (10 AGeV) and 100% magnetic field strength"; + } + else if (geoTag == "v25z") { fsdX = -90; fsdY = -70; fsdZ = 1756; @@ -212,7 +215,7 @@ void create_fsdgeo_firstRealistic(TString geoTag = "v23i") const Double_t plasticSizeZ = 0.5*wallThickness; // half-thickness const TString plasticUnitName = "unit_Plastic"; TGeoVolume* plasticUnit = new TGeoVolumeAssembly(plasticUnitName); - + std::cout << "Module array in one quadrant is " << std::ceil(holeRadius/smallCellSize) << "x" << std::ceil(holeRadius/smallCellSize) << " small module size is taken by hole in FSD," << nSmallCells_X << "x" << nSmallCells_Y << " small modules, " @@ -500,7 +503,7 @@ TGeoVolume* ConstructGeneralModule(const char* name, Double_t sizeXY, Double_t c shiftFiberDownRight->RegisterYourself(); TGeoTranslation *shiftFiberDownLeft = new TGeoTranslation(Form("shiftFiberDownLeft%s",suffix.Data()) ,-0.5*sizeXY+edgeMargin+bendRadius,-0.5*sizeXY+edgeMargin+bendRadius,0.5*cellThickness*(1.-fiberDepth)); shiftFiberDownLeft->RegisterYourself(); - + //fiber parts TGeoBBox *fiberBoxHorizontal = new TGeoBBox(Form("fiberBoxHorizontal%s",suffix.Data()), 0.5*sizeXY-edgeMargin-bendRadius, 0.5*fiberWidth, 0.5*fiberDepth*cellThickness); @@ -563,7 +566,7 @@ TGeoVolume* ConstructGeneralModule(const char* name, Double_t sizeXY, Double_t c std::pair<Bool_t, std::pair<Double_t,Double_t> > calc_x_trans_and_y_rot(Double_t zpos, TString geoTag) { const TString strPsdTube = "psd_tube"; // node of beampipe containing this string is what we need - TString strPipeFile = "${VMCWORKDIR}/geometry/pipe/pipe_v21" + geoTag(3,geoTag.Length()) + ".geo.root"; + TString strPipeFile = "${VMCWORKDIR}/geometry/pipe/pipe_v24" + geoTag(3,geoTag.Length()) + ".geo.root"; if(gSystem->ExpandPathName(strPipeFile)){ std::cout << __func__ << ": Error while expanding path to beampipe file!!" << std::endl; return std::make_pair(kFALSE,std::make_pair(0.,0.)); @@ -632,7 +635,7 @@ std::pair<Bool_t, std::pair<Double_t,Double_t> > calc_x_trans_and_y_rot(Double_t // | . // | / . /<-. // | / . / `. angle of rotation around y axis - // | / . / ; + // | / . / ; // ^ value of x_lt / x /-------- "this is 0 but in principle can be different" // | / . / tan(phi_y) = (x_FSD - x_lt "-x_glob") / (z_FSD - z_lt - z_glob) // | / . / from here we get: -- GitLab