kSPECIESN from AliESD: number of particle species identified by PHOS

diff --git a/PHOS/AliPHOSRecParticle.cxx b/PHOS/AliPHOSRecParticle.cxx
index 0ddfa7eec9da23128674df370a3f8c969f589860..df5777b8cd5809b78c365090a9770bd19328e5b7 100644 (file)
--- a/PHOS/AliPHOSRecParticle.cxx
+++ b/PHOS/AliPHOSRecParticle.cxx
@@ -36,8 +36,8 @@
   AliPHOSRecParticle::AliPHOSRecParticle(): fPHOSTrackSegment(0)  ,  fDebug( kFALSE )
 {
   // ctor
-  const Int_t nSPECIES = AliESDtrack::kSPECIES;
-  for(Int_t i = 0; i<nSPECIES+4 ; i++)
+  const Int_t nSPECIES = AliESDtrack::kSPECIESN;
+  for(Int_t i = 0; i<nSPECIES ; i++)
     fPID[i]=0.;
 }
 
@@ -72,8 +72,8 @@
   fPolarTheta  = rp.fPolarTheta;
   fPolarPhi    = rp.fPolarPhi;
   fParticlePDG = rp.fParticlePDG; 
-  const Int_t nSPECIES = AliESDtrack::kSPECIES;
-  for(Int_t i = 0; i<nSPECIES+4 ; i++)
+  const Int_t nSPECIES = AliESDtrack::kSPECIESN;
+  for(Int_t i = 0; i<nSPECIES ; i++)
     fPID[i]=rp.fPID[i];
 }
 
@@ -132,7 +132,6 @@ const Double_t * AliPHOSRecParticle::GetPID()
   // 8       K0L
 
  
-  const Int_t nSPECIES = AliESDtrack::kSPECIES;
   if (IsElectron()     ) fPID[0] = 1.0;
   if (IsChargedHadron()) {
     fPID[1] = 0.25;
@@ -153,10 +152,11 @@ const Double_t * AliPHOSRecParticle::GetPID()
     fPID[4] = 1.00;
   }
 
-  if (IsPhoton() || IsHardPhoton()) fPID[nSPECIES]  =1.0;
-  if (IsHardPi0())                  fPID[nSPECIES+1]=1.0;
-  if (IsFastNeutralHadron())        fPID[nSPECIES+2]=1.0;
-  if (IsSlowNeutralHadron())        fPID[nSPECIES+3]=1.0;
+  if (IsPhoton() || IsHardPhoton()) fPID[5] = 1.0;
+  if (IsHardPi0())                  fPID[6] = 1.0;
+  if (IsFastNeutralHadron())        fPID[7] = 1.0;
+  if (IsSlowNeutralHadron())        fPID[8] = 1.0;
 
+  if (IsElecCon()) fPID[9] = 1.0;
   return fPID;
 }

diff --git a/PHOS/AliPHOSRecParticle.h b/PHOS/AliPHOSRecParticle.h
index ea400c75d9dc4d2787cac35dd40e6e518b129f1e..7801503e50f12bd10e4b36efef9ad4505f871c92 100644 (file)
--- a/PHOS/AliPHOSRecParticle.h
+++ b/PHOS/AliPHOSRecParticle.h
@@ -49,7 +49,7 @@ class AliPHOSRecParticle : public AliPHOSFastRecParticle {
   Int_t fPHOSTrackSegment ; // pointer to the associated track segment in PHOS  
   Bool_t fDebug ; // to steer debug output
   TVector3 fPos ; // position in the global alice coordinate system 
-  Double_t fPID[AliESDtrack::kSPECIES+4] ; // PID probability densities
+  Double_t fPID[AliESDtrack::kSPECIESN] ; // PID probability densities
 
   ClassDef(AliPHOSRecParticle,3)  // Reconstructed Particle
 };