Store correct chi2/ndf for tpc constraint tracks

diff --git a/ANALYSIS/AliAnalysisTaskESDfilter.cxx b/ANALYSIS/AliAnalysisTaskESDfilter.cxx
index ea75a23e1c85c19b2e05138bdb11b5cddb1b817c..60132584cdd9cba181a987b660462966f3357941 100644 (file)
--- a/ANALYSIS/AliAnalysisTaskESDfilter.cxx
+++ b/ANALYSIS/AliAnalysisTaskESDfilter.cxx
@@ -968,18 +968,11 @@ void AliAnalysisTaskESDfilter::ConvertTPCOnlyTracks(const AliESDEvent& esd)
       AliExternalTrackParam exParam;\r
       // take the B-field from the ESD, no 3D fieldMap available at this point\r
       Bool_t relate = false;\r
-      relate = track->RelateToVertex(vtxSPD,esd.GetMagneticField(),kVeryBig,&exParam);\r
+      relate = track->RelateToVertexTPC(vtxSPD,esd.GetMagneticField(),kVeryBig,&exParam);\r
       if(!relate){\r
         delete track;\r
         continue;\r
       }\r
-      /*\r
-      Double_t  fSigmaConstrainedMax = 5.;\r
-      if(track->GetConstrainedChi2TPC() > fSigmaConstrainedMax*fSigmaConstrainedMax) {\r
-       delete track;\r
-       continue;\r
-      }\r
-      */\r
       track->Set(exParam.GetX(),exParam.GetAlpha(),exParam.GetParameter(),exParam.GetCovariance());\r
     }\r
     \r
@@ -1005,7 +998,13 @@ void AliAnalysisTaskESDfilter::ConvertTPCOnlyTracks(const AliESDEvent& esd)
                                                             selectInfo);\r
     aodTrack->SetTPCClusterMap(track->GetTPCClusterMap());\r
     aodTrack->SetTPCSharedMap (track->GetTPCSharedMap());\r
-    aodTrack->SetChi2perNDF(Chi2perNDF(track));\r
+    Float_t ndf = track->GetTPCNcls()+1 - 5 ;\r
+    if(ndf>0){\r
+      aodTrack->SetChi2perNDF(track->GetTPCchi2()/ndf);\r
+    }\r
+    else{\r
+      aodTrack->SetChi2perNDF(-1);\r
+    }\r
     aodTrack->SetFlags(track->GetStatus());\r
     aodTrack->SetTPCPointsF(track->GetTPCNclsF());\r
 \r