Add flag to apply shared cluster cut on AOD
authorMarco van Leeuwen <marco.van.leeuwen@cern.ch>
Fri, 6 Dec 2013 13:29:50 +0000 (14:29 +0100)
committerMarco van Leeuwen <marco.van.leeuwen@cern.ch>
Fri, 6 Dec 2013 13:29:50 +0000 (14:29 +0100)
JETAN/AliAnalysisTaskJetCluster.cxx
JETAN/AliAnalysisTaskJetCluster.h

index 270b614..6e5f04b 100644 (file)
@@ -119,6 +119,7 @@ AliAnalysisTaskJetCluster::AliAnalysisTaskJetCluster():
   fExternalWeight(1),
   fTrackEtaWindow(0.9),    
   fRequireITSRefit(0),
+  fApplySharedClusterCut(0),
   fRecEtaWindow(0.5),
   fTrackPtCut(0.),                                                     
   fJetOutputMinPt(0.150),
@@ -287,6 +288,7 @@ AliAnalysisTaskJetCluster::AliAnalysisTaskJetCluster(const char* name):
   fExternalWeight(1),    
   fTrackEtaWindow(0.9),    
   fRequireITSRefit(0),
+  fApplySharedClusterCut(0),
   fRecEtaWindow(0.5),
   fTrackPtCut(0.),                                                     
   fJetOutputMinPt(0.150),
@@ -1888,6 +1890,10 @@ Int_t  AliAnalysisTaskJetCluster::GetListOfTracks(TList *list,Int_t type){
          continue;
        }
         if(fRequireITSRefit){if((tr->GetStatus()&AliESDtrack::kITSrefit)==0)continue;}
+        if (fApplySharedClusterCut) {
+           Double_t frac = Double_t(tr->GetTPCnclsS()) /Double_t(tr->GetTPCncls());
+           if (frac > 0.4) continue;
+        } 
        if(TMath::Abs(tr->Eta())>fTrackEtaWindow){
          if(fDebug>10)Printf("%s:%d Not matching eta %d/%d",(char*)__FILE__,__LINE__,it,aod->GetNumberOfTracks());     
          continue;
index 7b21770..849f837 100644 (file)
@@ -63,6 +63,7 @@ class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE
     virtual void SetAODMCInput(Bool_t b){fUseAODMCInput = b;}
     virtual void SetEventSelection(Bool_t b){fEventSelection = b;}
     virtual void SetRequireITSRefit(Int_t i){fRequireITSRefit=i;}
+    virtual void SetSharedClusterCut(Int_t docut){fApplySharedClusterCut=docut;}
     virtual void SetRecEtaWindow(Float_t f){fRecEtaWindow = f;}
     virtual void SetTrackEtaWindow(Float_t f){fTrackEtaWindow = f;}
     virtual void SetTrackTypeGen(Int_t i){fTrackTypeGen = i;}
@@ -181,6 +182,7 @@ class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE
     Float_t       fExternalWeight;        // external weight
     Float_t       fTrackEtaWindow;        // eta window used for corraltion plots between rec and gen 
     Int_t         fRequireITSRefit;       // to select hybrids with ITS refit only
+    Int_t         fApplySharedClusterCut; // flag to apply shared cluster cut (needed for some AODs where this cut was not applied in the filtering)
     Float_t       fRecEtaWindow;          // eta window used for corraltion plots between rec and gen 
     Float_t       fTrackPtCut;            // minimum track pt to be accepted
     Float_t       fJetOutputMinPt;        // minimum p_t for jets to be written out
@@ -344,7 +346,7 @@ class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE
     TList *fHistList; //!leading tracks to be skipped in the randomized event Output list
    
 
-    ClassDef(AliAnalysisTaskJetCluster, 25) 
+    ClassDef(AliAnalysisTaskJetCluster, 26) 
 };
  
 #endif