Possibility for more selective PID (Giacomo)

diff --git a/PWG3/vertexingHF/AliRDHFCutsDplustoKpipi.cxx b/PWG3/vertexingHF/AliRDHFCutsDplustoKpipi.cxx
index cb7b882a4647e1b65fc8d28bcdb7b84e41cd4c1a..b1524ac7f9f40b7f83e2222265dce629f582ee37 100644 (file)
--- a/PWG3/vertexingHF/AliRDHFCutsDplustoKpipi.cxx
+++ b/PWG3/vertexingHF/AliRDHFCutsDplustoKpipi.cxx
@@ -34,7 +34,8 @@ ClassImp(AliRDHFCutsDplustoKpipi)
 
 //--------------------------------------------------------------------------
 AliRDHFCutsDplustoKpipi::AliRDHFCutsDplustoKpipi(const char* name) : 
-  AliRDHFCuts(name)
+AliRDHFCuts(name),
+fUseStrongPid(kFALSE)
 {
   //
   // Default Constructor
@@ -100,7 +101,9 @@ AliRDHFCutsDplustoKpipi::AliRDHFCutsDplustoKpipi(const char* name) :
 }
 //--------------------------------------------------------------------------
 AliRDHFCutsDplustoKpipi::AliRDHFCutsDplustoKpipi(const AliRDHFCutsDplustoKpipi &source) :
-  AliRDHFCuts(source)
+  AliRDHFCuts(source),
+  fUseStrongPid(source.fUseStrongPid)
+
 {
   //
   // Copy constructor
@@ -269,9 +272,11 @@ Int_t AliRDHFCutsDplustoKpipi::IsSelectedPID(AliAODRecoDecayHF *rd)
     if(isKaon<0) nNotKaons++;  
     if(sign==track->Charge()){//pions
       if(isPion<0)return 0;
+      if(rd->Pt()<2. && isPion<=0 && fUseStrongPid)return 0;
     }
       else{//kaons
        if(isKaon<0)return 0;
+       if(rd->Pt()<2. && isKaon<=0 && fUseStrongPid)return 0;
       }
     
       

diff --git a/PWG3/vertexingHF/AliRDHFCutsDplustoKpipi.h b/PWG3/vertexingHF/AliRDHFCutsDplustoKpipi.h
index 4624f21baa0560f3b73440134b9561501eac400c..718f212162c475f312724bd9b74ace248e5574ed 100644 (file)
--- a/PWG3/vertexingHF/AliRDHFCutsDplustoKpipi.h
+++ b/PWG3/vertexingHF/AliRDHFCutsDplustoKpipi.h
@@ -38,10 +38,14 @@ class AliRDHFCutsDplustoKpipi : public AliRDHFCuts
 
   Float_t GetMassCut(Int_t iPtBin=0) const { return (GetCuts() ? fCutsRD[GetGlobalIndex(0,iPtBin)] : 1.e6);}
   Float_t GetDCACut(Int_t iPtBin=0) const { return (GetCuts() ? fCutsRD[GetGlobalIndex(1,iPtBin)] : 1.e6);}
+  void SetUseStrongPid(Bool_t spid){fUseStrongPid=spid;}
+  Bool_t GetStrongPid(){return fUseStrongPid;}
   
  protected:
 
-  ClassDef(AliRDHFCutsDplustoKpipi,1);  // class for cuts on AOD reconstructed 
+ private:
+  Bool_t fUseStrongPid;//use strong pid
+  ClassDef(AliRDHFCutsDplustoKpipi,2);  // class for cuts on AOD reconstructed 
                                    // D+->Kpipi
 };