Usage of the error and shape parameterization

[u/mrichter/AliRoot.git] / TPC / AliTPCTracklet.cxx

diff --git a/TPC/AliTPCTracklet.cxx b/TPC/AliTPCTracklet.cxx
index e8b77d51ff86bc7778b8f74aaf111350d1bc25b8..9eddd5a84c3d1fdacedd6e31ff061c5194a6cc0e 100755 (executable)
--- a/TPC/AliTPCTracklet.cxx
+++ b/TPC/AliTPCTracklet.cxx
@@ -326,6 +326,9 @@ void AliTPCTracklet::FitKalman(const AliTPCseed *seed,Int_t sector) {
     if (!c) continue;
     Double_t r[3]={c->GetX(),c->GetY(),c->GetZ()};
     Double_t cov[3]={0.01,0.,0.01}; //TODO: correct error parametrisation
+    AliTPCseed::GetError(c, track,cov[0],cov[2]);
+    cov[0]*=cov[0];
+    cov[2]*=cov[2];
     if (!track->PropagateTo(r[0])) {
       isOK=kFALSE;
       break;
@@ -342,7 +345,10 @@ void AliTPCTracklet::FitKalman(const AliTPCseed *seed,Int_t sector) {
     AliTPCclusterMI *c=track->GetClusterPointer(i);
     if (!c) continue;
     Double_t r[3]={c->GetX(),c->GetY(),c->GetZ()};
-    Double_t cov[3]={0.01,0.,0.01}; //TODO: correct error parametrisation
+    Double_t cov[3]={0.01,0.,0.01}; 
+    AliTPCseed::GetError(c, track,cov[0],cov[2]);
+    cov[0]*=cov[0];
+    cov[2]*=cov[2];
     if (!track->PropagateTo(r[0])) {
       isOK=kFALSE;
       break;
@@ -359,7 +365,10 @@ void AliTPCTracklet::FitKalman(const AliTPCseed *seed,Int_t sector) {
     AliTPCclusterMI *c=track->GetClusterPointer(i);
     if (!c) continue;
     Double_t r[3]={c->GetX(),c->GetY(),c->GetZ()};
-    Double_t cov[3]={0.01,0.,0.01}; //TODO: correct error parametrisation
+    Double_t cov[3]={0.01,0.,0.01}; 
+    AliTPCseed::GetError(c, track,cov[0],cov[2]);
+    cov[0]*=cov[0];
+    cov[2]*=cov[2];
     if (!track->PropagateTo(r[0])) {
       isOK=kFALSE;
       break;