From f790fd5a15e0cb6a63ffe801966aa12c028b8372 Mon Sep 17 00:00:00 2001
From: ChunzhengLab <chunzheng.wang@icloud.com>
Date: Tue, 17 Jun 2025 17:34:15 +0200
Subject: [PATCH 1/5] move the energy deposition wrt. centre of response

---
 Detectors/Upgrades/ITS3/simulation/src/Digitizer.cxx | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/Detectors/Upgrades/ITS3/simulation/src/Digitizer.cxx b/Detectors/Upgrades/ITS3/simulation/src/Digitizer.cxx
index b91e17890a6d8..70425a1b4139f 100644
--- a/Detectors/Upgrades/ITS3/simulation/src/Digitizer.cxx
+++ b/Detectors/Upgrades/ITS3/simulation/src/Digitizer.cxx
@@ -90,7 +90,8 @@ void Digitizer::init()
     } else {
       LOGP(fatal, "ResponseFunction '{}' not implemented!", func);
     }
-    mSimRespIBShift = mSimRespIB->getDepthMax() - constants::silicon::thickness / 2.f;
+    mSimRespIB->computeCentreFromData();
+    mSimRespIBShift = mSimRespIB->getDepthMax() - constants::silicon::thickness / 2.f - mSimRespIB->getRespCentreDep();
     mSimRespOBShift = mSimRespOB->getDepthMax() - SegmentationOB::SensorLayerThickness / 2.f;
   }
 

From 6c8c9047a49f5f834331f5ae9ddfd76d4e74de22 Mon Sep 17 00:00:00 2001
From: ChunzhengLab <chunzheng.wang@icloud.com>
Date: Tue, 17 Jun 2025 19:22:26 +0200
Subject: [PATCH 2/5] fix the order of IB/OB for noise seeding

---
 .../Upgrades/ITS3/simulation/src/ChipDigitsContainer.cxx  | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/Detectors/Upgrades/ITS3/simulation/src/ChipDigitsContainer.cxx b/Detectors/Upgrades/ITS3/simulation/src/ChipDigitsContainer.cxx
index 102b15863683e..efe878536687d 100644
--- a/Detectors/Upgrades/ITS3/simulation/src/ChipDigitsContainer.cxx
+++ b/Detectors/Upgrades/ITS3/simulation/src/ChipDigitsContainer.cxx
@@ -31,12 +31,12 @@ void ChipDigitsContainer::addNoise(UInt_t rofMin, UInt_t rofMax, const o2::its3:
   int nel = 0;
 
   if (isIB()) {
-    // Inner barrel: use ITS3-specific noise interface with OB segmentation.
-    mean = params->getIBNoisePerPixel() * SegmentationOB::NPixels;
+    // Inner barrel: use ITS3-specific noise interface with IB segmentation.
+    mean = params->getIBNoisePerPixel() * SegmentationIB::NPixels;
     nel = static_cast<int>(params->getIBChargeThreshold() * 1.1);
   } else {
-    // Outer barrel: use base class noise interface with IB segmentation.
-    mean = params->getNoisePerPixel() * SegmentationIB::NPixels;
+    // Outer barrel: use base class noise interface with OB segmentation.
+    mean = params->getNoisePerPixel() * SegmentationOB::NPixels;
     nel = static_cast<int>(params->getChargeThreshold() * 1.1);
   }
 

From ccd75ae5c05c8e86e6f11f27c18ec77d11e47202 Mon Sep 17 00:00:00 2001
From: ChunzhengLab <chunzheng.wang@icloud.com>
Date: Tue, 17 Jun 2025 20:57:16 +0200
Subject: [PATCH 3/5] print out the resp centre in the CheckChipResponseFile

---
 .../ITS3/macros/test/CheckChipResponseFile.C  | 51 +++++++++++--------
 1 file changed, 31 insertions(+), 20 deletions(-)

diff --git a/Detectors/Upgrades/ITS3/macros/test/CheckChipResponseFile.C b/Detectors/Upgrades/ITS3/macros/test/CheckChipResponseFile.C
index 996a99d87ecbc..5cacf2eab9900 100644
--- a/Detectors/Upgrades/ITS3/macros/test/CheckChipResponseFile.C
+++ b/Detectors/Upgrades/ITS3/macros/test/CheckChipResponseFile.C
@@ -24,6 +24,7 @@
 
 #define ENABLE_UPGRADES
 #include "ITSMFTSimulation/AlpideSimResponse.h"
+#include "ITS3Simulation/ChipSimResponse.h"
 
 #include "ITS3Base/SegmentationMosaix.h"
 #include "fairlogger/Logger.h"
@@ -34,21 +35,21 @@ using SegmentationMosaix = o2::its3::SegmentationMosaix;
 double um2cm(double um) { return um * 1e-4; }
 double cm2um(double cm) { return cm * 1e+4; }
 
-o2::itsmft::AlpideSimResponse *mAlpSimResp0 = nullptr,
-                              *mAlpSimResp1 = nullptr,
-                              *mAptSimResp1 = nullptr;
+std::unique_ptr<o2::its3::ChipSimResponse> mAlpSimResp0, mAlpSimResp1, mAptSimResp1;
 
-o2::itsmft::AlpideSimResponse* loadResponse(const std::string& fileName, const std::string& respName)
+std::unique_ptr<o2::its3::ChipSimResponse> loadResponse(const std::string& fileName, const std::string& respName)
 {
   TFile* f = TFile::Open(fileName.data());
   if (!f) {
     std::cerr << fileName << " not found" << std::endl;
     return nullptr;
   }
-  auto resp = (o2::itsmft::AlpideSimResponse*)f->Get(respName.data());
-  if (!resp)
+  auto base = f->Get<o2::itsmft::AlpideSimResponse>(respName.c_str());
+  if (!base) {
     std::cerr << respName << " not found in " << fileName << std::endl;
-  return resp;
+    return nullptr;
+  }
+  return std::make_unique<o2::its3::ChipSimResponse>(base);
 }
 
 void LoadRespFunc()
@@ -56,18 +57,28 @@ void LoadRespFunc()
   std::string AptsFile = "$(O2_ROOT)/share/Detectors/Upgrades/ITS3/data/ITS3ChipResponseData/APTSResponseData.root";
   std::string AlpideFile = "$(O2_ROOT)/share/Detectors/ITSMFT/data/AlpideResponseData/AlpideResponseData.root";
 
+  std::cout<<"=====================\n";
+  LOGP(info, "ALPIDE Vbb=0V response");
   mAlpSimResp0 = loadResponse(AlpideFile, "response0"); // Vbb=0V
-  LOG(info) << "ALPIDE Vbb=0V response" << std::endl;
+  mAlpSimResp0->computeCentreFromData();
   mAlpSimResp0->print();
+  LOGP(info, "Response Centre {}", mAlpSimResp0->getRespCentreDep());
+  std::cout<<"=====================\n";
+  LOGP(info, "ALPIDE Vbb=-3V response");
   mAlpSimResp1 = loadResponse(AlpideFile, "response1"); // Vbb=-3V
-  LOG(info) << "ALPIDE Vbb=-3V response" << std::endl;
+  mAlpSimResp1->computeCentreFromData();
   mAlpSimResp1->print();
+  LOGP(info, "Response Centre {}", mAlpSimResp1->getRespCentreDep());
+  std::cout<<"=====================\n";
+  LOGP(info, "APTS response");
   mAptSimResp1 = loadResponse(AptsFile, "response1"); // APTS
-  LOG(info) << "APTS response" << std::endl;
+  mAptSimResp1->computeCentreFromData();
   mAptSimResp1->print();
+  LOGP(info, "Response Centre {}", mAptSimResp1->getRespCentreDep());
+  std::cout<<"=====================\n";
 }
 
