Setter for the maxpt for the tight PID selection (Andrea)

diff --git a/PWG3/vertexingHF/AliRDHFCutsD0toKpi.cxx b/PWG3/vertexingHF/AliRDHFCutsD0toKpi.cxx
index 1d89fef6808e68ed4fd921d6afde48d8df70644c..596e0e482ad2c11a8b5ade8cecb5ceab4fb63999 100644 (file)
--- a/PWG3/vertexingHF/AliRDHFCutsD0toKpi.cxx
+++ b/PWG3/vertexingHF/AliRDHFCutsD0toKpi.cxx
@@ -43,7 +43,8 @@ AliRDHFCuts(name),
 fUseSpecialCuts(kTRUE),
 fLowPt(kTRUE),
 fDefaultPID(kFALSE),
-fUseKF(kFALSE)
+fUseKF(kFALSE),
+fPtLowPID(2.)
 {
   //
   // Default Constructor
@@ -95,7 +96,8 @@ AliRDHFCutsD0toKpi::AliRDHFCutsD0toKpi(const AliRDHFCutsD0toKpi &source) :
   fUseSpecialCuts(source.fUseSpecialCuts),
   fLowPt(source.fLowPt),
   fDefaultPID(source.fDefaultPID),
-  fUseKF(source.fUseKF)
+  fUseKF(source.fUseKF),
+  fPtLowPID(source.fPtLowPID)
 {
   //
   // Copy constructor
@@ -114,6 +116,7 @@ AliRDHFCutsD0toKpi &AliRDHFCutsD0toKpi::operator=(const AliRDHFCutsD0toKpi &sour
   fUseSpecialCuts=source.fUseSpecialCuts;
   fLowPt=source.fLowPt;
   fDefaultPID=source.fDefaultPID;
+  fPtLowPID=source.fPtLowPID;
 
   return *this;
 }
@@ -727,7 +730,7 @@ Int_t AliRDHFCutsD0toKpi::IsSelectedPID(AliAODRecoDecayHF* d)
       else isD0D0barPID[1]=0;// not a D0bar if pi+ or if K+ excluded
     }
 
-    if(fLowPt && d->Pt()<2.){
+    if(fLowPt && d->Pt()<fPtLowPID){
      Double_t sigmaTPC[3]={3.,2.,0.};
      fPidHF->SetSigmaForTPC(sigmaTPC);
     // identify kaon
@@ -776,10 +779,10 @@ Int_t AliRDHFCutsD0toKpi::IsSelectedPID(AliAODRecoDecayHF* d)
     return 0;
   }
 
-  if(fLowPt && d->Pt()<2.){
-    if(fLowPt) fPidHF->SetSigmaForTPC(sigma_tmp);
+  if(fLowPt && d->Pt()<fPtLowPID){    
     if(combinedPID[0][0]<=0&&combinedPID[1][0]<=0){
       fWhyRejection=32;// reject cases where the Kaon is not identified
+      fPidHF->SetSigmaForTPC(sigma_tmp);
       return 0;
     }
   }

diff --git a/PWG3/vertexingHF/AliRDHFCutsD0toKpi.h b/PWG3/vertexingHF/AliRDHFCutsD0toKpi.h
index ad3cb34d99328bad233071d5bcb9d376bb7c1b57..d346e2127c5e9ddfc70017bd72ba0b7a6eb1f08e 100644 (file)
--- a/PWG3/vertexingHF/AliRDHFCutsD0toKpi.h
+++ b/PWG3/vertexingHF/AliRDHFCutsD0toKpi.h
@@ -49,10 +49,11 @@ class AliRDHFCutsD0toKpi : public AliRDHFCuts
   Int_t IsSelectedPIDdefault(AliAODRecoDecayHF *rd);
   Int_t IsSelectedSpecialCuts(AliAODRecoDecayHF *d) const;
   void SetUseSpecialCuts(Bool_t useSpecialCuts) {fUseSpecialCuts=useSpecialCuts;}
-  void SetLowPt(Bool_t lowpt) {fLowPt=lowpt;}
+  void SetLowPt(Bool_t lowpt,Double_t ptlow=2.) {fLowPt=lowpt;fPtLowPID=ptlow;}
   Bool_t GetUseSpecialCuts() const {return fUseSpecialCuts;}
   void SetUseDefaultPID(Bool_t defPID){fDefaultPID=defPID;}
   Bool_t GetIsUsedDefPID(){return fDefaultPID;}
+  Double_t GetPtForPIDtight()const {return fPtLowPID;}
   void SetUseKF(Bool_t useKF);
   Bool_t GetIsUsedKF() const {return fUseKF;}
   
@@ -64,8 +65,9 @@ class AliRDHFCutsD0toKpi : public AliRDHFCuts
   Bool_t fLowPt;           // flag to switch on/off different pid for low pt D0
   Bool_t fDefaultPID;      // flag to switch on/off the default pid
   Bool_t fUseKF;           // flag to switch on/off D0 selection via KF 
+  Double_t fPtLowPID;      // transverse momentum below which the strong PID is applied
 
-  ClassDef(AliRDHFCutsD0toKpi,5);  // class for cuts on AOD reconstructed D0->Kpi
+  ClassDef(AliRDHFCutsD0toKpi,6);  // class for cuts on AOD reconstructed D0->Kpi
 };
 
 #endif