Added upper p threshold for using TPC PID (Rossella)

diff --git a/PWG3/vertexingHF/AliAODPidHF.cxx b/PWG3/vertexingHF/AliAODPidHF.cxx
index cd48124a53eca3b81af20670f67f45049fbf4470..d17e7a18e3d233a4697c2e510ed66ad1310b1ae9 100644 (file)
--- a/PWG3/vertexingHF/AliAODPidHF.cxx
+++ b/PWG3/vertexingHF/AliAODPidHF.cxx
@@ -56,7 +56,8 @@ AliAODPidHF::AliAODPidHF():
   fMCLowEn2011(kFALSE),
   fPbPb(kFALSE),
   fTOFdecide(kFALSE),
-  fOldPid(kFALSE),
+  fOldPid(kTRUE),
+  fPtThresholdTPC(999999.),
   fPidResponse(0)
 {
  //
@@ -115,6 +116,7 @@ AliAODPidHF::AliAODPidHF(const AliAODPidHF& pid) :
   fPbPb(pid.fPbPb),
   fTOFdecide(pid.fTOFdecide),
   fOldPid(pid.fOldPid),
+  fPtThresholdTPC(pid.fPtThresholdTPC),
   fPidResponse(pid.fPidResponse)
   {
   
@@ -197,6 +199,7 @@ Int_t AliAODPidHF::ApplyPidTPCRaw(AliAODTrack *track,Int_t specie) const{
   
   Double_t dedx=pidObj->GetTPCsignal();
   Double_t mom = pidObj->GetTPCmomentum();
+  if(mom>fPtThresholdTPC) return 0;
   AliTPCPIDResponse tpcResponse;
   SetBetheBloch(tpcResponse); 
   UShort_t nTPCClus=pidObj->GetTPCsignalN();
@@ -515,6 +518,7 @@ Bool_t AliAODPidHF::TPCRawAsym(AliAODTrack* track,Int_t specie) const{
   if(!CheckStatus(track,"TPC")) return kFALSE;
   AliAODPid *pidObj = track->GetDetPid();
   Double_t mom = pidObj->GetTPCmomentum();
+  if(mom>fPtThresholdTPC) return 0;
   Double_t nsigma=999.;
   if(fOldPid){
   Double_t dedx=pidObj->GetTPCsignal();

diff --git a/PWG3/vertexingHF/AliAODPidHF.h b/PWG3/vertexingHF/AliAODPidHF.h
index 7bf55e46ad0b62c2f4e6129b43d9427b5391a3a6..ee27b2d3051fcc0c9a786783c39229a2b7eb2fc6 100644 (file)
--- a/PWG3/vertexingHF/AliAODPidHF.h
+++ b/PWG3/vertexingHF/AliAODPidHF.h
@@ -53,6 +53,7 @@ class AliAODPidHF : public AliAODPid{
  void SetPCompatTOF(Double_t pTOF){fPCompatTOF=pTOF;return;}
  void SetTOFdecide(Bool_t tOFdecide){fTOFdecide=tOFdecide;return;}
  void SetOldPid(Bool_t oldPid){fOldPid=oldPid;return;}
+ void SetPtThresholdTPC(Double_t ptThresholdTPC){fPtThresholdTPC=ptThresholdTPC;return;}
  void SetPidResponse(AliPIDResponse *pidResp) {fPidResponse=pidResp;return;}
  
  //Getters
@@ -78,6 +79,7 @@ class AliAODPidHF : public AliAODPid{
  Double_t GetnSigmaCompatTPC() const{return fnSigmaCompat[0];}
  Double_t GetnSigmaCompatTOF() const{return fnSigmaCompat[1];}
  Bool_t GetOldPid(){return fOldPid;}
+ Double_t GetPtThresholdTPC(){return fPtThresholdTPC;}
  AliPIDResponse *GetPidResponse() const{return fPidResponse;}
 
  Int_t RawSignalPID (AliAODTrack *track, TString detector) const;
@@ -132,11 +134,12 @@ class AliAODPidHF : public AliAODPid{
  Bool_t fPbPb; //  real data PbPb 
  Bool_t fTOFdecide; //  real data PbPb 
  Bool_t fOldPid; //  old PID method implemented
+ Double_t fPtThresholdTPC; //  pT threshold to use TPC PID
  AliPIDResponse *fPidResponse; //pid response
  
 
 
- ClassDef(AliAODPidHF,12) // AliAODPid for heavy flavor PID
+ ClassDef(AliAODPidHF,13) // AliAODPid for heavy flavor PID
 
 };