From 803dc399b238096dcf58224967bf3ba25a3efcd4 Mon Sep 17 00:00:00 2001
From: abercuci <abercuci@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Fri, 4 Nov 2011 08:42:47 +0000
Subject: [PATCH] update TRD code for PbPb 2011 - use new alignment/roads in
 tracking (Johannes/Dariusz/Alex) - use new ExB/v drift calibration in
 reconstruction (Raphaelle/Theo/Alex) - use chamber status in tracking and
 save it in ESD track (Raphaelle/Alex) - use same/opposite sign track
 selection in AttachClusters (Alex) - remove updates of tracklet reference
 after track update by TRD (Alex)

---
 .../Align/Data/Run151536_999999999_v3_s0.root | Bin 0 -> 17877 bytes
 .../TrkAttach/Run151536_999999999_v0_s1.root  | Bin 0 -> 18039 bytes
 TRD/AliTRDseedV1.cxx                          |  11 +++++---
 TRD/AliTRDseedV1.h                            |   3 +++
 TRD/AliTRDtrackV1.cxx                         |   5 +++-
 TRD/AliTRDtrackerV1.cxx                       |  11 ++++++--
 TRD/AliTRDtransform.cxx                       |  16 +++++++++--
 TRD/AliTRDtransform.h                         |   2 ++
 TRD/Cal/AliTRDCalTrkAttach.cxx                |  25 +++++++++---------
 TRD/Cal/AliTRDCalTrkAttach.h                  |   5 +++-
 TRD/Cal/AliTRDmakeTrkDB.C                     |   6 ++---
 11 files changed, 60 insertions(+), 24 deletions(-)
 create mode 100644 OCDB/TRD/Align/Data/Run151536_999999999_v3_s0.root
 create mode 100644 OCDB/TRD/Calib/TrkAttach/Run151536_999999999_v0_s1.root

