Removing obsolete mass hypothesis (M.Ivanov)
[u/mrichter/AliRoot.git] / TPC / AliTPCseed.cxx
index 3d802591f633af6103f7c0fea163737a1d4b8f27..bc323871acef6795d2629f1c08db9350259da4d9 100644 (file)
@@ -50,6 +50,7 @@ AliTPCseed::AliTPCseed():AliTPCtrack(){
   fCurrentCluster =0;
   fCurrentSigmaY2=0;
   fCurrentSigmaZ2=0;
+  fEsd =0;
   fCircular = 0;  // not curling track
 }
 AliTPCseed::AliTPCseed(const AliTPCseed &s):AliTPCtrack(s){
@@ -62,6 +63,7 @@ AliTPCseed::AliTPCseed(const AliTPCseed &s):AliTPCtrack(s){
   fPoints  = 0;
   fEPoints = 0;
   fCircular =0;
+  fEsd =0;
 }
 AliTPCseed::AliTPCseed(const AliTPCtrack &t):AliTPCtrack(t){
   //
@@ -73,7 +75,7 @@ AliTPCseed::AliTPCseed(const AliTPCtrack &t):AliTPCtrack(t){
   fRemoval =0;
   fSort =0;
   for (Int_t i=0;i<3;i++)   fKinkIndexes[i]=t.GetKinkIndex(i);
-
+  for (Int_t i=0;i<5;i++)   fTPCr[i]=0.2;
   for (Int_t i=0;i<160;i++) {
     fClusterPointer[i] = 0;
     Int_t index = t.GetClusterIndex(i);
@@ -93,6 +95,7 @@ AliTPCseed::AliTPCseed(const AliTPCtrack &t):AliTPCtrack(t){
   fCurrentSigmaY2=0;
   fCurrentSigmaZ2=0;
   fCircular =0;
+  fEsd =0;
 }
 
 AliTPCseed::AliTPCseed(UInt_t index,  const Double_t xx[5], const Double_t cc[15], 
@@ -124,6 +127,7 @@ AliTPCseed::AliTPCseed(UInt_t index,  const Double_t xx[5], const Double_t cc[15
   fCurrentCluster =0;
   fCurrentSigmaY2=0;
   fCurrentSigmaZ2=0;
+  fEsd =0;
 }
 
 AliTPCseed::~AliTPCseed(){
@@ -332,7 +336,7 @@ Int_t  AliTPCseed::GetProlongation(Double_t xk, Double_t &y, Double_t & z) const
 
 
 //_____________________________________________________________________________
-Double_t AliTPCseed::GetPredictedChi2(const AliTPCclusterMI *c) const 
+Double_t AliTPCseed::GetPredictedChi2(const AliCluster *c) const 
 {
   //-----------------------------------------------------------------
   // This function calculates a predicted chi2 increment.
@@ -392,7 +396,7 @@ Int_t AliTPCseed::Compare(const TObject *o) const {
 
 
 //_____________________________________________________________________________
-Int_t AliTPCseed::Update(const AliTPCclusterMI *c, Double_t chisq, UInt_t /*index*/) {
+Int_t AliTPCseed::Update(const AliCluster *c, Double_t chisq, UInt_t /*index*/) {
   //-----------------------------------------------------------------
   // This function associates a cluster with this track.
   //-----------------------------------------------------------------
@@ -450,7 +454,7 @@ Int_t AliTPCseed::Update(const AliTPCclusterMI *c, Double_t chisq, UInt_t /*inde
 
 
 //_____________________________________________________________________________
-void AliTPCseed::CookdEdx(Double_t low, Double_t up,Int_t i1, Int_t i2, Bool_t onlyused) {
+Float_t AliTPCseed::CookdEdx(Double_t low, Double_t up,Int_t i1, Int_t i2, Bool_t onlyused) {
   //-----------------------------------------------------------------
   // This funtion calculates dE/dX within the "low" and "up" cuts.
   //-----------------------------------------------------------------
@@ -660,7 +664,7 @@ void AliTPCseed::CookdEdx(Double_t low, Double_t up,Int_t i1, Int_t i2, Bool_t o
   }
   else{
     SetdEdx(0);
-    return;
+    return 0;
   }
   //  Float_t dedx1 =dedx;
   /*
@@ -677,27 +681,7 @@ void AliTPCseed::CookdEdx(Double_t low, Double_t up,Int_t i1, Int_t i2, Bool_t o
 
   
   SetdEdx(dedx);
-    
-  //mi deDX
-
-
-
-  //Very rough PID
-  Double_t p=TMath::Sqrt((1.+ GetTgl()*GetTgl())/(Get1Pt()*Get1Pt()));
-
-  if (p<0.6) {
-    if (dedx < 39.+ 12./(p+0.25)/(p+0.25)) { SetMass(0.13957); return;}
-    if (dedx < 39.+ 12./p/p) { SetMass(0.49368); return;}
-    SetMass(0.93827); return;
-  }
-
-  if (p<1.2) {
-    if (dedx < 39.+ 12./(p+0.25)/(p+0.25)) { SetMass(0.13957); return;}
-    SetMass(0.93827); return;
-  }
-
-  SetMass(0.13957); return;
-
+  return dedx;
 }
 Double_t AliTPCseed::Bethe(Double_t bg){
   //