Upper cut in track p for PID in D* cut class (Alessandro)

diff --git a/PWGHF/vertexingHF/AliRDHFCutsDStartoKpipi.cxx b/PWGHF/vertexingHF/AliRDHFCutsDStartoKpipi.cxx
index 26f1fb86d71304a966afd82ff2571af407152356..499d94f7c6ac88af2aa94a236a13064f78dfcf33 100644 (file)
--- a/PWGHF/vertexingHF/AliRDHFCutsDStartoKpipi.cxx
+++ b/PWGHF/vertexingHF/AliRDHFCutsDStartoKpipi.cxx
@@ -48,7 +48,8 @@ ClassImp(AliRDHFCutsDStartoKpipi)
 AliRDHFCutsDStartoKpipi::AliRDHFCutsDStartoKpipi(const char* name) : 
   AliRDHFCuts(name),
   fTrackCutsSoftPi(0),
-  fMaxPtPid(9999.)
+  fMaxPtPid(9999.),
+  fTPCflag(999.)
 {
   //
   // Default Constructor
@@ -116,7 +117,8 @@ AliRDHFCutsDStartoKpipi::AliRDHFCutsDStartoKpipi(const char* name) :
 AliRDHFCutsDStartoKpipi::AliRDHFCutsDStartoKpipi(const AliRDHFCutsDStartoKpipi &source) :
   AliRDHFCuts(source),
   fTrackCutsSoftPi(0),
-  fMaxPtPid(9999.)
+  fMaxPtPid(9999.),
+  fTPCflag(999.)
 {
   //
   // Copy constructor
@@ -523,6 +525,13 @@ Int_t AliRDHFCutsDStartoKpipi::SelectPID(AliAODTrack *track, Int_t type)
        if(type==3) isTOF=fPidHF->IsKaonRaw(track,"TOF");
       }
     }
+
+    //--------------------------------
+    // cut on high momentum in the TPC
+    //--------------------------------
+    Double_t pPIDcut = track->P();
+    if(pPIDcut>fTPCflag) isTPC=1;
+    
     isParticle = isTPC&&isTOF;
   }
   
@@ -541,6 +550,7 @@ Int_t AliRDHFCutsDStartoKpipi::SelectPID(AliAODTrack *track, Int_t type)
     isParticle = Bool_t(k==type);
   }
   
+
   return isParticle;
   
 }

diff --git a/PWGHF/vertexingHF/AliRDHFCutsDStartoKpipi.h b/PWGHF/vertexingHF/AliRDHFCutsDStartoKpipi.h
index b78234f9323e2430f7b7050df9088cf809b64381..c6bb9782ce791eb177864bd9ccd6a73a81037864 100644 (file)
--- a/PWGHF/vertexingHF/AliRDHFCutsDStartoKpipi.h
+++ b/PWGHF/vertexingHF/AliRDHFCutsDStartoKpipi.h
@@ -43,6 +43,8 @@ class AliRDHFCutsDStartoKpipi : public AliRDHFCuts
   virtual void SetStandardCutsPbPb2011();  
   void SetMaxPtPid(Float_t maxPt){fMaxPtPid = maxPt;}
 
+  void SetOffHighPtPIDinTPC(Float_t TPCrem =999.){fTPCflag = TPCrem;}
+
   void AddTrackCutsSoftPi(const AliESDtrackCuts *cuts) 
      {fTrackCutsSoftPi=new AliESDtrackCuts(*cuts); return;}
   virtual AliESDtrackCuts *GetTrackCutsSoftPi() const {return fTrackCutsSoftPi;}
@@ -51,8 +53,9 @@ class AliRDHFCutsDStartoKpipi : public AliRDHFCuts
 
   AliESDtrackCuts *fTrackCutsSoftPi; // cuts for soft pion (AOD converted to ESD on the flight!)
   Float_t fMaxPtPid; // maximum Dstar Pt for using PID
+  Float_t fTPCflag;   // 
 
-  ClassDef(AliRDHFCutsDStartoKpipi,5);  // class for cuts on AOD reconstructed D0->Kpipi
+  ClassDef(AliRDHFCutsDStartoKpipi,6);  // class for cuts on AOD reconstructed D0->Kpipi
 };
 
 #endif