From 92a3add6a9dc4d3609cf0a4e3c8ae0f0cec9b9ec Mon Sep 17 00:00:00 2001 From: miranov25 Date: Fri, 12 Jun 2020 15:21:52 +0200 Subject: [PATCH] ALIROOT-8439 - adding systematic error only after rescaling - keeping gold interface in ErrY2Z2 - needed by HLT --- TPC/TPCrec/AliTPCtracker.cxx | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/TPC/TPCrec/AliTPCtracker.cxx b/TPC/TPCrec/AliTPCtracker.cxx index 887e6f1da9f..0d3ccaf740e 100644 --- a/TPC/TPCrec/AliTPCtracker.cxx +++ b/TPC/TPCrec/AliTPCtracker.cxx @@ -311,8 +311,10 @@ Int_t AliTPCtracker::AcceptCluster(AliTPCseed * seed, AliTPCclusterMI * cluster) Double_t sy2=0, erry2LM=0, erry2HM=0, erry2Res=0,erry2OccudEdx=0, erry2Ratio=0;//ErrY2(seed,cluster); Double_t sz2=0, errz2LM=0, errz2HM=0, errz2Res=0,errz2OccudEdx=0, errz2Ratio=0;//ErrZ2(seed,cluster); ErrY2Z2(seed,cluster,sy2,sz2); - erry2Res=sy2; - errz2Res=sz2; + Float_t serry2=seed->GetErrorY2Syst(); + Float_t serrz2=seed->GetErrorZ2Syst(); + erry2Res=sy2-serry2; + errz2Res=sz2-serrz2; const Float_t kMinSigma2=0.02*0.02; if (AliTPCReconstructor::GetRecoParam()->GetUseClusterErrordEdxCorrection()||AliTPCReconstructor::GetRecoParam()->GetUseClusterErrordEdxMultCorrection()) { const Float_t kdEdxMIP=50., kClusterMIP=20.; @@ -372,8 +374,6 @@ Int_t AliTPCtracker::AcceptCluster(AliTPCseed * seed, AliTPCclusterMI * cluster) seed->SetErrorZ2(errz2Res); /// } // - Float_t serry2=seed->GetErrorY2Syst(); - Float_t serrz2=seed->GetErrorZ2Syst(); seed->SetErrorY2(erry2Res+serry2); seed->SetErrorZ2(errz2Res+serrz2); Double_t sdistancey2 = erry2Res+serry2+seed->GetSigmaY2(); @@ -950,8 +950,8 @@ void AliTPCtracker::ErrY2Z2(AliTPCseed* seed, const AliTPCclusterMI * cl, double // additional systematic error on the cluster double serry2=0,serrz2=0; AliTPCcalibDB::Instance()->GetTransform()->ErrY2Z2Syst(cl, tgp, seed->GetTgl(), serry2,serrz2); - // erry2 += serry2; // MI change -add error due distortion fluctuation later after rescaling - // errz2 += serrz2; // MI change -add error due distortion fluctuation later after rescaling + erry2 += serry2; + errz2 += serrz2; seed->SetErrorY2(erry2); seed->SetErrorZ2(errz2); seed->SetErrorY2Syst(serry2);