Upper cut in track p for PID in D* cut class (Alessandro)
authorfprino <fprino@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 22 Nov 2012 00:22:06 +0000 (00:22 +0000)
committerfprino <fprino@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 22 Nov 2012 00:22:06 +0000 (00:22 +0000)
PWGHF/vertexingHF/AliRDHFCutsDStartoKpipi.cxx
PWGHF/vertexingHF/AliRDHFCutsDStartoKpipi.h

index 26f1fb8..499d94f 100644 (file)
@@ -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;
   
 }
index b78234f..c6bb978 100644 (file)
@@ -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