Merge branch 'feature-movesplit'

[u/mrichter/AliRoot.git] / PWGPP / TRD / AliTRDcheckTRK.cxx

diff --git a/PWGPP/TRD/AliTRDcheckTRK.cxx b/PWGPP/TRD/AliTRDcheckTRK.cxx
index fb94a2ebbad103b35877e22062111b69b489abeb..11a731dd7489837f42ea4a3baf8c9200302400ed 100644 (file)
--- a/PWGPP/TRD/AliTRDcheckTRK.cxx
+++ b/PWGPP/TRD/AliTRDcheckTRK.cxx
@@ -212,7 +212,7 @@ Bool_t AliTRDcheckTRK::PropagateKalman(AliTRDtrackV1 &t, AliExternalTrackParam *
     return kFALSE;
   }
   if(ref->Pt()<1.e-3) return kFALSE;
-
+  Float_t prod(t.GetBz()*t.Charge());
 
   // Initialize TRD track to the reference
   AliTRDtrackV1 tt;
@@ -240,11 +240,9 @@ Bool_t AliTRDcheckTRK::PropagateKalman(AliTRDtrackV1 &t, AliExternalTrackParam *
     }
     if(HasTrkltRefit()){
       //      if(!tr->FitRobust(tt.Charge()>0.)) printf("W - AliTRDcheckTRK::PropagateKalman :: FitRobust() failed for Det[%03d]\n", det);
-      if(!tr->FitRobust(AliTRDgeometry::GetPadPlane(det))) printf("W - AliTRDcheckTRK::PropagateKalman :: FitRobust() failed for Det[%03d]\n", det);
-      else {
-       TGeoHMatrix *matrix = AliTRDgeometry::GetClusterMatrix(det);
-       if (matrix) tr->SetXYZ(matrix);
-      }
+      TGeoHMatrix *matrix = AliTRDgeometry::GetClusterMatrix(det);
+      if(!tr->FitRobust(AliTRDgeometry::GetPadPlane(det), matrix, tt.GetBz(), tt.Charge())) printf("W - AliTRDcheckTRK::PropagateKalman :: FitRobust() failed for Det[%03d]\n", det);
+      else tr->SetXYZ(matrix);
     }
     if(!AliTRDtrackerV1::PropagateToX(tt, tr->GetX0(), fgKalmanStep)) continue;
     if(!tt.GetTrackIn()) tt.SetTrackIn();