-std::vector<float> getCollectionSeediciencies(o2::itsmft::AlpideSimResponse* resp,
+std::vector<float> getCollectionSeediciencies(o2::its3::ChipSimResponse* resp,
                                               const std::vector<float>& depths)
 {
   std::vector<float> seed;
@@ -81,7 +92,7 @@ std::vector<float> getCollectionSeediciencies(o2::itsmft::AlpideSimResponse* res
   return seed;
 }
 
-std::vector<float> getShareValues(o2::itsmft::AlpideSimResponse* resp,
+std::vector<float> getShareValues(o2::its3::ChipSimResponse* resp,
                                   const std::vector<float>& depths)
 {
   std::vector<float> share;
@@ -103,7 +114,7 @@ std::vector<float> getShareValues(o2::itsmft::AlpideSimResponse* resp,
   return share;
 }
 
-std::vector<float> getEffValues(o2::itsmft::AlpideSimResponse* resp,
+std::vector<float> getEffValues(o2::its3::ChipSimResponse* resp,
                                 const std::vector<float>& depths)
 {
   std::vector<float> all;
@@ -135,7 +146,7 @@ void CheckChipResponseFile()
 
   int colors[] = {kOrange + 7, kRed + 1, kAzure + 4};
   struct RespInfo {
-    o2::itsmft::AlpideSimResponse* resp;
+    std::unique_ptr<o2::its3::ChipSimResponse>& resp;
     std::string title;
     int color;
   };
@@ -154,11 +165,11 @@ void CheckChipResponseFile()
   for (auto& r : responses) {
     if (!r.resp)
       continue;
-    auto seed = getCollectionSeediciencies(r.resp, vecDepth);
-    auto shr = getShareValues(r.resp, vecDepth);
-    auto all = getEffValues(r.resp, vecDepth);
+    auto seed = getCollectionSeediciencies(r.resp.get(), vecDepth);
+    auto shr = getShareValues(r.resp.get(), vecDepth);
+    auto all = getEffValues(r.resp.get(), vecDepth);
 
-    TGraph* grSeed = new TGraph(vecDepth.size(), vecDepth.data(), seed.data());
+    auto grSeed = new TGraph(vecDepth.size(), vecDepth.data(), seed.data());
     grSeed->SetTitle(Form("%s seed", r.title.c_str()));
     grSeed->SetLineColor(r.color);
     grSeed->SetLineWidth(2);
@@ -168,7 +179,7 @@ void CheckChipResponseFile()
     grSeed->Draw("SAME LP");
     leg->AddEntry(grSeed, Form("%s seed", r.title.c_str()), "lp");
 
-    TGraph* grShare = new TGraph(vecDepth.size(), vecDepth.data(), shr.data());
+    auto grShare = new TGraph(vecDepth.size(), vecDepth.data(), shr.data());
     grShare->SetLineColor(r.color);
     grShare->SetLineWidth(2);
     grShare->SetMarkerColor(r.color);
@@ -177,7 +188,7 @@ void CheckChipResponseFile()
     grShare->Draw("SAME LP");
     leg->AddEntry(grShare, Form("%s share", r.title.c_str()), "p");
 
-    TGraph* grEff = new TGraph(vecDepth.size(), vecDepth.data(), all.data());
+    auto grEff = new TGraph(vecDepth.size(), vecDepth.data(), all.data());
     grEff->SetLineColor(r.color);
     grEff->SetLineWidth(2);
     grEff->SetMarkerColor(r.color);

From fb4c608e42380577d3fdde6dfbaa32ef1e80c536 Mon Sep 17 00:00:00 2001
From: ChunzhengLab <chunzheng.wang@icloud.com>
Date: Thu, 19 Jun 2025 20:48:27 +0200
Subject: [PATCH 4/5] find the real centre

---
 ChipResponse.pdf                              | Bin 0 -> 23047 bytes
 .../ITS3/simulation/src/ChipSimResponse.cxx   |  31 +++++++++++++++---
 .../ITS3/simulation/src/Digitizer.cxx         |   3 +-
 core_dump_74974                               |   0
 4 files changed, 28 insertions(+), 6 deletions(-)
 create mode 100644 ChipResponse.pdf
 create mode 100644 core_dump_74974

diff --git a/ChipResponse.pdf b/ChipResponse.pdf
new file mode 100644
index 0000000000000000000000000000000000000000..0351c0af364bc2e7f77ff71398f4fe0524877b0f
GIT binary patch
literal 23047
zcmd741zeTM_dgCuH_{<oLTS18a)FDqq;!LH!$mrz8zdw|x=TV)S^=e7L{UKjl?J6t
zy8JyCWp`cuth=w*e*gQ|*GuQlGc#xAJ!ekLbDncvRZ4~%auJ4mermdZ02cxRgB;Cl
zafO9(c{IJ8EI~Yyrf#Noj@Gz5s;1VKt{@oDKpmGy#of&g;Q;*R{rVgEP|48(Xs_jJ
zdHN7nR20|J!2&rR^bg}ox>%aJIl6$@)Kye8LE4rst_Vj55dTG9=tVF;J1)=3bD&Df
z)XfsaCUpr4hVp^=A%ajaL;wQj1VdTDU{;_VK<`&gG!brgz!;J?2q$$*S0_gYSIdh|
z7FIxK1xqgvM;8m%ud`tP!&LeHFjZ9-M{^BJH;_IL027Eu)6&z;5SK^B(ZLOg5Qs+x
z0un$10F<C0__q=-Nbp+;27-XU75P9AB(}(w{2&Mt4rEaP1UUf^D8fOI6KsK^AP9n-
z1hN6pAHsiHhJYZa)PNBnr}Tg!Ag2_8F`%a;fkB`rG?9zrk@j?x)o?>X3jI1b0JNi>
zql<=<sktRcSQyx8S2zgzb!Ny8s-|vkmM#t;=+`a(eJ_!7_?|~j0O8SaH*-5-86@LC
zaCwwX?E!v}0ZNxmT`iGq00a^Cmag1a0G2|S19NgPceFq_Sc7<U5DpR!u87l$KWqnb
zmH&Vm^hau!9qlas9xWt_|9}?n&uRUW)t&IiKOhPFBS|?+yX%&22y@fFx4IJs{0Fr7
zeoBj5#mo-j>~8rtY5wuxe-$Wy$t?Uor}_7{@Z>!GgGCDbZ-W19p~zG659q@GNLSL)
z-34Lk@;4#<qk#GwkP7~U(m&ipqy+j0YlDD)M)Du*;Xjoy--(!SU~(hr{j)9nhoT1Z
zovhLDvNv<I`-9**k;(rAwV#L+bw_(shrh=bCu0Afko!)Q=$SfMNdc1B%+$^G4{JEl
zAOB>~?_|agkep~2<ZysK{ia=bq%5x^%q`VrFB|@5D<D<kSE4-M6!3}uHT+!#o#>=*
z^8)dxf*~O2H!Z0OhJv8qw4^GS7X<yLB~`(IX8filRl$59-fvn`70eIf{ic~!!2%%O
zZ>m`p3<vRk)6J@2K@jgZ<*W)wY~F9$Srw3|yr*;ljm&#W7Z5_cr*r|0%zH`~(8#b;
zx`0}Tozex=I_#7#pw?lhbO9EEozeyN7j{Y)*cKQ-7Z+(Fr~;tEP6<PKL9kQ807kG=
z!hp)>J0bkXz5V@EpO|{5Tf(Dc>gIy*)JIMZus)FX1Qc#~`t(HLKzIOa3XoamfB^`Q
z1WXm+F~Y&^<RKJL@_=n<{jKHY%Z{G<$dLqrg&FebYPeZC0cD=6Cj}UfqNRhin+-@1
z!XtyQ0}Kux89Ts8k+K9V1504TUEKg{(VhTTE|Oygs@bDC0aXP9dl@*9uUY04CY*bV
zQWdPws+vc~c7LLj*uI4LTJEhj%R46dx1L!daPGJSK8}%~!viIqNrVLVPz71yGeU2O
zD}!}Fl;VM?D8A-ot4Ox{O64~O{9XS16dCYE0pv~qf02cM4Ih9Ze#3|V2YjHv!iOaa
zf*S}XlibN<lHF}4f5nMSI5naA!A<pE?#&^AcZ{J8DR?hpz%--sE99gqdCE@3pGmCx
zXQDG>@X$EvZD_UL^7t2MTc)Wh?#56~(=cN;py8wBn~-1n2|&L~%Krh-DL(%yJP`gL
z@PY9E#z4IP2%t)~ZQS4xI6NYQP@ddvegF6bAlfNqs+g+V*6kh7)IP1{?q0I$zlv`V
z&B?*BCS~W!8qn21Q$AC6@!oMvRD2xC8BtVBl>P$p`#(YGcT@6zAOy&ruPg)<{w)@Q
zfPn)G2A`q?6u+?%vh?ks@&6L1zdC4A-wxV2;FzIlr#-{a$Y>|;Bk#c=HbI%GBWU;$
zi4#HlO`0P`!l@uR6*FMK>HCk6;QtlJ{9Pp4Fw_#BVK5Ppxzq*Mqt&AT2WRHxPr*P6
z)xQP<DNd1+1*yCKCCU4@*y7&<!xC3L9WP(p*7Y|qW_1ZVf5I8R`v(8(jQm$%!2d1I
zKqKSlQ|*~WV@HSn&Kt8LKZOIi*uUnD|C56P+$w+}e<?B`yn?^NLn^M4unmka<q^>N
zT#;Wb$@k1z$<7v4_iLO%G%SzT4u&5N5Y`XxWX>>o22n-dOQ+&Izo@Nu(bi%o?2Km~
zn~&ksLBj#xDNI`Vpi5nEQsu3a<trPd%9+PSZdEK@%b1Uh(-YtuCq5kPI@jKlI^`+S
zzSr38_|A1CAu(5W7;9^j(gX!p7vlF5F8kdv`(Ir4e*){TB{2m2FS8kjx{QqX1!6QV
zi`-j%mJVgsn0(k=sAj3bEV=2?;u}RKv4V1AWcmDQcozb>aNkmWliT>>^eI2V4XMKa
z19tp1Zi4@w2K&S5417nBXCMh8vQWNAso0Cs7FNm^X3HO@oA{f4ixQy5`3X))?)gic
zASV}YNSTf-{eLOcA^fmk9bAJr2=d@cd*n<pv?81mj8G#{-A%;pz=3RMxeOSDu7|$M
ziyS)`*Zy1-DiqNww>17qhP6R+SxfZRnpVu27nEHO`6xxOx)M6BpxM70Ub#T^=Dn3(
zqs(|%)<cy#8uXFqvCca56q8j%S;CZ%;miXz#{HErTdFS#qqk<hL^O1r9Uj;8iROwD
zYMp~)e!M}VA`<Hx<y*vjgieJLD<c#B6U2XMp8t1<L!rOaqr$2}4((w4Q7_zV&@Cy2
z{4&*C$ya%j@lT)L$|t;avC&^Do+sqWl;19!7~AI<cIT<g!-FTUDaTni`CA!Y9EzFx
ze_C~&dKSvJE+Lrn+R8jFFk7a&^@C6(Q8<C}(|Z~QDL%!R9^rS|A2@Ncg}q<1<hj0L
z`n-Sb@k=vbkERgEiEG&Clz1#=$yr?a_ubgYaNRq#?Vphx4LvSTv2wRcwsl{*?MX%U
ztWqQ9%VYyr$o*11gDQ!tFiWYqh?hZ0Q7`H9mdW?!rMZ&%XSRE^W(ehEJVob2vH_#V
z!p}1HC;RnFZ{a_)0Dr01{;mahY71lJF*1StPu<Ir6ynmK;PFe};N<%ML}2`%$qV?Y
zYXB6!c?J;3FQo-6S_JSh@xyj9r_M2Q=WXJjx=@^RQ(H>=yop0Qepl#S;1c*ZruWA}
zC_Y4ChEbyDm__Rfbd?SBLPv7_&wA-Ur?_!@f7q{vzBF!hF_#Mad_tk4h;GWq!excY
zS@CsVrA<6AO1L9a<&W5%T%Y{4v;{7JzFN=_<W<nWWJ2?ULHxh@Rs5&e0fleaK_S1?
z<-+l&lD3UiJkq?@L)=UpA3x&mH%a@cyG-v*{L%QAq|bT1L`7|3tPSHrS$oSdjkVbZ
zd)*J$f_R471ROKg<*PD^{Wq+~BE8*R7$ppS&iPSmOk+CbC0i(5!FwD0@Hy`YyDs~v
zm2{X=*1-?~$05bNv$5_8f$LWy1E|lUGMku8$zF*>H;=bQ6r~l%7pE0xu;?oBJ{OJ=
zelDyKX-Ur<V+z}zIcH}dt=DH3scSo~U=YOe?hW47L&4CE{03IMc!F1vu}F(6iq$yg
zCwufuC+F+c65!+fpRxB(U0(?FU$Y$AENyd9r0PNM`gHp*w+^GprDkwYr3Lp`T$Gd)
zVJ61oAkfZ=zK<8-B3@HTt(r$I$0t=vB3Fyi0|VjB5LTcOpRv1qDkW^RKtJ)0PA&xh
zlGFb)srgG4_jjbGy+uCqt`1E*6Szdk{>{iDm?yma6PEgqe4}3*qx}Dh3i$>GaF3Au
zr+WlK)Su$<OQ!~M@~ZemhXAF2OAP*<n*<2#myjsN3K6!6<4gJEm?;yFaPpg$2TBF8
zv8mn@M^(R-JGsAG9n1Fik=5An!G+c{9VXb_2Os8xQ{=J8bBKMlx=4o7oC~^2qTlRm
zVvg<5RJ{%Bh1Hl9urN7Xyd%<A!gy1VcIYfZWT5{p{GF!Xr;EWKbq@AzpU2NUx+#*A
zH$%5x!<g+c6HxC6ABgk^td{LI{*0fn+^_wg(@V~O?~eR$;S2rcRj{z45Mc;lV){jS
zpv|x-r2_#1AWY1T#OpR%w@fDYg|RfBZdZ>vl7<sW<vDpTA19M=-c6pNB%-_X$~2vv
zCldFC%e0ERO_g7N#0-sN4|`hzZ{a<C`-XSVni*N%?d?UFDCgI*#8nsF(<v{+oS!8#
zL2>-_Qz7=BIY)ovRs$mNE9Xcm26?59FAbPG42<lFV-ilpSXj;j5_f|y8jZ7Xlb<!(
zr!{goy9C0N)NRO#xnMx!PD4-JKvU4&0M<67%p}WaqOE~SI}G+FnweJ_q<(3Ar=ye@
z#ikX1^stqjP;9PWH71_+rsJI;<V`a56;`XC@Zo>v?fe}Ge{J=|h~S2*;)fn&n#JA{
z#$F|CqcAgNyspRg1n-@vUqV9;lW&4)YQ3COAQir8<P(qdb%AT(&!H0#S_UpobLM`q
ze<G}7NLstVrTDR<in+`9fT0{t^1ToDTpnzc-(<t-D0%x_vzDgOIX~VA!0_=7+l8NC
z`2SlE1oCVBA$;ozE*gYb`d}|gUtj6%$kGfYI$mW3u@HHp)XT1M5BmE~pHmYiS!R?l
zh=@bf;N3_E9e){LRBeqfTI2-=-o@x^ZgHtX?1S5dGNv~zUtW)LIjf5qDT5e<mdGcz
z5%1VU-M?b+Mm@=q8piukD@)uIM?YB46h|Uh+GO5e1uYBsS0lK~=AwOR4STvmoLdi(
z^Agr#@#YuF7216qFbcN)rK6vK{%beqKf4%#0$$HAUN=gq2`ZyY;P)Mn$40)LzRqAQ
z&LobG9)W*L{LEPv-+<9s+gsOX@3v2K#PmKPYc#Mp5?i-WdI9Ykj~>4|=8$V2N8jo-
zZ6M5rC^-X{ImiCV?lpzu7t%#v{rdsa7(&ncuX!bw=Q@NO5+{x4&MIeXJ{C5fDeN#B
zfK%a^XG>>zs?dE76}{U(m_h2Zkaaw!u_|MCmuNq{Sk%+e=7RCImc>YEpA;!@IYSs$
zl=U{7qp>^G8_wuTJ=$f^CD~<=blhu@)P1jLW%G&Jr}|5KlXd}5EhAPTk$0I7+vFH7
zO_}EnR0_NV(UqXI4Ve01oNqHdVAJz_MsZU-#%a)xLZrbmP8#m|){jDozSplF7mN~9
z9ZvY;E}y(s{!26U|LQX4*W!uqm;UT;Z<Kpxq5P^S`E(au(d`NI2?+4X#HrfQULo{U
zJ57{6_Cy7#qX5$7#%rR$pFr@tzjGQii42ECM$oxf0zq#lZ}5R|MN12WDey8Jq<`X{
z@(J>R;4ml<PAKW(==8N7X{PZ3F#<rK0vw8br0#eUu_p)w)cpR&9cgv|A()atKpil=
zD~RtTbl?wf_D?Fl?g{RxgUW}rpn+arxgH2;{KIQ@LuB7?ufRPsmkkxF&)2IRb6O)v
zThm_a>s3_5Fh-c5*0pm5hWULkm+XTKi?8iGO@wOXXu7NHr1C?~2v|X~ygsXG$;K(N
zi#;7B-0vUFcNEBZy}Z6Olk@n6*d;%n`EBElnJ>pHJhKeVU$&M@yLmnw!5x~9ci!uW
z?ZqE|_P9fQ<(T0(axmac%$>~D7`l#vv4Kg)<}aQHvpFA*uTv`RWrq(pR_Q8EW}DSv
z<g9Kv74c_BXImVGWY6%z0<N#4eR*>1k@tzaetJcc|Ahcc(tyaDL5KOePh?l+1=+OP
zMaX0uq{L#a*{1^&g+7cH;;MH^u<NXH%Ot83%7^*#poWD8LD|O=o~<qQUcD359DW?%
zm*}E;@7zwo#r`t7TN|?VL&5zHbcq}D*}S+oXDrR?x+D6_GNn0$_B+=`VnWmlM~T9b
z4QxNaldW|K<uGK8i34mba}{$}b=F3bs^wng7<O3;5X#S9ZWa-J;;ee_b<P>!e{l}_
zL4ucgP}MG>{pEFJLyUcYdq^Q&V*W*Kypx7f`|ImS63@rYFEg43MsJVrEd_YS9gpfW
z&yJ7#`?x~oHS|km-3K#QX4zt8KQRQ<`E60gdrylzqH$noYC4+Oom_5j{@}lF8eNg|
zXu+&0d=sZ&<-Oh(<LI^H#H(NIpJS>tYwd-;S+jp2^Y&tn_wcU$^RtV}hF?gn25SBF
z_FH!cUdv5*Nj2~vd}8>rEav|?PHeaTli0D-yXMr;(DqSb)*)K+Jnau9u7z%KjQNQ!
z-tAV_Bb9CwS1je(u?lT`onl1SkyzJ<V0!JFX~+xi8pMue7D~v&Z*q^XxO;iea;Zbc
zZpN{l2QYo;QzU_NC$w6al7OX>BJu1Cv(iUCJ<o)2u?i0*njP9aG+!{yz8l5xm7Se%
z-6KZqw*ob0oNs;XcVmI=n`0=P5aO_o=Sb2cZ7Og{{$xwX%Du`)S3_PxDb-w~HOZx(
z&6Im)Tutz%W{HQ+`njZ#esVT{At=yTwAsTq`$7NW>fq?LWnQ`H6$OF=m?*0)*PuHa
zu9IU7(<Vhg$>@ib!a{-NF0q$oJ6lE@8OD>@UsRtY(P()d2n{+4$!i2oMV>ROm$kj&
z*jJi^cvA&k^~s#=@huQD)<1gc=)Zno9MEtyW8A#=ae7+pc+WV%=U}8euKBT_`vLXQ
zB=gaiV&edjb59QB*5I_4WNyHYuiTB1=6O5pCB;cj`{W~if@J0GeiFfSFO=|iqY~}Z
zE;3`Nbm@ZE*m|fDSYS4(Xa18mG*H10M{%>xEVC%eg;;k!;?iA#-;yNKU61D*MTsrI
z>f7zPZI2RyIi9!{&wY@ZDMduJ&zv&h$?--pB%Fl7Yj+LS1y_kPnDbLu$&9xqE+&d=
zQeRSLmbuWHe8!cf8)RNDlkM9s5jWKPfMuVYR6m|9#Q&`eOnWLGraoXhSWvnsf7=A^
zU%PTi%DbVjy12Qvh3Nh~Mt@M5JEn+XKHIaT0}>%w7rg%G_r(~f3)$aoS!0X9U-p?1
zX5^{i&*3?-RHED{uo=aek<&xXY%wn6&WzVY3BTl6iX(Vi9VMK;iOh~XT@@u<)JlGs
zgC!?5T1teaqe79XYWwnYflY>Ym-3u#`tt8j$&AN)Bj5x1>#1}uKA{K`W~F0++H%Li
zE5ct)UhX9^`_B#KAF^d4s?9Sed+_aMrcxeHwHo{F^TqDw@yDCvTEAPOvW7#?(AMu4
zLLSz4mb-KiJ!~Meg%4PD!>8m3kUefB=`kZaI1jd%JN@2a*B$1aq~2QYJV>;~EMNa(
zbF=g-jfIr19(<9Hf`SX--kb!a*53PCv2b-Hi7nwism)Tlt#<QZq4<GPaK}ws$C-sh
znFDOYFM@&}!)|0TXCj<9>74Ii^1Ey)xXQQi7vgVSmZZhnVfFUbimf9#O49Rt%^(tO
z+W%CI2j;syXIfmWGaJFcR}o)ooc_X~XHEwoi#cy_l$PzE#2h==B2(s%z!Vjhk?m3w
zgOT64jhyFK0-qIvB6e;XPnIyGJDHW8O!Ly&6C^zIFP^yrtV`({;}c$yBMw=*yFoDq
zN9u#>v6i{&vAgRqpMA3ani0s?b?t7^(Vl~=Qp+AWZIeAi66U2#JL#}VN=vBO1rdWW
zWnNjkL^UP-P{75#ch99Ws29Ev*0xU3mX_^sI!*TE&E5ALl8vf8#(ts$1Enq<F>~dE
z3kmsJI!?LJSntU~;~`-7_vqFOPT+y8Lv9DrzZGK*uhvs9m^m*Lz$ASP5lqCQC+kl?
zR?g^JQy5JffM21{P8e7bEUb)5h;=KHo`wi+Hh1B5*_swB75J<K71y%NBIxAh)$bhk
zx?X&{5UCy0J+o#;qf}w#ifr-mCGUkRH?n0C7;@(Dym8U|Tv@$)#4ae!D=hb)-@dLl
z^F)PENb+ImqM%^V9nZ+LDtVphG>VvT=bo}@@~aOk1~g8nRClZh3dX0ihhb&Et}RI9
z<)$H;W-+=~V0542y0>k`r<QnD>A-YxQP%MN7ZMtZxoo$1`(nuU!Xa}yv58+Y_LVa(
zdprwb#RyC<0@?&n+LpAgU-CNP{LGO2&V{0|EdH34`ta;z0<ok?LKRm^$AgVBHNjW+
z)>7$|V`EJy1YXZP(IFH$-`m&7tfV}9gprT1k+XII*3zzFOD5r0tYunPXrv^O`*BNU
zg5~%PoAq+WC@b<g>)LItw0$9ano%CCoQLxZt)x^R#yv-4G4-N)UbyLaqJ*Sm(@+f~
z)aD#~Hs|eQrUV7EQ$ifb9;-GDDBtN33nV?z>#M07xcBPFW%Fr!(EU*1APt0XVNSSC
zA|h9CKyn(DPGxj3BrHoiW+iP}WhJR<Y6G3H+`uZI_Poe^Fg*>?#2SqZI-bx2no-v{
z@pD}9D#}HNng02D?5lSvH(z?9g!>cjR~COA!2BiIt7_*C(yh+xb$IWDX=R1e!m?!g
z;|pHuwtSw!SAnUID%_;3d6$+siO?IK?{tlmd{qhUAUfi&PsV$9NF#e#0x)mm>%_Bd
zQ`&-4;lXX1{dc^yZg{J0Wf^^Np<Cd6j1(Ug;AAI>4D;y-*6NCZo7S<?hLZhVBr0Uo
zH!^u_KHJf)TL3$|{$X0}gL69FdLn!fL0RL)AL|Ll6KeGQn7S!9Hv(UX#Eh8XnVY{V
zIQ5ikX>+RHamjSUgf8+uduO33e-`5Avlw1958jG+1v&;Fsi??2;80433uk1_&FcuB
zh^K`V^Vg^K<l49~zn^-+3*4!+G_prEmDl|9i;ZU^m;oV?VH_(*SsTJ=#HN>THy6y$
z*`+2R8BqAP^G4IT?WXmutoOU=(8uc-bR)azbR+vysNJ>b`o@n1i)f8y;_P`gt`B}R
zeSfpPSe?nv$GKB-MVkHbvNY7)PUmSEB>~5qoqO1$+xZQpY6-#iJX3UDcJ)sm>F{tk
z6FzGENIU<Ca)~j<p2wu;5v5a2lY@)Jdptq&_jvHCWq{r`+6#t&_s=xOC=S$RZH1jw
zH(M26Ob0ZUBwya6>s27>omF>S=aHhBO_?mjb`NY3Tf`$7TU?X%VZZFx(dbX{>26xh
zraBAyMQktj%k^Ze73ETA{hpjPj#|5|c(Gi1?)_}a*~IPeN|K!->@BS?IGlrChMs%Q
zUhG?*N1Wc2Y!*Iy&g?<?!wa+BTdiuun;X;no~p9`d(M7$c!u`7A7siJM82%qj`S#s
z?-)v9NJ+WtskW4a9m6NP!s)fCk*E+aV|;Pvd~M6fVwl{-Y;2-Yl&Yw%8Bvt*kiJo*
zs!>+AoGzPQ0|jER@ZdIe=lD=8noa~cMcxd^AnGNR?1M0w%D12P@Nj58mvpYdrOxea
zR6EHp+j!2h_quX<Fsjczd)UOKw~u)kCO6B!RyC{6CU<V%@>z_m@k4L!%f@95cJ;5E
z*|w@?U8YE1KQj#iUfLZ+^Wsn%KQbWPivv2y?AxaYj7T#*CwbF+SQq~>V#-L&VK!}V
ztJ~h7a?88VkBGY&FK-qf=ZiI;R$+ZKpO%4?Kl=T<Q|AwitV)4jbF|FpDKPI%b{yM#
zb^JRm6R*`2sb%W!@@Z`n-^Jv*w~_YvdQ1Mjrn0o;d=A?@=44-<dsRkkqgzkUt5?#`
zwOZITaWs0S`A>hK>)pJ}zT9kev&mDp-6H{_KWMBT^MreG&+4W>XLD2i(^T#SL}lV=
zatt+4W0|7YK%>*AJ^l9j(j&VXh`xb~8`oa=Q$_ul`{4=P=ZcHCqwwEt45USmxX7Mg
z+_O-SsYBen96|HeS<Rw|`>Nl#-TW~2BRJtz8O@}$M>b&FSF4XKG_#7>E<E+8j*_;}
z%P%vQe8)ZSB&^gaDr`tFGC&a%&^a`kOgJ+@vDH}DkQ=3_!IIVZUTGd5tfUrz(+zbj
z>8p{JbhzXW$`C4gAnMJ{8{4a>57X@~SbF64HtkVKGu^XITd2v+8y}}PV=oyi>61k~
zd<+;2X^%Mg2CsX|X3h5{7e_ouo;{~!_0D!sb@4csUpAQAwuE~8l|>(gwVb^p34gt^
z#-P575ru)KReJv3{fiRqY>d<c8!+jt^XKheazMGgj1)(F9zOpxK<1uIm2#dQxjfwG
zh}Aa~IEwm%`fac0l&s$6KNTo!)X~Rhy!*PE(!!H|(n>QAnj{)e8;>O^aJDD<HFl)v
z+?hcmPVvUhb=-b@S<TSuCE;idq1o2w7BUY)Lb;=37wc%t(%^?@yYWe)^wIiJ#y7L1
z80dQNHO{We(R823PqbOTcUwlBsr4y2&jV;U+Bvb5w=(RZTrJM9*WP#Wq{WgF3aJGV
zT+-DX#fZK<l^NWg=qb#|3byt<i-(i~s<ef0^Dcu63H6Efyi<*Wx2QUC&!)Y4uQs&M
zi5r%I(PhgdS?)(?%0-S@bxzz$6BL~qi^dbx<c*5Pr{*mjv2g>tzfv~5*jZYl$kaY6
zItq;dlYMhhDTot&?z6e44fMfIilcX$_!NO%0yXV6jm&G}l<pfM4^^WIK@~9SeXDco
zT^V2E{g>H!`6nNi<TR<x8Q<KGd`2Rvj;Ag;9&WAwV0A+|t+Ki-?^;6&e64#>s$Ll_
z_-cF}BVTN7?4Zm%cTrzOp#+b67Fe}+P$q?KSgK?Zzwa_nt~2m0!5E+3b-%PIC5SxO
z<j@XU$a<$Q&#G6;^aAgr*L*LECLTK8(eZAGcCM9onlSLeW5kOQxTa)f7jw<`h5mk{
z{VS{<Ja4<+mnEY9+UcJ-Ft^Oh_0Rb5zVW`<h=2TC(>zmc*zT%)ZjmQvdh?3(+uUSN
zflCJs`c(#Piq7^eYCN9_s-+>OF7~ySS0hI9uTPrR((@)EZqnWux4(L9_Nb(sj@#%X
z-0GQ5RkrdY*k_1-ydJ{vHWzR%gVGVDPq?3@6>-}ZJ^KVD;U)K~eY0G|H3HQzu-8rQ
zcmc5LbFAgTIxeN^yL`}C&$JgxT}9es8V@Yv5L~Spk4kFqSY8rHL)`pWU(Nm;znFc2
z%XySjEe#PB+?KhMD&yLVh`OYpt4>aymTYq<U#!@9TaC9fBf<ZvrXr@DnvjS!hu`_g
zHd&!)spERH?b}2MLt}Y%yG4bzpsm8T>Z7WrJDK%Og9hPm^{kRso3<q0Y%miJ?PnGv
zVu47O`>JhO_w-!}yB8nb5>!Q-zr>_F9iR;IEccLEuWw)C?3jWH4-uc08f@ghr;<TK
zN$LH7L{QF6fm8$APdORnZr0B?=mx9$%-?}eb<Sgmj&X@u9z0IdxIQyWgK0|lV93%d
z@>$1pmFo!ezNzeG+WsrhG9EAQ%siGL%x3wXx~2}g@iXjMn}Q8Igb*KL=8tDsP5h0$
z88T(({T8u>-1uj(ozvd*yltGtbtT1ka>Ut6CK<B7`@pSe;#CcdZuPDg9oCA`hI!I)
zeUt2>0;xJ{BeA#-yP7P@%Q5+zUAXVA`AaWz*C6aF8pmoWhmDM#(T0%^(B4n0qMG>N
zAvS33R1VSX$`4S1`iygd%azQ5%iJ>y9>ELy@P{qi99x#5%MWU7xGEakJn>tABEqgk
zWF5byaRn#PhtvK#O5KnMs*7VGWk1UOMdo$1(QEOK1;lr>^<WfpM%=!>ogFJcAFc@z
z>UGoxSyM4<Ys}~G5wB3zj3$zl1h|)gXKNy9S5caVJhynMgdLfs{o8itZax_<0xa4#
zv<cp&P-k>b)3t^yewLFGU^BHUvvnyKT?tvUeDOvP95|+2iYC6>b0{}+qcaA!hta|6
zx{kg+ONc;VzIu;%;|aRo?`#Tg+17psqoB&^+avAj<`pKw!+i}rM4>WndeUoMTJk)L
zp8$OWVFwBFP5bmFlc2~WUwS%7wYkF;RUDPl+ZL3lVxEdKJKd@I3BR~Dfr&_5)Vo0U
zbLfiw=AbxRj*L8MZ}+@o5r2tzBPllx@MH3DlGSrnDg!fk55|;X(sX%L>{Vca>dDi-
z6v2uGaJwnsn_;%?v20}6G-VgdI_fr(ijZAMjk_RDld}M;7`PW`BO1fQnjY7rGJyt`
zSX{o-Lz`ai0wX)r97kw~3y%wHAY9OhLU7*-#Y_6<?|I^+pAl6X?_nph5hX~>NsusX
zgtsqb&F8J5#5N1<6FhH`xhFie2cyaH6s3&IW@gh{mAct=1|H9aRZ%~1FDODZW<RH8
z(*2EXL?f>Vj2@89#MTwSJc3e1pI^LC6ju#eEoIzRuCWmgB&9aM<oD3yBOWNe+Sjq^
zwI5G4P!NG!9npl-%B_}<%B7fQ`!IcPG}cX9Kj~#)77%B@Oe~kz`I!8#XGn&M1xGrz
z)@fZDQ}e%$*k0+5nE@E^5sL{QR&k=8XVQYV0p}sLL^NOxVPsD;jL;HS&9O>VnHWU9
zNtN43nl@8S9g=VfYe6V$LWpxh`SUQV6I;672*}6U8<?TwG%1g1I$-e=pll%?A7)1L
zx4_Uwga#xhQ8Ql#SOYdQM_jJx^QWvO{XQ3yc*(|<qYllX@TDszPQ_ZN5`O8sVnLD+
zE*=g%n^i@k7%Zo`FP%sagdlWl9<7}eUXZ9}`jyf&=0-o)QjINFkLcAN6blB1imsQ^
z(koeWd8ospqGAvhB!-nH2Ntd96@Fx1;b^)|q#`hs_bCQMXOJ!lLC5xdTVKca`r0$C
z$fyrOukzPX2-~%Ahd&_Rh}-2{L=g#UKO+N~8B05HbZ(Kfwkel(1urLrU1-@>%b6f+
zx#h#^ZG-YiW&2%U8{pRhkI1=OwlShC=5?SmW0CYUfy)IipKWP)bOMht%?>Ok?aU6$
zmM7qn$XTKci{iSJ%x4_S-`$K_qP^zYj}fJVc#Xw2ewlDDa;acNZ?_Ji1+bi_>~<sb
z!C1WCYEP5ttODpjkH@xU<Jw2A_vT(?Obs9CuNg$4y*DSDa$fIxEOH^UUPQaPyh6lE
z3$FT1R$^NH_K<gc*N1b{X4$F9qm@R+RZFq!Z0~Dh>^!0uow`UDS9B-L4kMnJ&osnN
z(GSxL?Hk3~$Guv5DKPbMo65b8(58L4^5qVOy-)Pb(6idcLPp#OUBBgftM2DqefpdY
zU$O<>UWQz(5~?_Sw<S`H``-KzqnAyc`n`Fd9ATwj?6nMBhd4+$WqEE@A-DqLdN-WN
zqe7soG^-#==26$K+?D~YU5qo|3pmTw+OE3vhcT`Y%WL<YKc>F_K*QNUXh>%l1Np$%
zjiPV<!LIk#o+#RZ%j1q!l(JO<oWLzD!n5U47w(sLUdK5P)n%X@Gc-@ee@p4z<aIu;
zq^b5i*&yXzM7_fOl%V}P!K?C85`LS~dls{tBp6$22RFO74%kZKXkL&Q1@GIWPWV0L
z#6jJ1PBY>fTII?WRVH<5^tVRcTK8HHe^Qg5D@w!kT0xG1vc=VJFeyTOpN18Byzzwn
zPMi*TqYuvf^lXbi+w1qYMCr?f@OPBt7>L-p4aK6w_v=9QM$K`ceb349W|__P)VWd1
zRaL2+#fj=`xmlq7eZn@{Iev-F${t!`LDbF4&?gR#xO3A5z8`3NBv3cEloF?D$^=mS
z8f+E@w%T!X;}pnUoA{XmMaQU1hwP0}{d`}j<E1p`q3zc}Ki6Sua^%JdBtP3CE)Epk
zkX>K5bX6B;$ie7(XrRM{J9iVm#20Jx$Qjk|uFmRZ@iarC4NxvuJ>|Hrc$#3P3YDsr
zZMUC%cwkjs;fJV8y=(&yJ*&d)sZ=)DO~@-5iMiy3pZ2`H&VtXzjAf!Aj=-X}5r0~8
zlb9AgG4m=N0V+P-9p{8})bo^f#VZPD<UiHP^Gb+cG1ahDXe5upE~Jb{h>xt+!oj&S
zVI@u`jaoeO^wNiCRWZk^=Vjr*P_3g7Gz}a!cDZ?sJl=eZSR&oCFZ2p7C@}SwkwNm8
zn2w(F<<BI&c4CQOlboq;A%59?J-w7zvr3UR`TV_^t(3uO3r^;br9R2iOrBxZi4hVb
z<?Sw|9b!!M6s^zh7N9#N>wDoSHOtOwlfKX~<LN|&kC)J4!HhQ>OHtAtDN@d)#xZ!c
zZd`H1<YI_H96m!!R<YpcGIsUjK3m}nB}7T)VK{dNRaZ3rYK8BI)P21&(>*(MR4(Id
zjR{oA{T38l3<4uYn-027g4-pE=2WQ7@uqu<R;?S(O65cmvl#`+5ira%TM+p?Uydr!
zp65uDm(<4dNo7kKk62L-_S?V+EXE7hB2Z>Brtd^J;V*jIBE&@#BtQ(~)=8A0aXQm*
zY@O>Jk5Kj$El@YE%!vf~f?<m$&!5Op2l;APp>FU^y@djGMT>0Xl5>@@;6Psu%phN3
z@3-rf=sj4~LA_3}aRp)W9sx_@p5kr}MG!+DGgJE&qZK(8j147f^o{i<FFg7{U+|gv
zOeu}0;+<@_9D6~x;yx#TP!ZGwF}yH@W^|6}vBo+^fnwt%mu`TY1V7!Cn@(DU_TIvZ
z2HHlKpKls|fEG<;r4T=PO!XdR+ns8hOoEDG=OI`*o1BV4vHF6f>3mBH4$vAUx?G?^
z6Qv7_v?j{Q**gbdQvqq>uI`vzFY)!(jxr|SxrNkWj80&>okIeuK3h5YyC4=->eYQV
zh#IgoHD7~cVpmKu26WO3lIe}-*rR|MMDfu1mN{ThF@!R5v`b^i23!kTSE3GGUoUdN
z_HWw-`;(B!%f_MRlTwW*0`p>3J{-cC4(aJOW9+rLBqI+1iJ9k%zeKYvam$w&dX+fI
z+w`pz0FKm#d-SjkN$^ZY>s$QhlJF#FH=XbeC6dsMt(cGgQ<azZK=Rdf+}Ae;LxI+T
z8|xa3w}|qE2zJzsC3D1^*n`EJjuq8qfIT6m(95MNEdGXVz&t!QlpMr>o_Ha_$M7Bc
z$v%V`wQHjGVC9uOU3(2f+bhCH+w;bw!y>`EMHHf&G)eW5-`Q6K-PyO#owO4xuN>#4
zD!A!M?rsr&{_e2q6&zJxVps;g@g`px0}Wdfdv)%vCwP8mlRDP%#`_EHKKj#BOs(r_
zeZkY_EZX)F)79dmCcX`qg$X}Vfkl&aRec+JUM((Mh*Cj2oJEgw&d&<kl|qO3#_Ozo
z2%fC&6W72{9-1bA08DJMJIB1a3S;z5RznGOv^%4SozuHS+*FW=(PmwbExwCM(LwCX
zu@O7rBfI9;PwYE=PbA&FD{Shroxti=&Gh*Ql^4AohP|yphP?};k4oGI&lgYm1b%*W
z#j9<tkkI3Xi2TB<RC_RXFjkOdf_lYa@;RqgXD8xj*Q<WcIQeUt-4{sg)YS!I&zTZt
zB%U9{%(I_$w(3Yan|JB09zT11;IJ$<<(`CJS8o5qmy)-rA3I}n>56x*RoeBRC1~%K
zeI*ca1FqbQk&;-aeV4!_>&6AT8gy(3o$M$T>n00T90S3*RwL&ZA3gMleT6*^qCUlK
zvk^ItPG5OC<aenQk(ETnTyv3auk0w`OCVbvX}O1ZgIipHWP_VyUU(bsI9}7mwZ<^*
zuvKhtm9iq#P5N=^GZT$rH8-tLo(z@eqppp~^XLqKs+13vH|d4Y*2XTUlXtARiv^DA
z&vo)lEbO<`neFT@GG2GglQ!>Mac1W;+o1sL>80$h(5_sQ+vU-GCS%G(7~(7I`G#WW
zrSO+(IP6*)Mf=KI8udoeKIq%L(hc-U9JWA5x&yuffj+{B2X!K%a585fJP|VOVfvJh
z%yy1>TBvG1P+sN1?CHYWP&FUg2e^Unrk~!QG(`F;c(y7<#^TF+mH2mZ8<RA>T9&pj
zGKINUY7ihR6;|LU##m7)X?5*DK5D9n(ds1k_0BoKJl6&LoK-}zic6TJx*+3NENyPK
zz01a2_5_9p+*spTsRUEBEUC@;uJJ$2-@L8|yI@;hPTd~n!}<=!KWIc<fiJ#H{=yga
zVeRYu%G~WJeesj2SgBkFw=AvG5EK-U6?~;gm>kOw@u3yrqv*>H^(L(uy7`3AI<*0?
z@5GlO9OcZWuchW00g84im*oSO*QK@3Q*oX6JKOsr+V=NWYuLFeL~smHAARcbPDR~P
zMqXPd#x>icl#Iz_U|!x~5zT(tc+0yp8}+mDrYg!@wxeH`IDy@rr#vr3aN1BGr%)do
zvVvEohzJ5#)(ssPgOIn~7>`+7Ma6gEdPE2`pxaB7IV0l*A`{bMtuuK(>6(gj*$kvv
zQybUSXBb){&qO0AJRXDK8WQfeAEP+Mv8H(Kj78FIpu8_o%nK>e3@BDI4_b~v#{mT3
z1kF(7)A#dvD7lTV1+ZGDly9~Lze|AOtNJuRyDHc^XO?ShP;ARxuAEE#AeUqe_PM=;
zYExcZWoov8oW`KhbTAe8#0D3+!p>?Hw{z|+PdcIEm@?sd_Z@Sd=$wzhs}qt+R>UGH
z)t>0eOIoJ}-cfS;7PvQSHs*CbGL*97`gK<*l~)pev`0~RTj=8Y11oo<$F&I;;+#a3
zS^!7JCBMk2fmCpYA)$4m9&JJB0WS8?0M{3s%s~$+GGwPNpi_ci;HNB`0r9%*vI~lJ
z^yM{JIt0;3=MR<cBI7d?R{*F<?<(7MbeTyez4ocB`LZ<>Oh*n)aK%7=pst7^truFj
zY@LZ*!->Z<`5CVT>18o(xnW>BOitrk8$Eg>1iS@;O;CuAwU~YZ-Wh3~iT8f|I;3{(
z6-x`A?DcNo+z`}-#c}x<^{Sn)q`HcB;u=L$Q@cp{@cSM{BDl!LfM<m25%2;Bqj<+3
zYy9Ax@2Xs*$DJ`D_0=bFgF^N}d;Sf~+2$jhaAzz#GL<&M)TNA8NHPa+ZHvL(_mE`M
zOdoh{i~JIp_x+h?bRj(5x-@KfL@h*Rt@OIaT3uDPU9)#8?1I#;UnEbg!syytoFpku
zInZZRj)PV>ZIV)0;@t;RaZg6hC&j&9*VTT0jq24Tg)ZlyI<7X)y5C#D2@0xmUIxNv
z_eNXA2(p)-Q(=!3XynrKg$HZBX;{-Zo1OAp)?yYj!=z#F{mh+(N*ryk#J)Nl?K*ze
z=?l6me1tgstyT~74ck6>X4%(=#x-{aJu!Kth3m0)MSuT`E0*!c)NQmgLAoL`(MAy{
zy5$Y1(U&XeDOC(2TF&`~O@H7N?rI2G<>ridTcCS*-XKn+$TOx>qmF)-&~ajk?13zI
ztcps*S+$z6_b=%r7uDvDl9?qIo$E*WZuOK+a2`&Uq564aFkmb-8Q#R!i(9<)+1XDa
zGk2@*P2DqDwL1D7i==XAMDSZq19W$Sa`hUc$>7aJl~NlKqZ}`cHwDKA^~1~Z7UEkr
z-W(4M6xcLZBSgZR&Z^y~kIUNP(_W>zhN8n1^=4XrdeQ2Eu1GTeB#p}jT@ji^ytp78
zCbYMu20;fj_jTmUX`b9BzHc<y+v6F-Rr74;we#q^r412w`R21E_vwk2w`kDpQdXUV
zjQt|tOjpZJZ^dD!tTytkH5fPT9088LO2DmZdQi#zt>G4LSsD{m_z@$$g630j^x>5n
zBOCGp`Shi@I!$f9{kq1PH%@d7t>nqQ_q=LCw6!f4GI8Qb(r)q+3DY&MeT>s5ly!v7
z#NvA6aB|8Pk3PAP6ccApMUt<nz0N2z-A4Q1{#5ECQ{AVoH+VcRU2-J3;&v}dt0B}l
z<=u7&E6)A<Q`2#m4e_4G8QqSWooXPj(DvJ*^~`;APsae0t4v-v|6<`Fj6IIVkF#uH
z_PHkty6YiHW!*hlF<pZ<ljgemor~A++CS8KkSXu+CC9KS<p2|BtWhA=1hvKp`)F}`
z)JM`PWN|A5q}?z!_cCIr^k^5qV^R8vQ1>NKzsSlFd`}T-T$FMh8EVFfEh($^MXw<b
z8I*F?5bUWf?}RgT*%5Y5jUk+#DWyAZn`t_$Nehokccg$J3ttmi0gq$!Zw!5Q5o_;Q
zjO=3|sOXk_TB4jDA##qB)vP}boq*VfshVk<P}B^I*vB=Iyn{Z+_?3o?<WgGCjbR~M
zW3B4aosVnxBSco#rfXl8?kHMGEID6o#%dX(d+VP2Nn{olWs6Z^4aGB2Q65KW8OtQT
zV8TNv9E`fHXob4u{E3GAJnnFYMr1FTP736UHq4ubQPIt8<ccyJazzf+r}2;h)d#S!
zmz)=dc^1@L=#`C`685%dk7qF|0qe9<XJ3tl8sy_T#j9z|5Wl~3wr4C+qQ^*9fM70F
zeA}c~Q#WTCtwXKig9o?>`|`GpL4@r1w0aMH-$UGyY02dGbq_4Ww~O#EZ&Mu&)W+dg
z30Nb@Ct(hK(`X~S@fcM@-JsB==Oc#Fbysg-R9Qo@N2aI3;S+XfBO!Gy7M`DN>CW+i
zfcXV2t+y9pR0%NSCQM#@B)+YP$>`(zDJ7GM>y}7_wBijC=c9V`DgjH19{QN(ZX11|
zKS)M;0~-gw(wYg+tiRhwEO-pXhj$mE%CHJ1#Lhth;K5i*%B};M^|xKXyCcI$!?djl
zVEd~7tWVIm)DohF&$8E#j>{@Sg#0S8Pvv+>9cY_NU4#0Dq+o+DV2+!`iENjmU9jpm
z6DFCI>nl$M+JSw1;Z|+Wc0lkNCx{yhl?gzlbLma24{F~SJ}39fJI(ls8D=MyTcx<f
z)dFYK@zl)4swlud7pl#+({n=zs<HC!W7fhPyD)03p`12xO;dsSIClMJCU{XAJaDYR
z?5O*>jiM?{HvPfYH;e9S;m>55ZEIQp0LdijmyKnp)joB4GcGQIZ2>r4J@f&w6-JHL
z>z51~bx5kL&*E9F9(`=w$=Dfl_W3g8<m0lMJ@$C#$#Gl<mW5{-)#MgiKc`ASz`^qK
zEwRVh49@p@9;*0e9!z^3tcSnX36LP+JTEzj;&>7XCL8d1(z-U_%kH|4ScBhg%6k>D
z568&_VAw5@1+k~B#WxH3O1Hx1KGji#gu73Qzj)3#(^8aC;o`lly_De*?r5zC+zCEC
zve`laN^JU>A5=^|Rwrjkpn3PGd++&XW4YKl9qd<ey;cKTSA2~>R_^Iym(BPbREqf@
zZCQ)LR+-kcNbF#rP|wk9h0@8!aH&Y|(%pi+xGQogJD}IMa=hd7BWKUgv&Z|6#|!-d
z-UoX;0ejZn{k7Y^ey)X&tj9e+%D>D`a($I8F$_Oct9tkrP`Oq<A1Nmbsew1_MeDOo
z=W;lFd_UGcu^V>ul|&4W9@7eW8PDtq-tl>2%F`e4rmZ~3D1%2|>=K8HUfpVOytmC{
z6N5WlUYJT6&kkp2W$xC~+QR`S+bPG(PadSX&d%4NDrD3<m8I{I@e?+pF$mm?*#+z?
zn~!uDVdxas_4t}TVBIl@@|Y1@J!4e+xmN7XBlqhre&X{fUXe}iJ83-5Bg5mq7}uW9
zFl`n(*{}JIC&3&^OdC{sh9{J|{L8;!Ithr}`$=f&>32^MFh3N6e2xt3<oo^WA0XiL
zYl=%`#3={@B({{YbOrKNg7}a@sX%%t-K%;a2X{NWum6C7n9^@ismL!g)c+7PdP?c{
z%%8ud@dQFm-GMk(Al;b^!UY)pYgZmc)8ER!XXg6-gOA^)b{eSqcP96h5|GdB+bn;J
zT1C=5SvL^FdP4ZykKZYO2g*qv%9D(x$Us@-?13>A5J0}q-=H~(dOb}~`uoQ|xIE^j
zkc&Vf9b`((lN?rmNSz5h1D!k+6$SETO8oD3k`kxw{s6J3i=`DV7{tei3;yj7h@X!S
z$R}n6`mGN6@h|cZ$l<p-2#^fs`#LCu51HxkyLM2B0P>T-@9Mx{r~m-WcXdEt_@C>*
zU_Smoj|T?x3;b9o00%hhyYawaI1tPHr#j$ErXTu3Ajq6}-}eQ+$ppfe|5V3|j9~s#
z9gGk7h~rOnd;&kV<NtAO5cp5#0*C+DS5WYWIYYq$KhOmL1U}07elAeRj~GCq$mB=g
zwFCBu|A##SV!nUa4=5a&S?;^OVBV8Z`ajjde&iQk{vY><R}dMg|9xK=FzFxaV1hrk
zgCjHHeczW4g3Lzur#cuiqurnC;K;P2-`DX&ks0p(RLA$@I{4wp%zWRs6M!Nm+n?(A
z1b)~L0XQ;4-uLZ*sr;~)aM+JL4i^A&)%~fjAOx8U@%!-vVLxIZ2>*cwaC(8{gWvas
z{N%hrzyd(h!SCBaA+R6LG*E}kYIr&qHy2Zc9gsW}DY!Kd-j+!D0(>~&2#CC|A^%q#
htQ<iAvws!$u5P9-ZYPom_-36Sj(h&Rw2BPw{{!fHQ$YX#

literal 0
HcmV?d00001

diff --git a/Detectors/Upgrades/ITS3/simulation/src/ChipSimResponse.cxx b/Detectors/Upgrades/ITS3/simulation/src/ChipSimResponse.cxx
index 1c482983f0d0a..1c4f08438d4f2 100644
--- a/Detectors/Upgrades/ITS3/simulation/src/ChipSimResponse.cxx
+++ b/Detectors/Upgrades/ITS3/simulation/src/ChipSimResponse.cxx
@@ -49,14 +49,37 @@ void ChipSimResponse::computeCentreFromData()
     qVec.push_back(p.second);
   }
 
-  float intQ = 0.f, intZQ = 0.f;
+  struct BinInfo { float z0, z1, q0, q1, dq; };
+  std::vector<BinInfo> bins;
+  float totQ = 0.f;
   for (size_t i = 0; i + 1 < zVec.size(); ++i) {
     float z0 = zVec[i], z1 = zVec[i + 1];
     float q0 = qVec[i], q1 = qVec[i + 1];
     float dz = z1 - z0;
-    intQ += 0.5f * (q0 + q1) * dz;
-    intZQ += 0.5f * (z0 * q0 + z1 * q1) * dz;
+    float dq = 0.5f * (q0 + q1) * dz;
+    bins.push_back({z0, z1, q0, q1, dq});
+    totQ += dq;
   }
+  if (totQ <= 0.f) { mRespCentreDep = 0.f; return; }
 
-  mRespCentreDep = (intQ > 0.f) ? intZQ / intQ : 0.f;
+  float halfQ = 0.5f * totQ;
+  float cumQ = 0.f;
+  for (const auto& b : bins) {
+    if (cumQ + b.dq < halfQ) { cumQ += b.dq; continue; }
+    float qSlope = (b.q1 - b.q0) / (b.z1 - b.z0);
+    float dz = b.z1 - b.z0;
+    float A = qSlope * 0.5f;
+    float B = b.q0;
+    float C = cumQ - halfQ;
+    float disc = B * B - 4.f * A * C;
+    float x;
+    if (disc >= 0.f && std::abs(A) > 1.e-12f)
+      x = (-B + std::sqrt(disc)) / (2.f * A);
+    else
+      x = (halfQ - cumQ) / b.q0;
+    x = std::clamp(x, 0.f, dz);
+    mRespCentreDep = b.z0 + x;
+    return;
+  }
+  mRespCentreDep = mDptMax;
 }
diff --git a/Detectors/Upgrades/ITS3/simulation/src/Digitizer.cxx b/Detectors/Upgrades/ITS3/simulation/src/Digitizer.cxx
index 70425a1b4139f..b91e17890a6d8 100644
--- a/Detectors/Upgrades/ITS3/simulation/src/Digitizer.cxx
+++ b/Detectors/Upgrades/ITS3/simulation/src/Digitizer.cxx
@@ -90,8 +90,7 @@ void Digitizer::init()
     } else {
       LOGP(fatal, "ResponseFunction '{}' not implemented!", func);
     }
-    mSimRespIB->computeCentreFromData();
-    mSimRespIBShift = mSimRespIB->getDepthMax() - constants::silicon::thickness / 2.f - mSimRespIB->getRespCentreDep();
+    mSimRespIBShift = mSimRespIB->getDepthMax() - constants::silicon::thickness / 2.f;
     mSimRespOBShift = mSimRespOB->getDepthMax() - SegmentationOB::SensorLayerThickness / 2.f;
   }
 
diff --git a/core_dump_74974 b/core_dump_74974
new file mode 100644
index 0000000000000..e69de29bb2d1d

From 5d41f7639b2d653fd2ba37a9e55b7be4907befa8 Mon Sep 17 00:00:00 2001
From: ChunzhengLab <chunzheng.wang@icloud.com>
Date: Thu, 19 Jun 2025 22:33:42 +0200
Subject: [PATCH 5/5] should use the total eff but not only seed

---
 ChipResponse.pdf                              | Bin 23047 -> 0 bytes
 .../ITS3/macros/test/CheckChipResponseFile.C  |  25 +++---
 .../ITS3/simulation/src/ChipSimResponse.cxx   |  73 ++++++++++--------
 core_dump_74974                               |   0
 4 files changed, 53 insertions(+), 45 deletions(-)
 delete mode 100644 ChipResponse.pdf
 delete mode 100644 core_dump_74974

diff --git a/ChipResponse.pdf b/ChipResponse.pdf
deleted file mode 100644
index 0351c0af364bc2e7f77ff71398f4fe0524877b0f..0000000000000000000000000000000000000000
GIT binary patch
literal 0
HcmV?d00001

literal 23047
zcmd741zeTM_dgCuH_{<oLTS18a)FDqq;!LH!$mrz8zdw|x=TV)S^=e7L{UKjl?J6t
zy8JyCWp`cuth=w*e*gQ|*GuQlGc#xAJ!ekLbDncvRZ4~%auJ4mermdZ02cxRgB;Cl
zafO9(c{IJ8EI~Yyrf#Noj@Gz5s;1VKt{@oDKpmGy#of&g;Q;*R{rVgEP|48(Xs_jJ
zdHN7nR20|J!2&rR^bg}ox>%aJIl6$@)Kye8LE4rst_Vj55dTG9=tVF;J1)=3bD&Df
z)XfsaCUpr4hVp^=A%ajaL;wQj1VdTDU{;_VK<`&gG!brgz!;J?2q$$*S0_gYSIdh|
z7FIxK1xqgvM;8m%ud`tP!&LeHFjZ9-M{^BJH;_IL027Eu)6&z;5SK^B(ZLOg5Qs+x
z0un$10F<C0__q=-Nbp+;27-XU75P9AB(}(w{2&Mt4rEaP1UUf^D8fOI6KsK^AP9n-
z1hN6pAHsiHhJYZa)PNBnr}Tg!Ag2_8F`%a;fkB`rG?9zrk@j?x)o?>X3jI1b0JNi>
zql<=<sktRcSQyx8S2zgzb!Ny8s-|vkmM#t;=+`a(eJ_!7_?|~j0O8SaH*-5-86@LC
zaCwwX?E!v}0ZNxmT`iGq00a^Cmag1a0G2|S19NgPceFq_Sc7<U5DpR!u87l$KWqnb
zmH&Vm^hau!9qlas9xWt_|9}?n&uRUW)t&IiKOhPFBS|?+yX%&22y@fFx4IJs{0Fr7
zeoBj5#mo-j>~8rtY5wuxe-$Wy$t?Uor}_7{@Z>!GgGCDbZ-W19p~zG659q@GNLSL)
z-34Lk@;4#<qk#GwkP7~U(m&ipqy+j0YlDD)M)Du*;Xjoy--(!SU~(hr{j)9nhoT1Z
zovhLDvNv<I`-9**k;(rAwV#L+bw_(shrh=bCu0Afko!)Q=$SfMNdc1B%+$^G4{JEl
zAOB>~?_|agkep~2<ZysK{ia=bq%5x^%q`VrFB|@5D<D<kSE4-M6!3}uHT+!#o#>=*
z^8)dxf*~O2H!Z0OhJv8qw4^GS7X<yLB~`(IX8filRl$59-fvn`70eIf{ic~!!2%%O
zZ>m`p3<vRk)6J@2K@jgZ<*W)wY~F9$Srw3|yr*;ljm&#W7Z5_cr*r|0%zH`~(8#b;
zx`0}Tozex=I_#7#pw?lhbO9EEozeyN7j{Y)*cKQ-7Z+(Fr~;tEP6<PKL9kQ807kG=
z!hp)>J0bkXz5V@EpO|{5Tf(Dc>gIy*)JIMZus)FX1Qc#~`t(HLKzIOa3XoamfB^`Q
z1WXm+F~Y&^<RKJL@_=n<{jKHY%Z{G<$dLqrg&FebYPeZC0cD=6Cj}UfqNRhin+-@1
z!XtyQ0}Kux89Ts8k+K9V1504TUEKg{(VhTTE|Oygs@bDC0aXP9dl@*9uUY04CY*bV
zQWdPws+vc~c7LLj*uI4LTJEhj%R46dx1L!daPGJSK8}%~!viIqNrVLVPz71yGeU2O
zD}!}Fl;VM?D8A-ot4Ox{O64~O{9XS16dCYE0pv~qf02cM4Ih9Ze#3|V2YjHv!iOaa
zf*S}XlibN<lHF}4f5nMSI5naA!A<pE?#&^AcZ{J8DR?hpz%--sE99gqdCE@3pGmCx
zXQDG>@X$EvZD_UL^7t2MTc)Wh?#56~(=cN;py8wBn~-1n2|&L~%Krh-DL(%yJP`gL
z@PY9E#z4IP2%t)~ZQS4xI6NYQP@ddvegF6bAlfNqs+g+V*6kh7)IP1{?q0I$zlv`V
z&B?*BCS~W!8qn21Q$AC6@!oMvRD2xC8BtVBl>P$p`#(YGcT@6zAOy&ruPg)<{w)@Q
zfPn)G2A`q?6u+?%vh?ks@&6L1zdC4A-wxV2;FzIlr#-{a$Y>|;Bk#c=HbI%GBWU;$
zi4#HlO`0P`!l@uR6*FMK>HCk6;QtlJ{9Pp4Fw_#BVK5Ppxzq*Mqt&AT2WRHxPr*P6
z)xQP<DNd1+1*yCKCCU4@*y7&<!xC3L9WP(p*7Y|qW_1ZVf5I8R`v(8(jQm$%!2d1I
zKqKSlQ|*~WV@HSn&Kt8LKZOIi*uUnD|C56P+$w+}e<?B`yn?^NLn^M4unmka<q^>N
zT#;Wb$@k1z$<7v4_iLO%G%SzT4u&5N5Y`XxWX>>o22n-dOQ+&Izo@Nu(bi%o?2Km~
zn~&ksLBj#xDNI`Vpi5nEQsu3a<trPd%9+PSZdEK@%b1Uh(-YtuCq5kPI@jKlI^`+S
zzSr38_|A1CAu(5W7;9^j(gX!p7vlF5F8kdv`(Ir4e*){TB{2m2FS8kjx{QqX1!6QV
zi`-j%mJVgsn0(k=sAj3bEV=2?;u}RKv4V1AWcmDQcozb>aNkmWliT>>^eI2V4XMKa
z19tp1Zi4@w2K&S5417nBXCMh8vQWNAso0Cs7FNm^X3HO@oA{f4ixQy5`3X))?)gic
zASV}YNSTf-{eLOcA^fmk9bAJr2=d@cd*n<pv?81mj8G#{-A%;pz=3RMxeOSDu7|$M
ziyS)`*Zy1-DiqNww>17qhP6R+SxfZRnpVu27nEHO`6xxOx)M6BpxM70Ub#T^=Dn3(
zqs(|%)<cy#8uXFqvCca56q8j%S;CZ%;miXz#{HErTdFS#qqk<hL^O1r9Uj;8iROwD
zYMp~)e!M}VA`<Hx<y*vjgieJLD<c#B6U2XMp8t1<L!rOaqr$2}4((w4Q7_zV&@Cy2
z{4&*C$ya%j@lT)L$|t;avC&^Do+sqWl;19!7~AI<cIT<g!-FTUDaTni`CA!Y9EzFx
ze_C~&dKSvJE+Lrn+R8jFFk7a&^@C6(Q8<C}(|Z~QDL%!R9^rS|A2@Ncg}q<1<hj0L
z`n-Sb@k=vbkERgEiEG&Clz1#=$yr?a_ubgYaNRq#?Vphx4LvSTv2wRcwsl{*?MX%U
ztWqQ9%VYyr$o*11gDQ!tFiWYqh?hZ0Q7`H9mdW?!rMZ&%XSRE^W(ehEJVob2vH_#V
z!p}1HC;RnFZ{a_)0Dr01{;mahY71lJF*1StPu<Ir6ynmK;PFe};N<%ML}2`%$qV?Y
zYXB6!c?J;3FQo-6S_JSh@xyj9r_M2Q=WXJjx=@^RQ(H>=yop0Qepl#S;1c*ZruWA}
zC_Y4ChEbyDm__Rfbd?SBLPv7_&wA-Ur?_!@f7q{vzBF!hF_#Mad_tk4h;GWq!excY
zS@CsVrA<6AO1L9a<&W5%T%Y{4v;{7JzFN=_<W<nWWJ2?ULHxh@Rs5&e0fleaK_S1?
z<-+l&lD3UiJkq?@L)=UpA3x&mH%a@cyG-v*{L%QAq|bT1L`7|3tPSHrS$oSdjkVbZ
zd)*J$f_R471ROKg<*PD^{Wq+~BE8*R7$ppS&iPSmOk+CbC0i(5!FwD0@Hy`YyDs~v
zm2{X=*1-?~$05bNv$5_8f$LWy1E|lUGMku8$zF*>H;=bQ6r~l%7pE0xu;?oBJ{OJ=
zelDyKX-Ur<V+z}zIcH}dt=DH3scSo~U=YOe?hW47L&4CE{03IMc!F1vu}F(6iq$yg
zCwufuC+F+c65!+fpRxB(U0(?FU$Y$AENyd9r0PNM`gHp*w+^GprDkwYr3Lp`T$Gd)
zVJ61oAkfZ=zK<8-B3@HTt(r$I$0t=vB3Fyi0|VjB5LTcOpRv1qDkW^RKtJ)0PA&xh
zlGFb)srgG4_jjbGy+uCqt`1E*6Szdk{>{iDm?yma6PEgqe4}3*qx}Dh3i$>GaF3Au
zr+WlK)Su$<OQ!~M@~ZemhXAF2OAP*<n*<2#myjsN3K6!6<4gJEm?;yFaPpg$2TBF8
zv8mn@M^(R-JGsAG9n1Fik=5An!G+c{9VXb_2Os8xQ{=J8bBKMlx=4o7oC~^2qTlRm
zVvg<5RJ{%Bh1Hl9urN7Xyd%<A!gy1VcIYfZWT5{p{GF!Xr;EWKbq@AzpU2NUx+#*A
zH$%5x!<g+c6HxC6ABgk^td{LI{*0fn+^_wg(@V~O?~eR$;S2rcRj{z45Mc;lV){jS
zpv|x-r2_#1AWY1T#OpR%w@fDYg|RfBZdZ>vl7<sW<vDpTA19M=-c6pNB%-_X$~2vv
zCldFC%e0ERO_g7N#0-sN4|`hzZ{a<C`-XSVni*N%?d?UFDCgI*#8nsF(<v{+oS!8#
zL2>-_Qz7=BIY)ovRs$mNE9Xcm26?59FAbPG42<lFV-ilpSXj;j5_f|y8jZ7Xlb<!(
zr!{goy9C0N)NRO#xnMx!PD4-JKvU4&0M<67%p}WaqOE~SI}G+FnweJ_q<(3Ar=ye@
z#ikX1^stqjP;9PWH71_+rsJI;<V`a56;`XC@Zo>v?fe}Ge{J=|h~S2*;)fn&n#JA{
z#$F|CqcAgNyspRg1n-@vUqV9;lW&4)YQ3COAQir8<P(qdb%AT(&!H0#S_UpobLM`q
ze<G}7NLstVrTDR<in+`9fT0{t^1ToDTpnzc-(<t-D0%x_vzDgOIX~VA!0_=7+l8NC
z`2SlE1oCVBA$;ozE*gYb`d}|gUtj6%$kGfYI$mW3u@HHp)XT1M5BmE~pHmYiS!R?l
zh=@bf;N3_E9e){LRBeqfTI2-=-o@x^ZgHtX?1S5dGNv~zUtW)LIjf5qDT5e<mdGcz
z5%1VU-M?b+Mm@=q8piukD@)uIM?YB46h|Uh+GO5e1uYBsS0lK~=AwOR4STvmoLdi(
z^Agr#@#YuF7216qFbcN)rK6vK{%beqKf4%#0$$HAUN=gq2`ZyY;P)Mn$40)LzRqAQ
z&LobG9)W*L{LEPv-+<9s+gsOX@3v2K#PmKPYc#Mp5?i-WdI9Ykj~>4|=8$V2N8jo-
zZ6M5rC^-X{ImiCV?lpzu7t%#v{rdsa7(&ncuX!bw=Q@NO5+{x4&MIeXJ{C5fDeN#B
zfK%a^XG>>zs?dE76}{U(m_h2Zkaaw!u_|MCmuNq{Sk%+e=7RCImc>YEpA;!@IYSs$
zl=U{7qp>^G8_wuTJ=$f^CD~<=blhu@)P1jLW%G&Jr}|5KlXd}5EhAPTk$0I7+vFH7
zO_}EnR0_NV(UqXI4Ve01oNqHdVAJz_MsZU-#%a)xLZrbmP8#m|){jDozSplF7mN~9
z9ZvY;E}y(s{!26U|LQX4*W!uqm;UT;Z<Kpxq5P^S`E(au(d`NI2?+4X#HrfQULo{U
zJ57{6_Cy7#qX5$7#%rR$pFr@tzjGQii42ECM$oxf0zq#lZ}5R|MN12WDey8Jq<`X{
z@(J>R;4ml<PAKW(==8N7X{PZ3F#<rK0vw8br0#eUu_p)w)cpR&9cgv|A()atKpil=
zD~RtTbl?wf_D?Fl?g{RxgUW}rpn+arxgH2;{KIQ@LuB7?ufRPsmkkxF&)2IRb6O)v
zThm_a>s3_5Fh-c5*0pm5hWULkm+XTKi?8iGO@wOXXu7NHr1C?~2v|X~ygsXG$;K(N
zi#;7B-0vUFcNEBZy}Z6Olk@n6*d;%n`EBElnJ>pHJhKeVU$&M@yLmnw!5x~9ci!uW
z?ZqE|_P9fQ<(T0(axmac%$>~D7`l#vv4Kg)<}aQHvpFA*uTv`RWrq(pR_Q8EW}DSv
z<g9Kv74c_BXImVGWY6%z0<N#4eR*>1k@tzaetJcc|Ahcc(tyaDL5KOePh?l+1=+OP
zMaX0uq{L#a*{1^&g+7cH;;MH^u<NXH%Ot83%7^*#poWD8LD|O=o~<qQUcD359DW?%
zm*}E;@7zwo#r`t7TN|?VL&5zHbcq}D*}S+oXDrR?x+D6_GNn0$_B+=`VnWmlM~T9b
z4QxNaldW|K<uGK8i34mba}{$}b=F3bs^wng7<O3;5X#S9ZWa-J;;ee_b<P>!e{l}_
zL4ucgP}MG>{pEFJLyUcYdq^Q&V*W*Kypx7f`|ImS63@rYFEg43MsJVrEd_YS9gpfW
z&yJ7#`?x~oHS|km-3K#QX4zt8KQRQ<`E60gdrylzqH$noYC4+Oom_5j{@}lF8eNg|
zXu+&0d=sZ&<-Oh(<LI^H#H(NIpJS>tYwd-;S+jp2^Y&tn_wcU$^RtV}hF?gn25SBF
z_FH!cUdv5*Nj2~vd}8>rEav|?PHeaTli0D-yXMr;(DqSb)*)K+Jnau9u7z%KjQNQ!
z-tAV_Bb9CwS1je(u?lT`onl1SkyzJ<V0!JFX~+xi8pMue7D~v&Z*q^XxO;iea;Zbc
zZpN{l2QYo;QzU_NC$w6al7OX>BJu1Cv(iUCJ<o)2u?i0*njP9aG+!{yz8l5xm7Se%
z-6KZqw*ob0oNs;XcVmI=n`0=P5aO_o=Sb2cZ7Og{{$xwX%Du`)S3_PxDb-w~HOZx(
z&6Im)Tutz%W{HQ+`njZ#esVT{At=yTwAsTq`$7NW>fq?LWnQ`H6$OF=m?*0)*PuHa
zu9IU7(<Vhg$>@ib!a{-NF0q$oJ6lE@8OD>@UsRtY(P()d2n{+4$!i2oMV>ROm$kj&
z*jJi^cvA&k^~s#=@huQD)<1gc=)Zno9MEtyW8A#=ae7+pc+WV%=U}8euKBT_`vLXQ
zB=gaiV&edjb59QB*5I_4WNyHYuiTB1=6O5pCB;cj`{W~if@J0GeiFfSFO=|iqY~}Z
zE;3`Nbm@ZE*m|fDSYS4(Xa18mG*H10M{%>xEVC%eg;;k!;?iA#-;yNKU61D*MTsrI
z>f7zPZI2RyIi9!{&wY@ZDMduJ&zv&h$?--pB%Fl7Yj+LS1y_kPnDbLu$&9xqE+&d=
zQeRSLmbuWHe8!cf8)RNDlkM9s5jWKPfMuVYR6m|9#Q&`eOnWLGraoXhSWvnsf7=A^
zU%PTi%DbVjy12Qvh3Nh~Mt@M5JEn+XKHIaT0}>%w7rg%G_r(~f3)$aoS!0X9U-p?1
zX5^{i&*3?-RHED{uo=aek<&xXY%wn6&WzVY3BTl6iX(Vi9VMK;iOh~XT@@u<)JlGs
zgC!?5T1teaqe79XYWwnYflY>Ym-3u#`tt8j$&AN)Bj5x1>#1}uKA{K`W~F0++H%Li
zE5ct)UhX9^`_B#KAF^d4s?9Sed+_aMrcxeHwHo{F^TqDw@yDCvTEAPOvW7#?(AMu4
zLLSz4mb-KiJ!~Meg%4PD!>8m3kUefB=`kZaI1jd%JN@2a*B$1aq~2QYJV>;~EMNa(
zbF=g-jfIr19(<9Hf`SX--kb!a*53PCv2b-Hi7nwism)Tlt#<QZq4<GPaK}ws$C-sh
znFDOYFM@&}!)|0TXCj<9>74Ii^1Ey)xXQQi7vgVSmZZhnVfFUbimf9#O49Rt%^(tO
z+W%CI2j;syXIfmWGaJFcR}o)ooc_X~XHEwoi#cy_l$PzE#2h==B2(s%z!Vjhk?m3w
zgOT64jhyFK0-qIvB6e;XPnIyGJDHW8O!Ly&6C^zIFP^yrtV`({;}c$yBMw=*yFoDq
zN9u#>v6i{&vAgRqpMA3ani0s?b?t7^(Vl~=Qp+AWZIeAi66U2#JL#}VN=vBO1rdWW
zWnNjkL^UP-P{75#ch99Ws29Ev*0xU3mX_^sI!*TE&E5ALl8vf8#(ts$1Enq<F>~dE
z3kmsJI!?LJSntU~;~`-7_vqFOPT+y8Lv9DrzZGK*uhvs9m^m*Lz$ASP5lqCQC+kl?
zR?g^JQy5JffM21{P8e7bEUb)5h;=KHo`wi+Hh1B5*_swB75J<K71y%NBIxAh)$bhk
zx?X&{5UCy0J+o#;qf}w#ifr-mCGUkRH?n0C7;@(Dym8U|Tv@$)#4ae!D=hb)-@dLl
z^F)PENb+ImqM%^V9nZ+LDtVphG>VvT=bo}@@~aOk1~g8nRClZh3dX0ihhb&Et}RI9
z<)$H;W-+=~V0542y0>k`r<QnD>A-YxQP%MN7ZMtZxoo$1`(nuU!Xa}yv58+Y_LVa(
zdprwb#RyC<0@?&n+LpAgU-CNP{LGO2&V{0|EdH34`ta;z0<ok?LKRm^$AgVBHNjW+
z)>7$|V`EJy1YXZP(IFH$-`m&7tfV}9gprT1k+XII*3zzFOD5r0tYunPXrv^O`*BNU
zg5~%PoAq+WC@b<g>)LItw0$9ano%CCoQLxZt)x^R#yv-4G4-N)UbyLaqJ*Sm(@+f~
z)aD#~Hs|eQrUV7EQ$ifb9;-GDDBtN33nV?z>#M07xcBPFW%Fr!(EU*1APt0XVNSSC
zA|h9CKyn(DPGxj3BrHoiW+iP}WhJR<Y6G3H+`uZI_Poe^Fg*>?#2SqZI-bx2no-v{
z@pD}9D#}HNng02D?5lSvH(z?9g!>cjR~COA!2BiIt7_*C(yh+xb$IWDX=R1e!m?!g
z;|pHuwtSw!SAnUID%_;3d6$+siO?IK?{tlmd{qhUAUfi&PsV$9NF#e#0x)mm>%_Bd
zQ`&-4;lXX1{dc^yZg{J0Wf^^Np<Cd6j1(Ug;AAI>4D;y-*6NCZo7S<?hLZhVBr0Uo
zH!^u_KHJf)TL3$|{$X0}gL69FdLn!fL0RL)AL|Ll6KeGQn7S!9Hv(UX#Eh8XnVY{V
zIQ5ikX>+RHamjSUgf8+uduO33e-`5Avlw1958jG+1v&;Fsi??2;80433uk1_&FcuB
zh^K`V^Vg^K<l49~zn^-+3*4!+G_prEmDl|9i;ZU^m;oV?VH_(*SsTJ=#HN>THy6y$
z*`+2R8BqAP^G4IT?WXmutoOU=(8uc-bR)azbR+vysNJ>b`o@n1i)f8y;_P`gt`B}R
zeSfpPSe?nv$GKB-MVkHbvNY7)PUmSEB>~5qoqO1$+xZQpY6-#iJX3UDcJ)sm>F{tk
z6FzGENIU<Ca)~j<p2wu;5v5a2lY@)Jdptq&_jvHCWq{r`+6#t&_s=xOC=S$RZH1jw
zH(M26Ob0ZUBwya6>s27>omF>S=aHhBO_?mjb`NY3Tf`$7TU?X%VZZFx(dbX{>26xh
zraBAyMQktj%k^Ze73ETA{hpjPj#|5|c(Gi1?)_}a*~IPeN|K!->@BS?IGlrChMs%Q
zUhG?*N1Wc2Y!*Iy&g?<?!wa+BTdiuun;X;no~p9`d(M7$c!u`7A7siJM82%qj`S#s
z?-)v9NJ+WtskW4a9m6NP!s)fCk*E+aV|;Pvd~M6fVwl{-Y;2-Yl&Yw%8Bvt*kiJo*
zs!>+AoGzPQ0|jER@ZdIe=lD=8noa~cMcxd^AnGNR?1M0w%D12P@Nj58mvpYdrOxea
zR6EHp+j!2h_quX<Fsjczd)UOKw~u)kCO6B!RyC{6CU<V%@>z_m@k4L!%f@95cJ;5E
z*|w@?U8YE1KQj#iUfLZ+^Wsn%KQbWPivv2y?AxaYj7T#*CwbF+SQq~>V#-L&VK!}V
ztJ~h7a?88VkBGY&FK-qf=ZiI;R$+ZKpO%4?Kl=T<Q|AwitV)4jbF|FpDKPI%b{yM#
zb^JRm6R*`2sb%W!@@Z`n-^Jv*w~_YvdQ1Mjrn0o;d=A?@=44-<dsRkkqgzkUt5?#`
zwOZITaWs0S`A>hK>)pJ}zT9kev&mDp-6H{_KWMBT^MreG&+4W>XLD2i(^T#SL}lV=
zatt+4W0|7YK%>*AJ^l9j(j&VXh`xb~8`oa=Q$_ul`{4=P=ZcHCqwwEt45USmxX7Mg
z+_O-SsYBen96|HeS<Rw|`>Nl#-TW~2BRJtz8O@}$M>b&FSF4XKG_#7>E<E+8j*_;}
z%P%vQe8)ZSB&^gaDr`tFGC&a%&^a`kOgJ+@vDH}DkQ=3_!IIVZUTGd5tfUrz(+zbj
z>8p{JbhzXW$`C4gAnMJ{8{4a>57X@~SbF64HtkVKGu^XITd2v+8y}}PV=oyi>61k~
zd<+;2X^%Mg2CsX|X3h5{7e_ouo;{~!_0D!sb@4csUpAQAwuE~8l|>(gwVb^p34gt^
z#-P575ru)KReJv3{fiRqY>d<c8!+jt^XKheazMGgj1)(F9zOpxK<1uIm2#dQxjfwG
zh}Aa~IEwm%`fac0l&s$6KNTo!)X~Rhy!*PE(!!H|(n>QAnj{)e8;>O^aJDD<HFl)v
z+?hcmPVvUhb=-b@S<TSuCE;idq1o2w7BUY)Lb;=37wc%t(%^?@yYWe)^wIiJ#y7L1
z80dQNHO{We(R823PqbOTcUwlBsr4y2&jV;U+Bvb5w=(RZTrJM9*WP#Wq{WgF3aJGV
zT+-DX#fZK<l^NWg=qb#|3byt<i-(i~s<ef0^Dcu63H6Efyi<*Wx2QUC&!)Y4uQs&M
zi5r%I(PhgdS?)(?%0-S@bxzz$6BL~qi^dbx<c*5Pr{*mjv2g>tzfv~5*jZYl$kaY6
zItq;dlYMhhDTot&?z6e44fMfIilcX$_!NO%0yXV6jm&G}l<pfM4^^WIK@~9SeXDco
zT^V2E{g>H!`6nNi<TR<x8Q<KGd`2Rvj;Ag;9&WAwV0A+|t+Ki-?^;6&e64#>s$Ll_
z_-cF}BVTN7?4Zm%cTrzOp#+b67Fe}+P$q?KSgK?Zzwa_nt~2m0!5E+3b-%PIC5SxO
z<j@XU$a<$Q&#G6;^aAgr*L*LECLTK8(eZAGcCM9onlSLeW5kOQxTa)f7jw<`h5mk{
z{VS{<Ja4<+mnEY9+UcJ-Ft^Oh_0Rb5zVW`<h=2TC(>zmc*zT%)ZjmQvdh?3(+uUSN
zflCJs`c(#Piq7^eYCN9_s-+>OF7~ySS0hI9uTPrR((@)EZqnWux4(L9_Nb(sj@#%X
z-0GQ5RkrdY*k_1-ydJ{vHWzR%gVGVDPq?3@6>-}ZJ^KVD;U)K~eY0G|H3HQzu-8rQ
zcmc5LbFAgTIxeN^yL`}C&$JgxT}9es8V@Yv5L~Spk4kFqSY8rHL)`pWU(Nm;znFc2
z%XySjEe#PB+?KhMD&yLVh`OYpt4>aymTYq<U#!@9TaC9fBf<ZvrXr@DnvjS!hu`_g
zHd&!)spERH?b}2MLt}Y%yG4bzpsm8T>Z7WrJDK%Og9hPm^{kRso3<q0Y%miJ?PnGv
zVu47O`>JhO_w-!}yB8nb5>!Q-zr>_F9iR;IEccLEuWw)C?3jWH4-uc08f@ghr;<TK
zN$LH7L{QF6fm8$APdORnZr0B?=mx9$%-?}eb<Sgmj&X@u9z0IdxIQyWgK0|lV93%d
z@>$1pmFo!ezNzeG+WsrhG9EAQ%siGL%x3wXx~2}g@iXjMn}Q8Igb*KL=8tDsP5h0$
z88T(({T8u>-1uj(ozvd*yltGtbtT1ka>Ut6CK<B7`@pSe;#CcdZuPDg9oCA`hI!I)
zeUt2>0;xJ{BeA#-yP7P@%Q5+zUAXVA`AaWz*C6aF8pmoWhmDM#(T0%^(B4n0qMG>N
zAvS33R1VSX$`4S1`iygd%azQ5%iJ>y9>ELy@P{qi99x#5%MWU7xGEakJn>tABEqgk
zWF5byaRn#PhtvK#O5KnMs*7VGWk1UOMdo$1(QEOK1;lr>^<WfpM%=!>ogFJcAFc@z
z>UGoxSyM4<Ys}~G5wB3zj3$zl1h|)gXKNy9S5caVJhynMgdLfs{o8itZax_<0xa4#
zv<cp&P-k>b)3t^yewLFGU^BHUvvnyKT?tvUeDOvP95|+2iYC6>b0{}+qcaA!hta|6
zx{kg+ONc;VzIu;%;|aRo?`#Tg+17psqoB&^+avAj<`pKw!+i}rM4>WndeUoMTJk)L
zp8$OWVFwBFP5bmFlc2~WUwS%7wYkF;RUDPl+ZL3lVxEdKJKd@I3BR~Dfr&_5)Vo0U
zbLfiw=AbxRj*L8MZ}+@o5r2tzBPllx@MH3DlGSrnDg!fk55|;X(sX%L>{Vca>dDi-
z6v2uGaJwnsn_;%?v20}6G-VgdI_fr(ijZAMjk_RDld}M;7`PW`BO1fQnjY7rGJyt`
zSX{o-Lz`ai0wX)r97kw~3y%wHAY9OhLU7*-#Y_6<?|I^+pAl6X?_nph5hX~>NsusX
zgtsqb&F8J5#5N1<6FhH`xhFie2cyaH6s3&IW@gh{mAct=1|H9aRZ%~1FDODZW<RH8
z(*2EXL?f>Vj2@89#MTwSJc3e1pI^LC6ju#eEoIzRuCWmgB&9aM<oD3yBOWNe+Sjq^
zwI5G4P!NG!9npl-%B_}<%B7fQ`!IcPG}cX9Kj~#)77%B@Oe~kz`I!8#XGn&M1xGrz
z)@fZDQ}e%$*k0+5nE@E^5sL{QR&k=8XVQYV0p}sLL^NOxVPsD;jL;HS&9O>VnHWU9
zNtN43nl@8S9g=VfYe6V$LWpxh`SUQV6I;672*}6U8<?TwG%1g1I$-e=pll%?A7)1L
zx4_Uwga#xhQ8Ql#SOYdQM_jJx^QWvO{XQ3yc*(|<qYllX@TDszPQ_ZN5`O8sVnLD+
zE*=g%n^i@k7%Zo`FP%sagdlWl9<7}eUXZ9}`jyf&=0-o)QjINFkLcAN6blB1imsQ^
z(koeWd8ospqGAvhB!-nH2Ntd96@Fx1;b^)|q#`hs_bCQMXOJ!lLC5xdTVKca`r0$C
z$fyrOukzPX2-~%Ahd&_Rh}-2{L=g#UKO+N~8B05HbZ(Kfwkel(1urLrU1-@>%b6f+
zx#h#^ZG-YiW&2%U8{pRhkI1=OwlShC=5?SmW0CYUfy)IipKWP)bOMht%?>Ok?aU6$
zmM7qn$XTKci{iSJ%x4_S-`$K_qP^zYj}fJVc#Xw2ewlDDa;acNZ?_Ji1+bi_>~<sb
z!C1WCYEP5ttODpjkH@xU<Jw2A_vT(?Obs9CuNg$4y*DSDa$fIxEOH^UUPQaPyh6lE
z3$FT1R$^NH_K<gc*N1b{X4$F9qm@R+RZFq!Z0~Dh>^!0uow`UDS9B-L4kMnJ&osnN
z(GSxL?Hk3~$Guv5DKPbMo65b8(58L4^5qVOy-)Pb(6idcLPp#OUBBgftM2DqefpdY
zU$O<>UWQz(5~?_Sw<S`H``-KzqnAyc`n`Fd9ATwj?6nMBhd4+$WqEE@A-DqLdN-WN
zqe7soG^-#==26$K+?D~YU5qo|3pmTw+OE3vhcT`Y%WL<YKc>F_K*QNUXh>%l1Np$%
zjiPV<!LIk#o+#RZ%j1q!l(JO<oWLzD!n5U47w(sLUdK5P)n%X@Gc-@ee@p4z<aIu;
zq^b5i*&yXzM7_fOl%V}P!K?C85`LS~dls{tBp6$22RFO74%kZKXkL&Q1@GIWPWV0L
z#6jJ1PBY>fTII?WRVH<5^tVRcTK8HHe^Qg5D@w!kT0xG1vc=VJFeyTOpN18Byzzwn
zPMi*TqYuvf^lXbi+w1qYMCr?f@OPBt7>L-p4aK6w_v=9QM$K`ceb349W|__P)VWd1
zRaL2+#fj=`xmlq7eZn@{Iev-F${t!`LDbF4&?gR#xO3A5z8`3NBv3cEloF?D$^=mS
z8f+E@w%T!X;}pnUoA{XmMaQU1hwP0}{d`}j<E1p`q3zc}Ki6Sua^%JdBtP3CE)Epk
zkX>K5bX6B;$ie7(XrRM{J9iVm#20Jx$Qjk|uFmRZ@iarC4NxvuJ>|Hrc$#3P3YDsr
zZMUC%cwkjs;fJV8y=(&yJ*&d)sZ=)DO~@-5iMiy3pZ2`H&VtXzjAf!Aj=-X}5r0~8
zlb9AgG4m=N0V+P-9p{8})bo^f#VZPD<UiHP^Gb+cG1ahDXe5upE~Jb{h>xt+!oj&S
zVI@u`jaoeO^wNiCRWZk^=Vjr*P_3g7Gz}a!cDZ?sJl=eZSR&oCFZ2p7C@}SwkwNm8
zn2w(F<<BI&c4CQOlboq;A%59?J-w7zvr3UR`TV_^t(3uO3r^;br9R2iOrBxZi4hVb
z<?Sw|9b!!M6s^zh7N9#N>wDoSHOtOwlfKX~<LN|&kC)J4!HhQ>OHtAtDN@d)#xZ!c
zZd`H1<YI_H96m!!R<YpcGIsUjK3m}nB}7T)VK{dNRaZ3rYK8BI)P21&(>*(MR4(Id
zjR{oA{T38l3<4uYn-027g4-pE=2WQ7@uqu<R;?S(O65cmvl#`+5ira%TM+p?Uydr!
zp65uDm(<4dNo7kKk62L-_S?V+EXE7hB2Z>Brtd^J;V*jIBE&@#BtQ(~)=8A0aXQm*
zY@O>Jk5Kj$El@YE%!vf~f?<m$&!5Op2l;APp>FU^y@djGMT>0Xl5>@@;6Psu%phN3
z@3-rf=sj4~LA_3}aRp)W9sx_@p5kr}MG!+DGgJE&qZK(8j147f^o{i<FFg7{U+|gv
zOeu}0;+<@_9D6~x;yx#TP!ZGwF}yH@W^|6}vBo+^fnwt%mu`TY1V7!Cn@(DU_TIvZ
z2HHlKpKls|fEG<;r4T=PO!XdR+ns8hOoEDG=OI`*o1BV4vHF6f>3mBH4$vAUx?G?^
z6Qv7_v?j{Q**gbdQvqq>uI`vzFY)!(jxr|SxrNkWj80&>okIeuK3h5YyC4=->eYQV
zh#IgoHD7~cVpmKu26WO3lIe}-*rR|MMDfu1mN{ThF@!R5v`b^i23!kTSE3GGUoUdN
z_HWw-`;(B!%f_MRlTwW*0`p>3J{-cC4(aJOW9+rLBqI+1iJ9k%zeKYvam$w&dX+fI
z+w`pz0FKm#d-SjkN$^ZY>s$QhlJF#FH=XbeC6dsMt(cGgQ<azZK=Rdf+}Ae;LxI+T
z8|xa3w}|qE2zJzsC3D1^*n`EJjuq8qfIT6m(95MNEdGXVz&t!QlpMr>o_Ha_$M7Bc
z$v%V`wQHjGVC9uOU3(2f+bhCH+w;bw!y>`EMHHf&G)eW5-`Q6K-PyO#owO4xuN>#4
zD!A!M?rsr&{_e2q6&zJxVps;g@g`px0}Wdfdv)%vCwP8mlRDP%#`_EHKKj#BOs(r_
zeZkY_EZX)F)79dmCcX`qg$X}Vfkl&aRec+JUM((Mh*Cj2oJEgw&d&<kl|qO3#_Ozo
z2%fC&6W72{9-1bA08DJMJIB1a3S;z5RznGOv^%4SozuHS+*FW=(PmwbExwCM(LwCX
zu@O7rBfI9;PwYE=PbA&FD{Shroxti=&Gh*Ql^4AohP|yphP?};k4oGI&lgYm1b%*W
z#j9<tkkI3Xi2TB<RC_RXFjkOdf_lYa@;RqgXD8xj*Q<WcIQeUt-4{sg)YS!I&zTZt
zB%U9{%(I_$w(3Yan|JB09zT11;IJ$<<(`CJS8o5qmy)-rA3I}n>56x*RoeBRC1~%K
zeI*ca1FqbQk&;-aeV4!_>&6AT8gy(3o$M$T>n00T90S3*RwL&ZA3gMleT6*^qCUlK
zvk^ItPG5OC<aenQk(ETnTyv3auk0w`OCVbvX}O1ZgIipHWP_VyUU(bsI9}7mwZ<^*
zuvKhtm9iq#P5N=^GZT$rH8-tLo(z@eqppp~^XLqKs+13vH|d4Y*2XTUlXtARiv^DA
z&vo)lEbO<`neFT@GG2GglQ!>Mac1W;+o1sL>80$h(5_sQ+vU-GCS%G(7~(7I`G#WW
zrSO+(IP6*)Mf=KI8udoeKIq%L(hc-U9JWA5x&yuffj+{B2X!K%a585fJP|VOVfvJh
z%yy1>TBvG1P+sN1?CHYWP&FUg2e^Unrk~!QG(`F;c(y7<#^TF+mH2mZ8<RA>T9&pj
zGKINUY7ihR6;|LU##m7)X?5*DK5D9n(ds1k_0BoKJl6&LoK-}zic6TJx*+3NENyPK
zz01a2_5_9p+*spTsRUEBEUC@;uJJ$2-@L8|yI@;hPTd~n!}<=!KWIc<fiJ#H{=yga
zVeRYu%G~WJeesj2SgBkFw=AvG5EK-U6?~;gm>kOw@u3yrqv*>H^(L(uy7`3AI<*0?
z@5GlO9OcZWuchW00g84im*oSO*QK@3Q*oX6JKOsr+V=NWYuLFeL~smHAARcbPDR~P
zMqXPd#x>icl#Iz_U|!x~5zT(tc+0yp8}+mDrYg!@wxeH`IDy@rr#vr3aN1BGr%)do
zvVvEohzJ5#)(ssPgOIn~7>`+7Ma6gEdPE2`pxaB7IV0l*A`{bMtuuK(>6(gj*$kvv
zQybUSXBb){&qO0AJRXDK8WQfeAEP+Mv8H(Kj78FIpu8_o%nK>e3@BDI4_b~v#{mT3
z1kF(7)A#dvD7lTV1+ZGDly9~Lze|AOtNJuRyDHc^XO?ShP;ARxuAEE#AeUqe_PM=;
zYExcZWoov8oW`KhbTAe8#0D3+!p>?Hw{z|+PdcIEm@?sd_Z@Sd=$wzhs}qt+R>UGH
z)t>0eOIoJ}-cfS;7PvQSHs*CbGL*97`gK<*l~)pev`0~RTj=8Y11oo<$F&I;;+#a3
zS^!7JCBMk2fmCpYA)$4m9&JJB0WS8?0M{3s%s~$+GGwPNpi_ci;HNB`0r9%*vI~lJ
z^yM{JIt0;3=MR<cBI7d?R{*F<?<(7MbeTyez4ocB`LZ<>Oh*n)aK%7=pst7^truFj
zY@LZ*!->Z<`5CVT>18o(xnW>BOitrk8$Eg>1iS@;O;CuAwU~YZ-Wh3~iT8f|I;3{(
z6-x`A?DcNo+z`}-#c}x<^{Sn)q`HcB;u=L$Q@cp{@cSM{BDl!LfM<m25%2;Bqj<+3
zYy9Ax@2Xs*$DJ`D_0=bFgF^N}d;Sf~+2$jhaAzz#GL<&M)TNA8NHPa+ZHvL(_mE`M
zOdoh{i~JIp_x+h?bRj(5x-@KfL@h*Rt@OIaT3uDPU9)#8?1I#;UnEbg!syytoFpku
zInZZRj)PV>ZIV)0;@t;RaZg6hC&j&9*VTT0jq24Tg)ZlyI<7X)y5C#D2@0xmUIxNv
z_eNXA2(p)-Q(=!3XynrKg$HZBX;{-Zo1OAp)?yYj!=z#F{mh+(N*ryk#J)Nl?K*ze
z=?l6me1tgstyT~74ck6>X4%(=#x-{aJu!Kth3m0)MSuT`E0*!c)NQmgLAoL`(MAy{
zy5$Y1(U&XeDOC(2TF&`~O@H7N?rI2G<>ridTcCS*-XKn+$TOx>qmF)-&~ajk?13zI
ztcps*S+$z6_b=%r7uDvDl9?qIo$E*WZuOK+a2`&Uq564aFkmb-8Q#R!i(9<)+1XDa
zGk2@*P2DqDwL1D7i==XAMDSZq19W$Sa`hUc$>7aJl~NlKqZ}`cHwDKA^~1~Z7UEkr
z-W(4M6xcLZBSgZR&Z^y~kIUNP(_W>zhN8n1^=4XrdeQ2Eu1GTeB#p}jT@ji^ytp78
zCbYMu20;fj_jTmUX`b9BzHc<y+v6F-Rr74;we#q^r412w`R21E_vwk2w`kDpQdXUV
zjQt|tOjpZJZ^dD!tTytkH5fPT9088LO2DmZdQi#zt>G4LSsD{m_z@$$g630j^x>5n
zBOCGp`Shi@I!$f9{kq1PH%@d7t>nqQ_q=LCw6!f4GI8Qb(r)q+3DY&MeT>s5ly!v7
z#NvA6aB|8Pk3PAP6ccApMUt<nz0N2z-A4Q1{#5ECQ{AVoH+VcRU2-J3;&v}dt0B}l
z<=u7&E6)A<Q`2#m4e_4G8QqSWooXPj(DvJ*^~`;APsae0t4v-v|6<`Fj6IIVkF#uH
z_PHkty6YiHW!*hlF<pZ<ljgemor~A++CS8KkSXu+CC9KS<p2|BtWhA=1hvKp`)F}`
z)JM`PWN|A5q}?z!_cCIr^k^5qV^R8vQ1>NKzsSlFd`}T-T$FMh8EVFfEh($^MXw<b
z8I*F?5bUWf?}RgT*%5Y5jUk+#DWyAZn`t_$Nehokccg$J3ttmi0gq$!Zw!5Q5o_;Q
zjO=3|sOXk_TB4jDA##qB)vP}boq*VfshVk<P}B^I*vB=Iyn{Z+_?3o?<WgGCjbR~M
zW3B4aosVnxBSco#rfXl8?kHMGEID6o#%dX(d+VP2Nn{olWs6Z^4aGB2Q65KW8OtQT
zV8TNv9E`fHXob4u{E3GAJnnFYMr1FTP736UHq4ubQPIt8<ccyJazzf+r}2;h)d#S!
zmz)=dc^1@L=#`C`685%dk7qF|0qe9<XJ3tl8sy_T#j9z|5Wl~3wr4C+qQ^*9fM70F
zeA}c~Q#WTCtwXKig9o?>`|`GpL4@r1w0aMH-$UGyY02dGbq_4Ww~O#EZ&Mu&)W+dg
z30Nb@Ct(hK(`X~S@fcM@-JsB==Oc#Fbysg-R9Qo@N2aI3;S+XfBO!Gy7M`DN>CW+i
zfcXV2t+y9pR0%NSCQM#@B)+YP$>`(zDJ7GM>y}7_wBijC=c9V`DgjH19{QN(ZX11|
zKS)M;0~-gw(wYg+tiRhwEO-pXhj$mE%CHJ1#Lhth;K5i*%B};M^|xKXyCcI$!?djl
zVEd~7tWVIm)DohF&$8E#j>{@Sg#0S8Pvv+>9cY_NU4#0Dq+o+DV2+!`iENjmU9jpm
z6DFCI>nl$M+JSw1;Z|+Wc0lkNCx{yhl?gzlbLma24{F~SJ}39fJI(ls8D=MyTcx<f
z)dFYK@zl)4swlud7pl#+({n=zs<HC!W7fhPyD)03p`12xO;dsSIClMJCU{XAJaDYR
z?5O*>jiM?{HvPfYH;e9S;m>55ZEIQp0LdijmyKnp)joB4GcGQIZ2>r4J@f&w6-JHL
z>z51~bx5kL&*E9F9(`=w$=Dfl_W3g8<m0lMJ@$C#$#Gl<mW5{-)#MgiKc`ASz`^qK
zEwRVh49@p@9;*0e9!z^3tcSnX36LP+JTEzj;&>7XCL8d1(z-U_%kH|4ScBhg%6k>D
z568&_VAw5@1+k~B#WxH3O1Hx1KGji#gu73Qzj)3#(^8aC;o`lly_De*?r5zC+zCEC
zve`laN^JU>A5=^|Rwrjkpn3PGd++&XW4YKl9qd<ey;cKTSA2~>R_^Iym(BPbREqf@
zZCQ)LR+-kcNbF#rP|wk9h0@8!aH&Y|(%pi+xGQogJD}IMa=hd7BWKUgv&Z|6#|!-d
z-UoX;0ejZn{k7Y^ey)X&tj9e+%D>D`a($I8F$_Oct9tkrP`Oq<A1Nmbsew1_MeDOo
z=W;lFd_UGcu^V>ul|&4W9@7eW8PDtq-tl>2%F`e4rmZ~3D1%2|>=K8HUfpVOytmC{
z6N5WlUYJT6&kkp2W$xC~+QR`S+bPG(PadSX&d%4NDrD3<m8I{I@e?+pF$mm?*#+z?
zn~!uDVdxas_4t}TVBIl@@|Y1@J!4e+xmN7XBlqhre&X{fUXe}iJ83-5Bg5mq7}uW9
zFl`n(*{}JIC&3&^OdC{sh9{J|{L8;!Ithr}`$=f&>32^MFh3N6e2xt3<oo^WA0XiL
zYl=%`#3={@B({{YbOrKNg7}a@sX%%t-K%;a2X{NWum6C7n9^@ismL!g)c+7PdP?c{
z%%8ud@dQFm-GMk(Al;b^!UY)pYgZmc)8ER!XXg6-gOA^)b{eSqcP96h5|GdB+bn;J
zT1C=5SvL^FdP4ZykKZYO2g*qv%9D(x$Us@-?13>A5J0}q-=H~(dOb}~`uoQ|xIE^j
zkc&Vf9b`((lN?rmNSz5h1D!k+6$SETO8oD3k`kxw{s6J3i=`DV7{tei3;yj7h@X!S
z$R}n6`mGN6@h|cZ$l<p-2#^fs`#LCu51HxkyLM2B0P>T-@9Mx{r~m-WcXdEt_@C>*
zU_Smoj|T?x3;b9o00%hhyYawaI1tPHr#j$ErXTu3Ajq6}-}eQ+$ppfe|5V3|j9~s#
z9gGk7h~rOnd;&kV<NtAO5cp5#0*C+DS5WYWIYYq$KhOmL1U}07elAeRj~GCq$mB=g
zwFCBu|A##SV!nUa4=5a&S?;^OVBV8Z`ajjde&iQk{vY><R}dMg|9xK=FzFxaV1hrk
zgCjHHeczW4g3Lzur#cuiqurnC;K;P2-`DX&ks0p(RLA$@I{4wp%zWRs6M!Nm+n?(A
z1b)~L0XQ;4-uLZ*sr;~)aM+JL4i^A&)%~fjAOx8U@%!-vVLxIZ2>*cwaC(8{gWvas
z{N%hrzyd(h!SCBaA+R6LG*E}kYIr&qHy2Zc9gsW}DY!Kd-j+!D0(>~&2#CC|A^%q#
htQ<iAvws!$u5P9-ZYPom_-36Sj(h&Rw2BPw{{!fHQ$YX#

diff --git a/Detectors/Upgrades/ITS3/macros/test/CheckChipResponseFile.C b/Detectors/Upgrades/ITS3/macros/test/CheckChipResponseFile.C
index 5cacf2eab9900..32d5bad87ce21 100644
--- a/Detectors/Upgrades/ITS3/macros/test/CheckChipResponseFile.C
+++ b/Detectors/Upgrades/ITS3/macros/test/CheckChipResponseFile.C
@@ -57,25 +57,25 @@ void LoadRespFunc()
   std::string AptsFile = "$(O2_ROOT)/share/Detectors/Upgrades/ITS3/data/ITS3ChipResponseData/APTSResponseData.root";
   std::string AlpideFile = "$(O2_ROOT)/share/Detectors/ITSMFT/data/AlpideResponseData/AlpideResponseData.root";
 
-  std::cout<<"=====================\n";
+  std::cout << "=====================\n";
   LOGP(info, "ALPIDE Vbb=0V response");
   mAlpSimResp0 = loadResponse(AlpideFile, "response0"); // Vbb=0V
   mAlpSimResp0->computeCentreFromData();
   mAlpSimResp0->print();
   LOGP(info, "Response Centre {}", mAlpSimResp0->getRespCentreDep());
-  std::cout<<"=====================\n";
+  std::cout << "=====================\n";
   LOGP(info, "ALPIDE Vbb=-3V response");
   mAlpSimResp1 = loadResponse(AlpideFile, "response1"); // Vbb=-3V
   mAlpSimResp1->computeCentreFromData();
   mAlpSimResp1->print();
   LOGP(info, "Response Centre {}", mAlpSimResp1->getRespCentreDep());
-  std::cout<<"=====================\n";
+  std::cout << "=====================\n";
   LOGP(info, "APTS response");
   mAptSimResp1 = loadResponse(AptsFile, "response1"); // APTS
   mAptSimResp1->computeCentreFromData();
   mAptSimResp1->print();
   LOGP(info, "Response Centre {}", mAptSimResp1->getRespCentreDep());
-  std::cout<<"=====================\n";
+  std::cout << "=====================\n";
 }
 
 std::vector<float> getCollectionSeediciencies(o2::its3::ChipSimResponse* resp,
@@ -85,7 +85,7 @@ std::vector<float> getCollectionSeediciencies(o2::its3::ChipSimResponse* resp,
   bool flipRow = false, flipCol = false;
   for (auto depth : depths) {
     auto rspmat = resp->getResponse(0.0, 0.0,
-                                    um2cm(depth) + resp->getDepthMin() + 1.e-9,
+                                    um2cm(depth) + 1.e-9,
                                     flipRow, flipCol);
     seed.push_back(rspmat ? rspmat->getValue(2, 2) : 0.f);
   }
@@ -99,7 +99,7 @@ std::vector<float> getShareValues(o2::its3::ChipSimResponse* resp,
   bool flipRow = false, flipCol = false;
   for (auto depth : depths) {
     auto rspmat = resp->getResponse(0.0, 0.0,
-                                    um2cm(depth) + resp->getDepthMin() + 1.e-9,
+                                    um2cm(depth) + 1.e-9,
                                     flipRow, flipCol);
     float s = 0;
     int npix = resp->getNPix();
@@ -121,7 +121,7 @@ std::vector<float> getEffValues(o2::its3::ChipSimResponse* resp,
   bool flipRow = false, flipCol = false;
   for (auto depth : depths) {
     auto rspmat = resp->getResponse(0.0, 0.0,
-                                    um2cm(depth) + resp->getDepthMin() + 1.e-9,
+                                    um2cm(depth) + 1.e-9,
                                     flipRow, flipCol);
     float s = 0;
     int npix = resp->getNPix();
@@ -140,9 +140,12 @@ void CheckChipResponseFile()
   LoadRespFunc();
   LOG(info) << "Response function loaded" << std::endl;
 
-  std::vector<float> vecDepth(50);
-  for (int i = 0; i < 50; ++i)
-    vecDepth[i] = i;
+  std::vector<float> vecDepth;
+  int numPoints = 100;
+  for (int i = 0; i < numPoints; ++i) {
+    float value = -50 + i * (100.0f / (numPoints - 1));
+    vecDepth.push_back(value);
+  }
 
   int colors[] = {kOrange + 7, kRed + 1, kAzure + 4};
   struct RespInfo {
@@ -156,7 +159,7 @@ void CheckChipResponseFile()
     {mAlpSimResp1, "ALPIDE Vbb=-3V", colors[2]}};
 
   TCanvas* c1 = new TCanvas("c1", "c1", 800, 600);
-  TH1* frame = c1->DrawFrame(-1, -0.049, 50, 1.049);
+  TH1* frame = c1->DrawFrame(-50, -0.049, 50, 1.049);
   frame->SetTitle(";Depth(um);Charge Collection Seed / Share / Eff");
   TLegend* leg = new TLegend(0.15, 0.5, 0.4, 0.85);
   leg->SetFillStyle(0);
diff --git a/Detectors/Upgrades/ITS3/simulation/src/ChipSimResponse.cxx b/Detectors/Upgrades/ITS3/simulation/src/ChipSimResponse.cxx
index 1c4f08438d4f2..72b291fb0d653 100644
--- a/Detectors/Upgrades/ITS3/simulation/src/ChipSimResponse.cxx
+++ b/Detectors/Upgrades/ITS3/simulation/src/ChipSimResponse.cxx
@@ -25,61 +25,66 @@ void ChipSimResponse::initData(int tableNumber, std::string dataPath, const bool
 
 void ChipSimResponse::computeCentreFromData()
 {
-  std::vector<float> zVec, qVec;
   const int npix = o2::itsmft::AlpideRespSimMat::getNPix();
+  std::vector<float> zVec, effVec;
+  zVec.reserve(mNBinDpt);
+  effVec.reserve(mNBinDpt);
 
   for (int iz = 0; iz < mNBinDpt; ++iz) {
-    size_t bin = iz + mNBinDpt * (0 + mNBinRow * 0);
-    const auto& mat = mData[bin];
-    float val = mat.getValue(npix / 2, npix / 2);
-    float gz = mDptMin + iz / mStepInvDpt;
-    zVec.push_back(gz);
-    qVec.push_back(val);
+    int rev = mNBinDpt - 1 - iz;
+    float z = mDptMin + iz / mStepInvDpt;
+    float sum = 0.f;
+    const auto& mat = mData[rev];
+    for (int ix = 0; ix < npix; ++ix) {
+      for (int iy = 0; iy < npix; ++iy) {
+        sum += mat.getValue(ix, iy);
+      }
+    }
+    zVec.push_back(z);
+    effVec.push_back(sum);
   }
 
-  std::vector<std::pair<float, float>> zqPairs;
-  for (size_t i = 0; i < zVec.size(); ++i) {
-    zqPairs.emplace_back(zVec[i], qVec[i]);
-  }
-  std::sort(zqPairs.begin(), zqPairs.end());
-  zVec.clear();
-  qVec.clear();
-  for (auto& p : zqPairs) {
-    zVec.push_back(p.first);
-    qVec.push_back(p.second);
-  }
+  struct Bin {
+    float z0, z1, q0, q1, dq;
+  };
+  std::vector<Bin> bins;
+  bins.reserve(zVec.size() - 1);
 
-  struct BinInfo { float z0, z1, q0, q1, dq; };
-  std::vector<BinInfo> bins;
   float totQ = 0.f;
   for (size_t i = 0; i + 1 < zVec.size(); ++i) {
     float z0 = zVec[i], z1 = zVec[i + 1];
-    float q0 = qVec[i], q1 = qVec[i + 1];
-    float dz = z1 - z0;
-    float dq = 0.5f * (q0 + q1) * dz;
+    float q0 = effVec[i], q1 = effVec[i + 1];
+    float dq = 0.5f * (q0 + q1) * (z1 - z0);
     bins.push_back({z0, z1, q0, q1, dq});
     totQ += dq;
   }
-  if (totQ <= 0.f) { mRespCentreDep = 0.f; return; }
+
+  if (totQ <= 0.f) {
+    mRespCentreDep = mDptMin;
+    return;
+  }
 
   float halfQ = 0.5f * totQ;
   float cumQ = 0.f;
-  for (const auto& b : bins) {
-    if (cumQ + b.dq < halfQ) { cumQ += b.dq; continue; }
-    float qSlope = (b.q1 - b.q0) / (b.z1 - b.z0);
+  for (auto& b : bins) {
+    if (cumQ + b.dq < halfQ) {
+      cumQ += b.dq;
+      continue;
+    }
     float dz = b.z1 - b.z0;
-    float A = qSlope * 0.5f;
-    float B = b.q0;
-    float C = cumQ - halfQ;
-    float disc = B * B - 4.f * A * C;
+    float slope = (b.q1 - b.q0) / dz;
+    float disc = b.q0 * b.q0 - 2.f * slope * (cumQ - halfQ);
+
     float x;
-    if (disc >= 0.f && std::abs(A) > 1.e-12f)
-      x = (-B + std::sqrt(disc)) / (2.f * A);
-    else
+    if (disc >= 0.f && std::abs(slope) > 1e-6f) {
+      x = (-b.q0 + std::sqrt(disc)) / slope;
+    } else {
       x = (halfQ - cumQ) / b.q0;
+    }
     x = std::clamp(x, 0.f, dz);
     mRespCentreDep = b.z0 + x;
     return;
   }
+
   mRespCentreDep = mDptMax;
 }
diff --git a/core_dump_74974 b/core_dump_74974
deleted file mode 100644
index e69de29bb2d1d..0000000000000