From 401eff160cc8c87cdf8067cf73e6612bf3e4f861 Mon Sep 17 00:00:00 2001 From: masera Date: Fri, 13 Jun 2008 16:39:33 +0000 Subject: [PATCH] In AliITSRecoParam has been implemented the possibility of adding to the ITS cluster errors an additional contribution to account for large misalignments. This is needed for tracking the cosmics in ITS alone, and will be needed for tracking the cosmics in ITS+TPC. The additional error can be added as: AliITSRecoParam::SetClusterMisalError(XXXX); where XXXX is the error in cm (few mm might be needed). The error is then added for both y and z cluster coordinates. (A. Dainese) --- ITS/AliITSClusterParam.cxx | 17 +++++++++++++---- ITS/AliITSRecoParam.cxx | 1 + ITS/AliITSRecoParam.h | 8 +++++++- ITS/Calib/RecoParam/Run0_999999999_v0_s0.root | Bin 5092 -> 5145 bytes ITS/MakeITSRecoParam.C | 4 ++-- 5 files changed, 23 insertions(+), 7 deletions(-) diff --git a/ITS/AliITSClusterParam.cxx b/ITS/AliITSClusterParam.cxx index a52e20bcf4e..2c247421fc4 100644 --- a/ITS/AliITSClusterParam.cxx +++ b/ITS/AliITSClusterParam.cxx @@ -117,21 +117,30 @@ Int_t AliITSClusterParam::GetError(Int_t layer, // // Calculate cluster position error // + Int_t retval=0; switch(AliITSReconstructor::GetRecoParam()->GetClusterErrorsParam()) { case 0: - return GetErrorOrigRecPoint(cl,erry,errz); + retval = GetErrorOrigRecPoint(cl,erry,errz); break; case 1: - return GetErrorParamMI(layer,cl,tgl,tgphitr,expQ,erry,errz); + retval = GetErrorParamMI(layer,cl,tgl,tgphitr,expQ,erry,errz); break; case 2: - return GetErrorParamAngle(layer,cl,tgl,tgphitr,erry,errz); + retval = GetErrorParamAngle(layer,cl,tgl,tgphitr,erry,errz); break; default: - return GetErrorParamMI(layer,cl,tgl,tgphitr,expQ,erry,errz); + retval = GetErrorParamMI(layer,cl,tgl,tgphitr,expQ,erry,errz); break; } + // add error due to misalignment (to be improved) + Float_t errmisal2 = AliITSReconstructor::GetRecoParam()->GetClusterMisalError() + *AliITSReconstructor::GetRecoParam()->GetClusterMisalError(); + erry = TMath::Sqrt(erry*erry+errmisal2); + errz = TMath::Sqrt(errz*errz+errmisal2); + + return retval; + } //-------------------------------------------------------------------------- Int_t AliITSClusterParam::GetErrorOrigRecPoint(const AliITSRecPoint*cl, diff --git a/ITS/AliITSRecoParam.cxx b/ITS/AliITSRecoParam.cxx index 6c32454024f..afb817cca7d 100644 --- a/ITS/AliITSRecoParam.cxx +++ b/ITS/AliITSRecoParam.cxx @@ -98,6 +98,7 @@ fXPassDeadZoneHits(0), fUseTGeoInTracker(3), fAllowSharedClusters(kTRUE), fClusterErrorsParam(1), +fClusterMisalError(0.0), fComputePlaneEff(kFALSE), fHistoPlaneEff(kFALSE), fReadPlaneEffFromOCDB(kFALSE), diff --git a/ITS/AliITSRecoParam.h b/ITS/AliITSRecoParam.h index 9394562947a..aa9dff6a224 100644 --- a/ITS/AliITSRecoParam.h +++ b/ITS/AliITSRecoParam.h @@ -125,6 +125,10 @@ class AliITSRecoParam : public AliDetectorRecoParam void SetClusterErrorsParam(Int_t param=1) { fClusterErrorsParam=param; return; } Int_t GetClusterErrorsParam() const { return fClusterErrorsParam; } + void SetClusterMisalError(Float_t err=0.) { fClusterMisalError=err; return; } + Float_t GetClusterMisalError() const { return fClusterMisalError; } + + void SetUseAmplitudeInfo(Bool_t use=kTRUE) { for(Int_t i=0;i}+}V%-KJ?&jxoEkAk;{hc5u|u@e9Q*Z}}^AOL{& z^;JLMDu}M)9r1s6fCvEK*2}9Q{VLxa6-I*$1}5c^u8xCON?hUqz!~75Whn6+wiST; z-+fms0D#0$-O1HnM+olWW-sJs<74k_BjoFCydI!BX#Frr=;~s{(=1k0sxpR z*1N0t3k3jZ&i`*7{~Y(9c8%qqgn*ZS_ly0L5RpZSE-j0wnGX3X5#SOC2q5yYw*E8H zNgkotN&JC|Qs?pyF#MS=)=|^*iL-a8qOHCySo0hNtVro0IziHEg4h|uOrBU@67|`B zsn^pawEythy!U$<(i&_dkbDA4d3<6weNpcb76WX0N$F{AJ@5kKj()A!0TQ)h%$z*V zqxBI@`LQ|=O#NlrWt;0MqK?)Lctj+73!Z&knQ|FU$0)poNtYECHdH#wX^2blaMmh% zY^(J^(?xXH;ECy-b&a8<;6wzkT3IXwI5rDW)>~%1Zl?s>qR*g9o{KFFCXHug-}yhL)8ORIvq0#a`S{YF9SAbGiUEf?#v+!=Ed1 zyYI6S7L~mMJ|y)ZZlC)Pd3aUm?0leRx|kW1^zGugw!@o|dw{ZVx)7XZs+$Af9xyM+J$wy-1Y!K%0oWU6;7} zUIm&lRgJ0um)RxAORDRAfK@YDeOVz@H=_ErVEv6lL}JkUkyM8`F|2I`yccGya>+Iao=t~4 zzL*G2C(hAA?xc)B^F8Xt$KsgXH7gMrvCEtB>gx?16#m3NK+{WDU%3jdh`EHm(5NF+ zKX�A-HdR37>H%LMq88QvOF+tfH#M5Jt=BCRgRSa}1LH^1H^5=G)CHas8o*N2{A1 zTFzWKhF(3~QUwef33|0B`@%sham#UhQ3|B1G!Z7q9-F|mG?G`pMF7gLdT{;1D-t$u zWVx#lr8wqPaNEddK_)}*nS{N`^*vE4y9g`0u{pSvQ>ig|-90>4BLrTeE`v){tfv7#%;_n&9rnSBG^7u!GZ7!zeOWTGkQJy-JrnQb?uoB;7A*n zHv-8f@O8<9u>39F)Q9su!lx(O=WbOLB{QfcRIg$GMB9LVC$zfG9c91k^=!HPc9L%}_DHFI(Nq&sCmBL@Dkbh3lEu~zDBXYm;i8!7 zZxX@527y7!#|FZj^SfwJ-{Js)s-Q+9KegIQE%`PtQx02bvM%*U(u@6>b5e))=Lw%c zBR!6+&%b&V3_g>Lj8a}=$$ete<)W5(bjqeZzy}tItec3H5G8rX+L!>R)7r_E8AgVk zmaB&h%^~MqFl#W1CLe zS)W%n9UO0OXyF-g^CEr7c@cMUmS+`Wm~N91db97Xr1-5U?)*#Cr_8V1%FX&zLZ|n& zy$rv3-5m*?Z|wfFTrA&Db8S-)R=rk@KY!?7E| z^_ssg@gJ_vtBrLJP9>)e>pjde&Ki!>!P+JNlt~FU!;*e8lcy3v{)XItNESMRTIy6F zhFyq$m{bFoYr)J##B!L>Z>GDUr51L~ll+46@9hZ-L+{e-=HFA~RBj)c@sHKQzzwi{ z$%&@e^e$-mjL?A~*60=exzB78!2cpP02CnU;Puu?XJL$@#lZ3x1iB0yMZ4*>*euq| zA)KHu0RZ0y{&K2R@?_nlVB&?J<~I6;h6saC+X*dWNp-P|3q z0sR;A5+6b-1o95cZb4|~?=kqHkyz(5%nxK|5}1quS0Gc)7MW2z<5h8v;xyj4>QF7 zoPsV*`xKnk4!o&1KPHu1PYHeUR*uv_w`#X4i@d?w&eu{aqI*7S4`C0g3pQf-QlAyq zs&X%oHX0@MCQp)OkH|jAU|759aMU>HP$Ra1=yEM{XOnl6{nfLv&>yXvu3nN5!!Sh7 z(TNw`^+&pa^$ku?XI*M;O`m2dforqbYh*g#<+0n@bf;3;v*gnDA9{OeGxMI@PEt9g z(6CSJSg?z|#aYq#=!r2+lTriWtP_}2o+-Qcn&UZs;rt0SB&n>U*HAEJcVEKD(nRb7 zbxrkcp<2G{!E+BmW{tH|H^LmKj~};TG<${^7F&@~n*2*HOuS-ZUbnl8yzWcm9j9iljg zyD8)vnd4QB12-QK+9*|#*w~IiRV(|VElAPtTrC}l+nxTNLXpH~rBj!@9Zr$ACK9H# zRg%{)v_G^nmMl4Qu*h-mbh&<8p8RHM4XKlsmJlz+NJt8Rj91|AFISoOPZ;Sn6NrLA{*oi)r=U+^{zsae9rzf z7ntzd-bMkpnQzmQBXFvMWp`(9(KAaxxur$ze@KqZm}9j%GbpjT0d#SX4Tp7R41Dq7 zK7Oj}1Ya~h*Ddf_F1;03%8rj$so9$U%t`TaIXHQgL?75Ni3e?cv~rGnH*^P)BABc_ z8FRCnpRoL5^m{ZO5Q!{l`SBGAYd~=%KQM=VedRrrPdAu|&3??z1;XJP^j?rqCD^ zg-C&KQLIkv5JQB(wfcY%w#Q)t_7U+-Lx^wQRBaZ!g0oLcOmx;Y9r}EcOzjqJSzdL5 zMeB>#Jes!yUNoDRI1~6iVgfe5DnNQ(umqd{@_X~+5<`YgEjB!X_?7@x$>0}9mD?X2 zUt$91d@~@MP{una^7g|U^sd*Gufw}knP^eiUT8JO}l zWv=WVo)zN-{U_ht3N)HUu5qZcu7X(hVCP)NL|-o8K@&Wm3zKi=K{e|AxbPhFXHD#PW^E%wH*_d^!m=KMg|v^l`S_U;W;@mcx0x}~y7dGC>j!H)b{~2R0)EqN4w>W5Kjz`whoj^6#urdk-ZC z@=rS2i+gP-HZ&eXj4%qKPSk}a2oaKoA2|!{&GXW+M){JS1M9SUhCH9X(5vT61HD^K zwo`bNf1I8ojgV*kH0RGqL%lM$${v~4$~8}FwAhb~p)QNkFQp~Am~-wou`3NJ1jKLN zd-&aGaVD_9)wYJXT-lPirmrL8IlQPy8lHO1h2jR9+Vf2t9owl0N{P<)`en^w^cpMqqU+gT*9-^`0!R;l(yG}9Kq!Bj8xMyDWMcx1 z!Z!*{lh0?xT5^gdiG)5z}E`H-Kl0k zzIO-F)hp!Yrbq?vAF^4B@?Te00`pgoNdEnJ|J5T`2IIPn|jQe5% ce}NSL0seJ1{|o##p7iyrKR+GFUpy7?U+0{e9RL6T delta 4579 zcmai&bx@Si+xHickgf$5kd*H34uM^oMY@*;>4v3i0fCiTO1h;>5G4iaE&(M(LJ$y; z^g*9z=J(FL?>o<#`OdlKT<3h|y080>`##?s_H8yfuNNuAkunRa&zwa6GEjLj;4lI`Jv}HO@UPu0JsMH-G&fJaQO{*@$cLN z3jn}0P<3*JY4Sn6VP4*RZq^TzHJ^c=f|kCzfp!{7obvylf&(!B(jIj604jgb|5W&w zOFH;}7xH)gKLZ`~w<3(2f9Hk%R)l8pkhQ~8aCK5rR#b8K@e2GuV}J<+XAKQ-g8>LY z^R}|OpJ4ic6D|MY(EuSn%@`}C>{jCIF>*t7wB3#sBo_St~ zKA6;;6iljMg0#+zY@(Difm&>PkAXQfGz~n|$Hu@sT9~=)cukV^9i^4Rolk}2w2J_* zobJBt{3ab7GPY3JzjbFXHPBxUx2NjD3NGH(CrH7uYsR7V&o- zc`qD44(4MNn6eWRUU`bMMgDQW22QkoodMs5Pgj94tR2Uv#~qY!Hr+aK8f1b#B|^p1j0RL!4@h~InWk-(l)ykgL5u?o7YFVKSa56&r4Dpjl~?Fm z1(FmUVYP}jIusHi1$J2erBvc0tuz!?p7};`p@Vfw$S?*z;u37PGv4+_{2h%fBAze{yj4w_WLyNLgj=?6VM-IFux#|X z$@QMB=xuNlTUo^8^8Tj!=pEU80$Ny?`8Ip`WlsgzPaow0CsC}m5t2Fi)IA(Z2`C!= zI%}I=>DX}V*Tskc=F_NVmWm-99Vvs ze@)mE;}@D+J|<^BbDkv!O1qvesT39Ic;9?9F#3ri@w&HAer(*!QO++Yy`Wegz8DiK z(ep_SWRbIl7W$OV*D$~ygr@t7rU|`F3UyD{d5KspVm2`9+k-WGe6{pB1Tsj!C0-1X@8yr2dhzHuqFqERq@G$3KfT7rlLr) zcbX?FXg4*Xd95h5fhtut3c`|nE^JF6u)s@JvnC-lEs94X7;_1gzBp9AxNT4TsHJ76 z{kf*jWwAJ6UN)F#-{?p4VlmR8)fS#-^N|Ynv18UJX?!9Y*H*;Tama!d7{0 zFh>|U+h}i<8qj18$W?B zj)b?>z^({5PHhk@z6xQcgRB*{wHNg?STQ8?FM4=fH9oQP+r$wUqBD~|`wd*zrs7Gb zI_XK1&M}j~0p&$j6Ljq|llqX#Y(o{ppEMKnXh+w!C!{~M>^wW1AS;u&oo43OF#UCm ziEuaV4(rb?(&kyJ{b9P+@O{$X&Is7*jPLf&>F+8=tK&W`;9g>`d)oH2?+mJdR3b!b z5_Ke+<}Z${&@Fi+#M3YVx=#e10a*ltp!y~Apt?_vmV(oMY&|!=kLHN3rlz zH|1_4RR%Qotm*e~7te(Bw&w$4ypA={9(yX{_jBZOR*A;~DvPBZtJB@F!yL@+bac=B z6rV{Gu5r<3NL)K3S!7?0%>>DttN+3ppV%sSDO+U|W<&|Bs#*LNjM`s?ocoqbI_bey z3Fl}EXvEx9eq4VXZsV~5ve#4(=UtCkYr*VZj6J@sKI8<#5i(_Zp4QYJ#J832+ihOH z5NP1u4FI!t$Gl#59C99O4iVWFkE`kh`$1*L)A?u13dur?@uuF6)pdNT*TG|og38TR zUvt##n;T}>Ry*f<(<<>gMo70tBj)pbWfm1{-ucGVPM<@wpOi!2L7!TWJ#!mypzAJ@ znlHVjO;Zzdz_8SYrS6O+`c7U&11dhbp{iB zOFFL9R!eCbH z;t|Sty|LGJG7Y)H#_jw{9|Ih@g-ryr&)QMTabt%V_TO^H2#JN18>^q););x3X0*#U zHV6)K@IT$8VlVBDl|tNPb7AqML=f!1&B~AXnTz@fWxW+7qZ4?1T(Y@KC5Rb#S&1>z zFRKaP>9W!=g_@a;1uAK6l^)wWDy6^qbDIh?pS?P3>}?`YLgO%{b-XWAkIq4RyxJ zbrmng@NQUmE`>+etUBq%o_OMog!WZyFB>|2u^7V&xY$2Q$Yr{7rZQ?q;J2W=GeyB^ zSkLpdoUG=a-@eTo9jU}VL-yAtDrsjdud)QSuAAv`VCp_G-aAU_Jat%MlCwc{_fs@N z5X3am&tJX5G?yh>K?~ZiWY{aKvnVs45zf7CGEB+<1~J*vi1itwT5pm3rc9E;~I$#NDG;x(jg z$gD+twKMJxK!lZ|9i8HU9 z7&}xI@UOv?uAaum+{uG-LU0dUbmVPh-PH>x_F-kAa@(E8eo6V|@u-y0ofrj(W@Z@Od9;dDU_=5mLOi zXw(&)Nn(DxOK({sfxS07G+zpzJOiIEr{Jj-Le1o_yDR z=Y%<44o#}+#b#%=i=LedOPulp8mp=%*O&n5Kbzehe+tD5NF+j1nn#CUt>+0!^idED zSTIH#iC{$!-@yvT({F@wcSO%l@LX~R-fM>}#bNjJcUpc4qGG+g*XxH*g`88+VDa^? zgk@3og*^7_;^+->7)x|Hp!vP@$nT|XMSDA!dv00`zjaCO>Bu!tPA}&(+KiE$$Qj=#t-Z(n~Y-e5Sf0aAqTtXGchgDvBE1Rfum> z@q45xzJLFiTSKdNyeydy*Ynn3c1=O)yPzpd7Ou@7UJ>OkEJbD*=$a`uxVhUnjjWvUu4dP*P@^!rN)^&@tB3rO!_bveb5n~#=N=I#@j?UM*H;EHg7)RfJ6Wga(0 z&{GugyLayjDZ?xyHyplDVoPYB4@d5F>dP3I8s63ZUTnGLYkX|K{HQAH&gA}HHuL?) zz;jRVflN~b4Cgy;Ae^L|8mzJKI)Otbh%HAwC#O`EuWmapyz;^j0x(}4gyiCyr&Vw%kyFISpKP1$ zKyXG*mnFmv9rMr?v@K*9qm}x!J{C&hPbB%;F7*q0T+`f?2@VeDBC#<4Y~c#VA(qUI3iihMpK z#w7om^n?DmKidVc-FJtihVdSoFE*a-!B;5%=YNR(yF>r;KOP*42e;vW^FOE`Ha%Rp i|6gVWGMbD*IPf19_dnQw*3AEk#r+E<+x8cE1^6$y>}M+g diff --git a/ITS/MakeITSRecoParam.C b/ITS/MakeITSRecoParam.C index 7ff0188773b..52a3d6124ba 100644 --- a/ITS/MakeITSRecoParam.C +++ b/ITS/MakeITSRecoParam.C @@ -31,9 +31,9 @@ void MakeITSRecoParam(Int_t type=1) { return; break; } + itsRecoParam->SetClusterErrorsParam(2); + //itsRecoParam->SetClusterMisalError(1.0); // [cm] /* - itsRecoParam->SetUseTGeoInTracker(3); - itsRecoParam->SetClusterErrorsParam(0); itsRecoParam->SetFindV0s(kTRUE); itsRecoParam->SetAddVirtualClustersInDeadZone(kFALSE); itsRecoParam->SetLayerToSkip(0); -- 2.43.0