use new GetTPCncls call on AOD

diff --git a/PWG/EMCAL/AliEmcalPicoTrackMaker.cxx b/PWG/EMCAL/AliEmcalPicoTrackMaker.cxx
index ff992bb8968ae155fa40c2e6c2400f1cf51e4d41..3ce26bb9fdf50d360c9f53f3952fc734f6dd390b 100644 (file)
--- a/PWG/EMCAL/AliEmcalPicoTrackMaker.cxx
+++ b/PWG/EMCAL/AliEmcalPicoTrackMaker.cxx
@@ -34,6 +34,7 @@ AliEmcalPicoTrackMaker::AliEmcalPicoTrackMaker() :
   fIncludeNoITS(kTRUE),
   fUseNegativeLabels(kTRUE),
   fIsMC(kFALSE),
+  fCutMaxFractionSharedTPCClusters(0.4),
   fTracksIn(0),
   fTracksOut(0)
 {
@@ -59,6 +60,7 @@ AliEmcalPicoTrackMaker::AliEmcalPicoTrackMaker(const char *name) :
   fIncludeNoITS(kTRUE),
   fUseNegativeLabels(kTRUE),
   fIsMC(kFALSE),
+  fCutMaxFractionSharedTPCClusters(0.4),
   fTracksIn(0),
   fTracksOut(0)
 {
@@ -172,6 +174,11 @@ void AliEmcalPicoTrackMaker::UserExec(Option_t *)
          continue;
        }
       }
+      if (fCutMaxFractionSharedTPCClusters > 0) {
+       Double_t frac = Double_t(aodtrack->GetTPCnclsS()) / Double_t(aodtrack->GetTPCncls());
+       if (frac > fCutMaxFractionSharedTPCClusters) 
+         continue;
+      }
       if (TMath::Abs(track->GetTrackEtaOnEMCal()) < 0.75 && 
          track->GetTrackPhiOnEMCal() > 70 * TMath::DegToRad() &&
          track->GetTrackPhiOnEMCal() < 190 * TMath::DegToRad())

diff --git a/PWG/EMCAL/AliEmcalPicoTrackMaker.h b/PWG/EMCAL/AliEmcalPicoTrackMaker.h
index 9ef049380dbd7f708abb18412632ff8971d0f02b..11271d34991eb84446842c194b21b0128fcc2a8f 100644 (file)
--- a/PWG/EMCAL/AliEmcalPicoTrackMaker.h
+++ b/PWG/EMCAL/AliEmcalPicoTrackMaker.h
@@ -30,6 +30,7 @@ class AliEmcalPicoTrackMaker : public AliAnalysisTaskSE {
   void SetIncludeNoITS(Bool_t f)                         { fIncludeNoITS     = f   ; }
   void SetUseNegativeLabels(Bool_t f)                    { fUseNegativeLabels= f   ; }
   void SetMC(Bool_t a)                                   { fIsMC             = a   ; }
+  void SetCutMaxFractionSharedTPCClusters(Double_t c = 0.4)        { fCutMaxFractionSharedTPCClusters  = c ; }
 
  protected:
   Int_t              fAODfilterBits[2];     // AOD track filter bit map
@@ -46,6 +47,7 @@ class AliEmcalPicoTrackMaker : public AliAnalysisTaskSE {
   Bool_t             fIncludeNoITS;         // includes tracks with failed ITS refit
   Bool_t             fUseNegativeLabels;    // whether or not should use negative MC labels
   Bool_t             fIsMC;                 // whether it is a MC event or not
+  Double_t           fCutMaxFractionSharedTPCClusters;  // max fraction of shared TPC clusters
   TClonesArray      *fTracksIn;             //!track array in
   TClonesArray      *fTracksOut;            //!track array out