Correctly resetting clusters (Yu.Belikov)

diff --git a/TPC/AliTPCtrack.h b/TPC/AliTPCtrack.h
index c630f756da4c9f3e908a513b5d1edc3eb6c6b315..e05077a6bc92377ad3d0905cec51dc64da073a63 100644 (file)
--- a/TPC/AliTPCtrack.h
+++ b/TPC/AliTPCtrack.h
@@ -105,6 +105,7 @@ public:
   Int_t PropagateTo(Double_t xr,Double_t x0=28.94,Double_t rho=0.9e-3);
   Int_t Update(const AliCluster* c, Double_t chi2, UInt_t i);
   void ResetCovariance();
+  void ResetClusters() {SetNumberOfClusters(0); SetChi2(0.);}
   void UpdatePoints();//update points 
   Float_t* GetPoints() {return fPoints;}
   //

diff --git a/TPC/AliTPCtracker.cxx b/TPC/AliTPCtracker.cxx
index c6ffce2edcf2f67491d5b7d57640de1dd4d185d0..051d718b3528b6a25a2240037ccafe6e406c4589 100644 (file)
--- a/TPC/AliTPCtracker.cxx
+++ b/TPC/AliTPCtracker.cxx
@@ -634,6 +634,7 @@ Int_t AliTPCtracker::ReadSeeds(const TFile *inp) {
   Int_t n=(Int_t)seedTree->GetEntries();
   for (Int_t i=0; i<n; i++) {
      seedTree->GetEvent(i);
+     seed->ResetClusters();
      fSeeds->AddLast(new AliTPCseed(*seed));
   }
   
@@ -734,7 +735,7 @@ Int_t AliTPCtracker::RefitInward(AliESD* event) {
       if ( (status & AliESDtrack::kTRDrefit ) == 0 ) continue;
 
     AliTPCtrack* tpcTrack = new AliTPCtrack(*track);
-    AliTPCseed* seed = new AliTPCseed(*tpcTrack);
+    AliTPCseed* seed=new AliTPCseed(*tpcTrack); seed->ResetClusters(); 
 
     if ( (status & AliESDtrack::kTRDrefit) == 0 ) seed->ResetCovariance();
 
@@ -785,7 +786,7 @@ Int_t AliTPCtracker::PropagateBack(AliESD *event) {
        if ( (status & AliESDtrack::kITSout) == 0 ) continue;
 
     const AliTPCtrack t(*esd);
-    AliTPCseed s(t);
+    AliTPCseed s(t); s.ResetClusters();
 
     if ( (status & AliESDtrack::kITSout) == 0 ) s.ResetCovariance();
 

diff --git a/TPC/AliTPCtrackerMI.cxx b/TPC/AliTPCtrackerMI.cxx
index e9c528a4fefabdb5bff7b37e8724dfce539de0a2..5e63f8cd7e504c2bf1fc261bafcb90ba6f031c11 100644 (file)
--- a/TPC/AliTPCtrackerMI.cxx
+++ b/TPC/AliTPCtrackerMI.cxx
@@ -2761,7 +2761,7 @@ void AliTPCtrackerMI::ReadSeeds(AliESD *event, Int_t direction)
     ULong_t status=esd->GetStatus();
     if (!(status&AliESDtrack::kTPCin)) continue;
     AliTPCtrack t(*esd);
-    AliTPCseed *seed = new AliTPCseed(t/*,t.GetAlpha()*/);
+    AliTPCseed *seed=new AliTPCseed(t/*,t.GetAlpha()*/);seed->ResetClusters();
     for (Int_t ikink=0;ikink<3;ikink++) seed->GetKinkIndexes()[ikink] = esd->GetKinkIndex(ikink);
     if ((status==AliESDtrack::kTPCin)&&(direction==1)) seed->ResetCovariance(); 
     if ( direction ==2 &&(status & AliESDtrack::kTRDrefit) == 0 ) seed->ResetCovariance();
@@ -4551,6 +4551,7 @@ Int_t AliTPCtrackerMI::ReadSeeds(const TFile *inp) {
   Int_t n=(Int_t)seedTree->GetEntries();
   for (Int_t i=0; i<n; i++) {
      seedTree->GetEvent(i);
+     seed->ResetClusters();
      fSeeds->AddLast(new AliTPCseed(*seed/*,seed->GetAlpha()*/));
   }