]> git.uio.no Git - u/mrichter/AliRoot.git/commitdiff
re-introduce part on covariance matrix for real data
authorcblume <cblume@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 5 Aug 2009 16:41:01 +0000 (16:41 +0000)
committercblume <cblume@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 5 Aug 2009 16:41:01 +0000 (16:41 +0000)
TRD/AliTRDseedV1.cxx

index d84b11dca513bf45977f102314fd59e57d936a7b..1cd77b3c7006a0eac8dfb14bcd87a2a047b13358 100644 (file)
@@ -898,6 +898,7 @@ Bool_t      AliTRDseedV1::AttachClusters(AliTRDtrackingChamber *chamber, Bool_t tilt)
   Int_t t0 = 14;
   Int_t kClmin = Int_t(fReconstructor->GetRecoParam() ->GetFindableClusters()*AliTRDtrackerV1::GetNTimeBins());
 
+  Double_t sysCov[5]; fReconstructor->GetRecoParam()->GetSysCovMatrix(sysCov); 
   Double_t s2yTrk= fRefCov[0], 
            s2yCl = 0., 
            s2zCl = GetPadLength()*GetPadLength()/12., 
@@ -938,7 +939,7 @@ Bool_t      AliTRDseedV1::AttachClusters(AliTRDtrackingChamber *chamber, Bool_t tilt)
     // get standard cluster error corrected for tilt
     cp.SetLocalTimeBin(it);
     cp.SetSigmaY2(0.02, fDiffT, fExB, dx, -1./*zt*/, fYref[1]);
-    s2yCl = (cp.GetSigmaY2() + t2*s2zCl)/(1.+t2);
+    s2yCl = (cp.GetSigmaY2() + sysCov[0] + t2*s2zCl)/(1.+t2);
     // get estimated road
     kroady = 3.*TMath::Sqrt(12.*(s2yTrk + s2yCl));