Updated strong pid algo for Dplus (Giacomo)

diff --git a/PWGHF/vertexingHF/AliRDHFCutsDplustoKpipi.cxx b/PWGHF/vertexingHF/AliRDHFCutsDplustoKpipi.cxx
index 7be2c6a67378379c60e4a4f8b6fde6842263c3eb..173a59103db5177f605900717592266da6d0c54f 100644 (file)
--- a/PWGHF/vertexingHF/AliRDHFCutsDplustoKpipi.cxx
+++ b/PWGHF/vertexingHF/AliRDHFCutsDplustoKpipi.cxx
@@ -40,8 +40,8 @@ AliRDHFCutsDplustoKpipi::AliRDHFCutsDplustoKpipi(const char* name) :
 AliRDHFCuts(name),
   fUseStrongPid(0),
   fMaxPtStrongPid(0.),
-  fMaxPtStrongPidK(0.),
-  fMaxPtStrongPidpi(0.),
+  fMaxPStrongPidK(0.),
+  fMaxPStrongPidpi(0.),
   fUseImpParProdCorrCut(kFALSE)
 {
   //
@@ -125,8 +125,8 @@ AliRDHFCutsDplustoKpipi::AliRDHFCutsDplustoKpipi(const AliRDHFCutsDplustoKpipi &
   AliRDHFCuts(source),
   fUseStrongPid(source.fUseStrongPid),
   fMaxPtStrongPid(source.fMaxPtStrongPid),
-  fMaxPtStrongPidK(source.fMaxPtStrongPidK),
-  fMaxPtStrongPidpi(source.fMaxPtStrongPidpi),
+  fMaxPStrongPidK(source.fMaxPStrongPidK),
+  fMaxPStrongPidpi(source.fMaxPStrongPidpi),
   fUseImpParProdCorrCut(source.fUseImpParProdCorrCut)
 {
   //
@@ -146,8 +146,8 @@ AliRDHFCutsDplustoKpipi &AliRDHFCutsDplustoKpipi::operator=(const AliRDHFCutsDpl
 
   fUseStrongPid=source.fUseStrongPid;
   fMaxPtStrongPid=source.fMaxPtStrongPid;
-  fMaxPtStrongPidK=source.fMaxPtStrongPidK;
-  fMaxPtStrongPidpi=source.fMaxPtStrongPidpi;
+  fMaxPStrongPidK=source.fMaxPStrongPidK;
+  fMaxPStrongPidpi=source.fMaxPStrongPidpi;
   fUseImpParProdCorrCut=source.fUseImpParProdCorrCut;
 
   return *this;
@@ -325,11 +325,11 @@ Int_t AliRDHFCutsDplustoKpipi::IsSelectedPID(AliAODRecoDecayHF *rd)
     if(isKaon<0) nNotKaons++;  
     if(sign==track->Charge()){//pions
       if(isPion<0)return 0;
-      if(rd->Pt()<fMaxPtStrongPid && isPion<=0 && fUseStrongPid&2 && track->Pt()<fMaxPtStrongPidpi)return 0;
+      if(rd->Pt()<fMaxPtStrongPid && isPion<=0 && fUseStrongPid&2 && track->P()<fMaxPStrongPidpi)return 0;
     }
     else{//kaons
       if(isKaon<0)return 0;
-       if(rd->Pt()<fMaxPtStrongPid && isKaon<=0 && fUseStrongPid&1&& track->Pt()<fMaxPtStrongPidK)return 0;
+       if(rd->Pt()<fMaxPtStrongPid && isKaon<=0 && fUseStrongPid&1&& track->P()<fMaxPStrongPidK)return 0;
     }
   }
   

diff --git a/PWGHF/vertexingHF/AliRDHFCutsDplustoKpipi.h b/PWGHF/vertexingHF/AliRDHFCutsDplustoKpipi.h
index 38988e5180532aef31decba07c1f09d83e714fb4..6817f1e5efc4771885a7fd724bd98a563253ed0b 100644 (file)
--- a/PWGHF/vertexingHF/AliRDHFCutsDplustoKpipi.h
+++ b/PWGHF/vertexingHF/AliRDHFCutsDplustoKpipi.h
@@ -49,12 +49,12 @@ class AliRDHFCutsDplustoKpipi : public AliRDHFCuts
   Float_t GetDCACut(Int_t iPtBin=0) const { return (GetCuts() ? fCutsRD[GetGlobalIndex(1,iPtBin)] : 1.e6);}
   void SetUseStrongPid(Int_t spid){fUseStrongPid=spid;}
   void SetMaxPtStrongPid(Float_t spid){fMaxPtStrongPid=spid;}
-  void SetMaxPtStrongPidK(Float_t spid){fMaxPtStrongPidK=spid;}
-  void SetMaxPtStrongPidpi(Float_t spid){fMaxPtStrongPidpi=spid;}
+  void SetMaxPStrongPidK(Float_t spid){fMaxPStrongPidK=spid;}
+  void SetMaxPStrongPidpi(Float_t spid){fMaxPStrongPidpi=spid;}
   Int_t GetStrongPid() const {return fUseStrongPid;}
   Float_t GetMaxPtStrongPid() const {return fMaxPtStrongPid;}
-  Float_t GetMaxPtStrongPidK() const {return fMaxPtStrongPidK;}
-  Float_t GetMaxPtStrongPidpi() const {return fMaxPtStrongPidpi;}
+  Float_t GetMaxPtStrongPidK() const {return fMaxPStrongPidK;}
+  Float_t GetMaxPtStrongPidpi() const {return fMaxPStrongPidpi;}
   void SetUseImpParProdCorrCut(Bool_t use){
     fUseImpParProdCorrCut=use;
   }
@@ -67,11 +67,11 @@ class AliRDHFCutsDplustoKpipi : public AliRDHFCuts
  private:
   Int_t fUseStrongPid; //use strong pid 0 no,1 only for K,2 pi 3 both
   Float_t fMaxPtStrongPid;//Maximum pt of candidate to apply strong Pid
-  Float_t fMaxPtStrongPidK;//Maximum pt of track to apply strong Pid on K
-  Float_t fMaxPtStrongPidpi;//Maximum pt of track to apply strong Pid on pi
+  Float_t fMaxPStrongPidK;//Maximum P of track to apply strong Pid on K
+  Float_t fMaxPStrongPidpi;//Maximum P of track to apply strong Pid on pi
   Bool_t fUseImpParProdCorrCut; //switch for d0K*d0pi1 vs. d0K*d0pi2 cut
 
-  ClassDef(AliRDHFCutsDplustoKpipi,5);  // class for cuts on AOD reconstructed 
+  ClassDef(AliRDHFCutsDplustoKpipi,6);  // class for cuts on AOD reconstructed 
                                    // D+->Kpipi
 };