Add number of points used for dE/dx in TPC to AliAODpid

 https://savannah.cern.ch/bugs/?74184

 (M.v. Leeuwen)

diff --git a/ANALYSIS/AliAnalysisTaskESDfilter.cxx b/ANALYSIS/AliAnalysisTaskESDfilter.cxx
index 6bb3c9cfca59bce92bc4550a94925f8a6ead6497..65d0fc80807acad525043be0bcdf572477f7ba08 100644 (file)
--- a/ANALYSIS/AliAnalysisTaskESDfilter.cxx
+++ b/ANALYSIS/AliAnalysisTaskESDfilter.cxx
@@ -1304,6 +1304,7 @@ void AliAnalysisTaskESDfilter::SetDetectorRawSignals(AliAODPid *aodpid, AliESDtr
 \r
  aodpid->SetITSsignal(track->GetITSsignal());\r
  aodpid->SetTPCsignal(track->GetTPCsignal());\r
+ aodpid->SetTPCsignalN(track->GetTPCsignalN());\r
 \r
  //n TRD planes = 6\r
  Int_t nslices = track->GetNumberOfTRDslices()*6;\r

diff --git a/STEER/AliAODPid.h b/STEER/AliAODPid.h
index 8d8bde700b5fa07d1ad1e7adf3328968b066ad79..37bf8703287326996ca6c8703e5fd1492752f1c3 100644 (file)
--- a/STEER/AliAODPid.h
+++ b/STEER/AliAODPid.h
@@ -25,6 +25,7 @@ class AliAODPid : public TObject {
  //setters
   void      SetITSsignal(Double_t its)                         {fITSsignal=its;}
   void      SetTPCsignal(Double_t tpc)                         {fTPCsignal=tpc;}
+  void      SetTPCsignalN(UShort_t tpcN)                       {fTPCsignalN=tpcN;}
   void      SetTPCmomentum(Double_t tpcMom)                    {fTPCmomentum=tpcMom;}
   void      SetTRDsignal(Int_t nslices, Double_t* trdslices)   {fTRDnSlices = nslices; fTRDslices=trdslices;}
   void      SetTRDmomentum(Int_t nplane, Float_t trdMom)       {fTRDmomentum[nplane]=trdMom;}
@@ -36,6 +37,7 @@ class AliAODPid : public TObject {
 
   Double_t  GetITSsignal()       const {return  fITSsignal;}
   Double_t  GetTPCsignal()       const {return  fTPCsignal;}
+  UShort_t  GetTPCsignalN()      const {return  fTPCsignalN;}
   Double_t  GetTPCmomentum()     const {return  fTPCmomentum;}
   Int_t     GetTRDnSlices()      const {return  fTRDnSlices;}
   Double_t* GetTRDsignal()       const {return  fTRDslices;}
@@ -50,6 +52,7 @@ class AliAODPid : public TObject {
  private :
   Double32_t  fITSsignal;        //[0.,0.,10] detector raw signal
   Double32_t  fTPCsignal;        //[0.,0.,10] detector raw signal
+  UShort_t    fTPCsignalN;       // number of points used for TPC dE/dx
   Double_t    fTPCmomentum;      // momentum at the inner wall of TPC;
   Int_t       fTRDnSlices;       // N slices used for PID in the TRD
   Double32_t* fTRDslices;        //[fTRDnSlices]
@@ -62,7 +65,7 @@ class AliAODPid : public TObject {
   Double32_t  fEMCALMomentum[3]; // momentum of track
                                 // extrapolated to EMCAL surface
 
-  ClassDef(AliAODPid, 3);
+  ClassDef(AliAODPid, 4);
 };
 
 #endif

diff --git a/STEER/AliAODpidUtil.cxx b/STEER/AliAODpidUtil.cxx
index 62cfe6952caa5bb05b76dcde75598a59fc843e40..8e8b837b404190e86ddadb52f904de293a0024f3 100644 (file)
--- a/STEER/AliAODpidUtil.cxx
+++ b/STEER/AliAODpidUtil.cxx
@@ -72,11 +72,11 @@ void AliAODpidUtil::MakeTPCPID(AliAODTrack *track,Double_t *p) const
   //
 
   if ((track->GetStatus()&AliESDtrack::kTPCin )==0) return;
-  UShort_t nTPCClus=track->GetTPCNcls();
 
   Double_t mom = track->P();
   AliAODPid *pidObj = track->GetDetPid();
   if (pidObj) mom = pidObj->GetTPCmomentum();
+  UShort_t nTPCClus=pidObj->GetTPCsignalN();
    
   Double_t dedx=pidObj->GetTPCsignal(); 
   Bool_t mismatch=kTRUE;

diff --git a/STEER/AliAODpidUtil.h b/STEER/AliAODpidUtil.h
index bdbc8b2d2b001cc5e2da3b8c4ef238c85c6f54a0..a6ced09b1242eee92060332ecce5b953b873593c 100644 (file)
--- a/STEER/AliAODpidUtil.h
+++ b/STEER/AliAODpidUtil.h
@@ -56,9 +56,9 @@ private:
 
 inline Float_t AliAODpidUtil::NumberOfSigmasTPC(const AliAODTrack *track, AliPID::EParticleType type) const {
   
-  UShort_t nTPCClus=track->GetTPCNcls();
   Double_t mom = track->P();
   AliAODPid *pidObj = track->GetDetPid();
+  UShort_t nTPCClus=pidObj->GetTPCsignalN();
   if (pidObj)
     mom = pidObj->GetTPCmomentum();
   return fTPCResponse.GetNumberOfSigmas(mom,pidObj->GetTPCsignal(),nTPCClus,type);