diff --git a/OCDB/TRD/Align/Data/Run151536_999999999_v3_s0.root b/OCDB/TRD/Align/Data/Run151536_999999999_v3_s0.root
new file mode 100644
index 0000000000000000000000000000000000000000..e0079e64a7d86046ddd8dd74d470e37d166f4341
GIT binary patch
literal 17877
zcmd42bx>VRw<eB5@BqOrNN`DT=imen?(T9RxH|+VxF6g#L4qCJ9fG?%2e*Sv-tXRb
z=AC<|YQFjX_1jgwS9kB#>v`7dwX3?iYJCL)U0`6wMPXoI%wS+7Xa6+6KYh(V4fap-
zLimr45C#kk@+l0g0L_IkQuEu*RDuWE@jx)-&)0wJX@vR9khR!g)Etb$zbyYjgMmR*
zld`rmXHxqj0kE^Sbg(yfaAEo4>hO{4BNqp^A>ZGRp&N&xGdt_ozqtNCDZ?QCh4?R2
z7?_YhfL?#vzmQ>IBp&~7?g(7{#kzO;U+X#kVhv0D*Lt)+WW^=KBpqD7dj7Kyqa%;f
zBmgtr3*!NMn&By^ZV0_=k}s1eQjwMXWTr@*N<}7~^_eu8f7+4QkxCr5B#=hp?R#+|
zxla?7rsU#cVPY1aq$DFPPodzr^>(OmyI;*dXurl8yl!y&vS8Uh|DreFGS}c<W{j9j
zi2R7<%R5EE#q4i_<ct2BO<O0cZq3x(-0K~I06XiPWYEHvQ0jQW>xA|Js6~{!;S*O=
zqroMXbvRs@Q|Z?B5AdLHH<?9i*KI&H@@qBsNkgXr+P_}+jf3OxjmRg0Mv!lW-18gp
zVOb5$ZGy+5fv&lO{K_U|^gl5hYcbWdkae|@%PWKJU-^`cVk$d{wF9|+z}cE1&JChg
zet)M$&1Cn^Y)pClF|v*HIpGNu)C~&S0tHQjf-XVbDX;i=FNGU$fk9{psI1uZzeIEY
z4D!DrbhrPv$q)oML;@5%JxXQAK+!J5e+|rO{$BI90d_f+dyYqD<CJb33x>bNAANW(
zZ?^k)|6%`c-9KP~RR6H{pQhm@|03}Zr9T7Y$iE@~Qz=w?cuu5X;vXP?L89OZD*wC-
zG7S8uJAYXsr~mKd1(pBd^MAs5<o*{A_0oz;Ux<8D*yW7h{KWfP{KY(@o!|-ee~ZSS
z!MW(${|w?k`2Sl@b~BqVgQhJ10rD5*ALgZt-d#1F{}J#11&4+AUn2BxJ6TJQdg(9W
zuVJKOfrxMZxA=>C-Px|8=KnqPf2jO>=oSo>NcR6*BK*Ba|4!il-C!E0{Memn__s5`
z;r)LP+<%|@iI*(0{m&Hs59TY~|HsWR6#f6|(4RdP;vV~=)y*7AJNiGx|0&`B?Wq6t
z5Byg;c3tcEuRsbMseXFGcLl~S(Y=7ZTt{U-D^eHCY0^v5VbUYgCDKPy9J=(xZrR;J
zr)fS}Ju_BRg5;=v$$cuP#Ab#SH3wpz?95*Ay}WZ}Pb)zeavZkw)NZ-mqKj!yS?3{V
z2In$oYv)5}SsycDREp%de(8N`drDxgU9Ll}eJ+qP;(H&xiQ2~uLj7nOMf{R)s(sd5
zE>6=@mr}!0M^Z~tk5UN;+Cfx+<XqWFTA9qjO!~~yOsmX;Oqm2T*+?4ogb%W*G_v1|
z=)b5jWZ>yX(kP@9%cvGvPrFzrplgRwi6`d`O`6MO4rMZAmStLJ9%jn^Fw>4?R!{jL
zmr6VIeV6{L8eayfek_erO39$=uJxshX9ApdV7{n&F2m$InMgCqNGkR3#WI>S4BwOK
zKdaGZ;OIvb%4-zUt0r3yyO>YoU*ZqrAK@?IKjJ5#YKP{FsppkVs>|pO=`!e+>00X^
z>dKazsYf!ZCl`O#q_w987T6Uy6xbI4sqN{2C3YnaCH5six>HI>?n&-h?rAQB@-W_&
zo?neD15-b|P;Rd1;EStui;I2&xprK>^lZW5)Q!xenQx@9I;dD_p6>gQ8EX3H;_ruM
z(|}9BFyIKV1b75YzrTDx{Ql_u()-8v)5e#^!^TI(OU94J39yR+CZb%qM3WdY(u30U
z(xuW?(g)HqVP@iy?`PvJWO!(ZzQ@o@sZnRVUHvBUUX>n?F;pr~Ri@v}H1gwYnuR<M
z>#yP71^)%|TIYal_!Cm2eHa1c1-{O4SU1f7x2YGPzxu?B&@yZzydClaW-XbKK&M;4
zDa`*bCN^GM;cK_oG{f1H*yH15A6B%YyW%WTf|kX}JLl49Q`_M`YEjhD&BupyUGxH+
zsc%5760bsicAi^m+Wacd(8cxY*G$Y01RxQ{+PA^fvr5phDE5UK!FRrDv?xznwCY5^
zqq_`g!M=>JYWF|KO*D0~1PPtwT~re~jl73b6RX41YMqrjxi6~ZX4wN9^;^BBZjyPB
zHb{}~Jwj+U4|)+?Tlw5dbc}Ns3e*oT+0b`b3=g{5bhB9=@4W9jo<Df+cg%`B!>Zd>
z?<x$N8aev@YFWWY_j{bPvivbb=93~^Q*ICz_~uKoOymdnSf7ZwjPH0spz#N*YN6P!
z?&X_f?)IJyZ?aS-GMnn$3Ij0UoEoXekSB>y$tc&z@w1#SDSVag*=SXCcW39c+9}$!
zO|CuD{B!m(`{V{;C1+AdV55Q|o&s@tuJ^OgUbPK+eW?@&nuckB_17mDS5_piBG)ZW
zf?UY@O6|Fp2vXQ|ED=4kSx`uB`|SXNC6LSe_Fx@3*ES7euLg_sGdRS48#C`|vF>>D
z$pYPl6&GhEDKPLy>$^aMN5WA}Ifpc+%YmlXeZ1>AS_OUH$S$18dul5Zq+`D?$Gu7H
zlR{@x36!B-aBmvaMfBJSv4yF<Rj}cOkcTeBG6)aWn{(W}v2fiE6b8NDQ+Q{5_P#e`
z=k<xdd>&K(@NHs&x?_n4*~LYl9m3SLvFuA0v^Vlm`R=yLHei}<@^xll0~;fF9lv_C
zHt$;y^CZxU4>DXWa+%4c%C4-vXtEqJ2j4H^a3zMJk&R<OaCggLwew`OxEBF{P``iI
z#fk3nfazqQHUZNkwlmw}g6^irT{UG8_ViAPu!esE#uq=*mQ7F8Z!tSf(_kFN6*+X4
z+$k5T@`gb$H=?U8t=Qf*g<f&>L!GcN;XTgj*O#>qNyu&qQE_^`qBAY3$aTW#(dWKA
zYF^P@HUk)z>un|hZst{KiuE~A!C7MW^V}N&?(ya@0(f#I8|pqnD^kVm<7E|cZUUlT
z*SeoXXQ^1)-=o9#EP;gFA50w`zjtc4<DnNvV?k$&>KvdHBCkhE4ZDj;gzCK*DxBUg
z#Jg9hQ_aW9eq&}nZv_eh0#}fK#Nk&19&&UuxkE4QJVrV*+v%0kkQ!H|E(v`a@|2`2
zF#en{TIW8l3?O#k&%>dI^_|qNz2_H*6Pc&A8W(umin#kA(<~l7$$?awGT~5==d|{q
zaUm!59u#%s0nERUQm>RmKEd|yc|5$ZMt{O_BweK+3YZD-4qp{QuV=Gi#uarsQWcz7
z-2T>;!wf4B<xG#I0M_H6GR^Iw)|>L;LTc?HSsqJa?~Ld9MkeuaFYLUd6-|N!z17q1
zA{!MIE~y~iH_lyAI%&GnM|Z;y=)oC$-OO|8GfEPn9e!i6N^gjYGw!1fIIhmzoVgvu
zuS(r5xzn$a>cCoCzkK10)KtzkPr<(YDSO&`y!{OqkHXWjg)Dl&98JV1U$9AQWosbD
zCq^NW@m<FoDJ%*<Fh*6X+sg|N>jj~iaOASkCmS@RG9E)zx;(*;=w}-ZouYf!_o+7W
z2tGn8g0~N-NFs5PsMa%6_s30W7tGIy{g~42zwR-DaZ4K7M*RwWsys#9zG=Yodw719
zjIRlvn;&x{GBH|tICHHVkw53_+O;2YKv!dBb&O;06dgvMzc?gxG31XzEzwg3Uk)I<
z{W96OD!pAg9lYyWxVz#LWRJ&R7*`3stW%KfuzogHbMK?D9V^I(dTRU0+9AM)Ho6To
zDC>nM`Zpv_@mGO&8z8AHNJ6T!wk(r^Ef~#$pvrF1wj~#A++C%?dDb33P6N$>sj%)<
z!t2(;jv0B3V0&4&MAFlPmrCwNNY8QJPN~~4O=OhjubEzRyt_V4?n=pZ;M&(SDQTi5
z%j=#JG>}oP%h<Oz#K@WO9|IJ4%A|CL4E0)(rW>unYBOyx=zHU7D_dgQ4WDf;TO;#t
z$~tb&*e@@j&ByzvrRZta2^SSB<hxP8pyCMfawW@cDP)}~-!pc7AHrRgWnUIsSrtwY
zJieNfS}Ok@a5g2QJNLP&Tk|bpdBKtZykVXo@huA|@w3!Rg#0k3zL={S@aE7#?G?<^
zRO7l5+VxRx7;`9`p1fcPb!Y5Au$pg01phQq&IZE;lJDf^d}kxlxl!sA|C>P_sgkiK
z1k>u3P#=(E=QAA8z&dY|lHZc94gu+?#PmN;?sqdAb@r#Qjkx84bH@xkLs7QeK-I7H
z25*BEjslZic%s#j-44fD$Nb4UZ8Cautq_R~7w)|%GNg^Op3Ja@e}48$|ItJ*9Gq>~
z_831=kmFsPXv$%9qjHuM=XCIBc7=Oo_-P~F@GfZcP@SEtT6h{}w?_-|n8%}vzQCmH
zN+!kK8TR#t>`vv$oAs>ummYd>yPOT2()wk%r!-sGCDl6kL#v!Rg;6Ig)9F=|=nHFu
z{{&>hb9@|hJ@C{arMdDsqPsI(svT%U4{Q-+V?FCs_GON1(sqM|kFhBh8!E-}Ov>Gx
z?v`-v)(Zr(>0pV^V~o1z`qZ#5)#$X|LtffuuveZiA6^^k=@YPKp<&`9x0T=fjJmN-
z<K`VH___3s<tLv;fDHWnrnR<BXM8FjkigjwURNlQK={<X#GO~-*NKPyl!m!U#Hpu(
z)|Eb$9N(slX`Gr~7;2l1d2J*ei0f*X63AGdnR>0sn->0M@OU#)LE=W$XuCTxlhBA}
zk&79KQ7)h>*lUh&WHjt6X-c3&N0A9W=&vwdd|jk?J(=GS)r?`~SrQRuIk9_W8G<1T
z2lJ8{gdSz>m0howzTl&`nF~G!T-x%wE77#nG$xKSH1U|^?$9<!oZ}5cy5D|<%^|4Z
zMq0=dt~B4h*bDSUU+@9Qu|he^Z5uT<^L~`Fu9<kX$y;{E#VnFFTCy)#y*K_yhj0o4
zTIn|By_87t78EzbJm5nWhlXKyb>pReYxC7(vV0nx>Gucu@9VzieN(U^^Xv^PB~&r?
zH(^@L2zv7Bs=AaykC`c_@2UAVHv!Aq+|n!~;`Vgjrd1Dqf^%&-BC4%={*^C^oUxyC
zg*9zE=_Q=}=xZ5FY*(c%46c5mu}frD1_1;2nwvh3N#wu^;CUh2#qM-+>u8+*WG{4g
zVAAoB!?>;TVi)v(uwQUH%`hOCmPS+sTeujXSIi|48UB<!|M^-+5Z<8jQu!-BiBCRM
zBu9M#!b90r7b%xBV^HcZ^9@JEMf|vbL8Qhb;|&z)b{!mkWLfHGuCoWTWSC<>B;4|5
zw%PHbZbz4p-e$@(sUq}Iv?)AF%A3e>GBmN^7NK13Y-BraJw4Z8M&p@rI6M#jd%=~8
zFN7##E23o#!5%ldivln>zg$O3q2ePzJo<S?{kKuo-ZY&ykMdef3=k!uSM<3+T{m0B
z*H%;g{8r<k^)=C2shw#p-?zeTmeU|yxbo)YfJ5~7SdqTw@{?M4IyI&M=%Z0n{A<#g
zSG`gttJWh2#_l0x#TUGb&Wo7E*)d>zZW)qtbLYSr+xg>xe`v7872^q>Y3k{4^`e13
zj<LobOcOSza_A_{IDo3gh8Xm_!65HB3bJ__eBiy?)!YF~b`#o98~%WAO8M+zu1!eX
zMCj*CuNP?k^JwTr<fuP)jN?HcV3loW`ci)-<vTdvr03K^m`J_768CeDeqD8h8^QNW
zaL0i2T->iM-=oXQ2CH>{Pg?6;Gs+_-=clvyM3rn?hI$MlQ9mb?mF%yo!4%Ros>f6O
zXE%JT=>^|>RGg2NeuTl=CZwkMjZ653&#$1JxqrlFy%1ETdNUi2{rv2*vz+?uo0i8|
zE`9UEm*vUob(H^@1v~IIRD6y~*FFEGS4Zd3lvYQ?i|Pt<q`;?f8+|A(I50!%4%-34
z#b88Ra6g*yWN6v;BU$*RNV-W*Z|LQd=3&O2Zn1IAxKaZL{p&yoW)z`qSZ7yAyCO*t
z?8yt|^+B=7Qi%k9b=zX-oVD{8@jVuzg}O{|9^!-6Yb?$)+=Ivy`{85Jeud?6#4UxN
zn;+kS6I#`F<7_*X5`TIwH0!%Gte4`^Z=zFrp{PgNyozg~D;0l2Y<9u0%T@?I*C%$=
zU`N&JQ}90F{cnrM`&go3=f!7%u_>P|eu%{6cO}C$-<s2msH5ZK^zh5UlKf1#4$bm{
zUpH>qBY1!ZH`T`2{!pKXy!hWN-;VJ8{qFcw&^=;3YWBYLh`Ky**s!(4ZJyjMwNk0P
zJ}y?j1PUSLF5h7>)#SAiFuLUzafqr^#Ka?)Up1WJ7+@r$->bH=tc+vch~rJ2q@!N9
zE!vb$Xn*TUCv&$EXr1J@zOCd<3Z2Vqy(mj3{*YxbZo74kvpSyxB7jrT6FJsLVG`<u
zL7W-CZ1q;jj;^IhYm2)r%0WcG#4dODV)7H#4Z60k&PsL5ySH3mS%}$FonSf5D~NEJ
zO@tjY9#pzaAX)gLmb%)DJYe%3JSf}#h3IAoE7lS-{}<H0HU32r<H2YVTTv3~e`wo3
z|FdBbnXe1plxAl?N?NRb-^zF4V5;Qn)Z^4uOZy<fY`feYe0zxUvp*WeN*Ta)GgL*N
zdALI3LX_D9S4Tu0w58E&#11b{hHfQ#{PE#(k1+{EVS~0Fn(LWFKF)yh6no|_-_1g}
ztA2)CLde=%asPC^d{65FVk4u)P^NftQ#tR+E(QONVdFc!LOhJiL<TkHdT3%@jnL#i
z7SDDRXb=Rm^xS#4YZ|eu|EL+kue!=BILi!vwljWxmky6ep|bvgb2DEnZqM^F?f6hB
z%xeA5aTz0D&`^&HD26^J-qI`cy2!V%`PNig94eMTPSA0LVh}F(1F5OLbL}y7PfEsr
zzpehUTXQFPK3q~a39CB?nimdlHr2u*&&%Bi6D%^fO3899OKEZUn$;CI+<|M5_i&3K
z3nUfaZixT2ZLgz}HI7h4oGA7}R3hXKcz_OW6C4_FWz?8r#K9g35M$Y4-UGy*;2)b|
z-^+C1AIubnKk5|0{h%nPU;%Fin_f-wzEAaPDl|IZO$+`Z@@OFPlhvX>0d~eC7W7ui
z=}f`j*pf7YH?wi?79*F}Q^Y!2QG8;k|HxHJp_G+K$Wx0pAm{3IPHn?ir_|Gbf*+}6
zUEUu=0#p(4hSWizZ_po0s{|R>!g5HE+|SJT^Ix9*$`*JtfZ=DS+v@aK9bgzs*R_H+
zhBn&FR+Qai-@71Qq9lmQ=MM4P{Fe6jJ8LQ4=nG3NM+@*f{?Ia2?km6C8y6GQifzEL
z)oMV|oWt+92>r55iRZidyQ|ekXfD+CEE<D943;i2BVF#`Wr@z9kw5Q_S*XgFNpg1z
z+*{SkwFg(+2yyDKj`2QCC$7k0u9NLd`{bmKw2g1?xM!R8yGZ<7T>Ariw9jMqO08EM
z#V|qaul6%*tG2*Vyl<Cwqf*wf8>Pf2XC#;rZh!~=?Bu?f%hF4i=kDL3gRhcMrDjrb
z;jmF>J8(dkdu;bxS}UQKHWN)=3ZYaB!?m3jtcR&?`%gEBvsD+P5jP~{`TO5KHzT_%
zlWt69KC&HH34PsJyHhvnOg-c=xXx>tF(>A@EJ%h8cdk*QrRqSNwT<4-y-_=sZH#pC
zM_5j9B~{hTyLwgxLKSoGSfW#T`uPdEC-80d7=92|mThZ1J0il@AbU}ex0k_&+oKt(
zDpTO==w5pe#u%=r?jCoOAgMOlqS|8uz+}@Am|9_5BKZYvb(To6&jH$RqS)hA=RMBf
zy%VB<pc8n1{{p{!>%Hpc8=RAwRA}&xO?18jOLKQQ_eP0aZcqA!29GS;GGuGY+H5mj
zdl=@;z%SAVt5R0qI=(YU9QZncEi--D11MWZ;w@ADquY|YD`|`&6|6tEE%hs7T`8Be
ziXYPQia_L+;clM4PYX&m#mdzU&SllX;c%<K68pIXgN_^BWle=**QBRVq%>x`LTL1_
zt)Lt~t7lgtT7~|C{rk($z+(fI-aC1BXQ|@~a@(iE$T@Wxe{qNnKux+0;Rz01Mc+yA
zdu7b!M`)3iD1D2c^4^UH0R2d0=`<9BnmxXbxnv5h&*^>?jxDeCpxM#<*T_vK`o^7!
z*aDuo=aKuJ6oL*GqGkJ~4q<7mONum$7=*z4*ZwkQ-m|c~Wap`&e3B-FnmUtY%DLzJ
z#Wa4sfq_4f&}dxV{atEAl4b$;tLHaG3w)%R(gvcMcqMz;`mg@iZ4V0L7Z#(w_3K#%
zjl=jc6MCe%^I<;sYm8@Ld)6cOImkNY6M9w8#SEWVB0aE4$ix55?7*FFjpqZEp;QbV
z=vKIMBVdw1Gj#<0{D>j?!+YO@gxA)yw`wB(N>49yF<WI%FM8-l;g6iNBa}Lp$u%i-
zg3#DUoaE*nKKdFjMr3~JovAj?8T-ffDes`Bs!-u=tM>WAd7=XTX=0k~v*^d)XKCc0
z?FsosmP@{M_>Yk*&7O)~oxVYVkV-=(ThG*=3Gvp+t5KqE7NcLwqI;<)ee!g{>{w%0
zoWQIr4$iZxIF|xrq*C?S>`3u#^u`1V2uSMB>BU4H;&|<yDc+c%XY~XTpcu~gZP1!=
zUg$!Ez<(kev>7B2GQ2A)?O<(%$8Q<V^|?8-I70y(Pe0_{$e;|To?LMkI8X6+SnxN(
zu4=x1ZYLnO61uBQ_u`<9nX}c4!}AhBMQ$~<W}6gVP`BO6xy0Mu2d<HMYN}*=Q~9@6
z6F=uLOKASNf$@seTo}%SQhXmWe6^?Bnot>=gos17Idj_msUsgHNVc~+Yv}zTE`s_}
zlS+5?h~qBkDe5bY6WEfz<8616<o0PxS%Je$v^MDM9uHcHZ!n%;UP6@daoa6bzm!oI
zyqNf?s!uJR()%rkdfz^YZ7b8vk&F6nf76rD+%mcLnRmRp<R+d*RB!8@g!TgHh#0``
zgRhClY4jblykhd9)7z6*g#OhFaYpq9PhUpy&2#s$tLV{4ep9?d^yAPMlk_8fvFFa{
zhI#JJu|A})+3a&Y6u+eJpdc1v*AL9GOv2)9+_=yHaVC8Qya<zP<7|apdgXtuUnROJ
zn+7*PWeD36sXzIXTzV_hJpZhIhva;K%=`O9h0ObzGZ%h#9{lXgD8eT0<Z+be_U;tz
z{SGk;1JY5eqIMTK=b4ciWNO`aFq{K)^RD_^U>kJ}*Nhdt5|~BVRmxYkN_~IJEeA1|
zGD9u#IQ-He`>L)_xUi!lBh^I@o3j!n)Q{rA_b6KRJ8tAQ#;brQi5VampY1`9W!j$S
zo;gYq8Y6n}ah2~tXedVZ16kjVG<aq44K9*hTspzkNIwjE`?APs$dOj@%ywRd%^+mI
zNC$3y@`*seA!6ZM&d-$l`d6{L!>tD6Wq8h4HlI08GIJM@WZQZYld{N{TNSfd%%e=t
zdlK{SSn4@BLGZX%?OY(zt6{)=>MUB0t0dI<R<9XhhhxS@sfd8MhApwD>WC9*n+kKL
zUnK{Rf<&iDk3gvDw@_TI*n@UjHt$L7-1@$ExgWJT{cAJdo}G=jtAmN-q05`$ksz$N
z)T;vD=^ILrBgqko8Z8O9zz2dhk((v<BFrD>G-^bMf>x`b6LaRRCupMk%D$xOAWNVt
z`F=h0$5Ob8$ZgEa)C>#aE=Bl><>rmmN?+Yp7ClgkIEK}LPz>4@13jXw9~VxQWLBA1
zgTI(RS?|kr5{a;1PY>Kmo|{%ok=lPA>sQ#{A7_2>e`*)se$#i@Q%d365wIiNj@KT3
zy52#R9uc^2E_Au3hJIf?X)eeJxRT*JV-Gm8tuUtAhcJ~wR_?z9r3PI-+Q5Z%@SJ#?
zcuwQZmfdP(XshzT&w4g#)u|yZ0ik1VeR;gjuU;+`yPkb(T^kJ;Rt6>4UxFsVFH?($
z$EeqS4gP*Jq;DjV2p#opW+sESJS2=tD9V+3_QKs9a>wTC{6{-%^#{zTBBA?W+U(&?
zo7y4ji3^6y!>fjgyb|{MiNW7wtI>}IZO+BVGD^T8ZdZuntfHFb&4I&*!GY_Qboe*G
zj-Tj{1NZhZ`P14S?Qi^(aB9E3NoV<lgFdD${_5ll{juNArlKbu{hZ1$^3m9~$*~`c
z0FC|E20;Pdn5BeWbwVKkB`qWK{>(3o*h;-)j39bzE!*T_Efh;=qTZ=p;0a^*bn!RS
zIHT1IvW3P#{>v5<-L|B>;8KkpuR#n#KCfV{5@m;303-pbL%pp%Gw_mybxh!Yl)J7%
z4FI54q0Qbct17oeQ7N);N%{_Swvi}i@zKY|Dbeh)NI^^cDuNSA{qXr#siN{8csxr4
zf44_H(p6GN63mtnIP-=8P|9t?RNRGx**8nuAkXq#zcu?D+I~*s-eg)6U3<^UFI4QX
z`Pt4j`j6t?%Qxw^+6u!)M9s4ae3-=?G4JK7a=(DDX(Bhaorvo>kDSf4`?_VEC*haU
zd49GSazBCbV_aM0mD*GQxm%?A)5!Bf$vD5|^V9ZYg#SUu3)eUTx=(-Wd^nhhy=Xv|
zK#3htJ=Ok^D~E|)OJ<^vMg+r*bysd=<t=hz_`26R#>yLNy_9UQeF9=+Y}`Wi0v&q3
zb!Z3wokJ**w#H3DbZhP;^u+QfqF`;?w3?*tt6ZAKRkLhrC6b9pedZUlUO>^rlBc1e
z%d4sIYusIfZCA(ME7?M}a#Wh(Fl<hSK?#Y-FH|(l`Q=sKWHfOldv#qGClPetjjr?(
z8`xt@hxy3il7QA1mE+iKMfeFS1aszYn9-Fxx6r;FOzhg_n!?IGLU~5r%?#$BDEGGp
z{>r&2$T^=ihA5yeZ$B7TUaRYRA%Ch6Xk4cEL$<uxBUR%|MWMz*E+>7(LFQQcSltZ%
z%x?C;NWS*dkNvaqhlj6g7=$eWA&rat*cd74lEyS?uf}Q3W5uw0C+HQ#toqj_Se4G$
z{vRX17rnrAmuR>8MM?%qbe%fzbsuenI*bx`IO?Yv)rA7tAYzXDpC_1vj{WBDYF~Y{
z_9<65z@InhMdTF(E?pBFi9Ae&88)iF!->_iu99uX2GyK}-lskt!2tcriw02(;Y38h
zg+2vhX-hBd><L+@xQ5vzsj(jehb-Fi*4X@Elq%(S%ku<0G97FwfDCX~pks<d4i30C
zw2vjLW1s+paif7{7_kGwV{lgndd{w1mQSnu+tS!|E4}&n{)~ChRY;f9ayxi+Gx#M@
zV2=fV9w`|lm>#;qw6O|o9MVTs*)aD`Q;zQQ-f-n(JQh<Ti{>@DmpNauRc?835Vlvq
z7)<4tB3H)<b?33jh&aRw`px*g+FL0GFx&;xvY=>ia9vn25YLG8HCeo%8@?$xJ{f~`
zr}5+c;Ju7`^qF2f55yNR)9--&z-8Ufy^8*L!jYJWEd&?IvpC}0=8W71Hqi7iSAHkB
zD|mZ-)~rq`ct229KpsIkXIZ^cmg>HN8zXmEb|_2<iCjr~?P_3`dMZIuKSp=ab!Ol&
zTBk9ZtHl0#A&3zvEcg@r6#NPleo{c*lz#Y?v*#f$a`-(#@>mdXu$V{36YOF1q(Xir
z2!YQXUy0Faa6%HqBJ~y;B*L1eD?v^4elTuazE8aYWQqoj8#jq@c3cNp*lW{rXsxgl
zX80?SB1u)2pM0NvPRzW4Jf7py@K~3&_&SFl;(nxP5ii3w=)cvvUFgxwS<fcz6%WMK
zY5~67gm><b9GI`+e$VyMrBu+_`%uY25iNxK2xC!JnY4BL4emjMb!u*wg&op$wtuP<
z$*v(_Y5f2(qy*lvL1I9LU6nU2&!3C70xp-Pe$w{g?0~nwdy{2iUct9gC${x=k&rYM
zS+L?#jyP8>DBL}>dq>`l7ci}NW=?zvf2Y<C$lCuD6$Nzb@krOVl|a`BB-NR#-sRNB
zcHm9EPJhCbdO6q9mDNVgr>m~!)9!%#l&HV%wk!f0>m`++__Qw?h_jAZq0|%q;M8ub
z+Ct{jYi=QFUyB40Uj3}icZl@Ug1J{o(IB*jyDzh*VH=TF0R97x>a=lnaDs-dDk6HG
zuR6l5TBlwk+pA@#Wtwf>@qx>_j^k<Aa5#2}lyV!O3}}+}CMxXVdKp<z^s%NtBgsNU
zdA)Hel2w49nClduYn1`Eyn$@5sr@|YP=P`-g`5sH*xOdGwEF{Pm!w^{AN4!_qtAA-
zpBaMzY@A^vbZ)^<yUor|fkG}byiR^N?S?3vnb2>LC>(m1?Ty9hO2lpg`v-er)i)q=
zcdo5SJ!JUPqSk)6y`EN=(EDa_lk;IVb6wT(CQ6(Au;78ENPyS;d8)UA#iWSg!`mWn
z5gAXq16Pn6PVeDRg&*6ekT@QiOz!qOq#*_=O*+L_xSG3a$MYEtSo|RBcKj$0K)!F(
z;ELrf`xn;Xm&l)tD5(%&>KYSz&V`8K__a|G9&!n&U?3&PYzoabRrofbOd;pJjX=KW
z2IFB39P)$lm+J~kVu&svPc?pOd0!RgIXJjmVbm(#W+c;;wF+UHsoo{WMki>?sSEp1
z$~C@mwn?<i7i+LB{P~g52}1&bMhW8m5ve=J^P9icD{imn$KsL1$j&Rk8;(F4o^zMI
zgw}Co{A<IHb&4y_q=g52^yL_=;di4-FD{d{;bn$uVWgn);5N{)Y7J^7?1MvgM{5lQ
zT(U8U$rL0X!1-{j(RyVsS{)D@kHUbF#FhqJd-S33;P&e|7;lEbzS^k-<hFCj>@mMI
zG+VO_daHD`%5K5PoVQbP7n7lXBBWgrGG588Kl3aY_T!4r`Tghz)n<tIR8ig94wA(l
zfsVaAm-F~n51dE}P32(rFL$ElO_7Z^?FHU@f;g)aPyuKC_B@0$H=PXHXTg)VOb6{v
zwd@_i^n>b)Os)lT#0J-&E?8!ik^9){pBn_rMc+d%7|#-<?!CNgcIzK}m3w#QeXQ@_
zCh;S)GGj6Zy3gC}DDgK0wQVTO5&4u;da^(CeV0CEH;|OSDC^Zx_PgBrovcX$F#c7(
z`a5N~W9$mZFK;(jz0L0rpR(MXF4$ax3cnBWp<G{t0kPD-96>MNje?mqHhUS+E4e;7
z->p=HM)2&d<hbvUj6yFI?L0vuG7BGQ9<nWm*87p)>&lH^gPlTl&Iuv}2}S606N!aK
zv%(1Y3?Z6#Ey$p`L|-?8qu_i2sOAEjVrb1-%PJP18e$awJa<Whs%TeYQ42=L_>=1s
zlK@3c(^?Crpl7@Zi~lD~LDXi4GnXe$8O=K#Eo<^OX<~Z55b~C89^_}9;|XQsN9FKu
zvXw-YtS1W}{rQm)j^$a+IaRx|abaUE%|!x2x?j!-!{AC|Zr*v*Z@kj+9B+3(g*J}O
zJ3hQ@a>YCIYk0^8!N_wxxUe_GAprfGgoz742;v%*4IZ4H@V&>!wnOHK0SB%@WCTTf
z0ggXG1Q^tq0m~Ue)&`^kU1Zs-Mb)mBqMJ+Qo+88VL@#rqe{<<ecTvAtR`*(;2nw1^
zHSfX*rJN}4p7M!-=RM%<a2l$@+^E8JDSQm299$OpoxcAQ!vp^T<D9}BNz(Q+Z9|e`
z;FIf68NrI*v+K)S`-bBRp(W9a!u7DK5?jDf=x$Zh*kge7tIho9A_dB_zWiKW2joF-
z)6H7BJ52N`wOzfd4ZcAX?AY1&WJrrV^KV2bK~m6rska?N`0?*O^yf<BNX+(Pia^jI
z*q%=%rUZDa@MUgwXuXj<!9FijiG~9A8&}K5j3h!o0COxL(4D0oplejJs^<HI(J0k>
z>lmxwH5P6*pUct%wmhkHiitan>PHXk6GL2-E%x8hLM)L#*H*T5b-q1#{i#XW{|jNR
z(sl`|Ury35HzBw|W)FK;1ik#a=hVF+Xpxo|$KnHqp6mmE3~XGPw|?GV@kEs4Ii27u
z0m1yann2~LxJPIQm9b$?>Z4zMg-@<HB?u0Y%=?kg^6~tlwf-IL$gI)OLqt({wLVhW
z%`Ox(k`x9owd`P%>>q?_eZ$q5YDGEp!+o)u1Y}gdcynw(P;6}OtO7|Xh};-NYYoTL
z<apf^@i~?C8%*9osDScZ`OE(dk$tEAT&(EgxUHR&|77o!MrqL2kurc>d`@UHZkUNG
z?p;LxBVhnmHTJ6>5p%fqL#C{?2J1O3aCRm)9NG8=7<FDw%#ucIegEEIUr>GlvB{y-
z6Xi0$@_7uu0T#H_>n*o#eXo<oN0lrJFX=>jjE^A*KM74}(mzol<S2ajbySQ)UonMm
zv!n<Ah5chYSG(qGmVYCP@3J)bFq!@SD%|_go||$v!rvx5?A>4z&vjGjJsii`2bYfe
zBSegGq&B*<$e<Dff+LD*eK_2(aKShoslbfum*v3h#f=+CzGVwiNp6X<zMSaW3}mH(
z7?UmEr2Q{rv~h&zErSFd33{IK&c=yZQZp>^!={Ng6VIqE^&q{w&M$B~_E1h>Nk1QU
zIgqxDZUbSpVLXoZi~(j;7tnpZDTdu)(`5Rvk2`sFu6!mq@Ia}+>%`l+`CjH)=V*PE
z3zp;L>`<6)0(q#ww5i^fBJZN-1m7+taMN2?bEP){FuYjrR_cnDRD_3qm^~a^7O#^c
zN{rOpIzzs2vV)-Rt|0e|bFLeqi|aE7N5h7~)M72@DmHj))U5rEAjw#9e5m)1!+p|O
zue-2f9_^sDw`!J|djsefR0<QL7(R)k#m-Xg6?lu*5B6*`QS9MG$s%Ign-E6rS?)G)
zbSMnqu^`*Y97o6IHg*irXu-Fb{{%GPRYwo9;0Nx6iN{<061&S4MPPou*HZUaCTjDe
zaUcC`gwjihndr7m7V=>Gi;cK!5q3aB+qPyF6d*k_thOr1KnP4ZGbUoe5U0!XU2l(s
zttKsLn<JJ8+!+Y#s=^ZNtlBrs>KLKxntU$SL;Nl$`gnhPxtb9%dDwSjiF}vTZ8iQ=
zaS`*7l4K?x`m_o!afwyx5+F<HeJEd0(Z!?8K5?0Bip_S-YfDSBv4A-fOZO&qS2_sJ
z(F-fim;Bq}DDJc91*6f)iZlhcRcNiRGiu*^+z<2M$a*!#&<COPxvE_y(e!O{`qe#T
zaXMRquf8(@o=L>jU&&#W6R!UxNK1E*`)D<aET-s(`NukF*9#>`j%wcq3%;!m2*t}@
zLz<Z-zqR4h3Ty+n7J1t$gS@CBwiBY*qgmhw3_1)@aDJ{0b))cMH$dk`Gm2E$*L#ho
z)R!r8yMgV)k`P?>h?$a>6L;28{df23lqLj&$wYm;2Tgf1tI15s6#9Y2?qk_}t`h<D
z4`6Pvx2~vBQ(m1}IcxOmHQ>T{j0<~y%{5O)5eLsf<X)?b4!_@{5QMtUjQV9L(qQ;j
z6VP80*~QW#$zI2tcAV>uE2GU~Wc@=FYex+=Wi2T!@NEOy@1CTyfQ_nTP@LdzBt&T+
zR&@Oi%%6(Zl2^%1K~L#~dk=f^*ye>*25w)`^^VaffMmjOUsbIX-;uEaU=8#0P@3MO
zZyw2l`y7$kbMKp+uo2Z9ZkxWIjRk%z0rPV?D+YJBDutQudjw196z^w3bJ(??BgVJ|
zSd=Hpfk%sis6UPzYMxJaV1uK8an*bL$c&g;8hol5051*DL`&hg$W@jrCQbs8w;nLT
zcy3aum7z#E4X4>K3Yzo*rECv!?&V1I;jHD{(GeCG$sotFc3qUDErWgaHw`QLPbKU7
z4@X7a$oU1uy5@$#68fu3`5CQ7FPT$vbo4&wk{$U8*ESf0ZkA(3)y|%!A~^(<s#!2A
zB>PYYFD2L)5%+qs(R_e(za?fek!9MN@+-iw8?HTk_iO5pRXBt2ky4sRL3rWrTP#je
zS0fiIf62KEPHYRc2e%PACJ(O*Ji`GDZ^Dix-_>*!@A&&%ZwXBiBB=xSR_MN9;ARxU
zV_eitAT~56UBGp&hxFxp)Gk05aJy)!1C6#cY8GqV9&C?`|L%P$*i?hdocE#aQ^YJ>
zs7{X3#M%)8VuQQb7a1jEG0Kv&qfVUj$_k!N#vJsuRUd+BiL201cW+s1I)Xs?RM`27
z?<D}KP_D+^Gv733i5j$FSp=syaA4PEs=)pCts^PJ)bijIp06DkXWcdPtw6rt6y~V=
z`a8rJTByh6u~I{gIk2-+?npJXI{~DJ;1fY5aV^(;*DL$FhujVA%&C4i<M{ZG{v$UA
z1P7FOYso?-1ix56wHr5|kjjj!^T<E>5Pv9HMU63FL*EKlYLAzGmfMMIvAEi2dq2|)
zu!h&hHrp}a7w(WJZyTHAD8rf#u)ZVrl@-=R$Jsc=xT}f{d`iZu^@aKcQ^v{<xk71-
zdx7tUS1>W-_w$cc>-7+d!%z8=^WconM^RPg&41+XOzBm`n}0Z0s0^@8a>w%Vv^<{s
zs349_ogj63{w1Ksr=O8qOQWl@AUUVMs-3b9(qj*L%$&5(*%vi!fgJ@ty4^bUqSMQ4
z=-KIW!mn4>yFecJ5u#Y;`C6A+yDY;epvVgmfiq1R?%HZu!RFgfsoRpMfW=JXb||?w
zjH;1CL$n*V%f+6jNLda1)=r!)$suCTBy>qNsHZ~oa*w+xpe2+tcDhBniynoBMV&<H
z-7*%i(L;bZQa|)!Kzu@NE?#tRjlv<{LE3oXjF^|-3gL`(iqIcaj_~75iOtnP1j$ei
zQR-Eg7De2^DWZ*a5@Bd$R`O8M5)eOao)az>Cpy<mdbCQH8aTT@6)@)0LpfMmVx4@Y
z7~gbqJTEAzup!3xUgP~HdQjr0#M=_WplH*Z^wi*iy!OE7r{j?yUE4w7tkUGn@ZZIk
zWqkHO1X!8OO{Ckd!w_xLCnWG;l1>jQr2^ZrWOQR$Y!y-9-n&;cnNgG?vnA)Ng=?-#
zzp_64aN?Hcil*3HVe?Rc#e*ifo95G;|B|!1=;IsNa8?gr^(*G5^mjQfLb*Ba#63_l
zXl=xhDuM_=on*)hzbHL7PD7Clm)kuj7)X$Dm$@C6F{1-XSp(6gI><ylTrS@B1ov&Y
ze9kENwn`|YxZp3VS&LrCs|{urq(Qc0tU+A}Wr2B#TRpzd@e4JVa>gP$;|SR`GJoMX
zCB@EOo$7b?hi9gYkWsqR=+0LQu1O^qf=mD2B*^VM6Ra|9j;xjf?F#A38JXkb4_-m!
zA3Mf32v^Y@k{*8`?5STJs9sO`j(Cw`X_!=H)Ojz2EF|SS={-7-cO~JV9;~Stm(ul9
zmEHdIGIEkBKR&%ZNaYr@+uO47k&m7wFtEbs*W7GK&-_aCc+MWCD24-Ot4KB6B|NPJ
z;;;ljC-c6ED}tMS7-@Uf?Md2y+zW&r2?r|Or)b8Q-qW4EQ*a3=!OzEVBVPxQ6Mz(&
z1_{6MjiUVUO@+(Be0kDB+y-cbV9x4+kb8xyDIhWPA;uJIF22AjIEEOiJ&gz^z*Dna
z4h1q9z^mZ42G@H&2puF;=&CSHS^bjibQ9{Qyz$ExEaD6!!H<@qRH3}7s=ge#D&{OY
z+^`f3N`mXGxbx=a|EU`o<?pc?@rW1Cy}8W3Q;cb#YV9Wm?KNKDp25l&m<u+c4QImy
z#|17fj4S|3`q(3b!%b;mZf@=#Q=@o42M*&zy<J`{&ow(s%!(`4pfsGa=t(L<mEXpb
z`!3Livj<EfC;irr$Wa}FHqOL{=-jtDQ8Qf9jf&O5hJ{165nh%&kqw_{lqKFfJp9=o
zC0p~X@7N<)oVuvzD?$tWnQfj$(FW|Tr-bu9@ZRD1pe<QllKuR?L7J^SBTriKI$LYT
zOKX=u%jgF;!rpAzYLUEE*B8*etY<l0r3Qu7Z(g$C990v4-*<W_6$X(RhQ1|8h6c|Y
zWeOcAIXaR{sKBWZ;MS4rZVRY8l_>EmOOJOa#!?+sROn)<gV37YnmKClyB~Da!YKfx
z{hI<T>Ja0ESNKzp7yJ#^ZL{*>>Fjg2M_2|~2u#}IVeM(}ncw}&uu;!qAzdd>H%Na(
z6fT+xe0Jw;F8B*<&^&))3NtG5j$xZj9oZmFaK*>fjAyWQv-YSZmw`Fs=lJFR4&}ET
z8#eunQWR9P{IL>ZF5?&K7dj^2j^NS;F0?~|835~cPH!FjyV`q$%z+3aKF71h%)MSB
z)APN-=MDAKs@em`^PPiTy|fb?lc$v-(%xmc#b4E5!aYs-o-X$pI%l8q=Ie_;mR4Rt
z&t>&7cJxbFcD47gBFFJH^;U3UR3RBc>TOhSwqoOANHL<3t8u2;_J5=HkRGQ1Gb_3Q
zk40Bl;r#Chr*Rqu{bEV_dtf`Je>Ma#kpG&gK1F+kZm#g))5aOKnb^@4I19)KROzrK
zUE>#CJFhAd8z33x5Jo{<G2m-r-9lF$JrMM`8CqMQBT{Ld70GJu)9@je`vNsdR)IC#
zUg{OE{8$9b&F2;kEkOsj-M@sV;Mo^e9B+nbRXe>$@RQes<<6NwBwDkk=Re-KPmT#~
zV_U@9AU;>atywZ=*^;Gn^=2FzFz9;{5}UqH{>X@!o$Rbb5}O?`bR7E-x_SnuG7hn+
z`b{PZHeNLD`ADt&!CurIq=kVIJ8?WVj2ZGlhw*dPu_~o1$Mtmu@imxz@g%u+U+18s
zC&0a`SLc2!GJeI>8}|SwVUOs(JZF}|`1PT7jYG5ppEh<U-`Z&1?tEP79#dXF-+kEs
z?t{nA*YU4SLr&P1$E{|J&LYwDkM&=m%ccjJo}g=bM!fs;CwV>l8`68tFUgrX6Tv+9
zn(T2_RRPJ;AE=wTec%}gQX0P9P#_9NDW282F2N_FYRdcTRUr7RAtulChXinhr+2c|
z&}e0CJhcm_BjAJ>&OJQl-PLZ+xVc_b7!z^bq5}=z{i61=Wav2_HrxjmOjGr(Ly7t!
zJQH(E7l<_O{T*7^_^8{GKxP$ET|^dL)N1!XwG|`Ion@ZwMDmKOlU%#H@cH7vu$iOc
zNoc2lPRH1*3Z3P%toxoCRy!AxU_x|O&0b&Nb<MH#P(y-$w$jN(QGNx;vqyRNv%v%s
zPrRWwv@KQ)BfwRzkcy2Q-gV?w_RL-MK}gBsY0aGWT%tlsZ{82D5I?1Z-h?;E4%p;?
zJn%NY*2u7IE*ieH)Bqn}O%3ab)mAG@wGno)UchrRyd4mu<=?ASOT#Sa#Czi1KEA$p
z#Y0xef3d<8f<dWoc3Nw!3DcH5aIrmRON;GRo5J?1m3$|ibtKaOOhnV>_wuV-2NU7W
zwO@_lp>winzhb0=6!|+)$qv=p0HvGKJs+y$PRLRiZ&|gEoIvF)qzPWd*b2fb<y`#Q
zvkf^(8Z}cB1?Yx69npFDXT}kSseHt|rg(JB7dg<|>zE)QCY_F7c0PRKL*v7kVJgl#
zNP-%c=-o3Q`?XK{8^(=IQ7mp6M!pnvcZBuuvWzd}*~(i?BKs=on?eFPH03WKTZXew
zkQ!EY<qk7zr+%-v<mAIyN#bISef!Fq2y$Ma@6`0!6A=d|)nlpg0u|?_|6bJi%guwx
z9fg?xziOQSl>Fk|{;T9y^skcNzbc<$%n;S&t({%aRb9TC8{3<Im36QH;z$5pP3$}=
z?W`Sa&CUJ@XO}<qzlF#!w|IZ*fBUpXfFuJ%ufOnUKVLFAtnz;g#>@ye4^9yyKn&Vz
z$WN<Bu#?NBZ}pSZ{%TpOklK2E_5~-AYLPSp4sl8KZK0?a!xCl5%Cs{qC@9AZ6%JwW
zHC`93{t%L>|3z-8UIE?}IV*j2ck~H-=fm;A_l|^7IaV3nL3vMMm<udvZ6_62&;O~^
zbC}M@H~qB$m#~^{XKI{lE)zJcA#uaDY%V`;yl{78=3|HfS22N$mZVg|4ink`DPoPJ
zt4u3z#^$;e$64lk;4MZ5p!HJ-&16ZV;J${^TsfU@42##JK5yeSl{DX{&eEPzj#Au_
z6OFN`M(V^;GTc?E4h5v83SuXQC@tl8bf`+4*)m41^Z3SQlI_YT4}l62(kc`PsCv83
zAE#Usqy24E{D9x!RPZ(REvt0O7O%%bI|_erz64%Sm^|%o<76SyN0srg0$AGYhqpUe
z#)f~2V~-LZ9Dii7RWURKev^p2{~6>p@kv#C)#n{`uK8mhilYtUJy;fLDT?^f48J9p
zyV9y;Qj$YSf-Ucb)YU{=W~VDrzKUj>>{^>P;^gwGeC_HeWj^H!Fo6pzVr(DA7r|dg
zDr<z4sBl21c=CQ^C{epKzlAw?@}j+YwQyj|^S3uKz}Uz^qvc8clRQ2F$-V00F2Srt
zO*$2GNc&z-z>~3rfTP)f?a$^mX{4nlE-Mx^FB$#XZYz#U_zD)ymcSKPpEQp&A8h@3
z8`-*@q4oqT<XC;2>h&Kh3ANtYXGBywl0sqpb)!zavZaIm*tTqguZrGXfytd{pL9xa
zLB3-qp$)tCqX85<z-VX}DH6NsX<9VgqqNZjq@U9nRoUAtx4p>x;j>0cL{X|8(fd!2
zfkugU*fM=zxSs33R5@npE}ncQJ1<V46!lL}z_!K_ROpV<#@YLg*vz!Luh`mx`=)Kh
z5FftWkE&XI{&Rj^U1M9@Q%5%&-Tf}QUjr3YtpSd0$}y)L+}ZqXai71)?km~YD|qj<
zq4nG!T=^gZ>PwvxJ))gTZj@km7uRao^}u=`VUdLpWt))@R&3Ur8C8yNfK+k2sqc_o
z%9do?5OKiE$n5tuz3LuZ(L_%~N6$0F9#L(T^!*nYC5b84wyCmHxGhy!9t}Zc_!(bj
zA@^&?bhW1qy(%_8KEol8;q0r_l900|JK9;gE`__&?Yi_y4!2daJ<cj6{OPWMIaiH=
za}c{&uSoelz^V`BgqTx-6xjD>-1XCFmjQHC*&qZ%KU^)$>8%cy^4<l-$S!On+P85d
zwC$wLjJd&WaRe)~x8uP{o9b++k2|lHFSwPFR@Rk9sFoNd1G6ff0n(MRj%IY9FGDca
zFQ^KwKd1zw*QeH}<Pcs;4XUYqa&t?|F(Q4ZLYHEa5uL*oQQcQ~70OTK(V9Ci5>Ezz
zFSj80&W_}(Qn8hd8;bQu&DhW4Xw_oV^m@^52U+`Qe^!!6Tb=1gWcD>InrVMKA8!V|
z3wtV2b7s0wcIXz&Pw@WG#<vOkU)wGe-4|C2qKl)K%eZ9zqVIffPe&lZxLP{kXzeyO
z5K&V-AG$|xiW}1^9_c6fL4}I8tmNNyzYFI!;x=jR2G?cnj$P<=o@607l-}=d0>2V<
zUVoXyB=WGeouse?V3oC?w}J`y6*PrEBA{amn-HLifr2MEJMS8zB-OhP?ZuH|fo$9t
zgS;L{<9OlJJi8MgdT6}#uQZ1*-PO7*j=>Z~L`<3N9g)HHuEdQecx)u?1EaoU5SNW#
z(mQ^5TKUOsp`6+Gw=?&Xw0{<=UfqbL)JQhGy^q)089WcK4u(0axGtp;E41geYLsy2
zMK~t~<cMc$Sr8Lcu?qC9o1WUb=t7qk0{ImZ=i5bINSL~Y4ye(?6Q>{p+330fZa<nn
zZg3-onB|^YSr`m|QdK6697p=XB5ovvV5qxrPn9xV*~XgeD7NV#h6pBb67DDWnQ+v#
zI@RjZK2MapdeBRlpg<6tydxZi?Pv_Q0vjW0&=t!0Ki)IWAz3zO*6lGjagA$}!)tvZ
z?MFK}PLhUlj4jF#?z1@73Ky?*ppVzK8PF;2d>HCnzRm9UmEm%XN7rZSwuX4@5zUR2
z3BAo0+lBw2itgtccb1D+;fGaSbM*@fw}p4S{T0@V?g)IGbIh+|yRWytgv&e8l2G;3
z1^d49>4NvoF`D}Xd00O5vcF@lj<?XnCx#(Q$me5q(yRFG`7tMHn_QUe?|I@!>$x9J
z_;>nXYI@+UfWmsH>e-_MRyd3BWm_ly>JPCYjku^xe4IG^vxwdt+|x;yx8i*fb#01a
zLz)CW4irOodud5kIyYKhMYwXsEeXQE4T@%KMe1e>_{HnVJaZ|=3;)@;M+xds?c&Yx
z*jc|Gscr8Rra~q@sxZ?0d@>7e6f=|heACT^b~e26JQ)l;B<&oQpx?t{b!L#%6Tcz9
z3N!crG3Vte$#p3YD~~@^_@JV7d`o1=w!C-lr<{`R+FPxE`}=&M(@*!587}e}r`Ga3
zO><n^{KfnF`KjKEVygbBX@ot0_avD!t2*vOt8|`H6yN==b7r2AQ7w75`1x*&Ws`qT
zoV0NAKW_7hS#iK>y#;(Ki8;8k2cJp;JQD<5*^_)K2~Zobtw8dDAwX^6U^j39-EbXp
XG79+o6<~vdbT<H9;Ci*(bP59ijJ`@!

literal 0
HcmV?d00001

diff --git a/OCDB/TRD/Calib/TrkAttach/Run151536_999999999_v0_s1.root b/OCDB/TRD/Calib/TrkAttach/Run151536_999999999_v0_s1.root
new file mode 100644
index 0000000000000000000000000000000000000000..1dc60815eeb9687c63e0fca9e788c748c3b35ab9
GIT binary patch
literal 18039
zcmcF~V{m3c7iKcCZCjJPv2EM7Z5wZFO>EnoOl;e>jY+ceeOtS=|91auS9RS}_g0^K
zpHtQS^f`SWX9ovY5Re&I5D*Yk5D;+>5D>Aee|FTr4*IX}g8xs0*%<@`br1w}L3(%>
z@;Ay}I+n0{vr_NPzo-8%4HO^bKaR46N6U*q0{`pz-`hYyAXLSzY|ZommJW7i0DCJ-
zBRhbRv6-`pn~4>`&BYlYWNW4T?*jp<%4+g500j|I;eVzmK*Y$_${3*PY$N3AYGh&w
zP<FFtVq;=sVK?CVU(LXs(ZGd?!TCS%{(t)if%p%u|6&6HY5NCj_h0`nzJJkV|Gx!|
z@%}%taU=Zic$WXf1~l+L<Kg~+{|8vi-qqRb|Fl80WuaRcK?vPIJVBRoU|7|BzQ1Qo
z(<+!*C6%p2r%@324GN2e%Yl$#z)1>LNoo}_Wm+*&5)9OSb_FH#%oDFsUZiimm+hTT
zftQYpu}&L4g*Ez|15Tg3rf0957p{*Mft?vFlxYA7@?F5IaDM>uA7?{M1Jf9UCAA~>
zGwdoI>@iKPvuyU)=Q8IXC@9b$Uz$NNb$w4aFPnjN=APriBjdA+qJ63cyKO$2=Vxbx
zm>aEEZ8cZS{0>$n7DkZ{L5SZ*3EO$;fUAoS<mJ-~vhSOjgUhy)f?D76>X@*wxT&^?
z72<>dpVA$|*A3G5D*;hrF8|*bgm45bU|w)=D!xsxv8CY5q{uZe=%U2Xf8SDQVodSS
zl%Ca`y@%Dg*?(x4%@Cj=Z-wqeBnpyS`F0Wc=X`BYbb$RHmL=Y6WBB3I-xihwJ@j&4
zC{wp86nS@sgv7R#A9EB3I*>r=8;QqYVfh?(j{ee4(qm<ruo*B1`4E%_Awe~$968$&
zDELzom>)mhW%s6Lhpi2mNP3XQ9X`v6{1ty`R7xi4;N{0dD(BSO<==_GWZS4$LA)5Q
z3Go_*ef}1i@$j1U6lL&Yau*>WdL^(KE_Py%LCf;cZya;uqG|TCzu)JtPo_?ukx-rw
zVWuuWkZ%g_d{x4%&)&2k6g|Dk(EsnFGERqddyo9E)X^2)EgL7XC-pi=uCvc{^`aRh
zVBx{1t*O4h>0>MAq<#kZkd>cyaV8cz6l2HZ#$QB+&i6Me=bgqS4<Mi$cE=r_WKHG6
zfFaV*oQBw9K@0@SAT{H_``M97bC&wCQGi!eRl@Bg8M*z^sX|qqDAc-y!&p?p)8^vs
zwRi;#=yZGcg3kv76VV05?214<bT-kS8HzeOyZzAjMG%}F%2*Z%2mbo#|2=U+Q>aM3
zKT_Pf3?ls`uMZomN+KMACt*^uOJFIzjXdFW!q1w(K%MvG3%=Bo{Fm<BtKjj+(ssTl
z4?1_k0ax;lu>TZ;e=@krD0>$M+SoldbG7`Mt-&3PGFg%_o-+mU&<>>=U<o&AgVkr7
zD@&wN*+R6K4kjrULgzJvnS>pXp<p0}D>1#9E11y_K!DkqW%=fLocxN~hJ3{H#`?zW
ziT?e=`!;cG0BU0k(n^EM<k6;!MC?V>)hL10)YG-KN3|)Unc#~*b)G3fS}tL3`62A^
zxbQIp(LyrYemdRA6XZ@>n;~wy7umx2j`?MZ0a$_jK55nmp)`!~DSopp=<0KRfk$UY
z(&+4&L$*K(gyD5=vihX;ow;2@lhl4a0wC#KXbw%N60%sr@OttFCW7G(bAGzd@P=5A
zc_`_Micyz2C`eai2Ngb&N<^SBDb~BzS7j=oJ0GD!Qr1;R{2UM0`AgQ6U?fuj&5k5k
zIfSsMx!^B(hVgwLBSw^jX|+A-SqKKRL>KyES1eqp6!N6$g@?HDwo}0KLeTlJQ-Jmk
zi3Vyg3eI~oDJzQgmF7`32@em^?D{l489x~^%&pR-7g4sG+}3_V0O1bZyXDIW$rDN;
zV<(E0VPGT%KgF*;3PFc?+4-9B0-ctg#@_f=i<x%5wwJrM9W=xXkEC_X0f)7;IbE_7
z9ok;GbaH)@8c`fBe3NPH2J}6l?5m0pvc?%mJ%%eu=m0I6F(T><%Hr6KZGC2PHE0Q&
z!?KrHepw$_J(3F@bcVy`Z@AlR?4o6yy`r``ngI2yO3V4%Xbj&OUfp7V-w^&xejPa8
zDmTlFY%4dW7=js3F@w3GCiD>a6@<a2k*_vNlX`K6VVRzHudl#|L|@g@?3xr@ASPsL
z(fQ+RI~&f^JOFB&SCY4y$^6g#C?!_)!tG=70GqR>NmldW_9Sxk=(tE~sV{-gm+aU4
z;kEl{GYjlEjuyyK?Fr_%g>t_4oDui)?1bWgE)Er@5j2|uN0^<$78t^AZ+FQf@uJQ-
z^!p}F#Cr*qK&UDpxU%Tp35x!2K^oT^a+(H9USAepsl=;xayT>%S1gu(caRng#Z@Lh
zT87uR@h)N)gKS;9)8cwOL}@YNa&M*|1lnJwvpO6-!%1iM{?Z@}anpGP+b7k)g8sW*
zYBvb%YGqem_h&+Qvv>1&-l#zT;qp-)^XgfG&Bf8lJDhh**!iyyFsIwvjqhhBn-{*j
z+5m-a`tKyqG|ECc4K3VUt}>TR`%kk(g5P=qMiwqw+JZ6`4{uWs<lmH&5_Z~A(uVtP
z_Yy}7MvV8Z0{(_AJn)nTisuJi5PoU!V9s?nA}eMrkWy!W+2zH=rI{^6dLmRZ9h`xD
zGo3i4c_$A`b_IBQ%?T)aYiU%qqqo_=cSDy18(p)P2fazB4~V}O+@{$0iVK!RIpdog
z?tJ^}ZT*Ty@})^tV@`K4k6eQHrR!pvH3$^jZ#xv;3;9Qrd3OEs4KYk-v>|x3ZG^)Y
zi=%{q48^%S3oL@JfvU^nfbp@{5{0ig_GFfXey=0j5=OjY720C=6C@#RE?j+fci4H|
zNAQKEG~p;Cx#_6A;tnIrXA}4nAzacZ;CY!$<q%M_>A<&251lloPWv@Q0>&H`YAc=4
zC?fK}O9~f?*(jxMydpUO3EpbcyX{t*5dL*NOc$#u)F0w3M^C4GwhZUq$YuxPRXQO(
zK*<EWXRaa?eOT($-~U_4$Pv!fTC!CaeJ+q8lih=TVfZc72u!DRS<l?uQj)ZHyyNeW
zbI#z((5#Yz_n{wiqYy6z**en1l#Do0bl{NjdBTd8E|N;L?+Cu)EF@kv!U<TFD06a3
z5!sYFQ#4*wqx0!Jr(+(djQ~{)U;C8cYHTGyUNaz}<zknqX5cuX#lMg7>W(<BVD84s
z^W6s**H2;Ko2T^q>vLvMG-d2Cbj#{s;$dGBONucOcvIdFWN(QZ@4tz$Y+<Nc4!aa(
z)mgYgXgNWS)#m(U_d}n9W`OS@2NgB>`zz&98D;MEx@2hSmH3yS(cqR`G;*xfPPL&K
zXQ6y#D<Wc!AZuLx;yegyVIuyKjD8%+uX=3{i8&+D1C;Ej?M5g5J$7|hlW2!MZ)}A+
z)OZ5Bz>4?M=?fk_qt@cD`DiH-QC@x<1?p<S;+mfY2G5$L_pLp9S2d_W*`x8XX9N-Q
zH&&{gAw*#Jg)zPh0u`x*j-86SSLCQiBM(mAGRgg36WI{%9q?6rjYq_ZgS1*;)on)H
z(7rwMGjhXMk|UBJ7{Fo2CWv}@mI_z>1KIAfCqVEQbfM8oZ+Tj+IBWV#gK7jl(7iH8
zXTp}K&~B4^jX&~$j2O3poDyXYJMA%#k1oEvsbc!$9E<o|aLU9X*#@Y+q$SPC5mhdN
zIbZ_@qC~gC0$C?Fa<}l}%s<av=cv=moUR}PPZ=@o!wIbySkXCOf9$bhSo!UVt<j1(
zS?gQi_t%l0G-JDyhM*Y%<wF;+#-v5M-zwPf;a>ve{Hc+THXbf1v320ps6euxTAM2I
zZUuBT$Q6kb6%%bHk?0`DMfojmfoq#6O#aRGf({q%w0WGbDm(bKbR79PdVKfgn+0`q
zbC@Uj(ikFnK1ktgN*UTrv3Nu25QdDYzoRlk6!Lo{7K`HrJfuu{<3}`~O`lPpk(TN{
zb2~Y;r8K|`P2M1?*sNW$&#QARg+Fg4zvuYAlwhlCFJ`_bPNETt>C@J~e;%BwUK4QH
zD#n8AN)z_n_I{^f`Gvp64qaEXt0@yUPPx(k!w!o2xAWX?P!}IM-yMM?scmSt_)^IP
z9vX4eR{Dw>Gr{~aB!kK9Zr~x6?_P=#on#dLq#3pk0!%Soh^${MvceiuQ}_X-Zvx6Y
zsm>*$Q>Ja!BONC(yO;P;p$^I57*T=;y8sWz>EsveR(#SmRv@Cx6BFmsqY=`xdDBhO
z(Uy^Bec3sSsUx(Hk)zE=aB_NmVEuBchA>&Z1pWf7+oD*;!+n?1emgwBgGNIHweWx{
z1K&=sAU*nY0hny}kzD=dc~wYHv{#>a!K=!c5zno+2$xoW3Fvd2JFR^f^Qv6klf-Y%
z+-!wewu;!grfx^HbNbWCJmjZ*VXEAKB1yWG;j+!+4i5xa6mG9~Frq^Oz%FV7Uf?CY
ze8Iq*RYC#9(tvIuSVikH(GpoVu~`4ig>x?!IY5!q&rb%*Q2M&Jm>uO6*&A|!sS}Y1
zvs=p}^(Pd}SGnWU9tWqZ%bln3_x0Ve*t#3*Z$73uRFIa)Cyh#htDwTeX>>7Wp4Wgq
z|F_E?&@`^b`GEm@4E6H5n+IM|#EI!RnWu`a@RAH3=U6&CW0IdgzP=(6ce{k9_YUM?
zJXN|~Ef4(IJq=PCN}}vQ_|ds?nqnEl9+@BF+PCWV4!norG<e}3V7kYCRO_B(lgWbH
zb(q^Wzzye?tLMzS1kR1F+@*wCFg{)edY|z*odPa=`mMaMy@l7S!N;eCBaExp7qH1J
z7vyZ*E%L&aD<GgWq*jITiX>C2fn3c}y@&>{DR}E*Dwxy1Hp7=!hG?gc{ASaxw7t!s
zzvp?dZYRGGyWbK#`7$LC9#`Xya*&$!Rr2PH=wo<7>kPKymJqYvDcn#SCZT4=PNwuo
zjvq82Frt?fQ{<j7rCq}s3(N19H#N(BiAeC|+dyd*XqgoI0m6QBJX_2C*5z0iZf>z8
z!dLS~(IZsn`MH75=>&4&_x?B4{Sx;7FlQcu<@%JJqA*>HJ~nLR0`nOi^G*tINb$?*
z5^phr*7Pp2wXaDs5HEit3*iMyeTx=B$o*oxquA*=Ddmfi|2fG4s}_-nmGjoI5^dz!
zaC@QQf5qtehW}ZKq+qoBYS$5sQD!JucXoQmsdg%W!H3hGPYC6}K&WsjLFmALNL2VQ
zw`inC7CEeFRvL@c5-oRmnb1!KNgi4mK}>o0BUyK^{lH{_jZfv3*J2muf0L=QaD+HJ
z@^F2pIEryfzzN%kwm0i?r9K|@%k~swag5;HC2*TXO1`iQJb&b<`pdrLXm;XT_vf}h
zohw#=pRqg8h2Lnk->YqVrVZ;?Tdjv%2N!sY@kDxzcwIPv+hb&Pd|U8NWeYF*CR$L4
zo{cQVip_V^bhN^c5#RaDCuQOySLnEA2%sk)BUl>jy<PFtouDsx(%JI_QS9-7E@$r(
z=p%qPP}6ddKrr=f?N<(6?5=2*Gkmun*kSm;IQ&>iY$TlG0(n6tWJ-ZefyBFrhl`-2
z5!jvMPf+p!PKu7t+o2#tRGJ>E;mivk1XH-XP^7%ah&virzs|#mPuVQI#6=X2n%J2x
z#t+}SJX(3r2^rF;Tu}M+%(%OJe1&2U`T2BwqcAsFAYi(fFh9wf$7t(NrBsEG7us3`
zKjFHTS&%p=ZZ>-z{OdX0tE0A3k^g~>GMn^Ai4y(GI|5I^D2t}<i4XwYo|?yN*$~bb
zb$`Rog$d(+8;^1)26D89Zb=8X8QYnkN_RWY3yJ@oe?|QkJw*0w*T`m!d9*)@<>nXG
zF81`$LFDJ&zbSn#MeW%9@&Q8{`N_ckGvb95rd8)XEQ}tAY5(9E_(EMsvUf*6FANml
zj-m*yAhT{0C35j4C~UAE)Dg%OTAX{PJ4y2sTV!xx56|%iYVWkr*jK_U{pN0GPvzLx
z2TaHxmSQoe4ij2xY>KIF1g%aJ>Z(uWnQTU4rJBf_$1(b^0%WHhkz{2UWSSGFBQniM
zy6O63P?%JWZtPZQMb{mRSbITUT#0RVFS0`<^M=%Q$C4xtFYKQD(j1VrBDD?+3j|0d
z7O7#LSjfgB{Y(o?K2W$sS&EsQaZ$`z8i^1~N$&T#p9gDplpPTNH(f)iZG22wBo3<I
zLD*I4C4$Gl@Iwt4iQ>fy2;r~^H!E%z71iM&Pl+ekQxmN8I=_tj4c4H<4((rBlhBLi
z&;Wm0WeqzOo!2*qlUT>8I%dcHEOfVKuyw2vQ{_P1XN)E(WcPZkp@SA)G#A;Hd9`7R
zcf|XcZ)BfwW0Ebg&yudr36csHE7D#;mIVa(F1Ha1G>90ZOKMrfK+{%7I*1=3JCPTN
zc-hdC<|OuH6bys_AqNwt#py*0FDtd^=YxTIIV%Ws9D7vZ`fS`@Z^X`);VU5}lfXws
zdiuK|R8jZsUbg-nf}Lw!-tDx2Cr=&!!#m9HAB2F$)G#$^TZ--@>PL4w1Bjpv7hQ)K
zZDGafOf92j%v<mJlC#<os&*k&&s2#e2njC8&)L74vywak2N|Y>=Dm5bjtMh}JgoKV
z^u8$JTJL&JyUYfWdIltyFYV;;3ljz%ER!A>$MXq#(g3DlJ>83&57UT2-{9J{6wC}$
zPPnHZ9D^w@%meAAFL60|Yx$zJYHfBKl&_CZ<1Pf|BBS1CD22%@2v%KH8EF7I@Sj5*
z_2}><B(<et)O|$!w^mYg+QTBU{k>C2)9(9W!E+=ODo`ZOFMR!p*}TBSLz7xPWe5k`
zLx(oG!6^jprr*LZMa3Tw`2&SzGG26K<Mon-Sz_jYAmHg*`1rV)qQibsVE1N{8`tIK
z08n*L;(E%F5G0Qfhl361*Dz-pLLscZBH(k0jTE?lr#v)77{3@ipP*)HGXXB%dmod#
z>LrX|>=!hAm~Nh;_P}_yjonZ9hqh+50=b;gmO{j$kI`J|<HFap->+tyDLLfhyI{`Q
zJzu$T_`?8`+Jqvoo=fE%HzMDm&Z7a-YaW*^*rpD*v5j7Opl35g5@Ba$Led07e5OHA
z<U5(ZLy0F)k>IruGDEq7-3S(K$w$&azWH&2=5pLZ1c9oaOwxgmuy@)yuciO;wa(Yp
z!Blj!O%Jt+Yu6?YZi((+d;V1<1914^cIx4TT}o4HS3NKptpiW_Lq~k|_S5@Wm9Vz<
zm8{j|48mxN@aeQS?SzJOKZDlHML>9&TU#0;(HK;KeKeT&SuAoS#ic&AaFFUzlo|C0
z<=~5z0e_D?ZX?tLVjU9fFNIw>4FoIvo_)mRHKDibtm}bS_uLtY)P>u|GD;VhL~(3H
zyYlUhQ%3l{usIW1^iwm}j7=pj_#C#zaar!Ne7nbR71M8JOAHk=rS!ghJCD}noG$Wa
zN`9qkb=aWRWDIZ{FfGf9=;k^u7o2s3Q!k`m!%V#YX63$;-E}PP#x*BY^$cgze#tjU
z4bP$x3uOJznp1NZ`^S_P4Z)&1HHaNh$j)8b7r+EaQ$K1|6eg80<+DDWDQ79gib``&
z>I4urt-OH^waS8hdo`}g4j45HE7>o4sC@s(#gE$fR7Jnvl;gLskS*vG@|Z>aUf_hJ
zkUESnJ(A*~*dOF0?ogub%Z-{Nffy*|Z+BRj`hws!V7%VjGKS%7<_2qBIgDPZTCS-P
zNV2!Ea*~qhf6CoVJ7)$K7)g=w(p*d-Vty)C`y^?Ah)9Y5JRWf|+)e%!eWSN2Rl=Yk
zkn#onPPeA+kT^5&hsW+)hc9>e1r~kccg@JgJy|@*$zKz)ubB~jV&j0HEgZfm#(oQ%
z*>C>PP7HOcWa6g6De9jZ%p>;wvfmbGQ9wcLl{#~sF(*)Qh8?Ji7$NMmYMII1Pw<c#
zSpGT|jPO5Vd<y1RTfu+2P{3Nx!5BT~-j|z5Vc}D{-MH%7ks)J3<@C|lzQHcPUW$Gy
zs}E+eO8j~3KSyfd7iQcr4%m<S0reQYG>VnpARN6W-V*iMp!o4xKnyj6M!P!|B^+y^
zXx-&*=o!Y3?sDxbLme5|zeyEA<+`5lnf|hUe8xnL%PGW-WseYcCLQ91a#QdKbCmjV
zKi*+v{%kUN{L5Vcj75A(f>4JelQ^wOg5je(Rj7J}zv>CDli0k%0?<qYg-mfN?-w~q
zDB68E=4OaMFo~O?FW;mZPV+E0%QjfY$kddl;W>wcEifkWb9pwLDKZo?q+A{wE&X<Z
z*)3?!rm<faw91ft&+)q@o9_X%$=*ZYL16?wKX1nF>O=&}=(KNYEegnfDu6@kAvyxT
zKK-i-p0i%);C#@&<uz3iuQZcuOb=4+^6QZ(H#aD%B3PY)o#M68;9xEJ#7S}3Wxody
zO{sdR+urA)1eU_(z2(C;02dZQMS%(62{&PfVlr+#LPFv=S$vNq8%YQ)F0ydRM@dq2
zkfPaXMhCbWbIc^{39|;UlPQGQ>xN!-v6u}2%MD&X%ZNrutS{|9LdHF}bXzxWxx95m
z?P6!QEe|5fSGB}WzVJomw5Uuz2a)GiMLf)Dx5QYdZDi&k0HROFk9^cVgi0fYslL+6
zU!eK|c9W|aV%t^fUI^H5>82mAa|o2Hg{~fh#r?M9nD2XC+ra#;i9G)ds6_A9tLgck
zSn6mUhKusiz1%6X9;#mWL^ijJ@s?#$>9qaB&tO|oiKW-mO4Rt#u8<Y>uvp?K)+UHK
zaD(47wYRPJKBA6Wg}TE+S$Ov50sHJLkM9-jS~r`=?)W3J(}xf9Ur3U|#tXu@AxKI2
zoQ)m1&=Ra1Pg&@G;-Wmys)ySJ4oB%|tU^EF3gu0Kj^Owri_Z`AB@k1b_MllY#SFyf
z2M<lXNp#3;+~z?vdW?nony3waMEVLIXUn}6`~;#s|CYQEr8nPxyH)R_F5UtWMih*~
z^X*U6QcpjuEX5eZLJ^V@5|kIS`~9+790QTT9Mqs0LV$3wO`*OP6(mZ9(<y&z5e$kC
zCOYR0f4pQio|V6i81kA-q?9EuLE?*rQ<x{cAo*Q&^Tom^aB@Wbv@<8>w6>zvo-h8U
zkjsnnF5t;=U*O)+_jh1Kbgy0~Lt`fveynpko~~I~paj%B;t*4ihM133iy?EnPkK(5
zqAD22&kk&9SM{KWzUg78_aQ@O`mQ6HHdbtBtc1Ao-+^$5x_>ucJ{d`tk<W}uaxy;0
zQ?j;^sF)_Si%^RgO)k%GZE9h1CjvFA^u^?>rm0S~xl$mOQoRZqfE@$ICU}dlT8?Nq
zd<Cx6l*@0ynCRp*vqZj!U>&Vw>NwpT#PmaMLs4Zp3H^%|e~%l<$bWQIz9qQ>6>961
zNr)mcJThen<i&(+_siF8XXt=NavT}ecYjhPUDA@1ctt60sQrv}IR1BHtr2*gGsomV
zRrA;JL#^beDc@kf{0yPSJmCQ0(=yUxx$Vd`w6tY4J<CuxKRFdbZf<|Kbt`>yuSmR?
zCd5;N29=S$>U51!ynx9x9M`in-OUgKP5V0p_o?xUZS0OrtCNZSSaAB;j^EI9wl{;S
z`Yf;_&MBD(>_km#p|5zh3<DG?f5`pB!XK52wfzQu{X!cewFNDPnQFD*@j!wrQB{iC
zlhltaU=shnyyK~<u{Yr?nXuD{DfqFa)F74}N60yUJFy22QFy*l@*MT=Uz@*6HmKkw
zvB+%)Cf~s|=ExWOtPB?t;Tg0q+*(a7|6yv!{qO6&f0^}@lgy}M=w|+u8O&6nB7`+s
z?dI<)D+CCIoKDU$*2~o91fu8(xStUgHkFiWL<8NKq51+^ffTe27x4HxF?wB^T~yDd
zY~Ya*oqia-us1vcleIfV&_<8`_#5dTg3>_82_+6p&#b~8YdkB|dE8M_CFP3~m36sa
zko>BKtU)xN1ooTmO^jo`&nhD&=-56!)n=#OkTFq<Cf!NejH8n|srQ|@@9pO@EsJ3V
z6S<aHg`VK+9s+Wy2W}zilss`CZ=#kHzS9E+I&kz?69gQ12T+7(i09{cd<A~px5GgX
zCbRe}gM3aje)hYWp9G!{#CP?*7LR)rWdg>ZVmxh^^kzey7ZG*G-r|NWRIeWkh;~C!
zLm1s<>0w551&E214CJGle)-0uE_moWc9&K*SRs`k#U&9)-9gqqH{tMX0p%JTs~H~h
zKdjr%?z)b=gbAj!Ad@G<=0_&N$+$tSi_z&(F1G$d)z1w2E@mrcqKB`(SHjH1{U7>e
zEtz_Z(RRGA+uYzNfyG+lN+CEePV{a;FN|hgh<d!I;Arf7^wR4}jJ=ol&IlW#*j{Z{
zfn`!?@D6qSt$S(^@J~&7zF1rkB88WyCf<NYqUd+zWgi?!RW^(0G~b~=S<w+_--!`f
zE?=MULc*;o$9ipBKTtHpPJW|)#{e5%EY_v|0$VnEsd$#K?=o$oH{Y8f8!)hL*fzY;
z=4k$%uff=f4_Zt}h=MsAzPEft6_3C&x@`hz;fP##$H;JOoA}`<iKJRq&ewL?V*hj%
zXUrR=WS6}=%Dxa8z6iWHZcRuMYV=$dJm@E+XYukpYJqQV96CcPxeiJ~EE@w&dYXUp
zSmBoD>Ukpzp@sGCF6!FzSy;A=YF>~&fZ|#oIUhCl7xd<(tZ5iSx@BuCV0#dTb3bM0
z@IIm_zL6O3v>@CoR!XOLThwfBGO~OynY+k!<ptps0TbLC%YPIKw3eTk<3knlxD}eR
z58*RORal1C@&$7?(!8r9Fj?-MixfDqg#$AgVu}!*NV!=zV!}4Xid$mc@sI^bCuL2|
zU`5ILDYmX!^fMW=x*M#oKEX5Cq&@GH=;2X9&G>Uw-l#Chzg;HZ0j_Y|98ba{ZqG5i
zRAi^T7GKih5FsSJ^LSA0=!V0KHY?UoY(U4`CuPjJ6Gmz)B^|#E_?$@tk6$B@LNbEs
z*~Vy)as-Le5bY5*y~&b@@p`Zb?~fte@_6}163U<fjAV-*hePfa=6e}aLj#fDdwa*(
z+5&1Z<xVnw#9RkYjuKl#rm6{p+K_*mux^9GDQL9Ncy5D&oWBh8<2D5sogOg3B>Y9?
zk|Uekn}ino4`vkfr^FWF{WadqyX6_8oh`baX+@`Wus=&Kk*}{{kNk$y$*i3=DC;)T
zfbR{&<o?igd0E(q-9_2oNXnOGUhQ|8jUv3Zd4}F5fk8P_Y8l8IdY}4IP^CoR_>QlC
zVNnz)+s&2JNs=ml&G6kE0RNKrYpo_lQIFWRch|Ze1R1=$d+pMTo^{udk{(A<N6o)Q
zv#uc@^{D(?%m;wFA8SEPI)qC6(*8J4ZGw-PRe5#IX1V_dH+wNV6R0bRqRJo6eN8z_
z3OgC!9VLlL2HTHq?7sSn$`IS}PY=gfx~Gt=?^%m(Bt#z4*a6pXUgvIp#4TwIW0q!l
zzM1{nFl&gU6n(<-`g{}yE##LCb-{cBzVG^F`O#b-$<{7Cl1GzTdm3K-S?YLzX=_jq
zNE<z~W2_aym@PHy54D*~aFD>4-6|!Z*Ym(I7vn;WCet6d94T{9V=yN{RWbS`{KORP
zya#!1o*+#s#e9pYfF{*5#c0}g2$3Au^>3WWBTc%cXe%exg=^}~K{;j+Hr?puH{UY@
zS_*L1lfn_{2yt6}Z*VR2?961I`wCdRuQ0eR{$uJzA}`MC9X|Z8#5aL!Wr@(h0<=v9
zQ-h$kws!b2CeUny?VMfa%y69`MR|Wi&}W811FSUZpa_#`nW91bdTmV>S6*%iSIFUQ
zNxN?7)<k7Gissxvuc5dzr=b!<uha;KL<7iDqK<>*mVTW2m(H(uxjwJiR)Gp+Plbb(
zh!|+O?&_G-=gL_i1AT&mKGKto#5O#DCJ4Y64lDL4tGH@U(Vi%R?7V+DTl@#USYp;d
zmQ_9YP@MgJj9xz>VFJD9PL6Z1f};I$y%W;-V{oeTIYtQIioE6%li!uo^zIA7p@-zy
zpEd9AM=Yo6@o$l*y!V+;j7lPkH6q2&LGppWTj2uHp4p)v&y5;yVPy6B<wfZXYFk4;
zuFk;?WN(dTQk2w2yr?d(!z=96PwET_H#TCecOXqTNX?g4(|;K=T&zI0l-l=g5@-q7
z%4@~x;XVfE7>qaYh&z$SnCxfF-BXH1_oQY>T*hLTtdFdkG?*;KdCmXkH~ltP7}&&L
zI+Jjp9j<I0<#<P9b>u{gr+gJ`DSY~>p{goSrHU?6=fsAcR_`LPm<@?{m|oOHNe-<h
zK#G_5Bno+SRC-44Ad!Tjps~ZQ4GNtnXr)S{6BI5u-#ci>Qc&FTRSK|U8-&f!-FJ@+
z!Hx`EII?zdIz{|M+zPspDt+M_`&ti(d%wxON|$dcKrIDoh9S)Mv|p#JtXf&hCdZ#{
z8zNBL`_~#gDD+PNEqZKG(;HsLwg!7QS~WL2e=U6ek<&D6DDGK#aD~=&>TW$Z<O-!K
z#Ye2*2rRD9x(s=NlGMKGe(Rw}Xhc9v>hw+})+4``k9|!vAYvt*)c5_{9e)BMJ7mYp
zTJsZZ;UZ_)AUh}cE(8_5fbPa{&{M1aX<iunaTR06--X3jjsZ1}vlplH-~K#M-QV)%
z2)R*_^GD6ijHAwVX~Z+e&=t2$RQ^-k;A{=xCcmXg^85LT(lT<{orrpsi=+Nu8N9lK
za$VsyV-Rv<)h$faqG4OkZJjWUK{I{nE}55QPqT_*n4>%Z?!)<nf@5fdP)03g5uPbA
za8<{*C#UIAas2lEbYx@8lWlUJT0nu1I@)`+)G!^PWr%-ygC<+?UAM{S_dU33S&p*r
z!VqF=Y_2)pYS#&<Ca^K$`!pt;h*Tc+N;WKxv``-8X<lAJ)bj+LM6ZQF!1H)qKXPw~
zp~yj6Lw5~C#nN2h^W21S&vWI^zCUuz`|w*0O=37oXB!2775UwJ%hHC!#Dd{%VsGUf
zJ3J&U3<Rl7a=t*7BQwi3K9tSUO=QcXk8#DPT@UV)M6T`h+mQ5MhZpE!9!!VX%_Xd0
zqz9wc)^c<yvB4k8vai=Cz-#rJq^q-G$-wOcBWfEOZEPAhPe(=A2FsDF!e6pv;kew}
z1bT5mQYXdLacaoKtt8QqjXbiXzy&yYP;pUOYvW8+-CHEacYwzSg>aWGPXONq^`k9q
zjO^MDoXNK8RTLGK;qn}qcN`R%Jldp!5D8c!y4#fwt42T2QS<Gt+Vym&+Hro<{Kxvm
z;oAF9OTGwuv#kj_hIX4FE}Wcrz8=k9e`l!<C*|P|BZ<8Ki0R6Xk@2&M?C~%0pM&4(
za8{(os*C8!etw`#?<^&n{;J}xH-0?>SxBUdr2g`saZuU~V39f)ghJUyT>=azT4AeV
z?$8d0d_fsA&pw_vvf}XO9+nLs;V3=yFV6b<QlcD<VXyjbYXW`yzZd^zLVtGha16~d
z{F&X8$@0LB7ZSJ`JizYWnd(fybAUC`CpoV1CF2mV3v${3h|6rO^w&#ml@9SJ=L`HR
z^9Ga%xL&BmpR`yAt&op9`&Yxu*-v;hFTCwyUw-h>yZem@UhqH$T-Awx4v0)~Vfgn)
zASKnS#M8Pe3Zqsvit668CV+bx>NcGs7ed(dOvE@zuy%UHn7!Qo$UcP`7)XA8^&Af@
zBq#LEK9cL1Ktj9`S@apt;bW^qVY^)iIx%gp(%RnaZ*efAEX_Pl_fu#fw^(}j2<P%J
zKtE?c)JuIrdH#*>K|`VJg!D^_@ix@vJber8H_!0)Vgqi{f+j<{@oCVyqF4=C>~@?*
zNKS<i6l>%Fohy!;Js*$@^LdxWbdQ1d_nF~|8;GK11?AIOiV>%m8{Ftm7AQS~BhO&I
zGAKixLe7PU#DxO=`iYvwu?MD)Jo4WKdU<Pkod14&@ZX=qNm_OFs3$!~PpUuEdR<q6
zm->mGs3hAW()_G+myMJLbmth+)U6_#dwzV5B4m+hhRDiU&e0DkPYgPgT;Kqn-~9<=
z@c<O+voTPsu7@9B>PqVd?bv77_OrBLD|JAmvgG)oED3*>H@u3Aws1k}ePH<eb${`4
zXRY!gi_OZ0-*Ko+jbm&vzIF1~un_@r(^_=iNf2JRlh%?gZrS}rN9=FRe*Rk?x0Uv%
zyb%&NV+{o8nZ;8uR93Q+4Tu{9{|vFUsUcgXh&Yb|36%(CvURc)gbyP=R1EZ~cB%ZF
z<MgmBMIKE`PwB@nN9-7KBcJPXBf|yqzB%#wlk!e8lMt;0V)<8MJ8tjlEf;En70p0h
z=2f_EN$z;B?G}TOeLU=HuTFLo{TD|F8gO^>+%j``E|@l7q-~|gXnp`9j>UVrB6lmg
z3#7^BV6Xj_5gl;xAQZBaf25hs?b^NH4Q)-bJIeLQHQc<ie}LwF0=9y4q3QR#&ZV<U
z@-vM*=HAQx%+^1Qihn@>3Z$Lb(>n?frJMas=&3aK79@|J5OjP=A{eX{<g7T0iOO`w
zUTIL2Y3(K6t6BsbEyK*u#w-$6b$Uk)`|y*hb@lL8o`raI@^I~3h`Xzh46y4-Ozagi
z|0m}vgu;GsE}flVbA=9hx{sKxa2h7|KChF|=P~)Q_^@)e|6oXk*4E3M`mvI}RvPn0
z;K;;Nmv!^ZcnZsaXEZShi6y;F>lFsmy_UNfB=@0i+?cdfMdrcQSAq|d8STMtf*8Hb
z*T&-p-i1+UTJ^iek6-W!V+FI+!ZRP?s8!PV^`$`hxJGKjIiG!#qGw1AuGiHSF!i_z
zGsGSC;;=|19F)P5LnR0A?AI)RQQsT<6Q%d)SNgTkI?M9tzJ?5mTVhj`_c{o9kGySy
z>J}V~kc2DI)9?V4US@6iw8OnK17YOz4dGe3d=<Yt>IZ+k%j^u1A-7Hfh@J_=aLvby
za7z_?2{Mjg__zM__0$VVzM9=nG7)6B+oQ7Mh3*E)T%I`{QUa1+FXWk{iZGEUF<Nx&
zaX+5BjmOf`nLSyu)?HR=Uj)xDQa;eSZhHjg9KDw6M|eTtP&gVTd7rjGn)q8l;o0L2
z8{N9Lh|{TJ@=$JRX1c(NCP<qh?A!Wh?;8<hkO{VVwAmhZy)aIlSXNgyoC!Kl-D`E!
zTCvs<qS@@$4TMMju;(ISe9e*!C~9MS!g7!EpMR{y3C$!Hw-=BBjv>?5c7mShsFp+5
z$mj_SY;4OnH?JpGkTeBc)=0+kX{)8mPxq@0nblGJytW!P-G3TMrIxr<p&5eFtq2G(
z!9Wz%G+yl=C@!rXn{?_D(Oy{G*W~&Ky(5QHmw2<6V~tf7S4D@9BM2?)aqB7KuVIMM
zmQF&|Q%A^dQ`1qCDHr}xifq(j+{pVQwA^p_Tp9LguEmw1*GnbvxUP}$<=D0G*bOP0
z{Ia&)`WSvUrkEZ<#?0tXW@=K}!yfZkan<2~2|o}V4}bW{i0$6O=~HbBFODdXOlOW9
zxSb2Oq#$!wfqP5}Eo&)oFX=K}_gzRFIQEIYAG64ZZBbjG>0J8f5f!Sd{7^_2eU=%C
zDZY@P^AoX)TXt5uqb(lnDUXanZ#%J&1KXITdRwT@2ZrgvhJ%J6M_9J>BQLwLR3s&d
zO(E}h`vX$_TOCX=q_J#AN*5{XGcAzdw-yWmskl}ij;K**s6lxA9p^f7B|Blm8}I%q
zT`5!6D|+cOum-mBBAf&^FVm$8`tj25{L;q*^wAnJ8pn@VL?_3jL3pC=jQV&#(m&>p
z?!xhJAgHlth9&2?2gckHtwGIxtMFi(&$Y&wfr#lFRB*&wDT=^~*kyD!ytWP-=|QOz
zrrr2<mJhcK7QF@=e=p-Zm$Aayse982weg!KnqxpTsIepmVso>=1D(;V$}HebjT72C
zxEk&Xi&hw~$862BU*?o8_j;d0D{_^v-fq6PARbfQyTkqJr_g38VEJmpBIH!V9z6!7
z0~%~e$YJe+bkBxRMNsqE!m>hn>BRQApNX;Y`zq!ufU};X;BznX_;2lZ1=O>(^OX?U
znWXW`cDIPZTS6#o02<1zxU^kYarHS7$uDE}R@Zh2CA!mioM9i*%JfV+D<_gd{`HsT
z<4M-S;W2Q<(_4SNui~e+vTEd~5^t|c8`^Eoii_cCPhaBnl{t%Af;!|UUa@-G2fpn<
zsaNfw_ECYpqkzlUFW&k*JI67+diD^Urw9%@)1AJ>E~mQ8lc3j#sRkcR-71tZ4eq7-
z0d%F+j$DVDc)=8N5>R{kzXZzyd`eq-@u%0~p#~Y}VvkF>S0yFaKq;?y_0S0v(VKGi
z6o*F;rSV74U-7AA!_&q(%Mw9A{^!0SD*m6G_8<NuR|1w&O#)6?Hl2Ub>2b2B{mt2^
zYN;Oh;aZ=`BB9bgiaoqZVZ4_kT+xPsX}-MN`2!$x@Lmi6(!Bw-bgxigTQspXeb0z}
zn%t|S8e9X5if2bw@>Ns*lNGS@V3x#q7`xD#qUzcpE!dH2gT{j!?q=r+xK1zleQ_Wb
zn(4V<H=Gv7nA*j1zQ~1?vHg@ZApjD9Qs*Ep;*E5-SJcz6Y*Kc27a$8C*)|4z)HIoJ
ztB@CM^$ZBNGWbav|3Vx~@e}yB^)Qmsg18WTEBZ)JHJ184bAbJY6Zzm5C12G>G_qL`
zBhtcr>MH}6%&XrpC+glt$}jN9@9&qeSKWNl9g8RQx8HAqW3P6uI}RrAzc+w8%es%>
zbyCDMFLnVLEYS{+8^YvSr&?q`_TJiXDDJ$3JjTWdXPf2|7aHc@TF+|))F}>LKC?2#
zxXfan7|t_JwX?AjmYO*IzMA*i?3knG*(0OvQMVfHeitO>HWUO^RMeIP_Ev=TIk?2a
z{^eggH8%%+dZ48A?)>rEU2L8;3I20B3mQ0AY?)M;bq94aNHsyWzcB!~-RT8<BSGf9
zb^cdK@=x)=&-!1*gW!J@5C4&ufS5w4%38U&!mGGCn;F@eIZN4_JD`d>xEb4e5!+hX
z+nAaDySTXilRx;Ng0S-clRsSP*x2H)xbgE$==GA%C6$-ucvcU#iq4@*by@vQia3m7
zGKy3kO|6MAvvKzb=UKj<wl-HM{~RB-N?MkChn_8!(?*w^gVIK){)LMGyND3l8^~@M
z{Y4E6POX^FI8Q{sP6V<3L+-cB_2u=Xb^4|mdJCv`q}ToQ=e5_i%S?~WY$lDgQm`T`
ze*r!d^<C^up{WMg4D&=<n`VPv|3=+hy%1~M`MnUc6H!T%Ou172_xrR<)+>9Ew!?Va
zM8}Fb-^}LWjee2oD@0p-(qXcYYV@DU%-$gfgA(#-EGv)`iC>1ItQlfe>inzGzvT4|
zwDhOXYrFploH*KQK^~iH%oIYtOA5hbyK<+ph5Ki+W*eZN#)}uIurG9I6r9QbsRWO^
z(j-Kg+UJZZvcWMG|ItU>sb4m@(5_F4aDwq?u87l$=K`F<MgtoC7VIeqloWMlF|v05
zJrsB$*yIJQ`Xi+M*41LwJaNAsdu}ky2GKy>u}`GS%8w(O7sn*ynYUFPRc<@4xYs%B
z)CnUhI`BLc<<ewMZ;lvIym-hp-MLlH6<A~=WuRLPUIz&|r=8PGT3eSC{$M7=ODHD0
z<|4mKG|Cl(_$LLnlPE)IYGp;QPEgc@(>ur*(ta5O%mEe*iCw9$ui(-aIopL6)9>EK
zqBoR1GbFlv20z>G#1IuDBLun2m7n(6WXfF0V(^&#WGm?+HpQ{*<PG!N5a)KTrTA~*
z1yhLkpe6#lOj`f#8z#bzTTFVku%asVh3^IrWJaoIeQkR{Iou0YwsIu@-N>7aCFW*X
zfRLYK9>3*D?=8Z2z3{aG8G-^J@|+j<gujhulc>q5EbG>X2U#dykTe)tG0)kM@Y7T8
z?v$m!Xh6<5zksZ|;6&~QGJ7Gw|Du411RC2wcNnX>D9`ppSwzABt|R5(+%fdhk17iR
zZH_^wh1Pey>>4B4wM@C5U9rSNjyBX%DFW}|34zsUxfZ!bK|Uy07w8hk%4-%@nPEbN
zP@(On-zJX7*iO~o-X1$4xkj-9PmJVU$8ev&gX+~Bn_!JE`*lHPAQ4jCkSAfcysFSL
zA-oCrGK%Q&|K~-ob79Tud48`%o6Ewh(?&#(Y4%)wjMhZ+K9V}tM8+Z=wQZd5JGt}x
z?PLDoQ!c28uHuqw1j+JYX8Vn;IbHZDPr~f_Rl(!p6_R#SJdBpYskP)lVAD1k_-v|C
zlQo-<8tybwvl(CZOKnWYyp~)Z;9EB{JyrDFTEbL{mVBS?CN`c2OM}E}WQ1;4|JY`*
z=9h~%qOrpxSG0yHL}b5pZa$5~%R+R@1<vl>Yx@2(%T6Ywl-xW0b*Pj1F%RWeoTyb)
z|GdK2)0@UDb9-(eNp#7{de~z9Bb#cw-DWc9L`9Y72suKah8U4NYb3E{>Yy4F9dxZ`
za)_}i-0im(eh;Juhe(x^Bn$BV$FnDXOpR_GD?*tOv!mhVnD7x8c~hkX7`KLob~*}-
zSRs9!#=vM2>K~NmavV{jN;)*|Sq6e6<Jr@V@}OJ0$28XVHYt#w%$a`vzj|3b4v3UE
zVri?&VSVuZ3VN7~GC`3-VB4JZET-}s>knm{#DwW+{qshH?-W%=R_s<5<Z>XCp+jsw
zIR)g-QDbuZBpx_he<EPWWA9HC!_`S*r9k;}u7{aRvJEpO*h506xwQ11PiSJFiR+<c
zVUc5y%5Me@ihOeBEmHq>;cRW%2n@^FtIl?nZ6hPQ+jcZv!}%5dT!4WpyiMv~8AWb%
z`6?o{kKLSpuU>JOXItShUGX|n#zIBu_7?B&OX(sk;C1(SPd_m*+cJ%Rtl%^28ZsI$
z2bQ+f4R`Zz;x_%N9Rkm41!UsI{eDTB#_BPA^O~uEL|x&R5e*3YH`C25?h-+HY8MEx
zR_F~G;3ki5I3_WII_~-YbYny{h~rv1i^a@{4{;jLWNtUweW>5w-knucRKREu$^(s#
zWff+y<N&vB8CGN3E?VK*p1boM@U{2kLdQxDp#-;5EVzVwXO>Xm9bl+`0#x^8paa@=
zqe=d{1^p&{95KEl`LKmR?7M$y-NzGU`x_1n?*~htHlsdRMwW)I!1f(9Eeuhx9S4rJ
z8qn3nTWt0HJgaJ9MiRP6Yf?|4t_0L;L#tCx_1UUH`;d;W2tlwv7{X8e)*3Jco8M$&
z?9Nrp#M-s_f&bu5^6#_@yPnh~m`cfSP77wQkmY-%&=rw}o5LK@ZSaDqVWFNMLxLA+
z;^)e^#6Xq*#LuOlOc~Z=2tVYnAEnFz`SM~xl^1#-TWM`Y&q)`lSjZ}!+nhVUC{Emx
z!K}8@(m#{liU@q)ZNrk?p)pyrQ%S!Sv$r1XDI9IKLzo-pz%&4y{w<+2lKP)uac602
z#y)%@CRw*~xMLva>N&$FPBjVPE=uTcCQ7@_2xxkFT&0FaOpN;M*>UWPNR0&I6?13m
zw#?j?0FBHEpXdza4x3yUU9Uf^)BCI8!{%Q%SlqWC!56AOWS4~=b&B`88Hxd-QbY|{
z))<*gZ)Jo$X#Rh)0R<wHU|n|(pY5@%8TRV1-cxP7h%31Rf7rKMip>f+j7AOVri(8O
zd8}uaY&Y&X6Sgt7FM59g%0OpzsPC4r@yLrjp|Di{63Wmbh6hU2iF-TAo$MK@%t-y7
z4Wq%=N#fbWCZkkjAt;gDw+T`Ht(uO^`(y_5XM=JGKYexssb|Jnf}i)OfI0YfkwEKO
zue?wBGYF<XvyqAaiV>O!q+w;?uE3^}#SY-4+gZle{XoXK<U`(h)tB1{NX;1s4Q=mK
z(@f4Rx$>+>7X+=S1Fcih|A~3blF{K*uU<onw-7jB8X#c_j30#2r<x*GzEjefbdG#u
z<By=B3lUOJ`usx{f8^K6$MjsjVKbm@SR2dXg?0DTcdlm`rbyFi{-ONL#b_H4V#71E
z#*QgHLHcJ>Xs2IMEsz`BZUWr42{*3R?9`Yp&7Xt(f?yuLkYmX1H?rGCSN$)QR%UD0
zY<IdclnD+<8E}dWC9*4$KT+}Mp+uUL$D8a_)IVp&nv^aFNK|-rdD0%UFVe_2Ya)tP
zA*~=2sd^%(kR?hu;PBbO<IN}gyhU%B_C}(mw$?qR^StjLA1BYK@)pZ1a=rgXmX)ri
zKeO~2k@MZU<x8lYlPSf)3(s;5Li%vFj5enoqA<Q%+qCf|`E5T&|J?$oEsI-Djo+r?
z=Wdh?G3SI}u=jS`KLh0Jah3)J2>Z341qB3d!%ArRb;R6{5qndm9A6ag(I9`B75Yg$
zsyW(o@&5q*Hl6Te+Zl7UvM#)D_g5a7v@2u|$r@)x`qPsoO>-$G<b4<QxUbrbMaHfY
z&RTRqAUX}ZAn2I*VBTnNdcRw&syWbNW5s~CMxw83%tXalNlqDPBG*^Cie17%PA9r<
zj3e%90xEQv5t|Ht%L8(L*kHm^s24@fHqomvRcF%KZg(8u7{u~z?4E;-4M*6}I>J4T
zvOpM@;vS=lh+~;(?)aP*v9Y=sjldglkgcW#7j;0<0d(Fi(I1XNwh<R7cKa1S=$iP6
zhgpt6Hlh<_O0Rm>&yUiEu>}-1HNUoDe5frPO6$%j$X9R;2vB)3@$9(>C@ViQ-E4%2
zefoFj+vA70VvNJ$5Z}B*v=8s%V9*7gJ-70D+2?T^3F+7F!z;*_(}WmJXDqB9T8DLy
zWuE9t)@7VMao8jp`gsAX@-WRKg;Oa*j)Jc&_fDFM2AcT$QhMySU>A%kk%*fsAJgjf
zx)!G_9=cYiy$Y&=wH=E&*I9M5w`psGX;E-x4*h1n%%gb@;un7{s1>IHmQHMi`a0&|
zfW}pcP3LA~#k!2Q5UQqa<)@q}swNuzm?7qum|nwed}x?N(fv;9XQ+jGv<2VHlkfZJ
z$IUlnUth-H%pTtHfkMjzQzY!9M0x$~96Enx>f@?mqPko>Twq+Ja|@8O^t%P%ZkGq6
z^$Ka_$$aDLA@y-t0@e4MRw4h>#)uf6L4TfJOy@T*KDKo;ChEKi2ZP7u@KM(?!Y%((
zYG;&3FZqwypG*&=akp9iFMA#W#F6gG&}5n^6~bNAcyrT8lf@}mVP64e(8jT8;cbiX
z)evDuLW!?5RPM|R%^t#>n~LX9+Co*0IG>poZkGY8yL;#u3!L6EyWME(rBSIo?Hp6r
z=D%QZBtJXa<V)8?0G(RLdY2G?QzpGq@+x7zO&-D}@9pn{rJtXEb+YtPO{=@o*DMm0
z2F%BvJ~}W|LXtf{U`JE<aG{)u)0IT}8VMO+2V2xWJf@f+5#w>TUW&p_%)_CTx8hXp
z%gx`qifG>5SrUa`KG5p;|Lk@j{dHkJ2JuQsvq+y<G{M&s7(!Z&hJVp3Ur**p=n%iO
z<mscj`Or%_2W*4P;A7>!+b0SnM@0N0aN}BhKXE^)F8I;ZrK)l8vzsGV6_Xm1*17Qh
z^w9*Tth(*6t;v!C0i15Ii_)`-v~884LI{J1KBk2fBNwE>-Bt&;m-dgXcWLda>&Uk1
zHa$r!V74s0_sOUkR^IM6$J4fidtF5B4B6DKp(n`6zoPNz^(;YR&LyTv8eh1HO$%aK
zI$Ak){6#a?dB32Y{^QR^66na7wIr6j%!8)|Ip=Sm*OfP$kmi9`@bl8Gm0EjVEK@ff
zk+6;S?U2|NK!iT63@%7|sLu;WY?vor$hua7`bPL2<WuXOdON-E54O6{E(mqWfdwwB
z-$UfW*){;o@WrS<#*j1JD0^9)ooOUPOGCjVF{ng%TcLjqn_mW=z33Pdbyy<C_YiQ0
z({)rBOBr!@m`UDdQ^hrmjVZZI0q+wR7!w-@_iF%dwJnBBo4#{4D`1>efTqU#qfeoN
zcQW>Cjo+!7I|YkFTH`i?>EfhsqIs?FEsvmK%lpq>(-!Y4y_)&*KGts{Xr=}sY2{5o
zLO5a$hMT1vm6ZF`tIZKS61yujR8~ni*x~=!Y0>FnfxFdsecTSE=}&G`f2uE^wW+V;
zTeAIkUhd>OU7@EUDs$A2`AY}%J-8vHG%2D<%lrP~Rd;3`U$`*#+o>8(gXeQLDHrE#
zd#?L<yT@z&*0uVR(*$(6(xW1oOiXgtKl43Td*_gIeCFbyjK>#6(r$Ng7C&F0$u0ic
zFTiEmgSOD3u65BbG#^MLc2B?96nya3{VP-X{_lKIU?X~GL)@?J#rutTog#|sCBH@%
z@4w|+{C=VH9x1Cep};EXvr_xDx;&=mZ@JR~d_oJQjw%0M%_|dox?|A{d!r-CC*Q4H
zaXd@XQc7}Lc-B<W&70DAtG9N&`zRY{d?pd7Opf*Pj~fjvPQ6`0pBAnN<hiMT$X~JA
zuH$cdjq77indnt-rkROHM@bz0))3QS@aTGxHRI2XyKb%jl?kd2U8K4Lr7TL{u&=&R
zpSUh^`RaM4rwczS3m>`q$LY9U$DY=|ldk<z*(C3m7^J^*!y?bU0ovh}iV0%Jr`6=X
zlQh-%xw-V>>nHvlaoe{rozqlGwfb=HzOnSen<ts}HKeJ>=Uso}7?Lo*hTl~BU3}f^
zEq+oB`CR(*yFUeGeoHyXwLwzmwZxi7n*@JUsMVBO86V$oVzoDF>!+H_;_QsS(r&O$
zZC&jlKCiQk@pRyx(9~CfH#Zuo-j4hF_}VeeH(?=rtGxD4nmPZ<X7OkKbFIGqJv2Y;
z$7R7kXXeN?Jo1!1yngcRD~Wm8Q|mt3i(mS`ts}wtMf?kg|7W&buophMy@~C*`qx#h
zXJ_9~e(t>5cljn}zb8)51^&vqovE^%Qz^JatSs5~il)WEwEY)(p5(S_TvT1=zUWII
z(^J*uAtz?XC4R4}xU4VmM&skrT{F@iZBCjs@xIgpbv^;nLm%4Gmz3|{`^i_y>^uuw
ztl9Y^MY}6E{k@m`S*K;st6ggLIj3LQzc`RTwc>LNhud+D#(hQY)t_(AWx9F8?4e@e
zN*?9Y|9z)Cyq9=jmrC`v)jRgx3RUX;Fzep$cQel~Ec8&Z>|K8&abCvs6aSkI?}-fO
zGye6v_2r6h_r%ZjK0ByU&sqKH`)bciixb;o&xp5twK*a2gENogR_Si}MLwt3C!M@~
zaNn)pvm0(~xmCgW{lQP8S>KuOc<Lt!%jBFnw7y?ILA@nf;_Sm|g$4Xf<&SR!{l0Ye
zey@_m)IbfXo?~K3$8v7G?t0AzT(DbgSa&7*qqBXAsMDSmkw;HRZF+v!=zl<V|JKI&
z?1|@>*c9Hm&412K!aU5=xFb&Zy4>yKwpsBXVhVH)oO_wvo7?17v#bB3h<*II*H_a-
z4-}vIUI<Lj4&U}fF8$A#xO?Gb;Aqu$@QEho;4v)li6+34KEPvG^f=K3XcKVMjO3F;
efZDpi&gMdLHt4t);LsagoDFoeTY0(Z6b1l}Ex`-`

literal 0
HcmV?d00001

diff --git a/TRD/AliTRDseedV1.cxx b/TRD/AliTRDseedV1.cxx
index 3845362f913..4fc2f7cfba0 100644
--- a/TRD/AliTRDseedV1.cxx
+++ b/TRD/AliTRDseedV1.cxx
@@ -1018,6 +1018,8 @@ Bool_t  AliTRDseedV1::AttachClusters(AliTRDtrackingChamber *const chamber, Bool_
 // Parameters
 //  - chamber : pointer to tracking chamber container used to search the tracklet
 //  - tilt    : switch for tilt correction during road building [default true]
+//  - chgPos  : mark same[kFALSE] and opposite[kTRUE] sign tracks with respect to Bz field sign [default true]
+//  - ev      : event number for debug purposes [default = -1]
 // Output
 //  - true    : if tracklet found successfully. Failure can happend because of the following:
 //      -
@@ -1075,13 +1077,13 @@ Bool_t  AliTRDseedV1::AttachClusters(AliTRDtrackingChamber *const chamber, Bool_
   if(fkReconstructor->IsHLT()) cp.SetRPhiMethod(AliTRDcluster::kCOG);
   if(!IsCalibrated()) Calibrate();
 
-  Int_t kroadyShift(0);
+/*  Int_t kroadyShift(0);
   Float_t bz(AliTrackerBase::GetBz());
   if(TMath::Abs(bz)>2.){
     if(bz<0.) kroadyShift = chgPos ? +1 : -1;
     else kroadyShift = chgPos ? -1 : +1;
-  }
-  AliDebug(4, Form("\n       syTrk[cm]=%4.2f dydxTrk[deg]=%+6.2f rs[%d] Chg[%c] rY[cm]=%4.2f rZ[cm]=%5.2f TC[%c]", syRef, TMath::ATan(fYref[1])*TMath::RadToDeg(), kroadyShift, chgPos?'+':'-', kroady, kroadz, tilt?'y':'n'));
+  }*/
+  AliDebug(4, Form("\n       syTrk[cm]=%4.2f dydxTrk[deg]=%+6.2f Chg[%c] rY[cm]=%4.2f rZ[cm]=%5.2f TC[%c]", syRef, TMath::ATan(fYref[1])*TMath::RadToDeg(), chgPos?'+':'-', kroady, kroadz, tilt?'y':'n'));
   Double_t phiTrk(TMath::ATan(fYref[1])),
            thtTrk(TMath::ATan(fZref[1]));
 
@@ -2033,6 +2035,9 @@ Bool_t AliTRDseedV1::FitRobust(Bool_t chg)
     } else { 
       kScalePulls = attach->GetScaleCov();//*lyScaler;
     }
+    // Retrieve chamber status
+    SetChmbGood(calibration->IsChamberGood(fDet));
+    if(!IsChmbGood()) kScalePulls*=10.;
   }  
   Double_t xc[kNclusters], yc[kNclusters], sy[kNclusters];
   Int_t n(0),           // clusters used in fit 
diff --git a/TRD/AliTRDseedV1.h b/TRD/AliTRDseedV1.h
index 67350b7ea0d..d83427f1c05 100644
--- a/TRD/AliTRDseedV1.h
+++ b/TRD/AliTRDseedV1.h
@@ -66,6 +66,7 @@ public:
    ,kKink       = BIT(18) // kink prolongation tracklet
    ,kStandAlone = BIT(19) // tracklet build during stand alone track finding
    ,kPrimary    = BIT(20) // tracklet from a primary track candidate
+   ,kChmbGood   = BIT(21) // status of the detector from calibration view point
   };
 
   enum ETRDtrackletError { // up to 8 bits
@@ -96,6 +97,7 @@ public:
   void      Init(const AliRieman *fit);
   Bool_t    IsEqual(const TObject *inTracklet) const;
   Bool_t    IsCalibrated() const     { return TestBit(kCalib);}
+  Bool_t    IsChmbGood() const       { return TestBit(kChmbGood);}
   Bool_t    IsOwner() const          { return TestBit(kOwner);}
   Bool_t    IsKink() const           { return TestBit(kKink);}
   Bool_t    IsPrimary() const        { return TestBit(kPrimary);}
@@ -178,6 +180,7 @@ public:
   void      Reset(Option_t *opt="");
 
   void      SetC(Float_t c, Int_t typ=0) { fC[typ] = c;}
+  void      SetChmbGood(Bool_t k = kTRUE){ SetBit(kChmbGood, k);}
   void      SetChi2(Float_t chi2)    { fChi2 = chi2;}
   inline void SetCovRef(const Double_t *cov);
   void      SetErrorMsg(ETRDtrackletError err)  { SETBIT(fErrorMsg, err);}
diff --git a/TRD/AliTRDtrackV1.cxx b/TRD/AliTRDtrackV1.cxx
index 38148d07442..28f01efd64f 100644
--- a/TRD/AliTRDtrackV1.cxx
+++ b/TRD/AliTRDtrackV1.cxx
@@ -882,8 +882,11 @@ void AliTRDtrackV1::UpdateESDtrack(AliESDtrack *track)
     // first implementation store no. of time bins filled in tracklet (5bits  see "y" bits) and
     // no. of double clusters in case of pad row cross (4bits see "x" bits)
     // bit map for tracklet quality xxxxyyyyy
+    // 27.10.11 A.Bercuci
+    // add chamber status bit "z" bit 
+    // bit map for tracklet quality zxxxxyyyyy
     Int_t nCross(fTracklet[ip]->IsRowCross()?fTracklet[ip]->GetTBcross():0);
-    spd = Double_t(fTracklet[ip]->GetTBoccupancy() | (nCross<<5));
+    spd = Double_t(fTracklet[ip]->GetTBoccupancy() | (nCross<<5) | (fTracklet[ip]->IsChmbGood()<<9));
     track->SetTRDmomentum(p, ip, &spd);
     track->SetTRDslice(fTracklet[ip]->GetdQdl(), ip, 0); // Set Summed dEdx into the first slice
   }
diff --git a/TRD/AliTRDtrackerV1.cxx b/TRD/AliTRDtrackerV1.cxx
index e49f3787d2c..45b812697aa 100644
--- a/TRD/AliTRDtrackerV1.cxx
+++ b/TRD/AliTRDtrackerV1.cxx
@@ -875,7 +875,11 @@ Int_t AliTRDtrackerV1::FollowBackProlongation(AliTRDtrackV1 &t)
         AliDebug(4, "Failed Tracklet Init");
         break;
       }
-      if(!ptrTracklet->AttachClusters(chamber, kTRUE, t.Charge()>0?kTRUE:kFALSE, fEventInFile)){
+      // Select attachment base on track to B field sign not only track charge which is buggy
+      // mark kFALSE same sign tracks and kTRUE opposite sign tracks
+      // A.Bercuci 3.11.2011
+      Float_t prod(t.GetBz()*t.Charge());
+      if(!ptrTracklet->AttachClusters(chamber, kTRUE, prod<0.?kTRUE:kFALSE, fEventInFile)){
         t.SetStatus(AliTRDtrackV1::kNoAttach, ily);
         if(debugLevel>3){
           AliTRDseedV1 trackletCp(*ptrTracklet);
@@ -984,7 +988,10 @@ Int_t AliTRDtrackerV1::FollowBackProlongation(AliTRDtrackV1 &t)
   
 
     // register tracklet with the tracker and track
-    ptrTracklet->Update(&t);
+    // Save inside the tracklet the track parameters BEFORE track update.
+    // Commented out their overwriting AFTER track update
+    // A.Bercuci 3.11.2011
+    //ptrTracklet->Update(&t); 
     ptrTracklet = SetTracklet(ptrTracklet);
     Int_t index(fTracklets->GetEntriesFast()-1);
     t.SetTracklet(ptrTracklet, index);
diff --git a/TRD/AliTRDtransform.cxx b/TRD/AliTRDtransform.cxx
index 6be6039c34b..b17d918b98a 100644
--- a/TRD/AliTRDtransform.cxx
+++ b/TRD/AliTRDtransform.cxx
@@ -48,8 +48,10 @@ AliTRDtransform::AliTRDtransform()
   ,fCalPRFROC(0x0)
   ,fkCalVdriftDet(0x0)
   ,fkCalT0Det(0x0)
+  ,fkCalExBDet(0x0)
   ,fCalVdriftDetValue(0)
   ,fCalT0DetValue(0)
+  ,fCalExBDetValue(0)
   ,fSamplingFrequency(0)
   ,fPadPlane(0x0)
   ,fZShiftIdeal(0)
@@ -73,6 +75,7 @@ AliTRDtransform::AliTRDtransform()
   // Get the calibration objects for the global calibration
   fkCalVdriftDet     = fCalibration->GetVdriftDet();
   fkCalT0Det         = fCalibration->GetT0Det();
+  fkCalExBDet        = fCalibration->GetExBDet();
 
 }
 
@@ -87,8 +90,10 @@ AliTRDtransform::AliTRDtransform(Int_t det)
   ,fCalPRFROC(0x0)
   ,fkCalVdriftDet(0x0)
   ,fkCalT0Det(0x0)
+  ,fkCalExBDet(0x0)
   ,fCalVdriftDetValue(0)
   ,fCalT0DetValue(0)
+  ,fCalExBDetValue(0)
   ,fSamplingFrequency(0)
   ,fPadPlane(0x0)
   ,fZShiftIdeal(0)
@@ -112,6 +117,7 @@ AliTRDtransform::AliTRDtransform(Int_t det)
   // Get the calibration objects for the global calibration
   fkCalVdriftDet     = fCalibration->GetVdriftDet();
   fkCalT0Det         = fCalibration->GetT0Det();
+  fkCalExBDet        = fCalibration->GetExBDet();
 
   SetDetector(det);
 
@@ -128,8 +134,10 @@ AliTRDtransform::AliTRDtransform(const AliTRDtransform &t)
   ,fCalPRFROC(0x0)
   ,fkCalVdriftDet(0x0)
   ,fkCalT0Det(0x0)
+  ,fkCalExBDet(0x0)
   ,fCalVdriftDetValue(0)
   ,fCalT0DetValue(0)
+  ,fCalExBDetValue(0)
   ,fSamplingFrequency(0)
   ,fPadPlane(0x0)
   ,fZShiftIdeal(0)
@@ -151,7 +159,7 @@ AliTRDtransform::AliTRDtransform(const AliTRDtransform &t)
   }
   fkCalVdriftDet     = fCalibration->GetVdriftDet();
   fkCalT0Det         = fCalibration->GetT0Det();
-
+  fkCalExBDet        = fCalibration->GetExBDet();
 }
 
 //_____________________________________________________________________________
@@ -191,16 +199,19 @@ void AliTRDtransform::Copy(TObject &t) const
   if (fCalibration) {
     ((AliTRDtransform &) t).fkCalVdriftDet   = fCalibration->GetVdriftDet();
     ((AliTRDtransform &) t).fkCalT0Det       = fCalibration->GetT0Det();
+    ((AliTRDtransform &) t).fkCalExBDet      = fCalibration->GetExBDet();
   }
   else {
     ((AliTRDtransform &) t).fkCalVdriftDet   = 0;
     ((AliTRDtransform &) t).fkCalT0Det       = 0;
+    ((AliTRDtransform &) t).fkCalExBDet      = 0;
   }
   ((AliTRDtransform &) t).fCalVdriftROC      = 0x0;
   ((AliTRDtransform &) t).fCalT0ROC          = 0x0;
   ((AliTRDtransform &) t).fCalPRFROC         = 0x0;
   ((AliTRDtransform &) t).fCalVdriftDetValue = 0; 
   ((AliTRDtransform &) t).fCalT0DetValue     = 0;
+  ((AliTRDtransform &) t).fCalExBDetValue     = 0;
   ((AliTRDtransform &) t).fSamplingFrequency = 0;
   ((AliTRDtransform &) t).fPadPlane          = 0x0;
   ((AliTRDtransform &) t).fZShiftIdeal       = 0;
@@ -236,6 +247,7 @@ void AliTRDtransform::SetDetector(Int_t det)
   // Get the detector wise defined calibration values
   fCalVdriftDetValue = fkCalVdriftDet->GetValue(det);
   fCalT0DetValue     = fkCalT0Det->GetValue(det);
+  fCalExBDetValue    = fkCalExBDet->GetValue(det);
 
   // Shift needed to define Z-position relative to middle of chamber
   Int_t layer        = Geometry().GetLayer(det);
@@ -275,7 +287,7 @@ Bool_t AliTRDtransform::Transform(AliTRDcluster *c)
   Double_t t0  = fCalT0DetValue     + fCalT0ROC->GetValue(col,row);
   t0 /= fSamplingFrequency;
   // ExB correction
-  Double_t exb = AliTRDCommonParam::Instance()->GetOmegaTau(vd);
+  Double_t exb = fCalExBDetValue;//AliTRDCommonParam::Instance()->GetOmegaTau(vd);
 
   Float_t x = c->GetXloc(t0, vd);
 
diff --git a/TRD/AliTRDtransform.h b/TRD/AliTRDtransform.h
index b46e704ffe6..18f3b625b8a 100644
--- a/TRD/AliTRDtransform.h
+++ b/TRD/AliTRDtransform.h
@@ -53,8 +53,10 @@ class AliTRDtransform : public TObject {
   AliTRDCalROC       *fCalPRFROC;           //  Pad wise PRF calibration object
   const AliTRDCalDet *fkCalVdriftDet;       //  ROC wise Vdrift calibration object
   const AliTRDCalDet *fkCalT0Det;           //  ROC wise T0 calibration object
+  const AliTRDCalDet *fkCalExBDet;          //  ROC wise ExB calibration object
   Double_t            fCalVdriftDetValue;   //  ROC wise Vdrift calibration value
   Double_t            fCalT0DetValue;       //  ROC wise T0 calibration value
+  Double_t            fCalExBDetValue;      //  Det wise ExB calibration value
 
   Double_t            fSamplingFrequency;   //  ADC sampling frequency
 
diff --git a/TRD/Cal/AliTRDCalTrkAttach.cxx b/TRD/Cal/AliTRDCalTrkAttach.cxx
index 88dd3610abd..2a79169a667 100644
--- a/TRD/Cal/AliTRDCalTrkAttach.cxx
+++ b/TRD/Cal/AliTRDCalTrkAttach.cxx
@@ -34,8 +34,9 @@
 
 #include "AliLog.h"
 
-#include "AliTRDCalTrkAttach.h"
 #include "AliTRDcalibDB.h"
+#include "AliTRDgeometry.h"
+#include "AliTRDCalTrkAttach.h"
 
 ClassImp(AliTRDCalTrkAttach)
 
@@ -71,13 +72,13 @@ Double_t AliTRDCalTrkAttach::CookLikelihood(Bool_t chg, Int_t ly, Float_t pt, Fl
     return 0.;
   }
   //check likelihood array  
-  if (!fLike || fLike->GetEntries() != 6*2*kNcalib) {
+  if (!fLike || fLike->GetEntries() != AliTRDgeometry::kNlayer*kNcharge*kNcalib) {
     AliError("No usable AttachClusters calib object.");
     return 0.;
   }
-  
+  Int_t offset(ly*kNcalib*kNcharge); // offset for layer
   TGraphErrors *g(NULL);
-  if(!(g = (TGraphErrors*)fLike->At(ly*8+2*Int_t(kResPos)+Int_t(chg)))){
+  if(!(g = (TGraphErrors*)fLike->At(offset+Int_t(kResPos)*kNcharge+Int_t(chg)))){
     AliError("Failed retrieving AttachClusters graph.");
     return 0.;
   }
@@ -109,7 +110,7 @@ Double_t AliTRDCalTrkAttach::CookLikelihood(Bool_t chg, Int_t ly, Float_t pt, Fl
   dd      = (dyr - dym)/sym; dd*=dd;
   f[0] = TMath::Exp(-0.5*dd);
   // +++ process dphi
-  if(!(g = (TGraphErrors*)fLike->At(ly*8+2*Int_t(kResAng)+Int_t(chg)))){
+  if(!(g = (TGraphErrors*)fLike->At(offset+Int_t(kResAng)*kNcharge+Int_t(chg)))){
     AliError("Failed retrieving AttachClusters graph.");
     return 0.;
   }
@@ -119,7 +120,7 @@ Double_t AliTRDCalTrkAttach::CookLikelihood(Bool_t chg, Int_t ly, Float_t pt, Fl
   dd      = (dphi - dpm)/spm; dd*=dd;
   f[1] = TMath::Exp(-0.5*dd);
   // +++ process no of clusters
-  if(!(g = (TGraphErrors*)fLike->At(ly*8+2*Int_t(kNclMean)+Int_t(chg)))){
+  if(!(g = (TGraphErrors*)fLike->At(offset+Int_t(kNclMean)*kNcharge+Int_t(chg)))){
     AliError("Failed retrieving AttachClusters graph.");
     return 0.;
   }
@@ -129,7 +130,7 @@ Double_t AliTRDCalTrkAttach::CookLikelihood(Bool_t chg, Int_t ly, Float_t pt, Fl
  
   // process phi dependences
   // +++ process <s>/s
-  if(!(g = (TGraphErrors*)fLike->At(ly*8+2*Int_t(kSigma)+Int_t(chg)))){
+  if(!(g = (TGraphErrors*)fLike->At(offset+Int_t(kSigma)*kNcharge+Int_t(chg)))){
     AliError("Failed retrieving AttachClusters graph.");
     return 0.;
   }
@@ -156,7 +157,7 @@ Double_t AliTRDCalTrkAttach::CookLikelihood(Bool_t chg, Int_t ly, Float_t pt, Fl
   Double_t sm = boundary?y1:((phiTrk*(y0-y1) + (x0*y1-x1*y0))*invdx),
            ssm = 0.5*(g->GetErrorY(ip)+g->GetErrorY(jp));
   dd      = (sr - sm)/ssm; dd*=dd;
-  f[3] = TMath::Exp(-0.5*dd);
+  f[3] = 1.;//TMath::Exp(-0.5*dd);
 
   // Calculate likelihood
   Double_t length = f[0]*f[0]+f[1]*f[1]+f[2]*f[2]+f[3]*f[3];
@@ -188,12 +189,12 @@ Bool_t AliTRDCalTrkAttach::LoadReferences(const Char_t *file)
          cs[2] = {'n', 'p'};
 
   if(fLike) fLike->Clear();
-  else fLike = new TObjArray(6*2*kNcalib);
-  for(Int_t ily(0); ily<6; ily++){
+  else fLike = new TObjArray(AliTRDgeometry::kNlayer*kNcharge*kNcalib);
+  for(Int_t ily(0); ily<AliTRDgeometry::kNlayer; ily++){
     for(Int_t icalib(0); icalib<kNcalib; icalib++){
-      for(Int_t isgn(0); isgn<2; isgn++){
+      for(Int_t isgn(0); isgn<kNcharge; isgn++){
         if(!(g = (TGraphErrors*)gFile->Get(Form("%c%c%d", co[icalib], cs[isgn], ily)))) return kFALSE;
-        fLike->AddAt(g, ily*8+2*icalib+isgn);
+        fLike->AddAt(g, kNcharge*(ily*kNcalib+icalib)+isgn);
       }
     }
   }
diff --git a/TRD/Cal/AliTRDCalTrkAttach.h b/TRD/Cal/AliTRDCalTrkAttach.h
index 24a824c3367..1d786465077 100644
--- a/TRD/Cal/AliTRDCalTrkAttach.h
+++ b/TRD/Cal/AliTRDCalTrkAttach.h
@@ -19,13 +19,16 @@ class TObjArray;
 class AliTRDCalTrkAttach : public TNamed
 {
 public:
-  enum ETRDCalTrkAttachCalib {    
+  enum ETRDCalTrkAttachCalib {
     kResPos = 0    // relative position residual location
    ,kResAng        // angular residual location
    ,kSigma         // relative error location
    ,kNclMean       // mean no. of clusters/tracklet location
    ,kNcalib        // no. of calib objects
   };
+  enum ETRDCalTrkAttachParam {
+    kNcharge = 2
+  };
 
   AliTRDCalTrkAttach();
   virtual ~AliTRDCalTrkAttach();
diff --git a/TRD/Cal/AliTRDmakeTrkDB.C b/TRD/Cal/AliTRDmakeTrkDB.C
index 665222d1e58..50ba824dd66 100644
--- a/TRD/Cal/AliTRDmakeTrkDB.C
+++ b/TRD/Cal/AliTRDmakeTrkDB.C
@@ -10,11 +10,11 @@ void AliTRDmakeTrkDB(const Char_t *file)
   metaData->SetObjectClassName("TObjArray");
   metaData->SetResponsible("Alexandru Bercuci");
   metaData->SetBeamPeriod(1);
-  metaData->SetAliRootVersion("05-27-06b"); //root version
+  metaData->SetAliRootVersion("05-30-02"); //root version
   metaData->SetComment(
     "Likelihoods for Attach Cluster.\n"
-    " Tunned on Pb-Pb run 137161.");
-  AliCDBId id("TRD/Calib/TrkAttach", 95352, AliCDBRunRange::Infinity()); 
+    " Tunned on p-p run 159580.");
+  AliCDBId id("TRD/Calib/TrkAttach", 151536, AliCDBRunRange::Infinity());
   AliTRDCalTrkAttach attach;
   if(!attach.LoadReferences(file)) return;
 //   attach.SetNsgmDy(Int_t ns0, Int_t ns1);
-- 
2.43.0